第5章-解常微分方程.pptxVIP

第五章 解常微分方程计算物理基础本章内容§5.1 龙格-库塔法§5.5 刚性问题§5.9 物理应用§5.2 变步长的龙格-库塔法§5.6 事件§5.10 边值问题和本征值§5.11 指令bvp4c解边值问题与本征值问题（阅读）§5.3 常微分方程组的初值问题§5.7 误差（阅读）§5.4 用MATLAB解常微分方程§5.8 性能（阅读）解初值问题的指令ode23,ode45解刚性问题的指令ode23s解事件问题的指令events解边值问题的指令bvp4c初值问题的龙格-库塔法本征值问题的打靶法掌握解常微分方程的MATLAB指令21学习解常微分方程的数值方法学习指导§5.1 龙格-库塔法§5.1.1 基本思想1. 解常微分方程的数值算法?一阶微分方程初值问题： 取步长： 区间分成个子区间：§5.1.1 基本思想?求由中值定理：得平均斜率为若确定，就能从求所以各种方法的实质：确定不同精度的?§5.1.1 基本思想2. 欧拉公式点斜率值点斜率值?? 取，欧拉公式 若取欧拉公式改进为 ?由于未知，所以采用预报法来计算：先用欧拉法求得的预报值 =再用预报值计算处的斜率值§5.1.1 基本思想?得改进的欧拉公式?龙格-库塔法基本思想:多点的斜率值加权平均当作平均斜率，以提高精度。?§5.1.1 基本思想一般的显式龙格-库塔方法：???是常数确定原则：使局部截断误差关