Ressources libres - Lumières sur l’Univers
Entrée du siteSommaireGlossairePage pour l'impression<-->
- Mécanique céleste, Temps et Calendriers

En savoir plus: l'interpolation des éphémérides

ensavoirplusEn savoir plus

Une éphéméride est une table qui fournit, pour des dates particuières en principe équidistantes, les valeurs d'une ou de plusieurs fonctions qui décrivent l'évolution des positions des astres concernés, ou d'autres grandeurs relatives à ces astres. On appelle interpolation la technique de calcul qui permet d'évaluer les valeurs de ces fonctions pour des dates intermédiaires autres que celles qui figurent dans l'éphéméride.

Le cas le plus simple est celui où la fonction considérée varie assez lentement, ou est donnée avec assez peu de précision pour qu'aucun calcul ne soit nécessaire. On se contente alors d'une interpolation à vue (par exemple, magnitude de la planète Jupiter : aucun calcul n'est nécessaire, quand on examine l'éphéméride du chapitre 5, pour estimer que sa valeur est -1,8 le 2 juillet 2002).

Plus souvent, on devra employer l'interpolation linéaire et procéder comme suit. Soit :

t, l'instant pour lequel on veut calculer le fonction f ;

t1, l'instant immédiatement inférieur ou égal à t et figurant dans l'éphéméride, la valeur correspondante de f étant f1 ;

t2, l'instant immédiatement supérieur à t et figurant dans l'éphéméride, la valeur correspondante de f étant f2.

Alors, on calculera f par la formule :

f = f_1 + n*subsup(Delta;1;1)

n=(t-t_1)/(t_2-t_1) et subsup(Delta;1;1)=f_2-f_1

s'appellent respectivement facteur d'interpolation et différence première.

Dans certains cas, l'interpolation linéaire fournit un résultat trop imprécis. Pour s'en rendre compte, on forme le tableau suivant :

images/formule.png

où t0, t1, t2, t3, etc sont des instants figurant dans l'éphéméride, et f0, f1, f2, f3, etc les valeurs correspondantes de la fonction f données par cette éphéméride. Les quantités :

subsup(Delta;k;1)=f_(k+1)-f_k

subsup(Delta;k;2)=subsup(Delta;k+1;1)-subsup(Delta;k;1) {k=0,1,2,...}

subsup(Delta;k;3)=subsup(Delta;k+1;2)-subsup(Delta;k;2)

s'appellent respectivement différences premières, différences secondes, différences troisièmes, etc. L'interpolation linéaire est suffisante si les différences secondes sont inférieures à 4 unités du dernier ordre envisagé. Sinon (et à condition que les différences troisièmes n'excèdent pas 62 unités du dernier ordre envisagé, ce qui sera toujours le cas dans cet ouvrage), on emploie la formule de Bessel :

f=f_1+n*subsup(Delta;1;1)-(1/4)*n(1-n)*(subsup(Delta;0;2)+subsup(Delta;1;2))

où n est le même facteur d'interpolation que plus haut, et, compte tenu des expressions ci-dessus :

subsup(Delta;0;2)+subsup(Delta;1;2)=subsup(Delta;2;1)-subsup(Delta;0;1)

Page précédentePage suivante