wadec.m

这是偶自己的作品拉

X=imread('E:\matlab1\LENAcaise64.bmp','bmp');
Wavelet_Name='db1';
Wavelet_Scale=1;
[C,S]=wavedec2(double(X(:,:,1)),Wavelet_Scale,Wavelet_Name);
LL1=appcoef2(C,S,Wavelet_Name,Wavelet_Scale);
LH1=detcoef2('v',C,S,1);
HL1=detcoef2('h',C,S,1);
HH1=detcoef2('d',C,S,1);
K=0;
HL1=HL1*K;
LH1=LH1*K;
HL1=HL1*K;
C5= [LL1 HL1 LH1 HH1];
S1=S;
Y1=waverec2(C5,S1,Wavelet_Name);
X(:,:,1)=Y1;

[C1,S1]=wavedec2(double(X(:,:,2)),Wavelet_Scale,Wavelet_Name);
LL2=appcoef2(C1,S1,Wavelet_Name,Wavelet_Scale);
LH2=detcoef2('v',C1,S1,1);
HL2=detcoef2('h',C1,S1,1);
HH2=detcoef2('d',C1,S1,1);
K=0;
LL2=LL2*K;
LH2=LH2*K;
HL2=HL2*K;
C6= [LL2 HL2 LH2 HH2];
Y2=waverec2(C6,S1,Wavelet_Name);
X(:,:,2)=Y2;

[C2,S2]=wavedec2(double(X(:,:,3)),Wavelet_Scale,Wavelet_Name);
LL3=appcoef2(C2,S2,Wavelet_Name,Wavelet_Scale);
LH3=detcoef2('v',C2,S2,1);
HL3=detcoef2('h',C2,S2,1);
HH3=detcoef2('d',C2,S2,1);
K=0;
HL3=HL3*K;
LH3=LH3*K;
HL3=HL3*K;
C7= [LL3 HL3 LH3 HH3];
Y3=waverec2(C7,S1,Wavelet_Name);
X(:,:,3)=Y3;
figure(100);imshow(X);
RGB=cat(3,Y3,Y2,Y1);
figure(200);
imshow(RGB);