完整版肌電信號(hào)的識(shí)別

1、燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)燕山大學(xué)課程設(shè)計(jì)說明書題目：肌電信號(hào)分析及動(dòng)作識(shí)別學(xué)院（系）：電氣工程學(xué)院年級(jí)專業(yè)：10級(jí)儀表三班學(xué)號(hào)：學(xué)生姓名：指導(dǎo)教師：教師職稱：教授講師燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)電氣工程學(xué)院課程設(shè)計(jì)任務(wù)書基層教學(xué)單位：自動(dòng)化儀表系指導(dǎo)教師：謝平張淑清學(xué)號(hào)學(xué)生姓名（專業(yè)）班級(jí)10級(jí)儀表三班設(shè)計(jì)題目肌電信號(hào)分析及動(dòng)作識(shí)別設(shè)計(jì)技術(shù)參數(shù)1、獨(dú)立完成設(shè)計(jì)任務(wù)。2、編程，上機(jī)調(diào)試。3、連接硬件實(shí)驗(yàn)線路，實(shí)現(xiàn)所要求的功能。4、完成設(shè)計(jì)，提交課程設(shè)計(jì)報(bào)告。設(shè)計(jì)要求1、利用肌電傳感器采集人體肌電信號(hào)，利用Matlab軟件對(duì)肌電信號(hào)進(jìn)行時(shí)域和頻域（包括FFT,功率譜，倒譜）分析，計(jì)算均值、標(biāo)準(zhǔn)

2、差、方差、積分肌電值IEMG、均方根RMS等時(shí)域指標(biāo)參數(shù)和平均功率頻率MPF、中值頻率MF等頻域指標(biāo)。比較不同動(dòng)作下肌電信號(hào)特征，根據(jù)設(shè)定閾值給出動(dòng)作狀態(tài)，設(shè)計(jì)相應(yīng)的信號(hào)分析及顯示界面。2、基于肌電信號(hào)分析結(jié)果，通過串口發(fā)送命令給單片機(jī)系統(tǒng)，根據(jù)肌電信號(hào)動(dòng)作狀態(tài)控制相應(yīng)的數(shù)碼管或LED發(fā)光顯示。3、擴(kuò)展：也可通過動(dòng)作模式驅(qū)動(dòng)電機(jī)轉(zhuǎn)動(dòng)或其他控制輸出模塊。工作量軟件編程與硬件調(diào)試相結(jié)合，繪制設(shè)計(jì)流程圖，編制相應(yīng)軟件界面，實(shí)現(xiàn)單片機(jī)控制與信號(hào)處理任務(wù)的綜合應(yīng)用參考資料1）微型計(jì)算機(jī)控制系統(tǒng)賴壽宏，機(jī)械工業(yè)出版社（教材）2）單片機(jī)及應(yīng)用李大友，高等教育出版社（教材）3）信號(hào)處理原理及應(yīng)用謝平等機(jī)械工業(yè)

3、出版社（教材）4）«Matlab程序設(shè)計(jì)及其在信號(hào)處理中的應(yīng)用聶祥飛等西南交通大學(xué)出版社周次A周第二周應(yīng)完成內(nèi)容熟悉偉福單片機(jī)編程環(huán)境，調(diào)試單片機(jī)各基本功能模塊；熟悉matlab信號(hào)處理工具箱，信號(hào)處理系統(tǒng)基本功能模塊學(xué)習(xí)和調(diào)試單片機(jī)系統(tǒng)與信號(hào)處理系統(tǒng)綜合進(jìn)行硬件調(diào)試，撰寫課程設(shè)計(jì)報(bào)告指導(dǎo)教師簽字基層教學(xué)單位主任簽字2燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)第一章摘要2第二章系統(tǒng)總體設(shè)計(jì)方案3第三章肌電信號(hào)的時(shí)域參數(shù)處理及其分析4第四章肌電信號(hào)的頻域處理方法及其分析73.1 FFT分析73.2 功率譜分析83.3 倒譜分析93.4 平均功率頻率MPFffi中值頻率10第五章Matlab程序及GU

