二阶常微分方程解存在的问题

二阶常微分方程解的存在问题分析
摘要
本文首先介绍了二阶常系数齐次线性微分方程的一般解法——特征方程法及二阶常系数非齐次线性微分方程的待定系数法，然后又介绍了一些可降阶的微分方程类型。

接着，讨论了二阶变系数微分方程的幂级数解法并论述了如何利用变量代换法将某些变系数方程化为常系数方程。

另外，本文还介绍了求解初值问题的另一种方法——拉普拉斯变换法。

最后，给出了二阶微分方程的存在唯一性定理的证明以及它在科学研究、工程技术以及数学建模中解决实际问题的一些应用。

1.引言
1.1常微分方程的发展过程与研究途径
二阶线性微分方程是常微分方程中一类很重要的方程。

这不仅是因为其一般理论已经研究地比较清楚，而且还因为它是研究非线性微分方程的基础，在工程技术和自然科学中有着广泛的应用。

在科学研究、工程技术中，常常需要将某些实际问题转化为二阶常微分方程问题。

因此，研究不同类型的二阶常微分方程的求解方法及探讨其解的存在唯一性问题是十分重要的。

常微分方程已有悠久的历史，而且继续保持着进一步发展的活力，主要原因是它的根源深扎在各种实际问题之中。

牛顿最早采用数学方法研究二体问题，其中需要求解的运动方程就是常微分方程。

他把两个物体都理想化为质点，得到3个未知函数的3个二阶方程组，经简单计算证明，可化为平面问题，即两个未知函数的两个二阶微分方程组。

用现在叫做“首次积分”的办法，完全解决了它的求解问题。

17世纪就提出了弹性问题，这类问题导致悬链线方程、振动弦的方程等等。

20世纪30年代直至现在，是常微分方程各个领城迅速发展、形成各自相对独立的而又紧密联在一起的分支学科的时期。

1927－1945年间定性理论的研究主要是跟无线电技术联系在一起的。

第二次世界大战期间由于通讯等方面的要求越来越高，大大地激发了对无线电技术的研究，特别是非线性振动理论的研究得到了迅速的发展。

40年代后数学家们的注意力主要集中在抽象动力系统的拓扑特征, 如闭轨是否存在、结构是否稳定等, 对于二维系统已证明可以通过奇点及一些特殊的闭轨和集合来判断结构稳定性与否；而对于一般系统这个问题尚未解决。

在动力系统理论方面, 我国著名数学家廖山涛教授, 用从典范方程组到阻碍集一整套理论和方法, 解决了一系列主要问题, 特别是C’封闭引理的证明, 对结构稳定性的充要条件等方面都作出了主要贡献。

1.2问题的研究现状
在当代由电力网、城市交通网、自动运输网、数字通讯网、灵活批量生产网、复杂的工业系统、指令控制系统等提出大系统的数学模型是常微分方程组描述的。

对这些系统的稳定性研究, 引起了越来越多学者的兴趣, 但目前得到的成果仍然只是初步的。

常微分方程的概念、解法和其它理论很多，比如，方程和方程组的种类及解法、解的存在性和唯一性、奇解、定性理论等等。

下面就方程解的有关几点简述一下，以了解常微分方程的特点。

求通解在历史上曾作为微分方程的主要目标，一旦求出通解的表达式，就容易从中得到问题所需要的特解。

也可以由通解的表达式，了解对某些参数的依赖情况，便于参数取值适宜，使它对应的解具有所需要的性能，还有助于进行关于解的其他研究。

后来的发展表明，能够求出通解的情况不多，在实际应用中所需要的多是求满足某种指定条件的特解。

当然，通解是有助于研究解的属性的，但是人们已把研究重点转移到定解问题上来。

一个常微分方程是不是有特解呢？如果有，又有几个呢？这是微分方程论中一个基本的问题，数学家把它归纳成基本定理，叫做存在和唯一性定理。

存在和唯一性定理对于微分方程的求解是十分重要的。

由于大部分的常微分方程求不出十分精确的解，而只能得到近似解。

当然，这个近似解的精确程度是比较高的。

微分方程的近似解法（包括数值解法）具有十分重要的实际意义，而解的存在和唯一又是进行近似计算的前提。

