Эфемериды DE-405.

[Toth]

Если посмотреть на равноденствия

Равноде́нствие — астрономическое явление, когда центр Солнца в своём видимом движении по эклиптике пересекает небесный экватор
Отсюда - https://ru.wikipedia.org/wiki/Равноденствие
То есть, когда дельта = 0
Из IMCCE -

[filiasfox]

Равноде́нствие — астрономическое явление, когда центр Солнца в своём видимом движении по эклиптике пересекает небесный экватор
Отсюда - https://ru.wikipedia.org/wiki/Равноденствие
То есть, когда дельта = 0
Из IMCCE -

Увидел. Но вопрос относительно координаты Z остался.

[Toth]

А тип координат какой?
Я там подчеркнул , что ApparentOfDate, у вас наверное MeanJ2000.

[filiasfox]

Я там подчеркнул , что ApparentOfDate, у вас наверное MeanJ2000.

Поправил и увидел. Но что делать с Z, хотя и просто деление координат Y на косинус наклона оси потом даст растущую погрешность, тоже пока не ясно

[Toth]

Но вопрос относительно координаты Z остался.

А что Z - у вас по нему расхождение с IMCCE-405  всего 64 метра. Вот у Y - 490 м, но и это не так уж много.
Я не знаю , какая разница между 405 и 423, может быть примерно такая и есть.

[filiasfox]

А что Z - у вас по нему расхождение с IMCCE-405  всего 64 метра. Вот у Y - 490 м, но и это не так уж много.
Я не знаю , какая разница между 405 и 423, может быть примерно такая и есть.

Не 423. 422.
Но можно взять несколько дат для проверки.

[Toth]

Вот, NASA , правда там сейчас DE431 -

[filiasfox]

Вот, NASA , правда там сейчас DE431 -

Сверил. Похоже не на те данные смотрел, когда сверял. Нестыковка не более 200 метров по Y и Z.
Тогда возникает вопрос в отношении следующего. Опять смотрим 2019 год

2458563.583333333 = A.D. 2019-Mar-21 02:00:00.0000 TDB
 X = 1.489872660541371E+08 Y =-2.207569373830897E+05 Z =-9.577695313739940E+04

2458563.625000000 = A.D. 2019-Mar-21 03:00:00.0000 TDB
 X = 1.489891469888312E+08 Y =-1.220389968882262E+05 Z =-5.298085753371072E+04

2458563.666666667 = A.D. 2019-Mar-21 04:00:00.0000 TDB
 X = 1.489909509257654E+08 Y =-2.332097191772831E+04 Z =-1.018476824548427E+04

2458563.708333333 = A.D. 2019-Mar-21 05:00:00.0000 TDB
 X = 1.489926778663567E+08 Y = 7.539709111577214E+04 Z = 3.261129442939936E+04

2458563.750000000 = A.D. 2019-Mar-21 06:00:00.0000 TDB
 X = 1.489943278119982E+08 Y = 1.741151457915737E+05 Z = 7.540731019577320E+04

То есть как увязать время наступления равноденствия согласно этим насовским данным и тому, что есть?

[Toth]

То есть как увязать время наступления равноденствия согласно этим насовским данным и тому, что есть?

Так же, как я и говорил. MeanJ2000-->ApparentOfDate
То есть учесть прецессию, нутацию, аберрацию, световременную коррекцию.

[filiasfox]

Так же, как я и говорил. MeanJ2000-->ApparentOfDate
То есть учесть прецессию, нутацию, аберрацию, световременную коррекцию.

Честно говоря не уверен, что в сумме эти показатели дадут большую погрешность, что выходит например на промежутке 2010-2024 для весенних равноденствий. Пример ниже. Последний столбец - то, что было рассчитано (переход Z из минуса в плюс). Horizons (НАСА) дает значение, еще больше отходящее от требуемого. Может подскажете, что здесь может быть не так?

2010 20 17:32:13   14:58
2011 20 23:21:44   21:08
2012 20 05:14:25   03:21
2013 20 11:02:55   09:28
2014 20 16:57:05   15:41
2015 20 22:45:09   21:47
2016 20 04:30:11   03:51
2017 20 10:28:38   10:07
2018 20 16:15:27   16:13
2019 20 21:58:25   22:14
2020 20 03:50:36   04:25
2021 20 09:37:27   10:33
2022 20 15:33:23   16:51
2023 20 21:24:24   23:04
2024 20 03:06:21   05:08

[filiasfox]

То есть учесть прецессию, нутацию, аберрацию, световременную коррекцию.

Добавил расчет величин, которые вы указали. Данные если и поменялись, то не намного. Влияние прецессии, нутации, аберрации и световременной коррекции не дало нужного результата. Собственно результат мало отличается от того, что показал выше.

Тестовые данные и собственно данные, полученные из расчета (последний столбец - то, что получено из расчета). Примерно в 2019 году нестыковка почти исчезает. Потом снова растет в другую сторону. Чтобы это могло быть

