运动的身体中央,减少错误

3视图(30天)
大家好,
我已经写了代码采用数值描述的运动身体m在中央1米(e。g卫星围绕地球)。我已经完成了它的初始值设置粒子的轨迹是椭圆。
显然,需要能量和角动量守恒的所以我得到。然而,我希望尽可能减少误差,得到一个更准确的解决方案。
为能源,我设法得到一个很小的使用相对和绝对宽容。例如,RelTol AbsTol,一些参数的相对误差可以达到的最大值约为8% ~ 1 e -只有50期。然而,相对误差在角动量似乎没有受到任何影响我设置RelTol和AbsTol以及它达到最大值(这是周期性的,见下面的图)是稳定的维持在11.36%左右。
我已经试着玩的步长,它似乎没有影响。
你可以看到两张图片下面描述能量和角动量在50期。y轴描述了相对误差百分比(不)在每个参数和x轴描述时间的数量(50)。
有人有什么建议如何控制和减少角动量的相对误差?
提前谢谢。
6个评论
Meysam Mahooti
Meysam Mahooti 2017年4月15日
//www.tianjin-qmedu.com/matlabcentral/fileexchange/55533-satellite-orbits--models--methods-and-applications

登录置评。

接受的答案

大卫Goodmanson
大卫Goodmanson 2016年12月2日
嗨mms ^ 2,这一切都因为你一直使用的表达式角动量是不正确的。不包括乘以质量,角动量
L_nomass = v * x - vx * y
如果你跟踪表达式时期,其相对误差应该与能量。
3评论
FlyGuy
FlyGuy 2016年12月2日
别介意,我刚刚注意到詹姆斯的评论。
谢谢你们俩!

登录置评。

更多的答案(0)