二阶微分方程解法.pdf

《二阶微分方程解法.pdf》由会员分享，可在线阅读，更多相关《二阶微分方程解法.pdf（20页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、毕毕 业业 论论 文（设文（设 计）计）题题目：目：二阶常微分方程的解法二阶常微分方程的解法院（系、部）：院（系、部）：数学科学与应用学院数学科学与应用学院专专业：业：数学与应用数学数学与应用数学摘要摘要：本文主要是介绍了二阶常微分方程众多解法中的三种，分别为特征方程法，拉普拉斯变换法和常数变易法，研究并讨论了二阶常微分方程在特征方程法中特征方程根为实根，复根和重根的情形。我们选用了弹簧振子系统的振子运动，用这三种不同的方法来解决该问题。关键词关键词：二阶常微分方程；特征根法；常数变易法；拉普拉斯变换Abstract:Abstract:The main purpose of this pape

2、r is the second-order ordinary manydifferential equation solution of three,respectively as the characteristic equationmethod,Laplace transform method and variation of constants method,study anddiscuss the second-order often differential equation in the characteristic equationof the roots of the char

3、acteristic equation for real roots,complex roots and rootweight.We choose the spring oscillator the oscillator motion,these threedifferent methods to solve the problem.Keywords:Keywords:second order ordinary differential equation;Characteristic analysis;constant variation method;Laplasse transform1目

4、录1 绪论.31.1 二阶常微分方程的起源和发展史.31.2 二阶常微分方程的介绍.31.3 研究二阶常微分方程的目的与意义.42 二阶常系数常微分方程的几种解法.52.1 特征方程法.5 2.1.1 特征根是两个实根的情形.5 2.1.2 特征根有重根的情形.62.2 常数变易法.72.3 拉普拉斯变换法.93 二阶常微分方程解法的应用（分析例题）.113.1 特征方程法.113.2 常数变易法.133.3 拉普拉斯变换法.144 结论和启示.16谢 辞.18参考文献.1921 绪论1.1 二阶常微分方程的起源和发展史既然说到了微分方程，就不能不提到海王星的故事，它的发现是人类智慧的硕果，微

5、分方程在其中扮演了重要的角色，并且在其中也包含数学演绎法的作用。在发现了天王星之后，进行天文观测的人们发现它所处的位置总是和万有引力计算出的位置有些许不同，于是有人质疑万有引力定律的正确性。但也有一部分人认为，这也可能是天王星在受到一颗尚未发现的行星的吸引力才会造成的改变。不少人坚信这种假设是正确的，但却很少有人找到了正确的方法并付诸实践。而英国的一个学生亚当斯显然不是其中之一，他勇敢接受了这项任务，运用手头仅有的资料建立起微分方程，成功求出了海王星的位置与下次出现的时间。1843 年 10 月 21 日，满怀信心的亚当斯把结果寄给了天文台长艾利，但换来的却是质疑，艾利并不相信籍籍无名的他。然

6、而在两年后，另一名青年勒威耶也计算出了同样的数据，并把计算的结果给予了位于柏林天文台的助理员卡勒，在那个值得铭记的夜晚，卡勒在计算出的位置上发现了第七颗行星海王星。从二十世纪三十年代以来，常微分方程的研究像是走上了快车道，迅速发展并建立起了多个分支。19271945 年期间定性理论的主要研究是与无线电技术紧密联系在一起的。在第二次世界大战期间由于对通讯等方面的需求越来越高，极大地激发了对无线电技术的研究进展，尤其是对非线性振动理论的研究取得了迅速的发展。在四十年代之后各国大部分数学家们主要在研究对抽象动力系统的拓扑特征,例如闭轨的存在性、结构的稳定性等,对于二维系统来说，我们可以通过一些方法证

