非线性方第9讲.docVIP

第五章 常微分方程初值问题数值解法 §1 引言 常微分方程定解问题分两类： 第一类——初值问题 例如 （1） （2） 第二类——边值问题 例如 边值问题的数值解法本课程不介绍。 微分方程数值解法的目标 求未知函数在给定的离散点列处的近似值。 重点内容 一阶常微分方程初值问题。 然后，推广到常微分方程组、高阶常微分方程初值问题。 初值问题（1）的数值解法的基本特点： 算法采用递推的方式。 顺着点列的次序逐次计算点处未知函数的近似值（称为数值解），计算过程中要利用已经求出的近似函数值。 单步法与多步法 如果微分方程数值解的递推公式中只用到前面一点的值，称为单步法；用到前面多个点处的函数值，称为多步法。 符号 ——未知函数在点处的精确值 ——未知函数在点处的近似值 §2 简单的数值方法与基本概念 一、基于数值微分的求解公式 1、欧拉法 设有等距节点。 用一阶向前差商作为在处导数的近似值，得 用近似值替代精确值，并令 得计算未知函数近似值的递推公式： （4） 公式（4）称为欧拉公式。 根据公式（4），由，可递推得到的近似值。 欧拉法的解释 从分析的角度——用在处的一阶泰勒展开式作为其近似值。 从几何的角度——用曲线在处切线上的值近似。 欧拉法的评价 求解精度较差，且递推过程中会产生误差积累。 2、后退欧拉法 若用一阶向后差商近似导数，则得 即有 。 用近似值替代精确值，并令 （5） 公式（5）称为后退欧拉公式，这是一个隐式公式，一般需要用解非线性方程的方法求出。 后退欧拉公式的评价 精度较差，计算量大。 例1 取，求解初值问题 解 本题中， ，。 该微分方程为伯努利方程，注意到初值条件，其理论解为 0.1 0.2 0.3 0.4 0.5 1.0954 1.1832 1.2649 1.3416 1.4142 0.6 0.7 0.8 0.9 1.0 1.4832 1.5492 1.6125 1.6733 1.7321 （1）使用Euler法求解，即 得 …… 数值结果列表如下： 0.1 0.2 0.3 0.4 0.5 1.1000 1.1918 1.2774 1.3582 1.4351 0.6 0.7 0.8 0.9 1.0 1.5090 1.5803 1.6498 1.7178 1.7848 （2）后退Euler法求解，即 得 0.1 0.2 0.3 0.4 0.5 1.0907 1.1741 1.2512 1.3231 1.3902 0.6 0.7 0.8 0.9 1.0 1.4529 1.5114 1.5658 1.6160 1.7848 3．预测——校正的欧拉公式 一般情况下，从向后Euler公式中得到的显式表达式是困难的。为了避免解非线性方程，采取向前Euler公式和向后Euler公式相结合的方法，得到预测——校正Euler法。 预测： （6.1） 校正： （6.2） 求解过程为： 例2 用预测—校正法求解例1中的问题。 解 根据公式（6.1）、（6.2）算出各个节点处的函数值，列表如下： 0.1 0.2 0.3 0.4 0.5 1.1000 1.1827 1.2599 1.3323 1.4003 1.0918 1.1763 1.2546 1.3278 1.3964 0.6 0.7 0.8 0.9 1.0 1.4645 1.5249 1.5818 1.6352 1.6850 1.4609 1.5216 1.5786 1.6321 1.6819 二、基于数值积分的求解公式 1．梯形公式 因为积分 所以 （7） 将上述公式中的积分用梯形公式计算，即 并用近似值替代精确值，代入公式（7），得 （8） 公式（8）称为梯形公式，它也是一个隐式公式。 2．改进的Euler公式 一般情况下，梯形公式（8）也需要解非线性方程才能得到。为了避免解非线性方程，采取预测——校正的欧拉公式类似的方法，称这种方法为改进的Euler法。 预测： （9.1） 校正： （9.2） 求解过程为： 改进的欧拉法既简化了计算，又可以使求解精度较高。 例3 利用改进Euler法解初值问题 取步长 解 本例中， 利用求解初值问题的改进Euler公式 由得 同理： 本例中微分方程是一阶齐次线性微分方程，其理论解为 精确值为 例4 利用改进Euler法求解例1。 解 数值结果如下表： 0.1 0.2 0.3 0.4 0.5 1.1000 1.1873 1.2687 1.3454 1.4