问题二阶线性常微分方程边值问题的数值解法.docx

1、问题二阶线性常微分方程边值问题的数值解法【关键字】问题重庆三峡学院毕业论文论文题目：二阶线性常微分方程边值问题的数值解法专 业：数学与应用数学年 级： 2004级学 号： 0203作 者： 指导老师： 查中伟（教授）完成时间：2008年5月二阶线性常微分方程边值问题的数值解法摘 要 对于二阶线性常微分方程定解问题，大多数是不存在解析解的，有的方程既使存在解析解，然而解出其解也是很难办到的.尤其是工程计算中更需要的是数值解.因此,本文给出两种二阶线性常微分方程边值问题的数值解法.首先给出了利用差分法求其数值解，在介绍此方法的过程中，第一步构造了二阶线性常微分方程边值问题的差分格式（差分方程组），

2、然后论证了该边值问题的收敛性，最后利用二阶微分的四阶紧致差分公式对第一步构造的差分格式进行精度上的改进，得到了较好的结果.接着介绍了二阶线性常微分方程边值问题的第二种数值解法Taylor展开解法，该方法主要是先将边值问题转化为Fredholm积分方程，再经过数学处理即可得到关于近似解、近似解的一阶导数和近似解的二阶导数的线性方程组，最后利用Crammer法则解出了该二阶线性常微分方程边值问题的数值解.并且利用工程数学软件MATLAB，给出了计算机程序，使前面两种算法在计算机上得以实现.最后给出了具体实例，分别运用以上两种解法进行求解，对这两种方法的计算精度进行了对比分析.关键词 数值解；差分格

3、式；解的收敛性；MATLABNumerical Solution for Boundary Value Problems of the Second-order Linear ordinary Differential EquationsKai-min Cheng(Grade 2004, mathematics and applied mathematics, School of Mathematics and Computer Science, , Chongqing Wanzhou 404000)Abstract: There is no exact solution for the ma

4、jority of second-order linear ordinary differential equations solution problem. Some of them even exists the exact solution, but to solve its solution is a very difficult job. Especially, we need the numerical solution urgently in Engineering Mathematics. Based on this, this paper gives two kinds of

5、 numerical solution for Boundary Value Problems of the second-order linear ordinary differential equations. Firstly, this paper gives the difference method to solve its numerical solution. In the process of introducing this method, we construct differential format of the Boundary Value Problems in t

6、he first step. Then we demonstrate the convergence of the Boundary Value Problem. Finally, we improve the accuracy for the difference format constructed in first step by the four-order differential format Secondly, this paper gives method of expansion to solve its numerical solution. In this method,

7、 we first transform the boundary value problems into Fredholm integral equation, and then can get a group of linear equations with unknowns to the approximate solution, the first order derivative of the approximate solution and the second derivative of the approximate solution after mathematical tre

8、atment, and can solve it in Crammer rule. Thirdly, we write an algorithm program by using engineering mathematical software and make above two methods realized on the computer. Finally, this paper gives a specific example and solves it with above two methods separately. This paper also analyzes the

9、feasibility for their accuracy.Keywords: Numerical solutions; Differential format; The convergence of solutions; Matlab language0 引言当前，常微分方程的定解问题已经有很多重要结果，如解的存在性定理.在很多典型的常微分方程的解法上也有较大突破，同时也涌现出了一批较为经典的解法，如降阶法、积分变换法、变易常数法、特征方程法等方法.尽管如此，在数学领域中还存在着迄今为止还难以解出其解析解的微分方程，这就使得微分方程领域必然会产生一个新的微分方程分支微分方程数值解法.对于较

10、为简单的常微分方程，只需利用经典解法即可解出其解析解，边值问题也是如此.往往在实际工程中抽象出来的微分方程，其边值问题是相当复杂的，所以用求其解析解的方法来计算微分方程的边值问题往往是不适宜的，有时甚至是很难办到的.实际上，对于解微分方程，我们所要获取的或感兴趣的，往往只是一个或几个特定点上的数据，并且既使有的方程存在解析解，也并不意味着其一定能够表达成初等函数形式，如多项式函数、对数函数、指数函数、三角函数、反三角函数以及它们的有限组合形式.在利用数值计算方法理论的基础上，再辅之计算机若能很好的解决现实工程数学中的许多常微分方程的边值问题，这是非常好的.本文主要论述二阶线性常微分方程边值问题

