连续时间系统的时域分析报告

10101010..第二章连续时间系统的时域分析§引

言线性连续时间系统的时域分析，就是一个建立和求解线性微分方程的过程。一、建立数学模型主要应用《电路分析》课程中建立在KVL基础上的各种方法。线性时不变系统的微分方程的一般形式可以为：dndr()r(t)r(tr(t)nndtndtdmd(t)()()(tdtdtdt二、求解（时域解）1时域法将响应分为通解和特解两部分：1通解通过方程左边部分对应的特征方程所得到的特征频率，解得的系统的自然响应（或自由响应2特解：由激励项得到系统的受迫响应；............3代入初始条件，确定通解和特解中的待定系数。经典解法在激励信号形式简单时求解比较简单，但是激励信号形式比较复杂时求解就不容易，这时候很难确定特解的形式。2卷积法（或近代时域法，算子法）这种方法将响应分为两个部分，分别求解：1零输入响应：系统在没有输入激励的情况下，仅仅由系统的初始状态引起的响应r(t)

；2零状态响:状态为零（没有初始储能）的条件下，仅仅由输入信号引起的响应r(t)

。系统的零输入响应可以用经典法求解，在其中只有自然响应部分；系统的零状态响应也可以用经典法求解，但是用卷积积分法更加方便。借助于计算机数值计算以求出任意信号激励下的响数值解卷积法要求激励信号是一个有始信号，否则无..........pdtdtntLupdtdtntLuid..法确定初始状态。零输入响应与自然响应、零状态响应与受迫响应之间并不相等，具体对比见§2经典法在高等数学中已有详细介绍。本课程中重点介绍近代时域法。§2一、算子

系统微分方程的算子表示通过微分算子可以简化微分方程的表示。微分算子：令，，积分算子：

1p

())利用算子可以将电路中的电感和电容的伏安特性记为：LLdt

L1t1C可以感和容成为L和..........101010101

..的电阻，即感抗和容抗。利用算子可以将线性时不变系统的微分方程：dndr()r(t)r(tr(t)nndtndtdmd(t)()()(tdtdtdt表示为：nr(t)

n(t)pr(t)r(t)0

et)

(t)(t)t)按照代数运算法则，提取公因子，可以将上式简化为：p

pr(t)0(bpm

m)et)或进一步简化为：r()

(pm(

mm)10p)1

()定义：(p)则：

(m(n

mpp)(p)m1pn)Dp)10r(t)(p)(t)..........tt..注意上面只是微分方程的一种简单记法，并不代表能进行这样的计算。二、算子运算法则1

np

，其中m,n为任意常数。2pmpnpm,中同为任意正整数（或负整数3

1pp但是：

，1

1p

px

不一定等于

x

——微分和积分的次序不能交换：1ppdt

，但是：

即，1pxpxp2如果

p(t(t

,不一定能够推出............(t)t)

，只能得到

()yt)

，即等式两边的公共因子不能抵消。可见，大部分代数运算法则可以使用，但是有一些不能用。§2-3

系统的零输入响应零输入响应是下列齐次方程的解：()r(t)p

n

n

p

n

p)rt)1对它有两种解法：经典解法等效源法（或初始条件法）一、经典解法用经典法求解零输入响应由如下两步构成：1确定系统的自然频率：令D(p)=0，将p看成一个代数量，解得其n个特征根

1n

