二阶常微分方程边值问题的数值解法讲义.docVIP

摘 要 本文主要研究二阶常微分方程边值问题的数值解法。对线性边值问题，我们总结了两类常用的数值方法，即打靶法和有限差分方法，对每种方法都列出了详细的计算步骤和Matlab程序代码，通过具体的算例对这两类方法的优缺点进行了细致的比较。 关键字：常微分方程边值问题；打靶法；差分法； ABSTRACT This article mainly discusses the numerical methods for solving Second-Order boundary value problems for Ordinary Differential Equations. On the one hand, we review two types of commonly used numerical methods for linear boundary value problems, i.e. shooting method and finite difference method. For each method, we give both the exact calculating steps , we compare the advantages and disadvantages in detail of these two methods through a specific numerical example. Key words：Boundary-Value Problems for Ordinary Differential Equations；Shooting Method；Finite Difference Method； 目 录 第一章 引言 - 1 - 第二章 二阶线性常微分方程 - 3 - 2.1试射法（“打靶”法） - 4 - 2.1.1简单的试射法 - 4 - 2.1.2 基于叠加原理的试射法 - 5 - 2.2 有限差分法 - 11 - 2.2.1 有限差分逼近的相关概念 - 12 - 2.2.2 有限差分方程的建立 - 14 - 2.2.3 其他边值条件的有限差分方程 - 15 - 2.2.4 有限差分方程的解法 - 17 - 第三章 二阶非线性微分方程 - 22 - 3.1基于牛顿迭代法的打靶法 - 22 - 3.1.1 第一类边值条件推导 - 22 - 3.1.2 其他边值条件的推导 - 24 - 3.1.3 算法及程序代码 - 25 - 3.2 基于改进的牛顿迭代法的打靶法 - 31 - 3.2.1 算法的推导 - 31 - 3.2.2 算法及代码 - 32 - 第四章 改进算法的算例 - 38 - 第五章 总结 - 45 - 参考文献 - 46 - 致谢 - 47 - 第一章 引言 微分方程是现代数学中一个很重要的分支，从早期的微积分时代起，这个学科就成为了理论研究和实践应用的一个重要领域。在微分方程理论中，定解条件通常有两种提法：一种是给出了积分曲线在初始时刻的性态，相应的定解条件称为初值问题；另一种是给出了积分曲线首末两端的性态，这类条件则称为边界条件，相应的定解问题称为边值问题。 常微分方程边值问题在应用科学与工程技术中有着非常重要的应用，例如工程学、力学、天文学、经济学以及生物学等领域中的许多实际问题通常会归结为常微分方程边值问题[12]的求解。文献[9]给出了边值问题求解的方法，虽然求解常微分方程边值问题有很多解析方法可以求解，但这些方法只能用来求解一些特殊类型的方程，对从实际问题中提炼出来的微分方程往往不再适用，因而对常微分方程边值问题的数值方法的研究显得尤为重要。经典的数值方法主要有：试射法（打靶法）和有限差分法，见文献[2]。对于二阶线性边值问题，差分法的优点在于稳定性较好，但它的精度不高。而用打靶法求解线性问题时，解的精度较高，这是因为打靶法将边值问题的求解转化为相应的初值问题的求解，因而可以使用具有较高精度的Runge-Kutta法（见文献[1]），但是算法稳定性较差。 在本文中，我们首先总结了二阶线性边值问题的数值算法：打靶法、有限差分法。对每种方法都列出了详细的计算步骤和Matlab程序代码，通过具体的算例对这两类方法的优缺点进行了细致的比较。由于简单的打靶法过分依赖经验，我们考虑了基于线性叠加原理的打靶法，将线性边值问题转化为两个初值问题，并通过线性叠加得到原边值问题的解。 第二章 二阶线性常微分方程 二阶常微分方程一般可表示成如下的形式： ， (1) 边值条件有如下三类[9]： 第一类边值条件 ， (2) 第二类边