m0806.m

matlab的经典讲义

%这是第3题的MATLAB指令
%装载原始信号
load noismima.mat;
s=noismima(1:1000);
%绘制信号图形
subplot(3,1,1);
plot(s);
title('原始信号');

%默认阈值、用wdencmp函数去噪
[thr,sorh,keepapp,crit]=ddencmp('den','wp',s);

%用全局阈值选项进行去噪
[c,treed,perf0,perfl2]=...
    wpdencmp(s,sorh,4,'db3',crit,thr,keepapp);
subplot(3,1,2);
plot(c);
title('默认阈值去噪信号');

%调节适当阈值去噪
thr=thr+15;
[c1,treed,perf0,perfl2]=...
    wpdencmp(s,sorh,4,'db3',crit,thr,keepapp);
subplot(3,1,3);
plot(c1);
title('调节阈值的去信号');