系统建模与控制导论
上QQ阅读APP看本书,新人免费读10天
设备和账号都新为新人

2.2 部分分式展开

在用拉普拉斯变换解微分方程时,最后一步通常是求有理函数s的拉普拉斯逆变换[2],在之前的例子中,我们用了这样一个公式:

即有理函数被分解为两个已知逆变换的更简单的函数。这个分解过程被称为部分分式展开。我们将通过一系列例子说明如何做到这一点。

例10

将其改写为

于是

因此

这样

A=1

类似地,

于是

或者说

B=-1

于是我们有

通过拉普拉斯变换可得

因此

例11

我们要求它的拉普拉斯逆变换

首先,公式

s2+2s+5=0

的根为

我们称Fs)的分母的根为Fs)的极点,可以把它写作

从拉普拉斯变换表可知

确定σ=-1,ω=2,那么Fs)可以重新写为

用已有的拉普拉斯变换表可得

ft)=2e-tcos(2tust)+5e-tsin(2tust

我们可以用MATLAB来验证这个结果:

MATLAB的返回值应为

exp(-t)(2cos(2t)+5sin(2t))

我们需要处理复数,所以现在简单回顾一下复数运算。

题外话 复数回顾

c=a+jb

表示一个复数,其中ab是实数,c的共轭复数用c*表示,定义为

c*a-jb

注意

c**=(a-jb*=a+jb=c

c的大小定义为

它也等于c*的大小,即

我们也可以用极坐标的形式表示复数。定义

c≜tan-1ba

如图2-1所示,tan-1ba)在大多数计算机语言中表示为atan2(b,a),这样便于使c=a+bj的角度在正确的象限内。例如,令c1=-1+j,那么

c1=tan-1(1,-1)=atan 2(1,-1)=3π/4=2.3562

相反,如果我们考虑c2=1-j,那么

c2=tan-1(-1,1)=atan2(-1,1)=-π/4=-0.7854

我们可以把c写成极坐标形式,即

图2-1 c=a+jb=|c|ej∠cc*=a-jb=|c|e-j∠c

于是我们得到

也就是

|c*|=|c|

c*=-∠c

下面我们证明

类似地,

例12 (续)

我们用Fs)的复共轭极点对部分分式展开重做之前的例子。

那么

并且

因此

同样,

因此β1=1-2.5 j。类似地,

这里需要注意的一点是

永远都是这样。回到部分分式展开,我们已经证明了

结果可以用MATLAB进行检验:

MATLAB的返回值应为

现在我们将β1=1-2.5 j转换为极坐标形式(见图2-2)。

我们有

再次用MATLAB检验我们的答案:

图2-2 将1-2.5j转换为极坐标形式

MATLAB的返回值应为

β1写成极坐标形式:

在例11中我们得到

ft)=2e-tcos(2tust)+5e-tsin(2tust

所以结果应该为

2=2|β1|cos(∠β1)=2×2.6932cos(-1.19)

5=-2|β1|sin(∠β1)=-2×2.6932sin(-1.19)

此式成立(用MATLAB来检查)。

例13 多根问题

Fs)的部分分式展开为

首先乘ss+2)2得到

1=A0s+2)2+A1ss+2)+A2s

1=(A0+A1s2+(4A0+2A1+A2s+4A0

s的系数相等,我们可以得到

A0=1/4,A1=-1/4,A2=-1/2

因此

最终可得到[3]

例14

对于例13的部分分式展开,更简单的方法如下:

通过乘(s+2)2得到

然后设s=-2,可以得到A2=-1/2。于是可以得到

上式可变形为

或者

就像例13的结果一样:A0=1/4,A1=-1/4。

例15

我们可以用二次方程来解,即

s2+s+1=0

解得

于是Fs)可以写成

然后我们可以用式(2.16)对复共轭极点进行展开。然而,还存在另外一种使用式(2.17)的方法。可以写成[4]

通过乘以ss2+s+1)来消除分式,得到

或写成

解得

于是

例16 (再次)

让我们用“硬核”的方式重做前面式(2.16)的例子。有

于是

并且

同时

为了给出和例15中相同形式的答案,我们必须将β1转换为极坐标形式(见图2-3)。且

此外

图2-3 ∠β1=π/2+π/3

回顾之前我们的定义tan-1ba)与计算机语言命令atan2(b,a)相同。最终可得到

结果和式(2.18)中的一样。

2.2.1 非严格正则有理函数

在所有部分分式的例子中我们得到了严格正则有理函数。也就是说,我们有

假设考虑:

其中,Fs)是正则有理函数,因为deg{bs)}≤deg{as)},但不是严格正则有理函数,因为deg{bs)}=deg{as)}。为了使部分分式展开法有效,我们必须先将分子除以分母,如下所示:

在MATLAB中我们可以编写程序:

它的返回值应该为

例17

另一个例子是令,它不是正则有理函数。我们首先用长除法将s2+5s+4除以s3,可得

这样

接下来,我们对做部分分式展开,得到

最后得到

在MATLAB中我们可以编写程序:

它的返回值应该为

注意 在我们把拉普拉斯变换应用到物理系统中时,Fs)往往是严格正则有理函数。