clc;
clear;
% 加载信号
load noisdopp;
s = noisdopp;
subplot(321);plot(s);title('原始信号');
%%
[thr,sorh,deepapp,crit]=ddencmp('den','wp',s);
[x,~,~,~]=wpdencmp(s,sorh,3,'sym6',crit,thr,deepapp);
subplot(322);plot(x);title('去噪信号(3层)');
%%
[thr,sorh,deepapp,crit]=ddencmp('den','wp',s);
[x1,~,~,~]=wpdencmp(s,sorh,4,'sym6',crit,thr,deepapp);
subplot(323);plot(x1);title('去噪信号(4层)');
%%
[thr,sorh,deepapp,crit]=ddencmp('den','wp',s);
[x2,~,~,~]=wpdencmp(s,sorh,5,'sym6',crit,thr,deepapp);
subplot(324);plot(x2);title('去噪信号(5层)');
%%
[thr,sorh,deepapp,crit]=ddencmp('den','wp',s);
[x3,~,~,~]=wpdencmp(s,sorh,8,'sym6',crit,thr,deepapp);
subplot(325);plot(x3);title('去噪信号(8层)');