examp13_22.m

《小波分析理论与Matlab7实现例程》一书里面各章节的源码

clc;
clear;
% 产生heavy sine信号和噪声
init = 2055615866;
[xref,x] = wnoise(5,11,sqrt(7),init);
xref=xref(1:1000);
x=x(1:1000);
% 画出原始信号和噪声信号
subplot(221), plot(x),
title('原始信号')
subplot(222), plot(xref),
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('小波降噪信号')