第八章+常微分方程的数值解法.pptVIP

k1 k2 k* 公式中含有三个参数λ1，λ2和p，如果我们适当选取参数的值，可以使公式的局部截断误差为O(h3)。 对k1和k2作泰勒展开 代入（8-7）得 （*） 又 y (x)在xi处的二阶泰勒展开式为 当x = xi+1时， ，有 （**） （**） 比较（*）与（**）的系数即可发现， 要使公式（7-7）的局部截断误差满足 ，即要求公式具有二阶 精度只要下列条件成立即可。 （8-8） 满足条件（8-8）的一簇公式统称为 二阶龙格—库塔公式。 特别的，当 塔公式就成为改进欧拉公式。 时，龙格-库 　　改进欧拉公式就是以y(x)在xi和xi+1处 的斜率k1和k2的算术平均值作为y(x)在 [xi，xi+1]上的平均斜率k*来进行计算的。 若取 时，龙格-库塔 公式就称为变形的欧拉公式，其形式为 （8-9） 此处的 就是欧拉方法预报出的中点 处的近似解；而 等于中点的斜率值 则近似 ； 所以公式可以看作用中点斜率近似 代替平均斜率k*，因此，公式（8-9）也 称作中点公式。 第八章 常微分方程的数值解法 工程技术和自然科学中的许多实际问题，在数学上往往可以归结为求解微分线性方程的形式。高数中我们学过的解析方法主要是用在一些简单和特殊的微分方程求解中，而对于大量的一般形式的微分方程，不能用解析方法求出其精确解，而只能用数值方法求近似解。 在具体求解微分方程时，要有某种定解条件，微分方程与定解条件在一起称为定解问题。定解条件有两种：一种是给出积分曲线在初始点的状态（初始条件），相应问题称作初值问题；另一种是给出积分曲线首尾两端的状态（边界条件），相应的问题称做边值问题。 本章主要研究一阶常微分方程初值问题的数值解法。 一阶常微分方程的初值问题是指形如 的问题。当函数 f (x，y)满足 时初值问题（8-1）的解y = y(x)存在且唯 一，此条件称为李普希兹 (Lipschitz)条 件。L为正常数，称为Lipschitz常数。 （8-1） 我们讨论的数值解法基本思想如下： 常微分方程的数值解法期待建立求解y(xi)的近似值yi的递推格式，通过选取适当的初始值，求得y(x)在各离散点上的近似值——差分方程。 在区间[a，b]上插入一系列离散点 第二节 欧拉方法 欧拉法是解初值问题最常用方法， 我们用它来说明求解微分方程的基本 的技巧和概念。 欧拉法的计算格式是 (8-2) 这个递推公式可以在已知 yi 时推出 yi+1，下面给出公式的几种解释： （1）利用y (x)在x = xi处的一阶泰勒公 式，得 舍去二阶导数以后的余项，并以 yi和yi+1分别表示y(xi)和y(xi+1)的近似 值，就得到 （2）对于充分小的h，将微分方程中的 用差商 代替， 解出yi+1即可得到（8-2）式。h=xi-xi-1称 为步长；一般取等步长 h = (b-a)/n。 导数 并以yi和yi+1分别表示y(xi)和y(xi+1) 的近似值，就得到 如果采用左矩形公式近似计算上式中的积 分，并以yi和yi+1分别表示y(xi)和y(xi+1)的近似 值，也可得到 称为向前欧拉公式，也称为显式欧拉公式。 此公式对给定的y0，可以算出y1，然后一步一 步算出yi+1。 （8-2） （3）从xi到xi+1对式 积分，得 （8-3） 如果对（8-2）式中的积分采用右矩形公 式，将得到计算格式 称为向后欧拉公式，也称为隐式欧拉公 式。由于未知量yi+1含在右端，因而不能直接 计算，一般使用迭代法进行求解。 (8-4) 欧拉法的实质是将连续问题 化为在给定节点上的差分方程（离散问题）。 梯形公式 上述的两种公式是用矩形求积公式来计算 （8-3）式右端的积分而得到的，由于矩形公式 精度较低，所以导出的欧拉公式精度也很低。 如果用数值积分的梯形公式计算（8-3）式右端 的积分，则有 （8-5） 这就是梯形公式。梯形公式由于其右端含 有未知的yi+1，所以也是隐式方程。 以 f (x0，y0)为斜率过点(x0，y0)做一条 直线，它与直线x=x1的交点就是y1。依此 类推，yi+1是以f (xi，yi)为斜率过点（xi， yi）的直线与x=xi+1的交点，所以欧拉法也 称为欧拉折线法