7、明他的结构稳定性；而对于一般的系统来说这个问题依旧困扰着我们。在动力系统的研究方面,目前采用的办法是从典范方程组到阻碍集有详尽的理论指导,成功解决了一系列困扰人类多年的问题,其中最为突出的是 C封闭引理的证明,以及对结构稳定性的充要条件等方面都作出了杰出贡献。在当今社会，由于信息技术的飞速发展，大量的领域需要用到常微分方程组进行描述。前赴后继的杰出的学者们，为了各种稳定性及专业问题，终其一生都在研究常微分方程，也取得了不朽的成就，但依然有很大的疑问等着我们去解开。1.2 二阶常微分方程的介绍二阶微分方程在时间上大致与微积分同时产生。对于初学者来说，y f(x)这样的问题就是最简单的微分方程了。

8、二阶常系数线性微分方程是形如y py qy f(x)的微分方程。与其对应的二阶常系数齐次线性微分方程为y py qy 0，其中p,q是实常数。若函数y1和y2之比为常数，称y1和y2是线性相关的；若函数y1和y2之比不为常数，称y1和y2是线性无关的。31.3 研究二阶常微分方程的目的和意义就如上文所说，研究二阶常微分方程已经取得了不少成就，尤其是在对方程的求解方面。与此同时应用常微分方程的理论也经过大家的不懈努力，取得了丰硕的成果。但是对于进一步的发展所需，还是小巫见大巫，所以我们要更加努力的钻研，努力地完善这门学科的理论系统。在数学的发展历史中，数学分析占有非常重要的地位，我们大学学习的课

9、程也都是以数学分析作为基础，而微分方程正是数学分析的关键所在。同时它也发展出了数学分析中大部分思想以及理论。众所周知的，常微分方程自始至终都是人类用来探索自然变换，研究自身社会结构，工程问题以及大自然的生态结构的便利的道具。常微分方程由于与现实生活息息相关，所以对其的研究一直没有停止过，而且表现出欣欣向荣的活力。并且在多个学术领域中，常微分方程都占着决定性的作用，可以说常微分方程带领着人类的进步。而二阶常微分方程同样在常微分方程的整套理论中有着弥足轻重的地位，在各个研究领域中都有十分广泛的应用。42 二阶常系数常微分方程的几种解法2.1 特征方程法特征方程法中的特征方程，是为了对相对应的数学对

10、象进行深入的研究而人为引入的一些等式，当研究的对象改变时，它也会改变，这其中包括数列特征方程，微分方程特征方程和积分方程特征方程等等。d2xdx求微分方程2 p qx 0的通解。dtdt解：特征方程2 p q 0的根1，2,(1)若1，2是两个不相等的实根,那么上面这个微分方程就拥有两个实值解e1t,e2t,于是我们可以求得其通解为x c1e1tc2e2t（c1,c2为常数）.(2)若1，2相等,则该方程有二重根,因此方程的通解具有形状如x c1e1tc2e1t（c1,c2为常数）.(3)若1，2为共轭复根z a bi的情况,则该方程的通解具有形状x eat(c1sinbt c2cosbt)（

11、c1,c2为常数）.在数学中，许多公式与定理都需要进行证明,下面本文给出前两个解答的理论依据及证明过程。2.1.1 特征根是两个实根的情形设1，2为该特征方程的两个不等实根,那么我们可以得出与之对应的方程的两个解为e1t,e2t,我们确定这两个解在a t b上线性无关,所以它们能够组成该方程的基本解组。事实上,这时w(t)e1te2t1e1t2e2t e(12)t1112,而上面最后一个行列式是著名的范德蒙德（Vandermonde）行列式,该行列式与（2-1）相等。由于在之前假设了21,所以此行列式不等于零,从而w(t)0,于是得5到e1t,e2t线性无关,这就是要证明的结论。而该方程的通解

