数字信号处理第二章z变换与离散时间傅里叶变换DTFT课件.ppt

1、第二章 z变换和DTFT本章主要内容：本章主要内容：1、z变换的定义及收敛域变换的定义及收敛域 2、z变换的反变换变换的反变换 3、z变换的基本性质和定理变换的基本性质和定理 4、离散信号的、离散信号的DTFT 5、z变换与变换与DTFT的关系的关系 6、离散系统的、离散系统的z变换法描述变换法描述2.1 z变换的定义及收敛域变换的定义及收敛域 信号和系统的分析方法有两种：信号和系统的分析方法有两种：时域分析方法时域分析方法变换域分析方法变换域分析方法连续时间信号与系统连续时间信号与系统 LT FT离散时间信号与系统离散时间信号与系统 ZT FT 一、一、ZT的定义的定义),(:),()(21

2、zzXnxnnznxzX)()(z 是复变量，所在的复平面称为是复变量，所在的复平面称为z平面平面 二、二、ZT的收敛域的收敛域 对于任意给定序列对于任意给定序列x(n)，使其，使其z变换变换X(z)收敛的所有收敛的所有z值的集合称为值的集合称为X(z)的收敛域。的收敛域。级数收敛的充要条件是满足绝对可和级数收敛的充要条件是满足绝对可和()nnx n zM 1）有限长序列）有限长序列12()()0 x nnnnx nn其它21Z ()()nnn nX zx n z其 变换：Re zIm jz00Rocz 至少为：Roc-Region of convergence 除除0和和两点是否收敛与两点是

3、否收敛与n1和和n2取值情况取值情况有关外，整个有关外，整个z 平面均收敛。平面均收敛。11(1)111()()(1)(1)nnX zx n zx nzxz22(1)0122(0)(1)(1)()nnxzxzx nzx n zz0 0,021时，nnz0 0,021时，nnz0 0,021时，nn 如果如果n20，则收敛域不包括，则收敛域不包括点点 如果如果n10，则收敛域不包括，则收敛域不包括0点点 如果如果n10n2，收敛域不包括，收敛域不包括0、点点2）右边序列）右边序列11()()0 x nnnx nnn110:0:xxnRoc RznRoc Rz 当时，当时，Re zIm jz0 x

4、Rz 包括处10n 因果序列因果序列的的z变换必在变换必在处收敛处收敛在在处收敛的处收敛的z变换，变换，其序列必为其序列必为因果序列因果序列3）左边序列）左边序列220()()nnx nx nnn220:00:0 xxnRoczRnRoczR当时，当时，Re zIm jz0 xR20n 4）双边序列）双边序列n为任意值时皆有值:xxxxxxRRRocRRRoc RzR当时，当时，Re zIm jz0 xRxR10z()()()nnnnX zx n zx n z其 变换：Roc:0 xzR前式Roc:xRz 后式例例1znZT0,1收敛域应是整个收敛域应是整个z的闭平面的闭平面1 nnzn例例2

5、：求：求x(n)=RN(n)的的z变换及其收敛域变换及其收敛域Re zIm jz0X(z)=()=()nnNnnx n zRn z解：10=Nnnz2 1,.,1rjNzerN零点：01zN极点：（）阶:0Rocz 122111nnnnn nqqqq111Nzz21nq 时须满足11(1)NNzzz例例3：求：求x(n)=anu(n)的变换及其收敛域的变换及其收敛域Re zIm jz0a0X(z)=()=()=nnnnnnnnx n za u n za z解：0z 零点：za极点：:Rocza111 az11az当时Re zIm jz0aX(z)=()=(1)nnnnnx n za unz 解

6、：0z 零点：za极点：:Rocza111111a za zaz11a z当时11=nnnnnna zaz例例4：求：求x(n)=-anu(-n-1)的变换及其收敛域的变换及其收敛域10X(z)=()=nnnnnnnnnnnx n za za za z解：10=nnnnnna za z11nnnaza zaz11/azza 1011nnna zaz11azza例例5：求：求x(n)=a|n|，a为实数，求为实数，求ZT及其收敛域及其收敛域Re zIm jz0a1/a211(1)1()11(1)()azzaaX zazazaz za当时，0,z 零点：1,za a极点：:1/Rocaza1X()

7、az当时，无公共收敛域，不存在 给定给定z变换变换X(z)不能唯一地确定一个序列，不能唯一地确定一个序列，只有同时给出收敛域才能唯一确定。只有同时给出收敛域才能唯一确定。X(z)在收敛域内解析，不能有极点，故：在收敛域内解析，不能有极点，故：右边序列的的z变换收敛域一定在变换收敛域一定在模最模最大的有限极点所在圆的有限极点所在圆之外左边序列的的z变换收敛域一定在变换收敛域一定在模最模最小的有限极点所在圆的有限极点所在圆之内Re zIm jz0abcRe zIm jz0abcRe zIm jz0abcRe zIm jz0abc2.2 z反变换反变换 实质：求实质：求X(z)幂级数展开式幂级数展开

8、式 z反变换的求解方法：反变换的求解方法：围线积分法（留数法）围线积分法（留数法）部分分式法部分分式法 长除法长除法()()x nIZT X zz反变换反变换:从从X(z)中还原出原序列中还原出原序列x(n)()()()nnX zZT x nx n z1 1、围数积分法求解（留数法）围数积分法求解（留数法）若函数若函数X(z)zn-1在围数在围数C上连续，在上连续，在C以内有以内有K个极点个极点zk，而在，而在C以外有以外有M个极点个极点zm，则有：，则有：mzznkzzncnmkzzXsorzzXsdzzzXjnx)(Re)(Re)(21)(111Re zIm jz0 xRxRCRe ()(

9、)()rrz zrz zs F zzz F z1 1、围数积分法求解（留数法）围数积分法求解（留数法）根据复变函数理论，若函数根据复变函数理论，若函数X(z)在环状区域在环状区域 内是解析的，则内是解析的，则在此区域内在此区域内X(z)可展开成罗朗级数，即可展开成罗朗级数，即而而 其中围线其中围线c是在是在X(z)的环状的环状收敛域内环绕原点的一条收敛域内环绕原点的一条反时针方向的闭合单围线。反时针方向的闭合单围线。,0,xxxxRzRRR （）()nnxxnX zC zRzR11()2nncCX z zdzj Re zIm jz0 xRxRC0,1,2,n 若若F(z)在在c外外M个极点个极

10、点zm，且分母多项式，且分母多项式z的的阶次比分子多项式高二阶或二阶以上，则：阶次比分子多项式高二阶或二阶以上，则：11()()(,)2nxxcx nX z zdzcRRj 1()()nF zX z z()Re ()kz zkx ns F z()Re ()mz zmx ns F z 利用留数定理求围线积分，令利用留数定理求围线积分，令 若若F(z)在围线在围线c上连续，在上连续，在c内有内有K个极点个极点zk，则：，则：Re ()()()rrz zrz zs F zzz F z单阶极点的留数：单阶极点的留数：2()1/44(4)(1/4)zX zzzz例1：，求其z反变换Re zIm jz0C

11、41/4211()(,)2(4)(1/4)nxxczx nzdzcRRjzz 解：211()(4)(1/4)(4)(1/4)nnzzF zzzzzz其中：11()4nF zcz 当时在围线 内只有一阶极点14()Re ()zx ns F z1141()4(4)(1/4)nzzzzz415n11()(1)04nF zcznz 当时在围线 内有一阶极点和-阶极点4()Re ()zx ns F z 14441/4nzzzzz 2415ncz=4F(z)而围线 外只有一阶极点，且的分母多项式阶次高于分子多项式阶次两次以上244()(1)(2)1515nnx nu nun Re zIm jz0C41/4

12、2()4(4)(1/4)zX zzzz例2：，求其z反变换Re zIm jz0C41/4解：收敛域是圆的外部 lim()1X(z)z=zX z 又，即在处收敛()()00 x nx nn是一个因果序列，即，()x n是右边序列10()c(4)(1/4)0()0nznF zzzx n同样当时，由在 外无极点，且分母阶次比分子阶次高两阶以上，由围线外极点留数为 可得0n 当时1()(4)(1/4)nzF zzz144cz在围线 内有一阶极点，Re zIm jz0C41/441/4()Re ()Re ()zzx ns F zs F z111441(4)()114(4)()(4)()44nnzzzzz

13、zzzzz21(44)15nn21()(44)()15nnx nu n思考：n=0,1时，F(z)在围线c外也无极点，为何()0 x n 2 2、部分分式展开法求解、部分分式展开法求解IZTIZT：NMnrNkrkkikkknnzzCzzAzBzAzBzX01111)1(1)()()(常见序列的常见序列的ZT参见书参见书p.54页的表页的表2-1若函数若函数X(z)是是z的有理分式，可表示为：的有理分式，可表示为：利用部分分式的利用部分分式的z反变换和可以得到函数反变换和可以得到函数X(z)的的z反变换。反变换。()Re1,2,kkz zX zAskNrz用留数定理求系数：1125()2316

14、zX zzzz例：，求z反变换Re zIm jz032 23353123zzX zAReszzzz 112255516623zzzX zzzzzzz解：1252323X zAAzzzzz 12252123zzX zAReszzzz 1123X zzzz 1111231 21 3zzX zzzzz23z11()1nZT a u nzaaz11(1)1nZT a unzaaz 1112z2()nu n2z 111 3z3(1)nun 3z 231nnx nu nun 例例2 2 设设利用部分分式法求利用部分分式法求z z反变换。反变换。2|,)5.01)(21(1)(11zzzzX5.031234

15、)5.0)(2()(2zzzzzzzzX)()5.0(31234)(nunxnn解：解：3 3、幂级数展开法求解（长除法）、幂级数展开法求解（长除法）:一般一般X(z)是有理分式，可利用分子多项式除是有理分式，可利用分子多项式除分母多项式（长除法法）得到幂级数展开式，分母多项式（长除法法）得到幂级数展开式，从而得到从而得到x(n)。nnzxxzxznxzX1)1()0()1()()(根据收敛域判断根据收敛域判断x(n)的性质，在展开成相应的性质，在展开成相应的的z的幂级数的幂级数 将将X(z)X(z)的的 x(n)展成展成z的的 分子分母分子分母 按按z的的 因果序列因果序列 负幂级数负幂级数

16、 降幂排列降幂排列 左边序列左边序列 正幂级数正幂级数 升幂排列升幂排列xzRxzR例例1 1111 azzX)(az ROC1:)11 az111 az1 az221 zaaz22 za.2211zaaz111 az.2211zaaz,.,21aanx 长除法示例长除法示例解：由解：由RocRoc判定判定x(n)x(n)是因果序列，用长是因果序列，用长除法展成除法展成z z的负幂的负幂级数级数az ROC2:0,.,12aanx111 az.221zaza)11 az1za1 221zaaz 22za.221zazaza11 解：由解：由RocRoc判定判定x(n)x(n)是左边序列，是左边

