第五章常微分方程数值解法第一节欧拉法演示文稿.pptVIP

第五章常微分方程数值解法第一节欧拉法演示文稿 当前第1页\共有31页\编于星期三\9点 优选第五章常微分方程数值解法第一节欧拉法 当前第2页\共有31页\编于星期三\9点 解析解法：（常微分方程理论） 只能求解极少一类常微分方程；实际中给定的问题不一定是解析表达式，而是函数表，无法用解析解法。 如何求解 计算解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 节点间距 为步长，通常采用等距节点，即取 hi = h (常数)。 数值解法: 求解所有的常微分方程 步进式：根据已知的或已求出的节点上的函数值计算当前节点上的函数值，一步一步向前推进。因此只需建立由已知的或已求出的节点上的函数值求当前节点函数值的递推公式即可。 当前第3页\共有31页\编于星期三\9点 当前第4页\共有31页\编于星期三\9点 当前第5页\共有31页\编于星期三\9点 --------Euler’s Method §1 欧拉方法 /* Euler’s Method */ §1 Euler’s Method Taylor展开法 当前第6页\共有31页\编于星期三\9点 几何意义 亦称为欧拉折线法 /* Euler’s polygonal arc method*/ 几何直观是帮助我们寻找解决一个问题的思路的好办法哦 当前第7页\共有31页\编于星期三\9点 定义 　　　在假设 yn = y(xn)，即第 n 步计算是精确的前提下，考虑公式或方法本身带来的误差: Rn = y(xn+1) ? yn+1 , 称为局部截断误差 /* local truncation error */。 说明 显然，这种近似有一定误差， 而且步长越大，误差越大， 如何估计这种误差y(xn+1) ? yn+1 ？ §1 Euler’s Method 当前第8页\共有31页\编于星期三\9点 截断误差: 实际上，y(xn) ? yn， yn 也有误差，它对yn+1的误差也有影响，见下图。但这里不考虑此误差的影响，仅考虑方法或公式本身带来的误差，因此称为方法误差或截断误差。 局部截断误差的分析：由于假设yn = y(xn) ，即yn准确，因此分析局部截断误差时将y(xn+1) 和 yn+1都用点xn上的信息来表示，工具：Taylor展开。 ? 欧拉法的局部截断误差： Rn+1 的主项 /* leading term */ §1 Euler’s Method 当前第9页\共有31页\编于星期三\9点 定义 　　　若某算法的局部截断误差为O(hp+1)，则称该算法有p 阶精度。 ? 欧拉法具有 1 阶精度。 在xn点用一阶向前差商近似一阶导数 在第2章讨论牛顿插值公式时 介绍了差商的概念和性质， 各阶差商可以近似各阶导数，具有不同的精度， 且可以用函数值来表示。 上一章中数值微分的方法之一 就是用差商近似导数 Euler’s method §1 Euler’s Method 当前第10页\共有31页\编于星期三\9点 §1 Euler’s Method ? 欧拉公式的改进： 隐式欧拉法或后退欧拉法 /* implicit Euler method or backward Euler method*/ xn+1点向后差商近似导数 隐式或后退欧拉公式 当前第11页\共有31页\编于星期三\9点 由于未知数 yn+1 同时出现在等式的两边，故称为隐式 /* implicit */ 欧拉公式，而前者称为显式 /* explicit */ 欧拉公式。隐式公式不能直接求解，一般需要用Euler显式公式得到初值，然后用Euler隐式公式迭代求解。因此隐式公式较显式公式计算复杂，但稳定性好 收敛性 §1 Euler’s Method 当前第12页\共有31页\编于星期三\9点 见上图， 显然，这种近似也有一定误差， 如何估计这种误差y(xn+1) ? yn+1 ？ 方法同上，基于Taylor展开估计局部截断误差。 但是注意，隐式公式中右边含有f(xn+1 , yn +1 ) , 由于yn +1不准确，所以不能直接用 y (xn+1)代替f(xn+1 , yn +1 ) 设已知曲线上一点 Pn (xn , yn ),过该点作弦线，斜率为(xn+1 , yn +1 ) 点的方向场f(x,y)方向,若步长h充分小，可用弦线和垂线x=xn+1的交点近似曲线与垂线的交点。 几何意义 xn xn+1 Pn Pn+1 x y y(x) §1 Euler’s Method 当前第13页\共有31页\编于星期三\9点