拉普拉斯变换与常系数线性微分方程

                     

贡献者: JierPeter

预备知识 常系数齐次线性微分方程,拉普拉斯变换的性质

1. 拉普拉斯变换

导数定理

   拉普拉斯变换的导数定理(见拉普拉斯变换)使得我们可以把常系数线性微分方程化为代数方程,不论是否齐次。在常系数齐次线性微分方程中,我们通过把方程写为算子 ddt 的代数方程,实现了化微分方程为代数方程,但这种操作对非齐次方程没有用。拉普拉斯变换是另一种化为代数方程的思路,而且不拘泥于齐次与否。

   引用拉普拉斯变换的性质,将导数定理誊写如下:

定理 1 拉普拉斯变换的导数定理

   设 [0,+] 上有 k 次可导函数 f(t),其拉普拉斯变换为

(1)L(f(t))=F(s)=0+f(t)estdt ,

   则有

(2)L(ddtf(t))=sF(s)f(0) .

   由定理 1 ,容易计算得

(3)L(d2dt2f(t))=sL(ddtf(t))f(0)=s2F(s)sf(0)f(0) .

   以此类推,可得:

推论 1 

  

(4)L(dkdtkf(t))=snF(s)i,j>0i+j=n1sif(j)(0) .
其中 f(j) 表示 fj 次导函数。

函数乘以一个指数函数后的拉普拉斯变换

定理 2 

   对于 f(t),若设其拉普拉斯变换为 F(s),则

(5)L(eatf(t))=F(sa) .

常见函数的拉普拉斯变换

例 1 幂函数

(6)L(tn)=n!sn+1 .

   这一事实可以通过应用分部积分得证:

(7)L(tn)=0+tnestdt=1stnest|0++ns0+tn1estdt=ns0+tn1estdt=nsL(tn1) .

例 2 指数函数

(8)L(eat)=1sa .

例 3 三角函数

   利用 eix=cosx+isinx,将 a=ωi 代入式 8 ,分开实部和虚部,可得三角函数的拉普拉斯变换:

(9){L(cosωt)=ss2+ω2L(sinωt)=ωs2+ω2 .

   有了式 6 式 8 式 9 ,再结合定理 2 ,我们就可以计算出多数情况下的拉普拉斯变换了。

2. 微分方程

   拉普拉斯变换的一大实用之处,就是可以用来解常系数线性微分方程。我们先举一个例子来说明其解法。

例 4 拉普拉斯变换解非齐次方程

   考虑方程

(10)(d2dt2+ddt+1)f(t)=t .

   设 L(f(t))=F(s),那么根据推论 1 ,给式 10 左右同时进行拉普拉斯变换,得

(11)s2F(s)sf(0)f(0)+sF(s)f(0)+F(s)=1s2 ,
其中 f(0),f(0) 都是待定的初值。

   将式 11 整理为

(12)F(s)=1s2(s2+s+1)+sf(0)+f(0)+f(0)s2+s+1 .

   注意到

(13)s2+s+1=(s+12)2+34 ,
式 12 可以分解为
(14)F(s)=s(s+12)2+34s1s2+sf(0)+f(0)+f(0)(s+12)2+34 ,

   进一步整理为

(15)F(s)=1s21s+(f(0)+1)(s+12)(s+12)2+34+f(0)+12f(0)12(s+12)2+34 .

   式 15 的右边就是我们熟知的拉普拉斯变换的形式了:

(16){L(t)=1s2L(1)=1sL(e12tcos(32t))=s+12(s+12)2+34L(e12tsin(32t))=1(s+12)2+34 .

   根据拉普拉斯变换的线性性,求出式 16 中各式的逆变换,其线性组合就是 F(s) 的逆变换:

(17)f(t)=t1+(f(0)+1)e12tcos(32t)+33(f(0)+12f(0)12)e12tsin(32t) ,

   式 17 就是式 10 的通解。只需要把初值 f(0)f(0) 代进去即可得给定初值下的特解。

   例 4 中最麻烦的一步是从式 12 式 14 的因式分解。代数学的知识告诉我们,任何一个形如 1p(t) 的函数,其中 p(t) 是多项式,都可以化为 1s+a,1(s+a)2+b,s(s+a)2+b 的线性组合,而这三个式子的拉普拉斯逆变换都是已知的。因此,当常系数线性微分方程式 1 的右端的拉普拉斯变换是多项式的倒数的时候,我们总可以求方程两端的拉普拉斯变换,得到通解的拉普拉斯变换,并求逆变换以得到通解本身。

未完成:因式分解的技巧需引用。

                     

© 小时科技 保留一切权利