Трансформация между ICRF и Fixed в STK

Я строю пропагатор спутниковой орбиты в MATLAB и использую STK в качестве эталона.

Сравнивая пропагатор TwoBody с использованием RKN78 (MATLAB) с HPOP (степень 1, порядок 1), я получаю неплохие результаты в системе координат ICRF, но когда я преобразую результаты MATLAB в инерциальной системе отсчета в фиксированную систему отсчета Земли с учетом прецессии, нутацию и полярное движение Земли и сравнить результаты с HPOP (степень 1, порядок 1) в фиксированной системе отсчета, я получаю ошибку положения в несколько сотен метров за 1 период моделирования орбиты.

Я ищу более точный результат, так есть ли способ поиграть с преобразованием ICRF-Fixed в настройках STK?

AGI заявляет следующее в своей онлайн-документации для преобразования ICRF-Fixed.

«Земля реализует свой Фиксированный кадр из алгоритма преобразования между ней и ICRF. Преобразование включает в себя эффекты прецессии, нутации и вращения, а также отклонения полюсов и коррекции кадра».

Я не знаю, как не быть прямолинейным по этому поводу, поэтому я просто должен быть прямолинейным по этому поводу: сравнение одного метода распространения состояния (положения + скорости) с другим - игра дураков, если только вы не знаете истинное состояние. В группе, разработавшей ЗАО «Инженерная орбитальная динамика», был человек, занимавшийся выяснением истинного состояния различных транспортных средств с целью сравнения. Его задача была нелегкой, как и задача тех, кто пытался воспроизвести эти орбиты, используя различные методы распространения.
Ни его работа, ни работа разработчиков этих оценок не подлежат ответственности в ограниченном пространстве этого веб-сайта.
@DavidHammen Технически, воспринимая вопрос буквально, «есть ли способ» требует только простого ответа «Да» или «Нет», поэтому ответ, говорящий что-то в этом роде, будет действительным: «Да, есть способ, используемый АО , но это очень сложный подход по [причинам], поэтому его нельзя полностью указать здесь. Вот некоторые ресурсы, которые вы можете использовать, чтобы узнать больше об этом: [ресурсы]".
Из ваших ответов я могу сказать, что вы не поняли мой вопрос, поэтому позвольте мне уточнить мою точку зрения. Я спрашиваю, есть ли в графическом интерфейсе STK параметры, позволяющие пользователю контролировать способ выполнения преобразования кадра с инерционно-фиксированной фиксацией. Примером может служить флажок для включения или исключения нутации земли и прецессии из алгоритма преобразования.
Я вижу, что через некоторое время у этого поста есть комментарии, но нет ответа. Удалось ли вам продвинуться в этой теме? Есть ли у вас общий код на GitHub или аналогичном?

Ответы (1)

STK предоставляет различные способы отчетности TrueOfDate, MeanOfDateкоторые помогут вам сравнить промежуточные этапы в цепочке преобразования ICRF в ECEF. Вам также необходимо убедиться, что ваши константы дрейфа/движения полюса и UTC-UT1 совпадают между STK и вашим кодом Matlab. Несоответствие полярного движения легко приведет к ошибке в несколько сотен метров.

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