12、即可表示为x c1e1tc2e2t（其中c1,c2为任意数）.2.1.2 特征根为复根时讨论这个特征方程含有复根的情况,首先这个方程的系数为实数,同时也是常数，所以复根成对共轭出现。设1i是特征根之一,则2-i就是第二个特征根,根据这两个特征根我们就可以得到原方程的两个不同的实值解:e(i)t ei(cost isint),e(i)t eai(cost isint).根据定理,我们得到的复值解的实部和虚部与方程的解相同。那么由于这一对共轭d2xdx复根i对应于特征方程,可以求得方程2 p qx 0的两个实值解dtdteaicost,eaisint.2.1.3 特征根有重根的情形设特征方程有k重

13、根1则易得F(1)F(1)F(k1)(1)0,F(k)(1)0,先设1 0,表示特征方程有因子k,于是an an1 ank1 0,也就是该特征方程的形状为na1n1ankk 0,dnxdn1x而与其对应的方程Lxna1n1an1anx 0则变为dtdtdnydn1ydkya1n1ankk 0.dxndxdx易得它有k个解,t,t2,tk1,而它们是线性无关的。于是我们可以得出,特征方程中6的k重零根对应方程的k个彼此线性无关的解,t,t2,tk1。当这个k重根 0,我们作变量变换x ye1t,可注意到xmye1t可得mmm12m2m e1tymm1ym11y1y，2!L ye1tdny1tdn

14、1y1t,bb y e L y e1n1n1dtndt于是可将对应方程化为dnydn1ydn1yL1ynb1n1b1n1bny 0,dtdtdt而其中b1,b2,b3,bn仍为常数,其相应的特征方程为Gnb1n1bn1bn 0,直接进行计算易得F1e1t L e1t L1ete1t Ge1t,因此 F1 G,从而Fj(1)G,j 1,2,k.通过这样转化,问题就化为前面所讨论过的情况了。2.2 常数变易法接下来我们要说的是求解微分方程的一种极其重要的方法,常数变易法。我们在求解一阶线性微分方程时常常采用这种方法。这个办法的大体思想，是通过将常数N代入UX中来得到该方程的通解。看似简单，但拉格朗

15、日为之研究奋斗了11年，而我们所如今所采用的仅仅是他所得出的结论，证明过程省略了。它在非齐次线性微分方程和与其对应的齐次线性微分方程之间起着重要的链接作用。我们这里讨论的对一般二阶常微分方程的求解,要先得出该方程的一个特解,再用我们上面讲解的特征方程法求通解。d2xdx例题 求常微分方程2 p qx ft的通解。dtdt7d2xdx解方程2 p qx ft与其对应的齐次方程为dtdtd2xdx p qx 0,2dtdt其特征方程为2 p q 0.d2xdx由于2 p qx ft的通解根据理论可以得到就等于这个方程所对应的齐次dtdt线性微分方程的通解与求出的它本身的一个特解之和，而我们已经讨论

16、了二阶常系数齐次线性常微分方程的通解,所以只要再求出一个该方程本身的特解。2d xdx若为该方程的实根,则x e为方程2 p qx 0的解。根据常数变易法dtdttd2xdxt p qx f t的 其 中 一 个 解 为,代 入 原 方 程 并 进 行 化 简 得x c t e2dtdtct2 pct etft,这是即是关于ct的一阶线性微分方程,而其一个特解为c(t)e2pteptftdtdt,从而得到上面方程的一个特解为x ete2pteptftdtdt.讨论为该方程的复根的情况,可以设 a bia,bR,b 0，则x eatsinbt就是d2xdx原方程2 p qx fx的解，接下来用常

17、数变易法设其中的一个特解为dtdtd2xdxx cte sinbt，通过运用情形 1 解得方程2 p qx ft的一个特解为dtdtatp2atp2aef t esinbtdtatx e sinbtdt.2sin bt*x由于是特解,所以积分常量可以都取零。82.3 拉普拉斯变换法接下来要介绍的拉普拉斯变换法是一种积分变换法，又名为拉氏转换法。它是一个线性变换法，通过因数为实数tt 0，的函数中的实数转化为复数i。在某些情形下一个实变量函数在R中进行运算比较困难，若能将实变量函数进行拉普拉斯转换，并在复数域中进行各种运算，再使用拉普拉斯反变换来求得该方程或函数在实数域中的相应结果，是一种简便的

