第二章连续系统的时域分析

1、1第二章第二章 LTILTI连续系统的时域分析连续系统的时域分析 21 系统的微分算子方程与传输算系统的微分算子方程与传输算子子一、微分算子、积分算子与微分算子方程一、微分算子、积分算子与微分算子方程:引入如下算子：引入如下算子： 微分算子微分算子: tp dd 积分算子积分算子: tpp 1 d) (1 )()(dd)( tfptfttf 则：则：)()(dd)( )(tfptfttfnnnn )()(1d )(1 tfptfpft 2对于微分方程对于微分方程 )(4d)(d)(6d)(d5d)(d 22tfttftyttytty 算子形式算子形式)(4)()(6)(5)( 2tftfpty

2、typtyp 微分算子方程：微分算子方程： )() 4()() 65(2tfptypp 它是微分方程的一种表示，含义是在等式两边它是微分方程的一种表示，含义是在等式两边分别分别对变量对变量y(t)和和f(t)进行相应的微分运算进行相应的微分运算。形式上。形式上是代数方程的表示方法。可用来在时域中建立与是代数方程的表示方法。可用来在时域中建立与变换域变换域相一致的分析方法。相一致的分析方法。3微分算子的运算性质：微分算子的运算性质：性质性质1 1 以以p的正幂多项式出现的运算式，在形式的正幂多项式出现的运算式，在形式上可以像代数多项式那样进行展开和因式分解。上可以像代数多项式那样进行展开和因式分

3、解。)()4()()3)(2(tfptypp 性质性质2 2 设设A(p)和和B(p)是是p的正幂多项式，则的正幂多项式，则 )()()()()()(tfpApBtfpBpA 如：如：)()4()()2)(3(tfptypp 性质性质3 3 微分算子方程等号两边微分算子方程等号两边p的公因式不能的公因式不能 随便消去随便消去。 例如：例如：p y(t)= p f(t) y(t)= f(t)+c( (c为常数为常数) ) y(t)= f(t) 性质性质4 4 设设A(p)、B(p) 和和D(p)都是都是p的正幂多项式的正幂多项式4)()()()()()()()(tfpBpAtfpBpDpApD

4、)()()()()()()()(tfpBpAtfpDpBpDpA 但是但是 ：)(d)(dd)(1 tffttfppt )()()(d)(dd)(1 tfftfftfppt 例如：例如： 函数乘、除算子函数乘、除算子p的顺序不能随意颠倒，的顺序不能随意颠倒，对函对函数进行数进行“先除后乘先除后乘”算子算子p的运算的运算时，分式的分时，分式的分子与分母中子与分母中公共公共p算子算子( (或或p算式算式) )才允许消去才允许消去。5二、二、LTILTI连续系统的算子方程与系统的传输算子连续系统的算子方程与系统的传输算子 电路元件伏安关系电路元件伏安关系( (VAR) )的微分算子形式称为的微分算子

5、形式称为 算子模型算子模型，电压、电流比为，电压、电流比为算子感抗算子感抗和和算子容抗算子容抗 元件名称 电路符号 ui关系(VAR) VAR的算子形式 算子模型 电阻 电感 电容 电路元件的算子模型电路元件的算子模型 i(t) R)(tui(t) R)(tui(t)L)(tu)(tui(t)1/pC)(tui(t)Ci(t)pL)(tuttiLtu d)(d)( tdiCtu )(1)( )(1)(tipCtu )( )(tiRtu )( )(tiRtu )( )(tipLtu 6电路系统微分算子方程的建立方法电路系统微分算子方程的建立方法: : LpL;C 1/pC画出算子模型，按照电路理

6、论画出算子模型，按照电路理论中的列写方程方法列写。中的列写方程方法列写。例例1 1：电路如图电路如图( (a) )所示，激励为所示，激励为f(t)，响应为，响应为i2(t)。试列写其微分算子方程。试列写其微分算子方程。(a)1+f(t)-i153Fi22H4H1+ +f(t)- -i15 1 3pi22p4p(b)i1i2解：解：画出其画出其算子模型电路算子模型电路如如图图( (b) )所示。由所示。由回路回路法法可列出方程为可列出方程为 ：7 0)()5431()(31)()(31)()3121(2121tipptiptftiptipp 化简微分方程组化简微分方程组时要时要考察电路的阶数考察

7、电路的阶数以便确定以便确定公共因子是否可消去。公共因子是否可消去。化简后化简后所求微分算子方程为：所求微分算子方程为： )()() 27148( 3223tftippp 对于激励为对于激励为f(t)，响应为，响应为y(t)的的n阶阶LTI连续系统，连续系统，其微分算子方程为：其微分算子方程为：)()()()(01110111tfbpbpbpbtyapapapmmmmnnn 8将其在形式改写为将其在形式改写为)()()()(01110111tfpHtfapapapbpbpbpbtynnnmmmm )()()()()( 01110111pDpNapapapbpbpbpbtftypHnnnmmmm

8、式中：式中： 它代表了系统将激励转变为响应的作用，或它代表了系统将激励转变为响应的作用，或系统对输入的传输作用，故将系统对输入的传输作用，故将H(p)称为称为响应响应y y( (t t) )对激励对激励f f( (t t) )的传输算子的传输算子或或系统的传输算子系统的传输算子 系统传输算子与系统微分算子方程是对系统系统传输算子与系统微分算子方程是对系统的等价表示。它们之间可以可以转化。的等价表示。它们之间可以可以转化。922 LTI22 LTI连续系统的零输入响连续系统的零输入响应应 LTILTI的全响应可作如下分解：的全响应可作如下分解： 1、y(t) = 自由响应自由响应 + 强制响应；

9、强制响应； 2、y(t) = 零输入响应零输入响应yx(t) + 零状态响应零状态响应yf (t) 零输入响应零输入响应: :是指输入激励为零，仅由系统的是指输入激励为零，仅由系统的初始状态单独作用而产生的输出响应。初始状态单独作用而产生的输出响应。零状态响应零状态响应: :是指系统的初始状态为零，仅由是指系统的初始状态为零，仅由系统的输入激励单独作用而产生的输出响应。系统的输入激励单独作用而产生的输出响应。10二、通过系统微分算子方程求零输入响应二、通过系统微分算子方程求零输入响应零输入下零输入下LTI连续系统的微分算子方程为连续系统的微分算子方程为:1110 x()( )0,0nnnpap

10、a p a y tt要使上式成立，需满足要使上式成立，需满足D(p)=0（特征方程）（特征方程） 针对针对特征根特征根两种情况来求两种情况来求yx(t) 1特征根为特征根为n个单根个单根p1 , p2 , , pn ( (可为实根、虚可为实根、虚根或复根根或复根) ) 0 , eee)( 2121x tAAAtytptptpnn 将将yx(0-)、yx (0-)、yx(n-1)(0-)代入上式，确定代入上式，确定积分常数积分常数A1、A2、An 。 共轭复根时欧拉公式共轭复根时欧拉公式cos t = 0.5(ej t + e j t )及及sin t = j0.5(e j t ej t )化简

11、为三角化简为三角实函数实函数 112 2特征根含有重根特征根含有重根 设特征根设特征根p1为为r重根，其余特征根为单根，重根，其余特征根为单根，, , , , 2 1nrrppp 则则yx(t)的通解表达式为：的通解表达式为：0 ,ee )()( 1111 2321x tAAetAtAtAAtytptptpnrnrrr确定积分常数的方法同前。确定积分常数的方法同前。 123求解零输入响应求解零输入响应yx(t)的基本步骤：的基本步骤： ( (1) )通过微分算子方程得通过微分算子方程得D(p)求系统的特征根求系统的特征根; ; ( (2) )写出写出yx(t)的通解表达式的通解表达式; ; (

12、 (3) )由系统的由系统的0-状态值与状态值与0-瞬时的零输入系统求得瞬时的零输入系统求得初始条件初始条件yx(j )(0-), j=0, 1, 2, , n-1。(4) 将将0-初始条件代入初始条件代入yx(t)的通解表达式的通解表达式,求得积分求得积分常数常数A1, A2, , An 。( (5) ) 写出所得的解写出所得的解yx(t)，画出，画出yx(t)的波形。的波形。 13例：已知某线性时不变系统的输入输出方程为系统的初始状态为求系统的零输入响应。14例例2 电路如图电路如图( (a) )所示，已知所示，已知uC (0-) = 1V，iL(0-) = -1A，求，求t0时的零输入响

13、应时的零输入响应uCx(t)。1H12F CuCi 21R 42RLi CuCi 2 4LiP2P1解解 (1)画出算子模型电路画出算子模型电路, ,由节点法列出方程为由节点法列出方程为 0)()41212( tuPPcx15uC x (t), V0t, s4130.5 1化简可得化简可得 :0)()65(x2 tuppC解得特征根解得特征根: : p1=-2，p2=-3 0 ,ee)( 32 21x tAAtuttCV1A124(2)0-瞬时的等效电路瞬时的等效电路 sV1)0(1)0(21)1(21)0(x x x CCCiCui 343211212121AAAAAA代入初始条件代入初始条

14、件. 0 ,V34)( 3 2x teetuttC16作业：2-22-3 （3）（4）2-6 （a）（b）1723 LTI23 LTI连续系统的零状态响应连续系统的零状态响应 一、一、零状态响应零状态响应 零状态LTI连续系统H(p)(tf)(tyf)()()()()()(tfpDpNtfpHtyf )( )()()()(非非齐齐次次微微分分方方程程tfpNtypDf 一般情况下一般情况下零状态响应零状态响应可通过将可通过将f(t)分解为分解为更为简单的单元信号更为简单的单元信号，将各，将各单元激励下的响应进单元激励下的响应进行叠加行叠加来求解。来求解。18 任意信号任意信号f( (t) )用

15、一系列宽度为用一系列宽度为 的矩形窄脉冲近似。的矩形窄脉冲近似。kktGkftf)()()(信号的时域分解：信号的时域分解：19dtfktGkftfk)()( )(1)(lim)(0 任意信号任意信号f( (t) )可以可以分解为不同时刻出现的分解为不同时刻出现的受该时刻受该时刻f( (t) )加权的门信号的无穷级数和。加权的门信号的无穷级数和。任意信号可分解为无穷多个不同时刻出现的冲任意信号可分解为无穷多个不同时刻出现的冲激强度为该时刻函数值的冲激信号之和激强度为该时刻函数值的冲激信号之和 dtftf)()()(0 20零状态响应的求解过程零状态响应的求解过程零状态零状态LTI)(t )(t

16、h零状态零状态LTI)( t)( th零状态零状态LTI)()( tf)()( thf零状态零状态LTI dtftf)()()( dthftyf)()()( 冲激响应冲激响应时不变性时不变性齐次性齐次性叠加性叠加性21二、冲激响应h(t)h(t)定义定义： 零状态LTI H(p)(t )(th单位冲激响应：系统在单位冲激函数激励下的零状态响应，简称冲激响应，统一以符号h(t) 表示。22)()()()()()()(01110111tapapapbpbpbpbtpDpNtpHthnnnmmmm 据据D D( (p p) )的根的不同有理真分式的根的不同有理真分式H(p)可展开为不可展开为不同的部

17、分分式同的部分分式 1当当D D( (p p) ) 有有n个单特征根个单特征根p1 , p2 , , pn ( (可为实可为实根、虚根或复根根、虚根或复根) ) 23)()()()()()(21npppppppNpDpNpHnnjjppKppKppKppK 2211)()()()()(2211tppKtppKtppKtppKthnnjj ),()(tppKthjjj 令第令第j项为项为 )()()(tKthppjjj )()()(tKthpdttdhjjjj 24)()()(tKethpedttdhjtpjjtpjjj tpje )()(tKethdtdjtpjj tjttpjtKethdtd

18、j 0 0 )()( )()0()()(tKhethethjjtpjtpjjj t0冲激响应h(t)为为)(e)(e)(e)( 2 121tKtKtKthtpntptpn 252当当D D( (p p) )特征根有重根时：特征根有重根时：设设p1为为r重根，其余重根，其余(n-r)个为单根个为单根pj( (j=r+1, r+2, , n) )，则有理真分式，则有理真分式H(p)可展开为：可展开为：)()()()()()()(11nrrpppppppNpDpNpH nnrrrrrppKppKppKppKppK 11111112111)()(26重根相关的部分分式项的冲激响应重根相关的部分分式项的

19、冲激响应 rmttmrKthrmtpmrm , , 2 , 1 , )(e)!()(11 13 3、H H( (p p) )为某个关于为某个关于p pj j多项式时：多项式时：rjtpkthjrjj , , 2 , 1 , )()(1 rjtkthjrjj , , 2 , 1 , )()()(1 27求解单位冲激的步骤：求解单位冲激的步骤：（1）据算子微分方程求出传输算子）据算子微分方程求出传输算子H(p)（2）长除法化为长除法化为多项式与有理真分式之和。多项式与有理真分式之和。（3）有理真分式部分分式展开；）有理真分式部分分式展开；（4）据）据D(p)根的不同根的不同确定确定分式中的分式中的

20、系数；系数；（5）对照不同情况写出单位冲激响应。表）对照不同情况写出单位冲激响应。表2-2282379972)(2234 pppppppH例：求系统的单位冲激响应：例：求系统的单位冲激响应：29三、阶跃响应g(t)g(t)定义定义： 零状态LTI H(p)(t )(tg根据根据LTI系统的积分性，通过对冲激响应系统的积分性，通过对冲激响应h(t)进行积分而求得：进行积分而求得： ttdhtgdt00)()()()( 30作业：2-7 （2）（3）2-8 （b）31三三 卷积积分卷积积分 (*) )()()( dthfty (*) )()()()()(2121 dtfftftfty 定义定义:卷

21、积积分简称卷积积分简称卷积卷积. .32卷积积分上下限的确定是关键，讨论如下：卷积积分上下限的确定是关键，讨论如下：(1)若若f(t),h(t) 都为因果信号都为因果信号积分上下限积分上下限为为(0-, t)(*) )()()()()( 0 dthfthtftyt (2)若若f(t) 为因果信号为因果信号,h(t) 为无时限信号为无时限信号,积分上下积分上下限限为为(0-, )(*) )()()()()( 0 dthfthtfty (3)若若f(t) 为无时限信号为无时限信号,h(t) 为因果信号为因果信号,积分上下积分上下限限为为(- , t)(*) )()()()()( dthfthtft

22、yt (4)若若f(t), h(t)都为时限信号则都为时限信号则卷积后仍为时限信卷积后仍为时限信号，其号，其左边界为原两左边界之和，右边界为原两左边界为原两左边界之和，右边界为原两右边界之和右边界之和 33（1 1）卷积的运算规律）卷积的运算规律 据卷积的定义和积分的性质，可推知卷积有如据卷积的定义和积分的性质，可推知卷积有如下的运算规律下的运算规律 :1 1交换律交换律: : )()()()(1221tftftftf* 2 2分配律分配律: : )()()()()()()(3121321tftftftftftftf* 3 3结合律结合律 )()()()()()()()()(231321321

23、tftftftftftftftftf* 34（2）卷积的主要性质）卷积的主要性质1 1f( (t t) )与奇异信号的卷积与奇异信号的卷积(1)(1) f(t)* *(t)=f(t),即即f(t)与与(t)卷积等于卷积等于f(t)本本身身 (2)(2) f(t)* *(t)=f(t) ,即即f(t)与与(t)卷积等于卷积等于f(t)导数。导数。 (3)(3)()()()(1tfdfttft 2 2卷积的微分和积分：卷积的微分和积分：(1)(1) 积分积分 f1(t)* *f2(t) -1 = f1-1(t)* *f2(t)= f1(t)* *f2-1(t) 35(3)(3) 微分微分- -积分

24、积分: :f1(t)* *f2(t)=f1(t)* *f2-1(t)=f1-1(t)* *f2(t) (2)(2) 微分微分 f1(t)* *f2(t) = f1(t)* *f2(t)= f1(t)* *f2(t) 0)()()()(1221 dttffdttff3 3卷积时移：卷积时移：设设f1(t)* *f2(t)=y(t)，则：，则： f1(t)* *f2(t-t0)=f1(t-t0)* *f2(t)=y(t-t0) f1(t-t1)* *f2(t-t2)=y(t-t1-t2); 推论：推论：f(t-t1)* *(t-t2)=f(t-t1-t2) (t-t1)* *(t-t2)=(t-t

25、1-t2); 利用卷积性质求解较复杂的卷积利用卷积性质求解较复杂的卷积 （表表2-3）)()()()(212121tttrtttttttt 条件：条件：36例例7:例例3已知已知:)1()( 2)(1 tttf )2()()(2 ttttf 解：解：1.把信号写成标准的延时信号把信号写成标准的延时信号2.分配律写出各项卷积分配律写出各项卷积3.查查P45表表2-337012313y(t)t38 若若f1(t)，f2(t)收敛，收敛，将被卷积的一个信号尽量将被卷积的一个信号尽量化为化为冲激信号以及其延时冲激信号以及其延时，可使计算简化。，可使计算简化。)2()()1(2)(2)()( )()(1

26、2121 ddtttftftftftt )2()221()(21)1(2)(2 22 tttttt )3( 4) 1()2()4()( 222 tttttt 39例例8 8 试计算常数试计算常数K K与信号与信号f( (t t) )的卷积积分的卷积积分 解解 直接按卷积定义，可得直接按卷积定义，可得 ：)( )()()( 下下的的净净面面积积tfKKdfKtftfK 用微分用微分- -积分性质来求解将积分性质来求解将导致错误结果导致错误结果 0)(dd)( tdfKttfK 常数常数K 不收敛不收敛且任意信号且任意信号f(t)也并非一定也并非一定收敛。收敛。 40例例9 9 已知某系统的冲激响

27、应已知某系统的冲激响应h(t)=sint (t)，激，激励励f(t)的波形如图所示，试求系统的零状态响的波形如图所示，试求系统的零状态响应应yf(t)。可用微分可用微分- -积分性来求积分性来求)()cos1 ()(sin)( 01ttdtht 解：解： 系统的零状态响应求解系统的零状态响应求解f(t)0t24241)()sin()()cos1()( 02tttdtht )4()2(2)()( ttttf)4()2(2)()( ttttf)4()4sin()4( )2()2sin()2(2)()sin( )4()2(2)()()sin( )()( )()( )()()(21 ttttttttt

28、tttttttfthtfthtfthty*f”(t)0t24（1）（2）（1）42+-f(t)i(t)uc(t)+-p1/p例例10:图示电路图示电路,激励激励求求:零状态响应零状态响应uc(t)6()()( tttf)(11)(11)(2tfptfppptuc 解：解：列方程列方程+-f(t)i(t)uc(t)+-1H1F4311)()()(2 ptftupHc)()sin()(ttth )()sin(*)6()( )(*)()(ttttthtftuc ttdttdtt0sin*)6()()(sin*)6()( )6()6cos(1 )()cos1( tttt44图示电路，其输入电压图示电路

29、，其输入电压us(t)波形如图波形如图示，试用卷积积分法求零状态响应示，试用卷积积分法求零状态响应uc(t)0.1M 10Fu uc c( (t t) )u us s( (t t) )11)()()( ptutupHsc)()(tetht 解：解：u us( (t t)(V)(V)t t( (s s) )3210145)(*)3()1()1()()(tettttttutc detttttdtdt)(* )3()1()1()(0 )()1(*)3()1()(tetttt )3(1)()1(*)1()()3( tetetttt 利用卷积的性质利用卷积的性质46)3(1 )1(11)()1()3()

30、1( tetettetttt )3(1 )1(2)()1()3()1( tetettetttt 47四、系统全响应的求解方法：四、系统全响应的求解方法：（1）求单位冲激响应）求单位冲激响应h(t)（2）求卷积积分）求卷积积分 dthf)()( （3 3）求零输入响应）求零输入响应yX (t) 零状态响应零状态响应yf (t)（4）全响应：）全响应：)()()(xtytytyf 48例例11 图示电路已知图示电路已知i1(0-) = i2(0-) =1A， f1(t) = t (t)，f2(t) = (t)- (t-1)，求全响应，求全响应y(t) 。1 i1(t)+ +f1(t)- -+ +f2(t)- -1 1 + +y(t)- -i2(t)1H1H解：使用叠加原理求解解：使用叠加原理求解49211 21pppp