%使用wnoise函數產生無噪信號和加噪信號并顯示
snr=2.5;
init=5055615866;
[sref,s]=wnoise(2,10,snr,init);
sref=sref(1:1000);
s=s(1:1000);
subplot(2,2,1);
plot(s);
title('加噪信號');
%用db8小波對原始信號進行3層分解，并提取系數
wmane='db8';
[c,l]=wavedec(s,3,wmane);
a3=appcoef(c,l,wmane,3);
d3=detcoef(c,l,3);
d2=detcoef(c,l,2);
d1=detcoef(c,l,1);
%對信號進行強制性消噪處理并顯示結果
dd3=zeros(1,length(d3));
dd2=zeros(1,length(d2));
dd1=zeros(1,length(d1));
c1=[a3,dd3,dd2,dd1];
s1=waverec(c1,l,wmane);
subplot(2,2,2);
plot(s1);
title('強制消噪后的信號')
%用默認閾值對信號進行消噪并顯示結果。
[thr,sorth,keepapp]=ddencmp('den','wv',s);
s2=wdencmp('gbl',c,l,wmane,3,thr,sorth,keepapp);
subplot(2,2,3);
plot(s2);
title('默認閾值消噪后的信號');
%用給定的閾值進行消噪并顯示結果
alpha=2;
[thr2,nkeep]=wdcbm(c,l,alpha)
softd1=wthresh(d1,'s',2.4834);
softd2=wthresh(d1,'s',1.8841);
softd3=wthresh(d1,'s',1.1950);
c3=[a3,softd3,softd2,softd1];
s3=waverec(c3,l,wmane);
subplot(2,2,4);
plot(s3);
title('給定軟閾值消噪后的信號');