体的太阳能系统模拟器——为什么有innacuracies在x飞机而不是别人?

1视图(30天)

显示旧的评论

撒母耳Leeney 2020年12月9日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/688939-n-body-solar-system-simulator-why-are-there-innacuracies-in-the-x-plane-but-not-the-others

figupload.png

我写了一个基本体模拟器,模拟太阳系中行星的运动。

系统读取行星的位置和速度在一个给定的时间,测量每颗行星的时间为一个完整的轨道,然后比较了模拟的准确性对一个高度精确的模型。

我已经分析了数据的x, y和z的飞机。有趣的是,仿真在y和z平面高度精确,但大约10 x不准确的飞机。请看附呈的图。

知道这是为什么吗?我想也许是因为行星旋转的轴的角动量一些影响他们的轨道在x飞机吗?或者是暗物质的作用?

这里是一些代码,“模拟”。

                         为nk = 2:长度(LY)
                        
                         t = 0;
                        
                         dt = LY (nk) / num_it;
                        
                         G = 6.67 e-11;
                        
                         t_finish = LY (nk);%的一个轨道周期选定的星球
                        
                         N_t =圆(t_finish / dt);
                        
                         N =长度(p);%在模拟找到的对象数
                        
                         %在质心坐标系
                        
                         v = v -意思是((m * (1 1 1))。* v) /意味着(m);
                        
                         %找到accelleration
                        
                         一个= a_v (p m G, sft);
                        
                         %找到能源开始
                        
                         (柯、PE) = E (v p G、m);
                        
                         %主要模拟
                        
                         为n = 1: N_t
                        
                         %的踢漂移踢二阶技术使用。这个teqnique
                        
                         %用于保存能量相当有效。来源:https://en.wikipedia.org/wiki/Leapfrog_integration
                        
                         %的第一个1/2踢
                        
                         v = v + a * dt / 2;
                        
                         %漂移
                        
                         p = p + v * dt;
                        
                         %的新accellerations
                        
                         一个= a_v (p m G, sft);
                        
                         %的第二个1/2踢
                        
                         v = v + a * dt / 2;
                        
                         %的新时间
                        
                         t = t + dt;
                        
                         结束
                        
                         :psave (nk) = p (nk:);
                        
                         :vsave (nk) = v (nk:);
                        
                         结束