conv2.sci

小波分解源代码

function c=conv2(a,b)

//  Copyright Aldo I Maalouf
adims=ndims(a);
bdims=ndims(b);
if (adims>2)|(bdims>2),error('dimension of matrix should be less than 2');end

[ma,na] = size(a);
[mb,nb] = size(b);

tmpm=ma;
if ma>mb, tmp=mb;end

tmpn=na;
if na>nb, tmp=nb;end
c=zeros(ma+mb-1,na+nb-1);
for m=1:(ma+mb-1)
for n=1:(na+nb-1)

  c(m,n)=0;
  for j=1:tmpm
  for k=1:tmpn       
     if k>na|j>ma then c(m,n)=c(m,n);
     
     else if ((n-k+1)>nb)|(n-k+1)<1|(m-j+1)>mb|(m-j+1)<1,c(m,n)=c(m,n);
          else c(m,n)=c(m,n)+a(j,k)*b(m-j+1,n-k+1);
          end
      end
   end
   end
end
end

endfunction