《数值计算》课程教学课件（讲稿）第8章 常微分方程的数值解法（2/2）

m教材《数值分析》陈晓江主编武汉理工大学出版社Y参考书目《数值分析》李庆扬、王能超、易大义编清华大学出版社施普林格出版社《科学和工程计算基础》施妙根、顾丽珍编清华大学出版社上页下页返圆

上页 下页 返回 教材 《数值分析》 陈晓江主编 武汉理工大学出版社  参考书目 《数值分析》 李庆扬、王能超、易大义编 清华大学出版社 施普林格出版社 《科学和工程计算基础》 施妙根、顾丽珍编 清华大学出版社

第八章常微分方程的数值解法第一节问题的提出第二节欧拉法第三节龙格-库塔法第四节单步法的收敛性与稳定性第五节线性多步法第六节一阶方程组和高阶方程第七节边值问题的数值解法上页下页返圆

上页 下页 返回 第八章 常微分方程的数值解法 第一节 问题的提出 第二节 欧拉法 第三节 龙格-库塔法 第四节 单步法的收敛性与稳定性 第五节 线性多步法 第六节 一阶方程组和高阶方程 第七节 边值问题的数值解法

$4单步法的收敛性与稳定性一.收敛性定义若某算法对于任意固定的x=x,=xo+ih，当h>0(同时i→o)时有yi→(x;)，则称该算法是收敛的。y'=2y考察欧拉显式格式的收敛性。例：就初值问题y(0) = yo解： 该问题的精确解为 (x)= yoeax欧拉公式为yi+i=y;+hay;=(1+ah)ylim(1+ ah)/ah=e+h->0对任意固定的x=x;=ih，有J; = o(1 + 2h)*; /h上页= yo[(1+ 2h)/ahjai → oeai = y(x,)下页返圆

上页 下页 返回 一. 收敛性 定义 若某算法对于任意固定的 x = xi = x0 + i h，当 h0 ( 同时 i  ) 时有 yi  y( xi )，则称该算法是收敛的。 例：就初值问题 考察欧拉显式格式的收敛性。        0 y(0) y y y 解：该问题的精确解为 x y x y e  0 ( )  欧拉公式为 i i i i y y h y (1 h) y 1       0 y (1 h) y i i    对任意固定的x = xi = i h ，有 x h i i y y h / 0  (1   ) h e h h      1/ 0 lim(1 ) xi y e   0  §4 单步法的收敛性与稳定性 ( ) xi  y h xi y h    [(1   ) ] 1/ 0

二.稳定性y(x)=-30y(x) ,在区间[0,0.5]上的解。例：考察初值问题y(0)=1分别用欧拉显、隐式格式和改进的欧拉格式计算数值解J=e-30x节点X欧拉显式欧拉隐式精确解改进欧拉法0.01.00001.00001.00001.00000.12.5000x10-14.9787×10-2-2.00002.50000.26.2500x10-22.4788x10-34.00006.25000.31.5625x10-21.5626x1011.2341x10-4-8.00000.43.9063x10-33.9063x1016.1442x10-61.6000x1010.53.0590x10-79.7656x10-49.7656x101-3.2000×101上页下页返圆

上页 下页 返回 二. 稳定性 例：考察初值问题 在区间[0, 0.5]上的解。 分别用欧拉显、隐式格式和改进的欧拉格式计算数值解。        (0) 1 ( ) 30 ( ) y y x y x 0.0 0.1 0.2 0.3 0.4 0.5 节点 xi 欧拉显式 欧拉隐式 改进欧拉法 精确解 x y e 30  1.0000 2.0000 4.0000 8.0000 1.6000101 3.2000101 1.0000 2.5000101 6.2500102 1.5625102 3.9063103 9.7656104 1.0000 2.5000 6.2500 1.5626101 3.9063101 9.7656101 1.0000 4.9787102 2.4788103 1.2341104 6.1442106 3.0590107

定义若某算法在计算过的误差在以后的计入为复数。为保算中都逐步衰减，则证稳定性，通常取入的实部小于0一般分析时为简单起见，只广开王y'=ay当步长取为h时，将某算法应用于上式，并假设只在初值产生误差&。=Jy。-J。，则若此误差以后逐步衰减，就称该算法相对于h=ah绝对稳定，h的全体构成绝对稳定区域。我们称算法A比算法B稳定，就是指A的绝对稳定区域比B的大。上页下页返圆

上页 下页 返回 定义 若某算法在计算过程中任一步产生的误差在以后的计 算中都逐步衰减，则称该算法是绝对稳定的。 一般分析时为简单起见，只考虑试验方程 y   y λ为复数。为保 证稳定性，通常 取λ的实部小于0 当步长取为h 时，将某算法应用于上式，并假设只在初值 产生误差 ，则若此误差以后逐步衰减，就称该 算法相对于 绝对稳定， 的全体构成绝对稳定区域。 我们称算法A 比算法B 稳定，就是指A 的绝对稳定区域比B 的大。 0 0 0   y  y h   h h