[工学]第六章常微分方程初值问题的数值解法n.pptVIP

常微分方程初值问题的数值解法 举例： 本章要求 要求熟练掌握的内容： 欧拉格式、隐式欧拉格式、梯形格式和改进欧拉格式的基本公式，计算步骤，算法； 运用经典龙格－库塔法求常微分方程数值解。 要求掌握的内容： 龙格－库塔法的基本思想，计算格式的导出； 单步法局部截断误差及阶的定义和计算。 注： ? 龙格-库塔法的主要运算在于计算 Ki 的值，即计算 f 的值。Butcher 于1965年给出了计算量与可达到的最高精度阶数的关系： 7 5 3 可达到的最高精度 6 4 2 每步须算Ki 的个数 ?高于四阶时每步计算量增加较多，但精度提高不快，因此使用的比较少。 ?由于龙格-库塔法的导出基于泰勒展开，故精度主要受解函数的光滑性影响。对于光滑性不太好的解，最好采用低阶算法而将步长h 取小。 Lab 15. Runge-Kutta (Order Four) Use Runge-Kutta method of order four to approximate the solution of the initial-value problem , , and (1) You are supposed to write a function void RK