肌电信号的时域和频域分析

表面肌电信号一般只有微伏级电压,信号中往往夹带着低频(接近直流)和高频的干扰信号,真正有用的肌电信号大致在10Hz~500Hz之间。由上面的对比图可知，消除了噪声信号，有用信号得以保存。

其部分程序代码如下所示：

global a; global t; %global h;

[C L]=wavedec(a,3,'db5'); cA3=appcoef(C,L,'db5',3); cD1=detcoef(C,L,1);

13

cD2=detcoef(C,L,2); cD3=detcoef(C,L,3);

thr1=thselect(cD1,'rigrsure'); thr2=thselect(cD2,'rigrsure'); thr3=thselect(cD3,'rigrsure'); TR=[thr1,thr2,thr3]; SORH='s';

[XC,CXC,LXC,PERFO,PERF2]=wdencmp('lvd',a,... 'db5',3,TR,SORH); y1=fft(XC,20000); fs=2000; N=length(y1); mag1=abs(y1); f=(0:N-1)/N*fs;

plot(handles.axes5,f,mag1);

去噪前后的功率谱图如下所示：

14

部分程序代码如下： global a; global t; global s; s=fft(a,2000);

[C L]=wavedec(a,3,'db5'); cA3=appcoef(C,L,'db5',3); cD1=detcoef(C,L,1); cD2=detcoef(C,L,2); cD3=detcoef(C,L,3);

thr1=thselect(cD1,'rigrsure');

15

thr2=thselect(cD2,'rigrsure'); thr3=thselect(cD3,'rigrsure'); TR=[thr1,thr2,thr3]; SORH='s';

[XC,CXC,LXC,PERFO,PERF2]=wdencmp('lvd',a,... 'db5',3,TR,SORH); y1=fft(XC,2000); fs=2000; N=length(y1);

mag=abs(y1); f=(0:N-1)/N*fs; power2=(mag.^2)/2000; plot(handles.axes5,power2);

16