%裝載并顯示原始圖像
load flujet;
subplot(2,2,1);
image(X);
colormap(map);
title('原始圖像');
axis square;

%生成含噪圖像
init=2055615866;
randn('seed',init);
X1=X+20*randn(size(X));
subplot(2,2,2);
image(X1);
colormap(map);
title('含噪圖像');
axis square;

%用小波sym2對圖像X1進行一層小波包分解
T=wpdec2(X1,1,'sym2');

%設置閾值
thr=8.342;

%對圖像的小波包分解系數進行軟閾值量化
NT=wpthcoef(T,0,'s',thr);

%僅對低頻系數進行重構
X2=wprcoef(NT,1);

%畫出消噪后的圖像
subplot(2,2,3);
image(X2);
colormap(map);
title('消噪后的圖像');
axis square;