load sinsin 
% 打開原始圖像
% X中包含圖像信息
init=2055615866; randn('seed',init); 
x = X + 18*randn(size(X));
%產(chǎn)生噪聲圖像
% 用 wdencmp函數(shù)給圖像除噪
[thr,sorh,keepapp] = ddencmp('den','wv',x);
% 發(fā)現(xiàn)缺省值
xd = wdencmp('gbl',x,'sym4',2,thr,sorh,keepapp);
% 使用全局閾值參數(shù)除噪圖像
subplot(121)
imshow(x,map)
title('加噪圖像')
subplot(122)
imshow(xd,map)
title('除噪圖像')