Problem 12.33 Numerically integrate the Gauss planetary equations for a given perturbation p and set
of initial conditions.
Let us choose an earth orbit in which a constant tangential thrust is the perturbing force. As shown in
Equation 6.26, the perturbing acceleration in that case is
where T is the rocket thrust, m is the spacecraft mass, v is the velocity and v is the speed. For this problem
let the state vector in the geocentric equatorial frame at time t = 0 be
r
07000ˆ
I (km)
v01.1
r0
cos 45ˆ
Jcos 45ˆ
K
km s
Plot the variation of the six orbital elements
,
,
,
,
,
over the next two hours by integrating
Equations 12.84. Use MATLAB’s RKF integrator ode45.
Select the following spacecraft parameters:
Thrust = 0.5 kN.
Specific impulse = 300 s.
Initial mass = 2000 kg.