图解处理的小波分析matlab源码

%程序algorithm2D.m %利用Daubechies小波（m=2)分解与重构256*256图像，把函数值作为起始%系数 %试验图像 picture=256*ones(256); for i=1:256 picture(i,i)=0; picture(i,257-i)=0; end image(picture) title('原始图像') %分解和重构滤波器 g0=[0.68301;1.18301;0.31699;-0.18301]; k=[0;1;2;3]; g1=flipud(g0).*(-1).^k; h0=flipud(g0)/2; h1=flipud(g1)/2; %分解过程 %第一层分解 for k=1:256 s=[0;0;picture(:,k);0;0]; x=conv(s,h0); a=x(1:2:length(x));%向下采样 x=conv(s,h1); w=x(1:2:length(x));%向下采样 C（:,k)=[a;w]; end for k=1:256+8 s=rot90([0 0 C(k,:)0 0],3); x=conv(s,h0); a=x(1:2:length(x));%向下采样 x=conv(s,h1); w=x(1:2:length(x));%向下采样 CC(k,:)=rot90[a;w]); end LL=CC(1:132,1:132)； HL=CC(133:264,1:132); LH=CC(1:132,133:264); HH=CC(133:264,133:264); figure(2) axes('位置',[0.1 0.5 0.3 0.3]) image(LL) title('LL') axes('位置',[0.5 0.5 0.3 0.3]) image(LH) title('LH') axes('位置',[0.1 0.1 0.3 0.3]) image(HL) title('HL') axes('位置',[0.5 0.1 0.3 0.3]) image(HH) title('HH') clear C clear CC %第二层分解 for k=1:132 s=LL(:,k); x=conv(s,h0); a=x(1:2:length(x));%向下采样 x=conv(s,h1); w=x(1:2:length(x));%向下采样 C（:,k)=[a;w]; end for k=1:128+8 s=rot90([C(k,:),3); x=conv(s,h0); a=x(1:2:length(x));%向下采样 x=conv(s,h1); w=x(1:2:length(x));%向下采样 CC(k,:)=rot90[a;w]); end LL_LL=CC(1:68,1:68)； HL_LL=CC(69:136,1:68); LH_LL=CC(1:68,69:136); HH_LL=CC(69:136,69:136); clear C clear CC %重构过程 %第二层重构 s=[LL_LL LH_LL;HL_LL HH_LL]; for k=1:136 x=zeros(136,1); x(1:2:136)=rot90(s(k,1:68),3); y=zeros(136,1); y(1:2:136)=rot90(s(k,69:136),3); x=con(x,g0)+conv(y,g1); C(k,:)=rot90(x(4:length(x)-4)); end s=C; clear C for k=1:132 x=zeros(136,1); x(1:2:136)=s(1:68,k); y=zeros(136,1); y(1:2:136)=s(69:136,k); x=conv(x,g0)+conv(y,g1); C(k,:)=x(4:length(x)-4); end LL_rec=C; clear C %第一层重构 s=[LL_rec LH;HL HH]; for k=1:264 x=zeros(264,1); x(1:2:264)=rot90(s(k,1:132),3); y=zeros(264,1); y(1:2:264)=rot90(s(k,133:264),3); x=conv(x,g0)+conv(y,g1); C(k,:)=rot90(x(4:length(x)-4)); end s=C; clear C for k=1:260 x=zeros(264,1); x(1:2:264)=s(1:132,k); y=zeros(264,1); y(1:2:264)=s(133:264,k); x=conv(x,g0)+conv(y,g1); C(k,:)=x(4:length(x)-4); end picture_rec=C(3:258,3:258); figure(3) image(picture_rec) title('重构图像')