复根

南京师范大学泰州学院

毕业论文（设计）

（一六届）

题目：二阶常微分方程的解法

院（系、部）：数学科学与应用学院

专业：数学与应用数学

姓名：潘陆

学号08120146

指导教师：刘陆军

南京师范大学泰州学院教务处制

南京师范大学泰州学院本科毕业论文

1

摘要：本文主要是介绍了二阶常微分方程众多解法中的三种，分别为特征方程法，拉普

拉斯变换法和常数变易法，研究并讨论了二阶常微分方程在特征方程法中特征方程根为

实根，复根和重根的情形。我们选用了弹簧振子系统的振子运动，用这三种不同的方法

来解决该问题。

关键词：二阶常微分方程；特征根法；常数变易法；拉普拉斯变换

Abstract:Themainpurpoofthispaperisthecond-orderordinarymany

differentialequationsolutionofthree,respectivelyasthecharacteristicequation

method,Laplacetransformmethodandvariationofconstantsmethod,studyand

discussthecond-orderoftendifferentialequationinthecharacteristicequation

oftherootsofthecharacteristicequationforrealroots,complexrootsandroot

thespringoscillatortheoscillatormotion,thethree

differentmethodstosolvetheproblem.

Keywords:condorderordinarydifferentialequation;Characteristicanalysis;

constantvariationmethod;Laplastransform

南京师范大学泰州学院本科毕业论文

2

目录

1绪论.........................................................3

1.1二阶常微分方程的起源和发展史......................................3

1.2二阶常微分方程的介绍..............................................3

1.3研究二阶常微分方程的目的与意义....................................4

2二阶常系数常微分方程的几种解法...............................5

2.1特征方程法........................................................5

2.1.1特征根是两个实根的情形......................................5

2.1.2特征根有重根的情形..........................................6

2.2常数变易法........................................................7

2.3拉普拉斯变换法....................................................9

3二阶常微分方程解法的应用（分析例题）........................11

3.1特征方程法.......................................................11

3.2常数变易法.......................................................13

3.3拉普拉斯变换法...................................................14

4结论和启示..................................................16

谢辞.........................................................18

参考文献......................................................19

南京师范大学泰州学院本科毕业论文

3

1绪论

1.1二阶常微分方程的起源和发展史

既然说到了微分方程，就不能不提到海王星的故事，它的发现是人类智慧的硕果，

微分方程在其中扮演了重要的角色，并且在其中也包含数学演绎法的作用。在发现了天

王星之后，进行天文观测的人们发现它所处的位置总是和万有引力计算出的位置有些许

不同，于是有人质疑万有引力定律的正确性。但也有一部分人认为，这也可能是天王星

在受到一颗尚未发现的行星的吸引力才会造成的改变。不少人坚信这种假设是正确的，

但却很少有人找到了正确的方法并付诸实践。而英国的一个学生亚当斯显然不是其中之

一，他勇敢接受了这项任务，运用手头仅有的资料建立起微分方程，成功求出了海王星

的位置与下次出现的时间。1843年10月21日，满怀信心的亚当斯把结果寄给了天文台

长艾利，但换来的却是质疑，艾利并不相信籍籍无名的他。然而在两年后，另一名青年

勒威耶也计算出了同样的数据，并把计算的结果给予了位于柏林天文台的助理员卡勒，

在那个值得铭记的夜晚，卡勒在计算出的位置上发现了第七颗行星海王星。

从二十世纪三十年代以来，常微分方程的研究像是走上了快车道，迅速发展并建立

起了多个分支。

1927－1945年期间定性理论的主要研究是与无线电技术紧密联系在一起的。在第二

次世界大战期间由于对通讯等方面的需求越来越高，极大地激发了对无线电技术的研究

进展，尤其是对非线性振动理论的研究取得了迅速的发展。

在四十年代之后各国大部分数学家们主要在研究对抽象动力系统的拓扑特征,例

如闭轨的存在性、结构的稳定性等,对于二维系统来说，我们可以通过一些方法证明他

的结构稳定性；而对于一般的系统来说这个问题依旧困扰着我们。在动力系统的研究方

面,目前采用的办法是从典范方程组到阻碍集有详尽的理论指导,成功解决了一系列困

