所附的脚本使用MATLAB支持的子例程ode45、ode23或ode23s,并使用简单的不支持的子例程ode4(可在附件和文件交换的其他地方获得)来解决两体轨道矢量常微分方程。万博1manbetx所有这些子程序都使用龙格-库塔算法求解ODE。脚本将结果与精确解析解进行比较,精确解析解是一个适当的椭圆:因此我们可以比较子程序的准确性。我们发现所支持的子程序显示的误差范万博1manbetx围在0.1到0.5%之间;而ode4的峰值误差小于10^-7%,即小10^6倍。这些结果表明ode4可能是第一个尝试求解ODE的龙格-库塔子程序。附加的信息包含这些结果的图以及复制结果所需的所有脚本。
引用作为
约翰Keevil(2021)。ODE4给出了比ODE45、ODE23、ODE23更准确的结果(//www.tianjin-qmedu.com/matlabcentral/fileexchange/59044-ode4-gives-more-accurate-results-than-ode45-ode23-ode23s), MATLAB中央文件交换。检索.