拉普拉斯变换及逆变换
第十二章 拉普拉斯变换及逆变换
拉普拉斯(Laplace)变换是分析和求解常系数线性微分方程的一种简便的方法,而且在自动控制系统的分析和综合中也起着重要的作用。我们经常应用拉普拉斯变换进行电路的复频域分析。本章将扼要地介绍拉普拉斯变换(以下简称拉氏变换)的基本概念、主要性质、逆变换以及它在解常系数线性微分方程中的应用。
第一节 拉普拉斯变换
在代数中,直接计算
N=6.28⨯
5781
是很复杂的,而引用对数后,可先把上式变换为
⨯20⨯(1.164)
2
35
然后通过查常用对数表和反对数表,就可算得原来要求的数N。
这是一种把复杂运算转化为简单运算的做法,而拉氏变换则是另一种化繁为简的做法。
13
lgN=lg6.28+(lg5781-lg9.8+2lg20)+lg1.164
35
一、拉氏变换的基本概念
定义12.1 设函数f(t)当t≥0时有定义,若广义积分
⎰
+∞
f(t)e-ptdt在P的某一区域内
收敛,则此积分就确定了一个参量为P的函数,记作F(P),即
0 (12.1)
称(12.1)式为函数f(t)的拉氏变换式,用记号L[f(t)]=F(P)表示。函数F(P)称为f(t)
F(P)=
⎰
+∞
f(t)e-ptdt
的拉氏变换(Laplace) (或称为f(t)的象函数)。函数f(t)称为F(P)的拉氏逆变换(或称为,记作 F(P)象原函数)
L-1[F(P)]=f(t),即f(t)=L-1[F(P)]。
关于拉氏变换的定义,在这里做两点说明:
(1)在定义中,只要求f(t)在t≥0时有定义。为了研究拉氏变换性质的方便,以后总假定在t
(2)在较为深入的讨论中,拉氏变换式中的参数P是在复数范围内取值。为了方便起见,本章我们把P作为实数来讨论,这并不影响对拉氏变换性质的研究和应用。
(3)拉氏变换是将给定的函数通过广义积分转换成一个新的函数,它是一种积分变换。一般来说,在科学技术中遇到的函数,它的拉氏变换总是存在的。
例12.1 求斜坡函数f(t)=at (t≥0,a为常数)的拉氏变换。
解:L[at]=
⎰
+∞
ate-ptdt=-
a=0+
p
⎰
+∞0
a+∞a-pt+∞a+∞-pt-pt
td(e)=[-e]0+⎰edt ⎰0ppp0aa+∞
e-ptdt=[-2e-pt]0=2(p>0)
pp
1
二、单位脉冲函数及其拉氏变换
在研究线性电路在脉冲电动势作用后所产生的电流时,要涉及到我们要介绍的脉冲函数,在原来电流为零的电路中,某一瞬时(设为t=0)进入一单位电量的脉冲,现要确定电路上的电流i(t),以Q(t)表示上述电路中的电量,则
由于电流强度是电量对时间的变化率,即
⎧0,t≠0,Q(t)=⎨
⎩1,t=0.
dQ(t)Q(t+∆t)-Q(t)
=lim∆t→0dt∆t,
所以,当t≠0时,i(t)=0;当t=0时,
i(t)=i(0)=lim
Q(0+∆t)-Q(0)1
=lim(-)=∞
∆t→0∆t→0∆t∆t。
上式说明,在通常意义下的函数类中找不到一个函数能够用来表示上述电路的电流强
度.为此,引进一个新的函数,这个函数称为狄拉克函数。
定义12.2
t
⎪1⎪
设δε(t)=⎨,0≤t≤ε,当ε→0时,δε(t)的极限δ(t)=limδε(t)
ε→0
⎪ε
t>ε⎪⎩0,
称为狄拉克(Dirac)函数,简称为δ-函数。
⎧0,δ(t)=当t≠0时,δ(t)的值为0;当t=0时,δ(t)的值为无穷大,即⎨
⎩∞,
+∞ε1+∞
显然,对任何ε>0,有⎰δε(t)dt=⎰=1,所以⎰δ(t)dt=1。
-∞
t≠0t=0
。
ε
-∞
工程技术中,常将δ-函数称为单位脉冲函数,有些工程书上,将δ-函数用一个长度
等于1的有向线段来表示,这个线段的长度表示δ-函数的积分,叫做δ-函数的强度。
例12.2 求单位脉冲信号δ(t)的拉氏变换。 解:根据拉氏变换的定义,有
L[δ(t)]=
⎰
+∞0
δ(t)e
1
-pt
dt=
⎰
ε
(lim
1
ε→0
ε
)e
-pt
dt+lim
ε→0
⎰
+∞
ε
0⋅e
-pt
dt=lim
ε→0
⎰
ε
1
ε
e-ptdt
e-ptε11-e-pε1(1-e-pε)'1pe-pε
=lim[-]0=lim=lim=lim=1ε→0εppε→0εpε→0(ε)'pε→01,
即
L[δ(t)]=1。
⎧0,
例12.3 现有一单位阶跃输入u(t)=⎨
⎩1,
解:L[u(t)]=
t
,求其拉氏变换。
⎰
+∞
u(t)e
-pt
dt=⎰1 e-ptdt=[-
0at
+∞
1-pt+∞1
e]0=,(p>0)。 pp
例12.4 求指数函数f(t)=e(a为常数)的拉氏变换。 解:L[e]=
at
⎰
+∞
e e
at-pt
dt=⎰e-(p-a)tdt=
+∞
1
,(p>a),即 p-a
2
L[eat]=
类似可得L[sinωt]=
ω
p2+ω
1
(p>a)p-a
(p>0);L[cosωt]=2
p
(p>0)。 22
p+ω
三、拉氏变换的性质
拉氏变换有以下几个主要性质,利用这些性质,可以求一些较为复杂的函数的拉氏变换。
性质12.1 (线性性质) 若a1,a2是常数,且L[f1(t)]=F1(p),L[f2(t)]=F2(p),则
L[a1f1(t)+a2f2(t)]=a1L[f1(t)]+a2L[f2(t)]=a1F1(P)+a2F2(p) (12.2)
证明:
L[a1f1(t)+a2f2(t)]=
⎰
+∞0
[a1f1(t)+a2f2(t)]e
-pt
dt=a1
⎰
+∞0
f1(t)e
-pt
dt+a2
⎰
+∞0
f2(t)e-ptdt
=a1L[f1(t)]+a2L[f2(t)]=a1F1(p)+a2F2(p)
例12.5 求函数f(t)=解:
1
(1-e-at)的拉氏变换 a
1111111L[(1-e-at)]=L[1-e-at]={L[1]-L[e-at]}=-=aaaapp+ap(p+a) 性质12.2(平移性质) 若L[f(t)]=F[p],则
L[ef(t)]=F(p-a)(a为常数) (12.3) 证明:
位移性质表明:象原函数乘以e等于其象函数左右平移|a|个单位。
atat
L[ef(t)]=
at
⎰
+∞
ef(t)e
-at
at-pt
dt=
⎰
+∞
f(t)e-(p-a)tdt=F(p-a)
例12.6 求L[te],L[e解 因为L[t]=
at
sinωt]和L[e-atcosωt]。
1ωp,,,由位移性质即得 L[cosωt]=L[sinωt]=22222
p+ωpp+ω
1ω-at
L[teat]=,L[esinωt]=,
(p-a)2(p+a)2+ω2
p+a
L[e-atcosωt]=。22
(p+a)+ω
性质12.3(滞后性质) 若L[f(t)]=F[p],则
L[f(t-a)]=e-apF(p) (a>0) (12.4)
证明:
a=0,
在拉氏变换的定义说明中已指出,当t
L[f(t-a)]=
⎰
+∞
f(t-a)e
-pt
dt
⎰
a
f(t-a)e
-pt
dt+
⎰
+∞
f(t-a)e-ptdt
当t-aL[f(t-a)]=
⎰
+∞
f(τ)e
-p(τ+a)
dτ=e
3
-ap
⎰
+∞
f(τ)e-pτdτ=e-apF(p)
滞后性质指出:象函数乘以e-ap等于其象原函数的图形沿t轴向右平移a个单位。 由于函数f(t-a)是当t≥a时才有非零数值。故与f(t)相比,在时间上滞后了一个a值,正是这个道理,我们才称它为滞后性质.在实际应用中,为了突出“滞后”这一特点,常在f(t-a)这个函数上再乘u(t-a),所以滞后性质也表示为
例12.7 求L[u(t-a)]。 解:因为L[u(t)]=
L[u(t-a)f(t-a)]=e-apF(p)
1-ap1,由滞后性质得L[u(t-a)]=e。 ppa(t-τ)
例12.8 求L[eu(t-τ)]。
11ata(t-τ)
解:因为L[e]=,所以L[e,(p>a) u(t-τ)]=e-τp
p-ap-a⎧0,t2c,a≤t
解: f(t)可用单位阶梯函数表示为f(t)=cu(t)+cu(t-a)-2cu(t-3a),于是 L[f(t)]=L[cu(t)+cu(t-a)-2cu(t-3a)
cccc
=+e-ap-2e-3ap=(1+eap-2e-3ap)pppp ,
由拉氏变换定义来验证:
L[f(t)]=
=
⎰
a0
ce
-pt
dt+
⎰
3aa
2ce-ptdt
cc
(1-e-ap+2e-ap-2e-3ap)=(1+e-ap-2e-3ap)pp。
'
性质12.4(微分性质) 若L[f(t)]=F[p],并设f(t)在[0,+∞)上连续,f(t)为分
段连续,则
L[f'(t)]=pF(p)-f(0) (12.5)
证明:由拉氏变换定义及分部积分法,得
L[f'(t)]=
⎰
+∞0
f'(t)e
-pt
dt=[f(t)e
t→+∞
-pt+∞
]+P
⎰
+∞0
f(t)e-ptdt
,
可以证明,在L[f(t)]存在的条件下,必有limf(t)e
-pt
=0。因此,
微分性质表明:一个函数求导后取拉氏变换等于这个函数的拉氏变换乘以参数p,再减去函数的初始值。
应用上述结果,对二阶导数可以推得 同理,可得
L[f'(t)]=0-f(0)+pL[f(t)]=pF(p)-f(0)
L[f''(t)]=pL[f'(t)]-f'(0)=p{pF(p)-f(0)}-f'(0)=p2F(p)-{pf(0)+f'(0)}
L[f'''(t)]=p3F(p)-{p2f(0)+pf'(0)+f''(0)}
以此类推,可得
L[f(n)(t)]=pnF(p)-{pn-1f(0)+pn-2f'(0)+ f(n-1)(0)} (12.6)
由此可见,f(t)各阶导数的拉氏变换可以由p的乘方与象函数F[p]的代数式表示出
4
来.特别是当初值f(0)=f'(0)=f''(0)=f
(n-1)
(0)=0时,有更简单的结果
L[f(n)(t)]=pnF(p),(n=1,2, ) (12.7)
利用这个性质,可将f(t)的微分方程转化为F(p)的代数方程。 例12.10 利用微分性质求L[sinωt]和L[cosωt]。
解:令f(t)=sinωt,则f(t)=sinωtf(0)=(12.6)式,得 即
0,f'(=0)ωf,
"(=0)-ω2,ωsint由
L[-ω2sinωt]=L[f''(t)]=p2L[f(t)]-pf(0)-f'(0),
-ω2L[sinωt]=p2L[sinωt]-ω,
移项化简得
L[sinωt]=
利用上述结果,cosωt=
ω
p+ω
2
2
1
ω
(sinωt)'及(12.5)式,可得
1
L[cosωt]=L[
1
ω
(sinωt)']={p⋅
ω
L[(sinωt)']=-0}=
1
ω
{pL[sinωt]-sin0}
p
ωp2+ω2p2+ω2.
性质12.5(积分性质) 若L[f(t)]=F(p) (p≠0),且设f(t)连续,则
=
1ω
L[
证明:令ϕ(t)=
⎰
t0
f(x)dx]=
F(p)
p (12.8)
⎰
t
f(x)dt,显见ϕ(0)=0,且因ϕ'(t)=f(t),由微分性质,得
t
t
L[ϕ'(t)]=pL[ϕ(t)]-ϕ(0),而L[ϕ'(t)]=L[f(t)]=F(p),所以有
F(p)=pL[ϕ(t)]=pL[⎰f(x)dx],即L[⎰f(x)dx]=
1
F(p)。 p
积分性质表明:一个函数积分后再取拉氏变换,等于这个函数的象函数除以参数p。 例12.11 求L[t](n是正整数)。 解:因为
n
t=
⎰
t0
1dx,t=
2
⎰
t0
2xdx,t=
3
⎰
t0
3xdx
2
,…,
t=
n
⎰
t0
nxn-1dx
,
所以由(12.8)式即得
L[1]p1!
L[t]=L[1dx]===2,
0pppt
2L[t]2!2
L[t]=L[2xdx]==3,
0pp
⎰
t
⎰
L[t]=L[3
一般地,有
n
3
⎰
t0
3L[t2]3!
xdx]==4,
pp
2
L[t]=L[n
⎰
t0
x
n-1
nL[tn-1]n!dt]==n+1
pp
5
性质12.6 若L[f(t)]=F[p],则a>0时
L[f(at)]=
性质12.7 若L[f(t)]=F[p],则
1pF()
aa (12.9)
L[tnf(t)]=(-1)nF(n)(p) (12.10)
性质12.8 若L[f(t)]=F[p],且lim
t→0
f(t)
存在,则 t
+∞p
f(t)L[]=
t
例12.12 求L[tsinωt]。 解:因为L[sinωt]=
⎰
F(p)dp
(12.11)
ω
p+ω
2
2
,由(12.10)式可得
L[tsinωt]=(-1)
例12.13 求L[
sint
]。 t
1sint
解:因为L[sint]=2,而且lim=1,所以由(12.11)式可得
t→0p+1t
+∞
sint1π+∞L[]=dp=arctgp|=-arctgpp
pt2p2+1
+∞sintπ即⎰-ptdt=-arctgp。因此,当p=0时,得到一个广义积分的值
0t2
+∞sintπ
dt=
0t2
dω2pω
(2)=dpp+ω2(p2+ω2)2
⎰
⎰
这个结果用原来的广义积分的计算方法是得不到的。
现将拉氏变换的八个性质和在实际应用中常用的一些函数的象函数分别列表如下:
6
7
习题12.1
1.求下列函数的拉氏变换 (1)f(t)=e
-4t
(2)f(t)=t (3)f(t)=te
(4)f(t)=sin(ωt+ϕ)(ω,ϕ是常数) 2.求下列题中函数的拉氏变换 (1)3e
-4t
at
2
(2)5sin2t-3cost
(3)f(t)=⎨
⎧-1,⎩1,
0≤t≤4t≥4
(4)f(t)=⎨
⎧sint,⎩t,
0≤t≤πt≥π
8
0≤t≤2⎧0,⎪
2≤t
第二节 拉普拉斯逆变换
前面我们主要讨论了怎样由已知函数f(t)求它的象函数F(p)的问题.运算法的另一面是已知象函数F(p)要求它的象原函数f(t),这就是拉斯逆变换问题.在控制工程中,求拉
氏反变换的简便方法是利用拉氏变换表。同时把常用的拉氏变换的性质用逆变换形式一一列出.
性质12.9(先行性质)
L-1[a1F1(p)+a2F2(p)]=a1L-1[F1(p)]+a2L-1[F2(p)]=a1f1(t)+a2f2(t)。
性质12.10(平移性质) L[F(p-a)]=eL[F(p)]=ef(t)。 性质12.11(滞后性质) L[e例12.14 求F(p)=解:
-1
-ap-1
at
-1
at
F(p)]=f(t-a)⋅u(t-a)。
2p+3
的逆变换。 2
p-2p+5
2p+3-12(p-1)+5]=L[]22p-2p+5(p-1)+4
p-15-12
=2L-1[]+L[]22
2(p-1)+4(p-1)+4
p52
=2etL-1[2]+etL-1[2]
p+42p+4
55
=2etcos2t+etsin2t=et[2cos2t+sin2t]
22
f(t)=L-1[
在运用拉氏变换解决工程技术中的应有问题时,通常遇到的象函数常常是有理分式,对
于有理分式一般可采用部分分式方法将它分解为较为简单的分式之和,然后再利用拉氏变换表求出象原函数。
例12.15 求F(p)=
p+3
的逆变换。 32
p+4p+4p
解:先将F[p]分解为几个简单分式之和:
p+3p+3ABC
==++
pp+2(p+2)2, p3+4p2+4PP(p+2)2
313
用待定系数法求得A=,B=-,C=-,所以
424
331
p+3--F(p)=3=
p+4p2+4Ppp+2(p+2)2,
于是
f(t)=L-1[F(p)]=L-1[
313111
--]24p4p+22(p+2)
9
3-113111L[]-L-1[]-L-1[]2
p4p+22(p+2) 4331=-e-2t-te-2t
2 44
=
习题13.2
求下列题中函数的拉氏逆变换
24p
2. F(p)=2
p+16p-3
2p-81
3. F(p)=2 4. F(p)=
p+36p(p+1)(p+2)
p2p2+1
5.F(p)=3 6. F(p)= 22
p+6p+9pp(p-1)
1.
F(p)=
第三节 拉氏变换在电学中的应用
一、求解常微分方程
例12.16 求微分方程x'(t)+2x(t)=0满足初值条件x(0)=3的解。 解:第一步 对方程两边取拉氏变换,并设L[x(t)]=X(p):
L[x'(t)+2x(t)]=L[0], L[x'(t)]+2L[x(t)]=0,
pX(p)-x(0)+2X(p)=0。
将初始条件x(0)=3代入上式,得
这样,原来的微分方程经过拉氏变换后,就得到了一个象函数的代数方程。
第二步 解出X(p):X(p)=
(p+2)X(p)=3
3
p+2
-1
-1
第三步 求象函数的拉氏逆变换:x(t)=L[X(p)]=L[这样就得到了微分方程的解x(t)=3e
-2t
3
]=3e-2t p+2
-t
。
例12.17 有一个二阶动态电路满足微分方程y''-3y'+2y=2e,并且其初值条件
y(0)=2,y'(0)=-1,求其解。
解:对所给微分方程的两边分别作拉氏变换.设L[y(t)]=Y(p)=Y,则得
[p2Y-py(0)-y'(0)]-3[pY-y(0)]+2Y=
将初值条件y(0)=2,y'(0)=-1,代入,得到Y的代数方程
2p+1
(p2-3p+2)Y=
即
10
2
+2p-7p+1
2p2-5p-5
(p-3p+2)Y=
p+1
2
解出Y,得
将上式分解为部分分式
2p2-5p-5
Y=
(p+1)(p-2)(p-1)
再取拉氏逆变换,就得到满足所给初值条件的方程的特解为
17
4
Y=3+-3
p+1p-1p-2
17y(t)=e-t+4et-e2t
33
用拉氏变换还可以解常系数线性微分方程组。
二、电学应用举例
例12.18 求图示电路的输入运算阻抗Zin(s)
解:由串并联关系得
1⎫2s2+12s2+1⎛
=Zin(s) =2 s+⎪=22
s⎭s+2s+1s+1⎝
例12.19 求图(a)所示电路中的i(t)、uC(t)。
()()
(a) (b)
解:先画出运算电路如图(b)所示。由运算电路得
1
K3KK2s
I(s)===1++
10(s+1)s2+6s+10s+1s+3-js+3+j6+s+
s
其中
K1=I(s)(s+1)s=-1=
s
s2+6s+10
s=-1
1=-
5
K2=I(s)(s+3-js=-3+j=
s
s+1s+3+js=-3+j
=0.1-j0.7=
12
-81.87︒
12K3=I(s)(s+3+j)s=-3-j==0.1+j0.7=
则
⎡1⎤
i(t)=L-1[I(s)]=⎢-e-t+2e-3 tcos(t-81.87︒)⎥ ε(t) A
⎣5⎦UC(s)=
其中
K3K1K21010
I(s)= =++ss+1s+3-js+3+js+1s2+6s+10
K1 =UC(s)(s+1)s=-1=
10
s2+6s+10
=2
s=-1
K2 =UC(s)(s+3-j)s=-3+j=
10
s+1s+3+js=-3+j
=-1+j2=5
K3 =UC(s)(s+3+j)s=-3-j= -116.565︒
则
uC(t)=L-1[UC(s)]=[2e-t+25e-3 tcos(t+116.565︒)] ε(t)V
习题 12.3
1.求解一输入响应电路的微分方程。
2. 求图(a)所示电路中的回路电流i1和i2.。
di
+5i=10e-3t,i(0)=0dt
自测题
1. 求各函数的拉氏变换
(1)
0≤t≤1⎧0,
⎪
f(t)⎨-1,1≤t
⎪2,2≤t⎩
(3)
f(t)=8sin23t (4)f(t)=1+tet
13p+9 (2)F(p)=
22
p(p-1)p+2p+10
-p
2. 求各象函数的逆变换 (1)F(p)=
22e
(3)F(p)=5p-15p+7 (4)F(p)=
(p+1)(p-2)3-e-2pp
3.如图所示电路激励为i(t)=δ(t),响应为u1、u2。求阶跃响应S1(t)、S2(t)。
2-
4.求图示两电路的输入运算阻抗Zin(s)