Отличия GPS TLE и RINEX

У меня есть навигационный файл TLE и RINEX для GPS спутника PRN 20, на 01.04.17р. 19:56:56 и 20:00:00 соответственно. См. также Независимый формат обмена получателя, версия 3.03 .

Представляется, что значение долготы восходящего узла (Ω0) у обоих различается, тогда как все остальные параметры относительно схожи.

RINEX-файл:

M0 = -0.41907076212790001
e  =  0.0044262305600569997
ω  =  1.5754052384999999
i  =  0.92697366533970005
Ω0 = -0.4446281969834

TLE-файл:

M0 =  5.8389169820721856
e  =  0.0042077
ω  =  1.5770882387483363
i  =  0.92682044735729674
Ω0 =  2.7589292150900362

Что может быть причиной этого?

В то время как значения для M0близки по модулю 2π, значения для Ω0, по-видимому, отличаются только на π, или примерно на 180 °.

При проверке других наборов данных результат аналогичен, хотя разница между значениями различается.

ТЛЕ:

GPS BIIR-4  (PRN 20)    
1 26360U 00025A   17091.83120691 -.00000002  00000-0  00000-0 0  9995
2 26360  53.1029 158.0750 0042077  90.3605 334.5453  2.00565359123815

РИНЭКС:

20 17  4  1 20  0  0.0 4.665344022214D-04 1.591615728103D-12 0.000000000000D+00
5.600000000000D+01-1.231250000000D+01 5.538802141913D-09-4.190707621279D-01
-5.606561899185D-07 4.426230560057D-03 6.150454282761D-06 5.153646507263D+03
 5.904000000000D+05-6.705522537231D-08-4.446281969834D-01 3.725290298462D-08
 9.269736653397D-01 2.442500000000D+02 1.575405238500D+00-8.669646839876D-09
 2.714398779890D-11 0.000000000000D+00 1.942000000000D+03 0.000000000000D+00
 2.000000000000D+00 0.000000000000D+00-8.381903171539D-09 5.600000000000D+01
 5.899500000000D+05 4.000000000000D+00
Как предложение, вы проверили, что PRN20 TLE и RINEX - это один и тот же спутник? Смещение М 0 а также Ом 0 о которых вы сообщаете, типичны для разных спутников созвездия. Кроме того, ваш ярлык обуви TLE и RINEX PRN20, когда вы начали вопрос, спрашивая о PRN1 (может быть, опечатка?)
Где вы взяли файл RINEX?
Файл rinex взят с garner.ucsd.edu , TLE с celestrak.com/NORAD/elements/gps-ops.txt . @Julio Джулио, как я могу проверить, является ли PRN20 одним и тем же спутником для обоих файлов? Я знаю, что номера спутников различаются для файлов ГЛОНАСС rinex и tle, но я думал, что для GPS они одинаковы. Да, PRN1 — это опечатка.
@MartynaKasyr, на вашем месте я буду искать другой спутник в вашем TLE, который соответствует текущим значениям элементов орбиты вашего файла RINEX, чтобы едва подтвердить (или нет) несоответствие чисел PRN между TLE и RINEX. Здесь также есть открытый вопрос по этой теме space.stackexchange.com/questions/20887/…

Ответы (2)

Они находятся в разных системах отсчета. Набор элементов TLE определяется относительно системы отсчета истинного экватора, среднего равноденствия ( TEME ), а элементы в сообщении RINEX Nav определяются относительно центра Земли, фиксированной точки Земли (ECEF). Поскольку в обоих кадрах используется истинный экватор, другие элементы (например, наклон) различаются не так сильно, как расположение узла.

Они находятся в разных системах отсчета, но это еще не все. Даже после преобразования TEME в ECEF различия все равно будут, потому что вещи с одинаковыми именами на самом деле имеют очень разные значения. RINEX — это просто формат файла; на пользователей возлагается бремя понимания того, что сообщения GPS-навигации и TLE используют одни и те же имена полей для разных целей.

Спецификация космического сегмента/навигационного пользовательского интерфейса NAVSTAR GPS IS-GPS-200 (редакция M, от мая 2021 г.), раздел 20.3.3.4, «Подкадры 2 и 3», содержит 12 страниц инструкций о том, как должен работать пользовательский сегмент. расчеты эфемерид, которые не совсем совпадают с тем, как кто-либо еще использует эти термины. В частности, обратите внимание на эти утверждения (выделено мной):

