数值分析方法第七章.pptVIP

例 用欧拉法求初值问题 当h = 0.02时在区间[0, 0.1]上的数值解。 具体计算结果 表中y(xn)，是初值问题的真解 在xn上的值。 为近似值yn的误差。 从表中可以看出，随着n的增大，误差也在增大，所以说，欧拉法计算简便，对一些问题有较大的使用价值，但是，它的误差较大，所得的数值解精确度不高。 * ? 考虑一阶常微分方程的初值问题 /* Initial-Value Problem */: 其中f (x, y)为x, y的已知函数，y0为给定的初始值 第六章 常微分方程数值解 /* Numerical Methods for Ordinary Differential Equations */ 例如： 其解析解为： 只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解，而大量的微分方程问题很难得到其解析解。 这是微积分的发明者之一Leibniz在1686年曾经让当时 数学界人士求解的一阶微分方程式，吸引了许多数学家 的注意，大约经过150年的探索 到1838年, 刘维尔（Liouville)在理论上证明了这个微分方程不能用初等积分法求解，得借助于数值方法 只能依赖于数值方法去获得微分方程的数值解。 数值方法的基本思想是：在解的存在区间上取n + 1个节点 这里 hi可以不相等，但一般取成相等的，这时 在这些节点上采用离散化方法，（通常用数值积分、微分、泰勒展开等）将上述初值问题化成关于离散变量的相应问题。把这个相应问题的解yn作为y(xn)的近似值。这样求得的yn就是上述初值问题在节点xn上的数值解。一般说来，不同的离散化导致不同的方法。 ，i = 0,1, …, n-1称为由xi到xi+1的步长。 步进式： 根据已知的或已求出的节点上的函数值计算当前节点上的函数值，一步一步向前推进。 因此只需建立由已知的或已求出的节点上的函数值求当前节点函数值的递推公式即可。 §1 欧拉方法 /* Euler’s Method */ ? 欧拉公式： x0 x1 向前差商近似导数 记为 亦称为欧拉折线法 /* Euler’s polygonal arc method*/ 也称欧拉折线法. 用这条折线近似地代替曲线 由Euler法所得的折线明显偏离了积分曲线，可见此方法非常粗糙。 解：把 代入欧拉法计算公式，得 0.0021 0.923 0.9192 0.10 5 0.0018 0.9354 0.9336 0.08 4 0.0014 0.9503 0.9489 0.06 3 0.0005 0.9660 0.9650 0.04 2 0.0005 0.9825 0.9820 0.02 1 0 1.0000 1.0000 0 0 ?n = y(xn) - yn y(xn) yn xn n 定义 　　　在假设 yi = y(xi)，即第 i 步计算是精确的前提下，考虑的截断误差 Ri = y(xi+1) ? yi+1 称为局部截断误差 /* local truncation error */。 定义 　　　若某算法的局部截断误差为O(hp+1)，则称该算法有p 阶精度。 ? 欧拉法的局部截断误差： 欧拉法具有 1 阶精度。 Ri 的主项 /* leading term */ 衡量求解公式好坏的一个主要标准是求解公式的精度, 因此引入局部截断误差和阶数的概念。 Euler’s Method ? 欧拉公式的改进： ? 隐式欧拉法 /* implicit Euler method */ 向后差商近似导数 x0 x1 )) ( , ( ) ( 1 1 0 1 x y x f h y x y + ? ) 1 , ... , 0 ( ) , ( 1 1 1 - = + = + + + n i y x f h y y i i i i § Euler’s Method 由于未知数 yi+1 同时出现在等式的两边，不能直接得到，故称为隐式 /* implicit */ 欧拉公式， 而前者称为显式 /* explicit */ 欧拉公式。 一般先用显式计算一个初值，再迭代求解。 ? 隐式欧拉法的局部截断误差： 即隐式欧拉公式具有 1 阶精度。 Hey! Isn’t the leading term of the local truncation error of Euler’s method ? Seems that we can make a good use of it … § Euler’s Method ? 梯形公式 /* trapezoid for