Я строю пропагатор спутниковой орбиты в MATLAB и использую STK в качестве эталона.
Сравнивая пропагатор TwoBody с использованием RKN78 (MATLAB) с HPOP (степень 1, порядок 1), я получаю неплохие результаты в системе координат ICRF, но когда я преобразую результаты MATLAB в инерциальной системе отсчета в фиксированную систему отсчета Земли с учетом прецессии, нутацию и полярное движение Земли и сравнить результаты с HPOP (степень 1, порядок 1) в фиксированной системе отсчета, я получаю ошибку положения в несколько сотен метров за 1 период моделирования орбиты.
Я ищу более точный результат, так есть ли способ поиграть с преобразованием ICRF-Fixed в настройках STK?
AGI заявляет следующее в своей онлайн-документации для преобразования ICRF-Fixed.
«Земля реализует свой Фиксированный кадр из алгоритма преобразования между ней и ICRF. Преобразование включает в себя эффекты прецессии, нутации и вращения, а также отклонения полюсов и коррекции кадра».
STK предоставляет различные способы отчетности TrueOfDate
, MeanOfDate
которые помогут вам сравнить промежуточные этапы в цепочке преобразования ICRF в ECEF. Вам также необходимо убедиться, что ваши константы дрейфа/движения полюса и UTC-UT1 совпадают между STK и вашим кодом Matlab. Несоответствие полярного движения легко приведет к ошибке в несколько сотен метров.
В справке STK есть довольно хороший обзор « Системы отсчета центрального тела », который может вам помочь.
Дэвид Хаммен
Дэвид Хаммен
call2voyage
Ахмад Фарес
Мефитико