离散傅立叶变换-西南交通大学教务网.ppt

DSP 第三章 离散傅立叶变换(Discrete Fourier Transform ) 本章要点 离散傅立叶变换（Discrete Fourier Transform）DFT是有限长序列的一种傅立叶表示法。 时域T和频域ω皆离散的一种变换。 Fast Fourier Transform ,FFT是DFT的一种计算机算法实现。 本章内容 一、引言 傅立叶变换 傅立叶级数 序列的傅立叶变换 离散傅立叶变换 四种形式归纳 清议2009-11-01 ：向钱学森学习科学法则:28年前(注：1981年)，当我还是个懵懵懂懂的学生时，曾斗胆向内心崇敬的钱学森教授邮寄过一篇探讨系统科学应用于农业经济的万言文，恳求老先生百忙中拨冗赐教。没想到，一个多月后的一天，我居然收到了钱学森教授字迹工整、笔划略微颤抖的亲笔回信。当时的心情格外激动。28年后，当得知钱学森教授已驾鹤西去的消息时，我在沉痛悼念老先生一生科学救国、功勋万代之余，因当年承蒙教诲而感慨万千。 ?“任何学科的研究对象都是十分特定的”。这句话，是我从钱学森教授那封亲笔回信中感悟最深的。虽然经济学在方法论上无法与物理学相比，但这句话被我一直当作研究经济问题必须遵循的科学法则。? 打那以后，我渐渐深信，能够被冠以科学的学问，一定是杜绝了任何形式的空泛，对研究对象或问题本质的把握恰如其分，容不得丝毫马虎。 王大麻子2009-11-01:一个人拥有知识和爱国情怀之后，拥有良知最重要。。。正是那些有良知的思想家和科学家，使人类的星空灿烂无比，使人类能够达到更高的道德高度。 二、离散傅立叶级数（DFS） 周期复指数序列 三、离散傅立叶变换DFT 余数运算表达式 定 义 DFT的性质 1、DFT的线性 2、圆周移位 3、共轭对称性 圆周共轭对称 圆周共轭反对称 几种特例 4、DFT形式下的帕塞瓦定理 5、圆周卷积 图解 6、有限长序列的线性卷积与圆周卷积 其他 四、频域采样定理 采样Z变换 频域采样恢复 内插函数的特性 频响特性 续1 图示 总 结 总 结 五、利用DFT对连续时间信号的逼近（选讲） 图例1 图例2 图例3 续2 利用DFT计算连续时间信号的几个问题 混叠现象 例 频率泄漏 续1 栅栏效应 频率分辨力 分析 对连续信号xa(t)，采样频率fs，采样间隔Ts=1/fs 对实部： 0 N-1 N-1 k=0 称 Xep(k)为X(k)的共轭偶部(圆周共轭对称分量)。 对虚部同理可证： 称Xop(k)为X(k)的共轭奇部(圆周共轭反对称分量)。 0 N-1 N-1 k=0 1)当x(n)为实序列时，X(k)=Xep(k)有： 2)当x(n)为纯虚序列时，X(k)=Xop(k) 利用对称性只需计算X(0)～X(N/2-1)的值即可。 3)当x1(n)和x2(n)都是N点实序列时，构造新序列： 则： 因此，通过一次N点DFT运算完成了两个N点序列的DFT计算。 [证明]: 当x (n)＝y(n) 即序列x(n)在时域计算的能量与频域计算的能量相等。 设 x1(n)和x2(n)均为长度为N的有限长序列，且： 若： 则： N N 圆周卷积可看作延拓序列周期卷积后取主值区间而得 m m m m N-1 0 n N-1 n N 结果： 0 2 3 3 2 1 1 N-1 n N 在圆周上的操作图示如下： 离散线性系统输出响应(线性卷积) 可用圆周卷积代替线性卷积，减小运算量。 当x1(n)和x2(n)为有限长序列，其长度为： 对x1(n),x2(n)构造新序列，长度为L，两序列补零对齐 圆周卷积是线性卷积的周期延拓序列的主值序列。 为了不发生混叠，用圆周卷积代替线性卷积必须满足条件： L点DFT L点DFT L点IDFT 计算流程： 通常情况下，DFT可用FFT实现，而H(k)可离线计算好，上述算法比直接计算线性卷积快。 DFT的主要性质可参见教材 Page110,Table3-3：DFT的性质。 时域抽样: 对一个频带有限的信号,根据采样定理对其进行采样,所得采样序列的频谱是原带限连续信号频谱的周期延拓,主要满足奈奎斯特采样定理，采样信号的频谱不发生混叠，可完全不失真由采样序列恢复原信号。 频域抽样: 对有限长序列进行DFT所得X(k) 是序列傅氏变换的采样，故DFT就是频域抽样。这种频域采样的采样需要满足怎样的条件？是如何恢复连续频谱的？如何才能不失真的恢复呢？ 对X(z)在单位圆等距采样有： 当k取整数时，显然由于