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

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

于是

因此

这样
A=1
类似地,

于是

或者说
B=-1
于是我们有

通过拉普拉斯变换可得

因此

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

首先,公式
s2+2s+5=0
的根为

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

从拉普拉斯变换表可知

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

用已有的拉普拉斯变换表可得
f(t)=2e-tcos(2t)us(t)+5e-tsin(2t)us(t)
我们可以用MATLAB来验证这个结果:

MATLAB的返回值应为
exp(-t)(2cos(2t)+5sin(2t))
我们需要处理复数,所以现在简单回顾一下复数运算。
题外话 复数回顾
令
c=a+jb
表示一个复数,其中a和b是实数,。c的共轭复数用c*表示,定义为
c*≜a-jb
注意
(c*)*=(a-jb)*=a+jb=c
c的大小定义为

它也等于c*的大小,即

我们也可以用极坐标的形式表示复数。定义
∠c≜tan-1(b,a)
如图2-1所示,tan-1(b,a)在大多数计算机语言中表示为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∠c,c*=a-jb=|c|e-j∠c

于是我们得到

也就是
|c*|=|c|
∠c*=-∠c
下面我们证明:

类似地,

例12 (续)
我们用F(s)的复共轭极点对部分分式展开重做之前的例子。

那么

并且

因此

同样,

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

这里需要注意的一点是

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

结果可以用MATLAB进行检验:

MATLAB的返回值应为

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

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

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

MATLAB的返回值应为


把β1,写成极坐标形式:

在例11中我们得到
f(t)=2e-tcos(2t)us(t)+5e-tsin(2t)us(t)
所以结果应该为
2=2|β1|cos(∠β1)=2×2.6932cos(-1.19)
5=-2|β1|sin(∠β1)=-2×2.6932sin(-1.19)
此式成立(用MATLAB来检查)。
例13 多根问题
令

F(s)的部分分式展开为

首先乘s(s+2)2得到
1=A0(s+2)2+A1s(s+2)+A2s
或
1=(A0+A1)s2+(4A0+2A1+A2)s+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
解得

于是F(s)可以写成

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

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

或写成

解得

于是

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

于是

并且

同时

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

此外


图2-3 ∠β1=π/2+π/3
回顾之前我们的定义tan-1(b,a)与计算机语言命令atan2(b,a)相同。最终可得到

结果和式(2.18)中的一样。
2.2.1 非严格正则有理函数
在所有部分分式的例子中我们得到了严格正则有理函数。也就是说,我们有

假设考虑:

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

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

它的返回值应该为

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

这样

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

最后得到

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

它的返回值应该为

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