11、数值解的两种解法：差分法和Taylor展开法.考虑二阶线性常微分方程的边值问题 (1) 1 差分法1.1 差分格式的构造对于边值问题(1)将区间进等距划分，分点为：其中 称与为边界点，称为内部节点.在每个节点上将用差商近似表示.这里要求有相同的截断误差，以保证精度协调.对于内部节点，二阶导数用二阶中心差商表示，得 ， (2)一阶导数用一阶中心差商表示，得 ， (3)假设，则；于是方程(1)的差分方程： ， (4)其中，.称(4)为(1)的中心差分格式.1.2 差分格式(4)的收敛性作适当变换可以消除线性方程(1)中的一阶导数项.事实上，可令，再作，将之代入(1)得所以，不妨仅就缺少一阶导数项的

12、方程来讨论.对于边值问题： (5)这里假定，其对应的差分问题是： (6)下面讨论差分问题(6)的可解性.由于(6)式是关于变量的线性方程组，要证明它的解存在唯一，只要证明对应的齐次方程组只有零解.为此，引进文献5的极值原理：引理1(极值原理) 对于一组不全相等的数，记其中 如果，则的正的最大值只能是或；如果，则的负的最小值只能是或.证明 用反证法 考察的情形，设是正的最大值，即，且和中至少有一个小于，此时有:由于，故由上式推出，此与原假设矛盾.此外，的情形可类似地进行讨论，证毕.利用引理1的结论有：定理1 差分问题(6)的解存在并且是唯一的.证明 只要证明对应的其次方程组只有零解，由于这里，由

13、极值原理知，的正的最大值和负的最小值只能是或，而按边界条件，故所有全为零，证毕.下面运用极值原理论证差分方法的收敛性并估计误差.定理2 设是差分问题(6)的解，而是边值问题(6)的解在节点处的值.则截断误差有估计式： (7)式中 证明 由Taylor展开式，易得 (8)将(8)与(6)相减，知误差满足 (9)式中的一般不知道的，讨论下列差分问题 (9-1)式中.首先证明(9)和(9-1)两个差分问题的解存在下列关系： (10)事实上，由于，故有又 利用引理1知，即(10) 式成立.我们进一步考察 (11)这里，又，故由引理1（注意时就是），知，即，于是有 然而是容易求出的. 事实上，可以先求解

14、差分问题(9)所对应的边值问题得 容易验证是差分问题(11)的解，注意到在点达到最大值因此有估计式(7)，证毕.根据估计式(7)知，当时有，这表明差分问题(6)是收敛的.又因为(4)式是含有个未知数的线性方程组，方程的个数，要使方程组(4)有唯一解，还需要有两个边值条件，它们和(4)一起构成三对角方程组.通过分解，采用追赶法即可解出(4)(见3.1)1.3 差分格式(4)的改进从1.2的讨论可以得知，在1.1构造的差分格式是直接的中心差分格式，其截断误差是即有二阶精度.基于中心差分格式的分析，直接利用二阶微分的四阶紧致差分公式，我们得到了数值求解二阶线性常微分方程边值问题的一种四阶精度的差分格

15、式.1.3.1 改进后的差分格式的推导为了便于推导，可将方程(1)改写为 (12)为了使格式具有更高得精度，利用11中的四阶差分公式 (13)其中，为二阶中心差分算子 (14)可以表示、及等.将(13)代入(12)得: (15)即 (16)而(2)即为 (17)联立(16)(17)知 (18)即 (19)显然要使(19)具有四阶精度，必须对也进行四阶离散，于是利用Taylor公式知 (20)又由(12)可得 (21)于是有 (22)再将(2)、(3)、(21)及(22)代入(19)，经整理并略去高阶项可得 (23)其中 (24) (25) (26) (27)故差分方程(23)即为所要推导的四阶精度格式，其截断误差为.1.3.2 差分格式(23)稳定性分析现对改进后的差分格式(23)进行稳定性分析.首先引用文献12对差分方程的稳定性的定义.定义1(差分算子的稳定性定义) 如果对于充分小的网格步长，线性差分算子对任何离散函数均存在不依赖于的正常数，使得 (28)则称差分算子是稳定的. 为了论证差分格式(23)的稳定性，先引入文献13的一个引理.引理2 假设是正型差分算子，