work1.m

是個能做復立葉轉換並將矩陣轉換成實部和虛部的部份~相當的方便

A=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0;
   0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0; 
   0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0; 
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0; 
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];


Fourr=zeros(16);
Fouri=zeros(16);
ifour=zeros(16);
ifoui=zeros(16);
ifoui1=zeros(16);

for k=-7:8
   for l=-7:8
      for m=7:9
         for n=7:9
            sum=((k-1)*(m-1)+(l-1)*(n-1))/16;
            a=2*pi*sum;
            b=cos(a);
            c=sin(a);
            p=A(m,n)*b;
            q=-A(m,n)*c;
            Fourr(k+8,l+8)=Fourr(k+8,l+8)+p;
            Fouri(k+8,l+8)=Fouri(k+8,l+8)+q;
         end
      end   
   end   
end

for rk=1:16
   for rl=1:16
      for rm=-7:8
         for rn=-7:8
            rsum=(((rk-1)*(rm-1))+((rl-1)*(rn-1)))/16;
            ra=2*pi*rsum;
            rb=cos(ra);
            rc=sin(ra);
            w1=Fourr(rm+8,rn+8)*rb;
            w2=Fourr(rm+8,rn+8)*rc;
            w3=Fouri(rm+8,rn+8)*rb;
            w4=Fouri(rm+8,rn+8)*rc;
            rp=w1-w4;
            rq=w2+w3;
            ifour(rk,rl)=ifour(rk,rl)+rp;
            ifoui(rk,rl)=ifoui(rk,rl)+rq;
            ifoui1(rk,rl)=ifoui1(rk,rl)+rq;
         end
      end   
   end   
end
for i=1:16
   for j=1:16
      ifour(i,j)=abs(ifour(i,j)/256);
      ifoui(i,j)=ifoui(i,j)/65536;
      ifoui1(i,j)=fix(ifoui1(i,j)/256);
   end
end   
disp(Fourr);
disp(Fouri);
disp(ifour);
disp(ifoui);
disp(ifoui1);