20/03/2000 07h 35m 16s   07:32:09
20/03/2001 13h 30m 44s   13:48:25
20/03/2002 19h 16m 10s   19:53:43
21/03/2003 00h 59m 47s   01:57:16
20/03/2004 06h 48m 39s   08:07:04
20/03/2005 12h 33m 26s   14:11:09
20/03/2006 18h 25m 35s   20:24:50
21/03/2007 00h 07m 26s   02:25:50
20/03/2008 05h 48m 19s   08:27:42
20/03/2009 11h 43m 39s   14:43:22

[Toth]

Попробуйте вместо вектора положения выбрать наблюдения
PS Target Body - конечно Sun, центр. И т.д. , то что вам надо.

[filiasfox]

Попробуйте вместо вектора положения выбрать наблюдения
PS Target Body - конечно Sun, центр. И т.д. , то что вам надо.

Эффект тот же, что и писал выше. Равноденствие 2015 года: 20 марта 22:45:09
Результат ниже c Horizons (NASA).

$$SOE
 2015-Mar-20 00:00     23 55 46.70 -00 27 27.5
 2015-Mar-20 04:00     23 56 23.20 -00 23 30.2
 2015-Mar-20 08:00     23 56 59.70 -00 19 33.0
 2015-Mar-20 12:00     23 57 36.19 -00 15 35.7
 2015-Mar-20 16:00     23 58 12.67 -00 11 38.5
 2015-Mar-20 20:00     23 58 49.16 -00 07 41.3

 2015-Mar-21 00:00     23 59 25.64 -00 03 44.1
 2015-Mar-21 04:00     00 00 02.12 +00 00 13.1

 2015-Mar-21 08:00     00 00 38.59 +00 04 10.2
 2015-Mar-21 12:00     00 01 15.06 +00 08 07.3
 2015-Mar-21 16:00     00 01 51.52 +00 12 04.4
 2015-Mar-21 20:00     00 02 27.99 +00 16 01.5
 2015-Mar-22 00:00     00 03 04.45 +00 19 58.5
$$EOE

Результат с Miriade
Sun    2015-03-20T21:00:00.00    23 59 44.01453    -00 01 43.9107
Sun    2015-03-20T22:00:00.00    23 59 53.13484    -00 00 44.6172
Sun    2015-03-20T23:00:00.00    00 00 2.25494    +00 00 14.6749

ApparentOfDate, DE406

[filiasfox]

То есть учесть прецессию, нутацию, аберрацию, световременную коррекцию.

Спасибо. Вычисления стали лучше сходиться. Правда пришлось лезть в фортран. Но собственно световременная коррекция и корректный способ вычисления прецесии уменьшили расхождение до пяти минут. Буду смотреть дальше, что и как лучше пользовать

[Desert71]

Очевидно, что первым делом нужно вычислить номер нужного блока:
Пусть
MD0 = MJD(1999,12,24) – начальная юлианская дата файла ascp2000.405
MD – нужная дата, на которую необходимо вычислить эфемериды.
Interval = 32 (длина интервала в сутках, во время которого действителен отдельный блок коэффициентов), константа
Номер блока:
NumerBlock=  Int((MD - md0) / Interval) + 1
Начальная дата нужного блока:
MD1 = MD0 + (NumerBlock - 1) * Interval


Затем возникает необходимость вычислить длину и номер подинтервала:
subInterval = Interval / Nset_Body     длина подинтервала в сутках
где Nset_Body – количество подинтервалов для конкретного тела (третья строка группы 1050 в header.405)
Количество полных подинтервалов в блоке до нужной даты:
podint = Int((MD - MD1) / subInterval) 
Начальная дата нужного подинтервала:
Mdat = MD1 + podint * subInterval
Сдвиг начала чтения коэффициентов с учетом подинтервала:
IndBegin = 3 * Nkoef_Body * podint
Где Nkoef_Body – количество коэффициентов для данного тела (вторая строка группы 1050 в header.405).

Заморочно, да. Можно записать всё это покомпактнее, но так более понятно.

Для примера возьмем Землю и захотим вычислить ее барицентрические координаты на 03.01. 2016  5:30 UT.
MD0 =  51536
MD =  57390.2299558333 (это заданная дата + 32.184 сек + TAI_UTC )
NumerBlock  =  183
MD1 =  57360  (183-ий блок начинается с этой даты)
subinterval  =   16 (для Земли Nset_Body = 2)
podint  =  1 (столько подинтервалов нужно пропустить)
Mdat  =   57376 (нужны коэффициенты 2-го подинтервала, который начинается с этой даты)
Из 183-го блока читаем коэффициенты с 231 по 308 (положение коэффициентов для Земли, первая строка группы 1050 в header.405), всего их 78, они принадлежат двум подинтервалам по 16 дней каждый. Нам нужны коэффициенты второго подинтервала.
IndBegin  =    39
Первые 39 коэффициентов пропускаем, и читаем оставшиеся:  kX(от 1 до 13) для X, kY(от 1 до 13) для Y, kZ(от 1 до 13) для Z.

Получили три комплекта чисел по 13 штук в комплекте. Теперь дело за полиномами Чебышева.

Дай Бог тебе здоровья, жениха богатого... ;)ты сделал космическое описание, во всех смыслах. Ты меня просто спас!!!!