第三章离散傅立叶变换PPT讲稿.ppt

《第三章离散傅立叶变换PPT讲稿.ppt》由会员分享，可在线阅读，更多相关《第三章离散傅立叶变换PPT讲稿.ppt（48页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第三章离散傅立叶变换第1页，共48页，编辑于2022年，星期二离散傅立叶变换的几种可能形式离散傅立叶变换的几种可能形式傅立叶变换就是以时间为自变量的“信号”与以频率为自变量的“频谱”函数之间的一种变换关系，当自变量“时间”和“频率“取连续值或离散值时，就形成不同的形式的傅立叶变换对。第2页，共48页，编辑于2022年，星期二非周期的连续时间、连续频率非周期的连续时间、连续频率傅立叶变换傅立叶变换非周期连续时间信号x(t)和它的频谱密度函数X(j)构成的傅立叶变换对为正变换反变换以连续时间矩形脉冲为例：x(t)t(a)非周期连续时间函数x(t)X(j)(b)非周期连续频谱X(j)第3页，共48页

2、，编辑于2022年，星期二周期的连续时间、离散频率周期的连续时间、离散频率傅立叶级数傅立叶级数周期为T0的连续时间信号x(t)的傅立叶级数展开的系数为X(jk0)，构成的傅立叶变换对为：正变换反变换X(jk0)是以角频率0为间隔的离散函数，形成频域的离散频谱，0与时间信号的周期之间的关系为。傅立叶级数展开将连续时间周期函数分解为无穷多个角频率为0整数倍的谐波，k为各次谐波序号。T0-T0 x(t)(a)周期连续时间函数x(t)0X(jk0)(b)非周期离散时间函数X(jk0)图3.2周期连续时间函数及其傅立叶变换第4页，共48页，编辑于2022年，星期二非周期的离散时间、连续频率非周期的离散时

3、间、连续频率序列的傅立叶变换序列的傅立叶变换非周期离散时间信号的傅立叶变换就是序列的傅立叶变换，其变换对为正变换反变换式中是数字频率。如果序列x(n)是模拟信号x(t)经过抽样得到，抽样时间间隔为Ts，抽样频率为fS，抽样角频率为S=2/Ts，由于数字频率与模拟角频率之间的关系为=T，因此抽样数字频率S=STS，则上面的变换对也可写成:正变换反变换第5页，共48页，编辑于2022年，星期二仍以连续时间矩形脉冲为例:结果表明，时域的离散造成频域的周期延拓，而时域的非周期性对应与频域的连续性。tTSx(nT)(a)离散时间序列X(ejT)S-S0(b)序列的频谱图 离散时间序列及其傅立叶变换第6页

4、，共48页，编辑于2022年，星期二离散时间、离散频率离散时间、离散频率离散傅立叶变换离散傅立叶变换假如序列x(n)是模拟信号x(t)经过抽样得到，抽样时间间隔为Ts，则频率函数的周期为S=2/Ts；如果频率函数也是离散的，其抽样间隔为0，则时间函数的周期为=2/T。当时间函数序列一个周期内的抽样点数为N时，有 上式表明在频域中频谱函数的一个周期内的抽样点数也为N，即离散傅立叶变换的时间序列和频率序列的周期都是N，可以得到表示于一个周期内的常用的离散傅立叶变换对为正变换反变换第7页，共48页，编辑于2022年，星期二(a)周期离散时间序列tx(n)TST02T0-T0-2T000S-SX(k)

5、(b)周期离散时间序列的频谱图周期离散时间序列及其傅立叶变换第8页，共48页，编辑于2022年，星期二周期序列的离散傅立叶级（周期序列的离散傅立叶级（DFSDFS）周期序列周期序列一个周期为N的周期序列，对于所有n满足式中N为正整数。定义n=0到N-1的周期区间为的主值区间主值区间，而主值区间内的N个样本值组成的有限长序列称为的主值序列主值序列，即这一过程称为取主值序列。第9页，共48页，编辑于2022年，星期二对于一个有限长序列如将其以N为周期进行周期性延拓，则可得由于周期序列不是绝对可和的，无论z取任何值，其z变换都是不收敛的，即因此周期序列不能用z变换法或傅立叶变换来进行讨论。第10页，

