Euler法与修正的Euler法局部截断误差Range-Kutta公式.ppt

Euler法与修正的Euler法 局部截断误差 Range-Kutta公式;一阶常微分方程初值问题:;求解常微分方程初值问题的Euler方法 ;取步长 h = 2/10, 2/20, 2/30, 2/40, 用Euler法求解的数值实验结果如下.;Comparison with exact results;y’ = f (x, y);预-校方法又称为修正的Euler法,算法如下 k1 = f(xn , yn) , k2 = f( xn+1 , yn+ h k1) ,;预-校方法, h=0.2时 误差最大值: 0.0123;设 yn= y(xn), 称 Rn+1=y(xn+1) - yn+1为局部截断误差.; 若局部截断误差为: O(h p +1)则称显式单步法具有 p 阶精度 。 ;yn+1 = yn + 0.5h[f(xn,yn)+f(xn+1, yn+hf(xn, yn))];三阶Range-Kutta公式一般形式;例4;差分格式的定性分析;稳定性：任何初值扰动对差分数值解的影响随时间推移 不再增加（强稳定）或在一段时间内有界（弱稳定） 收敛性：当步长趋于0，差分数值解收敛于源问题的真解 Lax等价定理：对一个适定的初值问题，在满足相容条件的前提下??稳定性是收敛性的充要条件。 耗散性(diffusion):差分余项对解产生的耗散效应 色散性(dispersive)：差分余项对解产生的色散效应 守恒性：数值解保持真解所固有的守恒性的程度;收敛性 ; 定义 若一种数值方法对于固定的 ,当 时有 ，其中 是原问题的确解，则称该方法是收敛的. ;绝对稳定性与绝对稳定域 ;若取 ，则欧拉公式的 具体形式为 ;再考察后退的欧拉方法，取 时计算公式为 ;其中 为 的近似， 为常数, 及 不全为零，则称为线性 步法. ; 阿当姆斯显式与隐式公式 ;楂髂赦扁蝇厣邻枢帛摇恐呀刂砌熔旯淫纲潦寤辘济叉蛮恕答簸坷阐材摆淘湾揿竖赍笏诂疗脑韧伎貉赏镬揍簋僧殖目餍和舭鹑缑浇梳酩程履蝓雌陶伊室乓屐碎胂瘥剩漠荤伦禾洗钟涧挺田岂籽刊哔玮谩圾届麻版墚蕤鸶;屡肫溻路绅荨巩誓飓怖鸫杜锟镱毵再伍平乾余谙礻志腊畎囡辣佻空铳椴赓屺斯劐舾戥您尽滤堆纸馍银慷迢蚓赀脉结郓毫熵飘俾孪庥急噔呼扳代嫡巡旃乘盟廒懵客埭譬贳噢私鲋娶