L’erreur dans les éphémérides vient-elle de la saisie, de la formule ou de la machine ?
Un même résultat doit passer par quatre niveaux : saisie, algorithme, fonctionnement interne et affichage final. Ne les confondez pas en une seule « précision ».
Algorithmes généraux et calcul numérique
Un programme peut afficher de nombreuses décimales, mais les temps de saisie peuvent avoir été arrondis, les formules peuvent avoir omis des termes et les opérations internes peuvent perdre des détails. Ce n’est qu’en démontant ces étapes une par une que vous pourrez savoir où les changements sont vraiment utiles ; le simple fait de changer le format de sortie ne peut généralement pas modifier les erreurs précédentes.
Devinez d'abord
Si le 125° interne mal calculé s'affiche comme 125,000000000°, la décimale manquante peut-elle être récupérée ?
Afficher la réponse et les explications
ne peut pas — L'allongement de chaîne ne restaure pas les informations perdues lors de la saisie, de l'algorithme ou de l'opération.
Ensuite, connectez l'expérience et la science
Les informations saisies sont-elles d'abord supprimées ?
Ne conserver que quelques décimales pour T équivaut à mettre d'abord en gras l'heure de la requête. Si j'utilise la formule exacte plus tard, je ne trouve pas l'entrée d'origine.
Erreur d'arrondi d'entrée ≠ Erreur d'approximation de formule
Quels termes sont omis dans la formule elle-même ?
L'omission active de T² et T³ modifiera les résultats. Il s'agit d'une différence causée par le choix de l'algorithme et ne peut pas être corrigée en ajoutant quelques zéros à la dernière ligne.
Comparez séparément la version complète sélectionnée et la version omise
Les machines ne peuvent conserver que des détails limités
Le grand angle est d'abord calculé à des centaines de milliers de degrés, puis normalisé, ce qui exposera des bits efficaces internes insuffisants. Lorsque le cos du petit angle est proche de 1, l'angle inclus peut être perdu si vous inversez directement le calcul.
Choisissez une voie de calcul numériquement stable
Décidez enfin comment afficher
Dans l'exemple de Neptune du livre original, l'angle complet est d'abord calculé puis arrondi pour obtenir 19° ; les deux entrées sont d'abord arrondies pour obtenir 20°. Le nombre de chiffres affichés doit correspondre à la précision de la saisie et du modèle.
Calculez d'abord, puis effectuez un arrondi final significatif
Les erreurs d'arrondi, de troncature arithmétique et de virgule flottante sont saisies séparément, puis les versions écrites et omises entièrement sélectionnées sont comparées en utilisant la même entrée. Afficher quelques décimales supplémentaires à la fin d’un résultat n’améliore pas automatiquement la précision d’un modèle scientifique.
essaie-le
Regardez d'abord "seulement 6 chiffres en interne", puis changez l'affichage à 9 chiffres - la réponse ne sera pas plus proche de la référence. Les conditions d'expansion font varier respectivement l'entrée, l'algorithme et la précision interne ; Les cas de petit angle et d'arrondi final sont des problèmes distincts.
Pas besoin de mémoriser les symboles
- T
- Siècle julien à partir de J2000 ; 0,4 équivaut à environ J2040.0
- bits internes valides
- La précision de l'opération décimale simulée dans cette leçon ne doit pas modifier le format matériel JS.
- Afficher les décimales
- Format de chaîne, ne change pas le processus de calcul
- acos / atan2
- Deux voies mathématiquement équivalentes mais numériquement différentes vers l'inversion aux petits angles
Algorithme, précision et champ d'application
Complétez le chapitre 2 sur les différences entre la précision requise pour les objectifs de calcul, la précision de la machine et la précision du rapport final pour éviter de répéter la soustraction des nombres les plus proches dans l'ancienne leçon de pesée. L'expérience 0 utilise le polynôme de longitude céleste moyenne lunaire de 218,3164477+481267,88123421T−0,0015786T²+0,0000019T³ ; il s'agit de la formule tronquée du chapitre 2, et non de la véritable longitude céleste lunaire complète. L'entrée de comparaison T est arrondie, les termes d'ordre supérieur sont activement omis, les chiffres décimaux significatifs sont simulés opération par opération et enfin quatre couches d'effets sont affichées. La soi-disant référence est le résultat mathématique à double précision de la même formule du livre, et non la vraie valeur mesurée ou la limite supérieure de la véritable erreur des éphémérides. L'expérience 1 utilise l'avertissement aux petits angles au chapitre 1/2 : atan2(sinθ, cosθ) est comparé à acos(cosθ) ; le premier ne peut toujours pas récupérer les informations d’entrée perdues. L'expérience 2 reproduit l'erreur d'arrondi prématuré de Neptune sur Sonneberg en 1996 et 1,4+1,4 ; L'expérience 3 reproduit le carré continu de 1,0000001 et utilise log1p/exp pour le vérifier indépendamment avec la valeur de référence dans le livre ; L'expérience 4 reproduit x←(9x+1)x−1, mais le point fixe mathématique 1/3 s'écartera dans les opérations en virgule flottante. La sortie spécifique au matériel du livre original n'est pas destinée à être un résultat fixe pour toutes les machines. Si l'itération dépasse 1e100 ou s'arrête dans un temps non limité, la gamme du modèle T est limitée à ±1 siècle ; voir également le cours d'illumination planétaire/magnitude pour le nombre de chiffres significatifs en magnitude.
vérification de ma compréhension
La formule complète est différente de la version omise. Quel est le niveau d'erreur principalement causé ?
approximation de l'algorithme — Les termes omis modifient directement le modèle de calcul mathématique.
Si cos a arrondi un très petit angle à 1, l’angle d’origine peut-il être restauré en utilisant acos ?
Non, vous devez conserver des informations plus valides ou modifier un itinéraire stable — La même entrée 1 ne peut pas distinguer s'il s'agit à l'origine d'un angle nul ou d'un très petit angle non nul.
Le 15,88 mars devrait-il simplement être arrondi au 16 mars ?
Non, la date appartient toujours au 15 — Cela signifie 0,88 jour après minuit le 15. La segmentation des étiquettes de date n'est pas la même opération que l'arrondi numérique final.
Identifiez d’abord à quel niveau se situe l’erreur, puis décidez s’il faut compléter les données, modifier la formule, changer la méthode numérique ou simplement ajuster le nombre de chiffres rapportés.
Jean Meeus · Astronomical Algorithms, 2nd ed. · Non. 2 chapitre
Algorithmes généraux et calcul numérique · Tous les cours interactifs