f=[4,2;3,4];
[M,N]=size(f);
F=zeros(M,N);
j=sqrt(-1);
for u=1:M
for v=1:N
for h=1:M
for k=1:N
F(u,v)=F(u,v)+f(h,k)*exp(-2*pi*j*((u-1)*(h-1)/M+(v-1)*(k-1)/N));
end
end
end
end
F
输出结果:
f=[4,2;3,4];
[M,N]=size(f);
F=zeros(M,N);
j=sqrt(-1);
for u=1:M
for v=1:N
for h=1:M
for k=1:N
F(u,v)=F(u,v)+f(h,k)*exp(-2*pi*j*((u-1)*(h-1)/M+(v-1)*(k-1)/N));
end
end
end
end
F
输出结果: