数字信号处理复习总结-最终版_第1页
数字信号处理复习总结-最终版_第2页
数字信号处理复习总结-最终版_第3页
数字信号处理复习总结-最终版_第4页
数字信号处理复习总结-最终版_第5页
已阅读5页,还剩45页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

绪论:本章介绍数字信号处理课程的基本概念。0.1信号、系统与信号处理1.信号及其分类信号是信息的载体,以某种函数的形式传递信息。这个函数可以是时间域、频率域或其它域,但最基础的域是时域。分类:周期信号/非周期信号确定信号/随机信号能量信号/功率信号连续时间信号/离散时间信号/数字信号按自变量与函数值的取值形式不同分类:2.系统系统定义为处理(或变换)信号的物理设备,或者说,凡是能将信号加以变换以达到人们要求的各种设备都称为系统。3.信号处理信号处理即是用系统对信号进行某种加工。包括:滤波、分析、变换、综合、压缩、估计、识别等等。所谓“数字信号处理”,就是用数值计算的方法,完成对信号的处理。0.2数字信号处理系统的基本组成数字信号处理就是用数值计算的方法对信号进行变换和处理。不仅应用于数字化信号的处理,而且第1页共60页也可应用于模拟信号的处理。以下讨论模拟信号数字化处理系统框图。(1)前置滤波器将输入信号xa(t)中高于某一频率(称折叠频率,等于抽样频率的一半)的分量加以滤除。(2)A/D变换器在A/D变换器中每隔T秒(抽样周期)取出一次xa(t)的幅度,抽样后的信号称为离散信号。在A/D变换器中的保持电路中进一步变换为若干位码。(3)数字信号处理器(DSP)(4)D/A变换器按照预定要求,在处理器中将信号序列x(n)进行加工处理得到输出信号y(n)。由一个二进制码流产生一个阶梯波形,是形成模拟信号的第一步。5)模拟滤波器把阶梯波形平滑成预期的模拟信号;以滤除掉不需要的高频分量,生成所需的模拟信号ya(t)。0.3数字信号处理的特点1)灵活性。(2)高精度和高稳定性。(3)便于大规模集成。(4)对数字信号可以存储、运算、系统可以获得高性能指标。0.4数字信号处理基本学科分支数字信号处理(DSP)一般有两层含义,一层是广义的理解,为数字信号处理技术——DigitalSignalProcessing,另一层是狭义的理解,为数字信号处理器——DigitalSignalProcessor。0.5课程内容该课程在本科阶段主要介绍以傅里叶变换为基础的“经典”处理方法,包括:(1)离散傅里叶变换及其快速算法。(2)滤波理论(线性时不变离散时间系统,用于分离相加性组合的信号,要求信号频谱占据不同的频段)。在研究生阶段相应课程为“现代信号处理”(AdvancedSignalProcessing)。信号对象主要是随机信号,主要内容是自适应滤波(用于分离相加性组合的信号,但频谱占据同一频段)和现代谱估计。简答题:1.按自变量与函数值的取值形式是否连续信号可以分成哪四种类型?2.相对模拟信号处理,数字信号处理主要有哪些优点?3.数字信号处理系统的基本组成有哪些?第2页共60页第一章:本章概念较多,需要理解和识记的内容较多,学习时要注意。1.1离散时间信号1.离散时间信号的定义离散时间信号是指一个实数或复数的数字序列,它是整数自变量n的函数,表示为x(n)。一般由模拟信号等间隔采样得到:x(n)xatnTxa(nT)n。时域离散信号有三种表示方法:1)用集合符号表示2)用公式表示3)用图形表示2.几种基本离散时间信号(记住定义)1)单位采样序列2)单位阶跃序列第3页共60页3)矩形序列4)实指数序列5)正弦序列ω是正弦序列数字域的频率,单位是弧度。对连续信号中的正弦信号进行采样,可得正弦序列。设连续信号为,它的采样值为,因此(重点)这个式子具有一般性,它反映了由连续信号采样得到的离散序列,其数字频率与模拟频率的一般关系。另外需要说明的是,ω的单位为弧度,Ω的单位为弧度/秒。本书中,我们一律以ω表示数字域频率,而以Ω及f表示模拟域频率。例:已知采样频率FT=1000Hz,则序列x(n)=cos(0.4πn)对应的模拟频率为(400π)弧度/s。说明:本题旨在理解数字频率与模拟频率之间的关系:。FT(6)复指数序列复指数序列是以余弦序列为实部、正弦序列为虚部所构成的一个复数序列。(7)周期序列(重点)所有n存在一个最小的正整数N,满足:x(n)x(nN),则称序列x(n)是周期序列,周期为N。(注意:按此定义,模拟信号是周期信号,采用后的离散信号未必是周期的)例:正弦序列sin(0n)的周期性:2k当0N2k,k为整数时,sin[N0(nN)]sin(0n),即为周期性序列。周期0,式中,k、N限取整数,且k的取值要保证N是最小的正整数。可分几种情况讨论如下:(1)当2/0为整数时,只要k1,N2/0就为最小正整数,即周期为2/0。(2)当2/0不是整数,而是一个有理数时,设2/0P/Q,式中,P、Q是互为素数的整数(互为素数就是两个数没有公约数),取kQ,则NP,即周期为P。(3)当2/0是无理数时,则任何k皆不能使N为正整数,这时,正弦序列不是周期性的。例:X(n)=cos(0.4πn)的基本周期为(5)。[说明]基本周期的定义即计算公式:N2k,其中N和k均为整数,N为基本周期(使得N为最小整数时k取值)。本题ω=0.4π,代入上式得到:N5,k1。第4页共60页3.信号运算(1)加法:两个信号之和由同序号的序列值逐点对应相加得到。(2)乘法:两个信号之积由同序号的序列值逐点对应相乘得到。(3)移位:当,序列右移(称为延时);当,序列左移(称为超前)。(4)翻转:(5)尺度变换:或,其中M和N都是正整数。当时,序列是通过取x(n)的每第M个采样形成,这种运算称为下采样。对于序列,定义如下这种运算称为上采样。4.信号分解(重点)任一信号x(n)可表示成单位脉冲序列的移位加权和:简记为1.2时域离散系统时域离散系统定义x(n)T.y(n)y(n)Tx(n)线性系统(重点)判定公式:若y1(n)=T[x1(n)],y2(n)=T[x2(n)]则y(n)T[ax1(n)bx2(n)]ay1(n)by2(n)第5页共60页时不变系统(重点)判定公式:y(n)=T[x(n)]y(n-n0)=T[x(n-n0)]例:判断下列系统是否为线性、时不变系统。(重点)(1)y(n)x(n)2x(n1)3x(n2);(2)y(n)x2(n);解:y'(n)x(nn)2x(nn1)3x(nn2)(1)令:输入为x(nn0),输出为000n)x(nn)2x(nn1)3x(nn2)y'(n)y(n0000故该系统是时不变系统。y(n)T[ax1(n)bx2(n)]ax1(n)bx2(n)2(ax1(n1)bx2(n1))3(ax1(n2)bx2(n2))T[ax1(n)]ax1(n)2ax1(n1)3ax1(n2)T[bx2(n)]bx2(n)2bx2(n1)3bx2(n2)T[ax1(n)bx2(n)]aT[x1(n)]bT[x2(n)]故该系统是线性系统。(2)y(n)x2(n)令:输入为x(nn0),输出为y'(n)x2(nn0),因为y(nn0)x2(nn0)y'(n)故系统是时不变系统。又因为T[ax1(n)bx2(n)](ax1(n)bx2(n))2aT[x1(n)]bT[x2(n)]ax12(n)bx22(n)因此系统是非线性系统。线性时不变系统(LTI或者LSI系统)输入与输出之间关系(重点):h(n)T[(n)]第6页共60页y(n)x(m)(nm)my(n)T[x(m)(nm)]my(n)=x(m)h(nm)=x(n)*h(n)m重点:线性离不变系统的输出等于输入序列和该系统的单位脉冲响应的卷积【说明】离散时间LTI系统的单位冲激响应h(n)为系统对单位冲激序列δ(n)的零状态响应。单位冲激响应的概念非常重要。在时域,LTI系统可以由其单位冲激响应h(n)唯一确定,因此,我们常常用单位冲激响应描述LTI系统。在这种情况下,LTI系统的输入输出关系可以由卷积运算描述:y(n)=m

x(m)h(nm)=x(n)*h(n)物理意义:卷积和运算具有显式意义,即可以用来确定系统的输出。如果系统确定,则其单位冲激响应是唯一的。由此,可求系统对任意输入的响应。注意:计算卷积和的关键是求和区间的确定。因此,常常需要绘制序列x(m)和h(n-m)的图形。利用序列x(m)和h(n-m)的图形可助我们方便地确定求和区间。卷积的求解方法(重点):线性卷积是一种非常重要的一种运算,对它的求解,一般我们采用作图法。线性卷积满足交换律,设两序列长度分别是N和M,线性卷积后序列的长度为N+M-1。卷积的计算过程包括翻转、移位、相乘、相加四个过程。1)将和用和表示,画出和这两个序列;2)选择一个序列,并将其按时间翻转形成序列;3)将移位n,得到;4)将和相同m的序列值对应相乘后,再相加。例:设x(n)n,0≤n≤4,h(n)4,x(n)和h(n)如图1所示。求x(n)和h(n)的卷积y(n)。(重R(n)点)x(n)R4(n)4n1n012340123图1解方法一:用图解法求卷积和。(1)将x(n)和h(n)用x(m)和h(m)表示(图2中(a)、(b)图)。第7页共60页x(m)4R4(m)R4(m)mmm012340123-3-2-10(a)(b)(c)R4(1m)R4(2m)-2-101mmy(n)-101210(d)(e)R4(5m)mn01234501234567(f)(g)图2图解法求卷积过程(2)将h(m)进行反折,形成h(m)(图2中(c)图);将h(m)移位n,得到h(nm)(图2中(d)、(e)、(f)图)。(3)将x(m)和h(nm)相同m的序列值相乘,再相加,得到y(n)(图2中(g)图)。y(n)1,3,6,10,9,7,41≤n≤7再讨论解析法求线性卷积。用式y(n)x(m)h(nm)m求解上式首先要根据x(m)和h(nm)的非零值区间确定求和的上下限,x(m)的非零值区间为1≤m≤4,h(nm)的非零值区间为0≤nm≤3,或n3≤m≤n,由两个非零值区间可得n的取值区间为1≤n≤7,它们的乘积x(m)h(nm)的非零值区间应满足:1≤m≤4和n3≤m≤n因此当n1、n7时,y(n)0;nn(n1)y(n)m12当1≤n≤3时,m0;4(n1)(8n)y(n)m1当4≤n≤7时,2。mn3与图解法结果一致。y(n)用公式表示为第8页共60页n(n1)/21≤n≤3y(n)(n1)(8n)/24≤n≤70其他方法二:当序列x(n)和h(n)的长度分别为有限长N和M时,可采用“不进位乘法”求两序列线卷积。x(n)0,1,2,3,4h(n)1,1,1,1如图1所示:,y(n)0,1,3,6,10,9,7,4例:两线性时不变系统级联,其单位取样响应分别为h1(n)和h2(n),输入为x(n),求系统的输出y(n)。已知:x(n)u(n),h1(n)(n)(n4),h2(n)anu(n)。解:设第一个系统的输出为(n),则(n)x(n)h1(n)u(n)[(n)(n4)]u(n)u(n4)+(n1)+(n+(n)2)(n3)因而输出为y(n)(n)h2(n)[(n)(n1)(n2)(n3)]anu(n)anu(n)an1u(n1)an2u(n2)an3u(n3)系统因果性和稳定性的判定(重点)1)稳定系统:有界的输入产生的输出也有界的系统,即:若|x(n)|,则|y(n)|(记住!!)线性移不变系统是稳定系统的充要条件:n

