wav1.m

来自「Included are the files wav1.m, wav2.m, w」· M 代码 · 共 37 行

%wav1.m
%takes the wavelet transform of the input matrix
%and gives as an output the four submatrices
%daubechies wavelet coefficients are used as filter coefficients
function [out]=wav1(inp)
load wavcoef.mat;
[N1 N2]=size(inp);
%******************
rs=conv2(inp,Hll); 
rs1=rs(8:N1+7,8:N1+7);
rs=conv2(inp,Hlh);
rs2=rs(8:N1+7,8:N1+7);
rs=conv2(inp,Hhl);
rs3=rs(8:N1+7,8:N1+7);
rs=conv2(inp,Hhh);
rs4=rs(8:N1+7,8:N1+7);
%******************   
for k1=1:N1/2,
  rs11(k1,:)=rs1(2*k1-1,:);
  rs22(k1,:)=rs2(2*k1-1,:);
  rs33(k1,:)=rs3(2*k1-1,:);
  rs44(k1,:)=rs4(2*k1-1,:);
end;
for k1=1:N1/2,
  rs111(:,k1)=rs11(:,2*k1-1);
  rs222(:,k1)=rs22(:,2*k1-1);
  rs333(:,k1)=rs33(:,2*k1-1);
  rs444(:,k1)=rs44(:,2*k1-1);
end;
 
%******************
out(1:N1/2,1:N1/2)=rs111;
out(N1/2+1:N1,1:N1/2)=rs222;
out(1:N1/2,N1/2+1:N1)=rs333;
out(N1/2+1:N1,N1/2+1:N1)=rs444;
%******************