发布时间 : 星期五 文章肌电信号的时域和频域分析更新完毕开始阅读33c7743902020740bf1e9b18
表面肌电信号一般只有微伏级电压,信号中往往夹带着低频(接近直流)和高频的干扰信号,真正有用的肌电信号大致在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