第十二章 拉普拉斯变换及逆变换
拉普拉斯(Laplace)变换是分析和求解常系数线性微分方程的一种简便的方法,而且在自动控制系统的分析和综合中也起着重要的作用。我们经常应用拉普拉斯变换进行电路的复频域分析。本章将扼要地介绍拉普拉斯变换(以下简称拉氏变换)的基本概念、主要性质、逆变换以及它在解常系数线性微分方程中的应用。
第一节 拉普拉斯变换
在代数中,直接计算
N6.28357819.820(1.164)
235是很复杂的,而引用对数后,可先把上式变换为
然后通过查常用对数表和反对数表,就可算得原来要求的数N。
这是一种把复杂运算转化为简单运算的做法,而拉氏变换则是另一种化繁为简的做法。
13lgNlg6.28(lg5781lg9.82lg20)lg1.16435
一、拉氏变换的基本概念
定义12.1 设函数f(t)当t0时有定义,若广义积分
0f(t)eptdt在P的某一区域内
收敛,则此积分就确定了一个参量为P的函数,记作F(P),即
(12.1)
称(12.1)式为函数f(t)的拉氏变换式,用记号L[f(t)]F(P)表示。函数F(P)称为f(t)0F(P)f(t)eptdt的拉氏变换(Laplace) (或称为f(t)的象函数)。函数f(t)称为F(P)的拉氏逆变换(或称为
F(P)象原函数),记作
L1[F(P)]f(t),即f(t)L1[F(P)]。
关于拉氏变换的定义,在这里做两点说明:
(1)在定义中,只要求f(t)在t0时有定义。为了研究拉氏变换性质的方便,以后总假定在t0时,f(t)0。
(2)在较为深入的讨论中,拉氏变换式中的参数P是在复数范围内取值。为了方便起见,本章我们把P作为实数来讨论,这并不影响对拉氏变换性质的研究和应用。
(3)拉氏变换是将给定的函数通过广义积分转换成一个新的函数,它是一种积分变换。一般来说,在科学技术中遇到的函数,它的拉氏变换总是存在的。
例12.1 求斜坡函数f(t)at (t0,a为常数)的拉氏变换。
解:L[at]0ateptdtaaptaptpttd(e)[e]0edt 0ppp0- - -专业资料-
- - -
a0p0eptdt[aptae]0p2p2(p0)
二、单位脉冲函数及其拉氏变换
在研究线性电路在脉冲电动势作用后所产生的电流时,要涉及到我们要介绍的脉冲函
数,在原来电流为零的电路中,某一瞬时(设为t0)进入一单位电量的脉冲,现要确定电路上的电流i(t),以Q(t)表示上述电路中的电量,则
0,t0,Q(t)1,t0.
由于电流强度是电量对时间的变化率,即
dQ(t)Q(tt)Q(t)limt0dtt,
所以,当t0时,i(t)0;当t0时,
Q(0t)Q(0)1i(0)limlim()t0t0tt。
i(t)上式说明,在通常意义下的函数类中找不到一个函数能够用来表示上述电路的电流强
度.为此,引进一个新的函数,这个函数称为狄拉克函数。
定义12.2
t00,1设(t),0t,当0时,(t)的极限(t)lim(t)
0t0,称为狄拉克(Dirac)函数,简称为δ-函数。
t00,当t0时,(t)的值为0;当t0时,(t)的值为无穷大,即(t)。
t0,1显然,对任何0,有(t)dtdt1,所以(t)dt1。
0工程技术中,常将函数称为单位脉冲函数,有些工程书上,将函数用一个长度等于1的有向线段来表示,这个线段的长度表示函数的积分,叫做函数的强度。
例12.2 求单位脉冲信号(t)的拉氏变换。
解:根据拉氏变换的定义,有
L[(t)]0(t)e1ptdt0(lim10)eptdtlim00eptdtlim001eptdtept11ep1(1ep)1peplim[]0limlimlim10000ppp()p1,
即
L[(t)]1。
例12.3现有一单位阶跃输入u(t)0,1,t0t0,求其拉氏变换。
- - -专业资料-
- - -
解:L[u(t)]00u(t)eptdt1eptdt[0at1pt1e]0,(p0)。 pp例12.4求指数函数f(t)e(a为常数)的拉氏变换。 解:L[e]ateeatptdte(pa)tdt01,(pa),即 paL[eat]类似可得L[sint]p21(pa)pa
(p0);L[cost]2p(p0)。 22p三、拉氏变换的性质
拉氏变换有以下几个主要性质,利用这些性质,可以求一些较为复杂的函数的拉氏变换。
性质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)]eptdta10f1(t)eptdta20f2(t)eptdta1L[f1(t)]a2L[f2(t)]a1F1(p)a2F2(p)
1at例12.5求函数f(t)(1e)的拉氏变换
a解:
1111111L[(1eat)]L[1eat]{L[1]L[eat]}{}aaaappap(pa) 性质12.2(平移性质) 若L[f(t)]F[p],则
L[eatf(t)]F(pa)(a为常数) (12.3)
证明:
位移性质表明:象原函数乘以e等于其象函数左右平移|a|个单位。
00atL[ef(t)]atef(t)eatatptdtf(t)e(pa)tdtF(pa)例12.6 求L[te],L[eatsint]和L[eatcost]。
解 因为L[t]1pL[sint]L[cost],,,由位移性质即得 22222ppp1atL[teat],L[esint],(pa)2(pa)22pa。22(pa)性质12.3(滞后性质) 若L[f(t)]F[p],则
L[eatcost]证明:
L[f(ta)]eapF(p)(a0)(12.4)
- - -专业资料-
- - -
0aa=0,
在拉氏变换的定义说明中已指出,当t0时,f(t)0。因此,对于函数f(ta),
L[f(ta)]f(ta)eptdtf(ta)eptdtf(ta)eptdt当ta0(即ta)时,f(ta)0,所以上式右端的第一个积分为0,对于第二个积分,令ta,则
滞后性质指出:象函数乘以e等于其象原函数的图形沿t轴向右平移a个单位。 由于函数f(ta)是当ta时才有非零数值。故与f(t)相比,在时间上滞后了一个a00apL[f(ta)]f()ep(a)deapf()epdeapF(p)值,正是这个道理,我们才称它为滞后性质.在实际应用中,为了突出“滞后”这一特点,常在f(ta)这个函数上再乘u(ta),所以滞后性质也表示为
例12.7 求L[u(ta)]。 解:因为L[u(t)]L[u(ta)f(ta)]eapF(p)
1ap1,由滞后性质得L[u(ta)]e。 ppa(t)u(t)]。 例12.8 求L[e11ata(t)u(t)]ep解:因为L[e],所以L[e,(pa)
papa0,t0c,0taf(t)例12.9已知,求L[f(t)]。 2c,at3a0,t3a解:f(t)可用单位阶梯函数表示为f(t)cu(t)cu(ta)2cu(t3a),于是 L[f(t)]L[cu(t)cu(ta)2cu(t3a) cccceap2e3ap(1eap2e3ap)pppp,
由拉氏变换定义来验证:
L[f(t)]a0ceptdt3aa2ceptdt
cc(1eap2eap2e3ap)(1eap2e3ap)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)]0f(t)eptdt[f(t)ept]0Ppt0f(t)eptdt,
可以证明,在L[f(t)]存在的条件下,必有limf(t)et0。因此,
微分性质表明:一个函数求导后取拉氏变换等于这个函数的拉氏变换乘以参数p,再
- - -专业资料-
L[f(t)]0f(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){pn1f(0)pn2f(0)f(n1)(0)} (12.6)
由此可见,f(t)各阶导数的拉氏变换可以由p的乘方与象函数F[p]的代数式表示出来.特别是当初值f(0)f'(0)f''(0)f(n1)(0)0时,有更简单的结果
L[f(n)(t)]pnF(p),(n1,2,) (12.7)
利用这个性质,可将f(t)的微分方程转化为F(p)的代数方程。 例12.10 利用微分性质求L[sint]和L[cost]。
解:令f(t)sint,则f(t)sintf(0)0,f'(0),f\"(0)sint,由
2(12.6)式,得
L[2sint]L[f(t)]p2L[f(t)]pf(0)f(0),
即
2L[sint]p2L[sint],
移项化简得
L[sint]利用上述结果,costp22
1(sint)'及(12.5)式,可得
111L[cost]L[(sint)]L[(sint)]{pL[sint]sin0}pp22p22.
性质12.5(积分性质) 若L[f(t)]F(p)(p0),且设f(t)连续,则
tF(p)L[f(x)dx]0p(12.8)
{p0}1
证明:令(t)t0f(x)dt,显见(0)0,且因'(t)f(t),由微分性质,得
tL[(t)]pL[(t)](0),而L[(t)]L[f(t)]F(p),所以有
F(p)pL[(t)]pL[f(x)dx],即L[f(x)dx]0t01F(p)。 p积分性质表明:一个函数积分后再取拉氏变换,等于这个函数的象函数除以参数p。 例12.11 求L[t](n是正整数)。
n解:因为
tt01dx,t2t02xdx,t3t03xdx2,…,
tnt0nxn1dx,
所以由(12.8)式即得
- - -专业资料-
- - -
1L[1]p1!L[t]L[1dx]2,0pppt2L[t]2!2L[t]L[2xdx]3,0pptL[t]L[33t03L[t2]3!xdx]4,pp
2一般地,有
nL[tn1]n!L[t]L[nxdt]n10pp
性质12.6 若L[f(t)]F[p],则a0时
1pL[f(at)]F()aa (12.9)
性质12.7 若L[f(t)]F[p],则
ntn1L[tnf(t)](1)nF(n)(p) (12.10)
性质12.8若L[f(t)]F[p],且limt0f(t)存在,则 tf(t)L[]t例12.12 求L[tsint]。 解:因为L[sint]pF(p)dp (12.11)
p22,由(12.10)式可得
L[tsint](1)例12.13 求L[sint]。 t1sint解:因为L[sint]2,而且lim1,所以由(12.11)式可得
t0p1tsint1L[]dparctgp|arctgpp2pt2p1
sint即eptdtarctgp。因此,当p0时,得到一个广义积分的值
0t2sintdt0t2
d2p(2)dpp2(p22)2
这个结果用原来的广义积分的计算方法是得不到的。
现将拉氏变换的八个性质和在实际应用中常用的一些函数的象函数分别列表如下:
表12.1 拉氏变换的性质 序号 1 设L[f(t)]F(p) L[a1f1(t)a2f2(t)]a1L[f1(t)]a2L[f2(t)] - - -专业资料-
- - -
2 3 L[eatf(t)]F(pa) L[f(ta)u(ta)]eapF(p) (a>0) L[f'(t)]pF(p)f(0) 4 L[f(n)(t)]pnF(p)[pn1f(0)pn2f'(0)…f(n1)(0)] 5 L[t0f(x)dx]F(p)p 6 L[f(at)1pF()aa (a>0) 7 L[tnf(t)](1)nF(n)(p) f(t)L[]t8 pF(p)dp
表12.2 常用函数的拉斯变换表 序号 1 f(t) F(p) 1 (t) u(t) 2 1p1p2 3 t 4 t(n1,2,...) nn!pn1 5 eat 1pa 6 1eatap(pa) 1(pa)27 teat - - -专业资料-
- - -
8 tneat(n1,2,) n!(pa)n1 9 sint cost p22 10 pp22 11 sin(t) cos(t) psincosp22pcossinp222p(p22)2 12 13 tsint sinttcost 14 23(p22)2p22(p22)2 15 tcost 16 eatsint (pa)22pa(pa)22 17 eatcost 18 19 1(1cosat)2a 1p(p2a2) ee2tatbt ab(pa)(pb)) 20 1 1pp 21 t 1p - - -专业资料-
- - -
习题12.1
1.求下列函数的拉氏变换 (1)f(t)e4t
(2)f(t)t (3)f(t)te
at2(4)f(t)sin(t)(,是常数) 2.求下列题中函数的拉氏变换 (1)3e4t(2)5sin2t3cost
(3)f(t)1,1,0t4t4(4)f(t)sint,t,0tt
0t20,2t4(6)f(t)tneat (5)f(t)1,0,4t第二节 拉普拉斯逆变换
前面我们主要讨论了怎样由已知函数f(t)求它的象函数F(p)的问题.运算法的另一面是已知象函数F(p)要求它的象原函数f(t),这就是拉斯逆变换问题.在控制工程中,求拉氏反变换的简便方法是利用拉氏变换表。同时把常用的拉氏变换的性质用逆变换形式一一列
出.
性质12.9(先行性质)
L1[a1F1(p)a2F2(p)]a1L1[F1(p)]a2L1[F2(p)]a1f1(t)a2f2(t)。
1at1atL[F(pa)]eL[F(p)]ef(t)。 性质12.10(平移性质)
1apL性质12.11(滞后性质)[eF(p)]f(ta)u(ta)。
例12.14 求F(p)解:
2p3的逆变换。
p22p52p312(p1)5]L[]p22p5(p1)24
p15122L1[]L[]22(p1)42(p1)4
p522etL1[2]etL1[2]2p4p4 552etcos2tetsin2tet[2cos2tsin2t]22
f(t)L1[在运用拉氏变换解决工程技术中的应有问题时,通常遇到的象函数常常是有理分式,对
于有理分式一般可采用部分分式方法将它分解为较为简单的分式之和,然后再利用拉氏变换表求出象原函数。
- - -专业资料-
- - -
例12.15 求F(p)p3p34p24p的逆变换。
解:先将F[p]分解为几个简单分式之和:
p3p34p24Pp3P(p2)2ABCpp2(p2)2, 用待定系数法求得A34,B34,C12,所以
331F(p)p34p34p24Pp4p22(p2)2,
于是
f(t)L1[F(p)]L1[31314p4p2112(p2)2]
3L1[1]3L11114p4[p2]2L1[(p2)2]
3432t12t4e2te
习题13.2
求下列题中函数的拉氏逆变换
1.F(p)2p32.F(p)4pp216 3.F(p)2p81p236 4. F(p)p(p1)(p2)
5.F(p)p2p21p36p29p6. F(p)p(p1)2 第三节拉氏变换在电学中的应用
一、求解常微分方程
例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代入上式,得
(p2)X(p)3
这样,原来的微分方程经过拉氏变换后,就得到了一个象函数的代数方程。
- - -专业资料-
- - -
第二步解出X(p):X(p)3 p211第三步求象函数的拉氏逆变换:x(t)L[X(p)]L[这样就得到了微分方程的解x(t)3e2t3]3e2t p2t。
例12.17有一个二阶动态电路满足微分方程y''3y'2y2e,并且其初值条件
y(0)2,y'(0)1,求其解。
解:对所给微分方程的两边分别作拉氏变换.设L[y(t)]Y(p)Y,则得
[p2Ypy(0)y(0)]3[pYy(0)]2Y2p1将初值条件y(0)2,y'(0)1,代入,得到Y
的代数方程
(p23p2)Y2p12p7
即
(p23p2)Y2p25p5p1
解出Y,得
2p2Y5p5(p1)(p2)(p1)
将上式分解为部分分式
17Y341p13pp2
再取拉氏逆变换,就得到满足所给初值条件的方程的特解为
y(t)1et4et7e2t33
用拉氏变换还可以解常系数线性微分方程组。
二、电学应用举例
例12.18求图示电路的输入运算阻抗Zin(s)
解:由串并联关系得
Z(s) =2s12s212s21inss22s1s12 例12.19求图(a)所示电路中的i(t)、uC(t)。
- - -专业资料-
- - -
(a) (b)
解:先画出运算电路如图(b)所示。由运算电路得
1K3KK2sIss11
10s1s26s10s1s3js3j6ss其中
K1Iss1s1ss26s10s115K2Iss3js3j0.1j0.712ss1s3js3j
/ 81.8712/ 81.87K3Iss3js3j0.1j0.7则
1itL1Iset2e3 tcost81.87 ε(t) A
5UCs其中
K3K1K21010Is 2ss1s6s10s1s3js3jK1 UCss1s110s26s102
s1K2 UCss3js3j10
s1s3js3j1j25/ 116.565
K3 UCss3js3j5/ 116.565
则
uCtL1UCs[2et25e3 tcos(t116.565)] ε(t)V
- - -专业资料-
- - -
习题 12.3
1.求解一输入响应电路的微分方程。
didt5i10e3t,i(0)0
2. 求图(a)所示电路中的回路电流i1和i2.。
自测题
1. 求各函数的拉氏变换
0,0t(1)
f(t)11,1t2(2)f(t)5sin2t3cos2t
2,2t(3)
f(t)8sin23t(4)f(t)1tet
2. 求各象函数的逆变换 (1)F(p)1p(p1)2 (2)F(p)3p9 p22p10(3)F(p)5p215p7p (4)F(p)2ee2p(p1)(p2)3p
3.如图所示电路激励为i(t)(t),响应为u1、u2。求阶跃响应S1(t)、S2(t)。2Ω1/4F+1Ωu+12Hu2--
4.求图示两电路的输入运算阻抗Zin(s)
- - -专业资料-
- - -
-
- -专业资料-
因篇幅问题不能全部显示,请点此查看更多更全内容