。2确定零输入响应的形式解：1如果D(p)=0俱单根时，则可以确定其形式解为：..........tktk22..nr(CeeenCezi1nii

i其中

CC2

为待定常数。2如果有重根时设是一个k重，即

2k

，则形式解为：r(t)et2etkzi3k

1k

tk

tnti

i

1

Cei

i

ti

i3根据初始条件，确定待定系数：一般的初始条件为已知零时刻的响应及其各阶导数

r(0),rr''(0),...,r((0)

入形式解中就可以确定待定系数。当俱为单根时：r(0)12r'(0)1122

nr''(0)CC12n

nr

(n

(0)1

(n

1

(n

2

(n

n由上面的个方程就可以确定n个待定系数。或者记为矩阵形式：..........222122212..

rrr''(0)

1123123

1Cnn

r

(

(0)

n

nnnn

n

nn其它形式的初始条件，以及特征方程中有重根的情况下的待定系数也可以用相似的方法和过程解出。举例：例．

已知系统的转移算子及未加激励时的初始条件是：

pp

r

，求系统的零输入响应并指出其自然频率。解p

，

，

，

。所以零输入响应的形式为re2

............：12r'

e

e

，r12r'0c12解得：1得，r

，自然频率分别为12二、等效源法

。这种方法将初始条件看成是一个在的瞬间加上的激励源（阶跃或冲激后将系统的初始条件归为零，从而将求解零输入响应的问题转化为求解零状态响应的问题。这在后面求解零状态响应中一并处理。等效源法将在时域解法中用得不是很多，本课程将在Ch5中介绍其原理。小结：系统零输入响应的时域求解方法就是经典的齐次微分方程的解法。............习题：①.2.3(1)-图②.2.4(3)。******************************************以下内容涉及到系统零状态响应的求解过程，为了叙述的清楚起见我们这里先简单了解其求解的基本思路。系统零状态响应的求解过程求解零状态响应的基本思想：将任意信号分解为一系列“标准统一”的子信号之和（或积分求线性系统对各个子信号的响应；将各子信号的响应相叠加从而得到系统对激励信号的响应这其中利用到了线性系统的齐次性和叠加性。为了求解线性系统的零状态响应，必须解决以下几个问题：选取什么样的子信号？如何将信号分解为子信号的和或积分？3如何求系统对子信号的响应？............4如何求得最后的响应？在下面的各节中我们将就上面的问题一一进行讨论。其中，§2介绍了时域分解法中使用的子信号；§2介绍如何将任意信号分解为子信号之和；§2介绍如何求子信号的响应；§2-7～§2介绍如何通过叠加而求出系统的响应。§奇异函数本节解决的是时域法中子信号选取问题。子信号的选取对系统分析至关重要。为了利于分析，要求子信号具有：1完备性：任意函数（或绝大部分函数）都可以分解为该子信号的和，没有（或几乎没有）例外；简单性：容易求得系统对该子信号的响应；相似性：不同子信号的响应具有内在联系，可以类推。奇异函数是一种理想化的函数，这些函数或其各阶导数具有一个或多个间断点，在这些间断............点上的导数无法用一般方法确定。常用的有阶跃函数和冲激函数。1阶跃函数

(t)

(t0

t其中t>0。1任意函数乘以

(t)

以后其部分等于零，成为有始函数。在很多文献中，用u(t)表示阶跃函数。2冲激函数

()冲激函数的图形表示方法：位置，强度。..........tt..其中>0。1冲激函数有很多种定义方法。常见的有两种：1定义为

(t)

的导数：

)

dt

(t)显然，该函数只在为非零值，其它各处都为零；

(t)

和

)

互为微分和积分

()

)dt)

的几个特性：

)

；

)t

时。

)

——冲激函数是一个偶函数

t)f(t))f

，或：

t)(t)f(t)00

ft)

t)dtf(0)或：

f(t()f(t)0——该特性被称为冲激函数的取样特性2定义为分配函数这是冲激函数的另外一种定义方法，它通过该函数对另外一个函数的作用来定义这个函数，利用上面的冲激函数的抽样特性作为冲激函数的............定义，即：对于任意的函数

使满足公式

f(t

()f(t)0性质的函数被称为冲激函数。参见本章附录。这种定义在数学上比较严格，但是难于理解。冲激函数的推广：冲激偶……举例：例1．计算f

。解：由于

(t)t)(t)f()000

，则有：(3)t

例2．计算f

。..........1t101t10..解：由于

f(t)

()f

，则有：

cos

t)。例3．计算

dt4

。解：

2t)cos3cos42

。例4．计算sin02

；解

100

sin2

。例5．计算

e

。解：原式

=2+1=3t