|h(n)|(系统稳定的充分必要条件是系统的单位脉冲响应绝对可和)(记住!!)或:其系统函数H(z)的收敛域包含单位圆|z|=1(记住!!)第9页共60页2)因果系统:n0时刻的输出y(n0)只由n0时刻之前的输入x(n),nn0决定(记住!!)线性移不变系统是因果系统的充要条件:h(n)0,n0(记住!!)因果系统的单位脉冲响应必然是因果序列。(记住!!)或:其系统函数H(z)的收敛域在某圆外部:即:|z|>Rx(记住!!)3)稳定因果系统:同时满足上述两个条件的系统。线性移不变系统是因果稳定系统的充要条件:n

|h(n)|,h(n)0,n0(记住!!)或:H(z)的极点在单位圆内H(z)的收敛域满足:|z|Rx,Rx1(记住!!)例:判断线性时不变系统的因果性、稳定性,并给出依据。(重点)(1)y(n)1N1x(nk);Nk0nn0(2)y(n)x(k);knn0解:(1)只要N1,该系统就是因果系统,因为输出只与n时刻的和n时刻以前的输入有关。如果x(n)M,则y(n)M,因此系统是稳定系统。nn0(2)如果x(n)M,y(n)x(k)2n01M,因此系统是稳定的。系统是非因果的,knn0因为输出还和x(n)的将来值有关。注意:如果给出的是h(n),用上面要求记住的充要条件判断!例:设某线性时不变系统的单位取样响应为h(n)anu(n)(a为实数),分析系统的因果性和稳定性。(重点)解:讨论因果性:因为n0时,h(n)0,所以该系统是因果系统。讨论稳定性:1a1h(n)ann1aann0n0a1∵∴当a1时,系统是稳定的;否则,系统不稳定。例:设某线性时不变系统的单位取样响应为h(n)anu(n1)(a为实数),分析系统的因果性和稳定性。(重点)解:讨论因果性:因为n0时,h(n)0,所以该系统是非因果系统。讨论稳定性:第10页共60页n1111ah(n)ana()na1nnn1n1aa1∵∴当a1时,系统是稳定的;否则,系统不稳定。1.3线性常系数差分方程差分方程定义卷积和是一种LTI系统的数学模型,一般情况下,我们可以用差分方程描述LTI系统的输入输出NM关系。aky[nk]bkx[nk]k0k0差分方程给出了系统响应y[n]的内部关系。为得到y[n]的显式解,必须求解方程。差分方程求解○1经典法○2递推法○3变换域法(参见下章z域变换)(重点)例:设系统的差分方程为y(n)0.5y(n1)1.5x(n),输入序列为x(n)(n),求输出序列y(n)。解:一阶差分方程需一个初始条件。设初始条件为:y(1)0则y(0)0.5y(1)1.5x(0)1.5y(1)0.5y(0)1.5x(1)0.75y(2)0.5y(1)1.5x(2)0.375y(n)1.5(0.5)nu(n)设初始条件改为:y(1)1则y(0)0.5y(1)1.5x(0)2y(1)0.5y(0)1.5x(1)1y(2)0.5y(1)1.5x(2)0.5y(n)2(0.5)nu(n)该例表明,对于同一个差分方程和同一个输入信号,因为初始条件不同,得到的输出信号是不相同的。几点结论(重点)(1)对于实际系统,用递推解法求解,总是由初始条件向n>0的方向递推,是一个因果解。但对第11页共60页于差分方程,其本身也可以向n<0的方向递推,得到的是非因果解。因此差分方程本身不能确定该系统是因果系统还是非因果系统,还需要用初始条件进行限制。2)一个线性常系数差分方程描述的系统不一定是线性非时变系统,这和系统的初始状态有关。如果系统是因果的,一般在输入x(n)=0(n<n0)时,则输出y(n)=0(n<n0),系统是线性非时变系统。1.4模拟信号数字处理方法模拟信号数字处理框图xa(t):模拟信号输入预滤波:目的是限制带宽(一般使用低通滤波器)○1采样:将信号在时间上离散化A/DC:模/数转换○2量化:将信号在幅度上离散化(量化中幅度值=采样幅度值)○3编码:将幅度值表示成二进制位(条件fs2fc)数字信号处理:对信号进行运算处理D/AC:数/模转换(一般用采样保持电路实现:台阶状连续时间信号在采样时刻幅度发生跳变)平滑滤波:滤除信号中高频成分(低通滤波器),使信号变得平滑ya(t):输入信号经过处理后的输出信号2.连续信号的采样对连续信号进行理想采样,设采样脉冲,则采样输出在讨论理想采样后,信号频谱发生的变化时,可遵循下面的思路:1)由;2)由;3)根据频域卷积定理,由计算出。计算过程:第12页共60页1)2)周期信号可以用傅里叶级数展开,因此其中系数所以其傅里叶变换3)因此,采样后信号频谱产生周期延拓,周期为Ωs,同时幅度为原来的1/T倍。这是一个非常重要的性质,应熟练掌握。时域抽样定理(重点)一个限带模拟信号xa(t),若其频谱的最高频率为F0,对它进行等间隔抽样而得x(n),抽样周期为T,或抽样频率为Fs1/T;只有在抽样频率Fs2F0时,才可由xa(t)准确恢复x(n)。例:有一连续信号xa(t)cos(2ft),式中,f20Hz,2(1)求出xa(t)的周期。(2)用采样间隔T0.02s对xa(t)进行采样,试写出采样信号%xa(t)的表达式。(3)求出对应x%a(t)的时域离散信号(序列)x(n),并求出x(n)的周期。解:(1)xa(t)周期为T10.05sf第13页共60页^(tnT)fnT)(tnT)(T0.05s)(2)x(t)x(t)cos(2nn(3)x(n)的数字频率ω=0.8π,故225,因而周期N=5,所以x(n)=cos(0.8πn+π/2)0.82简答题:(重点)1.是不是任意连续信号离散后,都可从离散化后的信号恢复出原来的信号?为什么?2.一个连续时间信号经过理想采样以后,其频谱会产生怎样的变化?在什么条件下,频谱不会产生失真?3.说明时域采样定理的要点?4.离散信号频谱函数的一般特点是什么?5.画出模拟信号数字处理框图。并说明各部分的作用。名词解释:(重点)时域采样定理线性系统、时不变系统、稳定系统、因果系统第14页共60页第二章:本章涉及信号及系统的频域分析方法,概念较多,但很基础,学习时要注意。2.1序列的傅里叶变换的定义及性质1.定义DTFT是一个用来确定离散时间序列频谱的重要数学工具。物理意义:傅里叶变换是将对信号的时域分析转换为对其在频域的分析,便于研究问题。若序列满足绝对可和条件则其离散时间傅里叶变换(DiscreteTimeFourierTransform-DTFT:非周期序列的傅里叶变换)定义为X(ej)x[n]ejn------(记住!!)n反变换定义为:x[n]1X(ej)ejnd------2傅里叶变换对例:设x(n)R4(n),求其序列傅里叶变换。(重点)解第15页共60页jx(n)ejnjnX(e)DTFT[x(n)]RN(n)ennN1jNjN/2ejn1eeen01ejej/2eNsinN12ej2sin2当N4时

jN/2ejN/2j/2ej/2X(ej)sin2e3jsin/2(2-5)j2.1所示。X(e)的幅度和相位随变化曲线如图41arg[X(e)]20或X(ej)sin(4/2)sin(/2)x(n)1n0123arg[X(ej)]X(ej)42002图2.1R4(n)的幅度与相位曲线例:试求如下序列的傅里叶变换:(重点)(1)x1(n)(nn0)(2)x2(n)1(n1)(n)1(n1)22(3)x3(n)anu(n2),0a1(4)x4(n)u(n3)u(n4)解:X1(ej)(nn0)ejnejn0(1)nX2(ej)x2(n)ejn1ej11ej1jsin(2)n22X3(ej)anu(n2)ejnanejna2e2j,0a1(3)nn21aej第16页共60页333(4)X4(ej)u(n3)u(n4)ejnejnejnejnnn3n0n1e1e

4j

e1e

1ej31ejej=

j7sin7jej321sin22.性质1)周期性(重点):DTFT是关于ω的周期为2π的周期函数。X(ej)x(n)ej(2M)nX(ej(2M))M为整数nX1(ej)FT[x1(n)],X2(ej)FT[x2(n)],那么2)线性(重点):设FT[ax1(n)bx2(n)]aX1(ej)bX2(ej)3)时移特性(重点)4)频移特性5)时域卷积定理(重点)6)频域卷积定理7)帕斯瓦尔定理时域总能量等于频域一周期内总能量。幅度频谱为ω的偶函数,相位频谱为ω的奇函数。X(ejω)的实部为ω的偶函数,X(ejω)的虚部为ω的奇函数。对称关系的总结(重点):如果x[n]为复数序列,其DTFT为X(ejω),(a)x[n]实部的DTFT为X(ejω)的共轭对称部分-----------xre[n]Xcs(ej)1X(ej)X*(ej)2(b)x[n]虚部的DTFT为X(ejω)的反共轭对称部分-----------xim[n]Xca(ej)1X(ej)X*(ej)2第17页共60页(c)x[n]的共轭对称部分的DTFT为X(ejω)的实部-----------xcs[n]1x[n]x*[n]Xre(ej)2(d)x[n]的反共轭对称部分的DTFT为X(ejω)的虚部-----------xca[n]1x[n]x*[n]jXim(ej)2如果实序列x[n]的DTFT为X(ejω),(e)x[n]的偶对称部分的DTFT为X(ejω)的实部,-----------xev[n]1x[n]x[n]Xre(ej)2(f)x[n]的奇对称部分的DTFT为X(ejω)的虚部,-----------xod[n]1x[n]x[n]jXim(ej)2例:设系统的单位取样响应h(n)anu(n),0a1,输入序列为x(n)(n)2(n2),完成下面各题:(1)求出系统输出序列y(n);(2)分别求出x(n)、h(n)和y(n)的傅里叶变换。(重点)解:(1)y(n)h(n)*x(n)anu(n)*[(n)2(n2)]anu(n)2an2u(n2)X(ejw)[(n)2(n2)]ejwn12ej2wn(2)H(ejw)anu(n)ejwnanejwn11nn0aejwY(ejw)H(ejw)gX(ejw)12ej2w1aejw2.2时域离散信号的傅里叶变换与模拟信号傅里叶变换之间的关系:^X(ejT1Xa(jjks)2X(ejT)s2Fs)Xa(j)Tk式中T2.3序列的Z变换Z变换定义(重点)Z变换为离散时间信号与LTI系统分析的重要数学工具。给定一离散时间序列x(n),其z变换定义第18页共60页为:X(z)x(n)znRxzRx------(记住!!)k其中,zes,sj。z变换存在情况下的Z变量取值范围称为收敛域(ROC)。注意:Z变换+不同收敛域对应不同收敛域的不同序列唯一序列(Z变换+收敛域)(重点)例:求以下序列的Z变换及收敛域:(重点)1)2nu(n);(2)2nu(n1);(3)2n[u(n)u(n10)]解:(1)ZT[2nu(n)]2nu(n)zn2nzn11z1,z1nn0122ZT[2nu(n1)]2nu(n1)zn2nzn2nzn(2)nn1n12z1112z121z1,z29ZT[2nu(n)u(n10)]2nzn(3)n01210z10z121z1,0[说明]上题也可以改为求序列的傅立叶变换。可以利用X(ej)X(z)zej。2Z变换和DTFT之间的关系(重点)jDTFT为单位圆上的z变换。数学表达为:X(e)X(z)zej------记住并理解!序列特性与X(z)的收敛域ROC的关系。(重点)收敛区域要依据序列的性质而定。同时,也只有Z变换的收敛区域确定之后,才能由Z变换唯一地确定序列。一般来来说,序列的Z变换的收敛域在Z平面上的一环状区域:Rx|z|Rx总结:a.ROC不包含任何极点。b.有理z变换的收敛域ROC由其极点界定。第19页共60页c.对于有限长序列x[n],其z变换的收敛域ROC为整个z-平面,可能在z=0或z=∞除外。只有序列为(n)时,收敛域是整个Z平面。d.对于右边序列x[n],其z变换的收敛域ROC由其离原点最远的极点确定,其形式为zRx。e.对于左边序列x[n],其z变换的收敛域ROC由其离原点最近的极点确定,其形式为zRx。f.对于双边序列x[n],其z变换的收敛域ROC环状收敛域,,其形式为公共收敛域RxzRx。Z反变换(重点)常用序列的Z变换(重点--记住!!):Z[(n)]1,|z|01Z[u(n)]1z1,|z|1Z[anu(n)]11,|z||a|az1Z[bnu(n1)]11,|z||b|bz1逆变换x(n)1j?X(z)zn1dzx,C:收敛域内绕原点逆时针的一条闭合曲线c留数定理:x(n)[X(z)zn1在C内极点留数之和]留数辅助定理:x(n)[X(z)zn1在C外极点留数之和]第20页共60页利用部分分式展开:X(z)Ak,然后利用定义域及常用序列的Z变换求解。(重点)1akz1基本要求:用部分分式展开法求z反变换。(重点)例:假设X(z)11,收敛域ROC为0.3z0.5,则X(z)的z反变换0.5z110.3z11为((0.5)nu(n1)(0.3)nu(n))。(重点)说明:本题要求掌握序列的时域特性域z变换收敛域之间的对应关系。具体说,有限长序列的z变换的ROC是怎样的,右边序列的z变换的ROC是怎样的,因果序列的z变换的ROC是怎样的,左边序列的z变换的ROC是怎样的,反因果序列的z变换的ROC是怎样的。nu(n1)1ROC:z典型序列的z变换表达式是否记住了?1z1这两个典1nu(n)z1ROC:z1型z变换对,对求z变换或逆z变换非常重要。X(z)

