%裝載原始信號(hào)并圖示之
load noismima;
s=noismima(1:1000);
figure(1);
subplot(2,2,1);
plot(s);
xlabel('樣本序號(hào) n');
ylabel('幅值 A');
title('原始信號(hào)');

%采用默認(rèn)閾值、用wdencmp函數(shù)進(jìn)行消噪處理
[thr,sorh,keepapp,crit]=ddencmp('den','wp',s);

%用全局閾值選項(xiàng)進(jìn)行消噪處理
[c,treed,perf0,perfl2]= wpdencmp(s,sorh,3,'db2',crit,thr,keepapp);
subplot(2,2,3);
plot(c);
xlabel('樣本序號(hào) n');
ylabel('幅值 A');
title('默認(rèn)閾值消噪信號(hào)');

%根據(jù)前面的消噪效果，調(diào)節(jié)閾值大小進(jìn)行消噪
thr=thr+15;
[c1,treed,perf0,perfl2]= wpdencmp(s,sorh,3,'db2',crit,thr,keepapp);
subplot(2,2,4);
plot(c1);
xlabel('樣本序號(hào) n');
ylabel('幅值 A');
title('調(diào)節(jié)閾值后的消噪信號(hào)');