式中-N/2≤m≤N/2;
N=256.
求各位有的，给与指点，谢谢！

[此贴子已经被作者于2007-5-10 20:42:00编辑过]

不好意思,刚弄错了

[此贴子已经被作者于2007-5-10 20:46:16编辑过]

就做一个循环实现啊！！

下面是我编写DFT的程序,你可以找DFT的公式来看一下,然后对着程序你应该能编出你自己的程序!!
自己试一下吧~

for u=0:M-1
for v=0:N-1
for x=0:M-1
for y=0:N-1
fft=cos(2*pi*(((u*x)/M)+((v*y)/N)))-j*sin...
(2*pi*(((u*x)/M)+((v*y)/N)));
%%上一句是e指数项通过欧拉公式来实现
Fsingle = double(f(y+1,x+1))*fft;
F(u+1,v+1) = F(u+1,v+1)+Fsingle;
end
end
end
end