二阶常微分方程边值问题的数值解法毕业论文(36页).doc

《二阶常微分方程边值问题的数值解法毕业论文(36页).doc》由会员分享，可在线阅读，更多相关《二阶常微分方程边值问题的数值解法毕业论文(36页).doc（36页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、-二阶常微分方程边值问题的数值解法毕业论文-第 29 页二阶常微分方程边值问题的数值解法摘 要 求解微分方程数值解的方法是多种多样的，它本身已形成一个独立的研究方向，其要点是对微分方程定解问题进行离散化本文以研究二阶常微分方程边值问题的数值解法为目标，综合所学相关知识和二阶常微分方程的相关理论，通过对此类方程的数值解法的研究，系统的复习并进一步加深对二阶常微分方成的数值解法的理解，为下一步更加深入的学习和研究奠定基础对于二阶常微分方程的边值问题，我们总结了两种常用的数值方法：打靶法和有限差分法在本文中我们主要探讨关于有限差分法的数值解法构造差分格式主要有两种途径：基于数值积分的构造方法和基于T

2、aylor展开的构造方法后一种更为灵活，它在构造差分格式的同时还可以得到关于截断误差的估计在本文中对差分方法列出了详细的计算步骤和Matlab程序代码，通过具体的算例对这种方法的优缺点进行了细致的比较在第一章中，本文将系统地介绍二阶常微分方程和差分法的一些背景材料在第二章中，本文将通过Taylor展开分别求得二阶常微分方程边值问题数值解的差分格式在第三章中，在第二章的基础上利用Matlab求解具体算例，并进行误差分析关键词：常微分方程,边值问题，差分法，Taylor展开，数值解The Numerical Solutions ofSecond-Order Ordinary Differentia

3、l Equationswith the Boundary Value ProblemsABSTRACTThe numerical solutions for solving differential equations are various. It formed an independent research branch. The key point is the discretization of the definite solution problems of differential equations. The goal of this paper is the numerica

4、l methods for solving second-order ordinary differential equations with the boundary value problems. This paper introduces the mathematics knowledge with the theory of finite difference. Through solving the problems, reviewing what have been learned systematically and understanding the ideas and met

5、hods of the finite difference method in a deeper layer, we can establish a foundation for the future learning.For the second-order ordinary differential equations with the boundary value problems, we review two kinds of numerical methods commonly used for linear boundary value problems, i.e. shootin

6、g method and finite difference method. There are mainly two ways to create these finite difference methods: i.e. Taylor series expansion method and Numerical Integration. The later one is more flexible, because at the same time it can get the estimates of the truncation errors. We give the exact cal

7、culating steps and Matlab codes. Moreover, we compare the advantages and disadvantages in detail of these two methods through a specific numerical example. In the first chapter, we will introduce some backgrounds of the ordinary differential equations and the difference method. In the second chapter

8、, we will obtain difference schemes of the numerical solutions of the Second-Order ordinary differential equations with the boundary value problems through the Taylor expansion. In the third chapter, we using Matlab to solve the specific examples on the basis of the second chapter, and analyzing the

9、 errors. KEY WORDS: Ordinary Differential Equations, Boundary Value Problems, Finite Difference Method, Taylor Expansion, Numerical Solution毕业论文（设计）原创性声明本人所呈交的毕业论文（设计）是我在导师的指导下进行的研究工作及取得的研究成果。据我所知，除文中已经注明引用的内容外，本论文（设计）不包含其他个人已经发表或撰写过的研究成果。对本论文（设计）的研究做出重要贡献的个人和集体，均已在文中作了明确说明并表示谢意。 作者签名： 日期： 毕业论文（设计）授

10、权使用说明本论文（设计）作者完全了解*学院有关保留、使用毕业论文（设计）的规定，学校有权保留论文（设计）并向相关部门送交论文（设计）的电子版和纸质版。有权将论文（设计）用于非赢利目的的少量复制并允许论文（设计）进入学校图书馆被查阅。学校可以公布论文（设计）的全部或部分内容。保密的论文（设计）在解密后适用本规定。 作者签名： 指导教师签名： 日期： 日期： 注 意 事 项1.设计（论文）的内容包括：1）封面（按教务处制定的标准封面格式制作）2）原创性声明3）中文摘要（300字左右）、关键词4）外文摘要、关键词 5）目次页（附件不统一编入）6）论文主体部分：引言（或绪论）、正文、结论7）参考文献8