18、计算方法。拉普拉斯转换的运算步骤对于我们要求的解线性微分方程极其有效，它可以把微分方程转化为较易求解的代数方程来进行处理，使计算大幅度简化。在经典控制理论中，对于控制系统的分析但和综合等，都是立足在拉普拉斯转换的基础上的。我们引入拉普拉斯变换有一个主要优点，就是可以采用传递函数来代替常系数微分方程用以描述系统的特性。当我们要用图解来测定某些控制系统的整体特性并且分析该控制系统的运动过程等情况时，可以用上这种方法，一般来说比较简单。steftdt.由积分Fs0所定义的确定在复平面Re 上的复变数s的函数Fs,称其为函数ft的拉普拉斯转换,我们称ft为原函数,而将Fs称为像函数。拉普拉斯变换法主要

19、是通过借助拉普拉斯变换将常微分方程转化，使它成为复平面s的代数方程。并通过一些简单的代数运算,再利用拉普拉斯变换表,即可对微分方程进行求解。方法十分简单方便,在现实生活中也经常被采用。但是这个方法也存在着些许不足，主要是只适用于右端函数为原函数的情况。d2xdx求解方程22 x et,x1 x1 0.dtdt解首先使 t 1，将问题转化为d2xdx2 x et1,x0 x0 0,2dtdt再对新方程两边进行拉普拉斯变换,得到s2Xs 2sXs Xs11,s 1 e因此Xs11,s 13e查阅拉普拉斯变换表可得1x2e1,2从而得到xt1t 12et,29这即为所要求的解。当然,求解二阶或更高阶

20、的常微分方程我们能用的方法还有不少,这里由于篇幅我不能将其一一列出。但我们利用前面的一些结论就足以解决下面的几个物理问题了。103 二阶常微分方程解法的应用上面介绍了解二阶常微分方程常用的三种方法，为了更好的体现这三种方法的区别，我选取了一个比较简单的动力学方程来进行求解。通常，在对物理问题求解时，我们有三个步骤：第一步是对该问题进行彻底分析从而能做到对方程的初步建立并且对定解条件进行明确；第二步是对其解的性质进行探究或者求出方程来满足初始条件的特解；最后一步是定性分析对解，对原来的问题反过来进行解释，其中最为关键的步骤就是要将方程列出，而能列出方程的方法主要有两种，分别为微元分析法和瞬时变换

21、法。但是在研究阻尼运动的过程中，求解运动方程一直是令人头疼的问题。接下来我们就分别使用上述三个方法来求解这个动力学方程。3.1 特征方程法在我们要研究的弹簧振子系统中，测定物体的阻尼系数10.0s1，物体的质量为m 1.0kg,该弹簧所具备的劲度系数为k 75N m1,根据上述条件，假设质点由静止状态逐步开始运动，对弹簧振子的位移方程进行求解。解：根据牛顿的第二运动定律我们可以得出 kx cv ma1或d2xdxm2 c kx 0,2dtdt由于相对来说振动系统这是在这之前给定的，其中包含的常量为m,k,c，如果能够2确定k/m 0,c/m 2,那么上面的方程式就可以转变为：dd2 20 0,

22、32dtdt那么根据所得到的数据代入公式（3）我们就可以得到d2xdx22075x 0.4dtdt通过对以上公式的仔细观察和深入研究则可以得到对该方程进行求解能够使用特征值法，得出其特征方程可以表述为：2 2075 0,并且在该特征方程当中包含有不同的两个根1 15,2 5,这样相对应的公式4的两个根分别为1 e5t,2 e15t511那么根据公式5进行计算可以即可得到固有角频率数值为0k m 5 2，此时2阻尼系数值为10,也就是说20,于是方程5的解可以表述为tA,B数值）Ae5t Be15（初始条件可得.6通过公式6，我们可以看出振子所保持的属于一个非振动状态，在这样的背景下，我们要求的

