离散时间信号与系统的频域分析.ppt

《离散时间信号与系统的频域分析.ppt》由会员分享，可在线阅读，更多相关《离散时间信号与系统的频域分析.ppt（162页珍藏版）》请在三一办公上搜索。

1、第3章 离散时间信号与系统的频域分析,我们知道信号和系统的分析方法有两种，即时域分析方法和频率分析方法。在模拟领域中，信号一般用连续变量时间t的函数表示，系统则用微分方程描述。为了在频率域进行分析，用拉普拉斯变换和傅里叶变换将时间域函数转换到频率域。,时域离散信号和系统中，信号用序列表示，其自变量仅取整数，非整数时无定义，而系统则用差分方程描述。频域分析是用Z变换或傅里叶变换这一数学工具。其中傅里叶变换指的是序列的傅里叶变换，它和模拟域中的傅里叶变换是不一样的，但都是线性变换，很多性质是类似的。本章学习序列的傅里叶变换和Z变换，以及利用Z变换分析系统和信号频域特性。本章学习内容是本书也是数字信

2、号处理这一领域的基础。,3.1 序列的傅里叶变换3.2 序列的Z变换 3.3 Z变换的基本性质和定理3.4 逆Z变换3.5 Z变换、傅里叶变换、拉普拉斯变换的关系3.6系统函数与频率响应,补充内容:Fourier变换的几种可能形式,连续非周期 非周期连续傅里叶变换（FT）,周期连续 离散非周期傅里叶级数（FS）,离散非周期 周期连续序列的傅里叶变换,周期离散 离散周期离散傅里叶变换,连续时间、连续频率傅里叶变换(FT),时域连续函数造成频域是非周期的谱，而时域的非周期造成频域是连续的谱密度函数。,连续时间、离散频率傅里叶级数(FS),时域连续函数造成频域是非周期的谱，时域周期函数造成频域的离散

3、。,离散时间、连续频率序列的傅里叶变换,时域的离散化造成频域的周期延拓，而时域的非周期对应于频域的连续,离散时间、离散频率离散傅里叶变换,一个域的离散造成另一个域的周期延拓，因此离散周期序列的傅里叶变换的时域和频域都是离散的和周期的,四种傅里叶变换形式的归纳,3.1序列的傅里叶变换,序列的傅里叶变换的定义,理想采样信号的频谱是连续信号频谱的周期延拓,重复周期为s,式是周期函数 的傅里叶级数展开式，xa（nT）是傅里叶级数展开式的系数。,（）,序列傅里叶变换的定义,(3.1.4),为序列x(n)的傅里叶变换，可以用FT(Fourier Transform)缩写字母表示。,(3.1.3),上式即是

4、FT的逆变换。(3.1.4)和(3.1.5)式组成一对傅里叶变换公式。,(3.1.5),FT成立的充分必要条件是序列x(n)满足绝对可和的条件，即满足下式：,例3.1.1 设x(n)=RN(n)，求x(n)的FT,解：,设N=4，幅度与相位随变化曲线如图3.1所示。,(3.1.1),图 3.1 R4(n)的幅度与相位曲线,3.1.2 序列傅里叶变换的性质,那么,设,式中a,b为常数 2.时移与频移 设X(e j)=FTx(n)，那么,(3.1.6),(3.1.7),(3.1.8),1.线性,3、FT的周期性,(3.1.9),因此序列的傅里叶变换是频率的周期函数，周期是2。这样X(ej)可以展成

5、傅里叶级数，其实3.1.4式已经是傅里叶级数的形式，x(n)是其系数。,因为,所以,4.FT的对称性 在学习FT的对称性以前，先介绍什么是共轭对称与共轭反对称以及它们的性质。设序列xe(n)满足下式：xe(n)=x*e(-n)(3.1.10)则称xe(n)为共轭对称序列。为研究共轭对称序列具有什么性质，将xe(n)用其实部与虚部表示 xe(n)=xer(n)+jxei(n)将上式两边n用-n代替，并取共轭，得到 x*e(-n)=xer(-n)-jxei(-n),对比上面两公式，左边相等，因此得到 xer(n)=xer(-n)xei(n)=-xei(-n)(3.1.14)由上面两式得到共轭对称序

