program_13_19.m

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

% 装载噪声chirp信号
load noischir; x = noischir;
% 使用sym6对信号进行5层小波包分解
wname = 'sym6'; lev = 5;
tree = wpdec(x,lev,wname);
% 由第1层高频系数估计噪声的标准差，结点索引号为2
det1 = wpcoef(tree,2);
sigma = median(abs(det1))/0.6745;
% 使用wpbmpen选择阈值
alpha = 2;
thr = wpbmpen(tree,sigma,alpha)
% 上面的阈值，软阈值处理，保存低频，使用wpdencmp进行信号降噪
keepapp = 1;
xd = wpdencmp(tree,'s','nobest',thr,keepapp);
% 画出原始信号和降噪后的信号
figure(1)
subplot(211), plot(x),
title('原始信号')
subplot(212), plot(xd)
title('降噪后信号')