4、I11第六章系統(tǒng)整體調(diào)試及結(jié)果說明24第七章學(xué)習(xí)心得24參考文獻(xiàn)25燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)第一章摘要肌電信號(hào)是產(chǎn)生肌肉力的電信號(hào)根源，它是肌肉中很多運(yùn)動(dòng)單元?jiǎng)幼麟娢辉跁r(shí)間和空間上的疊加，反映了神經(jīng)，肌肉的功能狀態(tài)，在基礎(chǔ)醫(yī)學(xué)研究、臨床診斷和康復(fù)工程中有廣泛的應(yīng)用。其種類重要有兩種：一，臨床肌電圖檢查多采用針電極插入肌肉檢測(cè)肌電圖，其優(yōu)點(diǎn)是干擾小，定位性好，易識(shí)別，但由于它是一種有創(chuàng)傷的檢測(cè)方法，其應(yīng)用收到了一定的限制。二，表面肌電則是從人體皮膚表面通過電極記錄下來的神經(jīng)肌肉活動(dòng)時(shí)發(fā)放的生物電信號(hào)，屬于無創(chuàng)傷性，操作簡(jiǎn)單,病人易接受，有著廣泛的應(yīng)用前景。主要應(yīng)用領(lǐng)域有：一，仿生學(xué)。提出肌肉

5、生理模型來判別肌肉的動(dòng)作以來，電子假肢的研究進(jìn)入了新的發(fā)展時(shí)期，過去電子假肢的控制靠使用者人為開關(guān)和選擇運(yùn)動(dòng)模式來完成，現(xiàn)在則可通過檢測(cè)人體殘肢表面肌電信號(hào)提取出肢體的動(dòng)作特征，來自動(dòng)控制假肢運(yùn)動(dòng)，利用殘肢表面肌電信號(hào)的肌電假肢研制在國(guó)內(nèi)外都取得較大進(jìn)展。二，康復(fù)工程。如利用表面肌電信號(hào)提取出的特征作為功能性電刺激的控制信號(hào)，幫助癱瘓的肢體恢復(fù)運(yùn)動(dòng)功能。通過檢測(cè)表面肌電信號(hào)，并將其作為反饋信號(hào)提供給病人和醫(yī)生，便于進(jìn)行合理的治療和訓(xùn)練。三，運(yùn)動(dòng)醫(yī)學(xué)。表面肌電信號(hào)在運(yùn)動(dòng)醫(yī)學(xué)中也可發(fā)揮重要作用，通過檢測(cè)運(yùn)動(dòng)員運(yùn)動(dòng)時(shí)的表面肌電信號(hào)，及時(shí)反映出肌肉的疲勞和興奮狀態(tài)，有助于建立科學(xué)的訓(xùn)練方法。本次課程設(shè)

6、計(jì)的主要任務(wù)就是對(duì)微弱的肌電信號(hào)進(jìn)行時(shí)域和頻域的處理及分析，運(yùn)用數(shù)字處理及matlab的知識(shí)進(jìn)行“屈”和“伸”動(dòng)作識(shí)別。然后通過串口將數(shù)據(jù)發(fā)送到單片機(jī)下行微機(jī)進(jìn)行顯示。燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)第二章系統(tǒng)總體設(shè)計(jì)方案根據(jù)課程設(shè)計(jì)要求在上微機(jī)利用matble分析肌電信號(hào)并處理，基于肌電信號(hào)分析結(jié)果，通過串口發(fā)送命令給單片機(jī)系統(tǒng)，根據(jù)肌電信號(hào)動(dòng)作狀態(tài)控制相應(yīng)的數(shù)碼管顯示。并增加了擴(kuò)展模塊，通過動(dòng)作模式驅(qū)動(dòng)電機(jī)轉(zhuǎn)動(dòng)或其他控制輸出模塊。通過分析上位機(jī)matlab中對(duì)信號(hào)處理的結(jié)果，我們可以得到一系列的信號(hào)特征值，其中我們選取了具有代表意義積分肌電值來進(jìn)行處理，并給出對(duì)于屈伸動(dòng)作的閾值。通過對(duì)閾值的判

