第一范文网 - 专业文章范例文档资料分享平台

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

来源:用户分享 时间:2025/5/20 3:57:07 本文由loading 分享 下载这篇文档手机版
说明:文章内容仅供预览,部分内容可能不全,需要完整文档或者需要复制内容,请下载word后使用。下载word有问题请添加微信号:xxxxxxx或QQ:xxxxxx 处理(尽可能给您提供完整文档),感谢您的支持与谅解。

通过比较上图可以看出,第四幅图像的去噪效果比较明显,因次在本次设计中选择小波四层分解对信号进行去噪处理。

分别对20N 的肌电信号和50N的肌电信号进行小波去噪,去噪前后的效果图如下所示:

9

左边是20N小波去噪前后的图,右边是50N作用前后的效果图。通过前后图像的比较可以看出,当选用四层小波分解进行滤波时,已经能去掉大部分噪声的干扰,保存有用的信号。

小波去噪部分程序代码如下:

global a; global t; global s; global y; % global IR; % global SORH; M =length (a) ;

10

N=length(y); p=size(a); s=a(1:20000);

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

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

[XC,CXC,LXC,PERFO,PERF2]=wdencmp('lvd',a,... 'db5',4,TR,SORH); L=p(2); x=a; h=XC; F=0; M=0; for ii=1:L

m(ii)=(x(ii)-y(ii))^2; t(ii)=y(ii)^2; f(ii)=t(ii)/m(ii); F=F+f(ii); M=M+m(ii); end;

SNR=10*log10(F);

11

MSE=M/N; SM=SNR/MSE; % K=length(d); % t1=(0:K-1)/2000;

plot(handles.axes5,XC(1:20000));

小波去噪前后的幅频和相频图如下:

12

搜索更多关于: 肌电信号的时域和频域分析 的文档
肌电信号的时域和频域分析.doc 将本文的Word文档下载到电脑,方便复制、编辑、收藏和打印
本文链接:https://www.diyifanwen.net/c20qd0282rm0a6rj16zo8_3.html(转载请注明文章来源)
热门推荐
Copyright © 2012-2023 第一范文网 版权所有 免责声明 | 联系我们
声明 :本网站尊重并保护知识产权,根据《信息网络传播权保护条例》,如果我们转载的作品侵犯了您的权利,请在一个月内通知我们,我们会及时删除。
客服QQ:xxxxxx 邮箱:xxxxxx@qq.com
渝ICP备2023013149号
Top