扰人类多年的问题,其中最为突出的是C’封闭引理的证明,以及对结构稳定性的充要

条件等方面都作出了杰出贡献。

在当今社会，由于信息技术的飞速发展，大量的领域需要用到常微分方程组进行描

述。前赴后继的杰出的学者们，为了各种稳定性及专业问题，终其一生都在研究常微分

方程，也取得了不朽的成就，但依然有很大的疑问等着我们去解开。

1.2二阶常微分方程的介绍

二阶微分方程在时间上大致与微积分同时产生。对于初学者来说，)(xfy



这样

的问题就是最简单的微分方程了。二阶常系数线性微分方程是形如)(xfqyypy







的微分方程。与其对应的二阶常系数齐次线性微分方程为0







qyypy，其中qp,是

实常数。若函数

1

y和

2

y之比为常数，称

1

y和

2

y是线性相关的；若函数

1

y和

2

y之比不为

常数，称

1

y和

2

y是线性无关的。

南京师范大学泰州学院本科毕业论文

4

1.3研究二阶常微分方程的目的和意义

就如上文所说，研究二阶常微分方程已经取得了不少成就，尤其是在对方程的求解

方面。与此同时应用常微分方程的理论也经过大家的不懈努力，取得了丰硕的成果。但

是对于进一步的发展所需，还是小巫见大巫，所以我们要更加努力的钻研，努力地完善

这门学科的理论系统。在数学的发展历史中，数学分析占有非常重要的地位，我们大学

学习的课程也都是以数学分析作为基础，而微分方程正是数学分析的关键所在。同时它

也发展出了数学分析中大部分思想以及理论。众所周知的，常微分方程自始至终都是人

类用来探索自然变换，研究自身社会结构，工程问题以及大自然的生态结构的便利的道

具。常微分方程由于与现实生活息息相关，所以对其的研究一直没有停止过，而且表现

出欣欣向荣的活力。并且在多个学术领域中，常微分方程都占着决定性的作用，可以说

常微分方程带领着人类的进步。而二阶常微分方程同样在常微分方程的整套理论中有着

弥足轻重的地位，在各个研究领域中都有十分广泛的应用。

南京师范大学泰州学院本科毕业论文

5

2二阶常系数常微分方程的几种解法

2.1特征方程法

特征方程法中的特征方程，是为了对相对应的数学对象进行深入的研究而人为引入

的一些等式，当研究的对象改变时，它也会改变，这其中包括数列特征方程，微分方程

特征方程和积分方程特征方程等等。

求微分方程0

2

2

qx

dt

dx

p

dt

xd

的通解。

解：特征方程02qp的根

21

，,

(1)若

21

，

是两个不相等的实根,那么上面这个微分方程就拥有两个实值解

ttee21,,于是我们可以求得其通解为

ttececx21

21

（

21

,cc为常数）.

(2)若

21

，

相等,则该方程有二重根,因此方程的通解具有形状如

ttececx11

21

（

21

,cc为常数）.

(3)若

21

，为共轭复根biaz的情况,则该方程的通解具有形状

)cossin(

21

btcbtcexat（

21

,cc为常数）.

在数学中，许多公式与定理都需要进行证明,下面本文给出前两个解答的理论依据

及证明过程。

2.1.1特征根是两个实根的情形

设

21

，为该特征方程的两个不等实根,那么我们可以得出与之对应的方程的两个

解为ttee21,,

我们确定这两个解在bta上线性无关,所以它们能够组成该方程的基本解组。事

实上,这时

21

)(

21

11

)(21

21

21











t

tt

tt

e

ee

ee

tw,

而上面最后一个行列式是著名的范德蒙德（Vandermonde）行列式,该行列式与

）（

12

-相等。由于在之前假设了

12



,所以此行列式不等于零,从而0)(tw,于是得

南京师范大学泰州学院本科毕业论文

6

到ttee21,线性无关,这就是要证明的结论。而该方程的通解即可表示为

ttececx21

21



（其中

21

,cc为任意数）.

2.1.2特征根为复根时

讨论这个特征方程含有复根的情况,首先这个方程的系数为实数,同时也是常数，所

以复根成对共轭出现。设i

1

是特征根之一,则i-

2



就是第二个特征根,根据

这两个特征根我们就可以得到原方程的两个不同的实值解:

)sin(cos)(titeeiti,

)sin(cos)(titeeaiti.

根据定理,我们得到的复值解的实部和虚部与方程的解相同。那么由于这一对共轭

复根i对应于特征方程,可以求得方程0

2

2

qx

dt

dx

p

dt

xd

的两个实值解

teteaiaisin,cos.

2.1.3特征根有重根的情形

设特征方程有k重根

1

则易得

0)(,0)()()(

1

)(

1

)1(

11

kkFFFF,

先设

0

1

,表示特征方程有因子k

,于是

0

11



knnn

aaa,

也就是该特征方程的形状为

01

1





k

kn

nnaa,

而与其对应的方程0

1

1

1

1









xaa

dt

xd

a

dt

xd

xL

nn

n

n

n

n



则变为

0

1

1

1









k

k

kn

n

n

n

n

dx

yd

a

dx

yd

a

dx

yd



.

易得它有k个解,12,,,kttt

,而它们是线性无关的。于是我们可以得出,特征方程中

南京师范大学泰州学院本科毕业论文

7

的k重零根对应方程的k个彼此线性无关的解,12,,,kttt

。当这个k重根0

,我们作

变量变换tyex1



,可注意到

























yy

mm

ymyeyexm

mmm

t

m

t

m

1

2

2

1

1

1!2

1

11

，

可得

tt

n

n

n

n

n

teyLeyb

dt

yd

b

dt

yd

yeL111

1

1

1

1

























,

于是可将对应方程化为

0

1

1

1

1

1

11











yb

dt

yd

b

dt

yd

b

dt

yd

yL

n

n

n

n

n

n

n

,

而其中

n

bbbb,,,,

321

仍为常数,其相应的特征方程为

0

1

1

1







nn

nnbbbG,

直接进行计算易得

tt

t

tteGeeLeLeF1111

11





,

因此

GF

1

,

从而

kjGFj,,2,1,)(

1

.

通过这样转化,问题就化为前面所讨论过的情况了。

2.2常数变易法

接下来我们要说的是求解微分方程的一种极其重要的方法,常数变易法。我们在求

解一阶线性微分方程时常常采用这种方法。这个办法的大体思想，是通过将常数N代入

XU

中来得到该方程的通解。看似简单，但拉格朗日为之研究奋斗了11年，而我们所

如今所采用的仅仅是他所得出的结论，证明过程省略了。它在非齐次线性微分方程和与

其对应的齐次线性微分方程之间起着重要的链接作用。

我们这里讨论的对一般二阶常微分方程的求解,要先得出该方程的一个特解,再用

我们上面讲解的特征方程法求通解。

例题求常微分方程tfqx

dt

dx

p

dt

xd



2

2

的通解。

南京师范大学泰州学院本科毕业论文

8

解方程tfqx

dt

dx

p

dt

xd



2

2

与其对应的齐次方程为

0

2

2

qx

dt

dx

p

dt

xd

,

其特征方程为

02qp.

由于tfqx

dt

dx

p

dt

xd



2

2

的通解根据理论可以得到就等于这个方程所对应的齐次

线性微分方程的通解与求出的它本身的一个特解之和，而我们已经讨论了二阶常系数齐

次线性常微分方程的通解,所以只要再求出一个该方程本身的特解。

若



为该方程的实根,则tex为方程0

2

2

qx

dt

dx

p

dt

xd

的解。根据常数变易法

tfqx

dt

dx

p

dt

xd



2

2

的其中一个解为tetcx,代入原方程并进行化简得

tfetcptct



2,

这是即是关于tc的一阶线性微分方程,而其一个特解为

dtdttfeetctptp2)(,

从而得到上面方程的一个特解为

dtdttfeeextptpt2.

讨论



为该方程的复根的情况,可以设0,,bRbabia，则btexatsin就是

原方程xfqx

dt

dx

p

dt

xd



2

2

的解，接下来用常数变易法设其中的一个特解为

btetcxatsin，通过运用情形1解得方程tfqx

dt

dx

p

dt

xd



2

2

的一个特解为



dt

bt

btdtetfe

btex

tapap

at





2

22

sin

sin

sin

.

由于*x

是特解,所以积分常量可以都取零。

南京师范大学泰州学院本科毕业论文

9

2.3拉普拉斯变换法

接下来要介绍的拉普拉斯变换法是一种积分变换法，又名为拉氏转换法。它是一个

线性变换法，通过因数为实数0tt

，的函数中的实数转化为复数i。在某些情形下一

个实变量函数在R中进行运算比较困难，若能将实变量函数进行拉普拉斯转换，并在复

数域中进行各种运算，再使用拉普拉斯反变换来求得该方程或函数在实数域中的相应结

果，是一种简便的计算方法。拉普拉斯转换的运算步骤对于我们要求的解线性微分方程

极其有效，它可以把微分方程转化为较易求解的代数方程来进行处理，使计算大幅度简

化。在经典控制理论中，对于控制系统的分析但和综合等，都是立足在拉普拉斯转换的

基础上的。我们引入拉普拉斯变换有一个主要优点，就是可以采用传递函数来代替常系

数微分方程用以描述系统的特性。当我们要用图解来测定某些控制系统的整体特性并且

分析该控制系统的运动过程等情况时，可以用上这种方法，一般来说比较简单。

由积分dttfesFst





0

.

所定义的确定在复平面Re

上的复变数s的函数sF,称其为函数tf的拉普

拉斯转换,我们称tf为原函数,而将sF称为像函数。

拉普拉斯变换法主要是通过借助拉普拉斯变换将常微分方程转化，使它成为复平面

s

的代数方程。并通过一些简单的代数运算,再利用拉普拉斯变换表,即可对微分方程进

行求解。方法十分简单方便,在现实生活中也经常被采用。但是这个方法也存在着些许

不足，主要是只适用于右端函数为原函数的情况。

求解方程011,2

2

2

xxex

dt

dx

dt

xd

t.

解首先使1t

，将问题转化为

000,21

2

2

xxex

dt

dx

dt

xd

t,

再对新方程两边进行拉普拉斯变换,得到



es

sXssXsXs

1

1

1

22



,

因此



e

s

sX

1

1

1

3







,

查阅拉普拉斯变换表可得

12

2

1

ex,

从而得到

tettx21

2

1

,

南京师范大学泰州学院本科毕业论文

10

这即为所要求的解。

当然,求解二阶或更高阶的常微分方程我们能用的方法还有不少,这里由于篇幅我

不能将其一一列出。但我们利用前面的一些结论就足以解决下面的几个物理问题了。

南京师范大学泰州学院本科毕业论文

11

3二阶常微分方程解法的应用

上面介绍了解二阶常微分方程常用的三种方法，为了更好的体现这三种方法的区

别，我选取了一个比较简单的动力学方程来进行求解。通常，在对物理问题求解时，我

们有三个步骤：第一步是对该问题进行彻底分析从而能做到对方程的初步建立并且对定

解条件进行明确；第二步是对其解的性质进行探究或者求出方程来满足初始条件的特

解；最后一步是定性分析对解，对原来的问题反过来进行解释，其中最为关键的步骤就

是要将方程列出，而能列出方程的方法主要有两种，分别为微元分析法和瞬时变换法。

但是在研究阻尼运动的过程中，求解运动方程一直是令人头疼的问题。接下来我们就分

别使用上述三个方法来求解这个动力学方程。

3.1特征方程法

在我们要研究的弹簧振子系统中，测定物体的阻尼系数10.10s，物体的质量为

kgm0.1,该弹簧所具备的劲度系数为175mNk,根据上述条件，假设质点由静止状

态逐步开始运动，对弹簧振子的位移方程进行求解。

解：根据牛顿的第二运动定律我们可以得出

macvkx1

或

0

2

2

kx

dt

dx

c

dt

xd

m,2

由于相对来说振动系统这是在这之前给定的，其中包含的常量为ckm,,，如果能够

确定2/,/2

0

mcmk,那么上面的方程式就可以转变为：

022

0

2









dt

d

dt

d

,3

那么根据所得到的数据代入公式（3）我们就可以得到

07520

2

2

x

dt

dx

dt

xd

.4

通过对以上公式的仔细观察和深入研究则可以得到对该方程进行求解能够使用特

征值法，得出其特征方程可以表述为：075202

,并且在该特征方程当中包含有

不同的两个根

5,15

21

,这样相对应的公式4的两个根分别为

ttee15

2

5

1

,5

南京师范大学泰州学院本科毕业论文

12

那么根据公式5进行计算可以即可得到固有角频率数值为25

0

mk，此时

阻尼系数值为10

,也就是说2

0

2,于是方程5的解可以表述为

ttBeAe155（初始条件可得BA,数值）.6

通过公式6，我们可以看出振子所保持的属于一个非振动状态，在这样的背景下，

我们要求的质点也只是在原先的不平衡位置逐渐恢复到平衡状态当中，该质点并不具有

周期振动的特征。由于我们的关注点是在

0



情况下面,质点会呈现出逐渐衰减的振

动。然而由于会受到阻尼作用的影响，他不能保持自由振动系统的长久运作，振动会逐

渐衰减直至振动停止，如果我们要保持振动持续不停的状态的话，该质子就需要从外界

获得运动必要的能量，在学术界将这种由于受到外部的持续作用而产生振动的情况称为

强迫振动。

我们再举一个例子：假如在以上的振动系统当中振子由于受到某个外力

NtF30cos100的作用，在公式之中100

A

F表示的是该驱动力的幅度值，30

表

示的是该驱动力所具有的圆频率，f即是该驱动力所保持的频率。

解：由于在质点振动系统当中会受到驱动力的作用，那么我们就可以得到关于系统

振动的方程：

Fkx

dt

dx

c

dt

xd

m

2

2

,7

上述方程还可以表示为

tHx

dt

dx

dt

xd

30cos22

0

2

2

.8

m

F

HA表示的是在单位质量上该振子所受到的外力大小。而7和8这两个式子都

是强迫运动方程。我们观察该式子的本质发现，他们属于二阶非齐次常系数微分方程，

而根据定理我们对这个方程求解后所能得到的一般解同时也是该方程所能得到的一个

特解和与该非齐次方程相对应的齐次方程的一般解之和。我们在前面求出了其相对应的

自由振动方程的一般解，那么现在要做的就是找到8方程的一个特解，将其代入方程8

就能得到方程9：

tx

dt

dx

dt

xd

30cos1007520

2

2



,9

在这里可以假设方程9有着

tBtAx30cos30sin

1



这样的特解，将这个特别对方

南京师范大学泰州学院本科毕业论文

13

程9替代进去并且将其简化最终得到

ttBAtBA30cos430cos332430sin2433

,

通过比较同类项系数我们可以得到

555

44

555

32

BA，,这样我们就可以进一步得

出ttx30cos

555

44

30sin

555

32

1

,而根据前面所得到的结果，可以将原方程的通解表示出

来，为：

ttBeAetxtt30cos

555

44

30sin

555

32

)(155.

综上所述，题目中给我们的条件决定BA,的数值，之前的两项为该方程的瞬态解，

瞬态项对这整个系统所进行的自由衰减振动能够进行有效的描述，然而只能在运动的开

始阶段起作用，在经过长时间的运动之后，它起到的影响会随时间消逝并且在运动最后

完全消失。而之后两项所代表的稳态解，描述的是强迫运动的状态，由于幅值条件的固

定，所以我们称这样的状态为稳定状态。根据上面的公式，外力作用在质点上的时候，

整个系统的振动状态十分复杂，这时候的振动既包括了瞬态振动，也包括了稳定振动，

而这样的振动状态对于在强迫振动之中逐步建立起稳态振动的过程进行有效的描述。在

长时间的振动之后，瞬态振动终将消逝，这整个系统就会保持稳态振动的状态。

3.2常数变易法

从上面的分析之中我们可以了解到tex5即是属于特征方程075202

的实

根，于是我们就可以得到tex5为属于上面方程9当中的一个根，之后通过运用常数

变易法设置tetcx5，在这一过程中我们也可以得到方程其中一个解为x，将数值代

入到方程9当中并且对其进行简化可以得到

tetctct30cos100105



.

属于tc的一阶线性微分方程,并且得到在该方程当中一个特解为



1

5530cos

3

4

30sin

3

8

ctetetctt,

从而得出方程9的一个特解为（取

0

21

cc

）





























21

55530cos

3

4

30sin

3

8

cdtteteetxttt

南京师范大学泰州学院本科毕业论文

14

tt30cos

555

44

30sin

555

32

,

于是可得方程9的通解为

ttBeAetxtt30cos

555

44

30sin

555

32

155.

由上面的结论可知

Fkx

dt

dx

c

dt

xd

m

2

2

.10

将题目数据代入公式中可以得到

tx

dt

dx

dt

xd

2cos40020

2

2

.11

通过观察我们可以发现，在进行该类求解的过程中要使用常数变易法，首先就是要

求出公式11，而在之前的探索当中已经可以得到公式11的特征方程为

0400202.于是可以进一步的假设该特征方程的根为

i310，那么

tetxt310sin10

即为公式11的一个解。运用常数变易法可设为

tetctxt310sin10

.

这里套用情形1的解法，将

*()xt

代入12并进行化简得

tttx2cos

39604

99

2sin

39604

10

.

由于x是特解,所以积分常量可以为零。

3.3拉普拉斯变换法

依然使用上面的例子，通过牛顿第二运动定律可以得到

Fkx

dt

dx

c

dt

xd

m

2

2

,

将公式代入数据之后我们可以得到

tx

dt

dx

dt

xd

2cos40020

2

2



,12

根据条件，质点由静止状态开始运动，我们可以得到以下的公式

南京师范大学泰州学院本科毕业论文

15

0,0

0



dt

dx

x

t

,

对方程(12)进行拉普拉斯变换得到



4

40020

2

2





s

s

sXssXsXs,

即

40020

1

4

)(

22







sss

s

sX,

将上式右端分解为部分分式



439604

99

4

2

39604

10

22







s

s

s

sX

2

2

2

231010

10

39604

99

31010

310

118812

3101











s

s

s

.

通过拉普拉斯变换表可得

tttx2cos

39604

99

2sin

39604

10



tetett310cos

39604

99

310sin

118812

3101

1010

.

南京师范大学泰州学院本科毕业论文

16

4结论和启示

在本篇文章中，我们首先介绍了海王星的发现过程，并介绍了研究二阶常微分方程

的背景及研究现状，之后重点总结了解二阶常微分方程最具有代表性的三种不同的解

法，特征值法，常数变易法和拉普拉斯转换法，并且进行了一定程度的探究。

在特征方程法中，我们通过研究特征方程的特征根的情况，来判断原方程的根的情

况。若

21

，

相等,即该方程有二重根，这是第一种解的情况。而

21

，

为共轭复根

biaz时，为第二种情况，并给出了这两种情况的证明过程。

常数变易法则是使用了常数N代入xU

的方法，并且提出了定理，原方程的通解

就等于这个方程相对应的齐次线性微分方程的通解，同时也运用了特征方程法的思想。

最后我们介绍了比较简便的拉普拉斯变换法，它是一个线性变换法，通过因数为实

数0tt，的函数中的实数转化为复数i来进行运算。对某些情况下的实变量函数进行

拉普拉斯转换，使得我们能在复数域中对其进行各种运算，最后再使用拉普拉斯反变换，

得出结果，这是一种简便的计算方法。该运算步骤对我们这篇文章要求解的线性微分方

程是相对来说最好的办法，它可以把微分方程转化为较易求解的代数方程来进行处理，

使我们不用再进行繁琐的计算。

之后我们用这三种方法求解了弹簧振子的运动方程，以此来比较三种方法的相同点

与不同点。并且通过这个方程的解答，体现出了二阶常微分方程与我们的生活息息相关

的特性。

由于篇幅及本人专业知识的限制，无法进行更深层次的探究，比如求解二阶常微分

方程还有一种幂级数解法，但是它的解题过程非常繁琐，对计算的要求很高，计算量比

较大，还要考虑该函数是否解析以及幂级数在某个区间是否收敛等问题，所以在此不做

讨论。

二阶常微分方程研究的道路远不止此，在这个时代，二阶常微分方程被广泛应用于

网络，军事，医学等各种高科技领域，对二阶常微分方程的研究，是对目前社会发展的

推进。

类似于本篇中写出的弹簧振子方程，研究二阶常微分方程不仅能够推进自身学科的

进步，还能提升别的学科的发展速度，它的各个分支学科也与人类生活息息相关。

总而言之，在高速发展的现代，常微分方程的研究工作虽然仍有很长的路要走，但

由于其重要性不言而喻，所以我们要全身心的投入到研究中去。对于二阶常微分方程，

目前并没有研究出通用的求解办法，只能对少数特殊情况的方程进行求解。这无论是对

现实生活还是研究工作都是非常不利的，所以对常微分方程的求解问题，任重而道远。

我们唯有脚踏实地，艰苦奋斗，为了人类更好的明天而努力。

而通过本次论文的写作，我也明白了有志者事竟成这个道理，无论是开头讲的海王

星的发现，还是在写作过程中读到的各个伟人的事迹，都告诉我只要努力了，终会得到

回报，或许得到的回报有大有小，甚至为目标努力奋斗了许久却仍未得到回报，这时更

不能放弃，继续做下去还有希望，而如果放弃了得到的就一定是失败。然而，不相信自

己的人，没有努力的价值，所以，相信自己，勇敢拼搏，最终定能获得成功。

在对这篇作文的写作过程中，对物理学方程的研究遇到了一些麻烦，因为我的物理

南京师范大学泰州学院本科毕业论文

17

学知识比较薄弱。但是通过同学、老师的帮助，再在网络及学校图书馆翻阅了各种文献

之后，我还是成功的做成了。这让我想到了，在今后的教育工作中，我应该保持这种勇

于探索的精神，不能以非专业，没学过为借口来逃避，而是尽到一名教师应尽的责任。

南京师范大学泰州学院本科毕业论文

18

谢辞

时间如白驹过隙，还未来得及细细品味，大学生活就已经尾声，四年的努力与付出，

伴随这篇论文的结束，将要画上一个句号。

这边论文的设计在刘陆军老师认真负责的指导下成功完成，无论是选择课题，还是

搜寻资料进行写作，都离不开刘陆军老师的热情帮助。他不仅给予了我专业知识上的帮

助，还在精神层面上为我加油打气，他的关怀对于容易焦躁的我来说就像是一副定心丸，

使得我能够顺利的完成自己的工作。

快要毕业了，衷心的对我的各科老师表示感谢，他们教学认真负责，使得我能够更

好的掌握和运用专业知识，并在论文中体现出来，顺利完成毕业论文。

在论文写作过程中，我参考了很多优秀的学者的作品，在这里对每一位认真负责的

学者们表示感谢，你们是为人类进步铺路的工人。

对于我的各位同学和舍友，在我毕业论文设计的这段时间里，你们给了我很多帮助，

尤其是徐师鉴同学，帮助我一起校对了论文的格式，对于你们大家的鼎力相助，我表示

深深的感谢！

南京师范大学泰州学院本科毕业论文

19

参考文献

[1](瑞典)L.戈丁(g)著,胡作玄译.数学概观[M].科学出版社,2001：

121-147

[2]赵慈庚,朱鼎勋主编,大学数学自学指南[M].中国青年出版社,1984:74-91

[3]王高雄等编,常微分方程[M].高等教育出版社,1978:39-53

[4]李瑞遐,何志庆编著.微分方程数值方法[M].华东理工大学出版社,2005：41-58

[5]余德浩,汤华中编著.微分方程数值解法[M].科学出版社,2003：14-21

[6]胡燧林.一阶方程初值问题解的存在与唯一性定理的几点注记[J].韶关学院学报.

1988(02)：38-47

[7]弭鲁芳,纪在秀.论一类常微分方程解的最大存在区间[J].聊城大学学报(自然科

学版).2006(04):24-26+28

[8]赵慧娟,陈伟丽,赵晨霞,袁书娟.关于常微分方程初值问题数值解法的分析[J].中

国科教创新导刊.2011(08):83

[9]李孝诚,刘兆丽.常微分方程解题模式的构建[J].高等数学研究.2009(04)：96-99

[10]韦程东,高扬,陈志强.在常微分方程教学中融入数学建模思想的探索与实践[J].

数学的实践与认识2008(20)：228-233

[11]舒小保.一类二阶常微分方程边值问题的无穷多个解[J].系统科学与数学.

2008(01)：91-98

[12]《邯郸大学学报》（自然科学版）2000年第1期

[13]张鹏。二阶常系数线性微分方程讨论[J].工科教学，1994.10(4):30-32第五期

[14]同济大学数学教研室.高等数学[M].北京：高等教育出版社，1996,(4):387-394