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

%對圖像X用小波db5進行2層分解
[c,l]=wavedec2(X,2,'db5');
Csize=size(c);

%弱化不重要的分解系數(shù)
for i=1:Csize(2)
    if(c(i)>300)
        c(i)=2*c(i);
    else
        c(i)=0.5*c(i);
    end
end

%重構(gòu)圖像并顯示
X1=waverec2(c,l,'db5');
subplot(1,2,2);
image(X1);
colormap(map);
title('增強后圖像');