第五章常微分方程数值解法第一节欧拉法.ppt

第五章 常微分方程数值解 /* Numerical Methods for Ordinary Differential Equations */;解析解法：（常微分方程理论） 只能求解极少一类常微分方程；实际中给定的问题不一定是解析表达式，而是函数表，无法用解析解法。;--------Euler’s Method;几何意义;定义;截断误差: 实际上，y(xn) ? yn， yn 也有误差，它对yn+1的误差也有影响，见下图。但这里不考虑此误差的影响，仅考虑方法或公式本身带来的误差，因此称为方法误差或截断误差。 局部截断误差的分析：由于假设yn = y(xn) ，即yn准确，因此分析局部截断误差时将y(xn+1) 和 yn+1都用点xn上的信息来表示，工具：Taylor展开。;定义;§1 Euler’s Method;由于未知数 yn+1 同时出现在等式的两边，故称为隐式 /* implicit */ 欧拉公式，而前者称为显式 /* explicit */ 欧拉公式。隐式公式不能直接求解，一般需要用Euler显式公式得到初值，然后用Euler隐式公式迭代求解。因此隐式公式较显式公式计算复杂，但稳定性好; 见上图， 显然，这种近似也有一定误差， 如何估计这种误差y(xn+1) ? yn+1 ？ 方法同上，基于Taylor展开估计局部截断误差。 但是注意，隐式公式中右边含有f(xn+1 , yn +1 ) , 由于yn +1不准确，所以不能直接用 y (xn+1)代替f(xn+1 , yn +1 ) ;? 隐式欧拉法的局部截断误差：;§1 Euler’s Method;? 隐式欧拉法的局部截断误差：;比较尤拉显式公式和隐式公式及其局部截断误差; 若将这两种方法进行算术平均，即可消除误差 的主要部分/*leading term*/ 而获得更高的精度,称为梯形法;梯形法的迭代计算和收敛性;梯形法的简化计算 迭代计算量大，且难以预测迭代次数。为了控制计算量，通常只迭代一次就转入下一点的计算。用显式公式作预测，梯形公式作校正，得到如下预测校正系统，也称为改进尤拉法：;注：此法亦称为预测-校正法 /* predictor-corrector method */。可以证明该算法具有 2 阶精度，同时可以看到它是个单步递推格式，比隐式公式的迭代求解过程简单。后面将看到，它的稳定性高于显式欧拉法。;几何解释;令x=x1,得;? 中点欧拉公式 /* midpoint formula */;几何解释;; 算例 分别用显式Euler方法，梯形方法和预估－校正Euler方法解初值问题; 算例 分别用显式Euler方法，梯形方法和预估－校正Euler方法解初值问题;;例：确定下列公式：