如题所述
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;
stem(abs(X));
grid on