肌電信號分析及動作識別畢業(yè)設(shè)計(編輯修改稿)
2025-01-08 15:25 本頁面
【文章內(nèi)容簡介】 ?????? ???? ?? dPdttTES TTT )(2 1)(1 2 2 2lim ??t? 的功率譜密度為 ? ???P TEP TT ])([)(2lim ???? ??? 這樣的平均功率等于各個頻率分量（統(tǒng)計值）單獨貢獻(xiàn)出的 11 功率之連續(xù)和，是在頻率域上描述隨機過程統(tǒng)計特性的最主要數(shù)字特征。 12 倒譜分析 ： 信號的傅里葉變換譜經(jīng)對數(shù)運算后再進(jìn)行的傅里葉反變換。功率譜的對數(shù)值的逆傅氏變換稱為倒譜。 如序列倒譜 c? 定義如下： ? ?? ?? ?xFFc log? 1?? 其中 F 和 F1分別表示傅里葉變換和逆傅里葉變換。 13 第 5章 MATLAB軟件實現(xiàn) MATLAB簡介 MATLAB 是一種面向工程和科學(xué)計算的交互式計算軟件，它以矩陣運算為基礎(chǔ)，把計算、可視化、程序設(shè)計融合到了一個簡單易用的交互式工作環(huán)境中。同時由于 MATLAB是一個數(shù)據(jù)分 析和處理功能十分強大的工程實用軟件，它的信號處理與分析工具箱為語音信號分析提供了十分豐富的功能函數(shù) ,利用這些功能函數(shù)可以快捷而又方便地完成語音信號的處理和分析以及信號的可視化，使人機交互更加便捷。 MATLAB是解決工程技術(shù)問題的技術(shù)平臺。利用它能夠輕松完成復(fù)雜的數(shù)值計算，數(shù)據(jù)分析，符號計算和數(shù)據(jù)可視化等任務(wù)。 MATLAB軟件由主包和各類工具箱構(gòu)成。其中，主包基本是一個用 C/C++等語言編寫成的函數(shù)庫。該函數(shù)庫提供矩陣（或數(shù)組）的各種算法以及建立在此基礎(chǔ)上的各種應(yīng)用函數(shù)和一些相關(guān)的用戶有好操作界面。而工具箱 從深度和廣度上大大擴(kuò)展了 MATLAB 主包的功能和應(yīng)用領(lǐng)域。隨著自身的不斷完善和發(fā)展，MATLAB功能越來越強大，應(yīng)用也越來越廣泛。 隨著計算機技術(shù)和信息技術(shù)的發(fā)展，語音交互已經(jīng)成為人機交互的必要手段，而語音信號的采集和處理是人機交互的前提和基礎(chǔ)。聲卡是計算機對語音信號進(jìn)行加工的重要部件，它具有對信號濾波、放大、采樣保持、A/D 和 D/A轉(zhuǎn)換等功能。在 MATLAB 環(huán)境中 ,可以通過多種編程方法驅(qū)動聲卡，實現(xiàn)對語音信號的采集和回放。同時由于 MATLAB是一個數(shù)據(jù)分析和處理功能十分強大的工程實用軟件，它的信 號處理與分析工具箱為語音信號分析提供了十分豐富的功能函數(shù) ,利用這些功能函數(shù)可以快捷而又方便地完成語音信號的處理和分析以及信號的可視化，使人機交互更加便捷。 14 Matlab程序及 GUI function varargout = untitled20(varargin) gui_Singleton = 1。 gui_State = struct(39。gui_Name39。, mfilename, ... 39。gui_Singleton39。, gui_Singleton, ... 39。gui_OpeningF39。, @untitled20_OpeningF, ... 39。gui_OutputF39。, @untitled20_OutputF, ... 39。gui_LayoutF39。, [] , ... 39。gui_Callback39。, [])。 if nargin amp。amp。 ischar(varargin{1}) = str2func(varargin{1})。 end if nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。 else gui_mainf(gui_State, varargin{:})。 end function untitled20_OpeningF(hObject, eventdata, handles, varargin) = hObject。 guidata(hObject, handles)。 function varargout = untitled20_OutputF(hObject, eventdata, handles) varargout{1} = 。 function qu_Callback(hObject, eventdata, handles) a=load(39。39。)。 axes()。 15 plot(a)。 xlabel(39。時間 (s)39。)。 ylabel(39。被測變量 y39。)。 title(39。(屈 )原始信號 39。)。 grid on。 A=mean(a)。 % 平均值 B=std(a)。 % 標(biāo)準(zhǔn)差 C=Var(a)。 % 方差 D=mean(abs(a))。 % 積分機電值 E=sqrt(mean(a.^2))。%均方根有效值 RMS y=fft(a,10000)。 %做 10000點傅里葉變換 fs=1000。 N=length(y)。 mag=abs(y)。 f=(0:N1)/N*fs。 power=(mag.^2)/10000。%周期圖法功率譜 ss=sum(power)。 M2=*ss。 df=fs/N。 M1=*df*(sum(power(1:N1))+sum(power(2:N)))。 F=M1/M2。 G=M2/2。 A1=A*100000。 A2=round(A1)。 A3=A2/100000。 set(,39。string39。,num2str(A3))。 B1=B*100000。 B2=round(B1)。 B3=B2/100000
公司管理相關(guān)推薦
鄂ICP備17016276號-1