第二章z变换与离散时间傅里叶变换(DTFT)

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

2、(21zzXnxnnznxzX)()( 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其 变换：0Rocz 至少为： Re zIm jz0 除除0和和两点是否收敛与两点是否收敛与n1和和n2取值情况

3、取值情况有关外，整个有关外，整个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 xRz

4、包括处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的闭平面的闭平面

5、1 nnzn例例2：求：求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)

6、nnnnnx n za unz 解：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 zazazazza当时，0,z 零点：1,z

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

8、质：求实质：求X(z)幂级数展开式幂级数展开式 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)(R

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

10、Im jz0 xRxRC0, 1, 2,n 若若F(z)在在c外外M个极点个极点zm，且分母多项式，且分母多项式z的的阶次比分子多项式高二阶或二阶以上，则：阶次比分子多项式高二阶或二阶以上，则：11( )( )(,)2nxxcx nX z zdzcRRj1( )( )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单阶极点的留数：单阶

11、极点的留数：2( ) 1/44(4)(1/4)zX zzzz例1：，求其z反变换Re zIm jz0C41/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)而围线 外只有一阶极点，且的分母

12、多项式阶次高于分子多项式阶次两次以上244( )(1)(2)1515nnx nu nun Re zIm jz0C41/42( ) 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 j

13、z0C41/441/4( )Re ( )Re ( )zzx ns F zs F z111441(4)()114(4)()(4)()44nnzzzzzzzzzz21(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的有理分式，可表示为：的有理分式，可表示为： 利用部分分式的利

14、用部分分式的z反变换和可以得到函数反变换和可以得到函数X(z) 的的z反变换。反变换。( )Re1,2,kkz zX zAskNrz用留数定理求系数：1125( ) 2316zX zzzz例：，求z反变换Re zIm jz032 23353123zzX zAReszzzz 112255516623zzzX zzzzzzz解： 1252323X zAAzzzzz 12252123zzX zAReszzzz 1123X zzzz 111123121 3zzX zzzzz23z11( )1nZT a u nzaaz11(1)1nZT a unzaaz1112z2( )nu n2z 111 3z3(1

15、)nun 3z 231nnx nu nun 例例2 2 设设利用部分分式法求利用部分分式法求z z反变换。反变换。2|,)5 . 01)(21 (1)(11zzzzX5 . 031234)5 . 0)(2()(2zzzzzzzzX)()5 . 0(31234)(nunxnn解：解：3 3、幂级数展开法求解（长除法）、幂级数展开法求解（长除法）: : 一般一般X(z)是有理分式，可利用分子多项式除是有理分式，可利用分子多项式除分母多项式（长除法法）得到幂级数展开式，分母多项式（长除法法）得到幂级数展开式，从而得到从而得到x(n)。nnzxxzxznxzX1) 1 ()0() 1()()( 根据收

16、敛域判断根据收敛域判断x(n)的性质，在展开成相应的性质，在展开成相应的的z的幂级数的幂级数 将将X(z) X(z)的的 x(n) 展成展成z的的 分子分母分子分母 按按z的的 因果序列因果序列 负幂级数负幂级数 降幂排列降幂排列 左边序列左边序列 正幂级数正幂级数 升幂排列升幂排列xzRxzR例例1 1111 azzX)(az ROC1:)11 az111 az1 az221 zaaz22 za. 2211zaaz111 az. 2211zaaz,.,21aanx 长除法示例长除法示例解：由解：由RocRoc判定判定x(n)x(n)是因果序列，用长是因果序列，用长除法展成除法展成z z的负幂

17、的负幂级数级数az ROC2:0 ,.,12aanx111 az. 221zaza)11 az1za1 221zaaz 22za. 221zazaza11 解：由解：由RocRoc判定判定x(n)x(n)是左边序列，是左边序列，用长除法展成用长除法展成z z的正幂级数的正幂级数2( ) 1/44(4)(1/4)zX zzzz例：，求z反变换解：解：X(z)的的Roc为环状，故为环状，故x(n)是双边序列是双边序列 极点极点z=1/4对应右边序列，极点对应右边序列，极点z=4对应左边序列对应左边序列 先把先把X(z)展成部分分式展成部分分式161( )1515(4)()41/41/4X zzzz

