function  b = fft_2d(a) %正序輸出,逆序輸出 頻率抽取的遞歸算法N=length(a);if( N ~=2 )    for k=1:N/2        c(k)=a(k)+a(k+N/2);        c(k+N/2)=( a(k)-a(k+N/2) )*exp(-j*2*pi/N*(k-1));    end    b(1:N/2)=fft_2d(c(1:N/2));    b(N/2+1:N)=fft_2d(c(N/2+1:N));else    b(1)=a(1)+a(2);    b(2)=a(1)-a(2);end