6、列其实部是偶函数，而虚部是奇函数。类似地，可定义满足下式的称共轭反对称序列 xo(n)=-x*o(-n)(3.1.11),将xo(n)表示成实部与虚部如下式：xo(n)=xor(n)+jxoi(n)可以得到 xor(n)=-xor(-n)xoi(n)=xoi(-n)(3.15)即共轭反对称序列的实部是奇函数，而虚部是偶函数。,例 试分析x(n)=e jn的对称性 解：将x(n)的n用-n代替，再取共轭得到：x*(-n)=e jn 因此x(n)=x*(-n)，满足(3.10)式，x(n)是共轭对称序列，如展成实部与虚部，得到 x(n)=cosn+j sinn 由上式表明，共轭对称序列的实部确实是

7、偶函数，虚部是奇函数。,对于一般序列可用共轭对称与共轭反对称序列之和表示，即 x(n)=xe(n)+xo(n)（3.1.17）式中xe(n),xo(n)可以分别用原序列x(n)求出，将(3.17)式中的n用-n代替，再取共轭得到 x*(-n)=xe(n)-xo(n)（3.1.18）利用(3.1.17)和(3.1.18)两式，得到,(3.1.19),(3.1.20),利用上面两式，可以分别求出xe(n)和xo(n)。对于频域函数X(ej)也有和上面类似的概念和结论：X(ej)=Xe(ej)+Xo(ej)(3.1.21)式中Xe(ej)与Xo(ej)分别称为共轭对称部分和共轭反对称部分，它们满足

8、Xe(ej)=X*e(e-j)(3.1.22)Xo(ej)=-Xo*(e-j)(3.1.23)同样有下面公式满足：,(3.1.24),(3.1.25),(a)将序列x(n)分成实部xr(n)与虚部xi(n)x(n)=xr(n)+jxi(n)将上式进行FT，得到 X(e j)=Xe(e j)+Xo(e j),式中,上面两式中，xr(n)和xi(n)都是实数序列，容易证明Xe(ej)满足(3.1.22)式，有共轭对称性，它的实部是偶函数，虚部是奇函数。Xo(ej)满足(3.1.23)式，具有共轭反对称性质，其实部是奇函数，虚部是偶函数。最后得到结论：序列分成实部与虚部两部分，实部对称的FT具有共轭

9、对称性，虚部和j一起对应的FT具有共轭反对称性。,(b)将序列分成共轭对称部分xe(n)和共轭反对称部分xo(n)，即 x(n)=xe(n)+xo(n)（3.1.26）将(3.19)式和(3.20)式重定如下：,将上面两式分别进行FT，得到FTxe(n)=1/2X(ej)+X*(ej)=ReX(ej)=XR(ej)FTxo(n)=1/2X(ej)-X*(ej)=jImX(ej)=jXI(ej)因此对(3.26)式进行FT得到：X(ej)=XR(ej)+jXI(ej)(3.1.27),(3.1.27)式表示序列的共轭对称部分xe(n)对应着FT的实部XR(ej)，而序列的共轭反对称部分xo(n)

10、对应着FT的虚部。,因为h(n)是实序列，其FT只有共轭对称部分He(ej)，共轭反对称部分为零。H(ej)=He(ej)H(ej)=H*(e-j)因此实序列的FT的实部是偶函数，虚部是奇函数，用公式表示为 HR(ej)=HR(e-j)HI(ej)=-HI(e-j),按照(3.1.19)和(3.1.20)式得到 h(n)=he(n)+ho(n)he(n)=1/2h(n)+h(-n)ho(n)=1/2h(n)-h(-n)因为h(n)是实因果序列，按照上面两式he(n)和ho(n)可以用下式表示：,(3.1.28),(3.1.29),实因果序列h(n)分别用he(n)和ho(n)表示为 h(n)=