23、质点也只是在原先的不平衡位置逐渐恢复到平衡状态当中，该质点并不具有周期振动的特征。由于我们的关注点是在0情况下面,质点会呈现出逐渐衰减的振动。然而由于会受到阻尼作用的影响，他不能保持自由振动系统的长久运作，振动会逐渐衰减直至振动停止，如果我们要保持振动持续不停的状态的话，该质子就需要从外界获得运动必要的能量，在学术界将这种由于受到外部的持续作用而产生振动的情况称为强迫振动。我们再举一个例子：假如在以上的振动系统当中振子由于受到某个外力F 100cos30tN的作用，在公式之中FA100表示的是该驱动力的幅度值，30表示的是该驱动力所具有的圆频率，f即是该驱动力所保持的频率。解：由于在质点振动系

24、统当中会受到驱动力的作用，那么我们就可以得到关于系统振动的方程：d2xdxm2ckx F,7dtdt上述方程还可以表示为d2xdx20 x H cos30t.82 2dtdtFA表示的是在单位质量上该振子所受到的外力大小。而7和8这两个式子都m是强迫运动方程。我们观察该式子的本质发现，他们属于二阶非齐次常系数微分方程，而根据定理我们对这个方程求解后所能得到的一般解同时也是该方程所能得到的一个特解和与该非齐次方程相对应的齐次方程的一般解之和。我们在前面求出了其相对应的H 自由振动方程的一般解，那么现在要做的就是找到8方程的一个特解，将其代入方程8就能得到方程9：d2xdx22075x 100co

25、s30t,9dtdt在这里可以假设方程9有着x1 Asin30t Bcos30t这样的特解，将这个特别对方12程9替代进去并且将其简化最终得到33A24Bsin30t 24A33Bcos30t 4cos30t,通过比较同类项系数我们可以得到A 出x13244，B ,这样我们就可以进一步得5555553244sin30t cos30t,而根据前面所得到的结果，可以将原方程的通解表示出555555来，为：3244x(t)Ae5t Be15tsin30t cos30t.555555综上所述，题目中给我们的条件决定A,B的数值，之前的两项为该方程的瞬态解，瞬态项对这整个系统所进行的自由衰减振动能够进行

26、有效的描述，然而只能在运动的开始阶段起作用，在经过长时间的运动之后，它起到的影响会随时间消逝并且在运动最后完全消失。而之后两项所代表的稳态解，描述的是强迫运动的状态，由于幅值条件的固定，所以我们称这样的状态为稳定状态。根据上面的公式，外力作用在质点上的时候，整个系统的振动状态十分复杂，这时候的振动既包括了瞬态振动，也包括了稳定振动，而这样的振动状态对于在强迫振动之中逐步建立起稳态振动的过程进行有效的描述。在长时间的振动之后，瞬态振动终将消逝，这整个系统就会保持稳态振动的状态。3.2 常数变易法从上面的分析之中我们可以了解到x e5t即是属于特征方程2 2075 0的实根，于是我们就可以得到x

27、e5t为属于上面方程9当中的一个根，之后通过运用常数变易法设置x cte5t，在这一过程中我们也可以得到方程其中一个解为x，将数值代入到方程9当中并且对其进行简化可以得到ct10ct e5t100cos30t.属于ct的一阶线性微分方程,并且得到在该方程当中一个特解为84cte5tsin30t e5tcos30t c1,33从而得出方程9的一个特解为（取c1 c2 0）4 8xt e5te5tsin30t e5tcos30tdt1c23 3133244sin30t cos30t,555555于是可得方程9的通解为xt Ae5t Be15t3244sin30t cos30t.555555由上面的

