wname = 'db5';
%设置小波名
[Lo_D,Hi_D,Lo_R,Hi_R] = wfilters(wname);
% 计算与给定小波名相关的四个滤波器
subplot(221); stem(Lo_D);
title('分解低通滤波器');
subplot(222); stem(Hi_D);
title('分解高通滤波器');
subplot(223); stem(Lo_R);
title('重构低通滤波器');
subplot(224); stem(Hi_R);
title('重构高通滤波器');
load sinsin
% 打开原始图像
% X中包含图像信息
init=2055615866; randn('seed',init);
x = X + 18*randn(size(X));
%产生噪声图像
% 用 wdencmp函数给图像除噪
[thr,sorh,keepapp] = ddencmp('den','wv',x);
% 发现缺省值
xd = wdencmp('gbl',x,'sym4',2,thr,sorh,keepapp);
% 使用全局阈值参数除噪图像
subplot(121)
imshow(x,map)
title('加噪图像')
subplot(122)
imshow(xd,map)
title('除噪图像')
MATLAB图像处理—孙兆林编著—清华大学出版社