Небесно-механические этюды

[cross-track]

А эпохи ТЛЕ какие?

[Олег]

Что интересно

А эпохи ТЛЕ какие?

Что интересно - обе после коррекции.
28.08.2024 7:10
28.08.2024 10:25
Но между ними разница высот 3.45 км. То есть первый - до, второй после коррекции.
Как я понял, время эпохи - не то же что время формирования TLE (определения орбит ) . Например, в Целестраке сейчас TLE XMM-NEWTON с эпохой 06.09.2024  4:39 UTC.

[cross-track]

Что интересно

А эпохи ТЛЕ какие?

Что интересно - обе после коррекции.
28.08.2024 7:10
28.08.2024 10:25
Но между ними разница высот 3.45 км. То есть первый - до, второй после коррекции.
Как я понял, время эпохи - не то же что время формирования TLE (определения орбит ) . Например, в Целестраке сейчас TLE XMM-NEWTON с эпохой 06.09.2024  4:39 UTC.

Да, время эпохи - не то же что время фактического определения орбит. Но в данном случае разница между эпохой и временем формирования измеренной орбиты слишком велика (не менее 5.5 часов). Я такое встречаю впервые. 

Вполне возможно, что для МКС такое делается сознательно - публикуются заведомо не очень точные ТЛЕ с целью предотвращения каких-то зловредных замыслов.

Если исключить такую конспирологию, то можно допустить, что ТЛЕ хорошие, но из-за неспокойного Солнца (и сопутствующего геомагнитного поля) прогноз орбит даже на несколькочасовых промежутках может давать ощутимую неточность. Хотя в это тоже слабо верится.
И еще. Сама коррекция не точечная, она длилась около 18 минут. Так что разница между временем коррекции и вашим расчетным временем немного меньше.

Я попробую тоже найти точное время минимального сближения, если что-то будет интересное, напишу.

[Олег]

Если исключить такую конспирологию, то можно допустить, что ТЛЕ хорошие, но из-за неспокойного Солнца (и сопутствующего геомагнитного поля) прогноз орбит даже на несколькочасовых промежутках может давать ощутимую неточность.

Последняя буря была 28 августа, слабая. Хотя, наверное не только бури влияют.

Я посчитал ошибки от эпохи к эпохе по соседним TLE.
Берем TLE1 и TLE2 , считаем координаты на момент эпохи TLE2 , сравниваем. Потом  TLE2 и TLE3 , и т.д. 
Результат -

Спойлер

Ep	Ep-Ep, day	err,km
28.08.2024 0:53
28.08.2024 7:10	0.262	4.14
28.08.2024 10:25	0.136	181.81
28.08.2024 16:07	0.237	17.49
28.08.2024 16:42	0.024	1.06
29.08.2024 4:30	0.492	14.55
29.08.2024 9:29	0.207	0.69
29.08.2024 12:52	0.141	0.61
29.08.2024 15:20	0.103	0.70
29.08.2024 15:45	0.017	0.39
29.08.2024 17:15	0.062	2.08
30.08.2024 3:44	0.437	15.24
30.08.2024 10:24	0.278	0.45
30.08.2024 11:28	0.045	0.75
30.08.2024 12:05	0.026	0.19
30.08.2024 16:07	0.168	0.95
30.08.2024 16:46	0.027	0.25
31.08.2024 4:30	0.489	1.52
31.08.2024 12:51	0.348	1.89
31.08.2024 15:20	0.103	2.21
31.08.2024 15:47	0.019	0.79
01.09.2024 0:56	0.381	4.55
01.09.2024 4:08	0.134	0.16
01.09.2024 8:49	0.195	0.40
01.09.2024 11:28	0.110	0.63
01.09.2024 12:04	0.025	0.25
01.09.2024 14:33	0.104	1.11
01.09.2024 15:08	0.024	0.78
02.09.2024 1:45	0.442	3.06
02.09.2024 4:57	0.133	2.38
02.09.2024 9:32	0.191	0.29
02.09.2024 12:14	0.113	0.54
02.09.2024 12:53	0.027	0.36
02.09.2024 15:57	0.128	0.09
03.09.2024 1:01	0.378	0.37
03.09.2024 3:43	0.112	0.48
03.09.2024 5:45	0.085	0.37
03.09.2024 7:08	0.058	0.11
03.09.2024 11:27	0.180	0.42
03.09.2024 13:41	0.093	0.47
03.09.2024 22:17	0.358	0.36
03.09.2024 22:32	0.010	0.40
04.09.2024 0:10	0.068	2.21
04.09.2024 4:57	0.200	1.82
04.09.2024 9:41	0.197	0.30
04.09.2024 12:53	0.133	0.16
04.09.2024 13:46	0.037	0.45
04.09.2024 14:22	0.025	0.66
05.09.2024 3:42	0.556	8.89
05.09.2024 5:46	0.086	0.29
05.09.2024 10:18	0.189	0.81
05.09.2024 11:55	0.068	0.10

[свернуть]

[cross-track]

Отличная таблица! Не можете добавить график err,km =  f(Ep-Ep,hours)?

[Олег]

Не можете добавить график

Вот -
В 2 вариантах, ошибка в линейной и в логарифмической шкале. Зависимость есть, но довольно нечеткая.
Точку с ошибкой 181 км , от коррекции, удалил.