17、序列，用长除法展成用长除法展成z z的正幂级数的正幂级数2()1/44(4)(1/4)zX zzz z例：，求z反变换解：解：X(z)的的Roc为环状，故为环状，故x(n)是双边序列是双边序列 极点极点z=1/4对应右边序列，极点对应右边序列，极点z=4对应左边序列对应左边序列 先把先把X(z)展成部分分式展成部分分式161()1515(4)()41/41/4X zzzzzzz116()151/44zzX zzz22233416164 44 zzzzzzzz 23144zzz1114114161 141 146 zzzzz 12111416zz2123111()141544X zzzzzz 1

18、+16244()()(1)1515nnx nu nun201114154nnnnnnzz1 1、线性性、线性性2.3 Z变换的基本性质和定理变换的基本性质和定理)()()()(zbYzaXnbynax)()(zXzNnxN)()(azXnxan)()(zXdzdznnxR1R2R|a|RR2 2、序列的移位、序列的移位3 3、z z域尺度变换域尺度变换 （乘以指数序列）（乘以指数序列）4 4、z z域求导域求导 （序列线性加权）（序列线性加权）Z变换的基本性质（续）变换的基本性质（续）)(lim)0(zXxz)()1(lim)(1zXzxz)1()(zXnx5 5、翻褶序列、翻褶序列)()(z

19、Xnx1/RR6 6、共轭序、共轭序列列7 7、初值定理、初值定理8 8、终值定理、终值定理Z变换的基本性质（续变换的基本性质（续）)()()()(zYzXnynx9 9、有限项累加特性、有限项累加特性nmzXzzmxny0)(1)()(dvvHvzXjnhnxc)()(21)()(dvvvHvXjnhnxcn1)1()(21)()(ZTZT的主要性质参见书的主要性质参见书p.69p.69页的表页的表2-22-21010、序列的卷积和、序列的卷积和1111、序列乘法、序列乘法1212、帕塞瓦定理、帕塞瓦定理1LSI ()()(1)()()nnnh nb u nabu nx na u n例：已知

20、系统的单位抽样响应：，求系统输入的响应。()()()nzX zZT x nZT a u nzaza解：1()()()(1)nnH zZT h nZT b u nabu n1()(1)nnZT b u naZT bu n1 zzzaazzbzbzbzb()()()zY zX z H zzbzb()()*()()()ny nx nh nIZT Y zb u nRe zIm jz0ba2.4 序列序列ZT、连续信号、连续信号LT和和FT的关系的关系若：若：)()()()(jXtxsXtxaFTaaLTannsTastaaLTnaaenTxdtetxsXnTtnTxtx)()()()()()(连续信号

21、采样后的拉氏变换连续信号采样后的拉氏变换LT抽样序列：抽样序列：)()(nTxnxannznxzX)()(sTez 当当)()(|)(sXeXzXasTezsT两变换之间的关系，就是由复变量两变换之间的关系，就是由复变量s s平面到复平面到复变量变量z z平面的映射，其映射关系为平面的映射，其映射关系为zTsezsTln1,对比：对比：nnsTaaenTxsX)()(j js sj je ez z进一步讨论这一映射关系：进一步讨论这一映射关系：TereeereTTjTTjj,)(1sTez s平面到平面到z平面的平面的映射是映射是多值映射。T 辐射线辐射线=0 0T T平行直线平行直线=0 0

22、正实轴正实轴=0实轴实轴=0Z平面平面S平面平面:/TT:3/TT/3/TT:)()()1(sXzXa与kaksaezksaakTjsXTjksXTzXjksXTsXsT)2(1)(1|)()(1)()()()2(jXzXa与kaaTjezkTjjXTjXeXzXTj)2(1)()(|)(抽样序列在单位圆上的抽样序列在单位圆上的z z变换，就等于其理想抽样变换，就等于其理想抽样信号的傅里叶变换信号的傅里叶变换数字频率数字频率 表示表示z z平面的辐角，它和模拟角频率平面的辐角，它和模拟角频率 的的关系为关系为jez 在以后的讨论中，将用数字频率在以后的讨论中，将用数字频率 来作为来作为z z平

23、面上平面上单位圆的参数，即单位圆的参数，即ssfffT2所以说，数字频率是模拟角频率的归一化值，或所以说，数字频率是模拟角频率的归一化值，或是模拟频率对抽样频率的相对比值乘以是模拟频率对抽样频率的相对比值乘以2 2 kajezTkjXTeXzXj)2(1)(|)(2.5 离散信号的付氏变换离散信号的付氏变换DTFT一、一、DTFT的定义的定义变换对：变换对：)()(jDTFTeXnx njnjenxeX)()(deeXnxjnj)(21)(称为称为离散时间傅里叶变换（离散时间傅里叶变换（DTFT）。）。FT存在的充分必要条件是：存在的充分必要条件是：)()(nxenxjwn如果引入冲激函数，一

24、些绝对不可和的序列，如果引入冲激函数，一些绝对不可和的序列，如周期序列，其傅里叶变换可用冲激函数的如周期序列，其傅里叶变换可用冲激函数的形式表示出来。形式表示出来。二、比较二、比较ZT和和DTFT的定义：的定义：dweeXdzzzXjnxenxzXeXjwnjznnjwnezjj)(21)(21)()(|)()(1|1 利用利用ZT和和DTFT的关系可以有的关系可以有ZT计算计算DTFT。序列的傅里叶变换是序列的序列的傅里叶变换是序列的z变换在单位圆变换在单位圆上的值上的值notherwise,Nn,(n)Rx(n)N010111011zzzX(z)NNnnjjNjeeeX11)(NjjjjN

25、jNjNjjeNeeeeeeeX2)1(222222)2sin()2sin()()()(例例1、计算门序列的、计算门序列的DTFT)2sin()2sin()(NeXj2)1()(N (类似类似Sa(.)Sa(.)函数函数 )(线性相位线性相位)解：解：DTFT幅频特性：幅频特性：相频特性：相频特性：图示说明：图示说明：零极点图(N=8)Z平面11jj)(nRNn0N-11)(X022N=8N例例2 2、已知、已知 (),(),计算其计算其DTFTDTFT。,)(:)()(nuanfn1a)sin)cos1(111)(0jaaeeaeFjnjnnjFTDTZTazzFeFjjezezj111)(

26、)(由此可以得到由此可以得到FT的的幅频特性幅频特性和和相频特性相频特性aaeFjcos211)(2)cos1sin()(1aatg物理说明物理说明:若若 (语音信号处理中常用该指数语音信号处理中常用该指数 函数展宽单音信号的频谱函数展宽单音信号的频谱),),该信号该信号3db3db带宽带宽 (或或 )。具体求。具体求 解过程如下：解过程如下：令令 即即 可解出可解出kHz,f.as89940Hzf15rad.cc006021)F(e)F(ejjC021a)(aaC121cos2112rad.C0060Hz ffsc152)(nf)(jeFn022.a1121c三、三、FT与与DTFT的关系的

27、关系kaTajTkjXTjXeX)2(1|)()(kajkjXeX)2()(归一化归一化 利用利用FT与与DTFT关系计算下列序列的关系计算下列序列的 DTFT 1)()cos()()(30210nx;nnx;enxnj例：例：解：解：1）)(2)()(0110jXetxFTtj DTFTnjenx0)(1mjmeX)2(2)(01)()(cos000tFT)2()2(cos)(0002 mDTFTmmnnx)(2)(1)(33jXtxFT DTFTnx)(3mjmeX)2(2)(32）3）2.6 DTFT的一些性质的一些性质)()()()(22112211jjeFaeFanfanfa)()(

28、*jjeXeX1 1、线性性：、线性性：)(Re2)()()(jeeXnxnxnx)(Im2)()()(joeXjnxnxnx0)()(0jnjeeXnnx2 2、实序列：、实序列：实偶性：实偶性：实奇性：实奇性：3 3、时移特性：、时移特性：)()()(00jnjeXnxe)()(jeXddjnnx4 4、乘以指数序列、乘以指数序列 （调制性）（调制性）5 5、序列线性加权、序列线性加权)()(jeXnx6 6、序列翻褶、序列翻褶)()(jeXnx7 7、序列共、序列共轭轭8 8、卷积定理：、卷积定理：(时域时域)(频域频域)()()()(jjeYeXnynxdeYeXeYeXnynxjjj

29、j)()(21)()()()()(jnd)X(e(n)x2221deYeXnynxjjn)()(21)()(*DTFT的主要性质参见书的主要性质参见书p.78页的表页的表2-39 9、帕色伐尔定理：、帕色伐尔定理：(Parseval Theory)频域卷积在一周期内积分频域卷积在一周期内积分,称称周期卷积周期卷积。下面举例说明下面举例说明DTFT性质得使用。性质得使用。计算下列积分计算下列积分I的值。的值。jjd)be)(ae(I11111b,ajnjnbeu(n)baeu(n)a1111解：根据解：根据 debeaenubnuajnjjnn)1)(1(121)()(222200000ba|b

30、au(n)|bu(n)aInnmmnmnnn利用时域卷积定理有：利用时域卷积定理有：上式卷积上式卷积n=0时就是积分时就是积分I的值。的值。2.7 周期性序列的周期性序列的DTFT1、复指数序列的傅里叶变换、复指数序列的傅里叶变换00),2(20injieq复指数序列复指数序列ej 0n的傅里叶变换，是以的傅里叶变换，是以 0为中心，为中心，以以2 的整数倍为间距的一系列冲激函数，其积分的整数倍为间距的一系列冲激函数，其积分面积为面积为2 q思考，思考，DTFTcos(0n+ff、DTFT sin(0n+ff2、常数序列的傅里叶变换、常数序列的傅里叶变换iiiin)2(2)(1q常数序列的傅里

31、叶变换，是以常数序列的傅里叶变换，是以0为中心，以为中心，以2 的整的整数倍为间距的一系列冲激函数，其积分面积为数倍为间距的一系列冲激函数，其积分面积为2 3、周期为、周期为N的抽样序列串的傅里叶变换的抽样序列串的傅里叶变换kikNNiNn)2(2)(q周期为周期为N的周期性抽样序列，其傅里叶变换是频的周期性抽样序列，其傅里叶变换是频率在率在2/N的整数倍上的的整数倍上的一系列冲激函数之和，一系列冲激函数之和，这些冲激函数的积分面积为这些冲激函数的积分面积为2/N/N4、一般性的周期为、一般性的周期为N的周期性序列的傅里叶变换的周期性序列的傅里叶变换kkkNjkjkijkNkXNkNeXNkN

32、NeXnxkNNiNneXnx)2()(2)2()(2)2(2)()()2(2)()()(2iiiNnnxiNnxnx)()()()(1021021022)()()()()(NnnkNjNnnkNjNnkNnjkNjenxenxenxeXkXq周期性序列周期性序列 （周期为（周期为N）的傅里叶变换是）的傅里叶变换是一系一系列冲激函数串，其冲激函数的积分面积等于列冲激函数串，其冲激函数的积分面积等于 乘以乘以，而，而 是是x(n)的一个周期的一个周期的傅里叶的傅里叶变换变换X(ej)在频域中在频域中 2/N/N的整数倍的各抽样点的整数倍的各抽样点上的抽样值。上的抽样值。)(nx)(kX)(kX)

33、(nxq即：即：kkNkXNnxDTFT)2()(2)(1021020201020)(1)2()(1)2()(1)2()(221)(NkknNjNknjnjNknjkekXNdekNkXNdekNkXNdekNkXNnx 满足满足0 0 2/N从从00之前开始抽样；之前开始抽样；在在22之间结束抽样；之间结束抽样；此区间共有此区间共有N N个抽样值：个抽样值：0 0 k N1N1周期序列的周期序列的DFS正变换和反变换正变换和反变换21100()()()()NNjnknkNNnnX kDFS x nx n ex n W2110011()()()()NNjnknkNNkkx nIDFS X kX

34、 k eX k WNN周期序列的傅里叶级数（周期序列的傅里叶级数（DFS）2jNNWe其中：其中：2.8 Fourier变换的对称性质变换的对称性质共轭对称序列：共轭对称序列：*()()eex nxn*()()oox nxn()()()eox nx nx n共轭反对称序列：共轭反对称序列：任意序列可表示成任意序列可表示成xe(n)和和xo(n)之和之和:*1()()()2ex nx nxn*1()()()2ox nx nxn其中：其中：定义：定义：*1()()()()2jjjjeeX eXeX eXe*1()()()()2jjjjooX eX eX eX e其中：其中：()()()jjjeoX

35、 eXeXe()jX e同样，同样，x(n)的的Fourier变换变换 也可分解成：也可分解成：对称性质对称性质 序列序列 Fourier变换变换()()jx nX eRe()()jex nXeIm()()jojx nXe()Re()jex nX e()Im()jox njX e实数序列的对称性质实数序列的对称性质 序列序列 Fourier变换变换Re()()()jjex nX eX eIm()0()0jojx nX e()Re()jex nX e()Im()jox njX e*()()()jjjeX eX eX e实数序列的实数序列的Fourier变换满足共轭对称性变换满足共轭对称性Re()

36、Re()jjX eX eIm()Im()jjX eX e 实部是实部是的偶函数的偶函数虚部是虚部是的奇函数的奇函数()()jjX eX earg()arg()jjX eX e 幅度是幅度是的偶函数的偶函数幅角是幅角是的奇函数的奇函数2.9 离散系统的系统函数、系统的频率响应离散系统的系统函数、系统的频率响应LSI系统的系统的系统函数系统函数H(z)：单位抽样响应单位抽样响应h(n)的的z变换变换()()()()()nnY zH zZT h nh n zX z其中：其中：y(n)=x(n)*h(n)Y(z)=X(z)H(z)系统的系统的频率响应频率响应 ：()jH e()()()jjz eH e

37、H zDTFT h n 单位圆上的系统函数单位圆上的系统函数,单位抽样响应单位抽样响应h(n)的的DTFT1、若、若LSI系统为因果稳定系统系统为因果稳定系统稳定系统的系统函数稳定系统的系统函数H(z)的的Roc须包含单位圆，须包含单位圆，即频率响应存在且连续即频率响应存在且连续H(z)须从单位圆到须从单位圆到的整个的整个z域内收敛即系统域内收敛即系统函数函数H(z)的全部极点必须在单位圆内的全部极点必须在单位圆内xRz1 1）因果：）因果：2 2）稳定：）稳定：()nh n 序列序列h(n)绝对可和，即绝对可和，即()nnh n z 而而h(n)的的z变换的变换的Roc：1z 3 3）因果稳

38、定：）因果稳定：RocRoc：/4/4/6/60.2,0.2,0.4,2,2,1.5jjjjeeee例：一LSI系统的极点有：问什么情况下，系统为因果系统，什么情况下，系统为稳定系统Re zIm jz0140.2je40.2je0.41.562je62je2z 解：因果系统：0.41.5z稳定系统：2、系统函数与差分方程、系统函数与差分方程常系数线性差分方程：常系数线性差分方程：00()()NMkmkma y nkb x nm00()()NMkmkmkma z Y zb zX z101101(1)()()/()(1)MMmmmmmNNkkkkkb zc zH zY zX zKa zd z取取z

39、变换变换则系统函数则系统函数LSI311()(1)(2)()(1)483()()123y ny ny nx nx nx ny n例：已知离散系统的差分方程：（设系统初始状态为零)其中：为输入，为输出。）求系统函数，指出系统的零极点；）若该系统是因果稳定的，指出系统的收敛域；）求该因果稳定系统的单位抽样响应。z解：1）对差分方程两边取 变换：121311()()()()()483Y zz Y zz Y zX zz X z1112111111()33()3111()1114824zzY zH zX zzzzz111,0 ,324zz 零点：极点：系统函数：212z）由于系统为因果稳定系统，故收敛域

40、：Re zIm jz00.50.2511/3 111131131111241124zzH zzzzzz 121311112424zH zAAzzzzz 1121211103112324zzzH zAReszzzz3)H(z)h(n)对求z反变换即得单位抽样响应，用部分分式法 214141173114324zzzH zAReszzzz 10733()1124zzH zzz1:2-12Rocz 根据，查表得 10 17 1()323 4nnh nu n3、系统的频率响应的意义、系统的频率响应的意义1）LSI系统对复指数序列的稳态响应：系统对复指数序列的稳态响应：0()jnx nen 000()()

41、()()jn mjnjmmmy nh m eeh m e00()jnjeH e0()cos()x nAnf000()()cosarg()jjy nA H enH ef2）LSI系统对正弦序列的稳态响应系统对正弦序列的稳态响应输出同频输出同频 正弦序列正弦序列幅度受频率响应幅度幅度受频率响应幅度 加权加权相位为输入相位与系统相位响应之和相位为输入相位与系统相位响应之和()jH e03）LSI系统对任意输入序列的稳态响应系统对任意输入序列的稳态响应 ()()*()y nx nh n()()()jjjY eX eH e1()()()2jjj ny nH eX eed1()()2jj nx nX ee

42、d其中：其中：1()2jj nX eed微分增量（复指数）：微分增量（复指数）：4、频率响应的几何确定法、频率响应的几何确定法利用利用H(z)在在z平面上的零极点分布平面上的零极点分布1()11111(1)()()(1)()MMmmNMmmNNkkkkc zzcH zKKzd zzd()arg()11()()()()jMjmjj NMjjH emNjkkecH eKeH eeed频率响应：频率响应：则频率响应的则频率响应的mjjmmmcece kjjkkkdedl ef 11arg()arg()MNjmkmkH eKNMf令令幅角：幅角：11()MmjmNkkH eKl幅度：幅度：零点位置影响

43、凹谷点的位置与深度零点位置影响凹谷点的位置与深度 零点在单位圆上，谷点为零零点在单位圆上，谷点为零 零点趋向于单位圆，谷点趋向于零零点趋向于单位圆，谷点趋向于零 极点位置影响凸峰的位置和深度极点位置影响凸峰的位置和深度 极点趋向于单位圆，峰值趋向于无穷极点趋向于单位圆，峰值趋向于无穷 极点在单位圆外，系统不稳定极点在单位圆外，系统不稳定()()nh na u n ()()(1)1y nx nay naa例：设一阶系统的差分方程：，为实数 求系统的频率响应。1z()1()()1Y zH zzaX zaz解：两边求 变换，得 11()1(1cos)sinjjH eaeaja21/2()(12 co

44、s)sinarg()arctan1cosjjH eaaaH ea 幅度响应：相位响应：2()()(1)(2).y nx nax na x n例：设系统的差分方程(0a1)：110(1)()MMkkax nMa x nk1MM这就是个单元延时及个抽头加权后相加所组成的电路，常称之为横向滤波器，求其频率响应。21,2,.,1jiMizaeiM零点：，1110()()z1()01()MMMMMkkMkx nna zzaH za zzazzza解：令，两边取 变换()()01()0nnh nanMh nn当输入为，则输出为其它0(1)zMza极点：，阶，处零极点相消5、IIR系统和系统和FIR系统系统

45、无限长单位冲激响应（无限长单位冲激响应（IIR）系统：）系统：单位冲激响应单位冲激响应h(n)是无限长序列是无限长序列有限长单位冲激响应（有限长单位冲激响应（FIR）系统：）系统：单位冲激响应单位冲激响应h(n)是有限长序列是有限长序列0001()1MMmmmmmmNNkkkkkkb zb zH za za z0ka IIR系统：至少有一个系统：至少有一个0ka FIR系统：全部系统：全部0b全极点系统全极点系统（自回归系统，自回归系统，AR系统系统）：分子只有常数项：分子只有常数项0b零极点系统零极点系统（自回归滑动平均系统，自回归滑动平均系统，ARMA系统系统）:分子不止常数项分子不止常数项收敛域收敛域 内无极点，是全零点系统内无极点，是全零点系统0z（滑动平均系统，（滑动平均系统，MA系统）系统）00()()()MNmkmky nb x nma y nk0ka IIR系统：至少有一个系统：至少有一个有反馈环路，采用递归型结构有反馈环路，采用递归型结构0ka FIR系统：全部系统：全部无反馈环路，多采用非递归结构无反馈环路，多采用非递归结构Homework：P831(1)(2)(3)3(1)7 10 14 18