DFT离散傅里叶变换,按照公式编写程序就是了。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
>> n=1:50;
x=cos(n*(pi/6));
N=12;
X=zeros(1,N);
sum=0;
for k=1:N
for ii=1:N
sum=sum+x(ii)*exp(-i*2*pi*k/N);
end
X(k)=sum;sum=0;
end
stem(abs(X));
grid on
在Matlab中执行X=fft(x(n),N)语句,即可得到x(n)=cos(npi/6)的N点傅里叶变换后之后再画图即可。fft是DFT的快速形式 没有本质差别。