Code covered by the BSD License  

Highlights from
Optimal Impulsive Orbital Transfer

Optimal Impulsive Orbital Transfer



27 Nov 2012 (Updated )

MATLAB script for the solution of the one and two impulse orbit transfer between two Earth orbits.

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

% primer vector initialization

% 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 = stm12 \ (pvf' - stm11 * pvi');

Contact us