11、he(n)u+(n)(3.1.30)h(n)=ho(n)u+(n)+h(o)(n)(3.1.31),(3.1.32),例 x(n)=anu(n);0a1;求其偶函数xe(n)和奇函数xo(n)。解：x(n)=xe(n)+xo(n)按(3.1.3)式得到,按照(3.1.29)式得到,图 2.2.3 例图,5.时域卷积定理 设 y(n)=x(n)*h(n),Y(ej)=X(e j)H(e j)(3.1.33)证明,令k=n-m,该定理说明，两序列卷积的FT，服从相乘的关系。对于线性时不变系统输出的FT等于输入信号的FT乘以单位脉冲响应FT。因此求系统的输出信号，可以在时域用卷积公式计算，也可以在频

12、域按照(3.1.33)式，求出输出的FT，再作逆FT求出输出信号。,6.频域卷积定理 y(n)=x(n)h(n)(3.1.34),证明：,7.帕斯维尔(Parseval)定理,(3.1.35),帕斯维尔定理告诉我们，信号时域的总能量等于频域的总能量。要说明一下，这里频域总能量是指|X(e j)|2在一个周期中的积分再乘以1/(2)。最后，表3.1综合了FT的性质，这些性质在分析问题和实际应用中是很重要的。,表 3.1 序列傅里叶变换的性质,3.2 序列的Z变换,3.2.1 Z变换的定义 序列x(n)的Z变换定义为,(3.2.1),式中z是一个复变量，它所在的复平面称为z平面。注意在定义中，对n

13、求和是在之间求和，可以称为双边Z变换。还有一种称为单边Z变换的定义，如下式,(3.2.2),使(3.2.3)式成立，Z变量取值的域称为收敛域。一般收敛域用环状域表示,这种单边Z变换的求和限是从零到无限大，因此对于因果序列，用两种Z变换定义计算出的结果是一样的。本书中如不另外说明，均用双边Z变换对信号进行分析和变换。(3.2.1)式Z变换存在的条件是等号右边级数收敛，要求级数绝对可和，即,(3.2.3),图 3.2.1 Z变换的收敛域,常用的Z变换是一个有理函数，用两个多项式之比表示 分子多项式P(z)的根是X(z)的零点，分母多项式Q(z)的根是X(z)的极点。在极点处Z变换不存在，因此收敛域

14、中没有极点，收敛域总是用极点限定其边界。对比序列的傅里叶变换定义(3.1.4)式，很容易得到FT和ZT之间的关系，用下式表示：,(3.2.4),式中z=e j表示在z平面上r=1的圆，该圆称为单位圆。(3.2.4)式表明单位圆上的Z变换就是序列的傅里叶变换。如果已知序列的Z变换，可用(3.2.4)式，很方便的求出序列的FT，条件是收敛域中包含单位圆。例 3.2.1 x(n)=u(n)，求其Z变换。解：X(z)存在的条件是|z-1|1，,|z|1,由x(z)表达式表明，极点是z=1，单位圆上的Z变换不存在，或者说收敛域不包含单位圆。因此其傅里叶变换不存在，更不能用(3.2.4)式求FT。该序列的

15、FT不存在，但如果引进奇异函数()，其傅里叶变换可以表示出来。该例同时说明一个序列的傅里叶变换不存在，在一定收敛域内Z变换是存在的。,3.2.2 序列特性对收敛域的影响 序列的特性决定其Z变换收敛域，了解序列特性与收敛的一些一般关系，对使用Z变换是很有帮助的。1.有限长序列 如序列x(n)满足下式：x(n)n1nn2 x(n)=0 其它,即序列x(n)从n1到n2序列值不全为零，此范围之外序列值为零，这样的序列称为有限长序列。其Z变换为,设x(n)为有界序列，由于是有限项求和，除0与点是否收敛与n1、n2取值情况有关外，整个z平面均收敛。如果n10，则收敛域不包括z=0点；如果是因果序列（n0