zz例:已知z0.5z2,试求与X(z)对应的所有可能的序列x(n)。(重点)解:同一个Z变换函数,收敛域不同,对应的序列也不同。本题没有给定收敛域,所以必须先确定收敛域。X(z)有两个极点:z10.5,z22,因为收敛域总是以极点为边界,所以收敛域有以下三种情况:z0.5,0.5z2,z2,三种收敛域对应三种不同的原序列,分别讨论如下:(1)z0.5对应左边序列∴x(n)0.5nu(n1)2nu(n1)(2)0.5z2对应双边序列∴x(n)0.5nu(n)2nu(n1)(3)z2对应右边序列∴x(n)0.5nu(n)2nu(n)X(z)1(12z1)(10.5z1)z2,用部分分式展开法求逆Z变换。(重点)例:设解:先去掉z的负幂次,以便于求解,将X(z)的分子分母同乘以z2,得:z2X(z)0.5)(z2)(zX(z)zA1A2将等式两端同时除以z,得:z(z2)(z0.5)z2z0.5第21页共60页A1X(z)(zX(z)z4Res[,2]2)(z2)zzz2(z2)(z0.5)z23A2X(z)(z0.5)X(z)(z0.5)z1Res[,0.5]zz0.52)(z0.5)z23z(z4z1zX(z)23z0.5因而得:3zx(n)42nu(n)10.5nu(n)由收敛域知,x(n)为右边序列,得:33主要应用于单阶极点的序列。Z变换的性质○线性性质(重点)M(z)ZT[m(n)]mm1aX(z)bY(z)RzR○2序列的移位性质(重点)X(z)ZT[x(n)]RxzRxZT[x(nn0)]zn0X(z)RxzRx3序列乘以指数序列的性质(重点)○X(z)ZT[x(n)]RxzRxy(n)anx(n)a为常数Y(z)ZT[anx(n)]X(a1z)aRxzaRx○X(z)ZT[x(n)]xzx4序列乘以n的ZTRRZT[nx(n)]zdX(z)RxzRxdx○X(z)ZT[x(n)]xzx5复共轭序列的ZTRRZT[x*(n)]X*(z*)RxzRx○6初值定理X(z)ZT[x(n)]x(0)limX(z)z第22页共60页○终值定理limx(n)lim(z1)x(z)7zz1○8时域卷积定理(重点)设(n)x(n)*y(n)X(z)ZT[x(n)]RxzRxY(z)ZT[y(n)]RxzRx则W(z)ZT[(n)]X(z)Y(z)RwzRw○X(z)xzx9复卷积定理ZT[x(n)]RRZT[y(n)]Y(z)RyzRy(n)x(n)y(n)W(z)21j?cX()Y(z)dRxRyzRxRy○X(z)xzx10帕斯维尔定理ZT[x(n)]RRZT[y(n)]Y(z)RyzRyRxRy1,RxRy1那么x(n)y*(n)21cX()Y*(1*)1dnj?2.4离散时间系统的系统函数及频率响应系统函数定义(重点)一个线性时不变离散时间系统在时域中可以用它的单位取样响应h(n)来表征,即:y(n)x(n)h(n)对等式两边取Z变换并根据时域卷积定理,有:Y(z)X(z)H(z)Y(z)H(z)一般称H(z)为系统的系统函数(系统零状态响应的则:X(z)Z变换与输入的Z变换之比),它表征了系统的复频域特性。第23页共60页系统函数与差分方程的关系aky(nk)bkx(nk)(给定差分方程,能计算其系统函数,或给定系统函数,能计算得到k0k0差分方程。)(重点)频率响应(重点)频率响应是一个重要的概念,根据频率响应,可理解滤波。频率响应定义为系统单位冲激响应的DTFT:H(ej)h[n]ejnH(ejj)ejH(e)(重点)n其中,|H(ejω)|称为幅频响应,H(ej)称为相频响应。系统的频率响应是以2π为周期的ω的连续函数,这一点和连续系统的频率响应是不同的,学习时应加以注意。若h(n)为实数,则系统的幅度响应在区间内是偶对称的,而相位响应是奇对称的。注意:仅当稳定系统才有频率响应。频率响应H(ejω)可根据DTFT与z变换之间的关系简单得到:jjX(e)X(z)zejH(e)H(z)zej稳态响应的求解结论:对于LTI系统,如果输入为正弦序列x(n)=cos(ω00则输出响应y(n)必为相同形式的正弦序t+φ),列,但需在ω=ωjωH(ej)在ω=ω0的值进行移位,0的幅频响应|H(e)|进行加权,并通过相频响应即:y[n]=|H(ejω000j0)))|cos(ωt+φ+H(e例:假设实序列x[n]的DTFT记为X(ej),则其幅值X(ej)是关于ω的(偶函数)。说明:还记得反复强调的一句话,实序列的DTFT的幅度、实部是关于频率ω偶函数,而相位和虚部则是关于频率ω奇函数。例:对于一LTI离散时间系统其频率响应H(ej)1,如果系统输x(n)=cos(n),响应10.5ej3的稳态输出响应y(n)=(1.15cos(n0.52))。3H(ej)H(ej)ejH(ej)说明:将系统的频率响应写成幅度相位表达式:,则输出信号为:第24页共60页jjy[n]H(e3)cos(nH(e3))3。这里由于给出了H(ej)的具体表达式,所以需要分别计算jH(e3)j出和H(e3)之值。用系统函数极点分布分析系统的因果性和稳定性(重点)MY(z)bizii0(传输函数H(z)为系统的单位冲激响应h(n)的Z变换。)系统函数:H(z)MX(z)aizik01)稳定系统:有界的输入产生的输出也有界的系统,即:若|x(n)|,则|y(n)|线性移不变系统是稳定系统的充要条件:|h(n)|n或:其系统函数H(z)的收敛域包含单位圆|z|=1(牢记此结论!)2)因果系统:n0时刻的输出y(n0)只由n0时刻之前的输入x(n),nn0决定线性移不变系统是因果系统的充要条件:h(n)0,n0或:其系统函数H(z)的收敛域在某圆外部:即:|z|>Rx(牢记此结论!)3)稳定因果系统:同时满足上述两个条件的系统。线性移不变系统是因果稳定系统的充要条件:|h(n)|,h(n)0,n0n或:H(z)的极点在单位圆内H(z)的收敛域满足:|z|Rx,Rx1(牢记此结论!)例:.一因果LTI离散时间系统的传输函数H(z)1n10.5z1,则系统的单位冲激响应为(0.5u(n))。说明:根据传递函数求系统的单位冲激响应,其实就是将传递函数进行逆z变换,但要注意系统的因果性如何。1例:因果IIR离散时间LTI系统,其传输函数H(z)10.5z

