离散傅里变换(DFT)
第3章 离散傅里叶变换(DFT) 3.1 离散傅里叶变换的定义及物理意义 3.2 离散傅里叶变换的基本性质 3.3 频率域采样 3.4 DFT的应用举例 3.1 离散傅里叶变换的定义及物理意义 3.1.1 DFT (Discrete Fourier Transform)的定义 设x(n)是一个长度为M的有限长序列,定义x(n)的N点 离散傅里叶变换为: 证明IDFT[X(k)]的唯一性: 把(3.1.1)式代入(3.1.2)式有 例 3.1.1 x(n)=R4(n) ,求x(n)的4点、8点、16点DFT 解:设变换区间N=4,则 3.1.2 DFT与傅里叶变换(DTFT)和Z变换的关系 设序列x(n)长度为M,其Z变换和N(N≥M)点DFT分别为: 3.1.3 DFT的隐含周期性 前面定义的DFT变换对中,x(n)与X(k)均为有限长序列,但由于 的周期性,使离散傅里叶变换式中的x(n)与X(k)隐含周期性,且周期均为N。对任意整数m, 总有 任何周期为N的周期序列 都可以看作长度为N的有限长序列x(n)的周期延拓序列,而x(n)则是 的一个周期,即 DFT的矩阵方程表示 作 业 第三章《习题与上机题》 课本第105、106页 1、 2、(题目有误) 作 业 第三章《习题与上机题》 课本第106页 3、 3.4 DFT的应用举例 DFT的快速算法FFT的出现,使DFT在数字通信、 语言信号处理、图像处理、功率谱估计、仿真、系统分析、雷达理论、光学、医学、地震以及数值分析等各个领域都得到广泛应用。 各种应用的基础: (1)卷积、相关的具体计算 (2)以DFT作为连续傅里叶变换的近似 作 业 第三章《习题与上机题》 课本第108页 18、 19、 所以,上述分析方法不丢失信息,即可由 X(k) 恢复 Xa(jΩ) 或 xa(t)。 直接由分析结果 X(k) 看不到 Xa(jΩ) 的全部频谱特性,而只能看到 N 个离散采样点的谱线,这就是所谓的栅栏效应。 如果 xa(t) 持续时间无限长,上述分析中要进行截断处理,所以会产生所谓的截断效应,从而使谱分析产生误差。 理想低通滤波器的单位冲激响应 ha(t) 及其频响函数 Ha(f) 如图3.4.7(a)、(b)所示(图3.4.7(a)中只画出ha(t)所截取的一段)。图中, 图3.4.7 用DFT计算理想低通滤波器的频响曲线 采用DFT来分析 ha(t) 的频率响应特性: 由于 ha(t) 的持续时间为无穷长,因此要截取一段Tp,假设 Tp = 8 s,采样间隔 T=0.25 s (即采样频率Fs=4 Hz),采样点数N=Tp/T=32;频域采样间隔F=1/Tp=0.125 Hz; ha(t)为实信号,仅取正频率[0, Fs/2]频谱采样: 其中 图3.4.7 用DFT计算理想低通滤波器的频响曲线 式中, ,因此 (3.3.4b) 令 则 1 = - kN N W 上式 即为用 X(k) 表示 X(z) 的内插公式 φk(z)称为内插函数 将z=ejω代入(3.3.4a)式,并进行化简,可得 (3.3.7) 在数字滤波器的结构与设计中,我们将会看到,频域采样理论及有关公式可提供一种有用的滤波器结构和滤波器设计途径,(3.3.7)式有助于分析FIR滤波器频率采样设计法的逼近性能。 (3.3.8) 【例3.3.1】 长度为26的三角形序列x(n)如图3.3.1(b)所示。编写MATLAB程序验证频域采样理论。 解:先计算x(n)的32点DFT,得到其频谱函数X(ejω)在频率区间[0,2π]上等间隔32点采样X32(k),再对X32(k)隔点抽取,得到X(ejω)在频率区间[0,2π]上等间隔16点采样X16(k)。最后分别对X16(k)和X32(k)求IDFT, 得到: 绘制x16(n)和x32(n)波形图验证频域采样理论。 % 频域采样理论验证 M=26;