16、,x(n)=0），收敛域包括z=点。具体有限长序列的收敛域表示如下：,n10时，00时，0z,收敛域为除了0和 的整个 平面,例 求x(n)=RN(n)的Z变换及其收敛域 解：这是一个因果的有限长序列，因此收敛域为0z。但由结果的分母可以看出似乎z=1是X(z)的极点，但同时分子多项式在z=1时也有一个零点，极零点对消，X(z)在单位圆上仍存在，求RN(n)的FT，可将z=ej代入X(z)得到，其结果和例题中的结果(3.1.1)公式是相同的。,2.右序列 右序列是在nn1时，序列值不全为零，而其它nn1，序列值全为零。第一项为有限长序列，设n1-1，其收敛域为0|z|。第二项为因果序列，其收敛

17、域为Rx-|z|，Rx-是第二项最小的收敛半径。将两收敛域相与，其收敛域为Rx-|z|。如果是因果序列，收敛域定为Rx-|z|。,圆外为收敛域,例 求x(n)=anu(n)的Z变换及其收敛域 解：,在收敛域中必须满足|az-1|a|。3.左序列 左序列是在nn2时，序列值不全为零，而在nn2，序列值全为零的序列。左序列的Z变换表示为,如果n20，则收敛域为0|z|Rx+。,圆内为收敛域，若 则不包括z=0点,例 求x(n)=-anu(-n-1)的Z变换及其收敛域。,X(z)存在要求|a-1 z|1，即收敛域为|z|a|,4.双边序列 一个双边序列可以看作一个左序列和一个右序列之和，其Z变换表示

18、为,X(z)的收敛域是X1(z)和X2(z)收敛域的公共收敛区域。如果Rx+Rx-，其收敛域为Rx-|z|Rx+，这是一个环状域，如果Rx+Rx-，两个收敛域没有公共区域，X(z)没有收敛域，因此X(z)不存在。,例 3.2.5 x(n)=a|n|，a为实数，求x(n)的Z变换及其收敛域。解：,第一部分收敛域为|az|a|。如果|a|1，两部分的公共收敛域为|a|z|a|-1，其Z变换如下式：,|a|z|a|-1,如果|a|1，则无公共收敛域，因此X(z)不存在。当0a1时，x(n)的波形及X(z)的收敛域如图所示。,图 3.2.2 例图,例：,右边序列,例：,左边序列,收敛半径,圆内为收敛域

19、，若 则不包括z=0点,例：,双边序列,典型序列的Z变换,单位样值序列单位阶跃序列斜变序列指数序列正弦余弦序列,余弦序列的 Z 变换:,正弦序列的 Z 变换:,例,z变换X(z)及其收敛域才能唯一确定一个序列；X(z)的收敛域为Z平面以原点为中心的圆环；收敛域不包含任何极点（以极点为边界）；有限长序列的收敛域为整个平面（可能除去z=0和z=）；右边序列的收敛域在|z|=Rx-的圆外；左边序列的收敛域在|z|=Rx+的圆内；双边序列的收敛域在Rx-|z|Rx+的圆环。,总结,3.3 Z变换的基本性质和定理如果则有：,*即满足均匀性与叠加性；*收敛域为两者重叠部分。,1.线性,例已知,求其z变换。

20、,解：,2.序列的移位,如果则有：,证明,序列位移以后的Z变换的收敛域不发生变换，但是由于乘以因子，对有些单边序列在z=0和z=处可能有变化。,例 求序列x(n)=u(n)-u(n-3)的z变换。,例如 Z(n)=1,收敛域为整个z平面，即0|z|Z(n-1)=z-1,收敛域为0|z|Z(n+1)=z,收敛域为0|z|,3.Z域尺度变换(乘以指数序列),如果,，则,证明：,若a是实数，则表示收敛域的扩大或缩小，对零、极点而言只有幅度上的变换而没有相位上的变换；若a是模为1的复数，则表示零、极点只有相位上的变化，没有幅度上的变化，而且收敛域不变；若a是模不等于1的复数，则表示零、极点在相位与幅度