taa§信号的时域分解在近代时域法中使用的子信号是阶跃函数和冲激函数。本节讨论如何将信号分解为冲激函数和阶跃函数的和（或积分..........kk..一、任意有始函数表示为阶跃信号的和（积分）任意信号近似表示为阶跃函数的子信号为：f(f(0)(tf(1f(∶∶∶

(f(tffkf∶∶∶总和为：..........tt..f()()f()f()f(t)......a1ff(2)ff）)fknf(kff）)k令

d

(或)

，则：fkf

f(t)f(0)'(

这里，因为在t=0时发生状态的跳变，所以必须对左右的状态加以区分，这就引出了两个特殊的时刻：0——起始状态；0——初始状态；如果f(t)在=0处连续可导，上式可以简化为：f()

t

f

(

)d

这种分析方法在六十年代应用比较广泛，其时冲激函数没有得到应用。由于它需要计算函数的导数比较麻烦现在冲激函数应用以后，这种分解方法就很少应用了。............二、任意函数表示为冲激函数之和（积分）将任意函数近似表示为一系列矩性脉冲函数之和，定义：u(t)0

0则：..........

f(t(0)(t)0f(t)(u(t1f(t)f(2(t2…………f(t)f(k()tt..…………求和，可得：f(t)f(t)ft)ft)f(t)......01k(0)u()f(u(t......f(k(t)......f(0)

u()

f)

u(

......u()f()令

nd

f)()

u(t)，则

(t

(tf(t)f0

)d............这个公式实际上可以直接从冲激函数的定义或性质中推导出，但是上面的推导更利于观察其含义。这种分解不仅可以用于有始信号，也可以用于一般信号，这时候公式可以改写成为：f()

f这个公式更具有普遍性。公式中的积分上限也可以从t为这样不会影响结果。这时候公式为：f(t)f(0

)习题：①.2.5（。*******************************************§阶跃响应和冲激响应一、定义阶跃响应：系统对阶跃信号的零状态响应；冲激响应：系统对冲激信号的零状态响应；..........tt..

一般用

h()

表示系统的冲激响应，

rt

表示

阶跃响应；系统的冲激响应和阶跃响应之间有对应关系：(t)

dt

r(t

或

r()(0

)d

所以两者只要知道其一就可以了。如前所述，现在很少将信号分解为阶跃信号，所以一般没有必要求阶跃响应，只要求冲激响应就可以了。二、系统冲激响应的求解方法根据前面对信号的分解，有：f()

f

)d信号可以分解为多个冲激信号

f

)的积分。如果知道信号对

()

的响应，利用线性移不变系统的线性和移不变特性，就可以得到系统对任意子信号

f

)

的响应而就可以得到系统对整个信号的响应。求解系统冲激响应的方法有：1系统方程法：根据微分方程求解。书上页介绍了这种方法。............2系数平衡法：比较等式两边相同函数的系数，得到解答。书中的例题(2-3、2-4、2-5)都使用了这种方法。3初始条件法冲激激励转化0时刻的初始条件，然后利用零输入响应的求解方法求解。例题2-4中介绍了这种算法。4LT变换法利用拉普拉斯变换求解这种方法最简单。在后面Ch5中介绍。本节中重点介绍系统方程法。三、冲激响应的系统方程求解方法——系统方程法或Heaviside部分分式分解方法)1一阶系统的冲激响应的求解用算子表示的形式为：ht)

p

t)

，t)

，h()(t

)

，h'())

，微分方程两边同时乘以

e

，可以得到：'(t..........

t)t()dt()dt..=>

h'(e

ke

(t=>

dt

t)t)=>

0

d

0

)d=>

((0)

(t=>

ht(t

(意：零状态，

h

)或者简单记为：kh(t))p

t)2一般情况下，系统的特征根D(p)=0的根）无重根Npp)H(m1D(p)(p)n0

(ppp)1(p12n12npmpmm1一般使用Heaviside部分分式分解法，它将复杂系统变为许多个简单系统的和。1m<n时............借助于代数运算，通过部分分式求解，可以得到：kk)1p12n由此可以得到：h(tH())k2n)pp2nk1)pp2

)

kp

)e

(t)e

t

(t)

t

()keiii

)2当m=n时可以将H(p)分解为：kkk()pp12n(t)1p

))

(t

t

(t

t

(t

ki

i

(t)i其中b是m

3当m>n时............可以将H(p)分解为：Hpmpnk12pn则：

(t)C

()

())ke

1

()

t2

t)

n

()3一般情况下，系统的特征根D(p)=0的根）有重根假设，且1ll

l

nHp)1lp(2(kllnppll可以证明：

lk

)kte

t

()k

t22

()kt3et)p3..........12lt12lt..…………kt)p(

t

(则：tlh(t)kte(l

1

(tel

tl

(tel

tl

(t)en

tn

(t)有关和的情况，也可以通过相似的过程得到。例：，例题2-3。§2-7

卷积积分本节讨论如何通过冲激响应或阶跃响应求解系统对信号的响应。一、通过阶跃响应求解——杜阿美积分由§2-5节推导得到的公式：f(t)f(0)f'(

此时激励信号e的积分：

(t)可以分为一系列阶跃函数..........tttttttttttt..et)(0)

()

t

e'(

(

)而系统对阶跃信号的响应：

r()==>

(t

)rt

)

——时不变性==>

'(

)'(r(

)

——齐次性==>

e'(

(

)d'(

)rt

)d——叠加性

(t)

)

(t

)==>

0(0)r(t)0

)r(t

)所以：et)er()

e

)r(t

)d——杜阿美积分可见，如果得到了系统的阶跃响应，通过杜阿美积分，就可以计算出系统对任意连续可导的激励信号

(t)的响应

如果激励信号在t=0处可导，则上式为：et)

