Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2022. Т. 26, № 4. С. 738-763 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1888
EDN: DVJGRH
УДК 521.182
Сравнение элементов орбит больших планет, Луны и Солнца с использованием различных математических моделей на интервале времени с 1600 по 2200 гг.
А. Ф. Заусаев, М. А. Романюк
Самарский государственный технический университет,
Россия, 443100, Самара, ул. Молодогвардейская, 244.
Аннотация
Проведен анализ точности элементов орбит, полученных по координатам и компонентам скоростей, найденных с помощью коэффициентов многочленов Чебышева планетного каталога DE405. Для исследования на интервале времени с 1600 по 2200 гг. проведено сопоставление элементов орбит, найденных с помощью каталога DE405, а также полученных на основании численного интегрирования уравнений движения, основанных на взаимодействии движущихся материальных тел с окружающим пространством. На примере численного интегрирования уравнений движения Луны показано преимущество использования уравнений движения, основанных на взаимодействии движущихся материальных тел с окружающим пространством по сравнению с релятивистскими уравнениями. На основании сравнения элементов орбит Меркурия, полученных путем решения уравнений, основанных на взаимодействии движущихся материальных тел с окружающим пространством и найденных с помощью использования каталога DE405, показано, что на исследуемом интервале времени элементы орбит практически совпадают. Максимальное расхождение в средней аномалии на конце интервала интегрирования составляет менее 1'' (секунды дуги). Определены невязки вековых смещений перигелиев для Меркурия, Венеры, Земли + Луны и Марса, значения которых для DE405 соответственно равны: 43.08'', 8.4'', 3.83'' и 1.14''. Показано, что погрешности вековых смещений перигелиев планет Меркурия, Венеры, барицентра Земли + Луны и Марса, полученные
Математическое моделирование, численные методы и комплексы программ Научная статья
© Коллектив авторов, 2022 © СамГТУ, 2022 (составление, дизайн, макет)
Q ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Заусаев А. Ф., Романюк М. А. Сравнение элементов орбит больших планет, Луны и Солнца с использованием различных математических моделей на интервале времени с 1600 по 2200 гг. // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2022. Т. 26, № 4. С. 738-763. EDN: DVJGRH. DOI: 10.14498/vsgtu1888. Сведения об авторах
Анатолий Федорович Заусаев © https://orcid.org/0000-0002-5035-9615 доктор физико-математических наук; профессор; каф. прикладной математики и информатики; e-mail: zausaev_af@mail.ru
Мария Анатольевна Романюк А https://orcid.org/0000-0003-1352-6725 кандидат технических наук; доцент; каф. прикладной математики и информатики; e-mail: zausmasha@mail.ru
Сравнение элементов орбит больших планет, Луны и Солнца.
при использовании каталога ББ405, принимают следующие значения: 0'', 6.06'', 3.83'' и 1.08''. Для внешних планет: Юпитера, Сатурна, Урана, Нептуна и карликовой планеты Плутон на основании рассмотренных сравнений различных уравнений движения расхождений элементов орбит не обнаружено. На основании проведенных исследований показано, что использование гармонических координат в релятивистских уравнениях при создании каталога ББ405 оправдано только для Меркурия и внешних планет: Юпитера, Сатурна, Урана, Нептуна и карликовой планеты Плутон.
Ключевые слова: элементы орбит, численное интегрирование, дифференциальные уравнения движения.
Получение: 19 октября 2021 г. / Исправление: 15 октября 2022 г. / Принятие: 17 ноября 2022 г. / Публикация онлайн: 6 декабря 2022 г.
Исследование движения больших планет, Луны и Солнца сопряжено с большим объемом вычислений, точность которых зависит от выбора как математической модели движения, так и метода решения задачи. Математическая модель движения небесных тел, как правило, описывается системой дифференциальных уравнений второго порядка. В основе каждой математической модели лежат определенные физические предпосылки. Важную роль в механике играет понятие пространства.
Понятие об эфире исходит из глубокой древности. Две с половиной тысячи лет назад ученые древней Греции сформулировали и развили понятие об эфире. Эфирная концепция достигла своей кульминации в XIX веке, когда Максвелл, опираясь на созданную им модель эфира, получил фундаментальные уравнения электродинамики [1].
В теории тяготения Ньютона под пространством понимается плоское евклидово пространство. Механика Ньютона основана на законе всемирного тяготения и трех законах движения [2,3].
В основе релятивистской механики лежат два основных принципа: принцип эквивалентности и принцип общековариантности [4, 5]. Принцип эквивалентности говорит о том, что все физические процессы протекают одинаково в инерциальной системе, находящейся в однородном поле тяготения и в неинерциальной равномерно ускоренной системе. Из принципа общеко-вариантности следует, что уравнения должны иметь одну и ту же форму во всех без исключения системах отсчета. Опираясь на эти принципы, Эйнштейн пришел к выводу, что пространство событий общей теории относительности должно представлять собой псевдориманово четырехмерное пространство с метрикой [6]
йв2 = дар <1ха(!хР. (1)
Движение пробной частицы под действием тяготения — это свободное движение по инерции, происходящее по геодезическим линиям псевдориманова пространства, метрика которого формируется самими гравитационными массами.
С развитием квантовой теории появился термин — физический вакуум. Физический вакуум является особым видом материи, претендующим на первооснову мира. Физический вакуум стал предметом изучения физики бла-
годаря усилиям известных ученых: П. Дирака, Р. Фейнмана, Дж. Уилера, У. Лэмба, де Ситтера, Г. Казимира, Я. Б. Зельдовича [7-10]. Столь разнообразное представление от эфира до физического вакуума об окружающем нас пространстве говорит о сложности проблемы, стоящей перед наукой.
Следует отметить, что в настоящее время свойства окружающего нас пространства недостаточно изучены, поэтому при выводе дифференциальных уравнений движения материальных тел необходимо наделять его определенными динамическими свойствами.
Ранее нами получены дифференциальные уравнения движения небесных тел, основанные на взаимодействии материальных тел с окружающим пространством, которые имеют следующий вид [11-15]:
(й2Х (И2
<12У
£( ¥)
ЗаогГ2
0г
(И2
й2г (И2
£
г
Е
( Уг - У
А2 + Аг 3Д - ^ + 3(Д? - Г&)2
Заогг2
Ог
V Д УД2 +Дг 3 Д? - г0г + 3 (Д? - Г&)2
( 2? - 2 \
V Дг )
(2)
3аогг\
Ог
Д2 + Дг ^Д? - г0г + 3(Д? - г0г)2 1
где гог — эффективный радиус г-того тела; аог — соответствующее ускорение для г-того тела на расстоянии гог от центра массы; X, У, Z — барицентрические координаты возмущаемого тела; X?, Уг, ^г — барицентрические координаты возмущающих тел; Д? — взаимное расстояние между возмущаемым и возмущающим телом.
Дифференциальные уравнения движения в ньютоновой форме имеют следующий вид [1,2]:
^ X \—^ т о / X? — X ^
^ У <И2 <!2 г
г
Е *2
т?
т,
У^ к2тг
v д? ( - ¥
I Д?
г? - г ДО
(3)
где Д2 = (Хг - X)2 + (У? - У)2 + (г? - г)2; X, У, Я — барицентрические координаты возмущаемого тела; т?, X?, Уг, — массы и барицентрические координаты возмущающих тел.
Целью данной работы является оценка точности координат планет, представленных в планетном каталоге БЕ405, поэтому следует привести математическую модель в форме дифференциальных уравнений, которую автор использовал при создании БЕ405.
Дифференциальные уравнения движения в барицентрической системе координат с учетом ньютоновых и шварцшильдовских членов и формулы (1) имеют следующий вид [16]:
£
3=г
^ (rj - П)
1 -
2(£ + 7)^- № 2/3 - 1
г3
£
к=г
Пк
£ +Ч!) ■
+
2
2
с
с
Сравнение элементов орбит больших планет, Луны и Солнца.
+ (1+7)(?) -
v3 \ 2 2(1 + 7)
(n - rj)гг
+ Ъ? ( Гз - ri)j>j\ +
+ 72 Е ЩГ( ri - rJ) [(2 + 27)h - (1 + 27) fj] ( h - fj) +
j=i ij
+
3 + 47 Щj fj IT2
EM'j + y^
f.. Z-—/
ч'-^"' ij m=1
ym( fm П)
(4)
j=i
где Гг, т%, Гг — координаты, скорости и ускорения в барицентрической системе координат г-того тела; ^^ = k2mj; к2 — гравитационная постоянная; т^ — масса ^'-того тела; г^ = |rj — гг|; [ и 7 — релятивистские параметры, [ = 7 = 1; уг = |Тг1; с —скорость света.
При создании эфемерид Луны помимо гравитационных и релятивистских эффектов учитывалось влияние фигур Земли и Луны в математической модели. Ускорение Луны благодаря учету зональных и тессеральных гармоник в координатной системе £]( имеет вид [16]
'
{П1
!>( Щ У
п=1
(п + l)pn(sin^) 0
- cos вPn (sin^>)
+
+
П2 п
£(Щ)n £
п=1
m=1
-(п + 1)Pnm(sin^)[C, т sec ^Pnm(sin ф)[-С cos ^Pnm (sin ^>)[С,
nm cos т\ + Snm sin rnA] m sin mA + Snm cos mA] nm cos mA + Snm sin mA]
(5)
где у — гравитационная постоянная; г — расстояние между центрами масс двух тел; п\ и П2 — максимальные степени зональных и тессеральных гармоник несферичных тел соответственно; Рп^тр) — полином Лежандра степени п; РЩг(8\пр) —присоединенный полином Лежандра степени п и порядка т; Зп — зональные гармоники от несферичности тела; Спт, Зпт — коэффициенты тессеральных гармоник; р — широта притягиваемого тела в фиксированной системе координат £; А — восточная долгота притягиваемого тела. Вклад в инерциальное ускорение от несферичного тела возникает от взаимодействия ее собственной фигуры с внешней точечной массой, представленной в координатной системе £ ] (. Ось £ направлена вовне от несферичного тела к точечной массе. Ось ] направлена на восток, лежит в селенографической плоскости ХУ, перпендикулярна оси £. Ось ( направлена на север, образуя правую систему координат.
Кроме того, земные приливы оказывают на геоцентрическое ускорение Луны следующее воздействие [16]:
'm —
' lm
(1 + Vm )
x + y5 у - x5
(6)
где k2 — число Лява; щ — радиус Земли; rim — геоцентрическое расстояние Луны; x, у, z, —декартовы геоцентрические координаты Луны; уm — гравитационная постоянная, умноженная на массу Луны; щ — гравитационная постоянная, умноженная на массу Земли; 5 — фазовый угол.
3
2
2
3
2
5
3
В отличие от уравнений (4), дифференциальные уравнения (2) записаны в декартовой системе координат. В них отсутствуют понятия силы и массы, которые, по мнению Пуанкаре, не имеют точного определения [17].
Ранее проведенные исследования показали, что решение дифференциальных уравнений, основанных на взаимодействии движущихся материальных тел с окружающим пространством, обладает преимуществом по сравнению с решениями Ньтоновых и релятивистских уравнений [18]. Использование уравнений в форме Ньютона для численного интегрирования уравнений движения больших планет приводит к накоплению ошибок в координатах внутренних планет. При построении численной теории Луны автором DE405 наряду с уравнением (4) решаются уравнения (5) и (6), что значительно усложняет решение задачи [19].
В отличие от релятивистских уравнений (4), (5) и (6), нами получены координаты всех больших планет, Луны и Солнца на определенные моменты на интервале времени с помощью решения одной системы дифференциальных уравнений (2).
Координаты планет, Луны и Солнца хранятся в каталоге DE405 в форме коэффициентов полиномов Чебышева [20-22]. Для различных планет степени многочленов различны. Максимальная степень многочленов порядка тринадцать обеспечивает необходимую точность для нахождения координат внутренних планет и Луны. Коэффициенты полиномов Чебышева хранятся в каталоге DE405 с шагом тридцать два дня. Многочлены Чебышева обеспечивают гладкое приближение даже при использовании многочленов высоких порядков. Многочлен Чебышева порядка п (для |ж| ^ 1) определяется как [21] Тп(х) = cos(n ■ arccos х).
Из равенства cos(a + ft) + cos(a — ft) = 2 cos a cos ft следует, что
cos(n + 1)^> = 2 cos п<р cos <р — cos(n — 1)^>. (7)
Если p = arccos x, то после подстановки в (7) получим следующую рекуррентную формулу:
Гп+1(ж) = 2хТп(х) — Тп-г(х) при п ^ 1.
Многочлены Чебышева до 13-го порядка имеют следующий вид:
То(х) = 1;
Т1(х) = х;
Т2(х) = 2х2 — 1;
Т3 (х) = 4ж3 — 3х;
ТА(х) = 8ж4 — 8х2 + 1;
Т5(х) = 16ж5 — 20ж3 + 5х;
Тб(х) = 32ж6 — 48ж4 + 18ж2 — 1;
Тт(х) = 64ж7 — 112ж5 + 56ж3 — 7х;
Т8(х) = 128ж8 — 256ж6 + 160ж4 — 32ж2 + 1;
Тд(х) = 256ж9 — 576ж7 + 432ж5 — 120ж3 + 9х;
Tw(x) = 512ж10 — 1280ж8 + 1120ж6 — 400ж4 + 50ж2 — 1;
Тц(х) = 1024Ж11 — 2816ж9 + 2816ж7 — 1232ж5 + 220ж3 — 11ж;
Ti2(x) = 2048ж12 — 6144ж10 + 6912ж8 — 3584ж6 + 840ж4 — 72х2 + 1;
Т1 з (х) = 4096ж13 — 13352Ж11 + 16640ж9 — 9984ж7 + 2912ж5 — 364ж3 + 13ж.
Сравнение элементов орбит больших планет, Луны и Солнца...
Рассмотрим особенности многочлена Чебышева на отрезке [—1, 1]. Внутри этого отрезка Тп(х) обращается в нуль в п точках:
Тп(х) = 0 при х = ео8(^(к + 1/2)/п),
где к = 0,1,... ,п — 1.
При |х| ^ 1 значение |Тп(х)| не может превышать 1.
Чтобы аппроксимировать произвольную функцию /(х) в интервале [ а, Ь], заменим независимую переменную х нормированной переменной х в интервале [—1,1] соотношением
2х — (а + Ь)
х =--- при х е 1а, Ь\ ^ х е 1—1,11
о — а
и
х = х ■ - (Ь — а) + -(а + Ь) при х е [—1,1] ^ х е [а, Ь].
Функция /(х) теперь может быть представлена многочленами Чебышева вплоть до п-го порядка в виде
п
/(х) = Е (х). j=0
Коэффициенты Cj вычисляются следующим образом:
2 п
* = ^л Е /(хп+1)т х+1), + к=0
где хп+1 представляет к-тый корень Тп+1.
Коэффициенты многочленов Чебышева в каталоге БЕ405 для каждого 32-дневного интервала приведены на начало интервала, что соответствует значению х = —1 для интервала [—1,1]. В точке х = —1 многочлены Чебышева примут следующие значения: 1, —1, 1, —1, ..., (—1)г+1, где г = 0,1, 2, 3,....
Значения координат планет Луны и Солнца в точке х = —1 вычисляются с помощью следующего соотношения:
/(х) = Е (—1)^
1-
j=0
Для нахождения скоростей планет Луны и Солнца в точке х = —1 берутся производные от многочленов Чебышева, затем вычисляются коэффициенты путем подстановки х = —1 в полученные многочлены. Сумма произведений полученных коэффициентов на производные от многочленов, умноженных для каждого объекта на постоянный множитель, определяет скорость исследуемых объектов. Постоянные коэффициенты в точке х = —1 имеют следующие значения: и = 0, и = 1, И = —4, Из = 9, и = —16, и5 = 25, И = —36, и7 = 49, и8 = —64, Ид = 81, и 10 = —100, ип = 121, и 12 = —144, И13 = 169.
Компоненты скоростей планет, Луны и Солнца находятся из соотношений
п
М (х) = Сг £ с, и, (х), 3=о
где С? для каждой планеты имеет постоянное значение. Для Меркурия, Венеры, барицентра Земли + Луны, Марса, Юпитера, Сатурна, Урана, Нептуна, Плутона значения С? взяты из [23].
Для сравнения координат и элементов орбит планет, вычисленных путем решения уравнений (2) и с помощью каталога БЕ405, весь интервал с 9 декабря 1599 г. (ЛБ 2305424.5) по 31 декабря 2199 г. (ЛБ 2524592.5) был разбит на шесть подынтервалов с шагом приблизительно 100 лет, в начале каждого из которых производилось сравнение.
По начальным данным гелиоцентрических координат планет были вычислены элементы их орбит, которые представлены в табл. 1.
В табл. 2-5 представлены элементы орбит Меркурия, Венеры, Земли + Луны и Марса на различные моменты времени, полученные с использованием БЕ405, численного интегрирования дифференциальных уравнений (2) и путем решения уравнений движения в Ньютоновой форме (3). Элементы орбит в табл. 2-5 размещены в следующем порядке: в первых строках табл. 2-5 представлены элементы орбит Меркурия-Марса, найденные с помощью каталога БЕ405, во вторых строках — элементы орбит, полученные на основании решения уравнений (2), и в третьих строках — найденные с помощью решения уравнений в Ньютоновой форме (3).
В табл. 6-9 представлены разности элементов орбит Меркурия, Венеры, Земли + Луны и Марса, полученные с помощью решения уравнений (2), (3) и БЕ405.
Меркурий является ближайшей к Солнцу планетой. Солнце оказывает на его движение более существенное влияние по сравнению с другими большими планетами. В средине прошлого столетия французский математик и астроном Леверье доказал, что с помощью решения уравнений движения в Ньютоновой форме невозможно построить теорию движения Меркурия, согласованную с наблюдениями [24]. Создание общей теории относительности и разработанные на ее основе дифференциальные уравнения движения небесных тел позволили согласовать движение Меркурия с наблюдениями. Использование гармонической системы координат полностью обосновало невязки между наблюдениями и результатами Ньютоновой гравитационной теорией.
Как видно из табл. 6, сопоставление элементов орбит Меркурия на интервале времен с 1600 по 2200 гг., полученных с использованием решений уравнений (2) и БЕ405, указывает на практическое их совпадение. Максимальное расхождение имеет место в средней аномалии, что составляет -0.0002° (градуса). Данное различие составляет 0.72'', в то время как аргументы перигелиев различаются на 0.36''. Учитывая, что данное расхождение в элементах орбит Меркурия произошло через 350 лет от начального момента интегрирования, полученные результаты можно считать вполне удовлетворительными. По данным табл. 2 найдено вековое смещение перигелия Меркурия, равное 43.08''.
Орбита Венеры расположена дальше орбиты Меркурия, поэтому она подвержена значительно меньшему воздействию Солнца и является объектом
Таблица 1
Элементы орбит больших планет на дату 27 июня 1969 г. (юлианская дата 2440400.5)
[Orbital elements of the major planets as of June 27, 1969 (JD 2440400.5)]
Planets M (in deg.) a (in au) e ш (in deg.) П (in deg.) г (in deg.)
Mercury 287.7773 0.387099 0.205617 29.0367 48.3688 7.0068
Venus 195.5860 0.723330 0.006813 54.6870 76.7621 3.3949
Earth + Moon 173.6153 1.000002 0.016713 107.0781 355.6508 0.0039
Mars 299.3763 1.523648 0.093379 107.0781 49.6476 1.852301
Jupiter 174.2869 5.207621 0.047265 273.7981 100.4002 1.3049
Saturn 302.9131 9.522000 0.054053 340.1900 113.7871 2.4866
Uranus 10.5307 19.280834 0.051387 97.9204 74.0516 0.7733
Neptune 185.6073 30.176592 0.004906 280.7089 131.7369 1.7702
Pluto 331.3768 39.774492 0.253318 112.9821 110.2863 17.1349
Элементы орбит Меркурия [Elements of the orbits of the Mercury]
Current, date Data sources (calculated by) M (in deg.) a (in au) e to (in deg.) П (in deg.) г (in deg.)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 160.8368 160.8366 160.8507 0.387099 0.387099 0.387099 0.205541 0.205541 0.205541 27.9902 27.9903 28.0343 48.8311 48.8311 48.8312 7.0287 7.0287 7.0287
1699 12 28 JD 2341968.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 311.1034 311.1033 311.1135 0.387898 0.387898 0.387098 0.205569 0.205569 0.205569 28.2774 28.2774 28.3096 48.7065 48.7065 48.7065 7.0228 7.0228 7.0228
1800 1 17 JD 2378512.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 101.3785 101.3784 101.3849 0.387098 0.387098 0.387098 0.205596 0.205596 0.205596 28.5579 28.5580 28.5782 48.5815 48.5815 48.5815 7.0169 7.0169 7.0169
1899 12 4 JD 2414992.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 349.7370 349.7370 349.7397 0.387098 0.387098 0.387098 0.205625 0.205625 0.205625 28.8405 28.8405 28.8488 48.4567 48.4567 48.4567 7.0110 7.0110 7.0110
1999 12 24 JD 2451536.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 140.0110 140.0110 140.0098 0.387098 0.387098 0.387098 0.205630 0.205629 0.205629 29.1243 29.1243 29.1206 48.4305 48.4305 48.4305 7.0050 7.0050 7.0050
2099 12 11 JD 2488048.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 159.3294 159.3294 159.3245 0.387097 0.387098 0.387098 0.205653 0.205653 0.205652 29.4069 29.4069 29.3913 48.2052 48.2052 48.2052 6.9991 6.9991 6.9991
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 309.5923 309.5924 309.5836 0.387098 0.387098 0.387098 0.205667 0.205667 0.205667 29.6962 29.6961 29.6686 48.0794 48.0794 48.0794 6.9930 6.9930 6.9930
Элементы орбит Венеры [Elements of the orbits of the Venus]
Current, date Data sources (calculated by) M (in deg.) a (in au) e to (in deg.) П (in deg.) г (in deg.)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 306.3409 306.3347 306.3354 0.723344 0.723344 0.723344 0.006972 0.006972 0.006972 53.7427 53.7490 53.7514 77.7910 77.7910 77.7910 3.3976 3.3976 3.3976
1699 12 28 JD 2341968.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 174.8557 174.8510 174.8516 0.723336 0.723336 0.723336 0.006927 0.006927 0.006927 53.7582 53.7630 53.7646 77.5144 77.5145 77.5145 3.3969 3.3969 3.3969
1800 1 17 JD 2378512.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 42.8101 42.8069 42.8073 0.723325 0.723325 0.723325 0.006833 0.006833 0.006833 54.3392 54.3424 54.3434 77.2331 77.2331 77.2331 3.3963 3.3963 3.3963
1899 12 4 JD 2414992.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 168.7702 168.7689 168.7691 0.723329 0.723329 0.723329 0.006802 0.006802 0.006802 54.3748 54.3761 54.3765 76.9589 76.9589 76.9589 3.3956 3.3956 3.3956
1999 12 24 JD 2451536.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 36.4937 36.4941 36.4941 0.723327 0.723327 0.723327 0.006756 0.006756 0.006756 55.1887 55.1884 55.1881 76.6784 76.6784 76.6784 3.3946 3.3946 3.3946
2099 12 11 JD 2488048.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 213.8268 213.8293 213.8290 0.723329 0.723329 0.723329 0.006693 0.006693 0.006693 55.1227 55.1202 55.1194 76.4005 76.4005 76.4005 3.3937 3.3937 3.3937
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 81.7762 81.7799 81.7795 0.723328 0.723328 0.723328 0.006703 0.006703 0.006703 55.7055 55.7017 55.7002 76.1244 76.1244 76.1244 3.3928 3.3928 3.3928
-J
oo
Элементы орбит Земли + Луны
Elements of the orbits of the Earth & Moon]
Current, date Data sources (calculated by) M (in deg.) a (in au) e ш (in deg.) П (in deg.) г (in deg.)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 341.0999 341.0969 341.0971 0.999994 0.999994 0.999994 0.016843 0.016843 0.016843 105.9312 105.9354 105.9360 355.7385 355.7376 355.7376 0.0524 0.0524 0.0524
1699 12 28 JD 2341968.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 358.8414 358.8391 358.8393 1.000002 1.000002 1.000002 0.016853 0.016853 0.016853 106.4499 106.4533 106.4537 355.5786 355.5776 355.5777 0.0393 0.0393 0.0393
1800 1 17 JD 2378512.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 16.6266 16.6225 16.6252 1.000021 1.000021 1.000021 0.016815 0.016815 0.016815 107.1480 107.1504 107.1507 355.1933 355.1924 355.1924 0.0260 0.0260 0.0260
1899 12 4 JD 2414992.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 331.4083 331.4078 331.4078 1.000006 1.000006 1.000006 0.016773 0.016773 0.016773 107.7580 107.7595 107.7597 354.8222 354.8212 354.8212 0.0131 0.0131 0.0131
1999 12 24 JD 2451536.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 349.1741 349.1744 349.1744 1.000004 1.000004 1.000004 0.016709 0.016709 0.016709 323.9007 323.8567 323.8561 139.0106 139.0542 139.0548 0.0001 0.0001 0.0001
2099 12 11 JD 2488048.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 335.4455 335.4466 335.4465 1.000001 1.000001 1.000001 0.016665 0.016665 0.016665 288.7857 288.7852 288.7850 174.4148 174.4140 174.4141 0.0131 0.0131 0.0131
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 353.1145 353.1164 353.1163 0.999996 0.999996 0.999996 0.016626 0.016626 0.016626 289.2209 289.2197 289.2194 174.4119 174.4109 174.4110 0.0260 0.0260 0.0260
Элементы орбит Марса
Elements of the orbits of the Mars]
Current, date Data sources (calculated by) M (in deg.) a (in au) e ш (in deg.) П (in deg.) г (in deg.)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 129.2408 129.2394 129.2394 1.523749 1.523749 1.523749 0.093011 0.093011 0.093011 283.5284 283.5297 283.5298 50.7286 50.7286 50.7286 1.8819 1.8819 1.8819
1699 12 28 JD 2341968.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 198.9886 198.9876 198.9876 1.523632 1.523632 1.523632 0.093113 0.093011 0.093113 284.3246 284.3255 284.3255 50.4397 50.4398 50.4398 1.8739 1.8739 1.8739
1800 1 17 JD 2378512.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 268.8961 268.8955 268.8955 1.523698 1.523698 1.523698 0.093296 0.093296 0.093296 284.9704 284.9709 284.9710 50.1466 50.1466 50.1466 1.8658 1.8658 1.8658
1899 12 4 JD 2414992.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 305.1031 305.1028 305.1028 1.523683 1.523683 1.523683 0.093215 0.093215 0.093215 285.7719 285.7721 285.7722 49.8552 49.8552 49.8552 1.8578 1.8578 1.8578
1999 12 24 JD 2451536.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 14.9027 14.9029 14.9029 1.523666 1.523666 1.523666 0.093307 0.093307 0.093307 286.5371 286.5370 286.5370 49.5619 49.5619 49.5619 1.8499 1.8499 1.8499
2099 12 11 JD 2488048.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 67.9179 67.9185 67.9184 1.523690 1.523690 1.523690 0.093613 0.093612 0.093612 287.2980 287.2975 287.975 49.2578 49.2578 49.2578 1.8417 1.8417 1.8417
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 137.8190 137.8199 137.8189 1.523711 1.523711 1.523711 0.093507 0.093507 0.093507 287.9596 287.9588 287.9587 48.9631 48.9631 48.9631 1.8333 1.8333 1.8333
о
более легким для исследования движения. Однако, как следует из сравнения элементов орбит, полученных путем численного интегрирования дифференциальных уравнений (2) и с использованием БЕ405, приведенных в табл. 7, расхождение средних аномалий и аргументов перигелиев значительные. Как известно, величина невязки векового смещения перигелия планеты определяется по разности долгот перигелиев, найденных с помощью решения уравнений движения двумя различными методами. В нашем случае сравниваются долготы перигелиев, найденные на основе решения Ньютоновых уравнений (3) и полученные с помощью каталога БЕ405. Долгота перигелия определяется как сумма значений аргумента перигелия и долготы восходящего узла. Как видно из табл. 7, смещение долготы перигелия Венеры зависит главным образом от скорости движения аргумента перигелия, поскольку долготы восходящих узлов, найденные вышеуказанными методами, совпадают на всем рассматриваемом интервале интегрирования. Наибольшие различия средних аномалий и аргументов перигелиев наблюдаются вблизи концов интервала интегрирования. Например, 9 декабря 1599 г. значение средней аномалии, найденной с помощью решения уравнений (2) и полученной с использованием БЕ405, отличается от средней аномалии на -0.0062°, а аргументы перигелиев— на 0.0063°. На момент 31 декабря 2199 г. расхождение в средней аномалии составляет 0.0037°, в аргументе перигелия -0.0038°. Как видно из табл. 7, величина расхождения средних аномалий, найденная двумя методами, находится в прямой зависимости от вековых скоростей движения аргументов перигелиев. По разностям аргументов перигелиев, полученных путем решения уравнений (2) и Ньютоновых уравнений (3), а также найденных с помощью каталога БЕ405, определена невязка векового смещения перигелия Венеры и ее погрешность для БЕ405. Величина невязки векового смещения для Венеры равна 8.4'' с погрешностью 6.06''.
Как видно из табл. 4, барицентр Земли + Луны имеет малый угол наклона к эклиптике, поэтому изменения аргумента перигелия и долготы восходящего узла в процессе движения более значительны по сравнению с Венерой. В 1999 году наклонение орбиты барицентра Земли + Луны находилось вблизи нуля, принимая в некоторые моменты времени отрицательные значения. В этом случае нисходящий и восходящий узлы меняются местами, что наглядно представлено в табл. 4. Несмотря на значительные различия аргументов перигелиев и долгот восходящих узлов — на -0.0440° и 0.0436° (см. табл. 4), средние аномалии отличаются незначительно: на -0.0003°. На концах интервала интегрирования 9 декабря 1599 г. различие средних аномалий составляет -0.0030°, а в аргументе перигелия и долготе восходящего узла 0.0042° и -0.0009°. По данным табл. 8 найдена невязка векового смещения перигелия барицентра Земли + Луны и ее погрешность: 4.32'' — невязка, 3.83'' — погрешность.
Расхождение средних аномалий и аргументов перигелиев у Марса менее значительное по сравнению с Венерой и Землей (см. табл. 9). Максимальное различие средних аномалий и аргументов перигелиев, полученных путем решения уравнений (2) и с использованием БЕ405, имеет место 9 декабря 1599 г. и составляет -0.0014° и 0.0013° соответственно. Решением Ньютоновых уравнений с использованием БЕ405 найдены вековое смещение перигелия Марса и его погрешность: 1.14'' — смещение, 1.01'' — погрешность.
Таблица 6
Разности элементов орбит Меркурия, полученные с помощью DE405 и решения уравнений (2) и (3) [Differences in the elements of the orbits of the Mercury, obtained using the DE405 and solving the Eqs. (2) and the Eqs. (3)]
Current, date Differences (between the results obtained by) AM (in deg.) Да (in au) Де Дш (in deg.) ДП (in deg.) Дг (in deg.)
1599 12 9 JD 2305424.5 by the Eqs. (2) and the DE405 by the Eqs. (3) and the DE405 -0.0002 0.0139 0 0 0 0 -0.0001 0.0441 0 0.0001 0 0
2199 12 31 JD 2524592.5 by the Eqs. (2) and the DE405 by the Eqs. (3) and the DE405 0.0001 -0.0087 0 0 0 0 0.0001 -0.0276 0 0 0 0
Таблица 7
Разности элементов орбит Венеры, полученные с помощью DE405 и решения уравнений (2) и (3) [Differences in the elements of the orbits of the Venus, obtained using the DE405 and solving the Eqs. (2) and the Eqs. (3)]
Differences (between AM (in deg.) Да (in au) Дш (in deg.) ДП (in deg.) Дг (in deg.)
Current, date the results obtained by) Де
1599 12 9 by the Eqs. (2) and the DE405 -0.0062 0 0 0.0063 0 0
JD 2305424.5 by the Eqs. (3) and the DE405 -0.0055 0 0 0.0087 0 0
2199 12 31 by the Eqs. (2) and the DE405 0.0037 0 0 -0.0038 0 0
JD 2524592.5 by the Eqs. (3) and the DE405 0.0033 0 0 -0.0053 0 0
Таблица 8
Разности элементов орбит Земли + Луны, полученные с помощью DE405 и решения уравнений (2) и (3) [Differences in the elements of the orbits of the Earth & Moon, obtained using the DE405 and solving the Eqs. (2) and the Eqs. (3)]
Current, date Differences (between the results obtained by) AM (in deg.) Да (in au) Де Дш (in deg.) ДП (in deg.) Дг (in deg.)
1599 12 9 JD 2305424.5 by the Eqs. (2) and the DE405 by the Eqs. (3) and the DE405 -0.0030 -0.0028 0 0 0 0 0.0042 0.0048 -0.0009 -0.0009 0 0
2199 12 31 JD 2524592.5 by the Eqs. (2) and the DE405 by the Eqs. (3) and the DE405 0.0019 -0.0004 0 0 0 0 -0.0012 -0.0015 -0.0010 -0.0009 0 0
В табл. 10-14 приведены сравнения элементов орбит Юпитера, Сатурна, Урана, Нептуна и Плутона, полученные с использованием БЕ405 и численного интегрирования дифференциальных уравнений (2) и Ньютоновых уравнений движения (3).
На основании сопоставления элементов орбит можно заключить, что в пределах рассматриваемой точности различий в элементах орбит Юпитера, Сатурна, Урана, Нептуна и Плутона, полученных с использованием БЕ405 и численного интегрирования дифференциальных уравнений (2), не обнаружено. При сравнении элементов орбит полагалось следующее допущение: элементы орбит имеют одинаковую точность, если угловые элементы совпадают с точностью до секунды, а большие полуоси и эксцентриситеты — до шести значащих цифр.
В релятивистских уравнениях результаты вычислений существенным образом зависят от выбранной системы координат. Для уравнений (4) выбрана гармоническая система координат, которая оказалась наиболее подходящей для исследования движения Меркурия. Решение уравнений (4) позволило полностью согласовать величину векового движения перигелия Меркурия с наблюдениями.
Возникает вопрос: для всех ли больших планет при исследовании их движений использование гармонической системы координат является оправданным? Ответ на данный вопрос можно получить либо путем сравнения вычисленных положений планет с наблюдениями, либо путем сравнения координат планет, найденных с помощью решения более эффективным методом. В силу того, что величины вековых смещений перигелиев у Венеры, Земли и Марса малы, а их орбиты близки к круговым, решить данный вопрос путем сравнения с наблюдениями затруднительно. Кроме того, отсутствие на всем исследуемом интервале времени высокоточных наблюдений также препятствует решению данного вопроса.
Проверка эффективности применения уравнений (2) и (4) нами проводилась на примере исследования движения Луны с помощью решения этих уравнений. При исследовании движения Луны в БЕ405 наряду с решением уравнений (4) проводилось интегрирование уравнений (5) и (6). Полученные решения координат Луны согласованы как с оптическими, так и с радиолокационными наблюдениями.
В табл. 15 представлены геоцентрические координаты и компоненты скорости Луны, полученные с использованием БЕ405 и путем решения уравнений (2) и (4).
Как видно из табл. 15, максимальные расхождения координат, полученных путем решения уравнений (2) и (4), с данными, полученными с использованием БЕ405, достигаются на концах интервала интегрирования. При этом расхождения координат и компонент скоростей, найденные путем решения уравнений (2) и с использованием БЕ405, почти на порядок меньше по сравнению с расхождениями решений уравнений (4) и данными с использованием БЕ405.
При сравнении координат Луны, вычисленных с помощью решения уравнений (2) и (4), можно оценить лишь относительные погрешности решений. Сравнение геоцентрических расстояний Луны различными методами позволяет оценить абсолютные погрешности решений данными методами.
Разности элементов орбит Марса, полученные с помощью DE405 и решения уравнений (2) и (3) [Differences in the elements of the orbits of the Mars, obtained using the DE405 and solving the Eqs. (2) and the Eqs. (3)]
Differences (between AM (in deg.) Да (in au) Дш (in deg.) ДП (in deg.) Дг (in deg.)
Current, date the results obtained by) Де
1599 12 9 by the Eqs. (2) and the DE405 -0.0014 0 0 0.0013 0 0
JD 2305424.5 by the Eqs. (3) and the DE405 -0.0014 0 0 0.0014 0 0
2199 12 31 by the Eqs. (2) and the DE405 0.0009 0 0 -0.0008 0 0
JD 2524592.5 by the Eqs. (3) and the DE405 0.0009 0 0 -0.0009 0 0
Таблица 10
Элементы орбит Юпитера [Elements of the orbits of the Jupiter]
Current date Data sources (calculated by) M (in deg.) a (in au) e to (in deg.) П (in deg.) г (in deg.)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 119.7747 119.7746 119.7746 5.206547 5.206547 5.206547 0.047515 0.047515 0.047515 273.6536 273.6536 273.6536 99.7415 99.7415 99.7415 1.3112 1.3112 1.3112
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 330.1812 330.1812 330.1812 5.208676 5.208676 5.208676 0.049767 0.049767 0.049767 272.6198 272.6198 272.6198 100.8539 100.8539 100.8539 1.2993 1.2993 1.2993
Элементы орбит Сатурна [Elements of the orbits of the Saturn]
Current, date Data sources (calculated by) M (in deg.) a (in au) e ш (in deg.) П (in deg.) г (in deg.)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 105.4116 105.4116 105.4116 9.574720 9.574720 9.574720 0.057706 0.057706 0.057706 340.3124 340.3124 340.3124 114.6839 114.6839 114.6839 2.4793 2.4793 2.4793
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 245.9046 245.9046 245.9046 9.588370 9.588370 9.588370 0.053906 0.053906 0.053906 335.9770 335.9770 335.9770 113.1024 113.1024 113.1024 2.4927 2.4927 2.4927
Таблица 12 Элементы орбит Урана [Elements of the orbits of the Uranus]
Current, date Data sources (calculated by) M (in deg.) a (in au) e ш (in deg.) П (in deg.) г (in deg.)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 230.5747 230.5747 230.5747 19.185204 19.185204 19.185204 0.047213 0.047213 0.047213 95.1233 95.1233 95.1233 73.6074 73.6074 73.6074 0.7814 0.7814 0.7814
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 279.3215 279.3215 279.3215 19.183852 19.183851 19.183851 0.045924 0.045924 0.045924 96.7465 96.7465 96.7465 74.0277 74.0277 74.0277 0.7679 0.7679 0.7679
Элементы орбит Нептуна Elements of the orbits of the Neptune]
Current, date Data sources (calculated by) M (in deg.) a (in au) e ш (in deg.) П (in deg.) г (in deg.)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 61.7486 61.7487 61.7487 30.270547 30.270547 30.270547 0.009574 0.009574 0.009574 317.2798 317.2798 317.2798 131.7883 131.7883 131.7883 1.7671 1.7671 1.7671
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 353.1376 353.1376 353.1376 30.264791 30.264791 30.264791 0.014321 0.014321 0.014321 256.6947 256.6947 256.6947 131.7178 131.7178 131.7178 1.7702 1.7702 1.7702
Таблица 14
Элементы орбит Плутона [Elements of the orbits of the Pluto]
Current, date Data sources (calculated by) M (in deg.) a (in au) e ш (in deg.) П (in deg.) г (in deg.)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 154.6398 154.6399 154.6399 39.397879 39.397879 39.397880 0.249562 0.249562 0.249562 113.0311 113.0311 113.0311 110.5004 110.5004 110.5004 17.1365 17.1365 17.1365
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (3) 304.1102 304.1102 304.1102 39.257516 39.257516 39.257516 0.247022 0.247022 0.247022 114.7477 114.7477 114.7477 110.3397 110.3397 110.3397 17.1680 17.1680 17.1680
Геоцентрические координаты и компоненты скоростей Луны, полученные по ОЕ405 и путем решения уравнений (2) и (4)
[Geocentric coordinates and velocity components of the Moon calculated by the DE405, the Eqs. (2) and the Eqs. (4)]
Current date Data sources (calculated by) X (in au) Y (in au) Z (in au) Vx (in au/day) Vy (in au/day) V, (in au/day)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.002596600 -0.002596680 -0.002600840 0.000681096 0.000681446 0.000667855 0.000131635 0.000129649 0.000122528 -0.000155519 -0.000155389 -0.000152079 -0.000477071 -0.000477231 -0.000478121 -0.000251468 -0.000251209 -0.000251385
1699 12 28 JD 2341968.5 by the DE405 by the Eqs. (2) by the Eqs. (4) 0.002643269 0.002643220 0.002644275 -0.000298812 -0.000298181 -0.000287487 -0.000140145 -0.000141896 -0.000138329 0.0000430142 0.0000430006 0.0000405114 0.000536555 0.000536541 0.000536757 0.000178574 0.000178638 0.000178749
1800 1 17 JD 2378512.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.002440199 -0.00244029 -0.00243880 -0.000525672 -0.000524984 -0.000530840 -0.000134793 -0.000135636 -0.000138776 0.000142850 0.000142793 0.000144268 -0.000517844 -0.000517814 -0.000517424 -0.000278656 -0.000278747 -0.000278636
1899 12 4 JD 2414992.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.000198405 -0.000198560 -0.000195631 -0.002288453 -0.002288571 -0.002288897 -0.000976145 -0.000975855 -0.000975701 0.000607334 0.000607335 0.000607362 -0.000047648 -0.000047659 -0.000047016 -0.000040291 -0.000040248 -0.000040520
1999 12 24 JD 2451536.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.000808901 -0.000808941 -0.000807684 0.002094014 0.002093992 0.002094421 0.000849050 0.000848989 0.000849049 -0.000600679 -0.000600683 -0.000600792 -0.000198425 -0.000198425 -0.000198136 -0.000025796 -0.000025845 -0.000025728
2099 12 11 JD 2488048.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.001038446 -0.001037717 -0.001042765 -0.001965349 -0.001965712 -0.001963643 -0.001077336 -0.001077272 -0.001076232 0.000551769 0.000551867 0.000551261 -0.000245541 -0.000245301 -0.000246343 -0.000124050 -0.000124145 -0.000124718
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (4) 0.000388652 0.000386760 0.000396300 0.002506014 0.002505890 0.002504366 0.000898766 0.000899905 0.000899895 -0.000557717 -0.000557775 -0.000557422 0.000099983 0.000099700 0.000101601 0.000004597 0.000004257 0.000004940
Таблица 16
Геоцентрические расстояния Луны, полученные по ОЕ405 и путем решения уравнений (2) и (4)
[Geocentric distance of the Moon calculated by the DE405 and by solving the Eqs. (2) and the Eqs. (4)]
Current date Data sources (calculated by) X (in au) Y (in au) Z (in au) p (in au) \Ap\ (in km)
1599 12 9 JD 2305424.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.002596600 -0.002596680 -0.002600840 0.000681096 0.000681446 0.000667855 0.000131635 0.000129649 0.000122528 0.002687666 0.002687736 0.002688013 10.47 51.91
1699 12 28 JD 2341968.5 by the DE405 by the Eqs. (2) by the Eqs. (4) 0.002643269 0.002643220 0.002644275 -0.000298812 -0.000298181 -0.000287487 -0.000140145 -0.000141896 -0.000138329 0.002663794 0.002663768 0.002663447 3.89 51.91
1800 1 17 JD 2378512.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.002440199 -0.00244029 -0.00243880 -0.000525672 -0.000524984 -0.000530840 -0.000134793 -0.000135636 -0.000138776 0.002499815 0.002499804 0.002499760 1.64 8.23
1899 12 4 JD 2414992.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.000198405 -0.000198560 -0.000195631 -0.002288453 -0.002288571 -0.002288897 -0.000976145 -0.000975855 -0.000975701 0.002495845 0.002495852 0.002495886 1.04 6.13
1999 12 24 JD 2451536.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.000808901 -0.000808941 -0.000807684 0.002094014 0.002093992 0.002094421 0.000849050 0.000848989 0.000849049 0.002400021 0.002399999 0.002399965 3.29 8.38
2099 12 11 JD 2488048.5 by the DE405 by the Eqs. (2) by the Eqs. (4) -0.001038446 -0.001037717 -0.001042765 -0.001965349 -0.001965712 -0.001963643 -0.001077336 -0.001077272 -0.001076232 0.002470146 0.002470100 0.002470125 6.88 3.77
2199 12 31 JD 2524592.5 by the DE405 by the Eqs. (2) by the Eqs. (4) 0.000388652 0.000386760 0.000396300 0.002506014 0.002505890 0.002504366 0.000898766 0.000899905 0.000899895 0.002690527 0.002690519 0.002690486 1.20 6.13
В табл. 16 представлены разности геоцентрических расстояний Луны в километрах, полученные с помощью решений уравнений (2), (4) и DE405.
В первой строке табл. 16 представлены координаты и геоцентрическое расстояние Луны, найденное с помощью DE405, во второй и третьей строках — полученные с помощью решения уравнений (2) и (4). В последней колонке табл. 16 дается разность между геоцентрическими расстояниями Луны, найденными путем решения уравнений (2) и (4) и определенными с помощью планетного каталога DE405. Полагая, что координаты Луны, найденные с помощью DE405, более точные по сравнению с координатами, полученными с помощью решения уравнений (2) и (4), оценим разности геоцентрических расстояний между точным и приближенным решением. Тот метод, у которого данная разность меньше, будем считать более точным. Данное условие является необходимым, но не достаточным для оценки точности метода. Из сопоставления разностей геоцентрических расстояний Луны, найденных различными методами, следует, что максимальные расхождения разностей, полученных с помощью решения уравнений (2) и DE405, изменяются в пределах от 1.20 до 10.47 км, в то время как разности, найденные с помощью решения уравнений (4), изменяются в пределах от 3.77 до 51.91 км.
Как видно из табл. 16, решения, полученные с помощью уравнений (2), являются более точными в шести случаях из семи по сравнению с решениями уравнений (4). Следовательно, элементы орбит планет, найденные по координатам и скоростям, полученным путем решения уравнений (2), более точные по сравнению с элементами орбит, полученных на основании решения уравнений (4). Поскольку при создании DE405 движение больших планет рассчитывалось с помощью решения уравнений (4), путем сравнения элементов орбит планет, найденных с использованием решений уравнений (2) и (4), можно оценить точность полученных решений, представленных в каталоге DE405 в форме коэффициентов многочленов Чебышева.
На основании сравнения элементов орбит Венеры, барицентра Земли + Луны и Марса, полученных с помощью решения уравнений (2) и (4), найденные избыточные смещения вековых долгот перигелиев с использованием каталога DE405 составляют: 6.06'', 3.83'' и 1.01''.
В заключение отметим основные результаты проведенных исследований:
1) использование гармонической системы координат в релятивистских уравнениях движения оправдано для Меркурия и для внешних планет: Юпитера, Сатурна, Урана, Нептуна и карликовой планеты Плутон;
2) показано, что использование релятивистских уравнений при создании DE405 для исследования движения планет Венеры, барицентра Земли + Луны и Марса, приводит к избыточному вековому смещению долгот перигелиев этих планет;
3) величины погрешностей координат и компонент скоростей Венеры, барицентра Земли + Луны и Марса, найденных с использованием DE405, находятся в прямой зависимости от погрешностей величин вековых смещений долгот перигелиев этих планет;
4) решение, полученное с помощью Ньютоновых уравнений для планет Венеры, барицентра Земли + Луны и Марса, не уступает по точности решению, найденному путем решения уравнений (4).
Сравнение элементов орбит больших планет, Луны и Солнца.
Конкурирующие интересы. Заявляем, что в отношении авторства и публикации этой статьи конфликта интересов не имеем.
Авторский вклад и ответственность. Все авторы принимали участие в разработке концепции статьи и в написании рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.
Финансирование. Исследование выполнялось без финансирования.
Библиографический список
1. Максвелл Дж. К. Избранные сочинения по теории электромагнитного поля/ Классики естествознания. Т. 104. М.: Гостехиздат, 1952. 687 с.
2. Чеботарев Г. А. Аналитические и численные методы небесной механики. М., Л.: Наука, 1965. 368 с.; Chebotarev G. A. Analytical and Numerical Methods of Celestial Mechanics / Modern Analytic and Computational Methods in Science and Mathematics. vol. 9. New York: American Elsevier Publishing Co., Inc., 1967. xviii+331 pp.
3. Субботин М. Ф. Введение в теоретическую астрономию. М.: Наука, 1968. 800 с.
4. Богородский А. Ф. Всемирное тяготение. Киев: Наук. думка, 1971. 352 с.
5. Брумберг В. А. Релятивистская небесная механика. М.: Наука, 1972. 384 с.
6. Визгин В. П. Об открытии уравнений гравитационного поля Эйнштейном и Гильбертом (новые материалы)// УФН, 2001. Т. 171, №12. С. 1347-1363. EDN: MPJJEZ. DOI: https://doi.org/10.3367/UFNr.0171.200112d.1347.
7. Дирак П. Принципы квантовой механики. М.: Наука, 1979. 440 с.
8. Фейнман Р. Развитие пространственно-временной трактовки квантовой электродинамики// УФН, 1967. Т. 91, №1. С. 29-48. DOI: https://doi.org/10.3367/UFNr.0091. 196701c.0029.
9. Casimir H. B. G., Polder D. The influence of retardation on the London-van der Waals forces// Phys. Rev., 1948. vol.73, no. 4. pp. 360-372. DOI: https://doi.org/10.1103/ PhysRev.73.360.
10. Зельдович Я. Б. Теория вакуума, быть может, решает загадку космологии // УФН, 1981. Т. 133, №3. С. 479-503. DOI: https://doi.org/10.3367/UFNr.0133.198103c.0479.
11. Заусаев А. Ф. Теория движения п материальных тел, основанная на новом принципе взаимодействия// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2006. №43. С. 132-139. EDN: HUYDOB. DOI: https://doi.org/10.14498/vsgtu463. [Zausaev A. F. Theory of motion of n material bodies, based on a new interaction principle / Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2006. pp. 132-139].
12. Заусаев А. Ф., Заусаев А. А. Математическое моделирование орбитальной эволюции малых тел Солнечной системы. М.: Машиностроение-1, 2008. 250 с. EDN: QJUWVV.
13. Заусаев А. Ф. Исследование движения планет, Луны и Солнца, основанное на новом принципе взаимодействия // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2014. №3(36). С. 118-131. EDN: TLEURX. DOI: https://doi.org/10.14498/vsgtu1304.
14. Заусаев А. Ф. Сопоставление координат больших планет, Луны и Солнца, полученных на основе нового принципа взаимодействия и банка данных DE405 // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2016. Т. 20, №1. С. 121-148. EDN: WQPYBB. DOI: https://doi.org/10.14498/vsgtu1458.
15. Заусаев А. Ф., Романюк М. А. Численные методы в задачах математического моделирования движения небесных тел в Солнечной системе. Самара: СамГТУ, 2017. 265 с.
16. Newhall X. X., Standish E M., Williams J. G. DE 102: A numerically integrated ephemeris of the Moon and planets spanning forty-four centuries // Astron. Astrophys, 1983. vol. 125, no. 1. pp. 150-167.
17. Пуанкаре А. О науке. М.: Наука, 1983. 560 с.
18. Заусаев А. Ф., Романюк М. А. Сравнение различных математических моделей на примере решения уравнений движения больших планет и Луны // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2019. Т. 23, №1. С. 152-185. EDN: BGRYUJ. DOI:https:// doi.org/10.14498/vsgtu1663.
19. Standish E. M. JPL Planetary and Lunar Ephemerides: DE405/LE405. Interoffice memorandum: JPL IOM 312. F-98-048, 1998, August 26. 18 pp. ftp://ssd.jpl.nasa.gov/pub/ eph/planets/ioms/de405.iom.pdf.
20. Грибкова В. П. Эффективные методы равномерных приближений, основанные на полиномах Чебышева. М.: Спутник, 2017. 194 с.
21. Montenbruck O., Pfleger T. Astronomy on the Personal Computer. Berlin, Heidelberg: Springer, 2000. xv+300 pp. DOI: https://doi.org/10.1007/978-3-642-03436-7.
22. Хемминг Р. В. Численные методы для научных работников и инженеров. М.: Наука, 1972. 400 с.
23. Park R. S. JPL Planetary and Lunar Ephemerides, 2020. https://ssd.jpl.nasa.gov/ planets/eph_export.html.
24. Le Verrier U. J. Theorie du movement de Mercure / Annales de l'Observatoire imperial de Paris. vol. 5: Annales de l'Observatoire de Paris. Memoires. Paris: Mallet-Bachelier, 1859. 195 pp.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2022, vol. 26, no. 4, pp. 738-763 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1888
MSC: 70F15, 70M20, 65L05
Comparison of the orbital elements of major planets, the Moon and the Sun using various mathematical models on the time interval with 1600 to 2200
A. F. Zausaev, M. A. Romanyuk
Samara State Technical University,
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation.
Abstract
An analysis of the accuracy of the orbital elements obtained according to the coordinates and components of the velocities, found using the coefficients of the Chebyshev polynomials of the DE405 planetary catalog, is carried out. We compared the elements of orbital elements in the time interval from 1600 to 2200 years found using the DE405 catalog and obtained by numerical integration of the equations of motion based on the interaction of moving material bodies with the surrounding space. On the example of the numerical integration of the Moon motion equations, the advantage of using the equations of motion based on the interaction of moving material bodies with the surrounding space is shown in comparison with relativis-tic equations. Based on a comparison of the elements of Mercury's orbits, found by coordinates obtained by solving equations based on the interaction of moving material bodies with the surrounding space, and obtained using the DE405 catalog, it is shown that the orbital elements practically coincide on a given interval time. The maximum discrepancy in the mean anomaly at the end of the integration interval is less than 1'' (second). The discrepancies of the secular displacements of perihelions for Mercury, Venus, Earth + Moon and Mars were determined, the values of which for DE405 are respectively: 43.08'', 8.4'', 3.83'' and 1.14''. It is shown that the errors of the secular displacements of the perihelions of the planets Mercury, Venus, the barycenter of the Earth + Moon and Mars obtained using the DE405 catalog take the following values: 0'', 6.06'', 3.83'' and 1.08''. For the outer
Mathematical Modeling, Numerical Methods and Software Complexes Research Article
© Authors, 2022
© Samara State Technical University, 2022 (Compilation, Design, and Layout) Q ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Zausaev A. F., Romanyuk M. A. Comparison of the orbital elements of major planets, the Moon and the Sun using various mathematical models on the time interval with 1600 to 2200, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2022, vol. 26, no. 4, pp. 738-763. EDN: DVJGRH. DOI: 10.14498/vsgtu1888 (In Russian).
Authors' Details:
Anatoliy F. Zausaev © https://orcid.org/0000-0002-5035-9615
Dr. Phys. & Math. Sci.; Professor; Dept. of Applied Mathematics & Computer Science;
e-mail: zausaev_af@mail.ru
Mariya A. Romanyuk A https://orcid.org/0000-0003-1352-6725
Cand. Tech. Sci.; Associate Professor; Dept. of Applied Mathematics & Computer Science;
e-mail: zausmasha@mail.ru
planets: Jupiter, Saturn, Uranus, Neptune and the dwarf planet Pluto, on the basis of the considered comparisons of various equations of motion, no discrepancies in the orbital elements were found. Based on the studies carried out, it is shown that the use of harmonic coordinates in relativistic equations when creating the DE405 catalog is justified only for Mercury and the outer planets: Jupiter, Saturn, Uranus, Neptune and the dwarf planet Pluto.
Keywords: orbital elements, numerical integration, differential equation of motion.
Received: 19th October, 2021 / Revised: 15th October, 2022 / Accepted: 17th November, 2022 / First online: 6th December, 2022
Competing interests. We declare that we have no conflicts of interest with respect to the authorship and publication of this article.
Authors' contributions and responsibilities. Each author has participated in the development of the concept of the article and in the writing of the manuscript. The authors are absolutely responsible for submitting the final manuscript in print. Each author has approved the final version of the manuscript. Funding. The research has not received funding.
References
1. Maxwell J. C. Izbrannye sochineniia po teorii elektromagnitnogo polia [Selected Papers on Electromagnetic Field Theory]. Moscow, Gostekhizdat, 1952, 687 pp. (In Russian)
2. Chebotarev G. A. Analytical and Numerical Methods of Celestial Mechanics, Modern Analytic and Computational Methods in Science and Mathematics, vol. 9. New York, American Elsevier Publishing Co., Inc., 1967, xviii+331 pp.
3. Subbotin M. F. Vvedenie v teoreticheskuiu astronomiiu [Introduction to Theoretical Astronomy]. Moscow, Nauka, 1968, 800 pp. (In Russian)
4. Bogorodsky A. F. Vsemirnoe tiagotenie [Universal Gravitation]. Kiev, Nauk. Dumka, 1971, 352 pp. (In Russian)
5. Brumberg V. A. Reliativistskaia nebesnaia mekhanika [Relativistic Celestial Mechanics]. Moscow, Nauka, 1972, 384 pp. (In Russian)
6. Vizgin V. P. On the discovery of the gravitational field equations by Einstein and Hilbert: new materials, Phys. Usp., 2001, vol.44, no. 12, pp. 1283-1298. EDN: LGQDQR. DOI: https:// doi.org/10.1070/PU2001v044n12ABEH001014.
7. Dirak P. Printsipy kvantovoi mekhaniki [The Principles of Quantum Mechanics]. Moscow, Nauka, 1979, 440 pp. (In Russian)
8. Feynman R. P. The Development of the Space-Time View of Quantum Electrodynamics, Nobel Lecture, December 11, 1965. Preprint les Prix Nobel en 1965. Stockholm, The Nobel Foundation, 1966.
9. Casimir H. B. G., Polder D. The influence of retardation on the London-van der Waals forces, Phys. Rev., 1948, vol.73, no. 4, pp. 360-372. DOI: https://doi.org/10.1103/PhysRev.73. 360.
10. Zel'dovich Ya. B. Vacuum theory: a possible solution to the singularity problem of cosmology, Sov. Phys. Usp., 1981, vol.24, no. 3, pp. 216-230. DOI: https://doi.org/10.1070/ PU1981v024n03ABEH004772.
11. Zausaev A. F. Theory of motion of n material bodies, based on a new interaction principle, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2006, no. 43, pp. 132-139. EDN: HUYDOB. DOI: https://doi.org/ 10.14498/vsgtu463.
Comparison of the orbital elements of major planets, the Moon and the Sun...
12. Zausaev A. F., Zausaev A. A. Matematicheskoe modelirovanie orbital'noi evoliutsii malykh tel Solnechnoi sistemy [Mathematical Modelling of Orbital Evolution of Small Bodies of the Solar System]. Moscow, Mashinostroenie-1, 2008, 250 pp. (In Russian). EDN: QJUWVV.
13. Zausaev A. F. The investigation of the motion of planets, the Moon, and the Sun based on a new principle of interaction, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2014, no. 3(36), pp. 118-131 (In Russian). EDN: TLEURX. DOI: https://doi.org/10.14498/vsgtu1304.
14. Zausaev A. F. Comparison of the coordinates of the major planets, Moon, and Sun obtained based on a new principle of interaction and of the data bank DE405, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2016, vol. 20, no. 1, pp. 121-148. EDN: WQPYBB. DOI: https://doi.org/10.14498/vsgtu1458.
15. Zausaev A. F., Romanyuk M. A. Chislennye metody v zadachakh matematicheskogo mod-elirovaniia dvizheniia nebesnykh tel v Solnechnoi sisteme [Numerical Methods in the Problems of Mathematical Modeling of the Motion of Celestial Bodies in the Solar System]. Samara, Samara State Technical Univ., 2017, 265 pp. (In Russian)
16. Newhall X. X., Standish E M., Williams J. G. DE 102: A numerically integrated ephemeris of the Moon and planets spanning forty-four centuries, Astron. Astrophys., 1983, vol. 125, no. 1, pp. 150-167.
17. Poincaré J. H. O nauke [About Science]. Moscow, Nauka, 1983, 560 pp. (In Russian)
18. Zausaev A. F., Romanyuk M. A. Comparison of various mathematical models on the example of solving the equations of the movement of large planets and the Moon, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2019, vol.23, no. 1, pp. 152-185. EDN: BGRYUJ. DOI: https://doi.org/ 10.14498/vsgtu1663.
19. Standish E. M. JPL Planetary and Lunar Ephemerides, DE405/LE405. Interoffice memorandum: JPL IOM 312. F-98-048, 1998, August 26, 18 pp. ftp://ssd.jpl.nasa.gov/pub/ eph/planets/ioms/de405.iom.pdf.
20. Gribkova V. P. Effektivnye metody ravnomernykh priblizhenii, osnovannye na polinomakh Chebysheva [Efficient Methods of Uniform Approximations Based on Chebyshev Polynomials]. Moscow, Sputnik, 2017, 194 pp. (In Russian)
21. Montenbruck O., Pfleger T. Astronomy on the Personal Computer. Berlin, Heidelberg, Springer, 2000, xv+300 pp. DOI: https://doi.org/10.1007/978-3-642-03436-7.
22. Khemming R. V. Chislennye metody dlia nauchnykh rabotnikov i inzhenerov [Numerical Methods for Scientists and Engineers]. Moscow, Nauka, 1972, 400 pp. (In Russian)
23. Park R. S. JPL Planetary and Lunar Ephemerides, 2020. https://ssd.jpl.nasa.gov/ planets/eph_export .html.
24. Le Verrier U. J. Theorie du movement de Mercure, Annales de l'Observatoire imperial de Paris, vol. 5, Annales de l'Observatoire de Paris. Memoires. Paris, Mallet-Bachelier, 1859, 195 pp.