program_13_22.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"MATLAB在小波分析中的应用"

% 产生heavy sine信号和噪声
init = 2055615866;
[xref,x] = wnoise(5,11,7,init);
xref=xref(1:1000);
x=x(1:1000);
% 画出原始信号和噪声信号
subplot(221), plot(xref),
title('原始信号')
subplot(222), plot(x),
title('噪声信号')
% 使用wpdencmp函数进行信号降噪
n = length(x); 
thr = sqrt(2*log(n*log(n)/log(2))); 
xwpd = wpdencmp(x,'s',4,'sym4','sure',thr,1);
% 画出小波包降噪信号
subplot(223), plot(xwpd),
title('小波包降噪信号')
% 比较小波基降噪结果
xwd = wden(x,'rigrsure','s','one',4,'sym4');
% 画出小波降噪信号
subplot(224), plot(xwd),
title('小波降噪信号')