program_13_07.m

来自「小波变换源代码集 这是葛哲学编著《小波分析理论与MATLAB R2007实现》」· M 代码 · 共 21 行

% 装载bumps的噪声信号
load noisbump; x = noisbump;
% 使用sym6执行5层小波分解
wname = 'sym6'; lev = 5;
[c,l] = wavedec(x,lev,wname);
% 由第1层的高频系数估计噪声标准差
sigma = wnoisest(c,l,1);
% 信号降噪时，采用调整参数选择全局阈值
alpha = 2;
thr = wbmpen(c,l,sigma,alpha)
% 使用软阈值和保存的低频信号，进行信号降噪
keepapp = 1;
xd = wdencmp('gbl',c,l,wname,lev,thr,'s',keepapp);
% 画出原始信号和降噪后的信号
figure(1)
subplot(211), plot(x), title('原始信号')
subplot(212), plot(xd), title('降噪后的信号')