e(

通过变化积分变量，可以得到杜阿美积分的另外一种形式为：..........ttttttttt..et)e(0)r()

e'(t

)r(

)因为需要计算信号的导数，需要激励信号连续可导，所以这种方法目前不常用。二、通过冲激响应求解——卷积积分由§2-5节推导得到的公式：f(t)f0

)d此时，其激励信号e(t)可以分解一系列冲激函数的积分：et)

t

e

t

)d0系统对冲激信号的响应：

h(t)==>

)(t

)

——时不变性==>

e(

)e((t

——齐次性==>

)ht

)——叠加性所以：et)t0——卷积积分可见，如果得到了系统的冲激响应，通过卷积积分，就可以计算出系统对任意信号(t)的............应。与杜阿美积分相比，这里并不需要信号连续可导，所以其实用性大大优于杜阿美积分。

通过变化积分变量，同样可以得到卷积积分的另外一种形式为：et)

t0

et

以上公式的应用条件是：有始信号作用于因果系统。卷积积分有另外一种更加通用的形式是：et)

e

)h(

)d该公式的积分限有始信号作用于因果系统，与原公式的积分限等价。我们定义一种特殊的函数与函数之间的计算：卷积：()(tx(t则卷积积分可以表示为：r()(t)t)习题：①2.15；②2.16（*******************************************............§卷积及其性质一、卷积计算的几何解法et)

t0

e(t卷积积分的计算从几何上可以分为四个步骤：反褶——平移——相乘——>叠加（积分）............例题见P61-62。二、卷积计算的解析法根据卷积的定义求解。例题见。三、卷积积分表P60页四、卷积的性质卷积的计算不少类似于函数的乘法计算。它的很多性质与乘法运算性质相同，但是也有一些不同。通过这些性质，可以方便卷积的计算。1与乘法运算相同的性质：交换律：分配律：..........

u(t)t)()()tttttttttttt..(t)tt)(t(t)3结合律：())2与乘法运算不同的性质：卷积的微积分1微分：

()t)()(t(t)()

2积分：

u()(t)

)d

对上述两边求导数，得到：dt

)

)

(t)dt

(

)t)dt3多重微积分：u(t)3函数延时后的卷积

假设：

ut)()f(t)则：..........

(t)()f(t)2tt..五、几个特殊函数的卷积：1

f(t)

(t(t)或：

f(t)

(t)f(02

f(t)

'(t)f'()

(见分性质)或推广：

f()

(

(t)(n))003

f(t)

(tf(六、卷积积分的计算卷积积分的计算就是一般的定积分的计算。但是工程上遇到的卷积计算可能比较复杂。一般可以借助于图形帮助确定积分函数和积分边界。例：P69-71，例题。习题①2.19;②2.20（*******************************************............§2

线性系统响应的时域求解法一、近代时域法求解步骤求系统的转移算子H(p)求系统的零输入响应求解方法：经典法，等效源法如果系统的初始条件为零，则本步可以省略。3求系统的零状态响应求系统的冲激响应；通过卷积积分，求系统对激励信号的响应；如果积分难于计算，可以通过计算机数值积分计算。4将零输入响应与零状态响应相叠加，得到总响应。二、系统对指数激励信号(st(t)的响应若系统的特征方程无重根时，系统的全响应为：..........nnnntnnnnnnntnnn..r(t)r(t)(t)tzizsj

(t(j

tj

(t)

(tj

j)jjj

j

j

(ej)

(t零输响应

零状响应jj

j

j

jt

)je)jj自然应

受迫应jj

j

j

e

tj

)(est)受迫应自然应若系统的特征方程有重根时上式系统全响应的有关项则为：

rtCttk0响应信号按照其数学特性可以分为自然响应和受迫响应，也可以按照物理特性分为零输入响应和零状态响应。其中零输入响应与自然响应、零状态响应与受迫响应从表面上看相似，但是它们并不相同;1零输入响应是自然响应的一部分，但是自然响应还包括了零状态响应响应的一部分；2受迫响应是零状态响应的一部分，但零状态响应还包括自然响应的一部分；3总之，信号的响应包括三部分：............既是零输入，又是自然响应既是零状态，又有自然响应既是零状态，又有受迫响应。系统响应又有另外一种分法：瞬态响应：随时间增长而趋于零的部分；稳态响应：随时间增长而不趋向零的部分。对于稳定系统，自然响应必定属于瞬态响应，受迫响应则可能为瞬态响应，也可是稳态响应，具体情况视激励信号的形式而定。如果激励信号的指s系统的某个特征根相同，s

j

，则卷积项为：

j

j

j例1：电路如图所示，已