八章 常微分方程组的数值解.pptVIP

在工程和科学技术的实际问题中，常需求解微分方程，但常微分方程中往往只有少数较简单和典型的常微分方程（例如线性常系数常微分方程等）可求出其解析解,对于变系数常微分方程的解析求解就比较困难，而一般的非线性常微分方程的求解困难就更不用说了。大多数情况下，常微分方程只能用近似方法求解。这种近似解法可分为两大类：一类是近似解析法，如级数解法、逐次逼近法等;另一类是数值解法，它给出方程在一些离散点上的近 似值。 其中 x 是质量，m是离开平衡位o的距离，t为时间，c为弹簧系数。 在具体求解微分方程时，需具备某种定解条件，微分方程和定解条件合在一起组成定解问题。定解条件有两种：一种是给出积分曲线在初始点的状态，称为初始条件，相应的定解问题称为初值问题。另一类是给出积分曲线首尾两端的状态，称为边界条件，相应的定解问题称为边值问题。 我们现在讨论常微分方程的数值解法。先从最简单的一阶常微分方程的初值问题出发开始讨论。 由常微分方程理论可知：只要上式中的函数f(x,y)在区域 G={a≤x≤b,－∞＜y＜∞}内连续，且关于 y 满足Lipschitz条件，即存在与 x, y 无关的常数L，使 下面分析均假定满足上述条件。 这就是Euler公式（格式）。 利用它可由初值 出发逐步算出 。 这类形式的方法也称为差分方法。 定义：如果局部截断误差为 ，则这种数值算法的精度为p阶，故Euler格式的精度为一阶。 从几何意义上来看，如图， 当假定 为准确值，即在 的前提下来估计误差 ，这种截断误差称为局部截断误差。 由（2）、（3）知Euler公式在 处的局部截断误差为： 即为Euler格式（3）。 显式与隐式两类方法各有特点，使用显式算法远比隐式算法方便，但考虑数值稳定性等因素，人们常选用隐式算法。 隐试算法（6）常用迭代法来实现，而迭代过程实质上是逐步 显式化。 对右端利用左矩形公式可得 即 Euler格式 梯形公式也是隐式的，可用迭代法求解，与后退Euler方法一样，仍用Euler方法提供迭代初值，其迭代格式为： 为分析迭代过程的收敛性，将（12）与（8）相减得： （12） k=0，1，2，….. L为f(x,y)关于y的Lipschitz常数.如果选取h充分小使得 则 . 时有 这表明迭代过程（12）是收敛于（8）的解的。 由表可见，与精确解 相比，改进的Euler公式的精度较Euler 公式有明显的提高。 下面再看两步Euler公式（9），除了给出初值 外，还需要借助 其它单步法（如Euler公式，后退Euler公式及梯形公式等）再提供一个 Euler公式 改进的Euler公式 精确解 0 1 1 1 0.1 0.2 0.3 1 0.9000000 0.8100000 0.7290000 0.3486784 0.9050000 0.8190250 0.7412176 0.3685410 0.9048374 0.8187308 0.7408182 0.3678794 启动值 然后才能启动计算公式依次计算 用两步Euler公式与梯形公式相匹配，又可得到下面预测-校正 系统： （18） 两步法优美是由于它调用了两个节点上的信息，从而能以较少的 计算量获得较高的精度。 预测： 与改进的 Euler公式（13）（14）相比较易见（17）（18）的一个 突出特点是它的预测公式与校正公式具有同阶精度。据此可以比较 方便的估计截断误差，并基于这种估计，可以提供一种提高精度的 简易方法。 若预测公式（17）中的 和 都是准确的。即 则由两步Euler公式的截断误差公式（10）知： 而具体则有: