Computational Astronomy & Ephemeris Precision (天文計算與星曆精度)
At its core, QZSY is an applied astronomical system — its accuracy depends entirely on the precision of the planetary positions it uses. This module examines the mathematics behind planetary position calculation, the evolution of ephemeris precision from ancient China to modern Swiss Ephemeris, and the practical implications of computational accuracy for QZSY chart interpretation.
Fundamentals of Planetary Position Calculation
Calculating a planet's position at a given moment requires solving the equations of its orbital motion around the Sun (or around the Earth, for geocentric coordinates). The key parameters are:
- Orbital elements: Semi-major axis, eccentricity, inclination, longitude of ascending node, argument of perihelion, and mean anomaly at epoch.
- Perturbation corrections: Gravitational influences from other planets that deviate the orbit from a simple ellipse.
- Coordinate transformation: Converting heliocentric (Sun-centered) positions to geocentric (Earth-centered) positions as seen from the observer's location.
Chinese Classical Ephemeris Methods
Chinese astronomers developed increasingly sophisticated calculation methods over the centuries:
| Period | Method | Approximate Accuracy |
|---|---|---|
| Han Dynasty | San Tong Li (三統曆) — basic cycle-based | ~1-2° for outer planets |
| Tang Dynasty | Da Yan Li (大衍曆) — Indian-influenced refinement | ~0.5-1° for outer planets |
| Yuan Dynasty | Shou Shi Li (授時曆) — most advanced pre-modern Chinese | ~0.1-0.5° |
| Qing Dynasty | Shi Xian Li (時憲曆) — Jesuit-enhanced | ~0.05-0.1° |
| Modern | Swiss Ephemeris (JPL DE431) | ~0.001 arcsecond |
Precession and Ayanamsa
The precession of equinoxes (歲差, Sui Cha) is the slow westward drift of the equinox points relative to the fixed stars, completing a full cycle in approximately 25,772 years. This phenomenon creates the angular offset between the tropical (equinox-based) and sidereal (star-fixed) coordinate systems — the ayanamsa. In the current era, the ayanamsa is approximately 24 degrees, meaning a planet at 0° Aries in tropical coordinates is at approximately 6° Pisces in sidereal coordinates.
Chinese astronomers discovered precession independently — the astronomer Yu Xi (虞喜) documented it in 336 CE, approximately two centuries after Hipparchus' discovery in the West. Classical Chinese ephemerides incorporated various precession rates, some more accurate than others.
Error Analysis: Classical vs Modern
When reconstructing historical QZSY charts using modern ephemeris data, practitioners may find that classical-era calculations contained errors of several degrees — enough to place a planet in the wrong Lunar Mansion. This has implications for interpreting historical case studies: some classical 'validated' readings may have been based on incorrect planetary positions. Modern recalculation can both validate and challenge classical interpretations, contributing to the ongoing refinement of QZSY methodology.
Building a QZSY Ephemeris Engine
A modern QZSY ephemeris engine requires: (1) Swiss Ephemeris library integration (C library with Python/JavaScript wrappers), (2) Configurable ayanamsa for sidereal conversion, (3) Lunar mansion mapping tables (28 mansions with current-epoch star boundaries), (4) Ci/house calculation from Ascendant, (5) Four Remainders computation routines, and (6) Aspect and dignity assessment algorithms. The platform's apps/calc-engine Flask service implements this full pipeline.