《信号与系统——离散傅里叶变换讲义资料》由会员分享，可在线阅读，更多相关《信号与系统——离散傅里叶变换讲义资料（79页珍藏版）》请在金锄头文库上搜索。
1、1,第6章 离散傅里叶变换,6.3 离散时间系统的频域分析,6.2 离散时间傅里叶变换,6.4 离散傅里叶变换,6.5 信号频谱的数值计算,6.6 离散傅里叶变换的性质,6.7 快速傅里叶变换简介,6.1 引言,2,6.1 引 言,1、从连续信号频域分析看到： 从频域角度可以获得对LTI系统性质的更加深入的了解，使系统的分析与设计更加直观方便，与时域分析互补。,6.1 引言,3、本章学习注意： （1）与连续情况对应关系并找出相似之处和重要区别；,2、从第二章时域和第四章的复频域看到： 连续信号与系统和离散信号与系统之间可以通过抽样联系起来，二者在时域和复频域中均有对应关系。本章将会看到，二者在
2、频域之间也有对应关系。,3,（1）实际信号与计算机能处理的信号之间的矛盾；,实际信号的特点： 时域：连续时间信号，持续时间较长 频域：频谱连续 数字处理设备（计算机）的特点 存储空间有限：只能存储有限多数据（离散的数据点，有限长的时间范围） 表示空间有限：只能表示有限多的数值（取值在一定精度内，取值在一定范围内）,4、本章要解决的问题：,6.1 引言,（2）以理论分析为依据，以工程实现为目的。,5,6.2 离散时间傅里叶变换,6.2 离散时间傅里叶变换,连续非周期信号,，,FT,DTFT,IFT,IDTFT,连续,1、比较FT和DTFT,，,n 离散,连续,t 连续,离散非周期信号,6,6.2
3、 离散时间傅里叶变换,2、对DTFT的说明,频谱密度连续，,综合式IDTFT,序列离散,分析式DTFT,综合式IDTFT不是n的周期函数，,以 为周期 因为： n为整数,7,6.2 离散时间傅里叶变换,连续频谱密度,是积分式,有利因素：频谱密度可利用离散点的数据计算，为利用计算机提供了可能。,由于,周期,所以,不利因素：计算机无法直接处理和存储连续频谱，数字处理遇到困难。,积分范围是,是求和式, 由于,非周期 且非时限，,是离散序列,对无穷多项求和,所以,8, 若,对应dt,6.2 离散时间傅里叶变换,9,6.2 离散时间傅里叶变换,3、典型信号的DTFT,解：,10,解：,6.2 离散时间傅
4、里叶变换,解：,ROC边界在单位圆，信号不满足绝对可和的条件。但可以仿照连续时间信号情况，在变换中引入冲激函数。,由于离散时间信号的傅里叶变换是以 为周期的，考察下式给出的等间隔冲激频谱函数：,利用逆变换公式得,因此,11,解：由求解逆傅里叶变换的公式有,从图中可以看出，离散时间系统的理想低通滤波器的样值响应，与连续时间系统的理想低通滤波器的冲激响应类似，即在输入没有加入前就已有了响应。 说明离散时间系统的理想低通滤波器也是一个非因果系统。,6.2 离散时间傅里叶变换,例6-4 若离散时间系统的理想低通滤波器频率特性,如图所示，求它的逆傅里叶变换,（即单位样值响应）。,12,6.3 离散时间系
5、统的频域分析,6.3 离散时间系统的频域分析,1. 离散时间傅里叶变换的性质,2. 离散时间系统的频域分析,13,1. 离散时间傅里叶变换DTFT的性质,（1）周期性、连续性 周期性：是离散信号 的DFS、DTFT的共性， 有别于 的FS、FT。,注意与FT对应，与Z变换对应,6.3 离散时间系统的频域分析,1. 离散时间傅里叶变换的性质,14,（2）线性（所有线性变换的共性）,设 、 的傅里叶变换分别为 及,其中 、 为任意常数,则,6.3 离散时间系统的频域分析,1. 离散时间傅里叶变换的性质,（3）移位（与FT一致）,若 是傅里叶变换对，则有：,频域移位：,时域移位：,15,（4）时域线
6、性加权（频域微分）（与FT一致）,若 是傅里叶变换对，则,（5）反转与对称（与FT一致）,若 是傅里叶变换对，则,实偶,实奇,对称,为实序列,实偶,虚奇,1. 离散时间傅里叶变换的性质,6.3 离散时间系统的频域分析,16,（6）卷积定理（与FT同）,若 ， ，则,时域卷积：,频域卷积：,1. 离散时间傅里叶变换的性质,6.3 离散时间系统的频域分析,17,（7）帕斯瓦尔定理,若 是傅里叶变换对，则,即时域的全部信息量包含在频谱的一个周期内，所以只讨论频谱的一个周期就够了。,1. 离散时间傅里叶变换的性质,6.3 离散时间系统的频域分析,18,2. 离散时间系统的频域分析,6.3 离散时间系统
7、的频域分析,2. 离散时间系统的频域分析,时域：,频域：,输入信号的频谱函数 经系统后变为,在输入信号频谱给定的情况下，要想得到需要的输出频谱结构的过程，实际上是对H()进行设计的过程。在频域中通过输入输出信号的频谱可清晰地看到系统对信号每个分量的变换过程及对H()的要求。,19,6.3 离散时间系统的频域分析,2. 离散时间系统的频域分析,例6-5 求差分方程 所描述系统的频响函数H()。若输入为 ，求响应 yn,解：对差分方程两边取傅里叶变换，有,由频响函数的定义可知,由于,则,取逆变换得,20,6.4 离散傅里叶变换,6.4 离散傅里叶变换,1.离散傅里叶级数,2.离散傅里叶变换,21,
8、6.4 离散傅里叶变换,1. 离散傅里叶级数DFS,1. 离散傅里叶级数,DFS 离散周期信号,完备正交系,t 连续 k:(-,),级数,FS 连续周期信号,整周期抽样,（1）由FS引入DFS,22,系数,理论依据:,=,B=T,B=N,（2）对DFS的说明,一个周期主值,一个周期主值,6.4 离散傅里叶变换,1. 离散傅里叶级数DFS,23,一个周期主值,一个周期主值,周期N点,是周期函数，定义域k:,周期N点,重要意义：只要计算一个周期的N个点，即可得到全 域结果。,6.4 离散傅里叶变换,1. 离散傅里叶级数DFS,在一个周期内有N个谐波分量，第k个谐波分量为：,24,因为 离散， 所以
9、 周期,即：,故适合计算机工作：离散 有限数, DFS总是收敛，因为是有限数项的求和。,6.4 离散傅里叶变换,1. 离散傅里叶级数DFS,25,解：由于所给信号的数字频率为 ，则该信号的 周期N为,把余弦函数用指数函数表示,由于,于是,根据,n域与k域的周期相同,6.4 离散傅里叶变换,1. 离散傅里叶级数DFS,26,此信号的频谱是以N=6为周期的周期离散频谱。,可得,6.4 离散傅里叶变换,1. 离散傅里叶级数DFS,27,例6-7 图（a）所示序列的周期N =10，求其频谱。,解：,6.4 离散傅里叶变换,1. 离散傅里叶级数DFS,28, 计算机可以处理的数据形式,6.4 离散傅里叶
10、变换,2. 离散傅里叶变换DFT,2. 离散傅里叶变换,离散：数据离散地放在存储器的各个单元,有限：存储空间有限，计算速度有限,为了让计算机解决实际问题，必须要做的工作是：,（1）从理论研究到工程实际,29,6.4 离散傅里叶变换,时域 是信号（输入，输出）的原始形式，离散化为 ， ， （第五章已解决）。,（a）离散处理,频域 是系统设计的出发点，也应离散化为 （本节将解决）。,在时域对离散化后的序列截断或加窗，在频域对离散信号的频谱（周期）加窗即只取用一个周期。,（b）有限化处理（本节解决）,2. 离散傅里叶变换DFT,30,6.4 离散傅里叶变换, 已有的理论基础,时域 频域,2. 离散傅
11、里叶变换DFT,31,6.4 离散傅里叶变换, 分析离散傅里叶级数DFS（时频域均只需有限的离散数据）,离散,周期：频谱重复，一个周期内仅有N个点,（b） 对整个时域，只需要计算N个点，即可周期延拓到全时域 ；对整个频域，只需要计算N个点，即可周期延拓到全域 。,（a）对时域的每个第n点，只需要N个频域数据参与计算，对频域的每个第k点，只需要N个时域数据参与计算,2. 离散傅里叶变换DFT,32,6.4 离散傅里叶变换,（2）离散傅里叶变换DFT 从DFS到DFT (a) DFS的分析和综合式中实际参加计算的数据分别是：,它们是周期数据 在主值区间内的数据。对周期 延拓可得 ，对周期 延拓可得
12、 。这是在保证信息不损失条件下最少的数据量。,条件： 是时间有限序列,(b) 若 是非周期序列，且为时间有限（N点），即可把它看成 的主值数据，仅用这N个点即可计算出 对应的主值 。,也就是说：直接对非周期序列,2. 离散傅里叶变换DFT,33,6.4 离散傅里叶变换, 定义DFT,DFT,IDFT,注意： (a) 标度因子N换位，这是为了让正变换简单，与DTFT对应。,DTFT,IDTFT,(b) 是非周期N点时限信号，可看成周期信号的一个主周期。,(c) 是对有限长的 的频谱 等间隔 抽样并对主值区间加窗得到的。,2. 离散傅里叶变换DFT,34,展开成矩阵形式：,DFT,DTFT,6.4
13、 离散傅里叶变换,DFT可简写为：,(d) 频率间隔 对应的 是 的最小单位（基频分量），令 （对应z变换下的 ）,2. 离散傅里叶变换DFT,35,第6章 离散傅里叶变换,6.3 离散时间系统的频域分析,6.2 离散时间傅里叶变换,6.4 离散傅里叶变换,6.5 信号频谱的数值计算,6.6 离散傅里叶变换的性质,6.7 快速傅里叶变换简介,6.1 引言,36,6.5 信号频谱的数值计算,6.5 信号频谱的数值计算,1.周期信号的频谱分析,2.非周期信号的频谱分析,3.数据截断问题,37,1. 周期信号的频谱分析,6.5 信号频谱的数值计算,1. 周期信号的频谱分析,连续时间周期信号xT(t)
14、的样本xNn与xT(t)的离散频谱ck的关系：,对频带有限的xT(t)，若一个周期中抽样的离散点数N大于最高谐波次数km的二倍即：N 2km，则在忽略数值误差的情况下，ck可以精确计算：,38,1. 周期信号的频谱分析,6.5 信号频谱的数值计算,若xT(t)频谱无限分布，用DFT必然会出现频普混叠，从而带来误差。解决的方法：适当提高抽样频率，以减少频谱混叠的影响。一般取两种不同的抽样频率进行计算，当二者的计算频谱基本一致时可认为结果正确。,根据ck 求xT(t)在主值区间离散值：,39,2. 非周期信号的频谱分析,（1）DFT与DTFT的关系, 有限长,DFT 离散变化，共N个点,DTFT
15、连续变化，周期为,结论：DFT是对DTFT在频域 内取N个点等间隔抽样的结果，DFT的包络线即为DTFT。,2. 非周期信号的频谱分析,6.5 信号频谱的数值计算,40,DFT：,DTFT：,对DTFT连续的频谱离散化的结果 MATLAB实现：利用FFT计算 ， 用绘图语句stem绘出 的离散频谱图。,例6-9 非周期离散序列,连续的频谱，是DFT的包络线,MATLAB实现：利用FFT计算 ， 用绘图语句plot绘出 的包络线，即 的连续频谱图。,2. 非周期信号的频谱分析,6.5 信号频谱的数值计算,41,在分析信号频谱的时候，由于受到计算能力的影响，只能处理有限长的信号。这就必须截取时间函数的一个有限范围，即把观测到的信号限制在一定的时间间隔之内。换句话说，就是要取出信号的某一个时间段。这种过程就是截断数据的过程。这种截断过程相当于对信号进行加窗，即信号乘以窗函数 ，变成 的N点有限长序列，然后可以利用DFT计算DTFT，N的大小会影响结果的准确性，应视情况慎重选择。, 无限长 、,N越大，间隔越密，在 给定的情况下，可补零加大N。,（2）N与谱线间隔,DFT对DTFT离散化的频谱间隔,2. 非周期信号的频谱分析,6.5 信号频谱的数值计算, DFT与FT的关系 （a）关系,42,DFT与FT
