DSP离散傅里叶变换.pptx

《DSP离散傅里叶变换.pptx》由会员分享，可在线阅读，更多相关《DSP离散傅里叶变换.pptx（70页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、3.1.1 DFT的定义设x(n)是一个长度为M的有限长序列，则定义x(n)的N点离散傅里叶变换为:X(k)的离散傅里叶逆变换(IDFT)为：3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义 旋转因子：N为变换区间的长度，NMIk=0第1页/共70页IDFTX(k)唯一性的证明由于：所以，在变换区间上满足下式：IDFTX(k)=x(n),0nN-1 离散傅里叶逆变换是唯一的。3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义M为整数 第2页/共70页例 序列x(n)=R4(n)，求x(n)的8点和16点DFT。解:(1)设变换区间N=8，则：(2)设变换区间N=16，则3.1 3

2、.1 离散傅里叶变换的定义离散傅里叶变换的定义3结论：离散傅立叶变换(DFT)结果与变换区间长度N有关。第3页/共70页3.1.2 DFT和Z变换的关系设序列x(n)的长度为N，其Z变换和DFT分别为：比较上面二式可得关系式3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义e第4页/共70页DFT的物理意义：(1)x(n)的N点DFT 是x(n)的Z变换在单位圆上N点等间隔采样。(2)X(k)是x(n)的傅里叶变换X(ejw)在区间0,2上的N点等间隔采样,采样间隔为2/N。(3)变换区间长度N不同，变换结果不同，N确定后，X(k)与x(n)是一一对应的。(4)当N足够大时，|X(k)|

3、的包络可逼近|X(ejw)|曲线；(5)|X(k)|表示wk=2k/N频点的幅度谱线。3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义第5页/共70页3.1.3 DFT的隐含周期性 在DFT变换的定义对中，x(n)与X(k)均为有限长序列。(1)旋转因子WknN的周期性(周期为N)(2)X(k)隐含的周期性(周期为N)(3)序列x(n)隐含的周期性(周期为N)3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义K,m,N均为整数x(n+mN)=x(n)第6页/共70页任何周期为N的周期序列 都可以看作长度为N的有限长序列x(n)的周期延拓序列，而x(n)则是 的一个周期，即:一般定

4、义周期序列 中从n=0到N-1的第一个周期为 的主值区间，而主值区间上的序列称为 的主值序列。总结：是x(n)的周期延拓序列 x(n)是 的主值序列3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义 0nN-1 0nN-1第7页/共70页为了以后叙述方便，可用如下形式表示：(n)N表示n对N求余，即如果n=MN+n1，0n1N-1，M为整数，则：(n)N=n1例：设N5，则有：3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义x(n)N表示：x(n)以N为周期的周期延拓序列。第8页/共70页DFT和周期序列的DFS的关系设x(n)的长度为N，且 ，则周期序列 的离散傅立叶级数表示式

5、：上式中：说明：有限长序列x(n)的离散傅立叶变换X(k)，正好是x(n)的周期延拓序列 的离散傅立叶级数系数 的主值序列3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义注意：是一周期序列 第9页/共70页总结DFTFTZT单位圆上的N点等间隔采样0,2上的N点等间隔采样单位圆上的Z变换,Z=ejw =DFS =DFSx(n)N =X(k)N第10页/共70页例1：若N=5，x(n)=R4(n)，画出x(n)N图形。3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义nx(n)10 1 2 3 4nx(n)510 1 2 3 45 6 7 8 9-3-2-4-5第11页/共70页例

