Demonstrate the use of FFSim to analyze disturbance effects.

------------------------------------------------------------------------
See also FFSim, FFSimPlotter, Goals2DeltaElem
------------------------------------------------------------------------

Contents

%--------------------------------------------------------------------------
%   Copyright (c) 2009 Princeton Satellite Systems, Inc.
%   All rights reserved.
%--------------------------------------------------------------------------
%   Since version 8.
%--------------------------------------------------------------------------

if( ~exist('FormationFlying/Transformation/Goals2DeltaElem','file') )
  warning('This demo requires the Formation Flying Module.')
  return
end

Initialize simulation data structure

% Reference orbital elements [a, theta, i, q1, q2, W]
%  a        semi major axis
%  theta    latitude (perigee + true anomaly)
%  i        inclination
%  q1       e*cos(w), e is eccentricity, w is perigee
%  q2       e*sin(w)
%  W        right ascension
clear s;
s.el0       = [6928.14 0 0.617846555205993 0 0 0];

% Choose desired relative motion to be 250 m relative ellipse
s.goals.y0 = 0;      % center of ellipse in along-track direction
s.goals.aE = .25;    % size of ellipse
s.goals.beta = 0;    % phase on ellipse at equator crossing
s.goals.zInc = .1;   % cross-track oscillation (inclination difference)
s.goals.zLan = -.1;  % cross-track oscillation (right ascension difference)

% initialize relative state to achieve desired relative motion
s.dEl0      = Goals2DeltaElem( s.el0, s.goals );

% other simulation options
s.nOrbits   = 5;
s.nSPO      = 300;
s.distOpt   = [1 1 1];
s.area      = [1 .5];
s.mass      = [50 10];

Simulate

d = FFSim(s);

5.0 orbits remaining ... should be done in 0 hours, 2 minutes, 5 seconds
4.0 orbits remaining ... should be done in 0 hours, 0 minutes, 1 seconds
3.0 orbits remaining ... should be done in 0 hours, 0 minutes, 1 seconds
2.0 orbits remaining ... should be done in 0 hours, 0 minutes, 0 seconds
1.0 orbits remaining ... should be done in 0 hours, 0 minutes, 0 seconds

View results...

FFSimPlotter( d );


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

Published with MATLAB® R2019a