21、上都发生变化，相应的收敛域也会发生变化,4.序列的线性加权(Z域求导数),如果,，则,证明：,序列nx(n)的Z变换与X(z)的导数有关，且收敛域不发生变化。,5.共轭序列,如果,，则,证明：,6.翻褶序列,如果,，则,证明：,7.初值定理,证明：,8.终值定理,证明：,又由于只允许X(z)在z=1处可能有一阶极点，故因子（z-1)将抵消这一极点，因此(z-1)X(z)在上收敛。所以可取z 1的极限。,9.有限项累加特性,证明：,10.序列的卷积和(时域卷积定理),证明：,例已知网络的单位取样响应h(n)=anu(n),|a|1，网络输入序列x(n)=u(n)，求网络的输出序列y(n)。解：y

22、(n)=h(n)*x(n)求y(n)可用二种方法，一种直接求解线性卷积，另一种是用Z变换法。,由收敛域判定y(n)=0,n0。n0 y(n)=ResY(z)z n-1,1+ResY(z)z n-1,a,将y(n)表示为,11.序列相乘(Z域卷积定理),其中,C是在变量V平面上，H(z/v),X(v)公共收敛域内环原点的一条逆时针单封闭围线。,证明,由X(z)收敛域和Y(z)的收敛域，得到,例已知x(n)=u(n),y(n)=a|n|，若w(n)=x(n)y(n)，求W(z)=ZTw(n)解：,因此,W(z)收敛域为|a|z|；被积函数v平面上收敛域为max(|a|,0)|v|min(|a-1|

23、,|z|)，v平面上极点：a、a-1和z,c内极点z=a。,12.帕塞瓦尔定理(parseval),其中“*”表示复共轭，闭合积分围线C在公共收敛域内。,如果,则有:,证明 令 w(n)=x(n)y*(n)按照复卷积定理，得到,收敛域为R x-R y-|z|R x+R y+，按照假设，z=1在收敛域中，令z=1代入W(z)中。,如果x(n)和y(n)都满足绝对可和，即单位圆上收敛，在上式中令v=e j，得到,令x(n)=y(n)得到,上面得到的公式和在傅里叶变换中所讲的帕期维尔定理公式是相同的。上式还可以表示成下式：,3.4 逆Z变换一.定义：已知序列的Z变换及其收敛域，求序列称为逆Z变换,z

24、变换公式：,C为环形解析域内环绕原点的一条逆时针闭合单围线.,0,c,1.用留数定理求逆Z变换 如果X(z)zn-1在围线c内的极点用zk表示，根据留数定理,(3.4.6),式中 表示被积函数X(z)zn-1在极点z=zk的留数，逆Z变换则是围线c内所有的极点留数之和。如果zk是单阶极点，则根据留数定理,(3.4.7),F(z)在z平面上有N个极点，在收敛域内的封闭曲线c将z平面上极点分成两部分：一部分是c内极点，设有N1个极点，用z1k表示；另一部分是c外极点，有N2个，N=N1+N2，用z2k表示。根据留数辅助定理下式成立：,(3.4.9),注意(3.4.9)式成立的条件是F(z)的分母阶

25、次比分子阶次必须高二阶以上。设X(z)=P(z)/Q(z)，P(z)与Q(z)分别是M与N阶多项式。(3.4.9)式成立的条件是,N-M-n+12 因此要求 N-M-n1(3.4.10)如果(3.4.10)式满足,c圆内极点中有多阶极点，而c圆外极点没有多阶的，可以按照(3.4.9)式，改求c圆外极点留数之和，最后加一个负号。例 已知X(z)=(1-az-1)-1，|z|a，求其逆Z变换x(n)。,为了用留数定理求解，先找出F(z)的极点，极点有：z=a;当n0时z=0共二个极点，其中z=0极点和n的取值有关。n0时，n=0不是极点。n0时，z=0是一个n阶极点。因此分成n0和n0两种情况求x

