肌電信號分析及動作識別畢業(yè)設計(編輯修改稿)
2025-01-08 15:25 本頁面
【文章內容簡介】 ?????? ???? ?? dPdttTES TTT )(2 1)(1 2 2 2lim ??t? 的功率譜密度為 ? ???P TEP TT ])([)(2lim ???? ??? 這樣的平均功率等于各個頻率分量（統計值）單獨貢獻出的 11 功率之連續(xù)和，是在頻率域上描述隨機過程統計特性的最主要數字特征。 12 倒譜分析 ： 信號的傅里葉變換譜經對數運算后再進行的傅里葉反變換。功率譜的對數值的逆傅氏變換稱為倒譜。 如序列倒譜 c? 定義如下： ? ?? ?? ?xFFc log? 1?? 其中 F 和 F1分別表示傅里葉變換和逆傅里葉變換。 13 第 5章 MATLAB軟件實現 MATLAB簡介 MATLAB 是一種面向工程和科學計算的交互式計算軟件，它以矩陣運算為基礎，把計算、可視化、程序設計融合到了一個簡單易用的交互式工作環(huán)境中。同時由于 MATLAB是一個數據分 析和處理功能十分強大的工程實用軟件，它的信號處理與分析工具箱為語音信號分析提供了十分豐富的功能函數 ,利用這些功能函數可以快捷而又方便地完成語音信號的處理和分析以及信號的可視化，使人機交互更加便捷。 MATLAB是解決工程技術問題的技術平臺。利用它能夠輕松完成復雜的數值計算，數據分析，符號計算和數據可視化等任務。 MATLAB軟件由主包和各類工具箱構成。其中，主包基本是一個用 C/C++等語言編寫成的函數庫。該函數庫提供矩陣（或數組）的各種算法以及建立在此基礎上的各種應用函數和一些相關的用戶有好操作界面。而工具箱 從深度和廣度上大大擴展了 MATLAB 主包的功能和應用領域。隨著自身的不斷完善和發(fā)展，MATLAB功能越來越強大，應用也越來越廣泛。 隨著計算機技術和信息技術的發(fā)展，語音交互已經成為人機交互的必要手段，而語音信號的采集和處理是人機交互的前提和基礎。聲卡是計算機對語音信號進行加工的重要部件，它具有對信號濾波、放大、采樣保持、A/D 和 D/A轉換等功能。在 MATLAB 環(huán)境中 ,可以通過多種編程方法驅動聲卡，實現對語音信號的采集和回放。同時由于 MATLAB是一個數據分析和處理功能十分強大的工程實用軟件，它的信 號處理與分析工具箱為語音信號分析提供了十分豐富的功能函數 ,利用這些功能函數可以快捷而又方便地完成語音信號的處理和分析以及信號的可視化，使人機交互更加便捷。 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)。 % 標準差 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
公司管理相關推薦
鄂ICP備17016276號-1