《生医讯号处理概论.ppt》由会员分享，可在线阅读，更多相关《生医讯号处理概论.ppt（58页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、生医讯号处理概论 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望大綱Electromyography(EMG)簡介Page.03Electromyography(EMG)測量的方法Page.05Electromyography(EMG)處理的目的Page.Electromyography(EMG)處理的流程Page.Electromyography(EMG)處理的方法Page.Electromyography(EMG)處理的結果Page.報告結論Page.心得Pa
2、ge.Electromyography(EMG)簡介肌電圖檢查是利用神經及肌肉的電生理特性，以電流刺激神經，記錄其運動和感覺的反應，肌電圖包括了神經傳導檢查和針極肌電圖來作為幫助診斷神經肌肉疾病的一種檢查。它可以偵測神經病變或損傷的位置、性質與嚴重程度，是神經內科重要的檢查工具，以便協助醫師能夠正確做出臨床診斷、選擇治療方式、評估治療效果和推斷疾病的預後。針極肌電圖是將針極刺入肌肉，來記錄在各種電擊狀態下的電位活動，再經由多條肌肉的檢查來判定神經、肌肉病變的特性，部位及範圍和嚴重度。一般常用針極大致分為同軸針極及單極針極兩種，同軸針極的部分在材質上是比較耐用，在訊號干擾部分也比較少，相對的對患
3、者來說會比較有疼痛的感覺；在單極針的部份好處就是可記錄的面積大，對患者來講也比較不會有疼痛的感覺，但單極針也比較容易損壞而且干擾也大，也各有優缺點。傳統肌電圖&表面肌電圖傳統肌電圖是應用肌電儀記錄肌肉在放鬆和收縮時的生物電活動同時監聽聲音變化，結合神經傳導速度測定，可確定神經、肌肉功能狀態的檢查方法。通過病理狀態下，發生相應的肌電圖變化，協助診斷和鑒別診斷。表面肌電圖（surface electromyography，sEMG），又稱動態肌電圖（dynamic electromyography,DEMG），是從肌肉表面通過電極引導、記錄下來的神經肌肉系統活動時的生物電信號。它與肌肉的活動狀態和
4、功能狀態之間存在著不同程度的關聯性，因而能在一定程度上反映神經肌肉的活動。肌肉運動中產生的生物電通過兩個測量電極（相對於參考電極）產生電位差，差分放大器檢測到該信號後，經過放大、記錄後所得到的圖形，現代高檔的sEMG都是把放大的信號再轉化為數字信號，經過通訊系統傳輸給微機。微機中的分析軟件對所獲得的數據進行分析處理，從而完成測試評估等科研或臨床診斷任務。sEMG是一種簡單、無創、大範圍測量且容易被受試者接受的肌電活動。Electromyography(EMG)測量的方法表面導出法：把但及貼附在皮膚上導出電位的方法，可做神經傳導檢查和特出檢查。主要展示神經和肌肉接頭功能（需電刺激）。針電極法：把
5、針電極刺入肌肉測得某個運動單位中的肌肉放鬆時的自發電位和肌肉主動收縮時的運動單位電位變化，主要展示肌肉本身的功能狀態（不需任何電刺激）。表面導出法&針電極法表面導出法針電極法Electromyography(EMG)處理的目的偵測神經病變或損傷的位置、性質與嚴重程度，是神經內科重要的檢查工具，以便協助醫師能夠正確做出臨床診斷、選擇治療方式、評估治療效果和推斷疾病的預後。針極肌電圖是將針極刺入肌肉，來記錄在各種電擊狀態下的電位活動，再經由多條肌肉的檢查來判定神經、肌肉病變的特性，部位及範圍和嚴重度。可用於肌病、周圍神經病，神經終板疾病的診斷，可協助判別病變的部位、程度、預後。可以鑒別肌原性或神經
6、原性損害；區別周圍神經的髓鞘損害和軸索損害；確定神經損傷及壓迫的部位、程度及預後；判定神經吻合後功能恢復的情況；判定咀嚼肌的功能，膀胱括約肌、肛門括約肌的功能。肌電訊號量測與分析方法：肌電圖(electromyography,EMG)是肌肉或肌群生理條件的量測方法之一，這是因為神經傳導是靠電位的變化，將電極置於肌肉表面皮膚或肌肉內部，量測肌肉活動時，微弱的電流變化，並將其訊號放大且紀錄。表面電極(Surface Electrode)：在肌肉外緣的皮膚面直接貼上電極，以非侵入的方式(Non-Invasive)量測，將電極固定在預偵測的肌肉皮膚上即可，既方便又安全。雖然能記錄表層肌肉的活動，卻無法
7、選擇性的測量小肌肉群的運動訊號，故較適合用在測量表層大區域肌肉的訊號。電極貼片Electromyography(EMG)原圖Electromyography(EMG)頻譜分析EMG訊號與施力的關係：一般以肌電活動性的方式來評估施力的強度，而強度分析的指標有均方根值(root mean square,RMS)，平均振幅(MA)，皆可反應活動中，肌肉的施力程度，在不疲勞的狀況下，同一條肌肉肌電訊號越大則表示施力越大，在靜態等長收縮施力中更為明顯。另外以MVC(maximum voluntary contraction)最大肌力收縮量做為變數來評估肌電訊號，MVC指的是人在自由意志之下，肌肉收縮所得
8、到的最大收縮力量，與個人肌肉強度及疲勞狀態息息相關EMG訊號與肌肉疲勞的關係：肌肉疲勞與否和程度，訊號頻譜分析為一重要指標。頻譜分佈的變化與施力大小、速度大小以及疲勞度有關。在為疲勞之前，施力的變化會使頻譜分佈在較高之頻率下，疲勞之後，會使頻譜分佈在較低之頻帶，原因是當肌肉疲勞時，肌纖維間的傳導速度減慢。肌肉疲勞與頻譜分析圖矩形窗海明窗漢寧窗三角窗黑人窗FFT原圖FFT detending(1)FFT detending(2)矩形窗FFT%加窗FFT%矩形窗FFTfd=1/length(signalEMG)*fs;f=0:fd:(length(signalEMG)-1)*fd;rectwin_
9、wind_1_fft=fft(rectwin_wind_1)figure(15)subplot(211)plot(f,abs(rectwin_wind_1_fft)title(矩形窗FFT傅利葉頻譜)subplot(212)plot(f,angle(rectwin_wind_1_fft)title(矩形窗FFT傅利葉相位)海明窗FFT%海明窗FFTfd=1/length(signalEMG)*fs;f=0:fd:(length(signalEMG)-1)*fd;hamming_wind_1_fft=fft(hamming_wind_1)figure(16)subplot(211)plot(f,
10、abs(hamming_wind_1_fft)title(海明窗FFT傅利葉頻譜)subplot(212)plot(f,angle(hamming_wind_1_fft)title(海明窗FFT傅利葉相位)漢寧窗FFT%漢寧窗FFTfd=1/length(signalEMG)*fs;f=0:fd:(length(signalEMG)-1)*fd;hanning_wind_1_fft=fft(hanning_wind_1)figure(17)subplot(211)plot(f,abs(hanning_wind_1_fft)title(漢寧窗FFT傅利葉頻譜)subplot(212)plot(f
11、,angle(hanning_wind_1_fft)title(漢寧窗FFT傅利葉相位)三角窗FFT%三角窗FFTfd=1/length(signalEMG)*fs;f=0:fd:(length(signalEMG)-1)*fd;bartlett_wind_1_fft=fft(bartlett_wind_1)figure(18)subplot(211)plot(f,abs(bartlett_wind_1_fft)title(三角窗FFT傅利葉頻譜)subplot(212)plot(f,angle(bartlett_wind_1_fft)title(三角窗FFT傅利葉相位)布雷克曼窗FFT%布雷
12、克曼窗FFT fd=1/length(signalEMG)*fs;f=0:fd:(length(signalEMG)-1)*fd;blackman_wind_1_fft=fft(blackman_wind_1)figure(19)subplot(211)plot(f,abs(blackman_wind_1_fft)title(布雷克曼窗FFT傅利葉頻譜)subplot(212)plot(f,angle(blackman_wind_1_fft)title(布雷克曼窗FFT傅利葉相位)M點濾波%M點濾波M1=30;A1=ones(M1,1)/M1;filter_A1=filter(A1,1,EMG
13、3);M2=100;A2=ones(M2,1)/M2;filter_A2=filter(A2,1,EMG3);figure(20)subplot(211)plot(tc,EMG3,r,tc,filter_A1,b)title(M點平均濾波(M=30)xlabel(Frequency(Hz)ylabel(Magnitude)subplot(212)plot(tc,EMG3,r,tc,filter_A2,b)title(M點平均濾波(M=100)xlabel(Frequency(Hz)ylabel(Magnitude)M=30,100點比較M=30加窗FFTM=100 加窗FFTChebyshev
14、 I 程式碼Chebyshev I 圖片Chebyshev II 程式碼Chebyshev II 圖Butterworth 程式碼Butterworth 圖Elliptic 程式碼Elliptic 圖波德圖程式碼波德圖M點移動平均濾波器M=10 M=100同步平均同步平均PSD功率頻譜密度PSD功率頻譜密度圖週期圖+修正週期圖+修正週期圖週期圖取1024點週期圖窗型比較Pwelch取不同點75%50%Pwelch取不同點75%50%圖各窗比較50%各窗比較圖50%Yule walker AR階數比較10298Yule walker AR階數比較10298結論-訊號處理流程原始訊號去基線飄移高通
15、低通帶拒FFT強度相位頻譜FIR移動平均濾波器PSD功率頻譜密度非參數方法periodgram,welch參數方法Yule-walker AR模型結論使用布雷克曼窗,使得旁辦變少波形變得很漂亮以及使用M點移動平均濾波器,雖會將能量變小,但能夠抓出肌電圖波形大概樣子同步平均濾波器:用於週期性訊號較佳,以至於肌電圖不大適用心得這學期跟上學期比起來不一樣的是，上學期學的都是些基本的matlab使用方法，而這學期，我們是實際應用在這些生理訊號上面，其實不是說很難，只是上課一定要用心聽才會，而我們這組做的肌電圖，訊號量測出來後，其實不是很漂亮，可是經過我們長時間分析及濾波之後，終於也找到了比較適合肌電圖的濾波器。