拉普拉斯变换在求解微分方程中的应用

目 录

引言 ............................................................... 1 1 拉普拉斯变换以及性质 ............................................. 1 1.1 拉普拉斯变换的定义 ....................................................... 1 1.2 拉普拉斯变换的性质 ....................................................... 2 2 用拉普拉斯变换求解微分方程的一般步骤 ............................. 3 3 拉普拉斯变换在求解常微分方程中的应用 ............................. 4 3.1 初值问题与边值问题 ....................................................... 4 3.2 常系数与变系数常微分方程 ................................................ 5 3.3 含?函数的常微分方程 ..................................................... 6 3.4 常微分方程组 .............................................................. 7 3.5 拉普拉斯变换在求解非齐次微分方程特解中的应用 ........................ 7 3.6 拉普拉斯变换在求解高阶微分方程中的推广 .............................. 11 4 拉普拉斯变换在求解偏微分方程中的应用 ............................ 12 4.1 齐次与非齐次偏微分方程 ................................................. 12 4.2 有界与无界问题 .......................................................... 15 5 综合比较，归纳总结 .............................................. 19 结束语 ............................................................ 20 参考文献 .......................................................... 20 英文摘要 .......................................................... 21 致谢 .............................................................. 21

忻州师范学院物理系本科毕业论文（设计）

拉普拉斯变换在求解微分方程中的应用

物理系0801班 学 生 岳艳林 指导老师 韩新华

摘 要：拉普拉斯变换在求解微分方程中有非常重要的作用，本文首先介绍拉普拉斯

变换的定义及性质；其次给出拉普拉斯变换求解微分方程的一般步骤；然后重点举例拉普拉斯变换在求解常微分方程（初值问题与边值问题、常系数与变系数常微分方程、含?函数的常微分方程、常微分方程组、拉普拉斯变换在求解微分方程特解中的应用、拉普拉斯变换在求解高阶微分方程的推广）与典型偏微分方程（齐次与非齐次偏微分方程、有界与无界问题）中的应用举例；最后综合比较、归纳总结拉普拉斯变换在求解微分方程中的优势以及局限性。

关键词：拉普拉斯变换；拉普拉斯逆变换；常微分方程；偏微分方程；特解 引言

傅里叶变换和拉普拉斯变换是常用的积分变换，但对函数进行傅里叶变换时必须满足狄里希利和在???t???内绝对可积，但是在物理、无线电技术等实际应用中，许多以时间t为自变量的函数通常在t?0时不需要考虑或者没有意义，像这样的函数不能取傅里叶变换。为避免上述两个缺点，将函数进行适当改造，便产生了拉普拉斯变换[1]。

1 拉普拉斯变换以及性质 1.1 拉普拉斯变换的定义

??设函数f(t)当t?0时有定义，而且积分

?0f(t)e?stdt(s是一个复参量)在s的

??某一区域内收敛，则此积分所确定的函数可写为F(s)??0f(t)e?stdt.我们称上式

为函数f(t)的Laplace变换式.记为F(s)?L[f(t)],F(s)称为f(t)的Laplace变换（或称为象函数）.

若F(s)是f(t)的Laplace变换，则称f(t)为F(s)的Laplace逆变换（或称为象原函数），记为f(t)?L?1[F(s)][2].

Laplace变换的存在定理 若函数f(t)满足下列条件:

1

忻州师范学院物理系本科毕业论文（设计）

1?在t?0的任一有限区间上分段连续；

2?当t???时，f(t)的增长速度不超过某一指数函数，亦即存在常数M?0及c?0，使得f(t)?Mecｔ，0?t???成立（满足此条件的函数，称它的增大是不超过指数级的，c为它的增长指数）.

??则f(t)的Laplace变换F（s）=?f(t)e?stdt在半平面Re(s)?c上一定存在，

0右端的积分在Re(s)?c1?c的半平面内,F(s)为解析函数[2]. 1.2 拉普拉斯变换的性质

⑴线性性质 若?，?是常数，L[f1(t)]?F1(s), L[f2(t)]?F2(s),

则有L[?f1(t)??f2(t)]??L[f1(t)]+?L[f2(t)],

L?1[?F1(s)??F2(s)]??L?1[F1(s)]+?L?1[F2(s)].