18、zzz116( )151/44zzX zzz22233416164 44 zzzzzzzz 23144zzz1114114161 141 146 zzzzz 12111416zz2123111( )141544X zzzzzz 1+16244( )( )(1)1515nnx nu nun201114154nnnnnnzz1 1、线性性、线性性2.3 Z变换的基本性质和定理变换的基本性质和定理)()()()(zbYzaXnbynax)()(zXzNnxN)()(azXnxan)()(zXdzdznnxR1R2R|a|RR2 2、序列的移位、序列的移位3 3、z z域尺度变换域尺度变换 （乘以指数

19、序列）（乘以指数序列）4 4、 z z域求导域求导 （序列线性加权）（序列线性加权）Z变换的基本性质（续）变换的基本性质（续） )(lim)0(zXxz)() 1(lim)(1zXzxz)1()(zXnx5 5、翻褶序列、翻褶序列)()(zXnx1/RR6 6、共轭序、共轭序列列7 7、初值定理、初值定理8 8、终值定理、终值定理Z变换的基本性质（续变换的基本性质（续）)()()()(zYzXnynx9 9、有限项累加特性、有限项累加特性nmzXzzmxny0)(1)()(dvvHvzXjnhnxc)()(21)()(dvvvHvXjnhnxcn1)1()(21)()(ZTZT的主要性质参见书

20、的主要性质参见书p.69p.69页的表页的表2-22-21010、序列的卷积和、序列的卷积和1111、序列乘法、序列乘法1212、帕塞瓦定理、帕塞瓦定理1LSI ( )( )(1) ( )( )nnnh nb u nabu nx na u n例：已知系统的单位抽样响应：，求系统输入的响应。( ) ( )( ) 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( )( )* ( ) ( )

21、( )ny nx nh nIZT Y zb u nRe zIm jz0ba2.4 序列序列ZT、连续信号、连续信号LT和和FT的关系的关系若：若：)()()()(jXtxsXtxaFTaaLTannsTastaaLTnaaenTxdtetxsXnTtnTxtx)()()()()()(连续信号采样后的拉氏变换连续信号采样后的拉氏变换LT抽样序列：抽样序列：)()(nTxnxannznxzX)()(sTez 当当)()(| )(sXeXzXasTezsT两变换之间的关系，就是由复变量两变换之间的关系，就是由复变量s s平面到复平面到复变量变量z z平面的映射，其映射关系为平面的映射，其映射关系为z

22、TsezsTln1,对比：对比：nnsTaaenTxsX)()(j js sj je ez z进一步讨论这一映射关系：进一步讨论这一映射关系：TereeereTTjTTjj,)(1sTez s平面到平面到z平面的平面的映射是映射是多值映射。T 辐射线辐射线= =0 0T T平行直线平行直线 =0 0正实轴正实轴=0实轴实轴 =0Z平面平面S平面平面: :/TT: :3 /TT /3 /TT: : :)()() 1 (sXzXa与kaksaezksaakTjsXTjksXTzXjksXTsXsT)2(1)(1| )()(1)()()()2(jXzXa与kaaTjezkTjjXTjXeXzXTj)

