急！！！在线等，求解答：一个小波去噪的matlab程序，高手进

微信公开课小程序热门讨论「急！！！在线等，求解答：一个小波去噪的matlab程序，高手进」最新回复-%%%%%%%%%%%%%%%%%%心电信号降噪%%%%%%%%%%%%%%%BirgeMassart策略阈值降噪%基于小波变换的心电信号的降噪ecg=fopen(100.dat,r);%调用心电数据库r为只读，ecg是打开文件的识别符N=1201;%常数赋值，要读数据个数data=fread(ecg,N,int16);%从一个流中读N个数据，数据格式是int16,16进制整数data=data/10000;%数据缩小10000倍fclose(ecg);%关闭打开的文件x=data;%把数据转赋给x变量wavename=db5;%db5是小波名level=4;%4级分解[c,l]=wavedec(x,level,wavename);%4级小波分解，c保存各级分解系数，l是薄记矩阵，保存各级的系数的个数alpha=1.5;%1.5用于信号压缩，3用于降噪sorh=h;%为硬阈值[thr,nkeep]=wdcbm(c,l,alpha);%使用BirgéMassart策略计算一维小波分解或压缩的阈值thr和各级的系数个数nkeep[xc,cxc,lxc,perf0,perfl2]=wdencmp(lvd,c,l,wavename,level,thr,sorh);%小波压缩重构后的图像t1=0:0.004:(length(x)1)*0.004;%一行数据figure(4);%打开一个图形窗口subplot(211);%子图1plot(t1,x);%画图形title(从人体采集的原始的ECG信号);%加上子图名称subplot(212);%子图2plot(t1,xc);%画图形title(BirgeMassart策略阈值降噪后的ECG信号（wname=db5level=4）);%加上子图名称
供参考：lev=5;[c,l]=wavedec(x,lev,wname);sigma=wnoisest(c,l,1);alpha=2;thr1=wbmpen(c,l,sigma,alpha)[thr2,nkeep]=wdcbm(c,l,alpha)xd1=wdencmp(gbl,c,l,wname,lev,thr1,s,1);[xd2,cxd,lxd,perf0,perfl2]=wdencmp(lvd,c,l,wname,lev,thr2,h);[thr,sorh,keepapp]=ddencmp(den,wv,x)xd3=wdencmp(gbl,c,l,wname,lev,thr,s,1);subplot(411);plot(x);title(原始信号,fontsize,12);subplot(412);plot(xd1);title(使用penalty阈值降噪后信号,fontsize,12);subplot(413);plot(xd2);title(使用Birge-Massart阈值降噪后信号,fontsize,12);subplot(414);plot(xd3);title(使用缺省阈值降噪后信号,fontsize,12);s=[-1.580.420.460.78-0.490.59-1.3-1.42-0.16-1.47-1.350.36-0.44-0.141-0.5-0.2-0.06-0.60.42-1.520.510.76-1.50.16-1.29-0.65-1.480.6-1.65-0.55];[C,L]=wavedec(s,1,db3);ca1=wrcoef(a,C,L,db3,1);x1=ca1;[C,L]=wavedec(s,2,db3);ca2=wrcoef(a,C,L,db3,2);x2=ca2;[C,L]=wavedec(s,3,db3);ca3=wrcoef(a,C,L,db3,3);x3=ca3;[C,L]=wavedec(s,4,db3);ca4=wrcoef(a,C,L,db3,4);x4=ca4;cg=wrcoef(a,C,L,sym5,1);x5=cg;p=1:31;subplot(6,1,1);plot(p,s);ylabel(s);subplot(6,1,2);plot(p,x1);ylabel(ca1);subplot(6,1,3);plot(p,x2);ylabel(ca2);subplot(6,1,4);plot(p,x3);ylabel(ca3);subplot(6,1,5);plot(p,x4);ylabel(ca4)subplot(6,1,6);plot(p,x5);ylabel(ca5)%加入的重构，是不是你要的？ 更多有关「急！！！在线等，求解答：一个小波去噪的matlab程序，高手进」的疑问请扫码关注微信公开课+小程序！