Ah, en effet ! J'ai inversé les fichiers « a » et « k » lors de la mise en ligne.
Désolé pour cette erreur et la confusion.
Je viens de remettre les bons fichiers au bon endroit.
Et juste pour confirmer les valeurs :
- celles-ci sont également données dans les fichiers « octave/matlab » donnés à la fin de P9 ; cela peut être utile pour déboguer (chercher « Invariants » ou « LA_ ») ;
- la valeur de la projection sur k du moment d'inertie au temps t = 0 est bien : 1.530393213572621e-01
obtenu :
- soit par calcul dans le repère d'inertie, puis changement de repère :
IA_G = IG_G + M*d*d * [ [ 1 0 0 ]; [ 0 1 0 ]; [ 0 0 0 ]; ] IA_G = 0.0544870 0.0000000 0.0000000 0.0000000 0.0544870 0.0000000 0.0000000 0.0000000 0.0029452 S = 1.000000 0.000000 0.000000 0.000000 0.866025 0.500000 0.000000 -0.500000 0.866025 w_G = 0 0 60 LA_G = IA_G * w_G LA_G = 0.000000 0.000000 0.176715 LA_O =S' * LA_G LA_O = 0.000000 -0.088357 0.153039
- soit directement dans le repère principal :
IA_O = S' * IA_G * S IA_O = 0.0544870 0.0000000 0.0000000 0.0000000 0.0416016 0.0223182 0.0000000 0.0223182 0.0158307 w_O = S' * w_G w_O = 0.000000 -30.000000 51.961524 LA_O = IA_O * w_O LA_O = 0.000000 -0.088357 0.153039