6、2：已知长度为N的一个有限长序列x(n)，其N点DFT为X(k)。另一个长度为2N的序列 y(n)定义为：y(n)=x(0.5n)，n为偶数；0，n为奇数；试用X(k)表示y(n)的2N点离散傅立叶变换Y(k)。解：已知3.1 3.1 离散傅里叶变换的定义离散傅里叶变换的定义令：则：第12页/共70页3.2.1 线性性质如果x1(n)和x2(n)是两个有限长序列，长度分别为N1和N2 y(n)=ax1(n)+bx2(n)式中a、b为常数，取：N=maxN1,N2，则y(n)的N点DFT为 Y(k)=DFTy(n)=aX1(k)+bx2(k）,0kN-1 其中：X1(k)和X2(k)分别为x1(

7、n)和x2(n)的N点DFT。3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质第13页/共70页3.2.2 循环移位性质 1.序列的循环移位 设x(n)为有限长序列，长度为N，则x(n)的循环移位定义为 y(n)=x(n+m)NRN(N)(1)序列y(n)由x(n)以N为周期进行周期延拓而得到 (n)=x(n)N(2)再将 (n)左移m位，得到：(nm);(3)取 (nm)的主值区间得到有限长序列x(n)的循环移位y(n)3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质 0 N-1nn 0 N-1xxxx从左侧移出主值区的

8、序列值依次从右侧进入主值区从左侧移出主值区的序列值依次从右侧进入主值区n 0 N-1n 0 N-1x第14页/共70页2.时域循环移位定理 设 x(n)是 长 度 为 N的 有 限 长 序 列，y(n)为 x(n)的 循 环 移 位，即：y(n)=x(n+m)NRN(n)则:Y(k)=DFTy(n)=W-kmNX(k)其中:X(k)=DFTx(n),0kN-1 3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质提示：x(n)N和 均以N为周期，所以对其在任一周期上的求和结果相同令n+m=n，则有 证明：第15页/共70页3.频域循环移位定理如果：X(k)=DF

9、Tx(n),0kN-1 Y(k)=X(k+l)NRN(k)则 y(n)=IDFTY(k)=WnlNx(n)3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质3.2.3 循环卷积定理第16页/共70页1、时域循环卷积定理 有限长序列x1(n)和x2(n)，长度分别为N1和N2，N=max N1,N2。x1(n)和x2(n)的N点DFT分别为：X1(k)=DFTx1(n)X2(k)=DFTx2(n)如果：X(k)=X1(k)X2(k)3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质x2x1注意：对于x1(n)或x2(n)不足N点

10、，则分别在其尾部补零，使长度为N。则：x(n)第17页/共70页证明：直接对上式两边进行DFT令n-m=n，则有3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质x2第18页/共70页两个有限长序列循环卷积的过程：(1)上式中求和变量为m，n为参变量；(2)将x2(m)以N为周期作周期延拓得到x2(m)N；(3)翻转x2(m)N 形成x2(-m)N(4)对x2(-m)N进行循环移位x2(n-m)N，取主值序列，形成x2(n-m)N RN(m);(5)n=0,1,N-1时，x1(m)和x2(n-m)N R N(m)对应相乘，并对m在0N-1区间求和。3.2 3.

11、2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质x2第19页/共70页【例】：已知x1(n)=1，0n3；x2(n)=1，2n5；0，其它n;0，其它n;求y(n)=x1(n)x2(n)，循环卷积区间长度N为8。y(0)=x1(m)x2(-m)8R8(n)=1;y(1)=x1(m)x2(1-m)8R8(n)=0;y(2)=x1(m)x2(2-m)8R8(n)=1;y(3)=x1(m)x2(3-m)8R8(n)=2;y(4)=x1(m)x2(4-m)8R8(n)=3;y(5)=x1(m)x2(5-m)8R8(n)=4;Y(6)=x1(m)x2(6-m)8R8(n)=3;y

12、(7)=x1(m)x2(7-m)8R8(n)=2;3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质mx1(m)10 1 2 3 4mx2(m)11 204 53my(n)1014232 3 4 5 6 7第20页/共70页2、频域循环卷积定理如果：x(n)=x1(n)x2(n)则：3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质其中：X1(k)=DFTx1(n)X2(k)=DFTx2(n)0kN-1第21页/共70页证明：令：k-m=k，代入得到3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基

13、本性质第22页/共70页3.2.4 复共轭序列的DFT 设x*(n)是x(n)的复共轭序列，长度为N X(k)=DFTx(n)则:DFTx*(n)=X*(N-k),0kN-1 且:X(N)=X(0)证明：根据DFT的唯一性，只要证明上式右边等于左边即可。又由X(k)的隐含周期性有：X(N)=X(0)用同样的方法可以证明：DFTx*(N-n)=X*(k)3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质nn第23页/共70页3.2.5 DFT的共轭对称性序列的傅里叶的对称性是关于坐标原点的纵坐标的对称性，DFT的对称性关于N/2点的对称性。1、有限长共轭对称序列

14、和共轭反对称序列用xep(n)和xop(n)分别表示有限长共轭对称序列和共轭反对称序列，则二者满足如下定义式：xep(n)=x*ep(N-n),0nN-1 xop(n)=-x*op(N-n),0nN-1 当N为偶数时，将上式中的n换成N/2-n可得到 xep(N/2n)=x*ep(N/2n),0nN/2-1 xop(N/2n)=-x*op(N/2n),0nN/2-13.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质第24页/共70页3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质共轭对称序列示意图共轭反对称序列示意图第25页

15、/共70页2、任何一有限长序列都可表示成其共轭对称分量和共轭反对称分量之和 x(n)=xep(n)+xop(n),0nN-1 将上式中的n换成N-n，并取复共轭：x*(N-n)=x*ep(N-n)+x*op(N-n)=xep(n)-xop(n)由上两式可得：xep(n)=1/2x(n)+x*(N-n)xop(n)=1/2x(n)-x*(N-n)同理可以确定有限长序列X(k)的Xep(k)和Xop(k)Xep(k)=1/2X(k)+X*(N-k)；Xop(k)=1/2X(k)-X*(N-k)；3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质第26页/共70页3

16、、DFT的共轭对称性(1)如果x(n)=xr(n)+jxi(n)其中：xr(n)=Rex(n)=1/2x(n)+x*(n)jxi(n)=jImx(n)=1/2x(n)-x*(n)DFTxr(n)=1/2DFTx(n)+x*(n)=1/2X(k)+X*(N-k)=Xep(k)DFTjxi(n)=1/2DFTx(n)-x*(n)=1/2X(k)-X*(N-k)=Xop(k)3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质共轭对称分量共轭反对称分量第27页/共70页(2)如果x(n)=xep(n)+xop(n)，0nN-1 其中：xep(n)=1/2x(n)+x*

17、(N-n),x(n)的共轭对称分量xop(n)=1/2x(n)x*(N-n),x(n)的共轭反对称分量 那么：DFTxep(n)=1/2DFTx(n)+x*(N-n)=1/2X(k)+X*(k)=ReX(k)DFTxop(n)=1/2DFTx(n)-x*(N-n)=1/2X(k)-X*(k)=jImX(k)3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质虚部实部第28页/共70页3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质4、有限长实序列DFT的共轭对称性 设x(n)是长度为N的实序列，且X(k)=DFTx(n)，则

18、(1)X(k)共轭对称，即：X(k)=X*(N-k),0kN-1 (2)如果 x(n)=x(N-n)则：X(k)实偶对称，即：X(k)=X(N-k)(3)如果 x(n)=-x(N-n)则：X(k)纯虚奇对称，即：X(k)=-X(N-k)(4)有限长实序列DFT共轭对称性的应用 当N=偶数时，只需计算前N/2+1点的DFT;当N=奇数时，只需计算前(N+1)/2点的DFT;序列x(n)实偶对称序列x(n)实奇对称可减少运算量，提高运算效率第29页/共70页3.2 3.2 离散傅立叶变换（离散傅立叶变换（DFTDFT）的基本性质）的基本性质 通过计算一个N点DFT，可得到两个不同实序列的N点DFT

19、。设：x1(n)和 x2(n)为 两 个 实 序 列，构成新序列x(n)如 下 ：x(n)=x1(n)+jx2(n)对x(n)进行DFT 得到:X(k)=DFTx(n)=Xep(k)+Xop(k)Xep(k)=DFTx1(n)=1/2X(k)+X*(N-k)Xop(k)=DFTjx2(n)=1/2X(k)-X*(N-k)所以:X1(k)=DFTx1(n)=1/2X(k)+X*(N-k)X2(k)=DFTx2(n)=j1/2X(k)-X*(N-k)第30页/共70页3.3 3.3 频率域采样频率域采样时域采样定理 在一定条件下，时域离散采样信号可以恢复出原来的连续信号；问题 在频域进行离散采样，

20、得到的离散采样值能否恢复出原来的信号(或原频域函数)。条件是什么？内插公式？第31页/共70页3.3 3.3 频率域采样频率域采样 设任意序列x(n)的Z变换为：设：X(z)收敛域包含单位圆(即x(n)存在FT)。在单位圆上对X(z)等N点间隔采样，得到：序列x(n)的FT在区间0,2上的N点等间隔采样k=0k=2k=1k=3k=N-1设离散序列x(k)是长度为N的有限长序列xN(n)的DFT，即问题：xN(n)与原序列x(n)之间是怎样的关系？xN(n)=IDFTX(k)，0nN-1第32页/共70页3.3 3.3 频率域采样频率域采样DFT与DFS的关系:X(k)是xN(n)以N为周期的周

21、期延拓序列 的离散傅里叶级数系数 的主值序列，即:X(n)X(k)N=DFS X(n)X(n)为整数 其它m 因为：第33页/共70页3.3 3.3 频率域采样频率域采样由上面推导可得：结论：X(z)在单位圆上的N点等间隔采样X(k)的IDFT，为原序列x(n)以N为周期的周期延拓序列的主值序列。频域采样定理：假设 x(n)的长度为M，频域采样点数为N 若 N M，则xN(n)=IDFTX(k)=x(n)时域无混叠 若 NM，则xN(n)=IDFTX(k)x(n)产生时域混叠 故频率抽样(不失真)条件为:N Mr=x(n)N RN(n)第34页/共70页3.3 3.3 频率域采样频率域采样例：

22、已知x(n)=R4(n)，X(ejw)=FTx(n)，对X(ejw)在区间0,2进行6点的等间隔采样，求：X6(k)，k=0,1,.5 及相应的x6(n)=IDFTX6(k)，n=0,1,.5。解：0 3n1 k=0k=2k=1k=3k=4k=50 1 2 3 4 5 4.001.731.00K=X(k)=0-j1.734.0001.00j1.731.0013245n=0第35页/共70页3.3 3.3 频率域采样频率域采样直接由频域采样定理得:2,对X(ejw)在一个周期内进行3点采样,求 及相应的x3(n)=IDFTX3(k)，n=0、1、2。（时域无混叠）n 01 6-6 -339第36

23、页/共70页3.3 3.3 频率域采样频率域采样解：直接由频域采样定义得时域混叠 02-6n 1369 21-3第37页/共70页用频域采样用频域采样X(k)X(k)表示表示X(z)X(z)的内插公式和内插函数的内插公式和内插函数设序列x(n)长度为M，在频域02之间等间隔采样N点，NM，则有:得到N个采样点代入X(z)的表达式令：则：内插函数内插公式第38页/共70页3.3 3.3 频率域采样频率域采样当z=ej时，上面两式成为x(n)的傅里叶变换X(ej)的内插函数和内插公式，即：进一步化简可得：X(ejw)在每个采样点上的函数值等于原始采样点值X(k),而采样点间的函数值是由N个内插函数

24、 按采样值X(k)的加权线性组合。频率采样内插（恢复）第39页/共70页3.4 DFT3.4 DFT的应用举例的应用举例 DFT的快速算法FFT的出现，使DFT在数字通信、语言信号处理、图像处理、功率谱估计、仿真、系统分析、雷达理论、光学、医学、地震以及数值分析等各个领域都得到广泛应用。这些应用一般是以卷积和相关运算的具体处理为依据，或者以DFT作为连续傅里叶变换的近似为基础。本节主要内容 (1)用DFT计算卷积 (2)用DFT对连续信号和序列进行谱分析 第40页/共70页3.4.1 3.4.1 用用DFTDFT计算线性卷积计算线性卷积1.用DFT计算循环卷积已知：X1(k)=DFTx1(n)

25、，X2(k)=DFTx2(n)求：L点x1(n)x2(n)=?时域直接卷积法:频域间接法计算：由时域循环卷积定理知：Y(k)=X1(k)X2(k)对Y(k)进行L点IDFT得y(n)，即 y(n)=IDFTY(k)DFTDFTIDFTL点L点M点x1(n)N点x2(n)X2(k)X1(k)Lx2(n)x1(n)第41页/共70页3.4 DFT3.4 DFT的应用举例的应用举例2、线性卷积和循环卷积的关系设：x(n)和h(n)的长度分别为M和N两序列的线性卷积：两序列的循环卷积：yl(n)长度为N+M-1 LmaxN,M第42页/共70页3.4 DFT3.4 DFT的应用举例的应用举例结论：yc

26、(n)等于yl(n)以L为周期的周期延拓序列的主值序列。只要保证LN+M-1，即循环卷积长度L大于等于线性卷积长度(N+M-1)，yc(n)=yl(n)。第43页/共70页3.4 DFT3.4 DFT的应用举例的应用举例例：已知序列x(n)和h(n)如图所示，求 (1)y1(n)=x(n)*h(n)；(2)y2(n)=x(n)h(n)；(3)y3(n)=x(n)h(n)；(4)y4(n)=x(n)h(n)。解：(1)y1(n)=x(n)*h(n)=线性卷积长度为 N+M-1=4+5-1=8点长。根据yc(n)=yl(n)LRL(n)，可计算出各循环卷积值。(2)y2(n)=x(n)h(n)=x

27、(m)h(n-m)6 R6(n)=3,3,3,4,4,3;nh(n)0 1 2 3 41nx(n)01 2 3 4 515m=0第44页/共70页3.4 DFT3.4 DFT的应用举例的应用举例(3)y3(n)=x(n)h(n)=x(m)h(n-m)8=1,2,3,4,4,3,2,1;(4)y3(n)=x(n)h(n)=x(m)h(n-m)10=1,2,3,4,4,3,2,1,0,0;7m=09m=0nyl(n)13428点长234101 23 45 67ny2(n)343长6点4发生混迭3 301 23 45没有混迭y3(n)1342长8点2301 23 45 6714nny4(n)1342

28、长10点23没有混迭，补04101 23 45 67 8 9第45页/共70页3.3.用用DFTDFT计算线性卷积计算线性卷积(快速卷积）快速卷积）设:x(n)是长度为M点序列，h(n)为N点序列。线性卷积：yl(n)=x(n)h(n)，序列yl(n)长:N+M-1；循环卷积：yc(n)=x(n)h(n)，序列yc(n)长L；当 LN+M-1时，用DFT计算线性卷积yl(n)步骤如下：计算x(n)和h(n)长度为L的DFT:X(k)=DFTx(n)，H(k)=DFTh(n)计算 Yc(k)=X(k)H(k)，根据时域循环卷积定理；yl(n)=yc(n)=IDFTYc(k)图：用DFT计算线性卷

29、积 yl(n)框图M点x(n)N点h(n)yl(n)补LN个零点补LM个零点L点DFTL点DFTL点IDFTX(k)H(k)X(k)H(k)第46页/共70页4 4、重迭相加法计算无限长线性卷积、重迭相加法计算无限长线性卷积快速线性卷积法针对的是：两个长度相差并不大序列。对两个长度相差很大的序列，比如：x(n)序列长度很长，h(n)序列长度比较短，如何有效计算x(n)h(n)=？两种计算方法：重迭相加法和重迭保留法重迭相加法计算步骤：设h(n)为N点长，x(n)为无限长(1)将x(n)均匀分段，每段长度为M，用xk(n)表示第K段序列。即:xk(n)=x(n)RM(nkM)x(n),kM n

30、(K+1)M1；0 ,其它n;xk(n)=第47页/共70页3.4 DFT3.4 DFT的应用举例的应用举例(2)分段线性卷积：=每段长L=M+N1(3)对yk(n)后N-1个点和yk1(n)前N-1个点的幅度值重叠相加。第48页/共70页3.4.2 3.4.2 用用DFTDFT对信号进行谱分析对信号进行谱分析信号的谱分析：就是计算信号的傅里叶变换。连续信号与系统的傅里叶分析不便于直接用计算机进行计算，应用受到限制。DFT是一种时域和频域均离散化的变换，适合数值运算，成为分析离散信号和系统的有力工具。1.用DFT对连续信号进行谱分析工程中经遇到的连续信号xa(t)，其频谱函数Xa(j)也是连续

31、函数。先对xa(t)进行时域采样，得到时域离散信号x(n)=xa(nT)；对x(n)进行DFT，得到的X(k)是x(n)的傅里叶变换X(ejw)在区间0,2上的N点等间隔采样；x(n)和X(k)均是有限长序列；第49页/共70页 DFTDFT对对x xa a(t)(t)进行频谱分析进行频谱分析傅里叶变换理论 信号持续时间有限长，其频谱是无限宽。信号的频谱有限长，在时域中，该信号的持续时间无限长。上述两种情况，在时域或频域中进行采样，得到的序列都是无限长序列，不满足DFT的变换条件。采用的处理方法：在频域中用滤波器滤除高于折叠频率的高频分量，在时域中则是截取有限点进行DFT。结论：用DFT对连续

32、信号进行谱分析是一种近似的分析，近似程度与信号带宽、采样频率和截取的长度有关。第50页/共70页3.4 DFT3.4 DFT的应用举例的应用举例设连续信号xa(t)持续时间为Tp，最高频率为fc，如下图(a)所示。则xa(t)的傅里叶变换为：Tp第51页/共70页3.4 DFT3.4 DFT的应用举例的应用举例对xa(t)以采样频率fs=1/T2fc进行采样得：x(n)=Xa(nT)。设共采样N点，并对Xa(jf)作零阶近似(t=nT,dt=T)得:对 x(jf)在区间0,fs上等间隔采样N点，采样间隔为F，参数fs、Tp、N和F满足如下关系式：令f=KF，频域N点采样得:令X(jkF)=Xa

33、(k)，xa(nT)=x(n)，代入得函数值与区间长度T的乘积和F=fs/N=1/NT=1/Tp,FT=1/N第52页/共70页3.4 DFT3.4 DFT的应用举例的应用举例结论:(1)连续信号的频谱特性可以通过对连续信号采样，并进行DFT再乘以T的近似方法得到。(2)连续信号的时域采样信号可以通过对其频谱函数进行采样，并进行IDFT再乘以1/T的近似方法得到。误差现象：(1)分析的结果看不到xa(jf)的全部特性，只能看到N个离散采样点的谱特性，这就是栅栏效应。(2)如果持续时间无限长，分析时要进行截断处理，这样会产生频谱混叠和泄漏现象，使谱分析产生误差。第53页/共70页3.4 DFT3

34、.4 DFT的应用举例的应用举例【例】理想低通滤波器的单位冲激响应ha(t)及其频响函数Ha(if)如图所示。用DFT来分析ha(t)的频率响应特性。由于ha(t)的持续时间为无穷长，所以要截取一段Tp，假设Tp=8 s，采样间隔T=0.25 s，采样点数N=Tp/T=32。频域采样间隔F=1/NT=0.125 Hz。则H(k)=TDFTh(n),0k31，其中：h(n)=ha(nT)R32(n)整个频响有波动，高频部分误差较大第54页/共70页3.4 DFT3.4 DFT的应用举例的应用举例对连续信号进行谱分析主要关心的两个问题：谱分析的范围fc：受采样频率fs的限制，fc 2 fc 谱分辩

35、率:F=fs/N采样点数N的选择：N 2fc/F信号观察时间Tp的选择：Tp 1/F提高F:(1)如保持N不变，必须fs 降低，导致谱分析范围减小；(2)fs 不变，增加采样点数N，即增加Tp第55页/共70页例：对实信号进行谱分析，要求谱分辨率F10 Hz，信号最高频率fc=2.5kHz，试确定最小记录时间TPmin，最大的采样间隔Tmax，最少的采样点数Nmin。如果fc不变，要求谱分辨率增加一倍，最少的采样点数和最小的记录时间是多少?解：根据信号观察时间TP的选择原则：TP 1/F=1/10=0.1s 因为要求:fs2fc，最小的采样频率为2fc，所以:频率分辨率提高一倍，即：F=5 H

36、z TPmin=1/F=1/5=0.2sTmax=1/2fc=Nmin=2fc/FNmin=2fc/F观察时间增加一倍，采样点数增加了一倍第56页/共70页2.2.用用DFTDFT对序列进行谱分析对序列进行谱分析 单位圆上的Z变换就是序列傅里叶变换。X(ejw)是w的连续周期函数，对序列x(n)进行N点DFT，得到X(k)，X(k)是在区间0,2上的N点等间隔采样。序列x(n)的傅里叶变换可利用DFT来计算。第57页/共70页(2)(2)对周期序列的频谱分析对周期序列的频谱分析设序列 (n)=x(n+rN)是周期为N的周期序列，则其傅立叶变换为：周期序列的频谱结构可以用离散傅里叶级数系数 表示

37、 取 的主值序列 进行N点DFT，得到周期序列的频谱结构也可以用其主值序列的离散傅里叶变换X(k)来表示(分析)xx(n)第58页/共70页令：n=n+rN，r=0,1,m-1，n=0,1N-1,则 x xM(n)=(n)RM(n)即：M=mN，m为整数截取序列的长度M为 (n)的整数个周期x设：n=n+rN第59页/共70页3.4 DFT3.4 DFT的应用举例的应用举例周期序列的频谱结构也可以用xM(k)表示 分析:(1)只有在k=rm时，XM(rm)=m ，表示 (n)的r次谐波谱线，幅度扩大了m倍，在其它k值，XM(k)0。(2)X(r)与XM(rm)对应点的频率相等。(3)只要截取

38、(n)整数个周期进行DFT，就可得到它的频谱结构，达到谱分析的目的。k/m=整数k/m整数 X(r)xx第60页/共70页3.4 DFT3.4 DFT的应用举例的应用举例若事先不知道x(n)的周期，怎样进行频谱分析：先截取x(n)M点，则求xM(n)的DFT：xM(n)=x(n)RM(n)，XM(k)=DFTxM(n)，0k M-1;再截取x(n)2M点，则求x2M(n)的DFT：x2M(n)=x(n)R2M(n)，X2M(k)=DFTx2M(n)，0k 2M-1;将2次截取序列的频谱进行分析，是否满足误差要求，若不满足，应加大截取窗长度（增加M值），再将结果进行分析。第61页/共70页3.4

39、 DFT3.4 DFT的应用举例的应用举例16点DFT相当于在序列后补零第62页/共70页3.4 DFT3.4 DFT的应用举例的应用举例x(n)=cos n/416点相当于取周期序列的两个周期进行DFT第63页/共70页3.4 DFT3.4 DFT的应用举例的应用举例3用DFT进行谱分析的误差问题 DFT(实际中用FFT计算)可用来对连续信号和数字信号进行谱分析。在实际分析过程中，要对连续信号采样和截断，由此可能产生误差分析。(1)频谱混迭现象：原因:不满足时域采样定理 避免措施：采样频率fs2fc，以避免信号在w=处附近的混迭。具体方法是：采样时满足采样定理，采样前对信号进行预滤波，滤去信

40、号中频率高于fs/2的频率分量。第64页/共70页3.4 DFT3.4 DFT的应用举例的应用举例(2)栅栏效应：现象：N点DFT是在区间0,2上的N点等间隔采样，采样点之间的频谱函数值是不知道的，就好像从N+1个栅栏缝隙中观看信号的频谱特性，得到的是N个缝隙中看到的频谱函数值，这种现象称为栅栏效应。原因：对信号的频谱进行有限点采样。后果：栅栏效应可能漏掉(挡住)大的频谱分量 减少栅栏效应的措施：对原序列补0，增大N，以增加采样点；第65页/共70页3.4 DFT3.4 DFT的应用举例的应用举例(3)截断效应：原因:对序列x(n)截断所引起的。无限长序列x(n)截短成有限长序列y(n)，即

41、y(n)=x(n)RN(n)，则 Y(ejw)=FTy(n)=1/(2)X(ejw)*RN(ejw)=1/(2)X(ej)*RN(ej(w-)d，其中X(ejw)=FTx(n)RN(ejw)=FTRN(n)=e-jw(N-1)/2sin(wN/2)sin(w/2)=RN(w)ej(w)RN(w)w02NN2N矩形窗函数幅度谱主瓣旁瓣第66页/共70页3.4 DFT3.4 DFT的应用举例的应用举例例：x(n)=cos(w0n)，w0=/4，用DFT分析其频谱特性。解：序列的幅度谱X(ejw)=(w-/4-2l)+(w+/4-2l)加矩形窗截断后,Y(ejw)=1/2X(ejw)*RN(ejw)

42、，定性图如下可见，截断后的频谱Y(ejw)与原序列频谱X(ejw)存在差别表现为频谱泄漏：在上图中，原谱线是离散谱线，而截短后，原来的离散谱线向附近展宽，常称这种展宽为泄漏。使谱分辨率F降低。泄漏原因是截取的窗函数有限长。l=-w0-4-4Y(ejw)N/2加矩形窗后幅度谱 w0-4-4X(ejw)x(n)=cos(w0n)的频谱RN(w)w02NN2N矩形窗函数幅度频第67页/共70页3.4 DFT3.4 DFT的应用举例的应用举例谱间干扰：在主谱线两边形成很多旁瓣，引起不同频率分量间的干扰（简称谱间干扰），影响频谱分辨率F，旁瓣的信号很强时，可能湮没弱信号的主谱线，导致较大的偏差。上述两种现象都是由于截短序列引起的，统称截断效应。为了减小截短效应的影响，可采取以下措施。窗函数不变，增大采样点N值：使主瓣变窄（4/N），提高频率分辨率。但旁瓣个数，相对幅度大小不变，即谱间干扰不变。采样点N不变，改变窗函数：选用旁瓣小的窗函数，使旁瓣个数减少，相对幅度减小，谱间干扰减小。但旁瓣越小，其主瓣就越宽，从而使谱分辨率降低。谱分辨率与谱间干扰是一对矛盾体，要综合考虑和兼顾。第68页/共70页第三章 作业本章第一次作业习题 1（5）（7）习题 4习题 7本章第二次作业习题 14习题 15离散傅里叶变换(DFTDFT)第69页/共70页感谢您的欣赏！第70页/共70页