% 當前延拓模式是補零
% 裝載原始圖像
load sinsin 
% X 包含原始圖像
subplot(221),image(X),colormap(map);
title('原始圖像');
% 產生噪聲圖像
init=2055615866; randn('seed',init); 
x = X + 18*randn(size(X));
subplot(222),image(x),colormap(map);
title('噪聲圖像');
% 使用wdencmp進行圖像降噪
% 尋找默認值
[thr,sorh,keepapp] = ddencmp('den','wv',x);
% 使用全局閾值進行圖像降噪
xd = wdencmp('gbl',x,'sym4',2,thr,sorh,keepapp);
subplot(223),image(xd),colormap(map);
title('降噪后圖像');