Elliptical Orbit Plotting_28 Apr 2020

I am using MATLAB R2015a.

I have the following two equations with intial values:

wn = 4000*(pi./30); % rad/sec

T = (2*pi./wn);

t = (0:T./200:T);

Xc = 1.0; Xs = 0.50; Ys = 0.20; Yc = 1.40;

X = (Xc*cos(wn*t)+Xs*sin(wn*t)) = Mod_X.*cos(wn.*t-PHASE_X);

Y = (Yc*cos(wn*t)+Ys*sin(wn*t)) = Mod_Y.*cos(wn.*t-PHASE_Y);

For the ELLIPTICAL ORBIT:

Mod_X = sqrt(Xc.^2+Xs.^2); Mod_Y = sqrt(Yc.^2+Ys.^2);

PHASE_X = atan(Xs./Xc); PHASE_X_Degree = PHASE_X.*(180./pi);

PHASE_Y = (pi./2)- atan(Ys./Yc); PHASE_Y_Degree = PHASE_Y.*(180./pi);

I want to calculate the following: (Final Values are in BRACKETS)-

(i.) Semi-Major Axis (Value = 1.61245); (ii.) Semi-Minor Axis (Value = 0.806226); (iii.) Attitude Angle (Value = 56.3099)

And for the Two-Circular Orbits Data:

(i.) Forward Amplitude (Value = 1.20934); Phase Angle (Value = 362.875); (ii.) Backward Amplitude (Value = 0.403113); Phase Angle (Value = 119.745);

I am unable to get correct inclination of ELLIPTICAL ORBIT, two CIRCLES of FORWARD and BACKWARD Precession. For plotting i am using following command:

figure(1)

h = plot (X,Y,'-c');

set(h(1),'linewidth',2.5); xlabel('X, [m]'); ylabel('Y1, [m]'); grid on;

Help me in trhis regard.

Sandeep Soni