7、斷，使數(shù)碼管顯示積分肌電值，使點(diǎn)陣模塊顯示相應(yīng)的“屈”和“伸”字樣，使電機(jī)根據(jù)動(dòng)作進(jìn)行正反轉(zhuǎn)，蜂鳴器在“屈”動(dòng)作是發(fā)出鳴響。最后，將上微機(jī)的處理數(shù)據(jù)通過串口通信發(fā)送到下微機(jī)顯示，得到動(dòng)作識(shí)別的要求。燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)第三章肌電信號(hào)的時(shí)域參數(shù)處理及其分析(1)均值：對(duì)于一組隨機(jī)變量來說，均值是一個(gè)很重要的數(shù)值特征，用來描述一組變量的平均水平。具嚴(yán)格的數(shù)學(xué)定義非常簡(jiǎn)單，就是一個(gè)隨機(jī)變量關(guān)于概率測(cè)度的積分。因此，在此處，均值表示肌電信號(hào)的平均水平。公式如下：N(2)方差：方差是各個(gè)數(shù)據(jù)與平均數(shù)之差的平方的平均數(shù)。在數(shù)理統(tǒng)計(jì)和概率論中，方差用來度量隨機(jī)變量和其數(shù)學(xué)期望(即均值)之間的偏離程度

