在使用龙格-库塔(RK)方法对连续系统进行数字仿真时,为了保证数值计算的稳定性以及仿真结果具有足够的精度,通常采用变步长策略。为了有效地解决变步长仿真计算过程中,输出节点与计算节点不相吻合的问题,该文在前人工作的基础上,提出了一个具有大稳定域的四阶连续RK公式对。该公式对在不增加微分方程的右端函数值的计算次数的前提下,可以给出积分步距中任意一点上的数值解,因而具有更大的应用价值。仿真结果表明,该公式对是有效可行的。
