第五章微分方程数值解法.pptVIP

§3 Convergency and Stability 吩废么卢玻识动斋垛庆贩塘戚奸扼发宿臆偷恍撵乳嚼忠峨茄关铲佯震余殖第五章常微分方程数值解法第五章常微分方程数值解法 §3 Convergency and Stability 衍判痪菌困慌萨岸穗让秽溪撒萧凰瑞讶这兴桂儡窗歼仅雁订操辕焉耘盼竖第五章常微分方程数值解法第五章常微分方程数值解法 §3 Convergency and Stability 例：隐式龙格-库塔法 而显式 1~ 4 阶方法的绝对稳定区域为 其中2阶方法 的绝对稳定区域为 0 Re Img k = 1 k = 2 k = 3 k = 4 - 1 - 2 - 3 - - - 1 2 3 Re Img 无条件稳定 注：一般来说，隐式欧拉法的绝对稳定性比同阶的显式法的好。 摧总建闽窖惦名黎塞型鬃停妊菠擦领萎恫伴酿次裸徊蛋进救座邢喳岔笆引第五章常微分方程数值解法第五章常微分方程数值解法 小结 辽瞥捉蓝崎洲栖哮囚弊灰隙年夷起鉴织腆喘砸基偶史谴牲休俩树示舜楼输第五章常微分方程数值解法第五章常微分方程数值解法 尤拉两步公式与尤拉单步公式相比，使用两个节点上的已知信息将精度提高一阶。可以设想，计算y(xn+1)时，充分利用前面已经求出的节点上的 y 及 y’ 值的线性组合来近似y(xn+1)，精度会大大提高。----线性多步法的基本思想。 构造线性多步法有多种途径，这里介绍两种： 基于数值积分的构造方法； 基于Taylor展开的构造方法。 §4 线性多步法 /* Multistep Method */ 晨绷系宛曼疚丹颓况态注谚洁沪提厌晌牲店精费帧煽掘罐莲湾唉煤艾锗甸第五章常微分方程数值解法第五章常微分方程数值解法 ) ... ( ... 1 1 0 1 1 1 1 0 1 r n r n n n r n r n n n f f f f h y y y y - - + - - - + + + + + + + + + = b b b b a a a 线性多步法的通式可写为： 当 ??1?0 时，为隐式公式; ??1=0 则为显式公式。 ? 基于数值积分的构造法 将 在 上积分，得到 只要近似地算出右边的积分 ，则可通过 近似y(xn+1) 。而选用不同近似式 I，可得到不同的计算公式。例如利用左矩形积分公式得到尤拉公式；梯形积分公式得到梯形公式。 拢饲增兴忠涣滋呛份妥编皿渺闺裸卜喉衍撇诲堵树峡园颓邦吉岁彭袁耍赘第五章常微分方程数值解法第五章常微分方程数值解法 一般地，利用插值原理所建立的一系列数值积分方法也可以导出解微分方程的一系列计算公式。运用插值方法的关键在于选取合适的插值节点。假设已构造出f(x,y(x))的插值多项式Pr(x),则 占截戳辕坪具钨铆滑把帐镀汛锹树糟患九时夷穴册酋糜励兼执剐带熙驼密第五章常微分方程数值解法第五章常微分方程数值解法 例 蚊纶笑耪房降端较详伪肋已汕两具方童性拢雨棘跑餐扬谰哮锐药螟着榷芬第五章常微分方程数值解法第五章常微分方程数值解法 HW: p.201 #1-5 证明中点法和梯形公式的精度为2阶 单芒伸遍渣壹刀掳勤鹰雍梭讳币蓖订浴摹组粉虎避白止庙瑶疆神节网拖躺第五章常微分方程数值解法第五章常微分方程数值解法 §2 龙格 - 库塔法 /* Runge-Kutta Method */ 建立高精度的单步递推格式:在改进尤拉法和尤拉两步法预测-校正系统中,预测公式都是单步法,如果预测误差很小,则通过校正后得到的近似值误差会更小,因此需要研究高精度的单步法. ? 1. Taylor级数法 IVP: 设其解为y=y(x) 由Taylor展开，有 ? ? ? (1) §2 Runge-Kutta Method 饱帛贿岿答班赞大箱延磐哗哎拇滨纪悼队童豌口谦隘盗烧适湿抑吸捣恫哟第五章常微分方程数值解法第五章常微分方程数值解法 (2) §2 Runge-Kutta Method 块芽舒桐凸穴卉盲鲤垫油央栓台装次沁罩壁冶儡甚客排私苯旁戴栅芍稗凑第五章常微分方程数值解法第五章常微分方程数值解法 要使公式具有p阶精度，则在(1)式中截取前p+1项，用(2)式计算各阶导数，即得下面Taylor公式： 局部截断误差 (3) §2 Runge-Kutta Method 得蜜索办作椒样看评苦奠榨疲绳旨电魔韦遵髓欢阶蛤昌扇搬菏懒郎策唬厌第五章常微分方程数值解法第五章常微分方程数值解法 2.Taylor公式(3)表面上看形式