8、。在此處，研究信號(hào)的隨機(jī)變量和均值之間的偏離程度有著很重要的意義。其求解公式如下：X(y產(chǎn)(3)標(biāo)準(zhǔn)差：標(biāo)準(zhǔn)差也稱均方差，是各數(shù)據(jù)偏離平均數(shù)的距離的平均數(shù)，它是離均差平方和平均后的方根，用u表示。標(biāo)準(zhǔn)差是方差的算術(shù)平方根。標(biāo)準(zhǔn)差能反映一個(gè)數(shù)據(jù)集的離散程度。所以能很好的反映肌電信號(hào)的離散程度。求解公式如下：1強(qiáng)=萬£%八不二i(4)積分肌電值IEMG:積分肌電值就是對(duì)所有信號(hào)取絕對(duì)值后盡心均值的求解，由于對(duì)肌電信號(hào)直接求均值，均值近似為零，無法表征信號(hào)間的差異。若對(duì)肌電信號(hào)取絕對(duì)值后再進(jìn)行均值運(yùn)算后，均值恒大于零，因而可用于提取肌電信號(hào)的特征。公式如下：1V-1=討6叫Ni=Q燕山大學(xué)

9、數(shù)字信號(hào)處理課程設(shè)計(jì)(5)均方根RMS均方根就是一組數(shù)據(jù)的平方和除以數(shù)據(jù)的個(gè)數(shù)再開方，均方根是最理想的平方濾波方式的典型，讓濾波更平滑，更大限度的濾掉噪聲。因此，對(duì)肌電信號(hào)求均方根，可以濾除信號(hào)中的噪聲，使濾波后的信號(hào)更平滑、更明顯。公式如下：RMS.士ICTEMG2(f)dt(6)原始信號(hào)的時(shí)域及上述參數(shù)值的Matlab程序allclear;closea=load('qu.txt');/('shen.txt');N=10000;b=0:N-1;axis(0,10000,-1,1);plot(b,a);xlabel('時(shí)間(s)');ylabel

10、('被測(cè)變量y');gridon;fprintf('n數(shù)據(jù)基本信息：n')printf('均值=%7.5fn',mean(a);fprintf('標(biāo)準(zhǔn)差=%7.5fn',sqrt(var(a);fprintf('方差=%7.5fn',var(a);fprintf('積分肌電值IEMG=%7.5fn',mean(abs(a);fprintf('均方根有效值RMS=%7.5fn',sqrt(mean(a.A2);屈信號(hào)數(shù)據(jù)基本信息：均值=0.03502標(biāo)準(zhǔn)差=0.05775方差=0.00

11、334積分肌電值IEMG=0.05437均方根RMS=0.17246燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)屈J原始信方10(10200C300040005Q00600070008000900010000時(shí)間閭JO.-F“酬翩夏樟伸信號(hào)數(shù)據(jù)基本信息：均值=0.00337標(biāo)準(zhǔn)差=0.24421方差=0.05964積分肌電值IEMG=0.12826均方根RMS=0.19993c伸)原始信號(hào)時(shí)間15)燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)第四章肌電信號(hào)的頻域處理方法及其分析3.1 FFT分析：FFT是離散傅立葉變換的快速算法，可以將一個(gè)信號(hào)變換到頻域。有些信號(hào)在時(shí)域上是很難看出什么特征的，但是如果變換到頻域之后，就很容易

12、看出特征了。這就是很多信號(hào)分析采用FFT變換的原因。另外，F(xiàn)FT可以將一個(gè)信號(hào)的頻譜提取出來。采樣得到的數(shù)字信號(hào)FFT變換。N個(gè)采樣點(diǎn)，經(jīng)過FFT之后，就可以得到N個(gè)點(diǎn)的FFT結(jié)果。為了方便進(jìn)行FFT運(yùn)算，通常N取2的整數(shù)次方。假設(shè)采樣頻率為Fs,信號(hào)頻率F,采樣點(diǎn)數(shù)為No那么FFT之后結(jié)果就是一個(gè)為N點(diǎn)的復(fù)數(shù)。每一個(gè)點(diǎn)就對(duì)應(yīng)著一個(gè)頻率點(diǎn)。這個(gè)點(diǎn)的模值，就是該頻率值下的幅度特性。而每個(gè)點(diǎn)的相位就是在該頻率下的信號(hào)的相位。如果要要提高頻率分辨率，就需要增加采樣點(diǎn)數(shù)，也即采樣時(shí)間。頻率分辨率和采樣時(shí)問是倒數(shù)關(guān)系。Matlab實(shí)現(xiàn):clear;closea=load(y=fft(a,1024);fs

13、=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;figure;plot(f,mag);xlabel(ylabel(title(all;'qu.txt'）;/'shen.txt'%1024點(diǎn)傅立葉變換嫻頻譜頻率）；幅值）；'（屈/伸）肌電幅頻'）;figure;plot(f,angle(y);xlabel(ylabel(title(%相頻譜頻率）；相位）；'（屈/伸）肌電相頻'）;gridon;燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)3.2 功率譜分析：功率譜是信號(hào)或噪聲的自相關(guān)函數(shù)的傅里葉變換。如果一確

14、定信號(hào)平均功率為有限的，則該信號(hào)的自相關(guān)函數(shù)存在，如隨機(jī)信號(hào)或隨機(jī)噪聲是由二階隨機(jī)平穩(wěn)函數(shù)表示的，則其自相關(guān)函數(shù)存在。為了更好得描述能量信號(hào)、功率信號(hào)，我們引入能量譜密度和功率譜密度概念。能量譜密度、功率譜密度函數(shù)表示信號(hào)的能量、功率密度隨頻率變化的情況。通過研究功率譜密度，可以幫助了解信號(hào)的功率分布情況，確定信號(hào)的頻帶等。信號(hào)的功率譜即上述FFT分析后，幅頻值mag的平方再除以2得到。功率密度譜雖然描述了隨機(jī)信號(hào)的功率在各個(gè)不同頻率上的分布，但因?yàn)樗鼉H與幅度頻譜有關(guān)，沒有相位信息，所以從已知功率譜還難以完整地恢復(fù)原來的功率信號(hào)。燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)Matlab實(shí)現(xiàn)：a=loadfqu

15、.txt'y=fft(a,1024);fs=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;power1=(mag.A2)/2;plot(f,power1););/('shen.txt')%做1024點(diǎn)傅里葉變換%周期圖法求功率譜xlabel(ylabel(title('頻譜);'功率譜);,(屈/伸)肌電信號(hào)功率譜);500頻譜3000期百2UOQ與1000(屆)肌電信號(hào)功率譜4000(1$)肌電信號(hào)功率譜10地助日5oa頻譜100C10003.3 倒譜分析：倒譜是信號(hào)的傅里葉變換譜經(jīng)對(duì)數(shù)運(yùn)算后再進(jìn)行的傅里葉反變換

16、或者功率譜的對(duì)數(shù)值的逆傅氏變換稱為倒譜。倒頻譜函數(shù)C(q)(powercepstrum)其數(shù)學(xué)表達(dá)式為：C(q)=|IF(log(s(f)-2。其中，s(f)是信號(hào)s(t)的傅里葉變換，log()為取對(duì)數(shù)，IF為逆傅里葉變換。Matlab實(shí)現(xiàn)：前100點(diǎn)的倒譜變化比較明顯，所以取前100點(diǎn)波形圖。a=load('qu.txty=fft(a,1024);fs=1000;N=length(y););/('shen.txt')%ft1024點(diǎn)傅里葉變換燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)%求倒譜mag=abs(y);f=(0:N-1)/N*fs;w=rceps(a);n=1:100;

17、plot(n,w(1:100);xlabel('時(shí)間);ylabel('倒譜);title('(屈/伸)肌電信號(hào)倒譜);gridon;（屈）肌電信罵倒譜y=fft(a,1024);fs=1000;）;%10000點(diǎn)傅里葉變換103.4平均功率頻率MPF和中值頻率MF:平均功率頻率是總功率除以總時(shí)間中值頻率是各個(gè)時(shí)間段的功率的平均值。公式如下：Matlab實(shí)現(xiàn)：a=load('qu.txt燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)N=length(y);mag=abs(y);f=(0:N-1)/N*fs;power=(mag.A2)/2;ss=sum(power);M2=0.5

18、*ss;df=fs/N;M1=0.5*df*(sum(power(1:N-1)+sum(power(2:N);MPF=M1/M2;MF=M2/2;第五章Matlab程序及GUIMATLA是一種面向工程和科學(xué)計(jì)算的交互式計(jì)算軟件，它以矩陣運(yùn)算為基礎(chǔ)，把計(jì)算、可視化、程序設(shè)計(jì)融合到了一個(gè)簡(jiǎn)單易用的交互式工作環(huán)境中。同時(shí)由于MATLA是一個(gè)數(shù)據(jù)分析和處理功能十分強(qiáng)大的工程實(shí)用軟件，它的信號(hào)處理與分析工具箱為肌電信號(hào)分析提供了十分豐富的功能函數(shù)，利用這些功能函數(shù)可以快捷而又方便地完成肌電信號(hào)的處理和分析，再通過GUI界面設(shè)計(jì)，顯示信號(hào)處理后的結(jié)果，方便簡(jiǎn)明。打開GUI設(shè)計(jì)界面，根據(jù)要求添加相應(yīng)的按鈕，

19、并對(duì)按鈕進(jìn)行定義，調(diào)整界面整體布局，將界面運(yùn)行，會(huì)產(chǎn)生相應(yīng)的程序。在程序中，尋找相應(yīng)的按鈕定義，在定義下添加該按鈕要實(shí)現(xiàn)的功能函數(shù)。然后對(duì)整個(gè)程序進(jìn)行調(diào)試驗(yàn)證。得到可視化GUI界面及功能的實(shí)現(xiàn)。11燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)functionvarargout=untitled20(varargin)%UNTITLED20M-fileforuntitled20.fig%UNTITLED20,byitself,createsanewUNTITLED20orraisestheexisting%singleton*.%H=UNTITLED20returnsthehandletoanewUNTITLE

20、D20orthehandleto%theexistingsingleton*.%UNTITLED20('CALLBACK',hObject,eventData,handles,.)callsthelocal%functionnamedCALLBACKinUNTITLED20.Mwiththegiveninputarguments.%UNTITLED20('Property','Value',.)createsanewUNTITLED20orraisesthe%existingsingleton*.Startingfromtheleft,prope

21、rtyvaluepairsare%appliedtotheGUIbeforeuntitled20_OpeningFunctiongetscalled.An%unrecognizedpropertynameorinvalidvaluemakespropertyapplication%stop.Allinputsarepassedtountitled20_OpeningFcnviavarargin.%*SeeGUIOptionsonGUIDE'sToolsmenu.Choose"GUIallowsonlyone%instancetorun(singleton)".%Se

22、ealso:GUIDE,GUIDATA,GUIHANDLES%Edittheabovetexttomodifytheresponsetohelpuntitled20%LastModifiedbyGUIDEv2.501-Jul-201316:55:28%Begininitializationcode-DONOTEDIT12gui_Singleton=1;gui_State=struct('gui_Name''gui_Singleton''gui_OpeningFcn''gui_OutputFcn''gui_LayoutFcn

23、''gui_Callback'燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)mfilename,.,gui_Singleton,.,untitled20_OpeningFcn,untitled20_OutputFcn,，口,.,);ifnargin&&ischar(varargin1)gui_State.gui_Callback=str2func(varargin1);endvarargin:);ifnargoutvarargout1:nargout=gui_mainfcn(gui_State,elsegui_mainfcn(gui_State,varargin:);en

24、d%Endinitializationcode-DONOTEDIT%-Executesjustbeforeuntitled20ismadevisible.functionuntitled20_OpeningFcn(hObject,eventdata,handles,varargin)%Thisfunctionhasnooutputargs,seeOutputFcn.%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(

25、seeGUIDATA)%varargincommandlineargumentstountitled20(seeVARARGIN)%Choosedefaultcommandlineoutputforuntitled20handles.output=hObject;%Updatehandlesstructureguidata(hObject,handles);%UIWAITmakesuntitled20waitforuserresponse(seeUIRESUME)%uiwait(handles.figure1);%-Outputsfromthisfunctionarereturnedtothe

26、commandline.13燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)functionvarargout=untitled20_OutputFcn(hObject,handles)%varargoutcellarrayforreturningoutputargs(see%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversion%handlesstructurewithhandlesanduserdata(seeGUIDATA)%Getdefaultcommandlineoutputfromhandlesstructurevarar

27、gout1=handles.output;%-Executesonbuttonpressinqu.functionqu_Callback(hObject,eventdata,handles)%hObjecthandletoqu(seeGCBO)%eventdatareserved-tobedefinedinafutureversion%handlesstructurewithhandlesanduserdata(seeGUIDATA)a=load('qu.txt');axes(handles.axes1);plot(a);xlabel('時(shí)間(s)');ylab

28、el('被測(cè)變量y');title('(屈)原始信號(hào)');gridon;A=mean(a);%平均值B=std(a);%標(biāo)準(zhǔn)差C=Var(a);%方差D=mean(abs(a);%積分機(jī)電值E=sqrt(mean(a.A2);%勻方根有效值RMSy=fft(a,10000);%10000點(diǎn)傅里葉變換fs=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;power=(mag.A2)/10000;%率譜ss=sum(power);M2=0.5*ss;df=fs/N;eventdata,VARARGOUT);ofMATLABofM

29、ATLABM1=0.5*df*(sum(power(1:N-1)+sum(power(2:N);14燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)F=M1/M2;G=M2/2;%A,B,C,D,E,F=index_calculation(a);%設(shè)置顯示A1=A*100000;A2=round(A1);%-IfEnable%-Otherwise,='on',executesonmousepressin5pixelexecutesonmousepressin5pixelborderborder.oroverA3=A2/100000;%round函數(shù)是對(duì)數(shù)據(jù)取整，以此限制數(shù)據(jù)的有效位set(hand

30、les.pingjunzhi,'string',num2str(A3);B1=B*100000;B2=round(B1);B3=B2/100000;set(handles.biaozhuncha,'string',num2str(B3);C1=C*100000;C2=round(C1);C3=C2/100000;set(handles.fangcha,'string',num2str(C3);D1=D*100000;D2=round(D1);D3=D2/100000;set(handles.jifenjidianzhi,'string&#

31、39;,num2str(D3);E1=E*100000;E2=round(E1);E3=E2/100000;set(handles.junfanggen,'string',num2str(E3);F1=F*100000;F2=round(F1);F3=F2/100000;set(handles.pingjungonglv,'string',num2str(F3);G1=G*100000;G2=round(G1);G3=G2/100000;set(handles.zhongweipinglv,'string',num2str(G3);pingjun

32、zhi.15燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)functionpingjunzhi_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletopingjunzhi(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%-IfEnable='on',executesonmousepressin5pixelborder.%-Otherwise,executeson

33、mousepressin5pixelborderoroverbiaozhuncha.functionbiaozhuncha_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletobiaozhuncha(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%-IfEnable='on',executesonmousepressin5pixelborde

34、r.%-Otherwise,executesonmousepressin5pixelborderoroverjifenjidianzhi.functionjifenjidianzhi_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletojifenjidianzhi(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%-IfEnable='on',

35、executesonmousepressin5pixelborder.%-Otherwise,executesonmousepressin5pixelborderoroverpingjungonglv.functionpingjungonglv_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletopingjungonglv(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%-IfEnable%-Otherwise,fangcha.%handlesstructu

36、rewithhandlesanduserdata(seeGUIDATA)='on',executesonmousepressin5pixelborder.executesonmousepressin5pixelborderorover16燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)functionfangcha_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletofangcha(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructure

37、withhandlesanduserdata(seeGUIDATA)%-Executesonbuttonpressinshen.functionshen_Callback(hObject,eventdata,handles)%hObjecthandletoshen(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)a=load('shen.txt');axes(handles.axes2);plot(a)

38、;xlabel('時(shí)間(s)');ylabel('被測(cè)變量y');title('(伸)原始信號(hào));gridon;A=mean(a);%平均值B=std(a);%標(biāo)準(zhǔn)差C=Var(a);%方差D=mean(abs(a);%積分機(jī)電值E=sqrt(mean(a.A2);%勻方根有效值RMSy=fft(a,10000);%10000點(diǎn)傅里葉變換fs=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;power=(mag.A2)/10000;%周期圖法功率譜ss=sum(power);M2=0.5*ss;df=fs/N;M1=

39、0.5*df*(sum(power(1:N-1)+sum(power(2:N);F=M1/M2;G=M2/2;%A,B,C,D,E,F=index_calculation(a);%設(shè)置顯示A1=A*100000;17燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)A2=round(A1);A3=A2/100000;set(handles.pingjunzhi,'string',num2str(A3);B1=B*100000;B2=round(B1);B3=B2/100000;set(handles.biaozhuncha,'string',num2str(B3);C1=C*1000

40、00;C2=round(C1);C3=C2/100000;set(handles.fangcha,'string',num2str(C3);D1=D*100000;D2=round(D1);D3=D2/100000;set(handles.jifenjidianzhi,'string',num2str(D3);E1=E*100000;E2=round(E1);E3=E2/100000;set(handles.junfanggen,'string',num2str(E3);F1=F*100000;F2=round(F1);F3=F2/100000;

41、set(handles.pingjungonglv,'string',num2str(F3);G1=G*100000;G2=round(G1);G3=G2/100000;set(handles.zhongweipinglv,'string',num2str(G3);%-Executesonbuttonpressinqu1.functionqu1_Callback(hObject,eventdata,handles)MATLAB%hObjecthandletoqu1(seeGCBO)%eventdatareserved-tobedefinedinafutureve

42、rsionof%handlesstructurewithhandlesanduserdata(seeGUIDATA)a=load('qu.txt');y=fft(a,1024);%故1024點(diǎn)傅里葉變換fs=1000;N=length(y);mag=abs(y);18燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)f=(0:N-1)/N*fs;axes(handles.axes3);plot(f(1:512),mag(1:512);%故幅頻譜xlabel('頻率');ylabel('幅值');title('(屈)肌電幅頻譜)；axes(handles.axe

43、s4);plot(f,angle(y);%做相頻譜xlabel('頻率');ylabel('幅值');title('(屈)肌電相頻譜)；gridon;w=rceps(a);%t倒譜n=1:100;axes(handles.axes5);plot(n,w(1:100);xlabel('時(shí)間);ylabel('倒譜');title('(屈)肌電信號(hào)倒譜)；gridon;power1=(mag.A2)/2;%期圖法求功率譜axes(handles.axes6);plot(f,power1);xlabel('頻譜'

44、);ylabel('功率譜');title('(屈)肌電信號(hào)功率譜)；gridon;%-Executesonbuttonpressinshen1.functionshen1_Callback(hObject,eventdata,handles)%hObjecthandletoshen1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)a=load('shen.txt');y=ff

45、t(a,1024);%故1024點(diǎn)傅里葉變換fs=1000;N=length(y);mag=abs(y);f=(0:N-1)/N*fs;19燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)axes(handles.axes3);plot(f(1:512),mag(1:512);%故幅頻譜xlabel('頻率');ylabel('幅值');title('(伸)肌電幅頻譜)；axes(handles.axes4);plot(f,angle(y);%做相頻譜xlabel('頻率');ylabel('幅值');title('(伸)肌電相頻譜)

46、；gridonw=rceps(a);%t倒譜n=1:100;axes(handles.axes5);plot(n,w(1:100);xlabel('時(shí)間);ylabel('倒譜');title('(伸)肌電信號(hào)倒譜)；gridon;power1=(mag.A2)/10000;%周期圖法求功率譜axes(handles.axes6);plot(f,power1);xlabel('頻譜');ylabel('功率譜');title('(伸)肌電信號(hào)功率譜)；gridon;%-Executesonmousepressoveraxe

47、sbackground.functionaxes5_ButtonDownFcn(hObject,eventdata,handles)%hObjecthandletoaxes5(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)20燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)%-Executesonbuttonpressinqufa.functionqufa_Callback(hObject,eventdata,handles)%hObj

48、ecthandletoqufa(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalss%a=load('qu.txt');%D=mean(abs(a);%積分機(jī)電值%D1=D*100000;%D2=round(D1);%D3=dec2hex(D2);dtmf_series(01);dtmf_series(01);%-Executesonbuttonpressinshenfa.functions

49、henfa_Callback(hObject,eventdata,handles)%hObjecthandletoshenfa(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)globalss%a=load('shen.txt');%D=mean(abs(a);%積分機(jī)電值%D1=D*100000;%D2=round(D1);%D3=dec2hex(D2);dtmf_series(00);dtmf_series(00);functionUntitled_2_Callback(hObject,eventdata,handles)ofMATLAB%hObjecthandletoqu1(seeGCBO)%eventdatareserved-tobedefinedinafutureversion21燕山大學(xué)數(shù)字信號(hào)處理課程設(shè)計(jì)%handlesstructurewith