下面是我编写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