23、2(1)()(| )(抽样序列在单位圆上的抽样序列在单位圆上的z z变换，就等于其理想抽样变换，就等于其理想抽样信号的傅里叶变换信号的傅里叶变换数字频率数字频率 表示表示z z平面的辐角，它和模拟角频率平面的辐角，它和模拟角频率 的的关系为关系为jez 在以后的讨论中，将用数字频率在以后的讨论中，将用数字频率 来作为来作为z z平面上平面上单位圆的参数，即单位圆的参数，即ssfffT2所以说，数字频率是模拟角频率的归一化值，或所以说，数字频率是模拟角频率的归一化值，或是模拟频率对抽样频率的相对比值乘以是模拟频率对抽样频率的相对比值乘以2 2 kajezTkjXTeXzXj)2(1)(| )(2

24、.5 离散信号的付氏变换离散信号的付氏变换DTFT一、一、DTFT的定义的定义变换对：变换对：)()(jDTFTeXnx njnjenxeX)()(deeXnxjnj)(21)(称为称为离散时间傅里叶变换（离散时间傅里叶变换（DTFT）。）。FT存在的充分必要条件是：存在的充分必要条件是：)()(nxenxjwn如果引入冲激函数，一些绝对不可和的序列，如果引入冲激函数，一些绝对不可和的序列，如周期序列，其傅里叶变换可用冲激函数的如周期序列，其傅里叶变换可用冲激函数的形式表示出来。形式表示出来。二、比较二、比较ZT和和DTFT的定义：的定义：dweeXdzzzXjnxenxzXeXjwnjznn

25、jwnezjj)(21)(21)()(| )()(1|1 利用利用ZT和和DTFT的关系可以有的关系可以有ZT计算计算DTFT。 序列的傅里叶变换是序列的序列的傅里叶变换是序列的z变换在单位圆变换在单位圆上的值上的值notherwise,Nn,(n)Rx(n)N010111011zzzX(z)NNnnjjNjeeeX11)(NjjjjNjNjNjjeNeeeeeeeX2)1(222222)2sin()2sin()()()(例例1、计算门序列的、计算门序列的DTFT)2sin()2sin()(NeXj2) 1()(N ( (类似类似Sa(.)Sa(.)函数函数 ) )( (线性相位线性相位) )

26、 解：解：DTFT幅频特性：幅频特性：相频特性：相频特性：图示说明：图示说明：零极点图(N=8)Z平面11jj)(nRNn0N-11)(X022N=8N例例2 2、已知、已知 ( ),( ),计算其计算其DTFTDTFT。,)( :)()(nuanfn1a)sin)cos1 (111)(0jaaeeaeFjnjnnjFTDTZTazzFeFjjezezj111)()(由此可以得到由此可以得到FT的的幅频特性幅频特性和和相频特性相频特性aaeFjcos211)(2)cos1sin()(1aatg物理说明物理说明: : 若若 ( (语音信号处理中常用该指数语音信号处理中常用该指数 函数展宽单音信号

27、的频谱函数展宽单音信号的频谱) ,) ,该信号该信号3db3db带宽带宽 ( (或或 ) )。具体求。具体求 解过程如下：解过程如下： 令令 即即 可解出可解出kHz, f.as89940Hzf15rad.cc006021)F(e)F(ejjC021a)(aaC121cos2112rad.C0060Hz ffsc152)(nf)(jeFn022.a1121c三、三、FT与与DTFT的关系的关系kaTajTkjXTjXeX)2(1| )()(kajkjXeX)2()(归一化归一化 利用利用FT与与DTFT关系计算下列序列的关系计算下列序列的 DTFT 1)()cos()()(30210nx;nn

28、x;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)()(*jjeXeX1 1、线性性：、线性性：)(Re2)()()(jeeXnxnxnx)(Im2)()()(joeXjnxnxnx0)()(0jnjeeXnnx2 2、实序列：、实序列：实偶性：实偶性：实奇

29、性：实奇性：3 3、时移特性：、时移特性：)()()(00jnjeXnxe)()(jeXddjnnx4 4、乘以指数序列、乘以指数序列 （调制性）（调制性）5 5、序列线性加权、序列线性加权)()(jeXnx6 6、序列翻褶、序列翻褶)()(jeXnx7 7、序列共、序列共轭轭8 8、卷积定理：、卷积定理： ( (时域时域) ) ( (频域频域) )()()()(jjeYeXnynxdeYeXeYeXnynxjjjj)()(21)()()()()(jnd)X(e(n)x2221deYeXnynxjjn)()(21)()(*DTFT的主要性质参见书的主要性质参见书p.78页的表页的表2-39 9

30、、帕色伐尔定理：、帕色伐尔定理：(Parseval Theory)频域卷积在一周期内积分频域卷积在一周期内积分, ,称称周期卷积周期卷积。下面举例说明下面举例说明DTFT性质得使用。性质得使用。计算下列积分计算下列积分I的值。的值。jjd)be)(ae(I11111b,ajnjnbeu(n)baeu(n)a1111解：根据解：根据 debeaenubnuajnjjnn)1)(1 (121)()(222200000ba|bau(n)|bu(n)aInnmmnmnnn利用时域卷积定理有：利用时域卷积定理有：上式卷积上式卷积n=0时就是积分时就是积分I的值。的值。2.7 周期性序列的周期性序列的DT

31、FT1、复指数序列的傅里叶变换、复指数序列的傅里叶变换00),2(20injieq复指数序列复指数序列ej 0n的傅里叶变换，是以的傅里叶变换，是以 0为中心，为中心，以以2 的整数倍为间距的一系列冲激函数，其积分的整数倍为间距的一系列冲激函数，其积分面积为面积为2 q思考，思考，DTFTcos( 0n+ff、 DTFT sin( 0n+ff2、常数序列的傅里叶变换、常数序列的傅里叶变换iiiin)2(2)(1q常数序列的傅里叶变换，是以常数序列的傅里叶变换，是以0为中心，以为中心，以2 的整的整数倍为间距的一系列冲激函数，其积分面积为数倍为间距的一系列冲激函数，其积分面积为2 3、周期为、周

32、期为N的抽样序列串的傅里叶变换的抽样序列串的傅里叶变换kikNNiNn)2(2)(q周期为周期为N的周期性抽样序列，其傅里叶变换是频的周期性抽样序列，其傅里叶变换是频率在率在2 /N的整数倍上的的整数倍上的一系列冲激函数之和，一系列冲激函数之和，这些冲激函数的积分面积为这些冲激函数的积分面积为2/N/N4、一般性的周期为、一般性的周期为N的周期性序列的傅里叶变换的周期性序列的傅里叶变换kkkNjkjkijkNkXNkNeXNkNNeXnxkNNiNneXnx)2()(2)2()(2)2(2)()()2(2)()()(2iiiNnnxiNnxnx)()()()(1021021022)()()()

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

34、kekXNdekNkXNdekNkXNdekNkXNnx 满足满足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 k eX k WNN周期序列的傅里叶级数（周期序列的傅里叶级数（DFS）2jNNWe其中：其中：2.8 Fourier变换的对称性质变换的对称

35、性质共轭对称序列：共轭对称序列：*( )()eex nxn*( )()oox nxn ( )( )( )eox nx nx n共轭反对称序列：共轭反对称序列： 任意序列可表示成任意序列可表示成xe(n)和和xo(n)之和之和:*1( ) ( )()2ex nx nxn*1( ) ( )()2ox nx nxn其中：其中：定义：定义：*1()()()()2jjjjeeX eX eX eX e*1()()()()2jjjjooX eX eX eX e其中：其中：()()()jjjeoX eXeXe()jX e同样，同样，x(n)的的Fourier变换变换 也可分解成：也可分解成：对称性质对称性质

36、序列序列 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()Re()jjX eX eIm()Im()jjX eX e 实部是实部是的偶函数的偶

37、函数虚部是虚部是的奇函数的奇函数()()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 eH zDTFT h n 单位圆上的系统函数单位圆上的系统函数

38、,单位抽样响应单位抽样响应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）因果稳定：）因果稳定：RocRoc：/4/4/6/60.2,0

39、.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变换变换则系统函数则系统函数LSI311

40、( )(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 ）由于系统为因果稳定系统， 故收敛域： Re

41、 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(

42、)( )( )( )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 eed其中：其中：1()2jj nX eed微分增量（复指数）：微分增量（复指数）