<< Orbit properties Orbit properties CL_kp_E2M >>

celestlab >> Orbit properties > CL_csl_defWalker

CL_csl_defWalker

Definition of a Walker constellation

Calling Sequence

[oe, nums] = CL_csl_defWalker(type_oe, oe_ref, T, P, F [[, opts]])

Description

Parameters

type_oe:

(string) Type of orbital elements: "kep" or "cir"

oe_ref:

Reference orbital elements = orbital elements for satellite 1 in plane 1 (6x1)

T:

(integer) Total number of satellites (1x1)

P:

(integer) Total number of planes (1x1)

F:

(integer) Phasing parameter, between 0 and P-1 (1x1)

opts:

(struct, optional) Computation options, see above for details

oe:

Resulting orbital elements (6xN)

nums:

(struct) Structure containing the fields: n = sat. number, p = plane number, sp = sat. number in plane (1xN)

Authors

Examples

// Reference orbital elements ("cir" type)
cir0 = [7078.e3; 0; 0; 1.71; 0; 0];
T = 12; // 12 satellites in all
P = 3; // 3 planes
F = 1; // Gap between 2 sat in 2 adjacent planes = 1 * (2 * pi / 12)

// Orbital elements for all satellites
[cir, nums] = CL_csl_defWalker("cir", cir0, T, P, F);

// Plot
function Plot()
scf();
plot(CL_rMod(cir(5,:) * 180 / %pi, 0, 360), CL_rMod(cir(6,:) * 180 / %pi, 0, 360), "o");
e = CL_g_select(gce(), "Polyline");
e.mark_background = color("blue");
e.mark_size = 8;
a = gca();
a.data_bounds = [-10, -10; 370, 370];
a.tight_limits = "on";
a.x_ticks = CL_g_ticks(0 : 120 : 360);
a.y_ticks = CL_g_ticks(0 : 30 : 360);
a.sub_ticks = [0,0];
xtitle(msprintf("T=%d, P=%d, F=%d", T, P, F), "RAAN (deg)", "Argument of latitude (deg)");
CL_g_stdaxes(a);
endfunction

Plot();

Report an issue
<< Orbit properties Orbit properties CL_kp_E2M >>