软硬和改进阈值处理.docx

我目前接触了信号去噪方面, 可是在去噪时候试了网上的改进折中阈值函数, 发现所有的改进函数去噪信噪比都不如硬阈值,我这是那出错了,求好心人帮忙解答!!! (附) 程序 load data x=data; % 产生含噪信号 init=20556115866; randn('seed',init); nx=x+*randn(size(x)); s=nx'; %用 db1 小波对原始信号进行 3 层分解并提取系数[c,l]=wavedec(s,3,'db1'); ca3=appcoef(c,l,'db1',3); cd3=detcoef(c,l,3); cd2=detcoef(c,l,2); cd1=detcoef(c,l,1); sigma3=median(abs(cd3))/; sigma2=median(abs(cd2))/; sigma1=median(abs(cd1))/; thr3=sigma3*sqrt(2*log(2000)); thr2=sigma2*sqrt(2*log(2000)); thr1=sigma1*sqrt(2*log(2000)); % d1=wthresh(cd1,'h',thr1); ccd2=wthresh(cd2,'h',thr2); ccd3=wthresh(cd3,'h',thr3); >> c1=[d1]; s1=waverec(c1,l,'db1'); subplot(411);plot(s); title(' 原始信号') subplot(412);plot(s1); title(' 硬阈值去噪后信号') grid on >> % 进行软阈值处理>> % 进行软阈值处理 scd1=wthresh(cd1,'s',thr1); scd2=wthresh(cd2,'s',thr2); scd3=wthresh(cd3,'s',thr3); c2=[ca3 scd3 scd2 scd1]; s2=waverec(c2,l,'db