,则系统(稳定)。例:一FIR离散时间LTI系统总是(稳定)。说明:系统的稳定性如何判断?按照教材中的说法,就是系统传递函数的收敛域如果包括“单位圆”,则系统是稳定的。如果你熟悉了序列的z变换的ROC的性质,则此题不难回答。对于因果系统来说,其单位冲激响应为因果序列,故其z变换的ROC一定是某圆外部的整个区域。而这个圆就位于离原点最远的极点上,所以,对于因果系统,如果系统传递函数的全部极点都位于单位圆以内的话,则系统是稳定的。对于FIR系统,其单位冲激响应是一个有限长序列,其z变换的ROC为除了无穷远和原点之外的第25页共60页整个z平面,自然包括单位圆,所以FIR系统始终是稳定的。系统的频率特性可由系统函数零点及极点确定MMMbizi(1ziz1)(zzi)zMX(z)i0Ai1Ai1NNNakzk(1zkz1)(zzk)zNk0k1k1(式中,zk是极点,zi是零点;在极点处,序列x(n)的Z变换是不收敛的,因此收敛区域内不应包括极点。)系统函数H(z)的极点位置主要影响频响的峰值位置及尖锐程度,零点位置主要影响频响的谷点位置及形状。例:设一阶系统的差分方程为y(n)x(n)ay(n1),0a1,用几何法分析其幅频特性。(重点)解:对差分方程两边取Z变换,得:Y(z)X(z)az1Y(z)系统函数为:左所示:

