系统的数学模型—微分方程与传输算子
扫描二维码
随时随地手机看文章
不涉及任何数学变换,而直接在时间变量域内对系统进行分析,称为系统的时域分析。其方法有两种:时域经典法与时域卷积法。 时域经典法就是直接求解系统微分方程的方法。这种方法的优点是直观,物理概念清楚,缺点是求解过程冗繁,应用上也有局限性。所以在20世纪50年代以前,人们普遍喜欢采用变换域分析方法(例如拉普拉斯变换法),而较少采用时域经典法。20世纪50年代以后,由于δ(t)函数及计算机的普遍应用,时域卷积法得到了迅速发展,且不断成熟和完善,已成为系统分析的重要方法之一。时域分析法是各种变换域分析法的基础。 在本章中,首先建立系统的数学模型——微分方程,然后用经典法求系统的零输入响应,用时域卷积法求系统的零状态响应,再把零输入响应与零状态响应相加,即得系统的全响应。其思路与程序是: |
其次,将介绍:系统相当于一个微分方程;系统相当于一个传输算子H(p);系统相当于一个信号——冲激响应h(t)。对系统进行分析,就是研究激励信号f(t)与冲激响应信号h(t)之间的关系,这种关系就是卷积积分。 |
2-1 系统的数学模型——微分方程与传输算子 |
研究系统,首先要建立系统的数学模型——微分方程。建立电路系统微分方程的依据是电路的两种约束:拓扑约束(KCL,KVL)与元件约束(元件的时域伏安关系)。为了使读者容易理解和接受,我们采取从特殊到一般的方法来研究。 图2-1(a)所示为一含有三个独立动态元件的双网孔电路,其中 为激励, , 为响应。对两个网孔回路可列出KVL方程为 |
|
|||
上两式为含有两个待求变量 , 的联立微分积分方程。 为了得到只含有一个变量的微分方程, 须引用微分算子 ,即 |
|||
|
, ,…, |
|
|
在引入了微分算子 后,上述微分方程即可写 |
|||
|
|
||
即 |
|||
|
(2-1) |
||
根据式(2-1)可画出算子形式的电路模型,如图2-1(b)所示。将图2-1(a)与(b)对照, 可很容易地根据图2-1(a)画出图2-1(b),即将L改写成Lp,将C改写成 , 其余一切均不变。当画出了算子电路模型后,即可很容易地根据图2-1(b)算子电路模型列写出式(2-1)。 |
|||
给式(2-1)等号两端同时左乘以p,即得联立的微分方程,即 |
|||
|
|||
将已知数据代入上式,得 |
|||
|
(2-2) |
||
用行列式法从式(2-2)中可求得响应i1(t)为 |
|||
|
|||
注意,在上式的演算过程中,消去了分子与分母中的公因子p。这是因为所研究的电路是三阶的, 因而电路的微分方程也应是三阶的。但应注意,并不是在任何情况下分子与分母中的公因子都可消去。 有的情况可以消去,有的情况则不能消去,视具体情况而定。故有 |
|||
|
|||
即 |
|||
|
|||
即 |
|||
|
|||
上式即为待求变量为i1(t)的三阶常系数线性非齐次常微分方程。 方程等号左端为响应i1(t)及其各阶导数的线性组合, 等号右端为激励f(t)及其各阶导数的线性组合。 |
|||
利用同样的方法可求得i2(t)为 |
|||
|
|||
即 |
|||
|
|||
即 |
|||
|
|||
即 |
|||
|
|||
上式即为描述响应i2(t)与激励f(t)关系的微分方程。 |
|||
推广之,对于n阶系统,若设y(t)为响应变量, f(t)为激励,如图2-2所示,则系统微分方程的一般形式为 |
|||
|
(2-3) |
||
|
|||
用微分算子 表示则为 |
|||
|
|||
或写成 |
|||
|
|||
又可写成 |
|||
|
|||
式中 |
|
|
|
称为系统或微分方程式(2-3)的特征多项式; |
|||
|
|||
|
(2-4) |
||
H(p)称为响应y(t)对激励f(t)的传输算子或转移算子,它为p的两个实系数有理多项式之比, 其分母即为微分方程的特征多项式D(p)。H(p)描述了系统本身的特性,与系统的激励和响应无关。 |
|||
这里指出一点:字母p在本质上是一个微分算子,但从数学形式的角度,以后可以人为地把它看成是 一个变量(一般是复数)。这样,传输算子H(p)就是p的两个实系数有理多项式之比。 |
|||
例2-1 图2-3(a)所示电路。求响应u1(t),u2(t)对激励 的传输算子及u1(t),u2(t)分别对i(t)的微分方程。 |
|||
解 其算子形式的电路如图2-3(b)所示。对节点①,②列算子形式的KCL方程为 |
|||
|
|||
代入数据得 |
|||
|
|||
对上式各项同时左乘以p,并整理得 |
|||
|
|||
用行列式法联解得 |
|||
|
|||
|
|||
故得u1(t)对i(t),u2(t)对i(t)的传输算子分别为 |
|||
|
|||
|
|||
进而得u1(t),u2(t)分别对i(t)的微分方程为 |
|||
|
|||
|
|||
即 |
|||
|
|||
|
|||
可见,对不同的响应u1(t),u2(t),其特征多项式 都是相同的, 这就是系统特征多项式的不变性与相同性。 |