program_13_07.m

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

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