28、结论可知d2xdxm2 ckx F.10dtdt将题目数据代入公式中可以得到d2xdx220400 x cos2t.11dtdt通过观察我们可以发现，在进行该类求解的过程中要使用常数变易法，首先就是要求 出 公 式11，而 在 之 前 的 探 索 当 中 已 经 可 以 得 到 公 式11的 特 征 方 程 为220400 0.于是可以进一步的假设该特征方程的根为 103i，那么xt e10tsin10 3t即为公式11的一个解。运用常数变易法可设为xt cte10tsin10 3t.*x这里套用情形 1 的解法，将(t)代入12并进行化简得xt1099sin2tcos2t.396043960

29、4由于x是特解,所以积分常量可以为零。3.3 拉普拉斯变换法依然使用上面的例子，通过牛顿第二运动定律可以得到d2xdxm2ckx F,dtdt将公式代入数据之后我们可以得到d2xdx220400 x cos2t,12dtdt根据条件，质点由静止状态开始运动，我们可以得到以下的公式14xt0 0,dx 0,dts,s2 4对方程(12)进行拉普拉斯变换得到s2Xs 20sXs 400Xs即X(s)s1,22s 4 s 20s 400将上式右端分解为部分分式Xs10299s2239604 s 439604 s 4101 310 399s 1022.22118812s 10 10 339604s 1

30、0 10 3通过拉普拉斯变换表可得xt1099sin2tcos2t3960439604101 310t9910tesin10 3t ecos10 3t.11881239604154 结论和启示在本篇文章中，我们首先介绍了海王星的发现过程，并介绍了研究二阶常微分方程的背景及研究现状，之后重点总结了解二阶常微分方程最具有代表性的三种不同的解法，特征值法，常数变易法和拉普拉斯转换法，并且进行了一定程度的探究。在特征方程法中，我们通过研究特征方程的特征根的情况，来判断原方程的根的情况。若1，2相等,即该方程有二重根，这是第一种解的情况。而1，2为共轭复根z a bi时，为第二种情况，并给出了这两种情况

31、的证明过程。常数变易法则是使用了常数N代入Ux的方法，并且提出了定理，原方程的通解就等于这个方程相对应的齐次线性微分方程的通解，同时也运用了特征方程法的思想。最后我们介绍了比较简便的拉普拉斯变换法，它是一个线性变换法，通过因数为实数tt 0，的函数中的实数转化为复数i来进行运算。对某些情况下的实变量函数进行拉普拉斯转换，使得我们能在复数域中对其进行各种运算，最后再使用拉普拉斯反变换，得出结果，这是一种简便的计算方法。该运算步骤对我们这篇文章要求解的线性微分方程是相对来说最好的办法，它可以把微分方程转化为较易求解的代数方程来进行处理，使我们不用再进行繁琐的计算。之后我们用这三种方法求解了弹簧振子

32、的运动方程，以此来比较三种方法的相同点与不同点。并且通过这个方程的解答，体现出了二阶常微分方程与我们的生活息息相关的特性。由于篇幅及本人专业知识的限制，无法进行更深层次的探究，比如求解二阶常微分方程还有一种幂级数解法，但是它的解题过程非常繁琐，对计算的要求很高，计算量比较大，还要考虑该函数是否解析以及幂级数在某个区间是否收敛等问题，所以在此不做讨论。二阶常微分方程研究的道路远不止此，在这个时代，二阶常微分方程被广泛应用于网络，军事，医学等各种高科技领域，对二阶常微分方程的研究，是对目前社会发展的推进。类似于本篇中写出的弹簧振子方程，研究二阶常微分方程不仅能够推进自身学科的进步，还能提升别的学科

