MATLAB 离散傅里叶变换（DFT）、逆离散傅里叶变换（IDFT）、快速傅里叶变换（FFT）的实现

离散傅里叶变换（DFT）、逆离散傅里叶变换（IDFT）的实现

代码如下，其中xn为时序序列

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('幅度');

subplot(2,2,4);stem(0:5,abs(x),'filled');
axis([0,5,0,1.1*max(abs(x))]);
title('原信号x[n]经DFT与IDFT的信号');xlabel('时间n');ylabel('幅度');

function [Xk]=dft(xn,N)
n=0:N-1;
k=0:N-1;
Xk=xn*exp(-j*2*pi/N).^(n'*k);%n'为n的转置

function [xn]=idft(Xk,N)
n=0:N-1;
k