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

En savoir plus: approximation de Tchébychev (2)

ensavoirplusEn savoir plus

Soleil, Lune, Planètes

Soit t une date julienne appartenant à l'intervalle de temps t0, t0 + DT et y une coordonnée d'un astre. On calcule y(t) par la formule :

y=a_0+a_1*T_1*((x))+...+a_n*T_n*((x))

  • x est un paramètre compris entre -1 et +1 qui se déduit de t par la relation : x=-1+2*(t-t_0)/DT
  • a0, a1, ..., an sont les coefficients numériques publiés dans la Connaissance des Temps, valables pour l'intervalle de temps t0, t0 + DT contenant t :
  • T1(x), ..., Tn(x) sont des polynômes de Tchebychev, en x, que l'on calcule :
    • soit par la relation : (T_p)(x)=cos(p*arccos(x)) ;
    • soit par la relation de récurrence : (T_(P+1))(x)=2*x*(T_p)(x)-(T_(p-1))(x) Avec : (T_0)(x)=1; (T_1)(x)=x

Caractéristiques générales

Les éphémérides de l'année sont publiées corps par corps. Pour chaque corps, on compte un certain nombre de variables qui pour des intervalles de temps donnés, sont représentées par des tableaux donnant les coefficients de leurs développements en polynômes de Tchebychev ou sous forme de fonctions mixtes.

En haut de chaque page, on trouve :

  • le nom du corps concerné suivi de l'année de l'éphéméride ;
  • le type des variables représentées ;
  • éventuellement, les caractéristiques du système de coordonnées.

En bas de chaque page, on trouve les unités utilisées.

Les coefficients d'un tableau sont publiés en colonne. Pour chaque tableau, on précise :

  • le nom de la variable représentée ;
  • les dates de début et de fin de l'intervalle de temps correspondant au domaine de validité du tableau.
  • la longueur des intervalles de temps DT (en jours) figure en haut de chaque page ;
  • les coefficients sont notés 0, 1, ..., p,... n;
  • la valeur de contrôle du tableau, imprimée en italique, figure en tête de chaque tableau. C'est la valeur y(t0) de la variable y représentée par le tableau, pour t0 début de l'intervalle de temps t0, t0 + DT. Elle est égale à : y(t_0)=a_0-a_1+...+(-1)^n*a_n

La valeur de contrôle permet de vérifier que les coefficients ont été correctement entrés en mémoire sur une calculette ou un calculateur. De plus, l'ensemble des valeurs de contrôle des tableaux d'une même variable donne une idée de l'évolution de cette variable au cours de l'année.

Calcul d'une coordonnée

Pour calculer la valeur d'une coordonnée d'un astre pour une date t exprimée en UTC, on commence par :

  • effectuer, sur t, une correction déduite de la table 2 afin de se ramener à l'argument des éphémérides (identifié à TT) :
  • chercher, parmi les tableaux représentant la coordonnée, celui qui correspond à l'intervalle de temps t0, t0 + DT dans lequel se trouve t. Il faut insister sur le fait que la représentation, qu'elle soit en polynômes de Tchebychev ou sous forme de fonctions mixtes, n'est pas valable en dehors de l'intervalle t0, t0 + DT.

Le calcul se poursuit de la manière suivante :

Soleil, Lune, Planètes

  • On calcule x par la formule (2).
  • A partir des n + 1 coefficients du tableau a0, a1, ..., ap, ..., an (notés 0, 1,... p, ..., n), on calcule la coordonnée y(t) par la formule (1).
  • Les polynômes de Tchebychev tp(x) se calculent soit par la formule (3), soit par la relation de récurrence (4). Notons que, si on n'utilise que les p + 1 (p < n) premiers coefficients du développement en polynômes de Tchebychev, on commet, sur la représentation de y, une erreur au plus égale à la somme des valeurs absolues des n - p coefficients négligés.
  • Pour avoir le temps de passage du Soleil au méridien origine en un jour donné, numéroté j dans le mois, on applique la même méthode en prenant t = j + 0,5.
Page précédentePage suivante