26、(n)。n0 时，,n0时，增加z=0的n阶极点，不易求留数，采用留数辅助定理求解，检查(3.4.10)式是否满足，此处n0，只要N-N0，(3.4.10)式就满足。,图 2.5.4 例中n0时F(z)极点分布,例已知，求其逆变换x(n)。解：该例题没有给定收敛域，为求出唯一的原序列x(n)，必须先确定收敛域。分析X(z)，得到其极点分布如图所示。图中有二个极点z=a和z=a-1，这样收敛域有三种选法，它们是(1)|z|a-1|，对应的x(n)是右序列；(2)|a|z|a-1|，对应的x(n)是双边序列；(3)|z|a|，对应的x(n)是左序列。,图 2.5.5 例2.5.7 X(z)极点分布

27、图,下面按照收敛域的不同求其x(n)。(1)收敛域|z|a-1|,收敛域是因果的右序列，无须求n0时的x(n)。当n0时，围线积分c内有二个极点z=a和z=a-1，因此,最后表示成：x(n)=(an-a-n)u(n)。(2)收敛域|z|a|这种情况原序列是左序列，无须计算n0情况，当n0时，围线积分c内没有极点，因此x(n)=0。n0时，c内只有一个极点z=0，且是n阶极点，改求c外极点留数之和,最后将x(n)表示成 x(n)=(a-n-an)u(-n-1)(3)收敛域|a|z|a-1|这种情况对应的x(n)是双边序列。根据被积函数F(z)，按n0和n0两情况分别求x(n)。n0时，c内极点z

28、=a x(n)=ResF(z),a=an,n0时，c内极点有二个，其中z=0是n阶极点，改求c外极点留数，c外极点只有z=a-1，因此 x(n)=-ResF(z),a-1=a-n 最后将x(n)表示为 an n0 x(n)=x(n)=a|n|a-n n0,例,解,必然是因果序列，右边序列,2.幂级数展开法(长除法)因为 x(n)的Z变换为Z-1 的幂级数，即 所以在给定的收敛域内，把X(z)展为幂级数，其系数就是序列x(n)。如收敛域为|z|Rx+，x(n)为因果序列，则X(z)展成Z的负幂级数。若 收敛域|Z|Rx-,x(n)必为左边序列，主要展成 Z的正幂级数。,例 已知 用长除法求其逆Z

29、变换x(n)。解由收敛域判定这是一个右序列，用长除法将其展成负幂级数,1-az-1,例 已知求其逆Z变换x(n)。解：由收敛域判定，x(n)是左序列，用长除法将X(z)展成正幂级数,1-az-1,例 试用长除法求的z反变换。,解:收敛域为环状，极点z=1/4对应因果序 列，极点z=4对应左边序列(双边序列),*双边序列可分解为因果序列和左边序列。*应先展成部分分式再做除法。,3.部分分式法 有理式：数字和字符经有限次加、减、乘、除运算 所得的式子。有理分式：含字符的式子做分母的有理式，或两个多项 式的商。分子的次数低于分母时称为真分式。部分分式：把x的一个实系数的真分式分解成几个分式 的和，使

30、各分式具有 或 的形式，其中x2+Ax+B是实数范围内的不可约 多项式，而且k是正整数。这时称各分式为原分 式的“部分分式”。,对于大多数单阶极点的序列，常常用这种部分分式展开法求逆Z变换。设x(n)的Z变换X(z)是有理函数，分母多项式是N阶，分子多项式是M阶，将X(z)展成一些简单的常用的部分分式之和，通过查表(参考表3.2)求得各部分的逆变换，再相加即得到原序列x(n)。,通常，X(z)可表成有理分式形式：,因此，X(z)可以展成以下部分分式形式其中，MN时，才存在Bn；Zk为X(z)的各单极点，Zi为X(z)的一个r阶极点。而系数Ak，Ck分别为：,的z反变换。,例利用部分分式法，求,