因为如果解根本不存在，却要去近似地求它，问题本身是没有意义的；如果有解存在而不唯一，由于不知道要确定是哪一个解，却要去近似地确定它，问题也是不明确的。

解的存在唯一性定理保证了所要求的解的存在和唯一，因此它也是近似求解的前提和理论基础。

此外，我们将看到在定理的证明中还具体地提出了求近似解的途径，这就更增添了存在唯一性定理的实用意义。

由于种种条件的限制，实际测出的初始数据往往是不精确的，它只能近似地反映初始状态。

因此我们以它作为初值条件所得到的解是否能用做真正的解呢？这就产生了解对初值的连续依赖性问题，即当初值微小变动时，方程的解的变化是否也是很小呢？如果不然的话，这样所求得的解就失去了实用的意义，因为它可能与实际情况产生很大的误差。

在科学研究、工程技术中，常常需要将某些实际问题转化为二阶常微分方程问题，因此，研究不同类型的二阶常微分方程的求解方法及探讨其解的存在唯一性问题，是十分重要的。

1.3问题研究存在的不足与前景
现今对于二阶线性微分方程的研究已经取得了不少成就，尤其在二阶常系数线性微分方程的求解问题和解的存在唯一性定理等方面卓有成效。

二阶微分方程的解的存在唯一性定理不仅可判断解的存在唯一性，而且还有着广泛的应用。

而幂级数解法作为求解二阶变系数齐次线性微分方程的一种方法，其过程还是比较繁琐的，计算量偏大，且需要考虑函数是否解析，幂级数在某个区间是否收敛等。

另外，对于二阶变系数非齐次微分方程，目前还尚有通用的求解方法，只有一些特殊类型是可以求解的。

应该说，应用常微分方程理论已经取得了很大的成就，但是，它的现有理论也还远远不能满足需要，还有待于进一步的发展，使这门学科的理论更加完善。

