《《离散傅里叶变换》ppt课件》由会员分享,可在线阅读,更多相关《《离散傅里叶变换》ppt课件(82页珍藏版)》请在金锄头文库上搜索。
1、1,离散傅里叶变换,Discrete Fourier Transform,2,内容提要,离散傅里叶变换 (Discrete Fourier Transform,DFT)是时间函数是离散的,而且频谱函数也是离散的变换。 离散傅里叶变换定义 DFT物理意义 DFT基本性质 讨论频率取样理论。 DFT的应用,3,傅里叶变换的各种形式,连续时间、离散频率的傅里叶变换 对于周期为T的连续时间信号,可以采用傅里叶级数展开: 连续时间、连续频率的傅里叶变换 对于非周期的连续时间信号,可以进行傅里叶变换:,它在时域和频域都是连续的。,4,离散时间、连续频率的傅里叶变换 对于非周期的序列,其傅里叶变换在频域是以
2、2为周期的连续函数。,5,3.1 离散傅里叶变换(DFT) 设x(n)是一个长度为M的有限长序列, 则定义x(n)的N点离散傅里叶变换为,X(k)的离散傅里叶逆变换为,式中, , N称为DFT变换区间长度, NM, 通常称(3.1.1)式和(3.1.2)式为离散傅里叶变换对。 Note:有限长序列x(n)的DFT即X(k)仍是有限长序列。,6,例 3.1.1 x(n)=R4(n) ,求x(n)的8点和16点DFT. 解: 设变换区间N=8, 则,7,对长度为M的序列x(n),其Z变换,N点DFT,进行对比,可以看出,式中,,表示z平面单位圆上辐角,(k=0,1,N-1)的N个等间隔点。,3.1
3、.2 DFT与FT、Z变换的关系,8,说明:序列x(n)的N点DFT是其Z变换在单位圆上的N点等角距取样,如图3.4(a)。序列x(n)的DFT是其FT在区间0,2上的N点等间隔取样。如图3.4(b)。,9,3.2.3 DFT的隐含周期性,DFT变换对中,均为整数,所以式(3.1.1) 中, X(k)满足,同理可证明式(3.1.2) 中 x(n+mN)=x(n),10,任何周期为N的周期序列 都可以看做长度为N的有限长序列x(n)的周期延拓序列,而x(n)是 的一个周期,(3.1.5),(3.1.6),11,定义:,为叙述方便,将式(3.1.5)该写成,表示x(n)以N为周期的周期延拓序列,符
4、号(n)N表示n对模N的余数,即,这里k是商。,的主值区间:周期序列 中从n=0到N-1的范围,的主值序列:主值区间上的序列,12,由此对长度为N的序列x(n),且 ,则 的DFS为,结论:与DFT定义比较,可见有限长序列x(n)的DFT即X(k)是x(n)的周期延拓序列 的离散傅里叶级数系数 的主值序列。,例如,N=7, =x(n)7,则有,13,解:,因此得,X(0)=4.16114 X(1)=0.71063-j0.92558 X(2)=0.50746-j0.40597 X(3)=0.47017-j0.16987,X(4)=0.46235 X(5)= 0.47017+j0.16987 X(
5、6)= 0.50746+j0.40597 X(7)= 0.71063+j0.92558,Matlab实现 fft1.m,例3. 1 求有限长序列,的DFT,其中a=0.8,N=8。,14,关于离散傅里叶变换(DFT):,序列x(n)在时域是有限长的(长度为N),它的离散傅里叶变换X(k)也是离散、有限长的(长度也为N)。 n为时域变量,k为频域变量。 DFT的物理意义:序列x(n)的Z变换在单位圆上的等角距取样。序列傅里叶变换在区间0,2上的等间隔取样。 离散傅里叶变换(DFT)具有唯一性。 离散傅里叶变换与离散傅里叶级数没有本质区别,DFT实际上是离散傅里叶级数的主值,DFT也隐含有周期性。
6、,15,DFT隐含着周期性,因此在讨论DFT的性质时,常与DFS的概念联系起来,并把有限长序列看作周期序列的一个周期来处理。 设x1(n)和x2(n)的长度都为N,且它们对应的DFT分别为X1(k)和X2(k)。,1线性,设x3(n)=ax1(n)+bx2(n),a和b都为常数,则,若它们长度不等,取长度最大者,将短的序列通过补零加长,注意此时DFT与未补零的DFT不相等。,3.2 离散傅里叶变换的性质,16,2循环移位性质,一个长度为N的序列x(n)的循环移位定义为,循环移位分3步计算:,(1)将x(n)延拓成周期为N的周期序列 ; (2)将 移位得 或 x(n+m)N; (3)对 x(n+
7、m)N 取主值得 x(n+m)NRN(n)。 这个过程如下图所示。,a) 序列的循环移位:,17,从图中两虚线之间的主值序列的移位情况可以看出,当主值序列左移m个样本时,从右边会同时移进m个样本,而且好像是刚向左边移出的那些样本又从右边循环移了进来。因此取名“循环移位”。 显然,循环移位不同于线性移位,18,19,20,对长度为N的有限长序列x(n),其循环移位后序列y(n)的DFT为,证明:,b) 时域循环移位定理,令n+m=n,则有:,因为 以N为周期,上式中的求和区间任取一个周期即可,取主值区间为求和区间,得证。,21,若,则,c) 频域循环移位定理,22,3.2.3 循环卷积定理,长度
8、分别为N1和N2的有限长序列x1(n)和x2(n)的N点DFT分别为: ( N=max N1, N2 )。,则,由上式表示的卷积称为循环卷积,X1(k)=DFTx1(n) X2(k)=DFTx2(n),如果,X(k)=X1(k)X2(k),23,循环卷积的过程:,(1)周期延拓 x2(m)x2(m)N (2)折叠 x2(m)Nx2(-m)N (3) 移位和取主值 x2(-m)Nx2(n-m)NRN(m) (4)相乘 x2(n-m)NRN(m) x1(m) x2(n-m)NRN(m) (5)相加 summ0,1,N-1,循环反转序列,Note: 两个长度为N的序列的循环卷积长度仍为N,(与线性卷
9、积不同),记为:,24,25,26,循环卷积计算说明:x1(n)的N个值按顺时针方向均匀分布在内圆周上,x2(n)的N个值按反时针方向均匀分布在外圆周上,把内外圆周上对应的数值两两相乘,然后把乘积相加就得到y(0)。若将外圆周顺时针方向转动一格,将内外圆周上对应的数值两两相乘并把乘积相加,便得到y(1)。依次类推,可以得出y(n)的其它值。因此循环卷积也叫做圆卷积。,考虑到DFT关系的对偶性,可以证明,长为N的两序列之积的DFT等于它们的DFT的循环卷积除以N,即,频域循环卷积定理,27,28,3.2.4 复共轭序列的DFT,是长度为N的序列x(n)的复共轭序列,,则,且,类似,Note:对实
10、序列有,29,3.2.5 DFT的共轭对称性,1. 有限长共轭对称序列和共轭反对称序列,分别用xep(n)和xop(n) 表示有限长共轭对称序列和共轭反对称序列, 则二者满足如下定义式: xep(n)=x*ep(N-n), 0nN-1 (3.2.9) xop(n)=-x*op(N-n), 0nN-1 (3.2.10),当N为偶数时, 将上式中的n换成N/2-n可得到,30,图 3.2.3 共轭对称与共轭反对称序列示意图,(图中*表示对应点为序列取共轭后的值),31,如同任何实函数都可以分解成偶对称分量和奇对称分量一样, 任何有限长序列x(n)都可以表示成其共轭对称分量和共轭反对称分量之和, 即
11、 x(n)=xep(n)+xop(n), 0nN-1 (3.2.11) 将上式中的n换成N-n, 并取复共轭, 再将式(3.2.9) 和式(3.2.10) 代入得到 x*(N-n)=x*ep(N-n)+x*op(N-n)=xep(n)-xop(n) (3.2.12) xep(n)=1/2x(n)+x*(N-n) (3.2.13) xop(n)=1/2x(n)-x*(N-n) (3.2.14),32,2. DFT的共轭对称性,(a) 若将序列x(n)分成实部xr(n)与虚部xi(n),即 x(n)=xr(n)+jxi(n) 根据复共轭序列的DFT可得,再由DFT的线性性质可得,33,(b) 若将
12、序列x(n)分成共轭对称部分xep(n)与共轭反对称部分xop(n),即 x(n)=xep(n)+xop(n) 根据复共轭序列的DFT可得,因此,34,结论: 如果序列x(n)的DFT为X(k),则 x(n)的实部和虚部(包括j)的DFT分别为X(k)的共轭对称分量和共轭反对称分量; x(n)的共轭对称部分和共轭反对称部分的DFT分别为X(k)的实部和虚部乘以j.,35,有限长实序列DFT的共轭对称性: 对长度为N的实序列,X(k)=DFTx(n),则 X(k)共轭对称,即 若x(n)=x(N-n),则X(k)实偶对称,即 若x(n)=-x(N-n),则X(k)纯虚奇对称,即,36,这意味着,
13、对于时间有限信号,可以像频带有限信号进行时域采样而不丢失任何信息一样,可以在频域上进行采 样而不丢失任何信息。这正是傅里叶变换中时域和频域对偶关系的反映,这有着十分重要的意 义。DFT实现了频域离散化,开辟了在频域采用数字技术处理的新领域。 这使我们自然想到,对于任意一个频率特性,是否均能用频域采样的办法来逼近,这是一个很吸引人的问题,因为用频率采样来逼近,可使问题大大简化。因此我们要讨论频率采样的可行性以及所带来的误差。,3. 3 频率域采样,37,频率取样是指对序列的傅里叶变换或系统的频率特性进行取样 。,本节讨论在什么条件下能够用得到的频谱取样值无失真地恢复原信号或系统。,设任意长序列x
14、(n)绝对可和,其Z变换表示为,如果在单位圆上对X(z)进行等角距取样,取样点数为N,则得,根据DFT的定义,对X(k)求反变换,38,根据上面两式可得:,因为,所以,上式表明,在z平面的单位圆上对序列的Z变换进行等角距取样,将导致时间序列的周期延拓。这一结果与对连续时间信号取样导致频谱周期延拓类似。,现在我们来考察xp(n)与原序列x(n)的关系,看它如何才能代表原序列x(n)。,39,xp(n)是原非周期信号x(n)的周期延拓序列,因此xp(n)是一个周期序列,其主值为,在x(n)为有限长度M的情况下,如果取样点NM,那么x(n)周期延拓的结果不会产生混叠。这时,xp(n)的主值xN(n)
15、与原序列x(n)一样,因此xN(n)完全能代表原序列x(n),可由频域采样X(k)恢复x(n)。 如果NM,即延拓的周期N小余有限序列的长度M,则x(n)周期延拓后一定产生混叠,因而xN(n)不能无失真地代表原信号x(n)。在x(n)为无限长的情况下,对Z变换取样必然导致混叠失真,因此xN(n)不能代表原序列x(n)。,40,因此,对于长度为M的有限长序列,对Z变换取样即频率取样不失真的条件,是取样点数 N应等于或大于原序列的长度M,即NM。在NM时,Z变换的取样即DFT X(k),利用IDFT公式可由X(k)恢复原序列x(n),即,这就是频域采样定理。,41,对于有限长序列x(n),满足频域采样定理时,N点频域采样X(k)就足以不失真地代表序列的特性。因此,由此N个采样值X(k)应能完全地表达整个X(z)函数及频率特性X(ej)。即由N点X(k)可内插恢复出X(z)或X(ej)。,式中,所以,设序列x(n)长度为M,在频域02之间等间隔采样N点,N=M,则有,42,令,则,上式就是用X(z)在单位圆上的N个取样值X(k)表示X(z)的内插公式,内插函数为 。,因为,所以,4