6、共48页，编辑于2022年，星期二离散傅立叶级数离散傅立叶级数令，则DFS变换对可写成正变换 反变换 离散傅立叶级数表明是以N为周期的周期序列，其基波成分为，k次谐波成分为，为DFS的k次谐波分量的复系数。由于的周期性，当已知0N-1次谐波成分后，根据周期性就可以确定其余的谐波分量，因此，无论时域或频域中都只有N个序列值是独立的。第11页，共48页，编辑于2022年，星期二离散傅立叶级数的性质离散傅立叶级数的性质假定和是周期皆为N的两个离散周期序列，它们的DFS为、线性、线性式中为任意常数，可见由两个离散周期序列和线性组合成一个新的周期序列的DFS也是周期为N的离散周期序列。第12页，共48页

7、，编辑于2022年，星期二、移位特性、移位特性时域移位 频域移位 如果N,那么证明：第13页，共48页，编辑于2022年，星期二、时域卷积特性、时域卷积特性两个周期都为N的周期序列和，它们卷积的结果也是周期为N的周期序列，即 m的取值由0（N-1），因此称为周期卷积。05n000000555555mmmmmm111234图图 两个周期序列两个周期序列(N=6)的周期卷积过程的周期卷积过程第14页，共48页，编辑于2022年，星期二周期卷积与DFS的关系如下：设 若 则有 这就是时域卷积定理。这就是时域卷积定理。第15页，共48页，编辑于2022年，星期二证明：第16页，共48页，编辑于2022

8、年，星期二、频域卷积特性、频域卷积特性 对于时域周期序列的乘积，同样对应于频域的周期卷积。若 则第17页，共48页，编辑于2022年，星期二离散傅立叶变换离散傅立叶变换由于长度为N的有限长序列可以看作是周期是N的周期序列的一个周期，因此利用DFS计算周期序列的一个周期，就可以得到有限长序列的离散傅立叶变换设x(n)是长度为N的有限长序列，可以把它看作是周期为N的周期序列的一个主周期，而将看作是x(n)以N为周期进行周期延拓得到，即同理第18页，共48页，编辑于2022年，星期二离散傅立叶变换的正变换反变换第19页，共48页，编辑于2022年，星期二离散傅立叶变换的性质离散傅立叶变换的性质假定和

9、都是N点的有限长序列，有、线性、线性 若两个有限长序列和的线性组合为，则有 式中为任意常数。说明：（1）若和的长度均为N，则的长度为N；（2）若和的长度不等，的长度为N1，的长度为N2，则的长度为N=maxN1，N2，离散傅立叶变换的长度必须按N来计算。第20页，共48页，编辑于2022年，星期二、序列的圆周移位、序列的圆周移位有限长序列x(n)的圆周移位是以它的长度N为周期，将其延拓成周期序列 ，并将周期序列进行移位，然后取主值区间（n=0到N-1）上的序列值。因而一个有限长序列的右圆周移位定义为第21页，共48页，编辑于2022年，星期二x(n)x(n)Nx(n-2)Nx(n-2)NRN(

10、n)nnnn0000N-1N-1N-1N-1图3.6 序列的周期移位(N=6)第22页，共48页，编辑于2022年，星期二（）时域移位定理证明：由周期序列的时域移位性质由于有限长序列的DFT就是周期序列DFS在频域中的主值序列，有（）频域移位定理若则上式称为频率移位定理频率移位定理，也称为调制定理调制定理，此定理说明时域序列的调制等效于频域的圆周移位。第23页，共48页，编辑于2022年，星期二、共轭对称性、共轭对称性任一序列都可以表示成共轭对称分量和共轭反对称分量之和。周期序列的共轭对称分量和共轭反对称分量都是周期性的，周期仍为N，取出它们的主值序列就得到了有限长序列的相应的分量，分别称为圆

11、周共轭对称分量和圆周共轭反对称分量，公式推导如下：设有限长序列x(n)的长度为N，以N为周期的周期延拓序列为第24页，共48页，编辑于2022年，星期二则有同样可以证明则有限长序列的圆周共轭对称分量和圆周共轭反对称分量定义为由于满足，有第25页，共48页，编辑于2022年，星期二DFTDFT的一系列的对称性质：的一系列的对称性质：（）式中x*(n)是x(n)的共轭复序列。（）（）复序列实部的DFT等于序列DFT的圆周共轭对称部分，即（）复序列虚部乘j的DFT等于序列DFT的圆周共轭反对称部分，即（）若x(n)是实序列，则X(k)只有圆周共轭对称部分，即满足（）若x(n)是纯虚数序列，则X(k)