11、）致谢9）附录（对论文支持必要时）2.论文字数要求：理工类设计（论文）正文字数不少于1万字（不包括图纸、程序清单等），文科类论文正文字数不少于1.2万字。3.附件包括：任务书、开题报告、外文译文、译文原文（复印件）。4.文字、图表要求：1）文字通顺，语言流畅，书写字迹工整，打印字体及大小符合要求，无错别字，不准请他人代写2）工程设计类题目的图纸，要求部分用尺规绘制，部分用计算机绘制，所有图纸应符合国家技术标准规范。图表整洁，布局合理，文字注释必须使用工程字书写，不准用徒手画3）毕业论文须用A4单面打印，论文50页以上的双面打印4）图表应绘制于无格子的页面上5）软件工程类课题应有程序清单，并提供

12、电子文档5.装订顺序1）设计（论文）2）附件：按照任务书、开题报告、外文译文、译文原文（复印件）次序装订3）其它目 录前 言1第一章 二阶常微分方程2第二章 边值问题的数值解法72.1 有限差分逼近的相关概念72.2差分方程的建立82.3差分问题的可解性102.4差分方程的收敛性112.5 差分方程的稳定性122.6有限差分分方程的解法12第三章 具体算例15 3.1 二阶常微分方程算例的数值解15 3.2 算例结果分析25总 结26参考文献27附 录30前 言求解常微分方程边值问题在计算数学领域中一直占很重要的地位，但是常微分方程中仅有一些典型的方程能求出解析解，大部分是求不出解析解的因此常

13、微分方程数值解的研究具有重要的现实意义 用数值方法求解微分方程问题几乎是与微分方程同时出现的，其历史可以追溯到月一个半世纪前上个世纪中叶以后，由于微分方程本身的理论的深入发展，兼之电子计算机的问世，用数值方法求解微分方程问题的研究更进入了一个蓬勃发展的新局面求解常微分边值问题最有效的方法之一是有限差分法经典的有限差分法是利用差商代替导数（数值微分）或者差分插值（数值积分）的方法来构造差分格式为了构造具有较高截断误差的差分格式，近年来一些学者提出了利用样条函数或者参数样条函数的方法来近似代替未知函数通过配置的方法，构造出一些样条差分格式，但高阶数值微分公式和关于高次样条函数的高阶导数的计算都较为

14、困难，同时构造差分格式引起的计算量非常大，有的方法精度并不高，所以这些方法都不能很好地适应高阶微分方程 本文就二阶常微分方程边值问题，利用差分法求解数值解有限差分法是数值方法中最经典的方法这种方法发展较早，比较成熟，较多用于求解双曲型和抛物型问题用有限差分法近似求解常微分方程问题有多种多样的方法，并且也可以用不同的构造方法来建立这些有限差分法用Taylor级数展开方法是最常用的方法用Taylor法展开来建立差分格式，实际上等价于用差商来近似微商得到相应的差分格式第一章 二阶常微分方程 如果在一个微分方程中出现的未知函数只含一个自变量，这个方程就叫做常微分方程，也可以简单地叫做微分方程1 微分方

15、程差不多是和微积分同时先后产生的，苏格兰数学家耐普尔创立对数的时候，就讨论过微分方程的近似解牛顿在建立微积分的同时，对简单的微分方程用级数来求解后来瑞士数学家雅各布贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论常微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对常微分方程的发展产生了深刻的影响，当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具对常微分方程的研究1可分为以下几个阶段： 发展初期是对具体的常微分方程希望能用初等函数或超越函数表示其解，

16、属于“求通解”时代 莱布尼茨(Leibniz)曾专门研究利用变量变换解决一阶微分方程的求解问题，而欧拉(Euler)则试图用积分因子统一处理伯努利(Bernoulli)、里卡蒂（Riccati）微分方程就是在研究初等积分时提出后人以他们的名字命名的方程 早期的微分方程的求解热潮被刘维尔（Liouvile)于1841年证明里卡蒂方程不存在一般的初等解而中断 加上柯西（Cauchy)初值问题的提出，常微分方程从“求通解”转向“求定解”时代 首先是对常微分方程定解问题包括初值和边值问题的解的存在性2、唯一性等解的性质的研究 其次，针对线性微分方程，特别是二阶线性微分方程，通过专门定义一些特殊函数以求

