Code covered by the BSD License  

Highlights from
Bi-elliptic Transfer Between Coplanar Circular Orbits

Bi-elliptic Transfer Between Coplanar Circular Orbits

by

 

02 Nov 2012 (Updated )

Bi-elliptic orbital transfer between two coplanar circular Earth orbits.

pviniz (tof, r1, v1, dv1, dv2)
function pviniz (tof, r1, v1, dv1, dv2)

% primer vector initialization

% required by primer.m

% Orbital Mechanics with MATLAB

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

global mu pvi pvdi

% compute primer vector at first impulse
  
dv1m = norm(dv1);

pvi = dv1 / dv1m;

% compute primer vector at second impulse
  
dv2m = norm(dv2);

pvf = dv2 / dv2m;
 
% compute state transition matrix

[r2, v2, stm] = stm2(mu, tof, r1, v1);

% extract submatrices of state transition matrix

stm11(1:3, 1:3) = stm(1:3, 1:3);

stm12(1:3, 1:3) = stm(1:3, 4:6);

% compute initial value of primer derivative vector

pvdi = inv(stm12) * (pvf' - stm11 * pvi');

Contact us