Planet path using netwon's law with forward method

3 次查看（过去 30 天）
BabyMilo
BabyMilo 2012-2-18
编辑： Michael 2013-10-16
Given the initial position (ie x1,y1) and the velocity (ie vx,vy) of the planet.
How to plot a graph of the path of the planet around the center point (0,0) assuming the center is Earth, using time step of 10s for 500000s?
formula for forward method:
x(t+dt)=x(t)+v(t)*dt
v(t+dt)=v(t)+a(t)*dt
a=GM/r^2
G = 6.67e-11, M = 5.97e24
Many thanks.
  3 个评论
显示 1更早的评论
Walter Roberson
Walter Roberson 2012-2-18
Where is "r" defined?
If I recall correctly, this kind of forward prediction is unable to handle the orbit of Mercury, which has notable influence from drag of the rotational gravitational frame.
James Tursa
James Tursa 2012-2-19
@Walter: Well, OP has Euler's method listed which isn't going to do all that great with *any* orbit.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Gregory Vernon
Gregory Vernon 2012-2-19
You should look into velocity verlet methods.

类别

在 Help Center 和 File Exchange 中查找有关 Earth and Planetary Science 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by