PS Попробовал посчитать ошибку от смены ориентации МКС, и от этого изменения балл. коэффициента .
Площадь миделя считал мин. = 30*5=150 м кв, макс. = 30*100 = 3000 м кв.
Но получилось немного - порядка десятков метров за сутки. До сотни.

[cross-track]

Спасибо, отлично!

Из вашей таблицы и графика я выбрал "правильные эпохи" для ТЛЕ до и после коррекции (я понял, что для гарантии ТЛЕ до коррекции лучше брать днем ранее).
ТЛЕ_до        27 августа 17:30,
ТЛЕ_после   28 августа 16:42.

В результате получил:


Максимальное сближение - 28 августа 1:50, 0.7 км.

И без конспирологии обошлось)

[Олег]

О влиянии активности Солнца.
Выбрал за месяц с 5.08 по 5.09 за этот год и за 2019 ( минимум Солнца ). Удалил ошибки при коррекциях.
Пунктирные линии - регрессия.

[cross-track]

Не такая уж большая разница. Возможно дело в том, что минимум активности 2019 был не таким выраженным минимумом.

[Олег]

Решил подробнее посчитать ошибки SGP4 (TLE) как в первом посте
Но более длительные периоды, то есть определить, когда окончательно "протухнут" старые TLE. Это актуально для всяких секретных USA, у которых обычно TLE не первой свежести.
Взял один из самых стойких спутников - геодезический LARETS (27944/03042F). "Стойкий" - в смысле не подвержен труднопредсказуемым факторам - атм. торможению и давлению света.
Вот - про него про LARETS
Взял TLE с мая этого года. Поэтому макс. dt = 138 дней.
Результат - в виде графиков. Это один график в линейной, другой - то же самое в логарифмической шкале.
Штриховая линия - аппроксимация ошибка(км) er = 1.819*dt^1.519 . dt - в днях.
Ошибки наверное вызваны неучетом гармоник выше J4, лунно-солнечными возмущениями, и округлениями значений ЭО от разрядности TLE.

Видна какая-то периодичность непонятная. Если бы ошибка была близка к размеру орбиты, то понятно было бы. Типа отставание стрелочных часов на 23 час похоже на опережение в 1 час.

[Олег]

Про другие объекты. Взял ступени ракет, они не маневрируют. Баллист. коэф. ступеней примерно такие же, как у спутников с панелями.
Период поменьше - с 20 августа, макс. dt 27 дней.
В таблицах в правом столбце - dt(err=100),day - расчетное время, когда ошибка достигнет 100 км. Видно, что на низких орбитах предсказуемость хуже.

[ОАЯ]

Если ошибка возникает резко 12 марта, то это серьезно. 

[cross-track]

Решил подробнее посчитать ошибки SGP4 (TLE) как в первом посте
Но более длительные периоды, то есть определить, когда окончательно "протухнут" старые TLE. Это актуально для всяких секретных USA, у которых обычно TLE не первой свежести.
Взял один из самых стойких спутников - геодезический LARETS (27944/03042F). "Стойкий" - в смысле не подвержен труднопредсказуемым факторам - атм. торможению и давлению света.
Вот - про него про LARETS
Взял TLE с мая этого года. Поэтому макс. dt = 138 дней.
Результат - в виде графиков. Это один график в линейной, другой - то же самое в логарифмической шкале.
Штриховая линия - аппроксимация ошибка(км) er = 1.819*dt^1.519 . dt - в днях.
Ошибки наверное вызваны неучетом гармоник выше J4, лунно-солнечными возмущениями, и округлениями значений ЭО от разрядности TLE.

Видна какая-то периодичность непонятная. Если бы ошибка была близка к размеру орбиты, то понятно было бы. Типа отставание стрелочных часов на 23 час похоже на опережение в 1 час.

Атмосферное торможение для LARETS, конечно, очень мало (он теряет высоту с темпом 1-1.5 метра в день), однако на больших интервалах предсказания орбит это может сказаться довольно ощутимо. Нельзя ли сделать аналогичный расчет, но для B*=0 для всех TLE на этом интервале в 138 дней?

[cross-track]

Атмосферное торможение для LARETS, конечно, очень мало (он теряет высоту с темпом 1-1.5 метра в день), однако на больших интервалах предсказания орбит это может сказаться довольно ощутимо. Нельзя ли сделать аналогичный расчет, но для B*=0 для всех TLE на этом интервале в 138 дней?

Я сделал такую симуляцию для одного спутника с эпохой за 1 мая 2024, положив для него B*=0. Как видно из графика, ошибка на рассматриваемом интервале достигает 2300 км, хотя для этого же спутника с  B*=/=0 ошибка не превышала 400 км.

Симуляция для  B*=0


Симуляция для  B*=4.47*10^-5

[Олег]

Нельзя ли сделать аналогичный расчет, но для B*=0 для всех TLE

Вот, сравнение с исходным.

[cross-track]

Нельзя ли сделать аналогичный расчет, но для B*=0 для всех TLE

Вот, сравнение с исходным.

Из графиков (особенно первого) хорошо видно, что геодезический спутник "Ларец", несмотря на относительно высокую орбиту с минимальным атмосферным трением все же не может рассматриваться как модель без учета дрейфа в остатках атмосферы. И периодичность на графике, о которой вы писали в первом сообщении, может быть обусловлена корреляцией с солнечной и геомагнитной активностью, влияющей на атмосферу (если такая активность тоже менялась похожими "волнами").