常微分方程数值解ppt课件.pptxVIP

常微分方程分为 （1）初值问题 （2）边值问题;称 在区域D上对 满足Lipschitz条件 是指:;二、初值问题解的存在性;要计算出解函数 y(x) 在一系列节点 a = x0 x1… xn= b 处的近似值 采用离散化方法。;三、初值问题的离散化方法;§2 欧拉方法 /* Euler’s Method */;;? 欧拉法的局部截断误差：;例1: 用欧拉公式求解初值问题;;可用来检验近似解的准确程度。;? 欧拉公式的改进：;;一般先用显式计算一个初值，再迭代求解。;? 梯形公式 /*trapezoid formula */;? 中点欧拉公式 /* midpoint formula */;;? 改进欧拉法 /* modified Euler’s method */;注:此法亦称为预测-校正法 /* predictor-corrector method */;在实际计算时，可将欧拉法与梯形法则相结合， 计算公式为 ;改进欧拉法的截断误差;例2: ;;通过计算结果得比较可以看出，改进的Euler方法;§3 龙格 - 库塔法 /* Runge-Kutta Method */;? 考察改进的欧拉法，可以将其改写为：;将改进欧拉法推广为：;Step 2: 将 K2 代入第1式，得到;要求 ，则必须有：; 为获得更高的精度，应该如何进一步推广？;考虑一阶常微分方程初值问题;; 若 ;;n级显式Runge-Kutta方法;二阶Runge-Kutta方法;由此得;另一方面;改进的Euler方法;中点方法;二阶Heun方法;二级Runge-Kutta方法不超过二阶;因此局部截断误差只能达到;三级Runge-Kutta方法;记;又由于;因此要使局部截断误差为 ，必须 ;Kutta方法;三阶Heun方法;三级Runge-Kutta方法不超过三阶; 四级R-K方法;局部截断误差为O(h5);附注：;注： ? 龙格-库塔法的主要运算在于计算 的值，即计算 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系：;§4 单步方法的收敛性与稳定性 /* Convergency and Stability */;解：该问题的精确解为 ;? 稳定性 /* Stability */;定义;我们称算法A 比算法B 稳定，就是指 A 的绝对稳定 区域比 B 的大。;例：考察显式欧拉法;例：考察隐式欧拉法;例：隐式龙格-库塔法;而显式 1~ 4 阶方法的绝对稳定区域为