资源描述
第五章 z变换
5-1 概述
在电路理论中我们学过拉普拉斯变换,在那里把它作为连续时间傅里叶变换的一种推广,做这中推广的部分原因是由于拉普拉斯变换比傅里叶变换有着更广泛的适用范围,有许多信号,其傅里叶变换不存在,但却有拉普拉斯变换,比如一个不稳定的线性时不变系统的傅里叶变换不存在,但它的拉普拉斯变换却存在,运用拉普拉斯变换可以对系统的不稳定性作分析,从而找出使系统稳定的措施或找出系统不稳定的原因。
对于离散系统表述系统和信号的数学抽象是序列,其变量为离散变量,因此拉普拉斯变换已不适用。作为序列的傅里叶变换的推广就是z变换。作为一种重要的数学工具,它把描述离散系统的差分方程,变换成代数方程,使其求解过程得到简化。还可以利用系统函数的零、极点分布,定性分析系统的时域特性、频率响应、稳定性等,是离散系统分析的重要方法。Z变换在离散系统的作用与地位,与拉氏变化在连续时间系统相当。本章中要讨论z变换的定义、性质和它与傅里叶变换、拉普拉斯变换的关系。
5-2 z变换的定义及其收敛域
序列的傅里叶变换定义为:
(5-1)
式中是的复函数,变量是实数,也可以看成是复数变量的函数,这时就是复变函数,只不过“复数”变量只在虚轴上取值,现在若将这个“复数”延拓到实轴上取值,即,这时(5-1)式中的就变成,这样(5-1)式就应该写成
(5-2)
显然,上式表述在形式有点累赘,注意到本身也是一个复数变量,令,则(5-2)式就为
(5-3)
上式就是序列的z变换,通过前面的说明可以看到序列的z变换实际上就是序列傅里叶变换的推广,序列的z变换是复数变量z的函数,即是个复变函数。(5-3)式中对n的求和是从到,即是在时域坐标原点的两边范围内进行,所以(5-3)式定义的z变换通常称为双边z变换。z变换我们也用符号“Z”来表示,即
(5-4)
将(5-3)式展开有
(5-5)
这是一个以z为变量的幂级数,我们知道只有当一个幂级数收敛时,讨论这个幂级数的特性才有意义,也就是说序列的z变换存在,讨论序列的z变换才有意义。那么,序列是否存在?存在的条件是什么?这时我们学习z变换时首先要回答的问题。
序列的z变换是z平面上的一个函数,所谓z平面是指图5-1所示的平面,z平面的横轴表示复数变量z的实部,z平面的纵轴表示复数变量的虚部。在z平面中以坐标原点为圆心,以单位1为半径的圆称为z平面中的单位圆,它在z平面中的方程为。比较序列的傅里叶变换公式
图5-1z平面及单位圆
既有 (5-6)
这就是说,序列的傅里叶变换实际上就是序列z变换在z平面中单位圆上的取值。换句话说,序列的傅里叶变换是序列在z平面单位圆上的变换,z变换是将变换延拓到整个z平面中。在一次回到前面提出的问题,这种延拓是否存在,即(5-5)式所示的幂级数是否收敛。显然,(5-5)式的收敛性取决于离散信号和z的取值范围。如果序列给定,则(5-5)式收敛性取决于z的取值范围,我们称所有使(5-5)式幂级数(即序列的z变换)绝对收敛的z值的集合为序列z变换的收敛域。
例如:
则
这是个等比幂级数级数,为了求出幂级数的闭式解可以用以下极限方法求解,
上示方括号中的和式为等比级数求和,则有
当时,上式极限存在的充分必要条件是,即,也就是说,所有满足关系的z值都可使序列的z变换收敛,所以单位阶跃序列u(n)的收敛域为,对所有收敛域中的z值,阶跃序列z变换的闭式表达式为
5-3 序列z变换的基本特性
在以上的讨论中,我们强调了序列z变换收敛域的重要性,那么是否所有序列的z变换都可以找到收敛域呢?换句话说,是否所有序列的z变换都存在(既有意义)?下面就来讨论这个问题。
将复变量z表示成极坐标形式,则序列的z变换可以写成
因此序列的z变换可以看成是一个实指数序列乘以序列后得傅里叶变换。序列傅里叶变换一致收敛(存在)的条件是序列绝对可积,即
(5-7)
由此可见,由于序列乘上了实指数序列,即使序列的傅里叶变换不存在,但它的z变换却可能存在,这表明z变换适用的范围远比傅里叶变换宽,这也就是要引入z变换分析的原因之一。
前面我们已经提到序列z变换的收敛域与序列本身的特性有关,下面就来讨论序列z变换的基本特性。
1. 有限长序列(finite length frequence)
若序列的非零值点仅分布在有限整数集合上,则称此序列为有限长序列,即
n<n1或n>n2 (5-8)
这里。因为序列的非零值点仅分布在有限整数集合[n1, n2]上,所以这类序列的z变换为
(5-9)
要使这个和式收敛,在序列有界(即)的条件下,z变换的收敛域就取决于,的取值。现在分几种情况来讨论,
第一种情况:,在这种情况下,z变换可以写成
上式中第一个和式里的n<0,所以在处,,和式不收敛(发散),此外对所有z这个和式均收敛。上式中第二个和式里,所以在z=0处,,和式也不收敛(发散),此外对所有z值这个和式均收敛。综上所述,在第一种情况下,有限长序列的z变换收敛域为,即除了z=0和z=∞外,序列z变换在整个z平面上收敛。
第二种情况:,这时序列z变换为
这个和式中的n<0,所以z变换除了在处发散外,取任何z值均收敛,所以在第二种情况下,有限长序列的z变换收敛域为,即除了z=∞外,序列z变换在整个z平面上收敛。
第三种情况:,这时序列z变换为
这个和式中的,所以z变换除了在z=0发散外,取任何z值均收敛,所以在第三种情况下,有限长序列的z变换收敛域为,即除了z=0外,序列z变换在整个z平面上收敛。
第四种情况:,这是一种特殊情况,即,A为常数,序列的z变换为
它是一个常数,所以无论z取何值均收敛,即这时序列的z变换收敛于整个z平面。
2. 右边序列(right-sided frequence)
若序列的非零值点仅分布在某一点的右边,即
(5-10)
则此序列称为右边序列,其z变换为
(5-11)
设序列为有界序列,假定我们已知这个序列的z变换在处收敛,既有
(5-12)
现分两种情况讨论这个和式的收敛域,
第一种情况,,当时,
图5-2右边序列z变换的收敛域
即至少在的区域内是收敛的,这是一个圆外区域,其包含了处,如图5-2所示。
第二种情况,,这时序列的z变换可以写成
上式中第一个和式可以看成有限长序列的第二种情况,它的收敛域为除以外的整个z平面,而上式中的第二个和式就是右边序列的第一种情况,它的收敛域为一个圆外区域。综合这两个和式的收敛域就是序列z变换的收敛域,即,这是一个不包括处的圆外区域。
例:已知序列 ,求X(z)。
当 或
此例收敛域是以X(z)的极点1/3为半径的圆外。推论:在X(z)的封闭表示式中,若有多个极点,则右边序列的收敛区是以绝对值最大的极点为收敛半径的圆外。
3. 左边序列(left-sided frequence)
若序列的非零值点仅分布在某一点的左边,即
(5-13)
则称此序列为左边序列。其z变换为
(5-14)
用类似于右边序列的讨论,假定在处收敛,既有
第一种情况,,对于所有有
图5-3左边序列z变换的收敛域
由此可见,的收敛域是个圆内区域,且包含了z=0处,即,如图5-3所示。
第二种情况,,这时序列的z变换可以写成
上式中第一个和式就是左边序列的第一种情况,它的收敛域是个圆内区域,且包含了z=0处;第二个和式是有限长序列的第三种情况,它的收敛域为除z=0点外的整个z平面。综合上述分析结果可知这时序列z变换的收敛域为一个圆内区域,但不包含z=0点,即。
例:已知序列 ,求X(z)。
解
此例收敛域是以X(z)的极点b为半径的圆内。推论:在X(z)的封闭表示式中,若有多个极点,则左边序列的收敛区是以绝对值最小的极点为收敛半径的圆内。
4. 双边序列(two-sided frequence)
若序列的非零值点分布在整个整数集上,则此序列称为双边序列。双边序列的z变换刻意写为:
图5-4双边序列的z变幻的收敛域
上式中第一和式就是左边序列的第一种情况,它的收敛域为包括z平面原点的一个圆内区域,设为;第二个和式为右边序列的第一种情况,它的收敛域为包括无穷远处(z=∞)的一个圆外区域,设为。当时,设两个和式有公共的收敛区域,它是个圆环,这就是双边序列z变换的收敛域,如图5-4所示。如果,则左边序列与右边序列的z变换的收敛域没有公共区域,因此这时双边序列的z变换不存在。
例: 已知双边序列 ,c为实数,求X(z)。
解
n<0时,
n<0时,
5.典型序列的z变换
1.δ(n)
Z[δ(n)]=
2.u(n)
Z[u(n)]= |z-1|<1
= |z|>1
3.斜变序列nu(n)
Z[nu(n)]= |z-1|<1
可利用u(n)的z变换 = |z|>1
等式两边分别对z-1求导,得
两边各乘以z-1 |z|>1
4.指数序列
(1)anu(n) Z[anu(n)]= |z|>|a|
(2) -anu(-n-1) Z[-anu(-n-1)]= |z|<|a|
现举例说明序列的z变换,
已知,求序列的z变换。
解:
单位抽样序列的z变换是个常数,显然无论z取何值,z变换都收敛,因此它收敛域为整个z平面。
再来看指数序列的z变换,已知,它的z变换为:
若,则上式收敛,即当时,有
这个z变换的收敛域是个圆外区域,因为序列是个右边序列。
已知序列,求序列的z变换。
解: 这是个左边序列,有
要使上式收敛,必有,即
这个z变换的收敛域为圆内区域,与我们前面推导的结论是一致的。
常见序列的z变换见教课书177页表6.1所示。
5-4 z变换的性质
现在讨论z变换的性质,设
1. 线性特性
(5-15)
式中a,b均为常数,通常两个序列z变换叠加后其收敛域为这两个z变换收敛域的重叠部分,即
但当线性组合中出现零极点相消时,z变换的收敛域可能扩大,例如
而
这时收敛域为整个z平面,显然收敛域扩大了。
2. 时移特性
(5-16)
式中m是可正可负的整数,序列位移后的z变换除了在z=0或z=∞处收敛情况与原序列z变换收敛情况可能不同外,其它区域的收敛情况完全相同。注意这里时移特性指的是双边z变换的特性。
注[6]:除了位移特性外,其它各项特性,单边z变换都可以看成是因果序列的双边z变换,因而所有公式都同样适用。
对于单边z变换,即
则序列位移的z变换要考虑序列的初值,即当时有
左移特性
右移特性
3. z域微分特性
(5-17)
这是因为
上式两边对z求导有
所以
4. z域尺度变换特性
(5-18)
这时因为
尺度变换特性表明,序列乘上指数序列后,相当于z平面的尺度变化为
5. 时域卷积特性
(5-19)
这是个重要特性,它将时域中的卷积关系转换为z域中的乘积关系。
证明: 设,
则
证毕
注:如果两个序列在时域中是乘积的关系,则这两个序列在z域中是复数卷积的关系(不是一般的卷积),即
6. 初值定理
若为因果序列,即,则
这是因为
上式两边在时取极限,则右边除了项外,其它各项均为零,故有
7. 终值定理
若为因果序列,且处在处可以有一阶极点外,全部其它极点都在单位圆内,则
(5-20)
这是因为
考虑到因果序列的特性,上式可改写为
注意上式的收敛域,由于在单位圆上只有处可能有一阶极点,现在函数将抵消这个可能的极点,因此的收敛域必将包含单位圆,这就允许上式两边取极限,即
5-5 逆z变换
在复变函数理论中逆z变换的定义公式为:
(5-21)
这是个围线积分,C是在收敛域内,反时针包围z平面坐标原点的闭合曲线。逆z变换通常也可记为
因为(5-21)式的积分计算很麻烦,所以一般采用以下介绍的三种方法求解逆z变换。
1留数定理法
若在闭合曲线C内的所有极点的集合为,则的逆变换可由下式求出
(5-22)
式中为函数在极点上的留数,表示对整个极点集合上的留数求和。即等于函数在闭合曲线C内所有极点上留数的总和。
如果为单阶极点,则极点上的留数为
(5-23)
如果为N阶极点,则极点上的留数为
(5-25.)
例如:已知,求逆z变换。
解:
图5-5的收敛域及其零极点
注意的收敛域为一个圆外区域,且收敛域包含处,如图5-5所示。所以序列是右边序列的第一种情况,它的非零值点一定在坐标原点(n=0)的右边,既有
因此只需求出时的。
当时,函数在处有一个单阶极点,在处也有一个阶极点。在收敛域内画一条包含z平面坐标原点的围线C如图5-5所示,则在闭合曲线C内的极点只有一个单阶极点,因此
注意这是时的解。综合以上分析,有
图5-6的收敛域及其零极点
在求逆z变换时,不仅要考虑表达式,而且也要注意它的收敛域,相同的不同的收敛域完全可以导出不同的序列。
例如
此例的z变换表达式与上例完全相同,但二者的收敛域不同。此例中的收敛域是个圆内区域,如图5-6所示,这个圆内区域包含了处,所以这个z变换对应着左边序列的第一种情况,即它的非零值点分布在坐标原点(n=0)的左边,则有
现在只需讨论时的,因为
当n=0时,有
这个函数只有在z=a处有个一阶极点,它在闭合曲线C的外面。也就是说,在闭合曲线C内这个函数没有极点,因此它的留数为零,即
当n<0时,函数在闭合曲线内z=0处有个|n|阶极点,则有
综合以上分析有
以上两例充分说明了收敛域对z变换的重要性。
2幂级数法
如果是个右边序列或左边序列的第一种情况,则可将z变换写成幂级数形式:
(5-25)
或
(5-26)
(5-25)式称为降幂级数,对应于一个右边序列,即有;(5-26)式称为升幂级数,对应于一个左边序列,既有。
当为有理函数时,可以用长除法将其展开成幂级数形式,在某些情况下,也可以运用一些数学技巧将z变换展开成幂级数,例如
这个z变换的收敛域为一个圆外区域,它一定对应着一个右边序列,因此应将展开成降幂级数形式。查数学手册可知
令,则有
既有
所以
或者写成
现在介绍用长除法将z变换展开成幂级数,下面的讨论通过例子说明。
已知
图5-7序列的降幂除法
的收敛域为一个圆外区域,且包括无穷远处(),所以可以判定是个右边序列的第一种情况。
为求出右边序列必须用长除法将展开成降幂形式的级数,即将的分子、分母多项式都写成降幂形式,然后再做长除法,如图5-7所示,
即有
因此
这个结果与用留数定理法求出的结果一样。
现在再来求,的逆z变换。
因为的收敛域是个圆内区域,且包含了处,所以可以判定是左边序列的第一种情况。为了将左边序列的z变换展开成幂级数,必须将的分子、分母多项式按升幂形式排列,即
再做长除法如图5-8所示,则有
图5-8序列的升幂除法
所以
或写成
这个结果与前面留数定理求解结果一致。
以上两例对应的是单边序列的逆z变换,即或为右边序列的逆z变换,或为左边序列的逆z变换。如果对应的是双边序列,即的收敛域为圆环区域的情况下,又如何运用幂级数法求解逆z变换呢?我们通过讨论一个双边序列的例子来进行说明。
已知z变换
Re[z]
jIm[z]
Z=-2
Z=1
根据的收敛域可以判定这个z变换对应一个双边序列,将上式写成
图5-9双边序列的收敛域
由上式可以看出有两个极点,一个在圆环内,一个在圆环外,如图5-9所示。实际上,双边序列的z变换可以分成两部分
式中包含了在圆环内的所有极点,而不包含圆外的任何极点;与此相反,包含了圆环外的所有极点,而不包含圆环内的任何极点。这样分解后,的收敛域为一个圆外区域,它对应一个右边序列;的收敛域为一个圆内区域,它对应一个左边序列。对本例来说,可以这样分解
既有
注意和的分解不是唯一的,只要求包含在收敛域环内的所有极点,而包含在收敛域环外的所有极点。
因为仅包含收敛域环内的极点,所以它的收敛域必定是个圆外区域,因此它对应一个右边序列,我们可以用降幂除法求解逆z变换。即
因为仅包含收敛域环外的极点,所以它的收敛域必定是个圆内区域,因此它对应一个左边序列,我们可以用升幂除法求解逆z变换。即
根据z变换的线性特性可知的逆z变换为和的逆z变换之和,即
由上述例子可以看出,用幂级数法求解逆z变换时非常有效的,但也比较麻烦。
例3利用z变换的线性特性,将一个双边序列分解成两个单边序列之和,然后分别利用降幂长除法和升幂长除法求出逆z变换。由此可以得到启发,将一个复杂的有理函数展开成简单的有理函数之和,而这些有理函数的逆z变换又可以通过简单的方法(如查表)得到。根据这个思想我们可以导出逆z变换的另一种方法——部分分式展开法。
3部分分式法
z变换的部分分式展开法与拉普拉斯的部分分式展开法基本相似。通常情况下,部分分式展开法用于因果序列的逆z变换比较有效。其基本思想是:先对分解成部分分式,然后将每个部分分式再乘以z,这样就可以得到以z-1为变量的部分分式展开式,再通过一些简单的方法(通常为查数学手册)求出这些简单部分分式的逆z变换,将部分分式的逆z变换求和就是原z变换的逆z变换。
一个因果序列的z变换通常可以写成有理分式,即
如果的N个极点都是单阶极点,那么可以展开成部分分式:
式中是常数,可以由下式求出:
对应的逆z变换为
一般情况下,当具有高阶极点时部分分式法并不是一个有效的方法,所以本课程不再讲述。
例 已知 ,求x(n)。
解:1<|z|,是右边(因果)序列。
5-6 LTI系统分析
z变换的重要应用与拉普拉斯变换的应用是一样的,用于对LTI系统的分析与表征。拉普拉斯变换的卷积定理指出:
若
则
一个线性非时变系统,对任意输入的零状态响应可以用卷积求得,即
式中为系统的单位冲激响应。运用拉普拉斯卷积定理可以得到如下关系
式中,和分别为系统的输入,输出和单位冲激响应的拉普拉斯变换。在拉普拉斯变换域(通常由叫做s域)内,系统单位冲激响应的拉普拉斯变换称为系统函数(或转移函数)。根据上式关系,我们有
上式的含义是零状态系统的输入、输出信号的拉普拉斯变换之比就等于系统函数,而系统函数的拉普拉斯反变换就是系统的单位冲激响应,有时候我们可以利用上述关系来求系统的单位冲激响应。
例如已知一个LTI系统的微分方程为:
求系统的单位冲激响应和单位阶跃响应。
按照前面介绍过的方法,可以通过解微分方程方法来求解,这里要采用拉普拉斯变换求解。先求的系统单位冲激响应,对微分方程两边取拉普拉斯变换,得
即
所以系统函数为
可以用部分分式法将上式展开成部分分式,即
式中A和B为常数可以用下式求出:
所以
求系统函数的拉普拉斯反变换(查表P130)就得到系统的单位冲激响应,即
现在求系统的单位阶跃响应,令输入,它的拉普拉斯变换为,所以有
式中系数可以用下式求出
即
查表可得拉普拉斯反变换,即
与连续时间系统类似,LTI系统对任意输入,输出响应可以由卷积求出,即
式中为系统的单位抽样响应。根据z变换的卷积定理,有
式中,和分别是输入,输出和单位抽样响应的z变换。注意根据以上关系,我们仍可得到:
在z域中也被称为系统函数。上式的含义是零状态系统的输入,输出序列的z变换之比就等于系统函数。
对于离散时间系统我们运用z变换作为分析工具。例如已知一个LTI因果离散时间系统的差分方程为
求系统的单位抽样响应和单位阶跃响应。
先求系统的单位抽样响应。对差分方程两边求z变换得:
即
系统函数为(因为系统是因果的,所以的收敛域一定为一个圆外区域)
求逆z变换得系统的单位抽样响应
现在求系统的单位阶跃响应,零输入,它的z变换为,于是
因为和的收敛域都是圆外区域,所以的收敛域至少也是个圆外区域,即对应着一个右边序列。用部分分式法将展开
式中常数A和B可以有下式求出
即
于是系统的单位阶跃响应为
通过以上两个例子我们看到了拉普拉斯变换和z变换分别在求解连续时间系统的微分方程和离散时间系统差分方程的作用。其实拉普拉斯变换和z变换的作用远不止此,在前面我们已经介绍过系统的单位冲激响应反映了连续时间系统的时域特性,而单位抽样响应反映了离散时间系统的时域特性,它们都取决于系统的结构和参数。显然连续时间系统和离散时间系统的系统函数也取决于系统的结构和参数,它是系统的复频域描述函数。特定系统具有不同的特性(如高通、低通、带通等),但所有物理可实现系统都要满足稳定性和因果性。即系统的稳定性和因果性是物理系统的最基本特性。
在第一章中我们已经介绍了系统稳定性和因果性的概念,现在讨论系统满足因果性和稳定性的条件是什么。换句话说,系统的单位冲激响应或单位抽样响应和系统函数是如何反映系统的稳定性和因果性的。
因果性
概念描述:如果一个系统在任何时刻的输出之取决于现在的输入以及过去的输入,该系统就成为因果系统。
简单叙述:系统的输出响应不超前与输入激励。
数学描述:对于连续时间系统,如果,则有
对于离散时间系统,如果,则有
对于一个LTI连续时间系统,其因果性的充分必要条件是
式中是系统的单位冲激响应。这是一个右边函数,它的拉普拉斯变换为连续系统函数,一般说来,在这种情况下的收敛域为s域的某个右半平面。即一个因果系统的系统函数的收敛域为s域的某个右半平面,但相反的结论未必成立。进一步的约束条件是:对于一个具有有理系统函数的系统来说,系统的因果性就等效于收敛域位于最右边极点的右边的右半平面。
对于一个LTI离散时间系统,其因果性的充分必要条件是:
这是个右边序列,它的z变换为离散系统函数,对一个离散时间LTI系统当且仅当它的系统函数收敛域是个圆外区域,且包含无穷远处点,则该系统是因果的。
在这里仅对离散系统因果性充分必要条件进行证明:
充分性 若
则系统对任何输入的响应为
因为n<0时,h(n)=0,则当k>n时,h(n-k)=0,于是卷积和可以写成
上式说明,在计算y(n)时x(n)只在区间上取值,这就是说y(n)只与n以及以前的输入值有关,这就是系统的因果性。
必要性(用反证法)
假定一个因果系统的单位抽样响应h(n)在n<0时,至少有一个非零值点,即有
则系统对任意输入的响应为
在上式的第二项和式中k>n,即n-k<0,根据假定当n-k=i时,,这就是说计算y(n)时,至少与这项输入有关,显然这项输入对n这个时刻来说是未来的输入,这个结果与因果系统的定义矛盾,所以假定不能成立。即对因果系统来说,在n<0时,不可能存在非零值点。
稳定性
一个LTI连续系统稳定的充分必要条件就是它的单位冲激响应绝对可积,即
这等价于系统函数的收敛域要包含s域中的虚轴(),即一个LTI连续系统,当且仅当系统函数的收敛域包含s域中的虚轴()时,该系统是稳定的。
一个LTI离散系统稳定的充分必要条件是它的单位抽样响应绝对可和,即
这等价价于系统函数的收敛域要包含z平面中的单位圆(),即一个LTI离散系统当且仅当它的系统函数的收敛域包含z平面中的单位圆()时,该系统是稳定的。
5-7系统函数
可以用单位脉冲响应h(n)来表示LTI离散系统的输入输出关系
对应的 变换为
定义LTI离散系统输出 变换与输入 变换之比为系统函数
(5.7-1)
当x(n)=δ(n),H(z)=Y(z)。所以系统函数是系统单位脉冲响应h(n)的 变换。
(5.7-2)
N阶LTI离散系统的差分方程通常为
(5.7-3)
系统为零状态时,对两边取 变换,可得
解出
得到系统函数
(5.7-4)
式(6.7-4)是z-1的有理分式,其系数正是差分方程的系数,系统函数还可以分解为
(5.7-5)
式中,{ck}是H(z) 的零点;{dk}是H(z)的极点。由式(5.7-5)可见,除了系数A外,H(z)可由其零、极点确定。将零点{ck}与极点{dk}标在z平面上,可得到离散系统的零、极点图。
当离散系统的系统函数有原点以外的任意极点时,即式(5.7-5)中有dk≠0时,对应的单位脉冲响应h(n) 的时宽为无限,这样的系统称为无限冲激响应系统,简称IIR系统;当离散系统得系统函数只有原点处极点时,即式(5.7-5)中所有dk=0时,对应的单位脉冲响应h(n) 的时宽为有限,这样的系统称为有限冲激响应系统,简称FIR系统。FIR系统函数一般表示为
(5.7-6)
由于FIR系统可以具有线性相位,并且不存在系统稳定问题,因此得到越来越广泛的应用。
与连续系统相似,离散系统的特性可其零、极点分布密切相关,但将系统函数由有理分式形式分解为零、极点形式时,并不容易,而用MATLBA可以很方便地确定零、极点并作零、极点图。
例 已知某系统函数为
求其零、极点并绘出零、极点图。
解: MATLBA程序及结果如下:
b=[0.2 0.1 0.3 0.1 0.2];% 分子多项式系数
a=[1 -1.1 1.5 -0.7 0.3];%分母多项式系数
r1=roots(a)%求极点
r2=roots(b)%求零点
zplane(b,a)
答案
r1= 0.2367+0.8915i
0.2367-0.8915i
0.3133+0.5045i
0.3133-0.5045i
r2= -0.5000+0.8660i
-0.5000-0.8660i
0.2500+0.9682i
0.2500-0.9682i
离散系统的零、极点图如图5.7-1所示。
1
0.8
0.6
0.4
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
图5.7-1 例5.7-1系统的零、极点图
习 题
5.1求出下列序列的z变换
(a) (b)
5.2考虑序列
求出X(z),并画出X(z)的零极点。
5.3求出下列序列的z变换X(z),画出零极点图和收敛域:
(a)
(b)
(c)
5.4设
a>0
(a)画出a<1和a>1时的x[n]。
(b)求出X(s),画出零极点图以及a>1和a<1时的收敛域。
5.5求出下列序列的z变换及其相关的收敛域:
(a)
(b)
(c)
(d)
(e) ,a 、b为常数
5.6求出下列序列的z变换:
(a) (b)
5.7求出下面X(z)的反拉普拉斯变换:
,
5.7利用幂级数展开技术,求出下列X(z)的z反变换:
(a)
(b)
5. 8利用幂级数展开技术,求出下列X(z)的z反变换:
(a) (b)
5.9求出下面X(z)的z反变换:
5.10 一个离散时间LTI系统,如果输入x[n]为u[n],输出y[n]为 。
(a) 求出系统的冲激响应h[n]。
(b) 求出输入x[n]为 时的输出y[n]。
5.11求出下列序列x[n]的z变换:
(a)
(b)
5.12已知下列差分方程以及相关的输入和初始条件,确定输出y[n]。
(a)
(b)
151
展开阅读全文