33、的发展速度，它的各个分支学科也与人类生活息息相关。总而言之，在高速发展的现代，常微分方程的研究工作虽然仍有很长的路要走，但由于其重要性不言而喻，所以我们要全身心的投入到研究中去。对于二阶常微分方程，目前并没有研究出通用的求解办法，只能对少数特殊情况的方程进行求解。这无论是对现实生活还是研究工作都是非常不利的，所以对常微分方程的求解问题，任重而道远。我们唯有脚踏实地，艰苦奋斗，为了人类更好的明天而努力。而通过本次论文的写作，我也明白了有志者事竟成这个道理，无论是开头讲的海王星的发现，还是在写作过程中读到的各个伟人的事迹，都告诉我只要努力了，终会得到回报，或许得到的回报有大有小，甚至为目标努力奋斗

34、了许久却仍未得到回报，这时更不能放弃，继续做下去还有希望，而如果放弃了得到的就一定是失败。然而，不相信自己的人，没有努力的价值，所以，相信自己，勇敢拼搏，最终定能获得成功。在对这篇作文的写作过程中，对物理学方程的研究遇到了一些麻烦，因为我的物理16学知识比较薄弱。但是通过同学、老师的帮助，再在网络及学校图书馆翻阅了各种文献之后，我还是成功的做成了。这让我想到了，在今后的教育工作中，我应该保持这种勇于探索的精神，不能以非专业，没学过为借口来逃避，而是尽到一名教师应尽的责任。17谢 辞时间如白驹过隙，还未来得及细细品味，大学生活就已经尾声，四年的努力与付出，伴随这篇论文的结束，将要画上一个句号。这

35、边论文的设计在刘陆军老师认真负责的指导下成功完成，无论是选择课题，还是搜寻资料进行写作，都离不开刘陆军老师的热情帮助。他不仅给予了我专业知识上的帮助，还在精神层面上为我加油打气，他的关怀对于容易焦躁的我来说就像是一副定心丸，使得我能够顺利的完成自己的工作。快要毕业了，衷心的对我的各科老师表示感谢，他们教学认真负责，使得我能够更好的掌握和运用专业知识，并在论文中体现出来，顺利完成毕业论文。在论文写作过程中，我参考了很多优秀的学者的作品，在这里对每一位认真负责的学者们表示感谢，你们是为人类进步铺路的工人。对于我的各位同学和舍友，在我毕业论文设计的这段时间里，你们给了我很多帮助，尤其是徐师鉴同学，帮

36、助我一起校对了论文的格式，对于你们大家的鼎力相助，我表示深深的感谢！18参考文献1(瑞典)L.戈丁(Lars.Garding)著,胡作玄译.数学概观M.科学出版社,2001：121-1472 赵慈庚,朱鼎勋主编,大学数学自学指南M.中国青年出版社,1984:74-913 王高雄等编,常微分方程M.高等教育出版社,1978:39-534 李瑞遐,何志庆编著.微分方程数值方法M.华东理工大学出版社,2005：41-585 余德浩,汤华中编著.微分方程数值解法M.科学出版社,2003：14-216 胡燧林.一阶方程初值问题解的存在与唯一性定理的几点注记J.韶关学院学报.1988(02)：38-477

37、 弭鲁芳,纪在秀.论一类常微分方程解的最大存在区间J.聊城大学学报(自然科学版).2006(04):24-26+288 赵慧娟,陈伟丽,赵晨霞,袁书娟.关于常微分方程初值问题数值解法的分析J.中国科教创新导刊.2011(08):839 李孝诚,刘兆丽.常微分方程解题模式的构建J.高等数学研究.2009(04)：96-9910韦程东,高扬,陈志强.在常微分方程教学中融入数学建模思想的探索与实践J.数学的实践与认识 2008(20)：228-23311舒小保.一类二阶常微分方程边值问题的无穷多个解J.系统科学与数学.2008(01)：91-9812邯郸大学学报（自然科学版）2000 年第 1 期13张鹏。二阶常系数线性微分方程讨论J.工科教学，1994.10(4):30-32 第五期14同济大学数学教研室.高等数学M.北京：高等教育出版社，1996,(4):387-39419