龙格-库塔方法基本原理3.ppt

（1）用欧拉法计算 根据欧拉公式，将函数表达式及其初始值代入后，可得该系统仿真第一步的值： 数值积分公式应用 （2）用梯形法计算： 根据预报—校正公式，将函数表达式及其初始值代入后，可得仿真第一步的值。 用预报公式求起始值： 数值积分公式应用 再用校正公式得到系统仿真第一步的值： 数值积分公式应用 二阶龙格-库塔公式 （3）用二阶龙格—库塔法计算 根据公式先计算出两个系数，再计算仿真第一步的值： 数值积分公式应用 则系统仿真第一步的值为： 数值积分公式应用 四阶龙格—库塔（Runge—Kutta）法 （4）用四阶龙格—库塔公式计算 根据公式先计算出4个系数，再计算仿真第一步的值： 数值积分公式应用 数值积分公式应用 则系统仿真第一步的值为： 数值积分公式应用 从上述结果可以看出: 对于同一个系统进行仿真计算时，其值的精度是随着数值积分公式的变化而改变的，其中欧拉法计算精度最低，其次为梯形法和二阶龙格—库塔法，四阶龙格—库塔法计算精度最高。 数值积分公式应用 例2：用matlab演示 龙格-库塔法的误差估计 一个高精度的仿真方法必须将步长控制作为手段。实现步长控制涉及局部误差估计和步长控制策略两方面的问题。 龙格-库塔法的误差估计 RK方法的误差估计通常是设法找一个低一阶的R