17、解特殊方程，如贝塞尔（Bessel）函数、勒让德(Legendre)多项式等，这促成了微分方程与（复变）函数论结合产生微分方程解析理论 同时，由于天文计算的需要促进了常微分方程摄动理论以及小参数、幂级数等近似方法的研究 19世纪末，天体力学中的太阳系稳定性问题需要研究常微分方程解的大范围性态，从而使常微分方程的研究从“求定解问题”转向“求所有解”的新时代 首先，庞加莱（Poincare）创立了定性理论和方法研究常微分方程解的大范围性态由于希尔伯特（Hilbert）提出20世纪23个数学问题中关于极限环个数的第16问题，大大促进了定性理论的发展 另一方面李雅普诺夫（Lyapunov）提出的运动稳

18、定性理论，用于解决方程解的初值不影响原方程解的趋向问题，在天文、物理及工程技术中得到广泛应用，先后在前苏联、美国受到极大重视 同时，伯克霍夫（Birkhoff）在20世纪初在动力系统方面开辟了一个新领域，由于拓扑方法的渗入，20世纪50年代后经阿诺德（Arnold）、斯梅尔（Smale)等大数学家的参与而得到蓬勃发展 除定性、稳定性和动力系统理论外，还有非线性振动理论、摄动和奇异摄动理论及变换群理论在20世纪也得到了迅速的发展 20世纪六七十年代以后，常微分方程由于计算机技术的发展迎来了新的时期，从“求所有解”转入“求特殊解”时代，发现了具有新性质的解和方程，如混沌（解）、奇异吸引子及孤立子等

19、科技和数学界的重大发现是混沌、孤立子和分形，其中混沌，孤立子直接与微分方程有关洛伦茨在20世纪60年代发现了称为Lorenz方程的常微分方程，初始敏感的特征导致了混沌现象的发现引起了科学界的巨大震动，斯梅尔称之为“利用牛顿的定律推翻了牛顿决定论”孤立子本是物理上有重要意义的偏微分方程的新类型解，但它们往往对应于可积的哈密顿系统的常微分方程，从而引发了对停顿百年的常微分方程可积性的研究热潮常微分方程的研究领域其它学科或领域的结合而出现各种新的研究分支，如控制论、种群生态学、分支理论、泛函微分方程、脉冲微分方程、广义微分方程、时标微分方程等 常微分方程3属于数学分析的一支，是数学中与应用密切相关的

20、基础学科，其本身也在不断发展中，学好常微分方程基本理论与方法对进一步学习研究数学理论和实际应用非常重要 而常微分方程边值问题则是微分方程理论研究的一个基本问题，也是最为重要的课题之一 它在应用科学和工程领域有着非常重要的作用，例如工程学、力学、天文学、经济学以及生物学等领域中的许多实际问题通常会归结为常微分方程边值问题的求解 虽然常微分方程问题有许多解析方法可以求解，但这些方法只能求解一些特殊类型的方程，对从实际问题中提炼出来的微分方程往往不再适用，因而对常微分方程边值问题的数值方法的研究就显得尤为重要 经典的数值方法有: 试射法（打靶法）和有限差分法4而用打靶法求解线性问题时，解的精度较高，

21、这是因为打靶法将边值问题的求解方法转化为相应的初值问题的求解，因而可以使用具有较高精度的Runge_Kuta法，但是算法的稳定性较差常微分方程边值问题已被深入而广泛地研究，并取得了系统而深刻的结果科学和工程技术中有许多实际问题都可以转化为微分方程的求解问题，而大量的微分方程很难求出其解析解，因此，微分方程的数值解法的研究就显得具有重要意义边值问题主要研究微分方程的求解及解的性质，其补充条件由以自变量取某些值时未知函数及其导数的值而定，许多数学和物理问题都可以归结为微分方程边值问题解这类问题的基本思想和初等数学解方程的基本思想很相似，也就是要把研究的问题中已知函数和未知函数之间的关系找出来，从列

22、出的包含未知函数的一个或几个方程中去求得未知函数的表达式但是无论在方程的形式、求解的具体方法、求出解的性质等方面，都和初等数学中的解方程有许多不同的地方 常微分方程的概念、解法、和其它理论很多，比如，方程和方程组的种类及解法、解的存在性和唯一性、奇解、定性理论等等下面就方程解的有关几点简述一下，以了解常微分方程的特点求通解在历史上曾作为微分方程的主要目标，一旦求出通解的表达式，就容易从中得到问题所需要的特解也可以由通解的表达式，了解对某些参数的依赖情况，便于参数取值适宜，使它对应的解具有所需要的性能，还有助于进行关于解的其他研究 后来的发展表明，能够求出通解的情况不多，在实际应用中所需要的多是

