《数字与信号处理第1章.ppt》由会员分享,可在线阅读,更多相关《数字与信号处理第1章.ppt(262页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第1章 离散时间信号与系统 第1章 离散时间信号与系统 1.1离散时间信号离散时间信号序列序列1.2连续时间信号的采样连续时间信号的采样1.3离散时间系统时域分析离散时间系统时域分析1.4Z变换变换1.5拉氏变换、拉氏变换、傅氏变换与傅氏变换与Z变换变换1.6离散时间系统的频域分析(离散时间系统的频域分析(域和域和域)域)第1章 离散时间信号与系统 1.1离散时间信号离散时间信号序列序列离散时间信号只在离散时间上给出函数值,是时间上不连续的一个序列。它既可以是实数也可以是复数。一个离散时间信号是一个整数值变量n的函数,表示为x(n)或x(n)。尽管独立变量n不一定表示“时间”(例如,n可以表示
2、温度或距离),但x(n)一般被认为是时间的函数。因为离散时间信号x(n)对于非整数值n是没有定义的,所以一个实值离散时间信号序列可以用图形来描述,如图1-1所示。横轴虽为连续直线,但只在n为整数时才有意义。纵轴线段的长短代表各序列值的大小。第1章 离散时间信号与系统 图1-1离散时间信号x(n)的图形表示第1章 离散时间信号与系统 离散时间信号常常可以对模拟信号(如语音)进行等间隔采样而得到。例如,对于一个连续时间信号xa(t),以每秒fs=1/T个采样的速率采样而产生采样信号,它与xa(t)的关系如下:然而,并不是所有的离散时间信号都是这样获得的。一些信号可以认为是自然产生的离散时间序列,如
3、每日股票市场价格、人口统计数和仓库存量等。第1章 离散时间信号与系统 1.1.1序列的运算序列的运算1序列的移位序列的移位如图1-1所示的序列x(n),其移位序列w(n)为当m为正时,则x(n-m)是指序列x(n)逐项依次延时(右移)m位而给出的一个新序列;当m为负时,x(n-m)是指依次超前(左移)m位。图1-2显示了x(n)序列的延时序列w(n)=x(n-2),即m=2时的情况。第1章 离散时间信号与系统 图1-2图1-1序列x(n)的延时第1章 离散时间信号与系统 2序列的翻褶序列的翻褶如果序列为x(n),则x(-n)是以n=0的纵轴为对称轴将序列x(n)加以翻褶。x(n)及x(-n)如
4、图1-3(a)、(b)所示。图1-3序列的翻褶(a)x(n)序列;(b)x(-n)序列第1章 离散时间信号与系统 3序列的和序列的和两序列的和是指同序号n的序列值逐项对应相加而构成的一个新序列。和序列z(n)可表示为第1章 离散时间信号与系统 4序列的乘积序列的乘积两序列相乘是指同序号n的序列值逐项对应相乘。乘积序列f(n)可表示为第1章 离散时间信号与系统 5序列的标乘序列的标乘序列x(n)的标乘是指x(n)的每个序列值乘以常数c。标乘序列f(n)可表示为第1章 离散时间信号与系统 6累加累加设某序列为x(n),则x(n)的累加序列y(n)定义为它表示y(n)在某一个n0上的值y(n0)等于
5、在这一个n0上的x(n0)值与n0以前所有n上的x(n)之和。第1章 离散时间信号与系统 7差分运算差分运算前向差分x(n)=x(n+1)-x(n)后向差分x(n)=x(n)-x(n-1)由此得出x(n)=x(n-1)第1章 离散时间信号与系统 1.1.2几种常用序列几种常用序列1单位脉冲序列单位脉冲序列(n)这个序列只在n=0处有一个单位值1,其余点上皆为0,因此也称为“单位采样序列”。单位采样序列如图1-4所示。(1-1)第1章 离散时间信号与系统 图1-4(n)序列第1章 离散时间信号与系统 这是最常用、最重要的一种序列,它在离散时间系统中的作用,很类似于连续时间系统中的单位冲激函数(t
6、)。但是,在连续时间系统中,(t)是t=0点脉宽趋于零,幅值趋于无限大,面积为1的信号,是极限概念的信号,并非任何现实的信号。而离散时间系统中的(n),却完全是一个现实的序列,它的脉冲幅度是1,是一个有限值。第1章 离散时间信号与系统 2单位阶跃序列单位阶跃序列u(n)如图1-5所示。它很类似于连续时间信号与系统中的单位阶跃函数u(t)。(1-2)第1章 离散时间信号与系统 图1-5u(n)序列第1章 离散时间信号与系统(n)和u(n)间的关系为这就是u(n)的后向差分。而令n-m=k,代入此式可得这里就用到了累加的概念。(1-3)(1-4)(1-5)第1章 离散时间信号与系统 3矩形序列矩形
7、序列RN(n)(1-6)矩形序列RN(n)如图1-6所示。第1章 离散时间信号与系统 图1-6RN(n)序列第1章 离散时间信号与系统 RN(n)和(n)、u(n)的关系为:(1-7)(1-8)第1章 离散时间信号与系统 4实指数序列实指数序列式中,a为实数。当|a|1时,序列是发散的。a为负数时,序列是摆动的,如图1-7所示。第1章 离散时间信号与系统 图1-7指数序列(a)|a|1;(c)a=-|a|第1章 离散时间信号与系统 5正弦型序列正弦型序列x(n)=Asin(n0+)(1-10)式中:A为幅度;为起始相位;0为数字域的频率,它反映了序列变化的速率。0=0.1时,x(n)序列如图1
8、-8所示,该序列值每20个重复一次循环。第1章 离散时间信号与系统 图1-8正弦序列(0=0.1)第1章 离散时间信号与系统 6复指数序列复指数序列序列值为复数的序列称为复指数序列。复指数序列的每个值具有实部和虚部两部分。复指数序列是最常用的一种复序列:(1-11a)或或(1-11b)式中,0是复正弦的数字域频率。第1章 离散时间信号与系统 对第一种表示,序列的实部、虚部分别为如果用极坐标表示,则因此有:第1章 离散时间信号与系统 1.1.3序列的周期性序列的周期性如果对所有n存在一个最小的正整数N,满足(1-12)则称序列x(n)是周期性序列,周期为N。现在讨论上述正弦序列的周期性。由于则第
9、1章 离散时间信号与系统 若N0=2k,当k为正整数时,则这时的正弦序列就是周期性序列,其周期满足N=2k/0(N,k必须为整数)。可分几种情况讨论如下。(1)当2/0为正整数时,周期为2/0,见图1-8。(2)当2/0不是整数,而是一个有理数时(有理数可表示成分数),则式中,k,N为互素的整数,则为最小正整数,序列的周期为N。第1章 离散时间信号与系统(3)当2/0是无理数时,则任何k皆不能使N取正整数。这时,正弦序列不是周期性的。这和连续信号是不一样的。同样,指数为纯虚数的复指数序列的周期性与正弦序列的情况相同。下面,我们来进一步讨论,如果一个正弦型序列是由一个连续信号采样而得到的,那么,
10、采样时间间隔T和连续正弦信号的周期之间应该是什么关系才能使所得到的采样序列仍然是周期序列呢?第1章 离散时间信号与系统 设连续正弦信号xa(t)为这 一 信 号 的 频 率 为 f0,角 频 率 0=2f0,信 号 的 周 期 为T0=1/f0=2/0。如果对连续周期信号xa(t)进行采样,其采样时间间隔为T,采样后信号以x(n)表示,则有如果令0为数字域频率,满足第1章 离散时间信号与系统 式中,fs是采样频率。可以看出,0是一个相对频率,它是连续正弦信号的频率f0对采样频率fs的相对频率乘以2,或说是连续正弦信号的角频率0对采样频率fs的相对频率。用0代替0T,可得这就是我们上面讨论的正弦
11、型序列。第1章 离散时间信号与系统 下面我们来看2/0与T及T0的关系,从而讨论上面所述正弦型序列的周期性的条件意味着什么?这表明,若要2/0为整数,就表示连续正弦信号的周期T0应为采样时间间隔T的整数倍;若要2/0为有理数,就表示T0与T是互为互素的整数,且有(1-13)式中,k和N皆为正整数,从而有即N个采样间隔应等于k个连续正弦信号的周期。第1章 离散时间信号与系统 1.1.4用单位采样序列来表示任意序列用单位采样序列来表示任意序列用单位采样序列来表示任意序列对分析线性时不变系统(下面即将讨论)是很有用的。设x(m)是一个序列值的集合,其中的任意一个值x(n)可以表示成单位采样序列的移位
12、加权和,即(1-14)由于第1章 离散时间信号与系统 则因此,式(1-14)成立,这种表达式提供了一种信号分析工具。第1章 离散时间信号与系统 1.1.5序列的能量序列的能量序列x(n)的能量E定义为序列各采样样本的平方和,即(1-15)第1章 离散时间信号与系统 1.2连续时间信号的采样连续时间信号的采样在某些合理条件限制下,一个连续时间信号能用其采样序列来完全给予表示,连续时间信号的处理往往是通过对其采样得到的离散时间序列的处理来完成的。本节将详细讨论采样过程,包括信号采样后,信号的频谱将发生怎样的变换,信号内容会不会丢失,以及由离散信号恢复成连续信号应该具备哪些条件等。采样的这些性质对离
13、散信号和系统的分析都是十分重要的。要了解这些性质,让我们首先从采样过程的分析开始。第1章 离散时间信号与系统 采样器可以看成是一个电子开关,它的工作原理可由图1-9(a)来说明。设开关每隔T秒短暂地闭合一次,将连续信号接通,实现一次采样。如果开关每次闭合的时间为秒,那么采样器的输出将是一串周期为T,宽度为的脉冲。而脉冲的幅度却是重复着在这段时间内信号的幅度。如果以xa(t)代表输入的连续信号,如图1-9(b)所示,以xp(t)表示采样输出信号,它的结构如图1-9(d)所示。显然,这个过程可以把它看作是一个脉冲调幅过程。被调制的脉冲载波是一串周期为T、宽度为的矩形脉冲信号,如图1-9(c)所示,
14、并以p(t)表示,而调制信号就是输入的连续信号。因而有第1章 离散时间信号与系统 一般开关闭合时间都是很短的,而且越小,采样输出脉冲的幅度就越准确地反映输入信号在离散时间点上的瞬时值。当2h);(d)已采样信号频谱(s2h)第1章 离散时间信号与系统 如果信号的最高频谱h超过s/2,则各周期延拓分量产生频谱的交叠,称为混叠现象,如图1-10(d)所示。由于Xa(j)一般是复数,所以混叠也是复数相加。为了简明起见,在图1-10中我们将Xa(j)作为标量来处理。我们将采样频率之半(s/2)称为折叠频率,即它如同一面镜子,当信号频谱超过它时,就会被折叠回来,造成频谱的混叠。(1-30)第1章 离散时
15、间信号与系统 图1-11说明了在简单余弦信号情况下频谱混叠的情况。在图1-11(a)中,给出该余弦信号(1-31)的傅里叶变换Xa(j)。(1-32)图(b)是在0s/2时,的傅里叶变换。(d)和(e)则分别对应于0/T时低通滤波器输出的傅里叶变换,在没有混叠时((b)和(d)),恢复出的输出ya(t)为第1章 离散时间信号与系统 在有混叠时,则是(1-33)这就是说,作为采样和恢复的结果,高频信号cos0t已经被当作和低频信号cos(s-0)t是一样的东西被冒名顶替了。这个讨论就是奈奎斯特采样定理的基础。第1章 离散时间信号与系统 图1-11一个余弦信号采样中的混叠效果第1章 离散时间信号与
16、系统 图1-11一个余弦信号采样中的混叠效果第1章 离散时间信号与系统 由此得出结论:要想采样后能够不失真地还原出原信号,则采样频率必须大于两倍信号谱的最高频率(s2h),这就是奈奎斯特采样定理。即fs2fh频率h一般称为奈奎斯特频率,而频率2h称为奈奎斯特率。采样频率必须大于奈奎斯特率。在实际工作中,为了避免频谱混淆现象发生,采样频率总是选得比奈奎斯特频率更大些,例如选到(34)h。同时为了避免高于折叠频率的杂散频谱进入采样器造成频谱混淆,一般在采样器前加入一个保护性的前置低通滤波器,称为防混叠滤波器,其截止频率为s/2,以便滤除掉高于s/2的频率分量。第1章 离散时间信号与系统 同样方法,
17、可以证明(亦可代j=s到式(1-27),理想采样后,使信号的拉普拉斯变换在S平面上沿虚轴周期延拓。也就是说,在S平面虚轴上是周期函数。即有(1-34)式中:即分别是的双边拉普拉斯变换。第1章 离散时间信号与系统 1.2.3采样的恢复采样的恢复如果理想采样满足奈奎斯特定理,即模拟信号谱的最高频率小于折叠频率则采样后不会产生频谱混叠,由式1-27)知故将通过一个理想低通滤波器,这个理想低通滤波器应该只让基带频谱通过,因而其带宽应该等于折叠频率,它的特性如图1-12所示。第1章 离散时间信号与系统 图1-12采样的恢复第1章 离散时间信号与系统 采样信号通过这个滤波器后,就可滤出原模拟信号的频谱因此
18、,在输出端可以得到原模拟信号理想低通滤波器虽不可实现,但是在一定精度范围内,可用一个可实现的滤波器来逼近它。第1章 离散时间信号与系统 1.2.4由采样信号序列重构带限信号由采样信号序列重构带限信号理想低通滤波器的冲激响应为由与h(t)的卷积积分,即得理想低通滤波器的输出为第1章 离散时间信号与系统 这里h(t-nT)称为内插函数:(1-35)第1章 离散时间信号与系统 它的波形如图1-13所示,其特点为:在采样点nT上,函数值为1;其余采样点上,函数值都为零。图1-13内插函数第1章 离散时间信号与系统 由于ya(t)=xa(t),因此以上卷积结果也可以表示为(1-36)式(1-36)称为采
19、样内插公式,即信号的采样值xa(nT)经此公式而得到连续信号xa(t)。也就是说,xa(t)等于各xa(nT)乘上对应的内插函数的总和。在每一采样点上,只有该点所对应的内插函数不为零,这使得各采样点上信号值不变,而采样点之间的信号则由加权内插函数波形的延伸叠加而成,如图1-14所示。这个公式说明了,只要采样频率高于两倍信号最高频率,则整个连续信号就可以完全用它的采样值来代表,而不会丢掉任何信息。这就是奈奎斯特采样定理的意义。由上面讨论可看出采样内插公式只限于使用到限带(频带有限)信号上。第1章 离散时间信号与系统 图1-14采样内插恢复第1章 离散时间信号与系统 1.3离散时间系统时域分析离散
20、时间系统时域分析一个离散时间系统是将输入序列变换成输出序列的一种运算。若以T来表示这种运算,则一个离散时间系统可由图1-15来表示,即(1-37)离散时间系统中最重要、最常用的是“线性时不变系统”。第1章 离散时间信号与系统 图1-15离散时间系统第1章 离散时间信号与系统 1.3.1线性系统线性系统满足叠加原理的系统称为线性系统,即若某一输入是由N个信号的加权和组成,则输出就是系统对这几个信号中每一个的响应的同样加权和组成。如果系统在x(n)和x2(n)单独输入时的输出分别为y1(n)和y2(n)即:那么当且仅当式(1-38a)和式(1-38b)成立时,该系统是线性的(1-38a)第1章 离
21、散时间信号与系统 和(1-38b)式中,a为任意常数。上述第一个性质称为可加性,第二个称为齐次性或比例性。这两个性质合在一起就成为叠加原理,写成(1-39)式(1-39)对任意常数a1和a2都成立。该式还可推广到多个输入的叠加,即(1-40)式中,yk(n)就是系统对输入xk(n)的响应。在证明一个系统是线性系统时,必须证明此系统同时满足可加性和比例性,而且信号以及任何比例常数都可以是复数。第1章 离散时间信号与系统 例例1-1以下系统是否为线性系统:y(n)=2x(n)+3很容易证明这个系统不是线性的,因为此系统不满足叠加原理。证证很明显,在一般情况下所以此系统不满足叠加性,故不是线性系统。
22、第1章 离散时间信号与系统 同样可以证明,第1章 离散时间信号与系统 1.3.2时不变系统时不变系统系统的运算关系T在整个运算过程中不随时间(也即不随序列的延迟)而变化,这种系统称为时不变系统(或称移不变系统)。这个性质可用以下关系表达:若输入x(n)的输出为y(n),则将输入序列移动任意位后,其输出序列除了跟着移位外,数值应该保持不变,即若Tx(n)=y(n)则 Tx(n-m)=y(n-m)(m为任意整数)(1-41)满足以上关系的系统就称为时不变系统。第1章 离散时间信号与系统 例例1-2证明不是时不变系统。证证由于二者不相等,故不是时不变系统。同时具有线性和时不变性的离散时间系统称为线性
23、时不变(LTI)离散时间系统,简称LTI系统。除非特殊说明,本书都是研究LTI系统。第1章 离散时间信号与系统 1.3.3单位脉冲响应与系统的输入输出关系单位脉冲响应与系统的输入输出关系线性时不变系统可用它的单位脉冲响应来表征。单位脉冲响应是指输入为单位脉冲序列时系统的输出。一般用h(n)表示单位脉冲响应,即h(n)=T(n)有了h(n)我们就可以得到此线性时不变系统对任意输入的输出。下面讨论这个问题:设系统输入序列为x(n),输出序列为y(n)。从式(1-14)已经知道,任一序列x(n)可以写成(n)的移位加权和,即第1章 离散时间信号与系统 则系统的输出为由于系统是线性的,可利用叠加原理式
24、(1-40),则又由于系统的时不变性,式(1-41)对移位的单位脉冲的响应就是单位脉冲响应的移位。因此(1-42)第1章 离散时间信号与系统 如图1-16所示。上式称为序列x(n)与h(n)的离散卷积,为了同以后的圆周卷积相区别,离散卷积也称为“线性卷积”或“直接卷积”或简称“卷积”,并以“*”表示之。图1-16线性时不变系统第1章 离散时间信号与系统 图1-17离散卷积第1章 离散时间信号与系统 图1-17离散卷积第1章 离散时间信号与系统(1)翻褶:先在哑变量坐标m上作出x(m)和h(m),将h(m)以m=0的垂直轴为对称轴翻褶成h(-m)。(2)移位:将h(-m)移位n,即得h(n-m)
25、。当n为正整数时,右移n位;当n为负整数时,左移n位。(3)相乘:再将h(n-m)和x(m)的相同m值的对应点值相乘。(4)相加:把以上所有对应点的乘积累加起来,即得y(n)值。依上法,取n=,-2,-1,0,1,2,各值,即可得全部y(n)值。第1章 离散时间信号与系统 由式(1-42)不难看出,卷积与两序列的先后次序无关。证证令n-m=m代入式(1-42),然后再将m换成m,即得(1-43)因此第1章 离散时间信号与系统 1.3.4线性时不变系统的性质线性时不变系统的性质1交换律交换律由于卷积与两卷积序列的次序无关,即卷积服从交换律,故这就是说,如果把单位脉冲响应h(n)改作为输入,而把输
26、入x(n)改作为系统单位脉冲响应,则输出y(n)不变。(1-44)第1章 离散时间信号与系统 2结合律结合律可以证明卷积运算服从结合律,即这就是说,两个线性时不变系统级联后仍构成一个线性时不变系统,其单位脉冲响应为两系统单位脉冲响应的卷积,且线性时不变系统的单位脉冲响应与它们的级联次序无关,如图1-18所示。(1-45)第1章 离散时间信号与系统 3分配律分配律卷积也服从加法分配律:(1-46)也就是说,两个线性时不变系统的并联等效系统的单位脉冲响应等于两系统各自单位脉冲响应之和,如图1-19所示。以上三个性质,交换律前面已经证明了,另外两个性质由卷积的定义可以很容易加以证明。第1章 离散时间
27、信号与系统 图1-18具有相同单位脉冲响应的三个线性时不变系统第1章 离散时间信号与系统 图1-19线性时不变系统的并联组合及其等效系统第1章 离散时间信号与系统 1.3.5因果系统因果系统所谓因果系统,就是系统此时的输出y(n)只取决于此时,以及此时以前的输入,即x(n),x(n-1),x(n-2),。如果系统的输出y(n)还取决于x(n+1),x(n+2),,也即系统的输出还取决于未来的输入,这样在时间上就违背了因果关系,因而是非因果系统,也即不现实的系统。根据上述定义,可以知道,y(n)=nx(n)的系统是一个因果系统,而y(n)=x(n+2)+ax(n)的系统是非因果系统。第1章 离散
28、时间信号与系统 从式(1-43)卷积公式,我们可以看到线性时不变系统是因果系统的充分必要条件是 h(n)=0 n0(1-47)依照此定义,我们将n0,x(n)=0的序列称为因果序列,表示这个因果序列可以作为一个因果系统的单位脉冲响应。第1章 离散时间信号与系统 我们知道,许多重要的网络,如频率特性为理想矩形的理想低通滤波器以及理想微分器等都是非因果的不可实现的系统。但是数字信号处理往往是非实时的,即使是实时处理,也允许有很大延时。这是对于某一个输出y(n)来说,已有大量的“未来”输入x(n+1),x(n+2),,记录在存储器中可以被调用,因而可以很接近于实现这些非因果系统。也就是说,可以用具有
29、很大延时的因果系统去逼近非因果系统。这个概念在以后讲有限长单位脉冲响应滤波器设计时要常用到,这也是数字系统优于模拟系统的特点之一。因而数字系统可以比模拟系统更能获得接近理想的特性。第1章 离散时间信号与系统 1.3.6稳定系统稳定系统稳定系统是指有界输入产生有界输出(BIBO)的系统。如果对于输入序列x(n),存在一个不变的正有限值Bx,对于所有n值满足|x(n)|Bx(1-48)则称该输入序列是有界的。稳定性要求对于每个有界输入存在一个不变的正有限值By,对于所有n值,输出序列y(n)满足|y(n)|By(1-49)第1章 离散时间信号与系统 一个线性时不变系统是稳定系统的充分必要条件是单位
30、脉冲响应绝对可和,即(1-50)证证充分条件:若如果输入信号x(n)有界,即对于所有n皆有|x(n)|Bx,则即输出信号y(n)有界,故原条件是充分条件。第1章 离散时间信号与系统 必要条件:利用反证法。已知系统稳定,假设我们可以找到一个有界的输入输出y(n)在n=0这一点上的值为第1章 离散时间信号与系统 也即y(0)是无界的,这不符合稳定的条件,因而假设不成立。所以是稳定的必要条件。要证明一个系统不稳定,只需找一个特别的有界输入,如果此时能得到一个无界的输出,那么就一定能判定一个系统是不稳定的。但是要证明一个系统是稳定的,就不能只用某一个特定的输入作用来证明,而要利用在所有有界输入下都产生
31、有界输出的办法来证明系统的稳定性。第1章 离散时间信号与系统 显然,既满足稳定条件又满足因果条件的系统,即稳定的因果系统是最主要的系统。这种线性时不变系统的单位脉冲响应应该既是因果的(单边的)又是绝对可和的,即:这种稳定因果系统既是可实现的,又是稳定工作的,因而这种系统正是一切数字系统设计的目标。第1章 离散时间信号与系统 1.3.7常系数线性差分方程常系数线性差分方程连续时间线性时不变系统的输入输出关系常用常系数线性微分方程表示,而离散时间线性时不变系统的输入输出关系除了用式(1-42)表示外,常用以下形式的常系数线性差分方程表示,即(1-52)所谓常系数是指决定系统特征的a1,a2,aN,
32、b1,b2,bM都是常数。若系数中含有n,则称为“变系数”线性差分方程。差分方程的阶数等于未知序列(指y(n))变量序号的最高值与最低值之差。例如,式(1-52)即为N阶差分方程。第1章 离散时间信号与系统 所谓线性是指各y(n-k)以及各x(n-k)项都只有一次幂且不存在它们的相乘项(这和线性微分方程是一样的);否则就是非线性的。离散系统的差分方程表示法有两个主要的用途,一是从差分方程表达式比较容易直接得到系统的结构,二是便于求解系统的瞬态响应。求解常系数差分方程可以用离散时域求解法,也可以用变换域求解法。第1章 离散时间信号与系统 离散时域求解法有两种:(1)迭代法,此法较简单,但是只能得
33、到数值解,不易直接得到闭合形式(公式)解答。(2)卷积计算法,这用于系统起始状态为零时的求解。第1章 离散时间信号与系统 变换域求解法与连续时间系统的拉普拉斯变换法相类似,它采用Z变换方法来求解差分方程,这在实际使用上是简单而有效的。卷积方法,前面已经讨论过了,只要知道系统脉冲响应就能得知任意输入时的输出响应。Z变换方法将在后面几节中讨论。这里仅简单讨论离散时域的迭代解法。差分方程在给定的输入和给定的初始条件下,可用递推迭代的办法求系统的响应。如果输入是(n)这一特定情况,输出响应就是单位脉冲响应h(n)。例如,利用(n)只在n=0取值为1的特点,可用迭代法求出其单位脉冲响应h(0),h(1)
34、,h(n)值,下面举例说明。第1章 离散时间信号与系统 例例1-3常系数线性差分方程输入为x(n)=(n)初始条件为y(n)=0n0,可得n0时h(n)=y(n)=0,将式(1-53)改写为另一种递推关系y(n-1)=2y(n)-x(n)或y(n)=2y(n+1)-x(n+1)又利用已得出的结果h(n)=0(n0),则有:第1章 离散时间信号与系统 所以也可表示为这样的系统是非因果系统,而且是非稳定的。第1章 离散时间信号与系统 1.4.1Z变换的定义及收敛域变换的定义及收敛域1.Z变换的定义变换的定义一个离散序列x(n)的Z变换定义为1.4Z变变换换式中,z是一个复变量,它所在的复平面称为Z
35、平面。我们常用Zx(n)表示对序列x(n)进行Z变换,也即(1-54)(1-55)第1章 离散时间信号与系统 这种变换也称为双边Z变换,与此相应的单边Z变换的定义如下:这种单边Z变换的求和限是从零到无穷,因此对于因果序列,用两种Z变换定义计算出的结果是一样的。单边Z变换只有在少数几种情况下与双边Z变换有所区别。比如,需要考虑序列的起始条件,其他特性则都和双边Z变换相同。本书中如不另外说明,均用双边变换对信号进行分析和变换。(1-56)第1章 离散时间信号与系统 2.Z变换的收敛域变换的收敛域显然,只有当式(1-54)的幂级数收敛时,Z变换才有意义。对任意给定序列x(n),使其Z变换收敛的所有z
36、值的集合称为X(z)的收敛域。按照级数理论,式(1-54)的级数收敛的充分必要条件是满足绝对可和的条件,即要求(1-57)第1章 离散时间信号与系统 要满足此不等式,|z|值必须在一定范围之内才行,这个范围就是收敛域。一般收敛域用环状域表示,即Rx-|z|Rx+收敛域是分别以Rx-和Rx+为半径的两个圆所围成的环状域(图1-21中的斜线部分)。Rx-和Rx+称为收敛半径。当然Rx-可以小到零,R x+可以大到无穷大。常用的Z变换是一个有理函数,用两个多项式之比表示:第1章 离散时间信号与系统 1-21第1章 离散时间信号与系统 分子多项式P(z)的根是X(z)的零点,分母多项式Q(z)的根是X
37、(z)的极点。在极点处Z变换不存在,因此收敛域中没有极点,收敛域总是用极点限定其边界。Z平面上收敛域的位置,或者说Rx-及Rx+的大小和序列有着密切的关系,分别讨论如下。(1)有限长序列:序列x(n)只在有限区间n1nn2之内才具有非零的有限值,在此区间外,序列值皆为零。也即第1章 离散时间信号与系统 其Z变换为设x(n)为有界序列,由于X(z)是有限项级数之和,除0与两点是否收敛与n1、n2取值情况有关外,整个Z平面均收敛。如果n10,则收敛域不包括z=0点;如果是因果序列,收敛域包括点。具体有限长序列的收敛域表示如下:(1-58a)(1-58b)(1-58c)有时将开域(0,)称为“有限Z
38、平面”。第1章 离散时间信号与系统 例例1-5x(n)=(n),求此序列的Z变换及收敛域。解解这是n1=n2=0时有限长序列的特例,由于所以收敛域应是整个z的闭平面(0|z|),如图1-22所示。第1章 离散时间信号与系统 图1-22(n)的收敛域(全部Z平面)第1章 离散时间信号与系统 例例1-6求矩形序列x(n)=RN(n)的Z变换及其收敛域。解解这是一个有限项几何级数之和。因此第1章 离散时间信号与系统(2)右边序列:右边序列是指x(n)只在nn1时有值,在nn1时x(n)=0。其Z变换为(1-59)此式右端第一项为有限长序列的Z变换,按上面讨论可知,它的收敛域为有限Z平面;而第二项是z
39、的负幂级数,按照级数收敛的阿贝尔(N.Abel)定理可推知,存在一个收敛半径Rx-,级数在以原点为中心,以Rx-为半径的圆外任何点都绝对收敛。因此,综合此二项,只有二项都收敛时级数才收敛。所以,如果Rx-是收敛域的最小半径,则右边序列Z变换的收敛域为Rx-|z|右边序列及其收敛域如图1-23所示。第1章 离散时间信号与系统 图1-23右边序列及其收敛域(n10,|z|=除外)第1章 离散时间信号与系统 因果序列是最重要的一种右边序列,即n1=0的右边序列。也就是说,在n0时x(n)有值,n|a|这是一个无穷项的等比级数求和,只有在|az-1|a|处收敛如 图 1-24所 示。故 得 到 以 上
40、 闭 合 形 式 的 表 达 式,由 于 ,故在z=a处有一极点(用“”表示),在z=0处有一个零点(用“”表示),收敛域为极点所在圆|z|=|a|的外部。第1章 离散时间信号与系统 收敛域上函数必须是解析的,因此收敛域内不允许有极点存在。所以,右边序列的Z变换如果有N个有限极点z1,z2,zN存在,那么收敛域一定在模值为最大的这一个极点所在圆以外,也即对于因果序列,处也不能有极点。第1章 离散时间信号与系统 第1章 离散时间信号与系统(3)左边序列:左边序列是指在nn2时x(n)有值,而在nn2时x(n)=0,其Z变换为等式第二项是有限长序列的Z变换,收敛域为有限Z平面;第一项是正幂级数,按
41、阿贝尔定理,必存在收敛半径Rx+,级数在以原点为中心,以Rx+为半径的圆内任何点都绝对收敛。如果Rx+为收敛域的最大半径,则综合以上两项,左边序列Z变换的收敛域为如果n20,则式(1-61)右端不存在第二项,故收敛域应包括z=0,即|z|Rx+。第1章 离散时间信号与系统 例例1-8x(n)=-anu(-n-1),求其Z变换及收敛域。解解这是一个左边序列。其Z变换为此等比级数在|a-1z|1,即|z|Rx-;第二项为左边序列,其收敛域为|z|Rx+。如果Rx-Rx+,则存在公共收敛区域,X(z)有收敛域Rx-|z|Rx+,则无公共收敛区域,X(z)无收敛域,也即在Z平面的任何地方都没有有界的(
42、z)值,因此就不存在Z变换的解析式,这种变换就没有什么意义。第1章 离散时间信号与系统 例例1-9x(n)=a|n|,a为实数,求其Z变换及收敛域。解解这是一个双边序列,其Z变换为设第1章 离散时间信号与系统 若|a|1,则存在公共收敛域其序列及收敛域如图1-26所示。若|a|1,则无公共收敛域,因此也就不存在Z变换的封闭函数,这种序列如图1-27。序列两端都发散,显然这种序列是不现实的序列。第1章 离散时间信号与系统 图1-26双边序列及收敛域第1章 离散时间信号与系统 图1-27Z变换无收敛域的序列第1章 离散时间信号与系统 表表1-1几种序列的几种序列的Z变换变换第1章 离散时间信号与系
43、统 1.4.2Z反变换反变换已知函数X(z)及其收敛域,反过来求序列的变换称为Z反变换,表示为x(n)=Z-1X(z)Z反变换的一般公式为若(1-63)则(1-64)第1章 离散时间信号与系统 图1-28围线积分路径第1章 离散时间信号与系统 证证该积分路径c在半径为R的圆上,即 z=Rej Rx-RRx+则(1-65)第1章 离散时间信号与系统 这个积分公式(1-65)也称为柯西积分定律。因此或直接计算围线积分是比较麻烦的,实际上,求Z反变换时,往往可以不必直接计算围线积分。一般求Z反变换的常用方法有三种:围线积分法(留数法)、部分分式展开法和幂级数展开法。第1章 离散时间信号与系统 1.围
44、线积分法(留数法)围线积分法(留数法)这是求Z反变换的一种有用的分析方法。根据留数定理,若函数F(z)=X(z)zn-1在围线c以内有K个极点zk,而在c以外有M个极点zm(M、K为有限值),则有(1-66)或(1-67)第1章 离散时间信号与系统 ResX(z)zn-1,zk表示函数F(z)=X(z)zn-1在极点z=zk上的留数。式(1-66)表示函数F(z)沿围线c反时针方向的积分等于F(z)在围线c内部各极点的留数之和。式(1-67)说明,函数F(z)沿围线c顺时针方向的积分等于F(z)在围线c外部各极点的留数之和。由式(1-66)及式(1-67),可得(1-68)将式(1-66)及式
45、(1-67)分别代入式(1-64),可得:(1-69a)(1-69b)第1章 离散时间信号与系统 根据具体情况,既可以采用式(1-69a),也可以采用式(1-69b)。例如,如果当n大于某一值时,函数X(z)zn-1在围线的外部可能有多重极点,这时选c的外部极点计算留数就比较麻烦,而通常选c的内部极点求留数则较简单。如果当n小于某一值时,函数X(z)zn-1在围线的内部可能有多重极点,这时选用c外部的极点求留数就方便得多。第1章 离散时间信号与系统 现在来讨论如何求X(z)zn-1在任一极点zr处的留数。设zr是X(z)zn-1的单一(一阶)极点,则有(1-70)如果zr是X(z)zn-1的多
46、重极点,如l阶极点,则有(1-71)第1章 离散时间信号与系统 例例1-10已知求Z反变换。解解围线c以内包含极点a,如图1-29粗线所示。当n|a|第1章 离散时间信号与系统 式中,a是单阶极点。应用公式(1-70),则在z=0处有一个-n阶极点(n0),应用公式(1-71),则因此第1章 离散时间信号与系统 即这个指数因果序列是单阶极点的反变换,这个反变换是很典型的,在以下的部分分式中还要用到这个结果。实际上,由于收敛域在函数极点以外,并且包括点,因此可以知道该序列一定是因果序列。用留数法计算的结果也证实了这一点。所以,在具体应用留数法时,若能从收敛域判定序列是因果的,就可以不必考虑n0时
47、围线c内无极点;而n0时只在z=0处有一个-n阶极点。因此第1章 离散时间信号与系统 即上例中,在n0时,也可用围线外极点a的留数来求,见式(1-69b),则有即从收敛域在函数极点所在圆以内,就能判断序列是左边序列,计算出来结果也证实了这个结论。第1章 离散时间信号与系统 第1章 离散时间信号与系统 2.部分分式展开法部分分式展开法在实际应用中,一般X(z)是z的有理分式,可表示成X(z)=P(z)/Q(z),P(z)及Q(z)都是实系数多项式,且没有公因式。可将X(z)展开成部分分式的形式,然后利用表1-1的基本Z变换对的公式求各简单分式的Z反变换(注意收敛域),再将各个反变换相加起来,就得
48、到所求的x(n)。为了看出如何求得部分分式展开,假设X(z)可以表示成z-1的多项式之比,即(1-72)第1章 离散时间信号与系统 为了得到X(z)的部分分式,将上式进一步展开成以下形式:式中,ck是X(z)的非零零点,dk是X(z)的非零极点。如果Mmax|dk|,因此上式部分分式展开式中每一项都是一个因果序列的z函数,可以直接利用例1-10的结果,得(1-75)式中,系数Ak可利用留数定理求得(1-76)如果MN,且除一阶极点外,在z=di处还有s阶极点,则X(z)可展开成(1-77)第1章 离散时间信号与系统 式中,Bn可用长除法求得。Ak可由式(1-76)求出。系数Cm由下式得到(1-
49、78a)或(1-78b)展开式各项被确定后,再分别求右边各项的Z反变换,则原序列就是各项的反变换序列之和。第1章 离散时间信号与系统 例例1-12设试利用部分分式法求Z反变换。解解X(z)有两个极点,d1=2和d2=0.5,收敛域为|z|2,则X(z)的零极点如图1-31所示。由收敛域可知x(n)是一个右边序列。因为极点全部是一阶的,因此X(z)能表示为:由用式(1-76)求得系数为:第1章 离散时间信号与系统 因此X(z)为根据表1-1可得或表示为第1章 离散时间信号与系统 例例1-13在这个例子中要考虑例1-12中给出的X(z)所对应的全部可能序列。解解根据图1-31的零极点图和收敛域性质
50、,X(z)有三种不同的收敛域:(1)|z|2,如例1-12,情况1已经证明是一个右边序列。(2),情况2对应于一个左边序列。(3),情况3则对应于一个双边序列。第1章 离散时间信号与系统 第1章 离散时间信号与系统 因为X(z)的部分分式展开仅决定于X(z)的代数式,所以对所有三种情况都是一样的。针对X(z)的三种不同的收敛域,根据表1-1可得:情况1:情况2:情况3:第1章 离散时间信号与系统 3.幂级数展开法(长除法)幂级数展开法(长除法)因为x(n)的Z变换定义为z-1的幂级数,即所以只要在给定的收敛域内,把X(z)展成幂级数,则级数的系数就是序列x(n)。把X(z)展成幂级数的方法很多