小波变换matlab代码,matlab小波变换代码

clear;clc;

T=256;       %  图像维数

SUB_T=T/2;   %  子图维数

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%  1.调原始图像矩阵

load wbarb%  下载图像

f=X;         %  原始图像

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%  2.进行二维小波分解

l=wfilters('db4','l');    %  db10(消失矩为10)低通分解滤波器冲击响应(长度为20)

L=T-length(l);

l_zeros=[l,zeros(1,L)];    %  矩阵行数与输入图像一致，为2的整数幂

h=wfilters('db4','h');    %  db10(消失矩为10)高通分解滤波器冲击响应(长度为20)

h_zeros=[h,zeros(1,L)];    %  矩阵行数与输入图像一致，为2的整数幂

for i=1:T;   %  列变换

row(1:SUB_T,i)=dyaddown( ifft( fft(l_zeros).*fft(f(:,i)') ) ).';    %  圆周卷积FFT

row(SUB_T+1:T,i)=dyaddown( ifft( fft(h_zeros).*fft(f(:,i)') ) ).';  %  圆周卷积FFT

end;

for j&