Параметры эфемерид описывают орбиту в течение интервалов подбора кривой, описанных в разделе 20.3.4. В таблице 20-II дано определение параметров орбиты с использованием терминологии, типичной для кеплеровских параметров орбиты; однако следует отметить, что передаваемые значения параметров таковы, что они обеспечивают наилучшее согласование траектории в геоцентрических координатах с привязкой к Земле (ECEF) для каждого конкретного интервала соответствия. Пользователь не должен интерпретировать промежуточные значения координат как относящиеся к какой-либо общепринятой системе координат. Пользователь должен вычислить координаты ECEF положения фазового центра антенн КА, используя вариант уравнений, показанных в Таблице 20-IV. Параметры подкадров 2 и 3 кеплеровы по внешнему виду.; значения этих параметров, однако, получаются... путем аппроксимации кривой наименьших квадратов распространяющихся эфемерид фазового центра антенн КА

Фитнес-интервалы длятся три-четыре часа. Числа в навигационных сообщениях GPS дают значения, которые обеспечивают наилучшее соответствие всему этому интервалу, поэтому они не совпадают с оскулирующими элементами для орбит GPS в любую конкретную эпоху. Вы должны прочитать документы, чтобы узнать, как с ними обращаться, и вы должны использовать уравнения, которые они вам говорят (включая, например, особый способ итерации уравнения Кеплера для решения), или ответ, который вы получите, не будет иметь такая же статистика ошибок, что эти таблицы существуют, чтобы позволить вам достичь.

Эти предупреждения еще более верны для наборов двухстрочных элементов. Эти цифры еще дальше от оскулирования. Их называют средними элементами. Можно возразить, что это слово должно относиться к тому, как это делает GPS, но это слово не используется в орбитальной механике. «Среднее» обычно означает просто «среднее», но есть много способов что-то усреднить, и техническое использование этого термина в астродинамике является более конкретным и более общим, чем вы могли бы ожидать. По словам Герберта Уолтера (1967), «под средними элементами мы понимаем соприкасающиеся элементы, из которых вычтены короткопериодические и долгопериодические возмущения земного потенциала». Чтобы преобразовать средние элементы в соприкасающиеся, вы должны точно знать, какие возмущения были выбраны для вычитания, и вам нужно вычислить эти поправки и добавить их обратно, иначе вы получите неправильный ответ.

Элементы орбиты в TLE указаны в «форме Кодзаи», что означает использование точного списка возмущений, выбранных для удаления в Yoshihide Kozai, The Motion of a close earth Satellite , Astronomical Journal 64 (367-377) 1959. Например, значение быть помещенным в поле с надписью «большая полуось» не является обычным а (который в связанной статье называется а 0 ), а скорее "средний" а ¯ , что равно а 0 раз фактор

1 А 2 п 2 ( 1 3 2 грех 2 я ) ( 1 е 2 )

куда п знак равно а ( 1 е 2 ) а также А 2 - коэффициент гравитационного расширения, равный 3 2 Дж 2 р 2 , куда Дж 2 знакомая вещь и р - экваториальный радиус Земли.

Это, безусловно, самая простая такая формула. Эквивалентные выражения для аргумента перигея ( ю ¯ ) и прямое восхождение восходящего узла ( Ом ¯ ) запустить на несколько страниц. Обратите внимание, что преобразование из Козаи А 2 к ныне общему Дж 2 Я получил от Брауэра (1959) , который описывает другой и гораздо более сложный набор средних элементов и содержит таблицу, показывающую, как переводить девять различных форм коэффициентов гравитационного расширения. Статья Брауэра заканчивается утверждением:

Я теперь жалею, что представил к 2 , к 4 в моей статье 1946 г. Основная причина заключалась в том, что они дают особенно простую форму для выражения потенциала в экваториальной плоскости. Если бы я мог предвидеть рост интереса к предмету и замешательство, которому я способствовал, я бы выбрал... альтернативную форму, которую использовал Винти (1959) . Я намерен вернуться к этой форме и рекомендовать ее другим авторам.

Заявление Брауэра, по-видимому, сработало, поскольку астродинамика, по-видимому, твердо остановилась на форме, которую Винти приписывал Р.Х. Мерсон и Д. Г. Кинг-Хеле , Использование искусственных спутников для исследования гравитационного поля Земли: результаты спутника-2 (1957β), Nature 182 (640) 1958, знакомый

В ( р , θ ) знак равно грамм М р [ 1 п знак равно 2 Дж п ( р р ) п п п ( грех θ ) ]