二维离散傅里叶变换程序错误

S0704082

小弟对二维离散傅里叶变换不是很懂，请教这个程序如何编？如果使用FFT2这个函数？
clear
N=16;
m=1:N;
n=m;
x=1-exp(m/N+n/N);
y=zeros(N,N);
for m=1:N
    for n=1:N
        y(:,m)=y(:,m)+x(:,n)*exp(-j*2*pi*(m-1)*(n-1)/N);
    end
end
for m=1:N
    for n=1:N
        y(m,:)=y(m,:)+x(n,:)*exp(-j*2*pi*(m-1)*(n-1)/N);
    end
end
mesh(abs(y));