Начальные условия

Начальные условия для численного интегрирования осреднённых уравнений движения были получены с помощью данных эфемериды DE406/LE406.

На момент 1721423.5 юлианских дней (полночь 1 января 1 года нашей эры) значения барицентрических координат и скоростей семи больших планет, Солнца и центра масс Земля-Луна, заданные в системе стандартного экватора эпохи J2000.0, были преобразованы в оскулирующие кеплеровские элементы орбиты.

Оскулирующие элементы вычислены в гелиоцентрической эклиптической системе отсчёта. Положения эклиптики и точки весеннего равноденствия соответствуют опорной системе отсчёта J2000.0.

Оскулирующие элементы были преобразованы в совокупность сглаженных параметров. Сглаженные параметры послужили начальными условиями для численного интегрирования осреднённых уравнений движения - системы 48 дифференциальных уравнений первого порядка.

Численное интегрирование было выполнено методом Эверхарта 15 порядка. Был выбран постоянный шаг интегрирования, равный 100 годам. Использование меньшего шага приводило к увеличению времени расчётов, а оценки сглаженных элементов орбит планет с точностью до 15 знаков после запятой были равны этим же величинам, вычисленным с шагом 100 лет.

Здесь можно вернуться к возмущающей функции

Здесь можно вернуться на страницу Программные приложения