Y(z)1zH(z)1az1zaza,极点为za,零点为z0,如下图X(z)当0时,由于极点矢量长度最短,幅频特性出现峰值,随着的增加,幅度逐渐减小,当时,由于极点矢量长度最长,幅频特性出现谷值,随着的增加,幅度逐渐增大,直到2时,幅频特性出现峰值,如上图右所示。简答题:(重点)1.说明有限长序列、左边序列、右边序列、双边序列的概念和收敛域各是什么?2.说明系统频率响应的概念?系统的频率响应和系统函数是什么关系?(单位圆上(zej)的系统函数就是系统的频率响应)说明FIR系统为什么始终是稳定的?怎样在z域表示离散时间LTI系统?答案:传输函数H(z)表示离散时间LTI系统。第26页共60页第三章:DFT是为适应计算机分析傅里叶变换规定的一种专门运算,本章是数字信号处理课程的重点章节。前言信号处理中会遇到几种信号形式:(1)连续周期信号(2)连续非周期信号(3)离散非周期信号(4)离散周期信号(重点)各种信号在时域和频域之间总的来说都是傅里叶变换,但具体形式及应用是不同的。1.连续周期信号——傅里叶级数(FS)~连续周期信号x(t)可展开成傅里叶级数:~cnejn0tx(t)(*)n20~T0T0式中,,为x(t)的周期。cn1T0/2傅里叶级数的系数为:T0T0/2