12、只有圆周共轭反对称部分，即满足第26页，共48页，编辑于2022年，星期二（）例：设x1(n)和x2(n)都是实数序列，试求X1(k)和X2(k)解：先利用这两个实数序列构成复序列，有 又故同样故 因此可以用一次DFT计算出Y1(k)，然后用上面的公式计算出X1(k)和X2(k)。第27页，共48页，编辑于2022年，星期二例：试利用DFT的对称特性求和的DFT。解：设因为所以第28页，共48页，编辑于2022年，星期二而因为所以第29页，共48页，编辑于2022年，星期二、帕斯瓦尔（、帕斯瓦尔（ParsevalParseval）定理）定理证明：若y(n)=x(n),则即第30页，共48页，编

13、辑于2022年，星期二、圆周卷积、圆周卷积（）时域圆周卷积（）时域圆周卷积 设x1(n)和x2(n)都是N点的有限长序列，有若则此卷积过程与周期卷积和的过程是一致的，只不过这里要取结果的主值序列。公式中的只在0mN-1范围内取值，因而是圆周移位，因此这个卷积和称为圆周卷积和。圆周卷积和。第31页，共48页，编辑于2022年，星期二（）频域圆周卷积（）频域圆周卷积利用时域与频域的对称性，得到频域圆周卷积定理若则（）圆周相关定理（）圆周相关定理若则第32页，共48页，编辑于2022年，星期二（）用圆周卷积求线性卷积（）用圆周卷积求线性卷积 如果信号x(n)和单位抽样响应h(n)都是有限长序列，那么

14、是否能用圆周卷积的运算来代替线性卷积运算呢？下面就这个问题加以讨论：设x1(n)是N1点的有限长序列，x2(n)是N2点的有限长序列。x1(n)和x2(n)的线性卷积：x1(m)的非零区间为0mN1-1，x2(n-m)的非零区间为0n-mN1-1，将两个不等式相加，得到 0nN1+N2-2第33页，共48页，编辑于2022年，星期二 x1(n)和x2(n)的圆周卷积：假设x1(n)和x2(n)进行L圆周卷积,Lmax(N1,N2)，再讨论L等于何值时，圆周卷积才能代表线性卷积。将两个序列都补零为长度为L点的序列，即则第34页，共48页，编辑于2022年，星期二将任一序列（这里采用x2(n)）变

15、成L点周期延拓序列，即因此L L点的圆周卷积点的圆周卷积y(n)y(n)是线性卷积是线性卷积y yl l(n)(n)以以L L为周期的周期为周期的周期延拓序列的主值序列延拓序列的主值序列 结论：若若LNLN1 1+N+N2 2-1-1，则，则L L点圆周卷积能代表线性卷积。点圆周卷积能代表线性卷积。第35页，共48页，编辑于2022年，星期二图两个有限长序列的圆周卷积和线性卷积000nnn11234x1(n)y(n)=x1(n)x2(n)0n340n4y(n)=x1(n)x2(n)y(n)=x1(n)*x2(n)n134y(n)=x1(n)x2(n)0312(a)(b)(c)(d)(e)(f)

16、x2(n)第36页，共48页，编辑于2022年，星期二利用利用DFTDFT计算模拟信号的计算模拟信号的 傅立叶变换傅立叶变换(级数级数)对对对连续时间非周期信号的傅立叶对连续时间非周期信号的傅立叶 变换的变换的DFTDFT逼近逼近 连续时间非周期信号连续时间非周期信号x(t)x(t)的傅立叶变换对为的傅立叶变换对为用用DFTDFT方法计算这一变换对方法计算这一变换对:第37页，共48页，编辑于2022年，星期二、采样：对x(t)x(t)以T为间隔进行采样，即 由于 因此得到、截断：将序列x(nT)=x(n)截断成包含有N个抽样点的有限长序列，因此有由于时域抽样，抽样频率为fS=1/T，频域产生

17、以fS为周期的周期延拓，若频域为带限信号，则有可能不产生频域混迭，而成为连续周期频谱。、频域抽样：在频域的一个周期中取N个样点，每个样点间隔为F0，fSF0。频域抽样使频域的积分式变成求和式，而在时域就得到原来已经截断的离散时间序列的周期延拓，时域周期为F0。因此有 第38页，共48页，编辑于2022年，星期二得到一些参量关系：因此得到：第39页，共48页，编辑于2022年，星期二2 2对连续时间周期信号的对连续时间周期信号的 傅立叶级数的傅立叶级数的DFSDFS逼近逼近连续时间周期信号x(t)的傅立叶级数对为：T0为连续时间周期信号的周期.1、时域抽样x(n)=x(nT)=x(t)|t=nT

