%裝載采集的信號leleccum.mat
load leleccum;
%將信號中第2000到第3450個采樣點賦給s
indx=2000:3450;
s=leleccum(indx);
%畫出原始信號
subplot(2,2,1);
plot(s);
title('原始信號');
ylabel('幅值 A');
%用db1小波對原始信號進行3層分解并提取系數
[c,l]=wavedec(s,3,'db1');
a3=appcoef(c,l,'db1',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,'db1');
subplot(2,2,2);
plot(s1);
title('強制消噪后的信號');
xlabel('樣本序號 n');
ylabel('幅值 A');
%用默認閾值對信號進行消噪處理并圖示結果
%用ddencmp函數獲得信號的默認閾值
[thr,sorh,keepapp]=ddencmp('den','wv',s);
s2=wdencmp('gbl',c,l,'db1',3,thr,sorh,keepapp);
subplot(2,2,3);
plot(s2);
title('默認閾值消噪后的信號');
xlabel('樣本序號 n');
ylabel('幅值 A');
%用給定的軟閾值進行消噪處理
softd1=wthresh(d1,'s',1.465);
softd2=wthresh(d2,'s',1.823);
softd3=wthresh(d3,'s',2.768);
c2=[a3 softd3 softd2 softd1];
s3=waverec(c2,l,'db1');
subplot(2,2,4);
plot(s3);
title('給定軟閾值消噪后的信號');
xlabel('樣本序號 n');
ylabel('幅值 A');