~jn0tdtx(t)e幅度频谱是指各次谐波的振幅随频率的变化关系,即:cn~(n0)时域频域连续时间、周期非周期、离散频率2.连续非周期信号——傅里叶变换(FT)第27页共60页连续非周期信号x(t)的傅里叶变换为:X(j)x(t)ejtdt20T00)变成因为非周期可视为,则离散频谱间距T0,则X(j的连续函数。时域频域连续时间、非周期非周期、连续频率3.离散非周期信号——序列的傅里叶变换(DTFT)X(ej)x(n)ejnn如果把序列看成连续时间信号的采样,采样间隔为T,则数字频率和模拟角频率的关系为T,且x(n)x(nT),代入上式,得:X(ejT)x(nT)ejnTn时域频域离散时间、非周期周期、连续频率4.离散周期信号——离散傅里叶级数(DFS)~N的周期序列,即:~~rN)r为任意整数设x(n)是周期为x(n)x(n表3.1四种傅里叶变换形式的归纳一般规律:一个域的离散对应另一个域的周期延拓,一个域的连续必定对应另一个域的非周期。(重点)3.1离散傅里叶级数1.周期序列的离散傅里叶级数(DFS)说明:离散傅里叶级数系数,用DFS(DiscreteFourierSeries)表示。第28页共60页连续时间周期信号可以用傅里叶级数表示,离散周期序列也可以表示成傅里叶级数形式。周期为N的复指数序列的基频序列为次谐波序列为由于,即,因而,离散傅里叶级数的所有谐波成分中只有N个是独立的。因此在展开成离散傅里叶级数时,我们只能取N个独立的谐波分量,通常取k=0到(N-1),即*)式中,1/N是习惯上采用的常数,是k次谐波的系数。利用将(*)式两端同乘以,并对一个周期求和即由于所以也是一个以N为周期的周期序列。因此,时域离散周期序列的离散傅里叶级数在频域上仍然是一个周期序列。称为离散傅里叶级数系数,用DFS(DiscreteFourierSeries)表示。令,则其中,符号DFS[.]表示离散傅里叶级数正变换,IDFS[.]表示离散傅里叶级数反变换。例:设x(n)R4(n),将x(n)以N8~~为周期进行周期延拓,得到周期序列x(n),求x(n)的DFS。解:第29页共60页~7~2knX(k)=x(n)ejn03j4knen01ejk441ejk41ejk1ejk4jkjkejke2(e22)jkjkejke8(e88)j3ksin2ke8sin8k其幅度特性为:2.周期序列的傅里叶变换思路:由利用和DTFT的频移特性,可得第30页共60页傅里叶变换时域、频域对应关系:根据序列的傅里叶变换和离散傅里叶级数频域特性,再结合连续时间信号的傅里叶变换频域特性,我们可以得出傅里叶变换时、频域的一般对应关系:连续→非周期,离散→周期。这种对应关系很重要,要求熟记(重点)。3.2有限长序列的离散傅立叶变换(DFT)说明:(DiscreteFourierTransform,DFT离散傅里叶变换)定义(重点)N1x(n)WNkn,0≤k≤N1------(记住!!)X(k)DFT[x(n)]{DFS[x(nN)]}RN(k)n0x(n)IDFT[X(k)]{IDFS[X(kN)]}RN(n)1N1X(k)WNkn,0≤n≤N1------Nk0记住!j2其中,WNeN应当注意,虽然x(n)和X(k)都是长度为N得有限长序列,但他们分别是由周期序列xp(n)和Xp(k)截取其主周期得到的,周期为N的周期序列xp(n)可以看成长度为N的有限长序列x(n)周期延拓的结果。本质上是做DFS或IDFS,所以不能忘记它们的隐含周期性。尤其是涉及其位移特性时更要注意。(重点)DFT的隐含周期性:WNkWNkmNk,m为整数,N为自然数(重点)例:设x(n)R4(n),求x(n)的4点DFT。(重点)解:x(n)的4点离散傅里叶变换为:N1j2kn3j2j3sin(k)X(k)Nknkx(n)ee4e4n0n0sin(k)k0,1,2,34第31页共60页以N4为周期将x(n)延拓成周期序列,得~x(n):~N1~j2323X(k)x(n)eknejknejkn0n0其离散傅里叶级数为:

sin(k)sin(k)4例:设x(n)R4(n),求x(n)的8点DFT。(重点)解:x(n)的8点离散傅里叶变换为:N1j23j2j3X(k)knknkx(n)eNe8e8n0n0

sin(k)2sin(k)8k0,1,2,,7第32页共60页以N8为周期将x(n)延拓成周期序列,得~x(n):其离散傅里叶级数为:N1jkn3jknjksin(k)~~2232Ne8e8X(k)x(n)en0n0sin(k)8由例可见,离散傅里叶变换的结果与变换区间长度N的取值有关。离散傅立叶变换与DTFT、Z变换的关系(重点)X(k)X(j)|2kX(z)|Nze