18、设一个周期内的样点数为第40页，共48页，编辑于2022年，星期二因此得到用DFS(DFT)来逼近连续时间周期信号傅立叶级数对的公式:、将频域离散序列截断，截断长度等于一个周期（时域抽样造成的频域周期延拓的一个周期），有第41页，共48页，编辑于2022年，星期二利用利用DFTDFT计算模拟信号时计算模拟信号时 可能出现的问题可能出现的问题、频域的混迭失真及参数的选择、频域的混迭失真及参数的选择、截断效应、截断效应、栅栏效应、栅栏效应第42页，共48页，编辑于2022年，星期二、频域的混迭失真及参数的选择、频域的混迭失真及参数的选择）根据采样定理，只有当采样频率fS大于信号的最高频率fh两倍时

19、，才能避免频域混迭。即fS2 fh。也就是抽样间隔为T满足=1/fS 1/2fh。实际信号的持续时间都是有限的，从理论上来说，其频谱宽度是无限的，在工程上总是对信号先进行低通滤波预滤波或抗混迭滤波，限制高于的频率分量出现。）DFT得到的频率函数也是离散的，其频域抽样间隔为F0，即频率分辨力,T0=1/F0为最短信号记录长度。为了对全部信号进行采样，必须使抽样点数N满足条件第43页，共48页，编辑于2022年，星期二例:有一频谱分析用的FFT处理器,其抽样点数必须是2的整数幂,假定没有采用任何特殊的数据处理措施,已知给定的条件为:频率分辨率10Hz,信号最高频率4kHz.试确定以下参量：最小记录

20、长度；抽样点间的最大时间间隔（最小抽样频率）；在一个记录中最少点数解：最小记录长度 抽样点间的最大时间间隔 第44页，共48页，编辑于2022年，星期二、截断效应、截断效应在实际中遇到的序列x(n)，其长度往往是很长，甚至是无限长的，用DFT对其进行谱分析时，必须将它截断为长度为N的有限长序列，即根据频率卷积定理，有式中，其中部分称为主瓣。假设，则第45页，共48页，编辑于2022年，星期二|RN(ej)|X(ej)|Y(ej)|2/N-2/N/4-/40(a)RN(ej)的幅频曲线(b)X(ej)的幅频曲线(c)Y(ej)的幅频曲线第46页，共48页，编辑于2022年，星期二序列截断后的频谱

21、与原序列频谱有着明显的差别，这种差别对谱分析带来两方面的影响：1 1）频谱泄露）频谱泄露原序列x(n)的频谱是离散谱线，经截断后使每根谱线都带上一个辛格谱，就好象使谱线向两边延伸，通常将这种因时域上的截断导致频谱展宽称之为“泄露”，显然泄露使频谱变得模糊，分辨率降低。2 2）谱间干扰）谱间干扰因截断使在主谱线两边形成许多旁瓣，引起不同分量间的干扰，称之为谱间干扰，这不仅影响频谱分辨率，严重时强信号的旁瓣可能湮灭弱信号的主谱线，或者将强信号谱的旁瓣误认为是另一信号的谱线，从而形成假信号，使谱分析产生较大的偏差。截断效应是无法完全消除的，只能根据要求折衷选择有关参量。首先可以取更长的数据，也就是使

22、截断窗加宽，当然数据太长也必然会导致存储量和运算量增加，其次数据不要突然截断，也就是不要加矩形窗，而是缓慢截断，即加各种缓变的窗（如三角窗、升余弦窗等），使得窗谱的旁瓣能量更小，卷积后造成的泄露减小。第47页，共48页，编辑于2022年，星期二、栅栏效应、栅栏效应N点DFT是在频率区间0，2上对信号的频谱进行N点等间隔采样，得到的是若干个离散点X(k)，且它们只限制为基频F0的整数倍，这就好象在栅栏的一边通过缝隙看另一边的景象，只能在离散点的地方看到真实的景象，其余部分频谱成分被遮挡，所以称为栅栏效应。减小栅栏效应，可以在时域数据末端增加一些零值点，使一个周期内的点数增加，但是不改变原有的记录数据，即增加频域抽样点数N，频域抽样为，这样必然使谱线更密谱线更密，这样原来看不到的谱分量就可能看到了。第48页，共48页，编辑于2022年，星期二