请帮忙看一下这个程序

t12211017

clear
N1=128;dx=1;dy=1;N2=128;a=0.0081;b=0.74;g=9.81;u=8.5233;%y=zeros(50,1);
%for ii=1:50
y1=randn(N1,N2);
ff=zeros(N1,N2);
for v=1:N2
   for n=1:N2
ff(:,v)=ff(:,v)+(y1(:,n)*exp(-j*2*pi*(v-1)*(n-1)/(N2*dx)));
   end
end
%ff=ff.';
for u=1:N1
   for m=1:N1
ff(u,:)=ff(u,:)+(y1(m,:)*exp(-j*2*pi*(u-1)*(m-1)/(N1*dx)));
   end
end
z=ff;
tt=fft2(y1);
我是求得二维傅立叶变换，但是得到的结果和直接通过matlab自带的函数求得的不一样，不知道哪里出错了，还请指教