j2kNDFT的物理意义:X(k)为x(n)的傅里叶变换X(ej)在区间[0,2]上的等间隔采样。X(k)为第33页共60页X(z)在Z平面单位圆上的N点等间隔采样。时域分析记住结论:时域抽样对应频域的周期拓展,频率抽样对应时域的以周期N的周期拓展。y(n)x(nmN)这可以表述为如下公式:m3.3离散傅里叶变换的基本性质线性性质若y(n)ax1(n)bx2(n)则Y(k)DFT[y(n)]aX1(k)bX2(k)循环移位性质设x(n)是长度为M的有限长序列,则x(n)的N点循环移位定义为(NM):y(n)x((nm))NRN(n)循环移位的实现步骤:序列点数M不够时补零,补到所需点数Nx(n)(补充N-M个零点)将x(n)以N为周期延拓为周期序列~x(n)x((n))N移位~x(nm)x((nm))N取主值序列y(n)x((nm))NRN(n)循环卷积定理(重点)1)设序列h(n)和x(n)的长度分别为N和M。h(n)与x(n)的L点循环卷积定义为第34页共60页L1yc(n)[h(m)x((nm))L]RL(n)m0式中,L称为循环卷积区间长度,L≥max[N,M]。2)循环卷积矩阵y(0)cx(0)x(L1)x(L2)Lx(1)h(0)y(1)cx(1)x(0)x(L1)Lx(2)h(1)y(2)c=x(2)x(1)x(0)Lx(3)h(2)MMMMOMMy(L1)cx(L1)x(L2)x(L3)Lx(0)h(L1)特点:(1)第1行是序列{x(0),x(1),,x(L-1)}的循环倒相序列。注意,如果x(n)的长度M<L,则需要在x(n)末尾补L-M个零后,再形成第一行的循环倒相序列。(2)第1行以后的各行均是前一行向右循环移1位形成的。3)矩阵的各主对角线上的序列值均相等。循环卷积和线性卷积的区别线性卷积:翻折—>乘加—>移位:y(n)=x(n)*h(n)=∑h(k)x(n-k)循环卷积:补零—>周期延拓—>翻折—>循环移位—>对应值相加例:计算下面给出的两个长度为4的序列h(n)与x(n)的4点和8点循环卷积。(重点)yc(0)1432110yc(1)2143110yc(2)3214110解:按照循环卷积矩阵写出h(n)与x(n)的4点循环卷积矩阵形式为yc(3)4321110yc(0)1000043211yc(1)2100004313yc(2)3210000416yc(3)43210000110yc(4)0432100009yc(5)0043210007yc(6)0004321004h(n)与x(n)的8点循环卷积矩阵形式为yc(7)0000432100【补充】①计算h(n)与x(n)的线性卷积?②哪一种情况下计算的循环卷积结果就等于线性卷积?【说明】当循环卷积区间长度L大于等于y(n)=h(n)*x(n)的长度时,循环卷积结果就等于线性卷积。假设h(n)和x(n)都是有限长序列,长度分别是N和M。循环卷积等于线性卷积的条件是L≥N+M-1。(重点)3)时域循环卷积定理设h(n)和x(n)的长度分别为N和M,其L点循环卷积为L1yc(n)h(n)○x(n)h(m)x((nm))LRL(n)Lm0第35页共60页H(k)DFT[h(n)]L0kL1且DFT[x(n)]LX(k)则由DFT的循环卷积定理有Yc(k)DFT[yc(n)]LH(k)X(k)0kL14复共轭序列的DFT性质:设x*(n)是x(n)的复共轭序列,长度为N,X(k)DFT[x(n)]N,则DFT[x*(k)]NX*(Nk)0kN1例:给定一16-点实序列x(n),其16-点DFT记为X(k),已知X(13)=2+j3,则X*(3)=(2+j3)。说明:DFT的性质。实序列的DFT的共轭对称性:X(k)=X*(N-k),或X(N-k)=X*(k)。DFT的共轭对称性(重点)可总结出DFT的共轭对称性质:如果序列x(n)的DFT为X(k),则x(n)的实部和虚部(包括j)的DFT分别为X(k)的共轭对称分量和共轭反对称分量;而x(n)的共轭对称分量和共轭反对称分量的DFT分别为X(k)的实部和虚部乘以j。3.4频域采样定理离散傅里叶变换相当于信号傅里叶变换的等间隔采样,也就是说实现了频域的采样,便于计算机计算。那么是否任一序列都能用频域采样的方法去逼近呢?这是一个很吸引人的问题。我们考虑一个任意的绝对可和的序列x(n),它的z变换为如果对X(z)单位圆上进行等距离采样现在要问,这样采样以后,信息有没有损失?或者说,采样后所获得的有限长序列xN(n)能不能代表原序列x(n)。为了弄清这个问题,我们从周期序列开始第36页共60页由于所以也即是原非周期序列x(n)的周期延拓序列,其时域周期为频域采样点数N。在第一章我们看到,时域的采样造成频域的周期延拓,这里又对称的看到,频域采样同样造成时域的周期延拓。因此,如果序列x(n)不是有限长的,则时域周期延拓时,必然造成混叠现象,因而一定会产生误差。对于长度为M的有限长序列,只有当频域采样点数N大于或等于序列长度M时,才有即可由频域采样值X(k)恢复出原序列x(n),否则产生时域混叠现象,这就是所谓的频域采样定理。(重点)1zNN1X(k)X(z)1Wkz1内插公式:Nk0N3.5DFT的应用举例1.用DFT计算线性卷积(重点)用循环(周期)卷积计算有限长序列的线性卷积(重点)对周期要求:NN1N21(N1、N2分别为两个序列的长度)(记住!!)2.用DFT进行谱分析的误差问题(重点)(1)混叠现象利用DFT逼近连续时间信号的傅里叶变换,为避免混叠失真,按照抽样定理的要求,采样频率至少是信号最高频率的两倍。解决混叠问题的唯一方法是保证采样频率足够高。(2)截断效应任何带限信号都是非时限的,任何时限信号都是非带限的。实际问题中遇到的离散时间序列可能是非时限的、无限长序列,在对该序列利用DFT进行处理时,由于作DFT的点数总是有限的,因此就有一个必须将该序列截断的问题。序列截断的过程相当于给该序列乘上一个矩形窗口函数RN(n)。如果原来序列的频谱为,矩形窗函数的频谱为,则截断后有限长序列的频谱为第37页共60页截断后序列的频谱与原序列频谱必然有差别,这种差别对谱分析的影响主要表现在如下两个方面:①频谱泄露:由于矩形窗函数频谱的引入,使卷积后的频谱被展宽了,即的频谱“泄露”到其它频率处,称为频谱泄露。在进行DFT时,由于取无限个数据是不可能的,所以序列的时域截断是必然的,泄露是难以避免的。为了尽量减少泄露的影响,截断时要根据具体的情况,选择适当形状的窗函数,如汉宁窗或汉明窗等。②谱间干扰。在主谱线两边形成很多旁瓣,引起不同频率分量间的干扰(简称谱间干扰),特别是强信号谱的旁瓣可能湮没弱信号的主谱线,或者把强信号谱的旁瓣误认为是另一频率的信号的谱线,从而造成假信号,这样就会使谱分析产生较大偏差。(3)栅栏效应由于DFT是有限长序列的频谱等间隔采样所得到的样本值,这就相当于透过一个栅栏去观察原来信号的频谱,因此必然有一些地方被栅栏所遮挡,这些被遮挡的部分就是未被采样到的部分,这种现象称为栅栏效应。由于栅栏效应总是存在的,因而可能会使信号频率中某些较大的频率分量由于被“遮挡”而无法得到反映。此时,通常在有限长序列的尾部增补若干个零值,借以改变原序列的长度。这样对加长的序列作DFT时,由于点数增加就相当于调整了原来栅栏的间隙,可以使原来得不到反映的那些较大的频率分量落在采样点上而得到反映。产生原因说明:由傅里叶变换理论知道,若信号持续时间有限长,则其频谱无限宽;若信号的频谱有限宽,则其持续时间必然为无限长。(重点)如果用DFT分析连续信号的频谱,在对连续信号采样时,无法满足采样定理,那么就会出现频谱混叠现象。解决混叠问题的唯一方法是保证采样频率足够高。当连续信号无限长或很长时,在对连续信号采样时,采样点数太多以致无法存储和计算,需要将信号截断,这样将导致频谱的泄漏现象。为了尽量减少泄露的影响,截断时要根据具体的情况,选择适当形状的窗函数,如汉宁窗或汉明窗等。用DFT计算连续信号的频谱只能得到采样点上的频谱,而不能看到整个频谱,这种现象称为栅栏效应。可以通过增加点数,因为点数增加就相当于调整了原来栅栏的间隙,可以使原来得不到反映的那些较大的频率分量落在采样点上而得到反映。3.用DFT进行谱分析的参数选择问题(重点)对模拟信号频谱的采样间隔,称之为频率分辨率。(1)在已知信号的最高频率fc(即谱分析范围)时,为了避免频率混叠现象,要求采样频率Fs满足:Fs>2fc。(2)采样频率Fs,采样点数N,谱分辨率F=Fs/N,如果保持采样点数N不变,要提高频谱分辨率(减小F),就必须降低采样频率,采样频率的降低会引起谱分析范围变窄和频谱混叠失真。如维持Fs不变,为提高频率分辨率可以增加采样点数N。因为NT=Tp,T=Fs-1,只有增加对信号的观察时间Tp,才能增加N。3)采样点数N>2fc/F4)最小记录时间Tp≥1/F例:用DFT对实信号进行谱分析,要求频率分辨率F50Hz,信号最高频率为fmax1000Hz,试确定以下参数:(1)最小记录时间;(2)最大取样间隔;(3)最少采样点数;(4)若要求频率分辨率提高一倍,求最少采样点数。(重点)Tpmin10.02s解(1)Fmax第38页共60页Tmax10.5ms(2)2fmaxNmin2fmaxTpmin40(3)FmaxTmaxNmin2fmax200080(4)Fmax25简答题:(重点)一个序列的DFT与序列的傅里叶变换之间的关系是什么?2.序列的DTFT和序列的z变换间的关系是什么?序列的DFT和序列的Z变换间的关系是什么?有限长序列x(n)的长度为M,对其进行频域采样,不失真的条件是什么?两个有限长序列x1(n),0nM,x2(n),0nN,对它们进行线性卷积,结果用y(n)表示,y(n)的长度是多少?如果进行循环卷积,那么什么时候线性卷积和循环卷积的结果相等?用DFT进行谱分析带来哪些误差问题?采取什么措施可以减少这些误差?时域采样定理的要点是什么?频域采样定理的要点是什么?第39页共60页第四章:快速傅里叶变换并不是一种新的变换,它是离散傅里叶变换的一种快速算法。4.1直接计算DFT的问题及改进的途径直接计算DFT,需要次复数乘法,次复数加法。直接计算离散傅里叶变换,由于计算量近似正比于N2,显然对于很大的N值,直接计算离散傅里叶变换要求的算术运算量非常大。(重点)我们可以利用系数WNnk的特性来改善离散傅里叶变换的计算效率。1)的对称性2)的周期性利用的对称性和周期性,将大点数的DFT分解成若干个小点数的DFT,FFT正是基于这个基本思路发展起来的。(重点)说明:快速傅里叶变换FFT(FastFourierTransform)分类:按时间抽取(DIT)算法和按频率抽取(DIF)算法。第40页共60页4.2基2FFT的算法原理和FFT运算特点1)数据要求:N2M2)计算效率(乘法运算次数:1NM,加法计算次数:NM)(复数运算)2(DFT运算:乘法运算次数:N2,加法计算次数:N2)(复数运算)对于算法原理,要求能够看懂分解流图。1时域抽取法如下:(DecimationInTime,DIT–FFT)设序列x(n)长度为N,且满足N=2M,M为正整数。按n的奇偶把x(n)分解为两个N/2点的子序列:x(2r)x1(r)rN10,1,L,2,Nx(2r1)x2(r)r0,1,L12则x(n)的DFT为NN1N1122X(k)DFT[x(n)]x(n)WNnkx(2r)WN2rkx(2r1)WN(2r1)kn0r0r0N1N122x1(r)WN2rkWNkx2(r)WN2rkr0r0N1N122x1(r)WNrk/2WNkx2(r)WNrk/2r0r0X1(k)WNkX2(k)所以X(k)X1(k)WNkX2(k)k0,1,N/21N1N1X1(k)2x1(r)WNrk/22x1(r)WNrk/2DFT[x1(r)]N/2r0r0N1N1X2(k)2x2(r)WNrk/22x2(r)WNrk/2DFT[x2(r)]N/2r0r0将X(k)又可以写为X(k)X1(k)WNkX2(k)k0,1,L,N12XkNX1(k)WNkX2(k)k0,1,L,N122第41页共60页上式将N点DFT分解为两个N/2点的DFT运算,运算过程如下图示利用蝶形运算求解。运算量:由按时间抽取的FFT流图可见:每级都由N/2个蝶形单元构成,因此每级都需要N/2次复数乘法和N次复数加法。这样,M级运算共需要:mFNMNlog2N复数乘法:22复数加法:aFNMNlog2N而直接计算DFT需要:复数乘法:N2复数加法:N(N1)以乘法为例,对FFT算法与直接DFT算法的运算量进行比较:NN2Nlog2NN2/(Nlog2N)22864125.41281638444836.6102410485765120204.8可以看出:当N越大时,FFT算法的优越性越突出。DIT-FFT算法与DFT运算量的比较直接计算DFT与FFT算法的计算量之比为N22NN越大,FFT的优点越为明显Nlog2N2Nlog2说明:掌握给定点数的基2DIT-FFT蝶形图8点DFT的完整FFT流图:第42页共60页2频域抽样法(DecimationInFrequency,DIF–FFT)将长度为N=2M的序列x(n)前后对半分开,其N点DFT可表示为N1N12N1X(k)x(n)WNnkx(n)WNnkx(n)WNnkn0n0Nn2N1N1NN22nkx(n)WNnkxnWN2n0n02N12NWNNk/2WNnkx(n)xnk0,1,,N1n02按k的奇偶可将X(k)分为两部分取偶数时N12Nr0,1,L,NWN2nrX(2r)x(n)xn1n022N12Nx(n)xnWNnr/2n02k取奇数时X(2r1)N12x(n)n0

N1NWNn(2r1)r0,1,L,N12x(n)xnn022xnNWnWnr2NN/2N/21x1(n)x(n)xnNX(2r)rn2x1(n)WN/2令得到n0NN/21WNnx2(n)x(n)xnX(2r1)x2(n)WNrn/22n0注:DIT—FFT与DIF—FFT特点比较(重点)第43页共60页相同之处:1)DIF与DIT两种算

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论