31、解：,分别求出各部分分式的z反变换（可查 P37表3.2），然后相加即得X(z)的z反变换。,例已知，求逆Z变换。,解,因为收敛域为22。第二部分极点z=-3，收敛域应取|z|3。查表3.2得到 x(n)=2nu(n)+(-3)nu(-n-1),利用Z变换解差分方程 在第一章中介绍了差分方程的递推解法，下面介绍Z变换解法。这种方法将差分方程变成了代数方程，使求解过程简单。设N阶线性常系数差方程为,(2.5.30),1.求稳态解 如果输入序列x(n)是在n=0以前时加上的，n时刻的y(n)是稳态解，对(2.5.30)式求Z变换，得到,式中,(2.5.31),(2.5.32),2.求暂态解 对于N

32、阶差分方程，求暂态解必须已知N个初始条件。设x(n)是因果序列，即x(n)=0,n0，已知初始条件y(-1),y(-2)y(-N)。对(2.5.30)式进行Z变换时，注意这里要用单边Z变换。方程式的右边由于x(n)是因果序列，单边Z变换与双边Z变换是相同的。下面先求移位序列的单边Z变换。设,(2.5.33),按照(2.5.33)式对(2.5.30)式进行单边Z变换,(2.5.34),例已知差分方程y(n)=by(n-1)+x(n)，式中x(n)=anu(n),y(-1)=2，求y(n)。解：将已知差分方程进行Z变换,式中，,于是,收敛域为|z|max(|a|,|b|)，,式中第一项为零输入解，

