%裝載并顯示原始圖像
load flujet;
subplot(1,2,1);
image(X);colormap(map);
title('原始圖像');
%=====================================
%對圖像X用小波db3進行2層分解
[c,l]=wavedec2(X,2,'db3');
Csize=size(c);
%=====================================
%對分解系數作處理以突出所需部分并弱化不需要部分
for i=1:Csize(2)
    if(c(i)>300)
        c(i)=2*c(i);
    else
        c(i)=0.5*c(i);
    end
end
%=====================================
%重構圖像并顯示
X1=waverec2(c,l,'db3');
subplot(1,2,2);
image(X1);colormap(map);
title('增強圖像');