拉普拉斯变换及逆变换

第十二章 拉普拉斯变换及逆变换

拉普拉斯(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)连续,则

=

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)

第十二章 拉普拉斯变换及逆变换

拉普拉斯(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)连续,则

=

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)


相关内容

  • 阅读材料拉普拉斯变换法
  • 第十三章 拉普拉斯变换 一.教学基本要求 1.了解拉普拉斯变换的定义,会用拉普拉斯变换基本性质求象函数. 2.掌握求拉普拉斯反变换的部分分式展开法,. 基尔霍夫定律的运算形式. 运算阻抗和运算导纳.运算电路. 3.掌握应用拉普拉斯变换分析线性电路的方法和步骤. 二.教学重点与难点 教学重点:1. 拉 ...

  • Unit11 传递函数和拉普拉斯变换
  • Unit11 传递函数和拉普拉斯变换 传递函数 若已知图2-1所示线性系统的输入输出关系,该系统的特点也就知道了.在拉普拉斯域里输入输出关系叫做传递函数(TF 或者G 增益).通过定义,一个元件或系统的传递函数是输入的拉氏变换形式与输出的拉氏变换形式的配给量: 公式2-1 这个传递函数的定义要求系统 ...

  • 4-23 -高阶线性方程的例子.拉普拉斯变换
  • 4.3 高阶线性微分方程例子.拉普拉斯变换 (Examples of higher order Linear ODE.Laplace Transformation ) [教学内容] 1. 介绍简谐振动.节拍.共振.振幅反应等概念. 2. 介绍拉普拉斯变换的概念和性质.3. 介绍线性高阶微分方程的拉普 ...

  • 信号与线性系统
  • <信号与线性系统>课程教学大纲 课程编号:28121008 课程类别:学科基础课程 授课对象:信息工程.电子信息工程.通信工程等专业 指定教材:管致中,<信号与线性系统>(第4版),高等教育出版社,2004年 教学目的: <信号与线性系统>课程讨论确定信号经过线性 ...

  • 拉普拉斯变换的实际应用
  • 拉普拉斯变换的实际应用 在工程学上的应用 应用拉普拉斯变换解常变量齐次微分方程,可以将微分方程化为代数方程,使问题得以解决.在工程学上,拉普拉斯变换的重大意义在于:将一个信号从时域上,转换为复频域(s域)上来表示:在线性系统,控制自动化上都有广泛的应用. 拉氏变换在微分方程(组)初值问题中的应用 1 ...

  • 傅里叶变换在图像处理中的作用
  • 傅立叶变换在图像处理中的作用 (2011-05-21 20:01:34) 转载▼ 标签: 分类: 学习历程 杂谈 从现代数学的眼光来看,傅里叶变换是一种特殊的积分变换.它能将满足一定条件的某个函数表示成正弦基函数的线性组合或者积分.在不同的研究领域,傅里叶变换具有多种不同的变体形式,如连续傅里叶变换 ...

  • 信号与系统考试考点
  • 信号与系统的时域分析 1. 什么是LTI 系统?在时域中,我们如何表示系统?什么是系统的单位冲激响应? ◆ 系统的线性时不变性的证明与判断(书中例子1-14,1-16,1-17,1-18): ◆ 表示系统的时域数学模型:卷积表示,微分方程(连续时间系统),差分方程(离散时间 系统): ◆ 单位冲激响 ...

  • 拉普拉斯变换及其应用
  • 第17卷 第6期 2001年12月雁 北 师 范 学 院 学 报JO U RN AL O F Y AN BEI T EACHERS ' CO L L EGE V o l. 17No. 6Dec. 2001文章编号:1009-1939(2001) 06-0048-02 拉普拉斯变换及其应用 王振芳 ( ...

  • 试井分析中的离散数据拉普拉斯变换方法
  • 开发试采 天 然 气 工 业 2004年3月 试井分析中的离散数据拉普拉斯变换方法 张杰 李翠楠 张向阳 李治平 1 2 2 3 (1. 西南石油学院 2. 中国石油西南油气田分公司重庆气矿 3. 中国地质大学 北京) 张杰等. 试井分析中的离散数据拉普拉斯变换方法. 天然气工业, 2004; 24 ...

  • 拉普拉斯变换的基本定理
  • 本节介绍拉普拉斯变换(也称为拉氏变换)的基本性质,了解掌握了这些性质,可以更加方便地求解各种拉普拉斯正反变换. 一.线性定理 设 则: (式9-2-1) 式中 为常系数. 例9-2-1  求 . 和 的拉氏变换. 解: 同理: 二.微分定理 设 ,则: (式9-2-1) 同理可推广得到 的高阶导数的 ...