program_13_19.m

来自「基于小波分析理论研究的运用matlab R2007实现的程序代码」· M 代码 · 共 22 行

% 装载噪声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('降噪后信号')