数值分析第九章常微分方程的数值解.ppt

数值分析;第九章 常微分方程的数值解;很多科学技术和工程问题常用常微分方程的形式建立数学模型.;一、 Euler方法;重复上述方法，当 x=x2 时;解：;　　　在假设yn = y(xn)，即第i步计算是精确的前提下，考虑的截断误差Rn = y(xn+1) ? yn+1 称为局部截断误差.;Lipschitiz条件：若存在正数L,使得对一切;特别当n=m-1时，有;后退Euler公式; 由于未知数yn+1同时出现在等式的两边，不能直接得到，故称为隐式欧拉公式，而前者称为显式欧拉公式.;2 梯形公式和改进Euler方法;;具有2阶精度.;中点欧拉公式;;改进欧拉法;??? 用梯形公式求解初值问题(步长h=0.2);例 用改进欧拉法求解初值问题;由 y(1) = y0 =1 计算得;二、Runge-Kutta方法;首先希望能确定系数 ?1、?2、p，使得到的算法格式有2阶精度，即在 的前提假设下，使得 ;(2) 将 K2 代入第1式，得到;(3) 将yn+1与y( xn+1 )在xn点的泰勒展开作比较;所有满足上式的统称为2阶Runge-Kutta格式.;2 四阶Runge-Kutta方法;由于方程的个数少于未知量的个数，所以方程有无穷多个解，可以根据情况得到几种常用的解，即得到相应的四阶公式.;Gill公式;;例 用标准四阶Runge-Kutta法求初值问题;那么;解：;　　　 若某算法对于任意固定的x=xi =x0 + i h, 当 h?0 (同时i ? ?)时有yi ?y(xi )，则称该算法是收敛的. ;Euler方法的增量函数;例：考察欧拉显式格式的收敛性:; 设y(x)是微分方程初值问题的准确解，;单步法的收敛定理;在f(x,y)对y满足Lipschitz条件下, Euler法，改进Euler法和Runge-Kutta法的增量函数;考虑 在xn处的Taylor展式;例 设;所以;例：考察初值问题 在区间; 若一种数值解法仅在节点值yn上有大小为δ的扰动，于以后各节点值ym(mn)上，仅由δ所引起的扰动都不超过δ时，称该方法是稳定的.;反复应用可得;显式欧拉法:;隐式欧拉法(后退欧拉法)：;梯形公式：;龙格-库塔法：;例 设;解：;四、 线性多步法;特别的当 时, 上式为显式, 否则是隐式.;例 设 yn+1 = yn-1+2hf(xn, yn) 为求解常微分初值问题的线性二步法，试求该二步公式的局部截断误差主项，和精度.;代入可得;考虑 在xn处的Taylor展式;解得;考虑 在xn处的Taylor展式;所以;2 Adams外推公式;这里;而;i;3 Adams内插公式;i;4 预报-校正公式;预测-校正技术即保证了计算精度，又使隐式计算显式化，克服了隐式公式要反复迭代的困难.;例 用Adams外推公式和预报-校正公式求初值问题;k