2.常系数线性微分方程的解法
2.1 二阶常系数齐次线性微分方程的解法——特征方程法
若21,y y 是二阶常系数齐次线性微分方程
0'''=++qy py y ，其中q p ,均为常数 （2.1）
的两个线性无关的解，那么（2.1）的通解就可表示成
2211y C y C y +=（21,C C 为任意常数）
由此可知，只要找到方程（2.1）的两个线性无关的解，就能求出（2.1）的通解。

我们知道，当r 为常数时，函数rx e y =和它的各阶导数只相差一个常数。

因此，可以设想（2.1）有形如rx e y =的解，将rx e y =代入方程（2.1）得：
0)(2=++q pr r e rx
又0≠rx e ，则必有
02=++q pr r （2.2）
即如果rx e y =是（2.1）的解，则r 必满足方程（2.2）.
反之，若r 满足方程（2.2），则rx e y =就是（2.1）的一个特解。

我们称方程（2.2）是方程（2.1）的特征方程，它的根就称为特征根，且特征根
2
422,1q p p r -±-=. 下面根据特征根的不同情况分别进行讨论。

1）有两个不相等的实根)0(>∆:
2
421q p p r -+-=，2422q p p r ---= 易知x r e y 1=和x r e y 2=是方程（2.1）的两个线性无关的特解，则方程（2.1）的通解为：
x r x r e C e C y 2121+=；
2）有两个相等的实根)0(=∆:
2
21p r r -== 易知x r e y 1=是方程（2.1）的一个特解，设另一特解为x r e x C y 1)(2=，将2y 代入到（2.1）得：
0)(')2(''1211=+++++C q pr r C p r C
（2.3) 又2
21p r r -==，q p 42=，则可得0''=C ，不妨取x x C =)(，代入（2.3）得： x r xe y 12=，则方程（2.1）的通解为：
x r e x C C y 1)(21+= ；
3）有一对共轭复根)0(<∆：
βαi r +=1，βαi r -=2
易知x i e
y )(~1βα+=与x i e y )(~
2βα-=是方程（2.1）的两个线性无关的复值解。

而 )sin (cos )(x i x e e x x i ββαβα+=+，)sin (cos )(x i x e e x x i ββαβα-=-
若取
x e y y y x βαcos )(21~2~11=+=，x e y y y x βαsin )(2
1~2~12=-= 由解的叠加性知，21,y y 也是方程（2.1）的两个特解，又
常数≠==x x
e x e y y x x βββααcot sin cos 21， 于是，21,y y 就是方程（2.1）的两个线性无关的实值解。

从而方程（2.1）的通解为：
)sin cos (21x C x C e y x ββα+=。

下面举个例子进行简单说明。

例1：02'''=+-y y y
解：特征方程为
012-2=+λλ
特征根为121=，λ（二重），故所求同解为
)(21x C C e y x +=.
例2：0134'''=++y y y
特征方程为
01342=++λλ
特征根为i 3-221±=，λ,故所求同解为
)3sin 3cos (212x C x C e y x +=-
2.2 二阶常系数非齐次线性微分方程的解法
现在讨论二阶常系数非齐次线性微分方程
)('''x f qy py y =++ （2.4）
的求解问题。

这里q p ,是常数，)(x f 是连续函数。

我们可以由其对应的齐次线性微分方程(2.1)的通解出发，使用常数变易法求出（2.4）的特解。

因而，只要能求出（2.1）的特征根，（2.4）的求解问题就已经解决。

但是，这样的方法往往是比较繁琐的，而且必须经过积分运算。

事实上，只要求得方程（2.1）的通解，再求出该方程的一个特解，就可得出它的通解表达式。

下面，我们讨论当)(x f 是某些特殊形式的连续函数时，所适用的求解其特解的简便方法——待定系数法。

2.2.1类型Ⅰ:x n e x P x f λ)()(=
设)(x P n 是n 次多项式，即
n n n n n p x p x p x p x P ++++=--1110)( （1≥n ）
（1）当λ不是特征根时，（2.4）有形如x n e x Q x y λ)()(=的特解，其中)(x Q n 是关于x 的n 次待定的多项式，即
n n n n n q x q x q x q x Q ++++=--1110)(
（2）当λ是)1(≥k 重特征根时，（2.4）有形如x n k e x Q x x y λ)()(=的特解，其中)(x Q n 也是形如上述的n 次多项式。

其中)(x y 中)(x Q n 的系数可以由待定系数法求得。

例3：x e y y y -=++432'''
解：对应齐次方程的特征方程为
01322=++λλ 特征根为2
1--121==λλ，，齐次方程的通解为 x x e C e C y 2121--+=
由于0，1都不是特征根，故已知方程有形如
x Be A y +=1
的特解。

将上式代入已知方程，得
6
1,4-==B A 即
x e y 6
141-= 因而，所求通解为
x x x e e C e C y 6
142121-++=--. 2）2.2.2类型Ⅱ：x n m e x x P x x P x f λωω)sin )(cos )(()(+=
其中，）
（、x P x P n m )(分别为两个已知的关于x 的m 次和n 次多项式，ωλ,为常数。

由欧拉公式，得
i
e e x e e x x
i x i x i x i 2sin ,2cos ωωωωωω---=+=. 故)(x f 可以改写成
i e e e x P e e e x P x f x i x i x n x i x i x
m 2)(2)()(ωωλωωλ---++= x i n x i m e x P e
x P )(~)(~)()(ωλωλ-++= （2.5） 其中，)(),(~~x P x P n m 分别是m 次和n 次多项式。

可以看出，（2.5）式就相当于两个类型Ⅰ形状的函数相加。

由非齐次方程的叠加原
理，就可求出类型Ⅱ的特解了。

设有二阶非齐次方程
)()('''21x f x f qy py y +=++
（2.6）
且)(),(21x y x y 分别是方程
),('''1x f qy py y =++)('''2x f qy py y =++ 的解，则函数)()(21x y x y +是方程（2.6）的解。

根据叠加原理及类型Ⅰ讨论的结果，我们有
1）当ωλi ±不是特征根时，（2.4）有如下形式的特解
x i l x i l
e x p e x p x y )()2*()()1*()()()(ωλωλ-++=
即 ]sin )(cos )([)()2()1(x x p x x p e x y l l x ωωλ+=
（2.7） 2）当ωλi ±是)1(≥k 重特征根时，（2.4）有如下形式的特解 ])()([)()()2*()()1*(x i l x i l
k e x p e x p x x y ωλωλ-++= （2.8）
即 ]sin )(cos )([)()2()1(x x p x x p e x x y l l x k ωωλ+=
（2.9） 其中)(),(),(),()
2()1()2*()1*(x P x P x P x P l l l l 为两个待定多项式，),max(n m l =. 注意：当）
（、x P x P n m )(中有一个恒为零时，方程（2.4）仍具有形如（2.8）、（2.9）的特解。

即不能当0)(≡x P m 时，就令0)()1(=x P l ，而0)(≡x P n 时，就令0)()
2(=x P l . 例4：x y y y 2cos 44'''=++
解：特征方程为
024422=+=++）（λλλ， 它有二重特征根-2=λ。

另一方面，方程的非齐次项为)(212cos 22x i x i e e x -+=.由此可见，相应的i 2±=μ与特征根-2=λ是不相等的。

因此，我们可设方程有特解
x b x a y 2sin 2cos 1+=
其中常数a 和b 待定。

把它代入原方程，得出 x x a x b 2cos 2sin 82cos 8=-,
由此推知
8
1,0==b a . 所以，原方程的通解为
x e x C C y x 2sin 8
1)(221++=-. 3 二阶微分方程的降阶和幂级数解法
3.1 可将阶的一些方程类型
1.方程不显含未知函数y 和未知函数的一阶导数'y ，即)(''x f y = （3.1） 若令p y =',那么'''p dx
dp y ==，则方程（3.1）即降为关于p 的一阶微分方程 )('x f p =，
两边积分得：1)(C dx x f p +=⎰，两边再次积分，就能得到方程（3.1）的通解.
2. 方程不显含未知函数y ,即
)',(''y x f y = （3.2）
若令p y ='，则方程（3.2）就变为
),('p x f p =，
这是一个关于p x ,的一阶微分方程.
例5：0)1)(1('2''=--+y x xy
解：将方程化为
0112'''=-+-x
x y y 令z y =',则上式化为
x x 11)]-[ln(z 2'
+-= 两边积分得
C x x z +-=-2
ln )1ln(2
因此
121'2
+=-x xe C y
再积分一次的通解
2212
C x e C y x ++-=-
3. 方程不显含自变量x ，即
)',(''y y f y = （3.3）
若令p y ='，那么
dy
dp p dx dy dy dp dx dp y =⋅=='' 则方程（3.3）就变为
),(p y f dy
dp p = 这是一个关于p y ,的一阶微分方程. 例6：求解0)('22'''=--y y y yy
解：这是不显含x 的二阶方程。

易见0=y 为一解。

若0≠y ,方程两边同除以2y 得 0)('22
'''=--y y
y y y 令dy
dp p dx dy dy dp dx dp y p y =⋅==
=''',,则方程可化为 y y p dy dp +=
即有1)1
('=y p
所以上式通解为
Cy y p +=2
从而得到
)(C y y dx dy
+=
将变量分离，两边积分得
1ln C Cx C y y
+=+
化简得原方程通解为
11
1C Cx C
Cx e Ce y ++-=
1,C C 为任意常数，这是特解0=y 包含于上述通解中。

4.恰当导数方程型
二阶微分方程也可以表示成0)'',',,(=y y y x F 的形式。

若方程
0)'',',,(=y y y x F
（3.4） 的左端恰为某一函数0)',,(=y y x G 对x 的全导数，即
)'',',,()',,(y y y x F y y x G dx d
=
则称方程（3.4）为恰当导数方程。

于是，方程（3.4）可写成
0)',,(=y y x G dx d
则有
C y y x G =)',,(，（C 为任意常数）
这样就把原方程降为了一阶微分方程。