23、求满足某种指定条件的特解当然，通解是有助于研究解的属性的，但是人们已把研究重点转移到定解问题上来 一个常微分方程是不是有特解呢？如果有，又有几个呢？这是微分方程论中一个基本的问题，数学家把它归纳成基本定理，叫做存在和唯一性定理5因为如果没有解，而我们要去求解，那是没有意义的；如果有解而又不是唯一的，那又不好确定因此，存在和唯一性定理对于微分方程的求解是十分重要的 大部分的常微分方程求不出十分精确的解，而只能得到近似解当然，这个近似解的精确程度是比较高的另外还应该指出，用来描述物理过程的微分方程，以及由试验测定的初始条件也是近似的，这种近似之间的影响和变化还必须在理论上加以解决 边值问题及其求解

24、发展过程的介绍： 一般地说，n 阶微分方程的解含有 n个任意常数也就是说，微分方程的解中含有任意常数的个数和方程的解数相同，这种解叫做微分方程的通解通解构成一个函数族 如果根据实际问题要求出其中满足某种指定条件的解来，那么求这种解的问题叫做定解问题，对于一个常微分方程的满足定解条件的解叫做特解对于高阶微分方程可以引入新的未知函数，把它化为多个一阶微分方程组微分方程是联系着自变量、未知函数及其导数的关系式如果在微分方程中，自变量的个数只有一个，我们称这种微分方程为常微分方程5一般的n阶常微分方程具有形式那么，二阶常微分方程的形式为: 牛顿在研究天体力学和机械力学的时候，利用了微分方程这个工具，从

25、理论上得到了行星运动规律后来,法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量 微分方程的理论逐步完善的时候，利用它就可以精确地表述事物变化所遵循的基本规律，只要列出相应的微分方程，有了解方程的方法微分方程也就成了最有生命力的数学分支 现在，常微分方程在很多学科领域内有着重要的应用6，自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等这些问题都可以化为求常微分方程的解，或者化为研究解的性质的问题应该说，应用常微分方程理论已经取得了很大的成就，但

26、是，它的现有理论也还远远不能满足需要，还有待于进一步的发展，使这门学科的理论更加完善. 第二章 边值问题的数值解法 有限差分法是用于微分方程定解问题最广泛的数值方法，其基本思想是以差商近似代替导数，把微分方程化离散为差分方程组，并把相应的解作为微分方程定解条件的近似解.在本节我们讨论有关边值的问题，介绍两点边值问题的一种数值解法有限差分法7.2.1 有限差分逼近的相关概念设函数光滑，且，利用Taylor展开，可得 (2-1-1) (2-1-2)由(2-1-1)可以得到一阶导数的表达式 (2-1-3) 或者 (2-1-4)其中表示截断误差8项.因此，可得一阶导数的的差分近似表达式为 (2-1-5

27、) (2-1-6)由上述式子可知，差商(2-1-5)和(2-1-6)逼近微商为一阶，即为，为了得到更为精确地差分表达式，将(2-1-1)减(2-1-2)可得 (2-1-7)从而可以得到 (2-1-8)或者 (2-1-9)其中 (2-1-10)由此可知，差商逼近微商的精度为二阶，即为. (2-1-5)、(2-1-6)与(2-1-10)公式分别被称为逼近一阶微商的向前，向后和中心差分公式. 类似地，我们还可以给出二阶微商和高阶微商的差分近似表达式.例如将(2-1-1)和(2-1-2)相加可得 (2-1-10a)进而有 (2-1-10b)其中. 因此，二阶导数的差分近似表达式为 (2-1-10)2.

28、2差分方程的建立在具体求解微分方程时，必须附加某种定解条件.微分方程和定解条件一起组成定解问题，对于二阶微分方程通常有两种给法：一种是给出了积分曲线在初始时刻的性态，这类条件称为初始条件，相应的定解问题称为初值问题8；另一种是给出了积分曲线首末两端的性态，这类条件称为边界条件，相应的定解问题称为边值问题.下面我们将要讨论的就是边值问题的数值解法9.首先给出二阶常微分方程的及其边界条件，如下式： (2-2-1)下面将应用差分方法来解决这个问题.差分方法的关键，在于恰当的选取差商逼近微分方程中的导数，我们知道，逼近一阶导数可用向前差商，也可用向后差商或者中心差商.中心差商是向前差商和向后差商的算术

