Calling Planetary EphemerisSimulink Block From .m file. The resulting output is unclear.

3 views (last 30 days)
Tjd=juliandate(2000,0,0);
options = simset('SrcWorkspace','current');
[X, V] = sim('Planetary_Epmemeris',Tjd,options);
I am trying to get ephemeris data (position and velocity) for specified celestial bodies. The outputs that I am getting are not the three element cartisian vectors I assumed I would get for position and velocity. V: double and X: double. I have no idea what this output means, and the documentation gives no details on how to interpret it.

Answers (1)

Ryan G
Ryan G on 8 May 2014
Did you try starting from the example?
Implement the position and velocity for Saturn with respect to the Solar System barycenter for noon on January 1, 2000 using DE421 and AU units:
[position,velocity] = planetEphemeris([2451544.5 0.5],... 'SolarSystem','Saturn','421','AU') position = 6.3993 6.1720 2.2738 velocity = -0.0043 0.0035 0.0016
  1 Comment
jonathan
jonathan on 8 May 2014
Hey, look at that! I hadn't seen that there was a Matlab function for this, so I was useing SIM to call the block. Thanks! Any idea on the accuracy of these results?

Sign in to comment.

Categories

Find more on Gravitation, Cosmology & Astrophysics in Help Center and File Exchange

Products

Community Treasure Hunt

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

Start Hunting!