[理学]第九章常微分方程初值问题数值解法.docVIP

第九章 常微分方程初值问题数值解法 引 言 科学技术中常常需要求解常微分方程的定解问题．这类问题最简单的形式， 是本章将要着重考察的一阶方程的初值问题 我们只要，只要函数适当光滑—譬如关于满足利普希茨条件 . 理论上就可以保证初值问题，的解存在并且唯一. 虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程，实际问题中归结出来的微分方程主要靠数值解法. 所谓数值解法，就是寻求解在一系列离散节点 上的近似值相邻两个节点的间距 成为步长.今后如不特别说明，总是假定为定数，这时节点为 初值问题，的数值解法有个基本特点，它们都采取“步进式”，即求解过程顺着节点排列的次序一步一步地向前推进．描述这类算法，只要给出用已知信息计算的递推公式. 首先，要对方程离散化，建立求数值解的递推公式．一类是计算时只用到前一点的值，称为单步法．另一类是用到前面点的值称为步法．其次，要研究公式的局部截断误差和阶，数值解与精确解的误差估计及收敛性，还有递推公式的计算稳定性等问题． 简单的数值方法与基本概念 欧拉法与后退欧拉法 我们知道，在平面上，微分方程的解称为它的积分曲线．积分曲线上一点的切线斜率等于函数的值．如果按函数在平面上建立一个方向场，那么，积分曲线上每一点的切线方向均与方向场在该点方向相一致． 基于上述几何解释，我们从初始点出发，先依方向场在该点的方向推进到上一点，然后再从依方向场的方向推进到上一点，循此前进做出一条折线（图）． 一般地，设已做出该折线的顶点，过依方向场的方向再推进到，显然两个定点，的坐标有关系 即 这就是著名的欧拉公式．若初值已知，则依公式可逐步算出 例1 求解初值问题 ， 解 为便于进行比较，本章将用多种数值方法求解上述初值问题．这里先用欧拉方法，欧拉公式的具体形式为 取步长计算结果见表 表计算结果对比 初值问题有解，按这个解析式子算出的准确值同近似值一起列在表中，两者相比较可以看出欧拉方法的精度很差． 还可以通过几何直观来考察欧拉方法的精度．假设，即顶点落在积分曲线上，那么，按欧拉方法做出的切线便是过点的切线．从图形上看，这样定出顶点显著地偏离了原来的积分曲线，可见欧拉方法是相当粗糙的． 为了分析计算公式的精度，通常可采用泰勒展开将在处展开，则有 在的前提下，于是可得欧拉法的公式误差 ， 称为此方法的局部截断误差． 如果对方程从到积分，得 ． 右端积分用左矩形公式近似，再以代替，代替也得到，局部截断误差也是． 如果在中右端积分用右矩形公式近似，则得另一个式 称为后退得欧拉法． 后退的欧拉公式与欧拉公式有着本质的区别，后者是关于的一个直接的计算公式，这类公式称作是显式的；然而公式的右端含有未知的，它实际上是关于的一个函数方程，这类公式称作是隐式的． 显式与隐式两类方法各有特点．考虑到数值稳定性等其他因素，人们有时需要选用隐式方法，但使用显式算法远比隐式方便． 隐式方程通常用迭代法求解，而迭代过程的实质是逐步显示化． 设用欧拉公式 给出迭代初值，用它代入式的右端，使之转化为显示，直接计算得， 然后再用代入式，又有 ． 如此反复进行，得 ． 由于对满足利普希茨条件．由减得 由此可知，只要迭代法就收敛到解．关于后退欧拉法的公式误差，从积分公式看到它与欧拉法是相似的． 梯形方法 为得到比欧拉法精度高的计算公式，在等式右端积分中若用梯形求积公式近似，若用代替，代替，则得 称为梯形方法． 梯形方法是隐式单步法，可用迭代法求解．同后退的欧拉法一样，仍用欧拉方法提供迭代初值，则梯形法的迭代公式为 为了分析迭代过程的收敛性，将式与式相减，得 ， 于是有 式中为关于的利普希茨常数．如