29、平均.为逼近二阶导数，一般用二阶差商向前差商的向后差商(即向后差商的向前差商) (2-2-2) 设将积分区间划分为等分，步长，节点.差商替代相应的导数，可将边值问题(2-2-1)离散化得到下面的公式： (2-2-3)如果函数是非线性的，那么所归结出的差分方程也是非线性的，这时实际求解比较困难.如果所给方程(2-2-1)是如下形式的线性方程： (2-2-4) 则差分方程(2-2-2)相应的形式为 (2-2-5) 其中的下标表示在节点的取值.利用边界条件(2-2-3)消除式(2-2-5)中的和，整理得到关于的下列方程组： (2-2-6)这样归结出的方程组是所谓的三对角形的，即： (2-2-7)因为

30、它的稀疏矩阵仅在注对角线及其相邻的两条对角线上有非零元素.求解这种三对角形方程组，用所谓追赶法10特别有效.在后面我们将介绍.2.3差分问题的可解性 我们知道，通过自变量的适当变换可以消除线性方程(2-2-4)中的一阶导数项.下面仅就缺一阶导数项的方程来讨论这一问题.考察边值问题： (2-3-1)这里假定，其对应的差分问题是 (2-3-2)现在论证差分问题的可解性.由于式(2-3-2)是关于变量的线性方程组.要证明它的解的存在唯一，只要证明对应的齐次方程组只有零解.为此，我们要引进下述极值定理.定理1：对于一组不全相等的数 ，记 (2-3-3)假定 ，则的正的最大值只能是或；如果 ，则的负的最

31、小值只能是或.定理2：差分问题式(2-3-2)的解存在的并且唯一的.证明略，见参考文献11.2.4差分方程的收敛性对于任给如果数值解当(同时)时趋向准确解，则称该差分方程是收敛的.现在运用极值原理证明差分方程的收敛性并估计误差.定理3 设是差分问题(2-3-2)的解，而是边值问题(2-3-1)的解在节点的值，则截断误差有下列估计式： (2-4-1)证明略，可参见参考文献12.2.5 差分方程的稳定性 前面关于收敛性问题的讨论有个前提，必须假定差分方程的每一步计算都是准确的.事实情况并不是这样，差分方程的求解还会有计算误差，譬如由于数字舍入而引起的扰动.这类扰动在传播过程中会不会恶性增长，以至于

32、“淹没”了差分方程的“真解”，这就是差分方程的稳定性13问题.如果一种差分方法在节点值上大小为的扰动，导致以后各节点值 上产生的偏差均不超过，则称该方法是稳定的.在实际计算时，希望某一步的扰动在后面的计算中能够被控制，甚至是逐渐衰减的.也就是希望所应用的方法具有稳定性.边值问题的病态性(稳定性)与解得模式和边界条件之间的相互作用有关.对于良态的边值问题，解得增长和衰变模式要受到约束条件的控制.对于这些问题的研究在这里不作详细讨论.2.6有限差分分方程的解法下面我们来研究以下这些差分方程组的解法：将系数矩阵为对角占优的三对角方程组简记为 (2-6-1)其中 (2-6-2)其中，满足下列对角占优条

33、件：(1) (2-6-3)(3) 下面利用矩阵的直接三角分解法来推导求解三对角方程组(2-6-1)的计算公式.由系数矩阵的特点，可以将分解为两个三角阵的乘积，即其中为下三角阵，为单位上三角阵.下面说明这种分解是可能的 (2-6-4) 其中为待定系数.比较(2-5-4)两端即得 (2-6-5)由，得.下面用归纳法证明 (2-6-6)即从而由式(2-6-5)可求出.式(2-6-6)对是成立的.现设式(2-6-6)对成立，求证对亦成立.由归纳法假设，又由式(2-6-6)及的假设条件，有：也就是，由式(2-6-5)得到也就是说，由得假设条件完全确定了，实现了与的分解. 那么，求解等价于解两个三角方程组

