第3章 DFT及其FFT(12上).ppt

SCHOOL OF PHYSICS AND TECHNOLOGY N. N. U. 第 3 章离散傅里叶变换(DFT)及其快速算法(FFT) Discrete Fourier Transform (DFT) and Fast Fourier Transform (FFT) 本章主要内容 离散傅里叶变换（DFT）的定义 离散傅里叶变换的主要性质 频域采样 DFT的快速算法——快速傅里叶变换（FFT） DFT（FFT）应用举例 存在问题： 序列的傅里叶变换、Z变换是时域离散信号及系统分析与设计的重要数学工具；但变换结果均为连续函数，无法用计算机进行处理。 解决之道： 离散傅里叶变换（DFT）：对有限长时域离散信号的频谱进行等间隔采样，频域函数被离散化了，便于信号的计算机处理。 不过DFT运算量较大，快速离散傅里叶变换（FFT）算法是解决方案。 3.1 离散傅里叶变换（DFT） 的定义及物理意义 复习 周期序列的离散傅里叶级数（DFS） 3.1.1 DFT的定义 设x(n)是一个长度为M的有限长序列，x(n)的N点离散傅里叶变换（正变换－DFT）定义为： N称为DFT变换区间长度， N点离散傅里叶逆变换（IDFT）定义为 令 ，离散傅里叶变换与逆变换对可表示为： 逆变换的证明： 由于 于是 因此离散傅里叶逆变换是唯一的。 例3.1.1 分别计算序列的8点、16点DFT 解: 8点DFT 16点DFT 上图显示， 是 在频率区间 上的等间隔采样。 1. DFT与ZT、DTFT之间的关系 有限长序列 DFT与ZT、DTFT的关系： 是 在频率区间 上的等间隔采样。 DFT与DFS、ZT、DTFT之间的关系（续） 2. DFT与DFS之间的关系 有限长序列 周期序列：由 以N为周期进行周期延拓而成 周期序列的DFS 有限长序列的DFT DFT与DFS之间的关系: 有限长序列x(n)的DFT变换X(k)，就是x(n)的周期延拓序列　　 的DFS系数　　　的主值 区间序列 DFS与DTFT之间的关系: 周期延拓序列的频谱特性由傅里叶级数的系数 确定，幅度相差一个常数因子 DFT: 是 的主值区间序列，它表示了周期序列的频谱特性 3.1.3 DFT的矩阵方程表示 DFT的矩阵方程表示（续） IDFT的矩阵表示 3.1.4用MATLAB计算序列的DFT xn=[1 1 1 1 1 1 1 1]; %输入时域序列向量xn=R8(n) Xk32=fft(xn,32); %计算xn的32点DFT Xk64=fft(xn,64); %计算xn的64点DFT %以下为绘图部分 k=0:31;wk=2*k/32; %产生32点DFT对应的采样点频率(关于π归一化值) subplot(2,2,1);stem(wk,abs(Xk32),.); %绘制32点DFT的幅频特性图 title((a)32点DFT的幅频特性图);xlabel(ω/\pi);ylabel(幅度) subplot(2,2,3);stem(wk,angle(Xk32),.); %绘制32点DFT的相频特性图 title((b)32点DFT的相频特性图); xlabel(ω/\pi);ylabel(相位);axis([0,2,-3.5,3.5]) k=0:63;wk=2*k/64; %产生64点DFT对应的采样点频率(关于π归一化值) subplot(2,2,2);stem(wk,abs(Xk64),.); %绘制64点DFT的幅频特性图 title((c)64点DFT的幅频特性图);xlabel(ω/\pi);ylabel(幅度) subplot(2,2,4);stem(wk,angle(Xk64),.); %绘制64点DFT的相频特性图 title((d)64点DFT的相频特性图) xlabel(ω/\pi);ylabel(相位);axis([0,2,-3.5,3.5]) 小结 DFT 引入的目的 DFT 的定义 DFT与DFS、ZT、DTFT之间的关系 DFT、IDFT的计算 DFT的矩阵表示 3.2 离散傅里叶变换（DFT） 的主要性质 1. 线性性质 设x1(n)，x2(n)为有限长序列，长度分别为