Examples
pos = [[7000.e3; 0; 0], [0; 7000.e3; 0], [0; 0; 7000.e3]]; CL_fo_zonHarmAcc(pos, 2) CL_fo_zonHarmAcc(pos, 2:6) CL_fo_zonHarmAcc(pos, [2,4,6]) |  |  |
Acceleration due to gravity (zonal harmonics)
[acc] = CL_fo_zonHarmAcc(pos, [nz, er, mu, j1jn])
Acceleration due to zonal harmonics.
Notes:
- nz is an array of the zonal term numbers to consider.
- The origin for the position vector must be the central body.
- The coordinates frame can be any frame whose Z axis is the revolution axis of the perturbation field.
- If nz contains numbers in the range [2,6], explicit formulas are used (faster).
- If nz is equal to -1 or "all", all the harmonics are used.
See Force models for more details.
Position vector from the central body [m]. (3xN)
(optional, integer) Zonal term numbers to considerer (1xP) or "all" or -1
(optional) Equatorial radius [m]. Default value is %CL_eqRad. (1x1)
(optional) Gravitational constant [m^3/s^2]. Default value is %CL_mu. (1x1)
(optional) Zonal harmonics. Default value is %CL_j1jn.
Acceleration [m/s^2]. (3xN)
CNES - DCT/SB
D. Vallado - Fundamentals of Astrodynamics and Applications, 2nd edition
pos = [[7000.e3; 0; 0], [0; 7000.e3; 0], [0; 0; 7000.e3]]; CL_fo_zonHarmAcc(pos, 2) CL_fo_zonHarmAcc(pos, 2:6) CL_fo_zonHarmAcc(pos, [2,4,6]) | | |