2.2 拉普拉斯变换

拉普拉斯变换(简称为拉氏变换)是分析研究线性动态系统的常用数学工具。通过拉普拉斯变换可将时域的微分方程变换为复数域的代数方程,这不仅使运算方便,也使系统的分析大为简化。

2.2.1 拉普拉斯变换的定义

函数ft)的拉普拉斯变换定义为

式中,ft)为时间t的函数,且当t<0时,ft)=0;s为复变数;L为运算符号;Fs)为ft)的拉普拉斯变换。

例2-5】 已知阶跃函数ft),计算阶跃函数ft)的拉普拉斯变换Fs)。

式中,A是常数。

:依据拉普拉斯变换定义有

A=1时称为单位阶跃函数,即ft)=1(t),其拉普拉斯变换为

一般地,式(2-9)所示的阶跃函数常常写为ft)=A·1(t)。

例2-6】 已知正弦函数

式中,Aω是常数。试计算正弦函数ft)的拉普拉斯变换Fs)。

:依据拉普拉斯变换的定义有

由此可知,函数的拉普拉斯变换通常可利用拉普拉斯变换的定义求得,但在实际应用中,该计算过程较为烦琐。因此提供了另外一种计算方法,即通过查拉普拉斯变换对照表来得到函数的拉普拉斯变换。常用函数及对应的拉普拉斯变换见表2-1,利用此表可查出已知时间函数的拉普拉斯变换,或者已知拉普拉斯变换所对应的时间函数。

表2-1 拉普拉斯变换对照表

(续)

2.2.2 拉普拉斯变换定理

下面介绍几个在线性控制系统分析和计算中常用的函数和拉普拉斯变换定理。

1. 位移定理

有时间函数ft),其拉普拉斯变换为Fs)=L[ft)],则有

时间函数ft)经过a的平移后,所得时间函数ft-a)的拉普拉斯变换,相当于时间函数ft)的拉普拉斯变换Fs)与e-as相乘。

在时间函数eatft)的拉普拉斯变换中,用(s-a)去替换Fs)中的sa是实数或复数。式(2-11)非常适合求解e-atsinωt和e-atcosωt等类似函数的拉普拉斯变换。

例2-7】 已知脉动函数为

式中,A是常量。试计算该脉动函数的拉普拉斯变换。

:该脉动函数ft)是一个从t=0开始、幅度为A的阶跃函数,再叠加一个从t=t0开始、幅度为A的负的阶跃函数,即

ft)=A·1(t-A·1(t-t0

ft)的拉普拉斯变换为

例2-8】 已知函数e-atcosωt,计算其拉普拉斯变换。

:已知cosωt的拉普拉斯变换为

则e-atcosωt的拉普拉斯变换为

例2-9】 已知某信号ft)如图2-6所示,试求其拉普拉斯变换。

图2-6 信号ft)的形状示意图

:由图2-6可知,ft)由若干个阶跃信号叠加而成,其表达式为

ft)=a·1(t)+(b-a)·1(t-t1)-(b-c)·1(t-t2-c·1(t-t3

已知L[1(t)]=,由式(2-10)有

2. 线性性质

F1s)=L[f1t)],F2s)=L[f2t)],ab为常数,则有

3. 相似定理

Fs)=L[ft)],则有

式中,a为实常数。

例2-10】 计算函数ft)=2e-5t的拉普拉斯变换。

:该题可利用拉普拉斯变换的定义进行求解,也可利用ft)与e-at相乘来求解。设ft)=e-t,则ft)的拉普拉斯变换为

需要注意的是,拉普拉斯变换的积分下限在某些情况下是有所区别的。若函数ft)在t=0处包含脉冲函数,则其拉普拉斯变换的积分下限必须明确地指出是0+或者0-。这两种下限的函数ft)的拉普拉斯变换是不相同的。

4. 初值定理

若函数ft)及其一阶导数都可进行拉普拉斯变换,而且存在,则函数ft)的初值为

即原函数ft)在t→0+时的极限值,取决于其拉普拉斯Fs)在s→∞时的极限值。初值定理不需要限制sFs)的极点位置,故其适用于正弦函数sinωt和余弦函数cosωt

5. 终值定理

若函数ft)及其一阶导数都可进行拉普拉斯变换,存在,且除在原点处有唯一的极点外,sFs)在包含jω轴的右半s平面内是解析的,则函数ft)的终值为

即原函数ft)在t→∞时的极限值,取决于其拉普拉斯变换Fs)在s→0时的极限值。根据式(2-15)即可求得ft)在t→∞时的极限值。当ft)是正弦函数sinωt或余弦函数cosωt时,sFs)在s=±jω处有极点,并且不存在,故终值定理不适用于正弦函数sinωt和余弦函数cosωt;当t→∞时,ft)→∞,则不存在,那么终值定理也不适用于此类情况。

6. 微分定理

若函数ft)及其一阶导数都可进行拉普拉斯变换,那么的拉普拉斯变换为

式中,f(0)是ft)在t=0处的初始值。

