CL_mod_nutationMatrix — Nutation matrix and rotation vectors (IERS 1996)
[NUT[,OM[,OMP]]] = CL_mod_nutationMatrix(K,KP,KPP[,flag])
with EPS0 mean obliquity,DPSI longitude, DEPS obliquity
[EPS0 ; DPSI ; DEPS ] (3xN)
[EPS0P ; DPSIP ; DEPSP] (3xN)
[EPS0PP; DPSIPP; DEPSPP] (3xN)
(optional) "n" -> matrix only; "p" -> matrix and rotation vector; "s" -> matrix, rotation vector and its time derivative (default depends on number of outputs)
Nutation matrix (3x3xN)
(optional) Rotation vectors (3xN)
(optional) first time derivative of OM (3xN)
1 IERS Conventions (1996), Dennis D. McCarthy
2 Explanatory Supplement to the Astronomical Almanac, Seidelman (1992)
jj_tai = [19500:1:20500]; ss_tai = 43200*ones(jj_tai); [NUT,F,NUTP,FP,NUTPP,FPP]=CL_mod_nutationAngles(jj_tai,ss_tai,"s"); //calcul des angles de la nutation (dPsi,dEps) [eps0,eps0p,eps0pp]=CL_mod_meanObliquity(jj_tai,ss_tai,"s"); //calcul de l'obliquite moyenne K=[eps0; NUT]; Kp=[eps0p; NUTP]; Kpp=[eps0pp; NUTPP]; [NUT,OM,OMP] = CL_mod_nutationMatrix(K,Kp,Kpp,"s");