生医讯号处理概论.ppt

生醫訊號處理概論,報告主題：Electromyography(EMG-肌電圖)組長：姜禮翊組員：吳庭毅王振彥,大綱,Electromyography(EMG)簡介Page.03Electromyography(EMG)測量的方法Page.05Electromyography(EMG)處理的目的Page.Electromyography(EMG)處理的流程Page.Electromyography(EMG)處理的方法Page.Electromyography(EMG)處理的結果Page.報告結論Page.心得Page.,Electromyography(EMG)簡介,肌電圖檢查是利用神經及肌肉的電生理特性，以電流刺激神經，記錄其運動和感覺的反應，肌電圖包括了神經傳導檢查和針極肌電圖來作為幫助診斷神經肌肉疾病的一種檢查。它可以偵測神經病變或損傷的位置、性質與嚴重程度，是神經內科重要的檢查工具，以便協助醫師能夠正確做出臨床診斷、選擇治療方式、評估治療效果和推斷疾病的預後。針極肌電圖是將針極刺入肌肉，來記錄在各種電擊狀態下的電位活動，再經由多條肌肉的檢查來判定神經、肌肉病變的特性，部位及範圍和嚴重度。一般常用針極大致分為同軸針極及單極針極兩種，同軸針極的部分在材質上是比較耐用，在訊號干擾部分也比較少，相對的對患者來說會比較有疼痛的感覺；在單極針的部份好處就是可記錄的面積大，對患者來講也比較不會有疼痛的感覺，但單極針也比較容易損壞而且干擾也大，也各有優缺點。,傳統肌電圖f=0:fd:(length(signalEMG)-1)*fd;rectwin_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,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,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,%布雷克曼窗FFTfd=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,EMG3);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加窗FFT,M=100加窗FFT,ChebyshevI程式碼,ChebyshevI圖片,ChebyshevII程式碼,ChebyshevII圖,Butterworth程式碼,Butterworth圖,Elliptic程式碼,Elliptic圖,波德圖程式碼,波德圖,M點移動平均濾波器,M=10M=100,同步平均,同步平均,PSD功率頻譜密度,PSD功率頻譜密度圖,週期圖+修正,週期圖+修正,週期圖,週期圖取1024點,週期圖窗型比較,Pwelch取不同點75%50%,Pwelch取不同點75%50%圖,各窗比較50%,各窗比較圖50%,YulewalkerAR階數比較10298,YulewalkerAR階數比較10298,結論-訊號處理流程,原始訊號,去基線飄移,高通,低通,帶拒,FFT,強度相位頻譜,FIR移動平均濾波器,PSD功率頻譜密度,非參數方法periodgram,welch,參數方法Yule-walkerAR模型,結論,使用布雷克曼窗,使得旁辦變少波形變得很漂亮以及使用M點移動平均濾波器,雖會將能量變小,但能夠抓出肌電圖波形大概樣子同步平均濾波器:用於週期性訊號較佳,以至於肌電圖不大適用,心得,這學期