若函数ft)在t=0处不连续,f(0+)≠f(0-),则式(2-16)修正为

若函数ft)存在n阶导数,则ft)的n阶导数的拉普拉斯变换为

式中,f(0),f1(0),…,fn-1f(0)分别表示ft),,…,t=0处的值。

ft)及其各阶导数的初始值均为零,则ft)的n阶导数的拉普拉斯变换简化为snFs)。

例2-11】 已知有函数ft)=Acosωtt≥0,A是常量),计算其拉普拉斯变换。

:由查表法即可得知正弦函数sinωt的拉普拉斯变换。此例采用微分定理来推导。

查表2-1可知

依据微分定理,ft)的拉普拉斯变换计算如下:

7. 积分定理

若函数ft)及其一阶积分∫ft)dt都可进行拉普拉斯变换,则函数∫ft)dt的拉普拉斯变换为

式中,f-1(0)表示∫ft)dtt=0处的值。若f-1(0)=0,则∫ft)dt的拉普拉斯变换简化为Fs/s

若函数ft)在t=0处有脉冲函数,且f-1(0+)≠f-1(0-),则式(2-20)修正为

表2-2列出了拉普拉斯变换的常用定理和关系式。

表2-2 拉普拉斯变换常用定理和关系式

2.2.3 拉普拉斯反变换

由象函数Fs)求原函数ft)的过程称为拉普拉斯反变换,是用复变数表达式来推导其相应的时间函数表达式的数学运算。L-1代表拉普拉斯反变换,其计算公式为

式中,c是实常数,其值均大于Fs)的所有奇点的实部。积分的路线平行于jω轴,并且与jω轴的距离是c,以及该积分路线位于Fs)所有奇点的右侧。

从式(2-23)可知,利用拉普拉斯反变换的定义直接求解原函数是非常烦琐的过程。任何连续时间函数(即原函数)与其拉普拉斯变换(即象函数)是一一对应的,因此,在工程实践中常常采用查表法。对简单的象函数而言,可直接查表2-1求得其原函数;对复杂的象函数而言,通常先用部分分式展开法将复杂函数展开成简单函数的和,再查表2-1求其原函数。

象函数Fs)常常表示为两个s多项式的比,其形式如下:

式中,多项式系数a0a1,…,anb0b1,…,bm均为实常数;多项式Bs)的阶次不高于多项式As)的阶次即mnmn均为正整数。

为了将Fs)写成部分分式的形式,需要对Fs)的分母因式进行分解,则有

式中,s=-p1s=-p2,…,s=-pn称为Fs)的极点。下面分两种情况来介绍求拉普拉斯反变换的部分分式展开法。

1.As=0无重根情况

此时,Fs)的所有极点彼此不相等,即p1p2≠…≠pnFs)可展开为n个简单的部分分式之和,每个部分分式都以As)的一个因式作为分母,即

式中,ck为待定常数,k=1,2,…,n,称为Fs)在极点-pk处的留数,可计算如下:

通过查表2-1、表2-2以及拉普拉斯变换的线性性质,可求得原函数

例2-12】 求Fs)=的原函数ft)。

:将Fs)分解为

依据式(2-27),有

因此,可求得原函数

例2-13】 求Fs)=的原函数ft)。

:由于Fs)分子的阶次大于分母的阶次,对其整理有

式(2-29)中右边最后一项的拉普拉斯反变换见例2-12,因此只需求式中的第一项和第二项的拉普拉斯反变换。查表2-1、表2-2可知,L[δt)]=1,,故Fs)的原函数ft)为

Fs)有共轭复数极点p1p2时,其展开式为

式中,留数ckk=3,4,…,n)计算由式(2-27)确定。c1c2计算如下:

p 1是共轭复数,故式(2-31)的两边也是复数。若极点p1p2是共轭复数,则留数c1c2也是共轭复数。

2.As=0有重根情况

As)=0有r个重根p1,则Fs)可写为

式中,p1Fs)的重极点,pr+1,…,pnFs)的(n-r)个互不相等的极点;brbr-1,…,b1cr+1,…,cn为待定留数,其中,cr+1,…,cn按式(2-27)计算,但brbr-1,…,b1计算如下:

因此,原函数ft)为

例2-14】 求Fs)=的原函数ft)。

:由题可知,分母As)=0有4个根,即二重根s1=s2=-1,以及单根s3=0,s4=-3。将Fs)展为部分分式,则有

按式(2-32)计算得

按式(2-27)计算得

则原函数ft)为

例2-15】 求Fs)=的原函数ft)。

:由题可知,Fs)的分母多项式包含有一对共轭复数根。将Fs)先分解并整理如下:

若将一对共轭复数根直接采用留数法分解后再求解比较烦琐,因此本例采用查表法。经查表2-1和表2-2可得

综上所述,部分分式展开法的特点是将象函数Fs)展开成部分分式,使得Fs)的每一项都是简单函数,那么经查询拉普拉斯变换表或者常用函数的拉普拉斯变换,即可快速求得Fs)的拉普拉斯反变换ft)。