MATLAB 离散傅里叶变换及其逆变换和快速傅里叶变换(DFT, IDFT, FFT)的实现

 clc;clear;

    
 xn=[7,6,5,4,3,2];
    
 Xk=dft(xn,6);
    
 x=idft(Xk,6);
    
  
    
 subplot(2,2,1);stem(0:5,abs(Xk),'filled');
    
 axis([0,5,0,1.1*max(abs(Xk))]);
    
 title('x[n]经DFT后的幅度');xlabel('频率w');ylabel('幅度');
    
  
    
 subplot(2,2,2);stem(0:5,angle(Xk),'filled');
    
 axis([0,5,1.1*min(angle(Xk)),1.1*max(angle(Xk))]);
    
 title('x[n]经DFT后的相角');xlabel('频率w');ylabel('相角');
    
  
    
 subplot(2,2,3);stem(0:5,xn,'filled');
    
 axis([0,5,0,1.1*max(xn)]);
    
 title('原信号x[n]');xlabel('时间n');ylabel('幅度')