线性常微分方程的若干初等解法探讨 优秀毕业论文.doc

线性常微分方程的若干初等解法探讨作 者：xxx指导教师：葛玉丽 摘要：介绍求解常微分方程的几种初等解法，如常数变易法，积分因子法，拉普拉斯变换法等，在学习过程中，通过对不同类型的方程解法，揭示了常微分方程的求解规律，从而找到最优解法.关键词：常数变易法；积分因子；特征根法；拉普拉斯变换 0引言常微分方程是数学分析或基础数学的一个组成部分，在整个数学大厦中占据着重要位置，学好常微分方程基本理论与方法对进一步学习研究数学理论和实际应用均非常重要，对于常微分方程的初等解法，既是常微分方程理论中有自身特色的部分，也与实际问题密切相关；恰当对初等解法进行归类，能正确而又敏捷地判断一个给定的方程属于何种类型，从而能按照所介绍的方法进行分解.1 一阶常微分方程的求解方法1.1方程能解出1.1.1变量分离方程形如的方程称为变量分离方程.分别是的连续函数.例1 .解 将变量分离 得 ;两边积分得: ;因而通解为： （为任意常数）.这是一种相当简洁的解法，是最基本的解法，对于比较复杂的方程，需经过一系列变换，最后利用变量分离求解.1.1.2 常数变易法对于一阶线性齐次方程 它的通解为 从此出发，将通解中的任意常数换成待定函数，假设 （1）为一阶线性非齐次方程 （2）的解，为了确定，将（1）代入的左边，得到 从而得到 ， 即 积分后得到 ，其中为任意常数把代入（1）中，得到方程（2）的通解为 例2 解方程：.解 方程变形为 令，则 ;代入变形方程为： ;利用常数变易法，其中;则它的通解为 ;代回原来的变量，得到 ;即原方程的通解为 ;此外，方程还有解 .常数变易法实际上也是一种变量替换法，虽然用其来解一阶非齐次线性微分方程时和变量代换法并无原则区别，但将它推广到解高阶线性微分方程和线性微分方程组时就显出了它的优越性，变易常数思想是解微分方程的重要数学思想，对非线性方程（如贝努利方程，黎卡提方程）也可使用常数变易法求解，并且常数变易法在数学分析中有很多应用，比如求解中值问题及存在性问题，祥见文献1.1.3 积分因子法把一阶线性微分方程 (1）改写为如下的对称形式： （2），一般而言，（2）不是恰当方程，但以因子 乘（2）两侧，得到方程：，即 它是恰当方程，由此可直接积分，得到 这样就求出了方程的通解 （3） 为任意常数，其中为积分因子，一般情况下，积分因子是很难寻求的，只有在很特殊的情况下才很容易求得.例3 求解.解 因为;则方程不是全微分方程，若把原方程改写为可以看出积分因子，因为上式两端同乘以，有;即 从而得到方程的通积分 ，或.此解法，目的明确，方法自然，学生很容易接受，逐步改变了一上来就直接用任意常数变易法求解一阶线性微分方程的方法，取而代之是按上述方法一步步求解，这一过程使我们顺利掌握了一阶线性微分方程的通解，同时更容易理解任意常数变易法，这样从不同角度，用不同方法解决了同一问题，更能深刻的体会到任意常数变易法的巧妙之处.1.2 方程不能解出这时把看作是的函数，再看是否能解出，成为方程可用以上方法求解；但对于不能显性表示为或或的方程，可分为两类：1.2.1 方程能就（或）解出（或）这时令（或）把问题转化为求解关于与（或）之间的一阶方程 （或），再利用以上方法，求得通解为 （或）则它与（或）一起构成原方程的通解的参数形式.例4 研究克莱洛（claivaut）方程 （1）.解 令代入原方程 假定两次可微且;两端对求导，得 取 则;代入（1）得到通解 取，则即（2）;由于，则（2）中第一式存在隐函数，代入第二式就得到一个解，则这个解也可以由联立方程 来表达.故克莱洛方程除了通解之外，还有一个由所决定的解.例5 求解.解 令，代入原方程 ;两边同时对x求导，则,则 ,则当时，；当时，则，为任意常数,则得到方程参数形式的通解， ；且当时，也是方程的解.总结：由于此方程的形式与前面所分析的类型不一致，可以先观察所给的方程的形式，利用变量代换的思想，经过一系列变换，化为我们最熟悉的形式.1.2.2方程不能就，或解出对于形如或的方程，引入参数，将方程表示为参数形式，再注意到关系式，就将问题转化为求解关于（或）与的一阶方程，且其导数（或）已表示为的已知函数，最后的工作就是求积分的问题.例6 求解.解 令，则原方程可化为：，则 ，；由于 ，则 ，两边同时积分，则；则原方程的通解为，.例7 .解 令，代入原方程为；则；由 ，则，；即 ,两边同时积分：；则原方程的通解为 ，.以上总结了一阶常微分方程的几种解法，熟悉各种类型方程的解法，正确而又敏捷地判断一个给定的方程属于何种类型，从而按照所介绍的方法进行求解，这是最基本的要求.但是我们所遇到的方程未必都恰好是所介绍过的方程类型，因此要注意学习解题的技巧，善于根据方程的特点，引进恰当的变换，将方程化为能求解的新类型，从而求解；一阶微分方程的求解有众多方法，技巧性很强，想进一步详细了解可参考常微分方程手册.2高阶常微分方程的求解方法高阶常系数线性微分方程的一般形式是 （1）其中为常数，为连续函数；依据常系数线性微分方程的通解结构理论，知方程（1）的通解可表示成该方程的一个特解与其对应的齐次方程的通解之和. 方程（1）对应的齐次方程，由于它具有线性结构，一般采用euler待定指数函数法可以得到通解，因而非齐次方程（1）的通解的计算只需寻到它的一个特解即可；有关特解的计算方法较多，如常数变易法，待定系数法，积分法等，因此接下来介绍线性微分方程的求解方法的几种归类.2.1常数变易法例8已知齐次线性微分方程的基本解组，求下列方程对应的非齐次线性微分方程的通解：.解 应用常数变易法，令,将它代入方程，则可得: 解得: ;由此 ；则原方程的通解为 .总结：利用一阶常微分方程的常数变易法的思想，推广到高阶常微分方程，关键是找出决定的方程组，从而求出高阶方程的通解. 由此可知，常数变易法一般用于给定非齐次线性微分方程特解的方程，这种方法简洁明了，但是比较局限，是最基本的解法.2.2特征根法主要是利用把微分方程的求解问题化为代数方程的求根问题的思想.我们知道简单的一阶方程，其中为常数，它有特解，由于与都是常系数线性齐次方程，因而猜想方程也有形如的解，其中是待定常数，为了确定出使为的解的，先将它代入方程中，实际上有，其中称为特征多项式.则为方程的解的充要条件是，即应是方程的根.下面分两种情况讨论：特征根互异：首先，假设有个互异的实根，这时，依上述讨论，方程有个特解，则函数为方程的通解，其中为任意常数.例9求方程的通解. 解 特征方程为，故特征根为，因而基本解组为,故所求通解为，其中为任意常数.特征根有重根：设是重特征根（），由上述讨论知，是的一个解，但这时由于互异的特征根的个数小于，故相应地线性无关的解的个数也小于，要得到通解，这些特解是不够的，对应于，除解外还应补上哪些解呢？ 先来研究二阶常系数方程并设，特征方程为，特征根为，即；易见，为二重特征根，因而，首先有特解；现在求已知方程的和线性无关的另一个特解，由知，取,则另一特解可取为，即当是二重特征根时，二阶方程除了有解之外，还有与它线性无关的另一个特解.根据以上讨论，对于一般的情形，我们有如下的定理：如果方程有两两互异的特征根，它们的重数分别为，且，则与它们对应的方程的特解是.例10 求方程的通解. 解 特征方程是故特征根是 ,则它们对应的解为： ，故所求通解为： ，其中为任意常数.总结：欧拉待定指数函数法，即特征根法，在高阶常微分方程中占据了十分重要的位置，要熟练掌握不同类型的解法，从而对于给定的方程能游刃有余.2.3 阶常系数线性非齐次方程解法对于形如的解法，它的通解等于其对应的齐次方程的通解与它本身的一个特解之和.2.3.1 比较系数法（待定系数法）下面分两种类型讨论:设，其中及为实常数.当不是特征根时，有形如的特解，其中当是（）重特征根时，有形如的特解，其中，对于中的的系数，则可以由待定系数法求得.例11求方程的通解解 先求对应齐次方程的通解，其特征方程是;故特征根为 从而，对应齐次线性方程通解为;由于不是特征根，因而已知方程有形如的特解.为确定将它代入原方程中，由于,故 .比较上式等号两端的同次幂系数，可得 ,故已知方程特解为，则原方程的通解为 .例12 求方程.解 由于则故齐次方程通解为： ，由于为二重特征根,故有 ,故 ,则原方程的通解为 .设，其中为常数，而是带实系数的多项式，其中一个的次数为，一个的次数不超过，则有形如的特解.其中为特征方程的根的重数，而均为特定的带实系数的次数不高于的的多项式根据欧拉公式，有则再利用迭加原理,于是有两种形式：（1） 如果不是特征根，则特解具有形式其中是系数待定的次多项式.（2）如果是重特征根，则特解应具有形状.例13 求解方程.解 先求对应的齐次方程，我们有,故特征根为；由于迭加原理，则原方程可化为(1)对于，由于是特征根，故方程具有形如的特解，现将上式代入，则;则的通解为.(2)对于，由于不是特征根，故方程具有形如的特解.现将上式代入，则,则的通解为.故原方程的通解为.总结：比较系数法用于方程右端是某些基本函数的情况，常见的有：多项式，指数函数，正弦（或余弦）函数以及它们的某种乘积组合，然后根据的前面所归纳的类型，从而求出方程的特解，进而求出通解.2.3.2 拉普拉斯变换 它无需求出已知方程的通解，而是直接求出它的特解来，从而在运算上得到很大简化，这一方法的基本思想是：先通过拉普拉斯变换将已知方程化为代数方程，求出代数方程的解，再通过逆拉普拉斯变换，便可得到所求初值问题的解.由积分所定义的确定于复平面上的复变数的函数称为的拉普拉斯变换，其中与有定义，且满足不等式，这里m,为某两个正常数，这时为原函数，而称为像函数.例14 求函数的拉普拉斯变换.解 .例15 解方程.解 由于，从而则 ，故 ，由于 ，故所求初值解为 .当然，方法本身也有一定的局限性，它要求所考察的微分方程的右端函数必须是原函数，否则方法就不适用了，关于拉普拉斯变换的一般概念及基本性质，请参阅有关书籍.2.3.3幂级数解法幂级数解法待定的是级数的系数，因而通常计算较大，其实幂级数解法适用二阶以上的高阶齐次线性微分方程与非齐次线性微分方程，也能求其特解或通解. 二阶线性方程.在近代物理学以及工程技术中有着很广泛的应用，其中幂级数解法不但对于求解方程有意义，而且还由此引出了很多新的超越函数，在理论上是很重要的.下来给出两个定理，若要了解定理证明过程，可参考有关书籍.定理1 如果在某点的邻域内解析，即它们可展成的幂级数，且，则的解在的邻域内也能展开成为的幂级数.定理2 如果在的邻域内解析，而为的重零点，是的不低于重的零点，（若），是的不低于重的零点，（若），则方程至少有一个形如的广义幂级数解，其中r为某一实数. 若要了解幂级数的详细解法可以参考常微分方程，这里不做具体分析. 总之，不同的方法用于不同类型的方程，这是应用之时必须特别注意之点.参 考 文 献1 朱思铭，王寿松等.常微分方程m.北京：高等教育出版社.2006（3）:126-129.2 汤光宋，余复民.应用交换变量位置法解两类一阶常微分方程j.兰州工业 高等专科学校学报.1996,(1):20-25. 3 焦洪田.一阶非线性微分方程的常数变易法j.雁北师范学院学报.1999（6）：44-45. 4 周斌.常数变易法在数学分析中的应用j.内江师范学院学报.2003,18（4）:56-58.5 曹玉平.一阶线性变系数微分方程组的矩阵解法j.河北理工学院学报.2005（2）. 6 阮炯.差分方程和常微分方程m.上海：复旦大学出版社.2007:39-90.7 黄雪燕.常微分方程的化归思想j.长春师范学院学报.2007,26(04):24-26. 8 李鸿祥.对“关于的特解”一文的意见j.高等数学研究. 2007，9（03）:55.9 刘林平.常系数线性微分方程的拉普拉斯变换解法j.内蒙古农业大学学报.2006,27（04）:157-159.10 黄启昌，任永泰,陈秀东等.常微分方程m.东北师范大学数学系微分方程 教研室编.人民教育出版社.2008:173-180. ordinary differential equation of elementary method of classificationguo zhen -zhenabstract: describing several elementary solutions to solve ordinary differential equations, such as the constant variation, integrating factor method, la