34、与，先后求出与，从而得到以下解三对角方程组的追赶式公式14：步1 计算的递推公式步2 解： ；步3 解： .将计算系数及的过程称为追的过程，将计算方程组的解的过程称为赶的过程.追赶法公式实际上就是把Gauss消去法用到求解三对角方程组上的结果.这时由于特别简单，因此使得求解的计算公式也非常简单，计算量较小.第三章 具体算例 3.1 二阶常微分方程算例的数值解例1：应用差分格式计算边值问题 (3-1)该定解问题的精确解为 由下面三个表格中可以看出，当步长变小时，误差也随之变小当步长为0.01时，误差的精度达到当步长缩小为原来的一半时，误差精度变为而当步长变为0.002时，误差精度就达到下面的表格

35、依次表示当步长分别取0.01、0.005及0.002时部分节点的数值解，精确解和误差值。x 数值解 精确解|数值解-精确解|0.1 0.2 0.3 0.40.50.60.70.80.91.0100521064163231.0408137039811691.0941772788901941.1735130190552281.2840258452347541.4333273942110011.6323114112108441.8964738519904672.2479012517352121.010050167084168 1.040810774192388 1.094174283705210 1

36、.173510870991810 1.284025416687741 1.433329414560340 1.632316219955379 1.896480879304952 2.2479079866764720.19393321546434810.29297887811363670.2995184983589638 0.2148063417362067 0.4285470129872238 0.20203493396042660.48087445352340550.7027314484320257 0.6734941259800564表3-1 算例3-1部分节点处数值解、精确解和误差的绝对

37、值（）x数值解精确解|精确解-数值解|0.10.20.30.40.50.60.70.80.91.0025034002819671.0100506518584101.022755671856521 1.0408115065211241.064495228119752 1.0941750323230961.130319790985197 1.1735114077723231.2244604328788051.0100501670841681.0408107741923881.0941742837052101.1735108709918101.2840254166877411.43332941456

38、03401.6323162199553791.8964808793049522.2479079866764720.2726761714200170.4847742423041980.6376920747008800.7323287360971880.7692018919591700.7486178856819240.6709111863845150.5367805127054480.347756889507167表3-2 算例3-1部分节点处数值解、精确解和误差的绝对值（）x数值解精确解|数值解-精确解|0.10.20.30.40.50.60.70.80.91.0004000986425581

39、.0016013163685441.0036065389636691.0064205888949161.0100502446453321.0145042679347701.0197934389791891.0259305999878251.0329307071432121.0004000800106681.0016012806829401.0036064877830041.0064205237606611.010050167084168 1.014504179460763 1.0197933410959711.0259304941903821.0329305949196555.35475607

40、691543 1.06855014619586 1.59922919204178 2.12751829712504 2.65342295618226 3.17694863385259 3.69810076433080 4.21688472335147 4.73330589757792表3-3 算例3-1部分节点处数值解、精确解和误差的绝对值（）图3-1给出了精确解曲线和步长为0.1的数值解曲线图3-6给出了取不同步长的误差图曲线图3-1 精确解曲线图3-2 步长0.01数值解曲线图3-3 步长为时的误差曲线图3-4步长为时的误差曲线图3-5 步长为时的误差曲线图3-6 不同步长时的误差曲线例2

41、 考虑下面方程 (3-2)解：可解出方程的精确解即解析解为，对方程进行Matlab编程，下面给出部分节点的数值解，精确解及误差图3-4给出了精确解曲线和步长为0.1的数值解曲线图3-12给出了取不同步长的误差图曲线图3-5，3-6给出了步长为0.005，0.002，的数值解曲线. 在图3-7中表示出该方程精确解的曲线. 在图3-8中画出了步长为0.01时的数值解曲线. 图3-9到图3-11分别表示了步长为0.01，0.005，0.002时的误差曲线图 . x 数值解 精确解|数值解-精确解|0.1 0.20.30.40.5 0.60.70.80.90.110523074942274 0.244

42、291156225608 0.404971599857503 0.596745965162951 0.824377623797421 1.093287899469460 1.4096417854116101.780444412427270 2.2136495810241000.1105170918075650.2442805516320340.4049576422728010.596729879056508 0.8243606353500641.0932712802343101.4096268952293301.780432742793970 2.2136428000412600.598313470928651.06045935741594 1.39575847021289 1.60861064432138 1.698844735653631.66192351596095 1.48901822767389 1.166963329768400.67809828405707表3-4 算例3