⑵微分性质 若L[f(t)]?F(s),则有L[f'(t)]?sF(s)?f(0). 高阶推广 若L[f(t)]?F(s),则有L[f??(t)]?s2F(s)?sf(0)?f'(0).

一般，L[fn(t)]?snF(s)?sn?1f(0)?sn?2f'(0)?????sf(n?2)(0)?f(n?1)(0).

⑶积分性质 若L[f(t)]?F(s),则L[?t0f(t)dt]?1L[F(s)]. s⑷位移性质 若L[f(t)]?F(s)，则L[eatf(t)]?F(s?a)(Re(s?a)?c). ⑸延迟性质 若L[f(t)]?F(s),又t?0时f(t)=0,

则对于任一非负实数?，

有L[f(t??)]?e?s?F(s),或L?1[e?s?F(s)]?f(t??)[2].

⑹相似性性质 若L[f(t)]?F(s),则L[f(at)]?1sF(). aa⑺卷积性质 若L[f1(t)]?F1(s),L[f2(t)]?F2(s)，

则L[f1(t)1?f2(t)]?F1(s)F2(s)，

其中f1(t)1?f2(t)??f1(?)f2(t??)d?称为f1(t)与f2(t)的卷积[3].

0t由于从定义以及性质求拉普拉斯变换或拉普拉斯逆变换困难且复杂，在控制工程中，常常通过查阅已编好的“拉氏变换对照表”来实现。拉氏变换对照表列

2

忻州师范学院物理系本科毕业论文（设计）

出了工程上常用的时间函数及其对应的拉氏变换，可以根据该表查找原函数的象函数，或者从象函数查找原函数。对于表中不能找到的形式，可以把它展开成部分分式，再求拉普拉斯变换或拉普拉斯逆变换。以下是本文将用到的几种常用的拉普拉斯变换函数对[3]:

原函数 象函数 原函数 象函数 1 1 s1 s??tn(n为整数) 1sinat tn!sn?1 e?t arctana ssin?t ? 22s???s2??2 cos?t s 22s??s s2??2sh?t ch?t tsin?t 2?s (s2??2)21 tcos?t a2ts2??2 222(s??)1?as es?(t) erfc() 表一：拉普拉斯变换函数表

2 用拉普拉斯变换求解微分方程的一般步骤

像其他方法求解微分方程一样，应用拉普拉斯变换求解微分方程也有规范的步骤，其一般步骤[4]如下：

1、根据自变量的变化范围和方程及其定解条件的具体情况来决定对哪一个自变量进行拉普拉斯变换，然后对线性微分方程中每一项取拉普拉斯变换，使微分方程变为s的代数方程；

2、解象函数的代数方程，得到有关变量的拉普拉斯变换表达式，即象函数； 3、对象函数取拉普拉斯逆变换，得到微分方程的时域解。 流程图法[5]如下：

3

忻州师范学院物理系本科毕业论文（设计）

原函数 微分方程的解

取拉普拉斯逆变换

象函数

解代数方程 取拉普拉斯变换

微分方程

象函数的代数方程

图一：拉普拉斯变换求解微分方程的流程图

拉普拉斯变换在物理和工程等领域有着广泛的应用，通过拉普拉斯变换，可以方便地对线性控制系统进行分析、研究，可以对一些级数进行求和，还可以求解微分方程[1]。接下来重点讨论拉普拉斯变换在求解微分方程中的应用。

3 拉普拉斯变换在求解常微分方程中的应用 3.1 初值问题与边值问题

例：求解初值问题y''?4y'?3y?e?t,y(0)?y'(0)?1[2]. 解：设Y(s)?L[y(t)],对方程两边同时取拉普拉斯变换，有

[s2Y(s)?sy(0)?y'(0)]?4[sY(s)?y(0)]?3Y(s)?1, s?11, s?1结合初始条件，有[s2Y(s)?s?1]?4[sY(s)?1]?3Y(s)?s2?6s?6711131整理展开成部分分式，有Y(s)?. ??????(s?1)2(s?3)4s?12(s?1)24s?3?1?1111]?e?t,可知L[]?e?t,L[]?e?3t. s??s?1s?3?1n!由拉普拉斯变换函数表L[n?1]?tn,并结合位移性质L[e??tf(t)]?F(s??)，

s由拉普拉斯变换函数表L[?1可知L[?11?t]?te， 2(s?1)对方程两边同时求反演，整理可得方程的解为

7131y(t)?L?1[Y(s)]?e?t?te?t?e?3t?[(7?2t)e?t?3e?3t]。

4244

4

忻州师范学院物理系本科毕业论文（设计）

例：求解边值问y''?y?0,y(0)?0,y'(2?)?1[2]. 解：设Y(s)?L[y(t)],对方程两边同时取拉普拉斯变换，有

[s2Y(s)?sy(0)?y'(0)]?Y(s)?0,

结合初始条件，有[s2Y(s)?y'(0)]?Y(s)?0,

y'(0)111整理展开成部分分式，有Y(s)?2?y'(0)(?),

2s?1s?1s?1由拉普拉斯变换函数表L[?1?1?1111]?e?t,可知L[]?et,L[]?e?t. s??s?1s?1对方程两边同时求反演，整理可得方程的解为

y(t)?L?1[Y(s)]?1'y(0)(et?e?t)?y'(0)sinht. 21,

sinh2?为了确定y'(0)，将条件y(2?)?1代入上式可得y'(0)?所以，方程的解为y(t)?sinht.

sinh2?3.2 常系数与变系数常微分方程

例：求解常系数微分方程y''?2y'?y?0,y(0)?0,y(1)?2[2]. 解：设Y(s)?L[y(t)],对方程两边同时取拉普拉斯变换，有

[s2Y(s)?sy(0)?y'(0)]?2[sY(s)?y(0)]?Y(s)?0, 结合初始条件，有[s2Y(s)?y'(0)]?2[sY(s)]?Y(s)?0,

y'(0)y'(0)整理展开成部分分式，有Y(s)?2?, 2s?2s?1(s?1)由拉普拉斯变换函数表L[可知L[?1?1n!sn?1]?tn,并结合位移性质L[e??tf(t)]?F(s??),

1]?tet. 2(s?1)对方程两边同时求反演，整理可得方程的解为y(t)?L?1[Y(s)]?y'(0)tet, 为了确定y'(0)，将条件y(1)?2代入上式可得y'(0)?所以，方程的解为y(t)?L?1[Y(s)]?2, e2tte?2tet?1. e5

忻州师范学院物理系本科毕业论文（设计）

例：求解变系数微分方程

[2]. ty''?2y'?ty?0,y(0)?1,y'(0)?c0,(c0为常数）解：设Y(s)?L[y(t)],对方程两边同时取拉普拉斯变换，L[ty'']?2L[y']?L[ty]?0,

'''即L[ty]?L[y]?4L[ty]?0,

d2d[sY(s)?sy(0)?y'(0)]?2[sY(s)?y(0)]?Y(s)?0, dsdsdd两边积分可得[?2sY(s)?s2Y(s)?y(0)]?2[sY(s)?y(0)]?Y(s)?0,

dsdsdd结合初始条件，有[?2sY(s)?s2Y(s)?1]?2[sY(s)?1]?Y(s)?0,

dsdsd1整理可得Y(s)??2,

dss?1亦即?两边积分可得Y(s)??arctans?c,

1 ?arctans?arctan，

s??22sa11由拉普拉斯变换函数表L?1[arctan]?sinat,可知L?1[arctans]?sint.

stt1对方程两边同时求反演，可得方程的解为y(t)?L?1[Y(s)]?sint.

t欲求待定系数c,可利用limY(s)?0，所以从c?，Y(s)?3.3 含?函数的常微分方程

例：质量为m的物体挂在弹簧系数为k的弹簧一端，当物体在t?0时在x方向受到冲击力f(t)?A?(t)(t),其中A为常数。若物体自静止平衡位置x?0处开始运动，求该物体的运动规律x(t)[2]. 解：根据牛顿定律，有mx''?f(t)?kx,

其中?kx由胡克定律所得，是使物体回到平衡位置的弹簧的恢复力。所以，

'''mx?kx?f(t)(t?0),且x(0)?x(0)?0. 物体运动的微分方程为

??这是二阶常系数非齐次微分方程，对方程两边取拉普拉斯变换，设

L[x(t)]?X(s),L[f(t)]?L[A?(t)]?A,并考虑到初始条件，则得

ms2X(s)?kX(s)?A, 如果记?0?

2A1k. ,有X(s)??2ms??02m6

忻州师范学院物理系本科毕业论文（设计）

由拉普拉斯变换函数表L?1[?s2??]?sin?t,可知L?1[21s??02]?21?0sin?0t.

对方程两边同时取反演，从而方程的解为x(t)?Asin?0t. m?0A,角频率是?0,称?0m?0可见，在冲击力作用下，运动为一正弦振动，振幅是为该系统的自然频率（或称固有频率）。 3.4 常微分方程组 例：求解三维常微分方程组

?x''?x?y?z?0,?[2] '''''?x?y?y?z?0,x(0)?1,y(0)?z(0)?x(0)?y(0)?z(0)?0?x?y?z''?z?0,?解：设X(s)?L[x(t)],Y(s)?L[y(t)],Z(s)?L[z(t)],对方程组的两个方程两边分别取拉普拉斯变换并结合初始条件，有

?(s2?1)X(s)?Y(s)?Z(s)?0,?2?X(s)?(s?1)Y(s)?Z(s)?0 ?X(s)?Y(s)?(s2?1)Z(s)?0.?解该方程组，整理展开成部分分式，有

?s32s1sX(s)?????,?222233(s?1)(s?2)s?2s?1? ??s1s1s?Y(s)?Z(s)??????.2222?33(s?1)(s?2)s?2s?1?取其逆变换，可得原方程组的解

21?x(t)?cosh(2t)?cost,??33 ??y(t)?z(t)??1cosh(2t)?1cost.?33?3.5 拉普拉斯变换在求解非齐次微分方程特解中的应用

形如y(n)?a1y(n?1)?????an?1y'?any?f(x)的方程称为n阶常系数非齐次线性微分方程，这里a1,a2,???,an?1,an为常数，f(x)为连续函数。我们平时用到的

7

忻州师范学院物理系本科毕业论文（设计）

f(x)主要有三种形式：f(x)?e?x,

f(x)?e?xp(x)(其中p(x)?p1x?p2x2???pnxn),

f(x)?sin?x、f(x)?cos?x[6].

该非齐次微分方程的解即该非齐次微分方程的特解与对应的齐次微分方程的通解。对于该方程的通解可用多种方法求特解，如：比较系数法、常数变易法、算子法等。下面将用拉普拉斯变换法求解该方程的特解。

设Y(s)?L[y(t)],F(s)?L[f(x)],为求特解令初始条件为零，对方程两边同时取拉普拉斯变换，得到Y(s)?式分别作介绍。 （1）f(x)?e?x 此时，Y(s)?1(s??)(s?a1snn?1F(s)，下面结合f(x)的三种形nn?1s?a1s???an?1s?an???an?1s?an)，

ABsn?1?Csn?2???D对其进行部分分式分解，令Y(s)?， ?nn?1s??(s?a1s???an?1s?an)则该齐次微分方程特解的形式与自由项f(x)有关，也就是说与变换项

A有关；s??Bsn?1?Csn?2???D对应的齐次微分方程的通解由n决定，只要该项分母中不n?1(s?a1s???an?1s?an)含有特解因子s??，则特解只取决于若sn?a1sn?1???an?1s?an则A?(s??)Y?(x)?s??As??[7]

。

?0,

1???an?1s?an)(s?a1snn?1s??，

即相应的拉普拉斯变换特解为Y?(x)?11?[ns??(s?a1sn?1???an?1s?an)s??].

对方程两边同时求反演，整理可得原微分方程的特解为y?(t)?L?1[Y?(s)]. 例：求解常系数线性齐次方程y''?y'?e2x的特解。

8

忻州师范学院物理系本科毕业论文（设计）

解：设Y(s)?L[y(t)],令初始条件为零， 对方程两边同时取拉普拉斯变换，有(s2?s)Y(s)?整理展开成部分分式，有Y(s)?此时(s2?s)Y?(x)?1. s?21ABs?C??, 22(s?2)(s?s)s?2s?ss?2?0,则

s??]?11?[ns??(s?a1sn?1???an?1s?an)11?[2s?2s?s11]??, s?22s?2对方程两边同时求反演，整理可得原微分方程的特解为

111y?(t)?L?1[Y?(s)]?L?1[?]?e2x

2s?22若sn?a1sn?1???an?1s?ans???[(s??)ms??](sn?m?b1sn?m?1?bn?m)?0,

A'B'sn?m?1?C'sn?m?2???D'令Y(s)?， ?(s??)m?1(sn?m?b1sn?m?1???bn?m)同理，相应的拉普拉斯变换特解为

Y?(x)?11?[(s??)m?1(sn?m?b1sn?m?1???bn?m)s??].

例：求解常系数线性齐次方程y'''?5y''?8y'?4y?e2x的特解。 解：设Y(s)?L[y(t)],令初始条件为零，

对方程两边同时取拉普拉斯变换，有(s3?5s2?8s?4)Y(s)?则Y(s)?11?, 322(s?2)(s?5s?8s?4)(s?2)(s?2)(s?1)1. s?2此时s3?5s2?8s?4s?2?0,

A'A's?B'?, 令Y(s)?3(s?1)(s?2)则相应的拉普拉斯变换特解为

Y?(x)?11?[(s??)m?1(sn?m?b1sn?m?1???bn?m)s??]?11.[(s?2)3s?1s?2]?1, 3(s?2)

9

忻州师范学院物理系本科毕业论文（设计）

对方程两边同时求反演，整理可得原微分方程的特解为

y?(t)?L?1[Y?(s)]?L?1[1122x]?xe. 2(s?2)3（2）f(x)?e?xp(x)(其中p(x)?p1x?p2x2???pnxn). 例：求微分方程y''?5y'?6y?xe2x的特解。 解：设Y(s)?L[y(t)],令初始条件为零，

对方程两边同时取拉普拉斯变换，有(s2?5s?6)Y(s)?11?, 222(s?2)(s?5s?6)(s?2)(s?2)(s?3)1, 2(s?2)则Y(s)?此时s2?5s?6s?2?0, 令Y(s)?As?BCs?D?, 22(s?2)s?5s?6As?B?Y(s)?(s?2)2??1?ss2?3s?2s2?4s??41s2?5s?61?s?,s?2s?4s??42?12?ss?4s??42?1?s(2?s)(1?s)s2?4s??4

相应的拉普拉斯变换特解为Y?(s)?As?B?2(s?2)11?s11????,

(s?2)2s?2(s?2)2(s?2)3对方程两边同时求反演，整理可得原微分方程的特解为

y?(t)?L?1[Y?(s)]?L?1[?1112x12x2x?]?(?xe?xe)??xe(1?x).

22(s?2)2(s?2)3（3）f(x)?sin?x、f(x)?cos?x

例：求解微分方程y''?4y'?5y?sin2x的特解[7]。 解：设Y(s)?L[y(t)],令初始条件为零，

对方程两边同时取拉普拉斯变换，有(s2?4s?5)Y(s)?令Y(s)?2, s2?4As?BCs?D?, 22s?4s?4s?510

忻州师范学院物理系本科毕业论文（设计）

As?B?Y(s)(s2?4)??2(4s?1)16s2?1s2??42s2?4s?52(4s?1)?,?65s??42?24s?1s??42?2(4s?1)(4s?1)(4s?1)s2??4

2(4s?1)?B1相应的拉普拉斯变换特解为Y?(s)?As???(8?22s?4?65(s?4)65s2?), 22s?4s?4对方程两边同时求反演，整理可得原微分方程的特解为

y?(t)?L?1[Y?(s)]?L?1[8?s21?]??(8cos2x?sin2x). 22s?4s?4653.6 拉普拉斯变换在求解高阶微分方程中的推广

对于n阶常系数线性齐次微分方程y(n)?a1y(n?1)?????an?1y'?any?0满足以下两个引理[8]：

引理1 n阶常系数线性齐次方程的解（积分曲线）具有平移不变性。也就是说，若y=y(x)为n阶常系数线性齐次方程的一个解，则对任意的常数c,y?y(x?c)也是n阶常系数线性齐次方程的解。

引理2 若y?y(x,x0,y0)为n阶常系数线性齐次方程的一个解，y?y(x,x0,y0)经平移后变为y?y(x?x0,0,y0),则y?y(x?x0,0,y0)也是n阶常系数线性齐次方程的解。

下面给出利用拉普拉斯变换方法求解三阶常系数线性齐次方程

y(3)?py''?qy'?ry?0满足在任意点的初始条件 y(x0)?y0,y'(x0)?y0,y''(x0)?y0的解。

设方程的解为y?y(x,x0,y0)?y(x?x0,0,y0),这样，我们便将初值点平移到了x?x0?0点，于是可用如下的拉普拉斯变换方法求解该初值问题。 令y(t)?y(x?x0,0,y0)(其中t?x?x0),

12y(0)?y0,y'(0)?y0,y''(0)?y0,y(3）(0)?y0.

设Y(s)?L[y(t)],对方程两边同时取拉普拉斯变换，得到L[y(3)?py''?qy'?ry]?0, 由拉普拉斯变换的导数性质L[f'(t)]?sF(s)?f(0)以及

高阶导数推广L[fn(t)]?snF(s)?sn?1f(0)?sn?2f'(0)?????sf(n?2)(0)?f(n?1)(0)可得，

'(2)(3)11

忻州师范学院物理系本科毕业论文（设计）

[s3Y(s)?s2y(0)?sy'(0)?y''(0)]?p[s2Y(s)?sy(0)?y'(0)]?q[sY(s)?y(0)]?rY(s)?0. 结合初始条件，有

[s3Y(s)?s2y0?sy0?y0]?p[s2Y(s)?sy0?y0]?q[sY(s)?y0]?rY(s)?0. 整理可得Y(s)?1s3?ps2?qs?r[(s2?ps?q)y0?(s?p)y0?y0].

12121对上式两边同时取拉普拉斯逆变换，可得

L?1[Y(s)]?L?1[1s3?ps2?qs?r{(s2?ps?q)y0?(s?p)y0?y0}].

12进行变量还原，便得到所求初值问题的解为

y?y(x,x0,y0)?y(x?x0,0,y0)?y(t)?y(x?x0).

例：求解二阶常系数线性齐次方程y''?y?0，该方程满足初始条件

y()?1,y'()??1[8] 44解：首先转化初值条件y?y(x,???4,1)?y(x??4,0,1)?y(t)(其中t?x??4).

设Y(s)?L[y(t)],对方程两边同时取拉普拉斯变换，得到L[y''?y]?0, 即[s2Y(s)?s?1]?Y(s)?0. 整理成部分分式，有Y(s)?由拉普拉斯变换函数表L[?1s?1s1??. 222s?1s?1s?1?1ss]?cos?t,L[]?cost, 可知222s??s?1?1?1?1]?sin?t,L[]?sint, 由拉普拉斯变换函数表L[2可知

s??2s2?1对方程两边同时求反演，整理可得方程的解为y(t)?L?1[Y(s)]?cost?sint. 变量还原，得到原初值问题的解为

y?y(x,

?4,1)?y(x??4,0,1)?y(t)?cost?sint?cos(x??4)?sin(x??4)?2cosx.

4 拉普拉斯变换在求解偏微分方程中的应用 4.1 齐次与非齐次偏微分方程 例：求解齐次偏微分方程

12

忻州师范学院物理系本科毕业论文（设计）

??2u2?xy,(x?0,y???),??x?y??2?uy?0?x,??ux?0?3y.?[2] ?解：对该定解问题关于y取拉普拉斯变换，并利用微分性质及初始条件可得

L[u(x,y)]?U(x,s),

L[?u]?sU(x,s)?u(x,0)?sU?x2, ?y?sdU?2x, dx?2u??u?u?uL[]?L[()]?sL[]??x?y?y?x?x?xy?0x2L[xy]?2,

s2L[u]x?0?Ux?0?3. s2这样，原定解问题转化为含参数s的一阶常系数线性非齐次微分方程的边值问题：

?dUx2s?2x?2,??dxs ?3?U?.2x?0?s?dUx2dUx2方程s?2x?2可转化为s?2x?2

dxdxssx3x2解此微分方程，可得其通解为U?3??c,其中c为常数。

s3s为了确定常数c，将边界条件Ux?0?33c?. 代入上式，可得22ssx3x23?2. 所以，U(x,s)?3?ss3s2?1x1由拉普拉斯变换函数表L[]?1,可知L[]?x2.

ss?1

13

忻州师范学院物理系本科毕业论文（设计）

x3x32?13y,L[2]?3y. 由拉普拉斯变换函数表L[n?1]?t,可知L[3]?ss23s?1n!n?1方程两边取反演，从而原定解问题的解为

x3y2u(x,y)?L[U(x,s)]??3y?x2.

6?1例：求解非齐次偏微分方程

2??2u2?u(x?0,t?0),?2?a2?g,(g为常数），?t?x??u?u?0,?0,?t?0?tt?0??u?0.?x?0[2] ?解：对该问题关于t取拉普拉斯变换，，并利用微分性质及初始条件可得

L[u(x,t)]?U(x,s),

?2u?uL[2]?s2U(x,s)?sut?0??t?t?s2U,

t?oL[g]?g,s

?2u?2d2L[2]?2L[u(x,t)]?2U,?x?xdx

L[ux?0]?Ux?0?0.

这样，原定解问题转化为含参数s的二阶常系数线性非齐次微分方程的边值问题：

?d2U121g?sU??,?222sdxaa??U?0,limU?0.s???x?0

d2U121gd2U11g方程2?2sU??2可转化为2?2s2U?2,

dxaasdxaas解此微分方程，可得其通解为U(x,s)?c1e为了确定常数c1,c2,将边界条件U

x?0sxa?c2es?xa?g,其中c1,c2为常数。 3s?0,limU?0代入上式，

s??14

忻州师范学院物理系本科毕业论文（设计）

可得c1?0,c2??g, s3sx?xggg?asa所以，U(x,s)?3(1?e)?3?3e.

sss由拉普拉斯变换函数表L[由拉普拉斯变换函数表L[?1?1?1n!ssn!]?tn,可知L[n?1n?1?1gg2]?t.

2s3]?tn,并结合延迟定理L?1[e?st0F(s)]?f(t?t0),

g?asgxx可知L[3e]?(t?)2u(t?).

2aasx方程两边取反演，从而原定解问题的解为

gg?sggxxu(x,t)?L[U(x,s)]?L[3?3ea]?t2?(t?)2u(t?).

22aass?1?1x（或）

x?g2t,t?;??2a???gt2?g(t?x)2,t?x.?2aa ?24.2 有界与无界问题 例：求解有界偏微分方程

2??2u2?u?2?a2,(0?x?l,t?0),?x??t??ux?0?0,ux?l??(t),??u?0,?u?0.t?0??tt?0?[2]

解：对该定解问题关于t取拉普拉斯变换，记

L[u(x,t)]?U(x,s),

?2u?uL[2]?s2U?sux?o??t?t?2ud2UL[2]?, ?xdx2t?0?s2U,

15

忻州师范学院物理系本科毕业论文（设计）

L[ux?o]?UL[ux?l?Ux?0?0,

x?l??(s).

这样，原定解问题转化为含参数s的二阶常系数线性齐次微分方程的边值问题：

?d2Us2 ?2?2U?0,a?dx?U?x?0?0,Ux?l??(s).该方程的通解为U(x,s)?c1esxa?c2ex?0s?xa,其中c1,c2是常数。

为确定常数c1,c2，将边界条件U将边界条件U因此c1??c2?从而

U(x,s)??(s)??e??(s)???s?(l?x)a?0,代入上式，可得c1?c2?0,即c1??c2;

slax?l??(s)代入上式，可得?(s)?c1e?c2es?la.

?(s)e?eslas?la.

esxasla?es?xas?la??(s)(esxasla?es?xas?la)(e)(e.s?la?e?e?3sla))e?e?e?(e?e?es?(3l?x)as?la?3slas?(l?x)a?e?s?(3l?x)a

1?e4lsa1?e4lsa为了求U(x,s)的拉普拉斯逆变换，注意到分母为1?e期为

?4lsa,所以逆变换u(x,t)是周

4l?(s)的关于ｌ的周期函数。根据周期函数的拉普拉斯变换式，其中表4l?sa1?ea4l明?(t)是以为周期的周期函数，即

aL[?(t)]??(s)1?e?4lsa?11?e?4lsa?4la0?(?)e?s?d?,

由拉普拉斯变换函数表L?1[?(s)1?e?4lsa]??(t),

并结合延迟定理L?1[e?st0F(s)]?f(t?t0), 可知L[?1?(s)1?e?4lsa?e?l?xsa]??(t?l?xl?x)u(t?). aa16

忻州师范学院物理系本科毕业论文（设计）

同理可知

L[?1?(s)1?e?4lsa?e?l?xsa]??(t?l?xl?x)u(t?).aa3l?x3l?x)u(t?).aaL[?1?(s)1?e?4lsa?e?3l?xsa]??(t?L[?1?(s)1?e?4lsa?e?3l?xsa]??(t?3l?x3l?x)u(t?). aa方程两边取反演，从而原定解问题的解为

u(x,t)?L?1[U(x,s)]??(t?l?xl?xl?xl?x)u(t?)??(t?)u(t?)aaaa

3l?x3l?x3l?x3l?x??(t?)u(t?)??(t?)u(t?).aaaa其中u(a)为单位阶跃函数，

?0,a?0,即u(a)??

1,a?0.?例：求解无界偏微分方程

2??u2?u??a2?hu,(h为常数），（x?0,t?0),?t?x???ux?0?u0(常数），??ut?0?0.?[2] ?解：对该问题关于t取拉普拉斯变换，记

L[u(x,t)]?U(x,s),

L[?u]?sU(x,s)?u?tx?0?sU,

?2u?2d2UL[2]?2L[u(x,t)]?2, ?x?xdxL[ux?0]?U?u0. sx?0这样，原定界问题转化为含参数s的二阶常系数线性齐次微分方程的边值问题：

17

忻州师范学院物理系本科毕业论文（设计）

?d2Us?h?2U?0,??dx2a ?u0?U?,limU?0(为自然定解条件）.x?0?x??s?解此微分方程可得通解为

U(x,s)?c1es?hxa?c2e?s?hxa,其中c1,c2为常数。

x?0为确定常数c1,c2，将边界条件U?u0u代入上式，可得c1?c2?0;

ss将边界条件limU?0代入上式，可得c1?0.

x??因此，c2?u0. ss?hxau?所以，U(x,s)?0es?1.

s?hxau0?从而U(x,t)?L[U(x,s）]?L[es?1]，

u1由拉普拉斯变换函数表L?1[]?1，可知L?1[0]?u0。

ss1a2)?由拉普拉斯变换函数表L?1[e?as]?erfc(s2t?xs1?ax2]?erfc()?可知L[es2at??1???a2te??d?，

2???x2ate??d?.

2如果令f(t)?2?'???xe??d?,显然f(0)?0,

22atxs1?a由导数性质L[f(t)]?sF(s)?f(0)可知f(t)?L[s?e],

s'?1亦即L?1[e?xsa]?f'(t)??2?e?(x2at)2?dxx()?edt2at2at?t?x24a2t,

由位移性质L[e??tf(t)]?F(s??), 可知L[e?1s?h?xa]?x2at?te?x24a2t?e?ht?x2at?te?(x24a2t?ht).

18

忻州师范学院物理系本科毕业论文（设计）

由卷积定理L[f1(t)1?f2(t)]?F1(s)F2(s),

?u可得U(x,t)?L[0]?L?1[es?1s?hxa],

令??x2a??1,最后可得该定解问题的解为

s?hxax24a2tuu(x,t)?L[0]?L?1[es??t?]?u0?x2at?td??e?(?ht)u0(?)x2a(t??)?(t??)0e?[x24a2(t??)?h(t??)]2u0????xe?(??2hx24a2?2

)d?.2at

5 综合比较，归纳总结

从以上的例题可以看出，用拉普拉斯变换方法求解微分方程有如下的优缺点

[1~13]

：

⑴拉普拉斯变换对像函数要求比傅里叶变换弱，其使用面更宽。但拉普拉斯变换像其他变换一样都有其局限性，只有满足其存在定理时才可以使用拉普拉斯变换。而在微分方程的一般解法中，并没有任何限制；

⑵用拉普拉斯变换方法求解微分方程，由于同时考虑初始条件，求出的结果便是需要的特解。而微分方程的一般解法中，先求通解，再考虑初始条件确定任意常数，从而求出特解的过程比较复杂；

⑶零初始条件、零边界条件使得拉普拉斯变换方法求解微分方程更加简单。而在微分方程的一般解法中，不会因此而有任何简化；

⑷用拉普拉斯变换求解微分方程，对于自变量是零的初始条件，求其特解是非常方便的。但微分方程的一般解法并没有简化；

⑸用拉普拉斯变换方法求解微分方程，对方程的系数可变与否、对区域有界与否、对方程和边界条件齐次与否并无特殊关系。而在微分方程的一般解法中，会遇到很多困难；

⑹用拉普拉斯变换方法求解微分方程组，可以在不知道其余未知函数的情况下单独求出某一个未知函数。但在微分方程的一般解法中通常是不可能的；

⑺拉普拉斯变换可以使解n个自变量偏微分方程的问题，转化为解n?1个自变量的微分方程的问题，逐次使用拉普拉斯变换，自变量会逐个减少，有时还可

19

忻州师范学院物理系本科毕业论文（设计）

将解n个自变量偏微分方程的问题最终转化为解一个常微分方程的问题，比微分方程的一般解法更为简单、直接；

⑻比较系数法和常数变易法只需进行代数运算和积分运算，要求相对较低。相比之下，算子法要先将方程化为算子形式然后利用算子的性质进行分解，对初学者而言要求相对较高，然而算子法却具备比较系数法和常数变易法无法具备的应用条件，有适应面广、计算量小、准确度高、简单易行的特点。 结束语

通过列举拉普拉斯变换在求解微分方程中的应用，可以看出拉普拉斯变换是一种特别成功的数学方法，求解微分方程的步骤比较明确、规律性比较强、思路清晰且容易掌握。灵活使用拉普拉斯变换，可以巧妙地推出一些复杂问题的答案，便于学生理解进而提高教学质量。

参考文献

[1] 李高翔.拉普拉斯变换在微分方程组求解中的应用[J] .高等函授学报,2009，22(3):22-24.

[2] 张元林.工程数学积分变换（第四版）[M].北京:高等教育出版社,2003：68-138. [3] 梁昆淼.数学物理方法(第三版)[M].北京:高等教育出版社,1998：120-121.

[4] 黄会芸.拉普拉斯变换在高等数学中的应用[J].潍坊教育学院学报,2009,22(4):44-45. [5] 全生寅.论解N阶常微分方程的Laplace变换法[J].青海大学学报,2000(5):61-62. [6] 李曼生,陈莉.拉普拉斯变换在求解微分方程中的应用[J].广西右江民族师专学报,2006(3):5-8.

[7] 张刁民.拉普拉斯变换求二阶常系数非齐次微分方程的特解[J].河南教育学院学报,2005(1):27-28.

[8] 李连忠,何乐亮.拉普拉斯变换应用的一个推广[J].山东师范大学学报,2007(1):148. [9] 姜立新.Laplace变换的应用研究[J].枣庄学院学报,2010,27(2).37-40.

[10] 唐妍霞.利用Laplace变换求解一维波动方程的定解问题[J].河北北方学院学报,2010(3):16-19.

[11] 王振芳.拉普拉斯变换及其应用[J.]雁北师范学院学报,2001(6):48-49.

[12] 谢小良.基于Laplace变换下微分方程的解法及应用[J].湖南城市学院学报（自然科学）.2003,24(3).85-86.

20

忻州师范学院物理系本科毕业论文（设计）

[13] 杨芳,吴小欢.n阶常系数非齐次线性微分方程特解的求解方法[J].广西师范学院学报,2009(04):97-100.

Application of Laplace Transform to General Solutions of Differential Equations

Department of Physics 0801 Student Yanlin Yue

Tutor Xinhua Han

Abstract: Through to the Laplace transform in solving ordinary differential equation, the typical application of partial differential equation, for example, comprehensive comparison, summarizes the Laplace transform in solving differential equations and the advantage of the limitations.

Key words: Laplace transform; Laplace inverse transform; ordinary differential equation; partial differential equation; particular integral 致谢

感谢我的导师韩新华老师，她渊博的专业知识，言谨的治学态度，精益求精的工作作风，诲人不倦的高尚师德，严以律己、宽以待人的崇高风范，朴实无华、平易近人的人格魅力对我影响深远。不仅使我树立了远大的学术目标、掌握了基本的研究方法，还使我明白了许多待人接物与为人处事的道理。本论文从选题到完成，每一步都是在导师的细心指导下完成的，倾注了导师大量的心血。在此，谨向导师表示崇高的敬意和衷心的感谢！在论文即将完成之际，我的心情无法平静，本论文顺利完成，还有许多可敬的师长、同学、朋友给了我无言的帮助，在这里请接受我诚挚的谢意！

21

本文来源：https://www.bwwdw.com/article/53bf.html