例7：02)(2'''=++x y yy
解：这是一个不显含x 二阶方程。

将原方程写为
C x yy =+2'')( 积分一次得
1'C x yy +-= 即
1'
222
12C x y =+）（ 两边积分得
213
26
2C x C x y ++-= 因而，原方程的通解为
213
2
223
C x C x y ++-=
5.关于未知函数及其各阶导数都是齐次的方程
方程0)'',',,(=y y y x F 关于未知函数及其各阶导数都是齐次的是指
)'',',,(y y y x F 满足
)'',',,()'',',,(y y y x F t ty ty ty x F k ≡.
作变换⎰
=zdx
e y （z 是新未知函数），则有
⎰=zdx ze dx
dy ,⎰+=zdx e dx dz z dx y d )(222
代入到（3.4）中，有
0))(,,())(,,(22=⎰+=⎰+
⎰
zdx
k zdx zdx
e dx
dz z z x F e dx dz z ze x F 因为方程0)'',',,(=y y y x F 关于未知函数及其各阶导数都是齐次的，约去非零公因子
⎰
zdx
e ，得到
0))(,,(2=+
dx
dz
z z x F 上式经整理后可化为0)',,(=z z x f 的形式，这就是关于新未知函数z 的一阶微分方程。

注意:若0<y ，则可作变换⎰-=zdx
e y 。

实际问题中，我们作变换⎰
±=zdx
e y 后，还
要考虑0=y 是不是方程的解。

例8：求解方程
02)23()(3)(32'22'23'3''4=+++-+-y y x y x xy y x y x y x 解：这是左端关于''',,,y y y x 的三次齐次方程。

令
τue y e x x ==, 则
)(22
2
2τ
ττ
τ
d du d u d
e dx y d u d du dx dy +=+=- 代入原方程，消去公因子τ3e ,得到
0)(3
2
2=--τ
ττd du d du d u d (*) 令
p d du
=τ
(**) 则
du
dp
p d u d =2
2τ (***) 将（**），（***）代入（*），得到
0)1(
2=--p du
dp
p 因此
0=p 或
21p du
dp
+=
由
21p du
dp
+= 得
)tan(1C u p +=
即
)tan(1C u d du
+=τ
再积分可得
τe C C u 21)sin(=+ 即原方程的通解为
x C x C xar y 12)sin(-=
由0=p 得C u =,即Cx y =也是解，但此解包含在上述通解中。

3.2 二阶线性微分方程的幂级数解法
二阶线性微分方程
0)(')('')(210=++y x p y x p y x p
(3.5)
在近代物理学以及工程技术中有着广泛的应用，但是，当它的系数)(),(),(210x p x p x p 不为常数时，它的解往往不能用“有限形式”表示出来。

而幂级数解法就解决了这个问题，它不但对于求解方程有意义，而且由此引出了很多新的超越函数，在理论上具有很重要的地位。

定理 1 如果)(),(),(210x p x p x p 在某点0x 的邻域内解析，即它们可以展成)(0x x -的幂级数，且0)(00≠x p ，则（3.5）的解在0x 的邻域内也能展成)(0x x -的幂级数
.)(00n n n x x a y -=∑∞
=
（3.6）
定理2 如果)(),(),(210x p x p x p 在某点0x 的邻域内解析，而0x 是)(0x p 的s 重零点，是)(1x p 的不低于1-s 重的零点（若1>s ），是)(2x p 的不低于2-s 重的零点（若2>s ），
则方程（3.7）至少有一个形如
n n n
r
x x a
x x y )()
(00
0--=∑∞
=
（3.7）
的广义幂级数解，其中r 是某一常数。

3.3 二阶变系数线性微分方程的常系数化
对二阶变系数齐次线性微分方程
0)(')(''=++y x q y x p y
（3.8）
（其中)(),(x q x p 均为连续函数） 作变换z x f y )(=，则有
')()(''z x f z x f y +=，'')(')('2)(''''z x f z x f z x f y ++=
代入到（3.8）中，得
0))(')(''('))('2(''=+++++z f x q f x p f z f x p f fz
（3.9）
不妨令'z 的系数等于零，即
0)('2=+f x p f
从而
⎰
=-
dx x p e
f )(2
1
则
f x p x p f f x p f )]('2
1
)(41['',)(21'2-=-
= 代入到方程中，整理得
0)(''=+z x Q z （).('2
1
)(41)()(2x p x p x q x Q --
=） 当)(x Q 取某些特殊的函数时。

我们有：
1）2)(x
C
x Q =
（C 为常数），方程（3.9）可化为欧拉方程。

2）C x Q =)(（C 为常数），方程（3.9）可化为常系数线性方程。

4.拉普拉斯变换
我们已经知道二阶常系数线性方程
)('''x f qy py y =++ （4.1）
的通解结构和求解方法，但是，在实际问题中往往还要求（4.1）的满足初始条件
0000')(',)(y x y y x y ==的解。

我们当然可以先求出（4.1）的通解，然后由初始条件确定其中的任意常数。

此外，还有另外一种方法可以求解初值问题，即拉普拉斯（Laplace ）变换法.因为它无需求出已知方程的通解，而是直接求出它的特解来，从而在运算上得到很大简化。

1拉普拉斯变换的定义
设函数)(t f 在区间),0[+∞上有定义，如果含参量s 的无穷积分
⎰⎰
-+∞
∞→-=T
st T st
dt t f e dt t f e
)(lim )(
对s 的某一取值范围是收敛的，则称
⎰+∞
-=0)()(dt t f e s F st
（4.2）
为函数)(t f 的拉普拉斯变换，)(t f 称为原函数，)(s F 称为象函数，并且记为
)()]([s F t f =
2一些特殊函数的拉普拉斯变换 1）)0(Re 1
]1[>=s s
2）)0(Re 1
][2>=s s t
3）)0Re ,(!
][1>=+s n s n t n n 是正整数
4）)Re (Re 1
][a s a s e at >-=
5）)Re Re ,()
(!
][1
a s n a s n e t n at n >-=
+是正整数 6）)0(Re ][sin 2
2>+=s s t ωω
ω
7）)0(Re ][cos 2
2>+=
s s s
t ωω
3拉普拉斯变换的基本性质
1）线性性质：设函数)(1t f ，)(2t f 满足定理3的条件，则在它们的象函数共同的定义域上，有
)]([)]([)]()([22112211t f C t f C t f C t f C +=+
其中21,C C 为任意常数。

2）原函数的微分性质：如果)(),(''),(')(t f t f t f n 均满足定理3的条件，则
),0()]([)]('[f t f s t f -=
).0()0(')0()]([)]([)1(21-------=n n n n n f f s f s t f s t f
3）象函数的微分性质：如果)()]([s F t f = ，则
)],([)(t tf s F ds d -= )].([)1()(t f t s F ds
d n n n n
-= 4）如果)]([)(t f s F =，则).()]([a s F t f e at -=
5 二阶微分方程的存在唯一性
5.1 存在唯一性定理
如果在二阶微分方程
)',,(''y y x f y =
（5.1）
中，令1'y y =，则'''1y y =，它就可化为方程组
⎪⎪⎩⎪⎪⎨
⎧==)
,,(111
y y x f dx
dy y dx
dy
（5.2）
我们称（5.2）为一阶微分方程组。

从而，要讨论二阶微分方程的初值问题的存在唯一性，就只需讨论一阶微分方程组的初值问题的存在唯一性。

令⎪⎪⎭⎫
⎝⎛=)()()(1x y x y x Y ，⎪⎪⎭⎫ ⎝⎛=),,(),(11y y x f y Y x F ，
并定义：
⎪⎪⎪
⎪⎭
⎫ ⎝⎛=dx dy dx dy dx x dY 1)(，⎰⎰⎰⎪
⎪⎪⎭⎫ ⎝⎛=x x x x x x dx y y x f dx y dx Y x F 000),,(),(11， 则（5.2）可记成向量形式
),(Y x F dx
dY
= （5.3）
初始条件100100)(,)(y x y y x y ==可记为
00)(Y x Y =，其中⎪⎪⎭
⎫
⎝⎛=1000y y Y
则二阶微分方程
⎩⎨
⎧===0
000')(',)()
',,(''y x y y x y y y x f y （5.4）
的初值问题就可记为
⎪⎩⎪
⎨⎧==0
0)(),(Y x Y Y x F dx
dY
（5.5）
此外，我们把二维向量⎪⎪⎭
⎫
⎝⎛=)()()(1x y x y x Y 的范数||||Y 定义为
||||||||1y y Y +=.
下面，我们给出初值问题（5.5）的解的存在与唯一性定理。

定理3 如果函数),(Y x F 在三维空间的区域
b Y Y a x x R ≤-≤-||||,|:|00
上满足：
1）连续；
2）关于Y 满足李普希兹()Lipschitz
条件，即存在0>L ，使对于R 上任意两点),(),,(21Y x Y x ，有
||||||),(),(||2121Y Y L Y x F Y x F -≤-，
则初值问题（5.5）的解在区间[]h x h x I +-=00,上存在唯一，其中
||),(||max ),,
min(),(Y x F M M
b
a h R Y x ∈==.
类似于一阶微分方程的初值问题的存在唯一性定理的证明，下面来简单证明一下定理3.
引理：如果函数),(Y x F 在三维空间的区域
b Y Y a x x R ≤-≤-||||,|:|00
上连续，则初值问题（5.5）的解)()Y(x x Φ=,[]h x h x I x +-=∈00,,与积分方程
⎰+=x
x dx x Y x F Y x Y 0
))(,()(0
（5.6）
在区间[]h x h x I +-=00,上的连续解等价，其中),
min(M
b
a h =，||),(||max ),(Y x F M R Y x ∈=.
由引理我们知道，要证明定理3，只要证明积分方程（5.6）的连续解在区间
[]h x h x I +-=00,上存在唯一就行了。

1 存在性的证明
下面用皮卡()Picard 逐次逼近法来证明积分方程(6)的连续解的存在性，可分三个步骤进行。

（1）构造区间I 上的逐次近似的连续向量函数列{})(x Y n . 令00)(Y x Y =，构造毕卡逐次逼近向量函数序列如下：
⎪⎩⎪
⎨
⎧=+==-⎰),2,1())(,()()(100
00 n d Y F Y x Y Y x Y n x x n
ξξξ 向量函数 )(x Y n 称为 （5.5）的第n 次近似解。

用数学归纳法可以证明：
b Mh x x M d M d Y F Y x Y x
x n x x n ≤≤-=≤≤-⎰⎰-|||||||))(,(|||||)(||0100
ξξξξ
即曲线)(x Y Y n =未越出区域R ，保证了逐次逼近可以一直进行下去。

（2）证明函数序列{})(x Y n 在区间I 上一致收敛。

考虑向量函数项级数
+-++-+-+-)]()([)]()([)]()([)(112010x Y x Y x Y x Y x Y x Y x Y n n （5.7） 它的部分和是
)()]()([)]()([)]()([)()(1120101x Y x Y x Y x Y x Y x Y x Y x Y x S n n n n =+-++-+-+=-+ 所以，要说明函数序列{})(x Y n 在区间I 上一致收敛，只需证明级数（5.7）在区间I 上一致收敛。

ξξξd Y F Y x Y x
x ))(,()(0010⎰=-
∴|||||))(,(|||||)(||00010
x x M d Y F Y x Y x
x -≤≤-⎰ξξξ
由数学归纳法，我们可以得到：
)!
1(||||)()(||1
01+-=-++n x x M L x Y x Y n n
n n
而h x x ≤-||0,易于看出级数（5.7）每一项的绝对值都不会超过正项级数
的对应项。

上面的级数显然是收敛的。

从而，级数（5.7）在区间I 上一致收敛。

设其和函数为)(x Φ，从而函数序列{})(x Y n 在区间I 上一致收敛于)(x Φ。

由于)(x Y n 在区间I 上是连续的，因而)(x Φ也是连续的。

（3）证明)(lim )(x Y x n n ∞
→=Φ是积分方程（5.6）的解。

对ξξξd Y F Y x Y n x
x n ))(,()(100
-⎰+=两边取极限，得
⎰-∞→+=x
x n n d Y F Y x 0
))(,(lim )(10ξξξΦ
要证)(lim )(x Y x n n ∞
→=Φ是积分方程（5.6）的解，只需证
⎰⎰=-∞→x
x x x n n d F d Y F 0
))(,())(,(lim 1ξξξξξξΦ
+++++-!
!2||1
20n h M L h LM Mh Y n n
{})(x Y n 在区间I 上一致收敛，
,0,0>∃>∀∴N ε使N n >∀时，有Lh
x x Y n ε
<
-||)()(||Φ.
|||))(,())(,(|||||))(,())(,(||0
11ξξξξξξξξξξξd F Y F d F d Y F x
x n x
x x
x n ⎰⎰⎰-≤---ΦΦ
εε
ξε
ξξξ≤-=≤-≤⎰
⎰-|||||||)()(|||010
0x x Lh
L
d Lh
L d Y L x
x x
x n Φ
⎰⎰=∴-∞→x x x
x n n d F d Y F 0
))(,())(,(lim 1ξξξξξξΦ
则)(lim )(x Y x n n ∞
→=Φ是积分方程（5.6）的解。

2 唯一性的证明
设)(~
x Φ也是积分方程（6）的解，且满足.00~
)(Y x =Φ则有
,)(,()(~
0~
ξξξd F Y x x
x ΦΦ⎰+=
于是
ξξξξξξξξd L d F F x x x
x x
x ||)()(|||||))](,())(,([||||)()(||0
~
~
~
⎰⎰-≤-=-ΦΦΦΦΦΦ
.||,||)()(|||0~
h x x d L O x
x ≤--+=⎰ξξξΦΦ
由Bellman 不等式得：
.||,0||)()(||0~
h x x x x ≤-≡-ΦΦ
得出矛盾。

因此，(5.6)在h x x ≤-||0的解唯一。

综上，（5.5）的存在唯一性定理得证。

结论
关于二阶线性微分方程的研究已经取得了不少成就，尤其在二阶常系数线性微分方程的求解问题和解的存在唯一性定理等方面卓有成效。

二阶微分方程的解的存在唯一性定理不仅可判断解的存在唯一性，而且还有着广泛的应用。

而幂级数解法作为求解二阶变系数齐次线性微分方程的一种方法，其过程还是比较繁琐的，计算量偏大，且需要考虑函数是否解析，幂级数在某个区间是否收敛等。

另外，对于二阶变系数非齐次微分方
真诚为您提供优质参考资料，若有不当之处，请指正。

程，目前还尚有通用的求解方法，只有一些特殊类型是可以求解的，还有待于进一步的发展和研究。

参考文献
[1] 吴美捷，浅谈待定系数法[A]. 当代经理人，2006,(03).
[2]朱思铭，王寿松，王高雄等. 常微分方程[M].北京：高等教育出版社, 2006 .
[3] 吴志坚,吴筱坚, 傅里叶变换与拉普拉斯变换[J], 石油大学学报(自然科学版),1996,(10)
[4]丁同仁,李承治,常微分方程教程[M]．北京：高等教育出版社，2004.
[5]朱乃明，李虹莉.常微分方程[M]．重庆:西南师范大学出版社,2005.
[6] 陈兆华,费仁允, 欧拉公式的证明与应用[J].数学通报,2005.
[7]都长清，焦宝聪，焦炳照.常微分方程[M]．北京：首都师范大学出版社，2001.
[8]黄启昌.常微分方程[M]．北京：高等教育出版社，1982.
[9]王克，潘家齐.常微分方程学习指导书[M]．北京：高等教育出版社，2007.
[10]田巍，李奇.二阶常系数非齐次线性微分方程特解的特征根公式法[J].高师理科学刊,2007,(06).
[11]刘培进.二阶常系数线性非齐次微分方程的公式解法[J].山东师范大学学报(自然科学版)2002,(03).
[12]孙梅娟, 倪致祥.可线性常系数化的二阶常微分方程[J].阜阳师范学院学报(自然科学版),2010,(02).
[13]石正华.浅谈二阶变系数齐次微分方程的求解问题[J].南昌教育学院学报,2012,(01).
[14]李录苹,王通.关于几类二阶微分方程的解法[J].雁北师范学报,2006,(02).
[15]李永利, 桑改莲. 一类二阶变系数齐次微分方程通解的求法[J].高等数学研究, 2006,(03).
[16]胡劲松等.一种二阶变系数线性微分方程的求解方法[J].重庆工商大学学报(自然科学版),2005,(03).
[17]陈湘.对常系数非齐次线性微分方程的一种讲授法[J].高等数学研究，2010,(03)
[18]李中平.用观察法求二阶变系数齐线性方程的非零特解[J].高等数学研究,2010,(03).
[19]Slimane Benaicha. Periodic Boundary Value Problem for Second Order Ordinary Differential Equations. Applied Mathematical Sciences, Vol. 3, 2009, no. 6, 267 – 276.
[20]L. F. Shampine. Numerical Solution of Ordinary Equations. Chapman and Hall, 1994.
[21]R.P. Agarwal. Boundary value problems for high order differential equations, World Scientic, Singapore, 1986.
21 / 21。