Sun-Synchronous Orbit Demo

Calculate the synchronous inclination given the other elements See Vallado, 2nd Edition, Figure 11-11

------------------------------------------------------------------------------

%	 Copyright (c) 2012 Princeton Satellite Systems, Inc. All rights reserved.
% ------------------------------------------------------------------------------

e = [0 0.1 0.2 0.3 0.4 0.5];
a = 6500:100:12000;

i = zeros(length(e),length(a));
for k = 1:length(e)
  i(k,:) = SunSyncInclination( a, e(k) );
end

i = i*180/pi;

Plot2D(i(1,:),a,'Inclination (deg)','Semimajor Axis (km)','Sun-Synchronous Elements')

c = 'bgrcmy';
hold on
for k = 2:length(e)
  plot(i(k,:),a,c(k))
end

legend(num2str(e'))

axis([94 121 6500 13000])

%--------------------------------------
% PSS internal file version information
%--------------------------------------