33、第二项为零状态解。,3.5 Z变换、拉氏变换、傅氏变换的关系,一.Z变换与拉氏变换的关系1.理想抽样信号的拉氏变换设 为连续信号，为其理想抽样信号，则,序列x(n)的z变换为，考虑到,显然，当 时，序列x(n)的 z 变换就等于理想抽样信号的拉氏变换。,2.Z变换与拉氏变换的关系(S、Z平面映射关系）S平面用直角坐标表示为：Z平面用极坐标表示为：又由于 所以有：,因此，;这就是说，Z的模只与S的实部相对应,Z的相角只与S虚部相对应。,=0,即S平面的虚轴 r=1,即Z平面单位圆；,0,即S的左半平面 r1,即Z的单位圆内；,0,即S的右半平面 r1,即Z的单位圆外。,(1).r与的关系,=0，

34、S平面的实轴，=0，Z平面正实轴；=0(常数)，S:平行实轴的直线，=0T,Z:始于 原点的射线；S:宽 的水平条带，整个z平面.,0,jImZ,ReZ,(2).与的关系（=T）,序列的傅氏变换与拉氏变换（双边）的关系,序列的傅氏变换可以看作是其拉氏变换（双边）在虚轴上的特例,二.Z变换和傅氏变换的关系,连续信号经理想抽样后,其频谱产生周期延拓，即 我们知道,傅氏变换是拉氏变换在虚轴S=j 的特例,因而映射到Z平面上为单位圆。因此,这就是说,（抽样）序列在单位圆上的Z变换,就等 于理想抽样信号傅氏变换。用数字频率作为Z平面的单位圆的参数，表示Z平面的辐角，且。,所以，序列在单位圆上的Z变换为序

35、列的傅氏变换。,系统函数,系统函数与系统的频率响应,本节将以系统函数和传输函数为核心来研究系统的变换域分析方法，它们分别是h(n)的Z变换和傅立叶变换。,1、系统函数：若系统单位脉冲响应为h(n)，则线性时不变离散系统零状态响应的输入输出关系为：,称H(z)为线性时不变离散系统的系统函数，它是单位脉冲响应的Z变换，即：,2、系统的频率响应（传输函数）,传输函数,若给系统输入单频率的复信号，则系统的输出为：,物理意义,结论：当输入为一个单频率的信号时，输出亦为同一频率的信号，但其幅度与相位都因为 的加权而发生了变化，且 的值是随频率的变化而变化的。,系统的因果性与稳定性,系统函数,H(z)有2个

36、极点，和，给定的收敛域为，包括无穷远点，故系统为因果系统。但收敛域不包括单位圆，因此系统是不稳定的。,解：,系统函数,若将收敛域改为，这时，收敛域包括单位圆，但不包括无穷远点，此时系统稳定但非因果。实际上这时系统的单位脉冲响应为，显然不是因果的。,该例表明:同一个系统函数，如果收敛域不同，系统的特性是完全不同的。由于任何物理可实现系统都必定是因果的，对于这种非因果但稳定的系统，有时可采用将单位脉冲响应截取一段后保存在存储器中，通过延时使之变成因果系统来近似实现。,系统函数的零极点与频率响应,极零点,对其进行因式分解得：,H(z)的零点,H(z)的极点,系统的频率响应,对于稳定系统，其极点应全部

37、位于单位圆内，或单位圆包括在收敛域内，其傅立叶变换存在。将 代入H(z)，得到系统的频率响应：,极零点分布与系统的频率响应,系统函数的零极点与频率响应,系统函数的零极点与频率响应,系统的频率响应特性,系统函数的零极点与频率响应,由幅频特性可知：当频率点变到极点附近时，Pk就变小，就会在该极点附近的频率出现峰值，极点越接近单位圆，峰值就越尖锐；同样，当频率点变到零点附近时，Qk就变小，就会在该零点附近的频率出现低谷，当零点在单位圆上时，该零点就是传输零点。可见在单位圆附近的零极点对系统的幅频特性有较大的影响。零点可在单位圆内外，但极点只能在单位圆内，否则系统将不稳定。而系统的相位响应对幅度特性没

38、有影响。,例 设一阶系统的差分方程为:,解:对差分方程两边取Z变换:,，a为实数,求系统的频率响应。,这是一因果系统，其单位抽样响应为而频率响应为:幅度响应为:相位响应为:,零极点分布情况,0,0,-1,0,a,1,系统的分类,IIR系统和FIR系统的定义,根据离散时间系统的单位脉冲响应h(n)在时域中的长度可将其分为两种类型：当h(n)的长度为无限长时称为无限长脉冲响应系统，简称为IIR系统。当h(n)的长度为有限长时称为有限长脉冲响应系统，简称为FIR系统。,通常可以根据系统函数的零极点来判断系统是IIR系统还是FIR系统。,系统的分类,无限长单位脉冲响应（IIR）系统,可见输出不但与输入

39、有关，还与以前的输出及其加权值有关，即系统中存在着输出对输入的反馈回路。这种结构常被称作为递归结构，在求解差分方程时需采用迭代的方法。,有限长单位脉冲响应（FIR）系统,系统的分类,对于FIR系统，它的h(n)为有限长，若已知输入x(n)，可通过卷积直接算出输出y(n)。例如假定h(n)取值范围为0nN-1 则:,最小相位系统,全通系统与最小相位系统,定义:零极点都在单位圆内的系统称为最小相位系统；系统函数的零点都在单位圆外的系统称为最大相位系统；单位圆内外都有零点的系统称为混合相位系统。,设系统具有M个零点，单位圆内有mi个，单位圆外有mo个，有N个极点，单位圆内有ni个，单位圆外有no个。对稳定系统no=0，N=ni。当频率由零变到2时，稳定系统的相位改变量为：,当系统函数的所有零点都在单位圆内时mo=0，当由0变到2时：最小相位系统的相位改变量为0。,全通系统与最小相位系统,全通系统,信号通过全通系统时只有相位上的变化，而无幅度上的变化，因此全通系统是一个纯相位调节系统。,复习,3.1 序列的傅里叶变换3.2 序列的Z变换 3.3 Z变换的基本性质和定理3.4 逆Z变换3.5 Z变换、傅里叶变换、拉普拉斯变换的关系3.6系统函数与频率响应,