Orbit Propagator

Plot a circular, elliptic, and hyperbolic trajectory.

Code:

  Code:

clc;
clear;
close all;
theta = [0:pi/90:2*pi];
subplot(3,1,1)
rc= 2./(1-0*cos(theta));
polar(theta,rc), title('Circular Trajectory - Eccentricity = 0 ')
subplot(3,1,2)
re= 2./(1-0.5*cos(theta));
polar(theta,re), title('Elliptical Trajectory - Eccentricity = 0.5 ')
subplot(3,1,3)
rh= 2./(1-1.5*cos(theta));
polar(theta,rh), title('Hyperbolic Trajectory - Eccentricity = 1.5 ')

Plot:

Orbit Propagator.jpg (48.11 KiB) Viewed 442 times