Почему кажется, что элементы орбиты LRO постоянно колеблются - за исключением недавнего времени?

Я нанес на карту некоторые элементы орбиты лунного разведывательного орбитального аппарата из базы данных JPL Horizons и вижу постоянные колебания. Период большой полуоси и параметры, связанные с эксцентриситетом, составляют около 27,25 дня (см. Последний график), что соответствует периоду обращения около 27,32 дня. Однако колебания наклона колеблются с удвоенной частотой.

Почему элементы орбиты LRO постоянно колеблются?

Кроме того, являются ли внезапные эпизоды очень постоянного эксцентриситета «реальными» или просто артефактами сращивания/сшивания/склеивания разных симуляций вместе? Даже когда эксцентриситет кажется постоянным, наблюдаются колебания наклона и, в меньшей степени, большой полуоси.

edit: я только что добавил некоторую информацию об истории различных траекторий, которые были объединены. «Период спокойной эксцентричности» начинается с середины 2016-Oct-21, но заканчивается в ту же печально известную дату, что и этот отрезок.2016-Dec-07558day_20160907_01.bsp_V0.2

Но помните, я спрашиваю о самих шевелениях не меньше, чем об их отсутствии.

введите описание изображения здесь

SPACECRAFT TRAJECTORY: 
  Updated irregularly (on Horizons) or by request. 
  Concatenated historical (reconstructed) trajectories are from PDS. 

  Trajectory name                 Start (TDB)         Stop (TDB)
  ---------------------------  -----------------  -----------------
  Reconstructed trajectory     2009 Jun 18 22:16  2016 Sep 15 00:01
  558day_20160907_01.bsp_V0.2  2016 Sep 15 00:01  2016 Dec 07 00:01 predict
  558day_20161207_01.bsp_V0.2  2016 Dec 07 00:01  2017 Jan 04 00:01 predict
  558day_20170104_01.bsp_V0.1  2017 Jan 04 00:01  2018 Jul 16 00:01 predict
  558day_20170216_01.bsp_V0.1  2017 Feb 16 00:01  2018 Aug 28 00:01 predict

введите описание изображения здесь

введите описание изображения здесь

Внизу: график зависимости эксцентриситета от времени (дни) для двух интервалов, сдвинутых на 327 дней, показывает разницу в 12 колебаний. Извлеченный период составляет 27,25, что близко к периоду обращения Луны около 27,32 дня.

введите описание изображения здесь

Похоже активировалась активная стабилизация с помощью RCS.
Мои деньги на прерывание телеметрии, которое сделало бы чистую синусоиду в данных наклона артефактом фильтра Калмана.
@Schlusstein База данных JPL Horizons содержит расчетные орбиты, результаты численного моделирования траекторий космических кораблей на основе моделей гравитации и других сил. Учтите, что данные распространяются на 2018 год — для получения телеметрии из будущего не использовались машины времени! Теперь я возьму эти деньги - сколько ты сказал? :)
@uhoh Да, и будущие данные выглядят точно так же, как данные, о которых спрашивают, и это большая часть того, как я пришел к своему выводу.
@Schlusstein нет фильтра Калмана. Я не могу понять смысл вашего комментария. Как вы можете заключить, что фильтр Калмана существует? Они обычно используются для зашумленных измеренных данных , чтобы улучшить интерпретацию измерений при наличии значительного количества экспериментального шума. Это точное численное моделирование, я ожидаю, что любой шум будет выше пятнадцатого знака после запятой или около того.
@uhoh Какие входные данные использовались для выполнения этих симуляций? Где вы нашли данные, которые вы построили?
@uhoh Я читаю о ГОРИЗОНАХ. Это не просто численное моделирование. Данные, полученные из фактических (предположительно отфильтрованных) измерений, с использованием численного моделирования для заполнения пробелов в данных. Это, конечно, намного точнее, чем численное моделирование, поскольку позволяет скорректировать накопленную ошибку и негравитационные орбитальные возмущения. ssd.jpl.nasa.gov/pub/ssd/Horizons_doc.pdf Я также рекомендую вам узнать больше о фильтрации Калмана, так как числовое интегрирование может быть легко встроено.
@uhoh И да, я действительно посмотрел на вопрос. Однако за прошедшие четыре дня я забыл, что это веб-инструмент.
@Schlusstein, давайте придерживаться этого вопроса в том виде, в котором он был задан. Согласно выходным данным Horizons, запуск 2016 Sep 15 00:01этих данных является чисто прогнозным, если я правильно понимаю. Таким образом, все, что показано в первом наборе графиков, начиная с 2016.7, является строго численным интегрированием (релятивистски верная версия) Ф знак равно м а . Разве это не правильно?
@uhoh Я не уверен, что это правда. Вы не заметите, что области странного поведения заканчиваются точно на датах окончания этих строк: «558day_20160907_01.bsp_V0.2 2016 15 сентября 00:01 2016 07 декабря 00:01 прогнозировать». Я считаю, что это вряд ли совпадение. Возможно, конечные даты — это то место, где заканчиваются фазы предсказания. Однако я также не вижу плоских участков в необработанных данных.
@Schlusstein Обсуждение с вами очень полезно - мы делали одно и то же в одно и то же время! Пока вы писали свой комментарий, я обновил вопрос той же информацией. Я посмотрел выходной файл Horizons, сохраненный на моем диске - в каждой строке сразу после JD и календарной даты первое число с плавающей запятой - это эксцентриситет, и он резко начинает зависать примерно на 3,23 или 3,24E-02 на время флэта. спотовые даты (от JD 2457682.5до 2457729.5).
Что вы имеете в виду, говоря: «Я спрашиваю о самих шевелениях столько же, сколько и об их отсутствии»? Относится ли это к периодическому поведению значений наклона в течение времени, когда эксцентриситет практически фиксирован?
@Schlusstein Меня так же интересует, например, почему наклон LRO колеблется с амплитудой 1 градус каждые 14 дней, как и то, почему это колебание внезапно падает до 0,1 градуса на некоторое время, а затем начинается снова. Я спрашиваю об этих экскурсиях не меньше, чем об их внезапной остановке и начале.

Ответы (2)

JPL HORIZONS показывает орбиты, интерполированные из реальных данных, а не из чистого моделирования. Я подозреваю, что здесь произошло то, что нет орбитальных данных для трех участков данных, которые не демонстрируют колебания. Мы точно знаем, что для одного из этих разделов нет необработанных данных, потому что он находится в будущем. Я разговаривал с экспертом по контролю, и он сказал, что мои подозрения, вероятно, верны.

Когда вы используете веб-инструмент, он говорит, что существует несколько рядов данных, и все, кроме первого, имеют слово «прогноз» в конце строки. Даты в конце этих рядов данных соответствуют возвращенным колебаниям данных. Я подозреваю, что новая серия данных начинается, когда контакт с космическим кораблем возобновляется.

Я считаю, что период колебаний ~ 14 дней соответствует половине периода обращения Луны. В этом случае основной движущей силой такого поведения будет гравитационное возмущение Земли. Тот факт, что эти осцилляции происходят довольно заметно в отсутствие реальных данных, указывает на то, что модель хорошо справляется с учетом этого возмущения. Что касается того, почему амплитуда сильно падает, я не могу сказать точно, так как недостаточно знаю используемую модель, но могу предположить, что это связано с неучтенной неравномерной плотностью Луны и негравитационными возмущениями, связанными с его ориентация относительно солнца, например световое давление и выделение газа.

Если вы посмотрите на данные об эксцентриситете, вы увидите колебание с периодом ~ 14 дней поверх колебания с периодом ~ 27 дней с некоторыми рассеянными высшими гармониками. В данных по перицентру и апоапсису также наблюдаются сильные колебания периода ~27 дней. Было бы интересно применить ко всему этому преобразование Фурье. Вы могли бы посмотреть на относительные величины возмущения по периодам и, возможно, попытаться выделить некоторое поведение с периодом в один год.

«Интерполированный» — не совсем правильное слово, может быть лучше «подогнанный». Обычно единственные высокоточные данные, которые мы получаем от космических аппаратов, вращающихся вокруг других тел, — это доплеровская задержка. Специально закодированный/модулированный сигнал передается из сети дальнего космоса (DSN) на космический корабль, который ретранслирует его обратно на Землю на другой частоте, но тщательно привязанной по фазе к входящей несущей. Доплеровский сдвиг и абсолютная задержка дают относительную скорость и расстояние по отношению к наземной станции, но это не истинные векторы состояния относительно тела, находящегося на орбите.
Моделирование используется для получения смоделированных результатов типа доплеровского запаздывания, и они итеративно корректируются до тех пор, пока они не смогут воспроизвести измеренные данные доплеровского запаздывания. Иногда также могут быть включены данные космического корабля (камеры, акселерометры), но обычно они не обеспечивают достаточную точность, чтобы соответствовать данным доплеровского запаздывания. Три наземные станции DSN днем ​​и ночью заняты наблюдением за космическими кораблями. Вот несколько крутых данных доплеровского запаздывания с космического корабля вокруг Луны, но в данном случае это пассивный радар.

Основываясь на некоторых очень полезных ответах Джона Джорджини из JPL и небольшом количестве дополнительного чтения, я могу резюмировать следующее:

Фактический расчет орбит космического корабля, сочетающий в себе как эффекты гравитации, так и другие силы в Солнечной системе, а также орбитальные маневры самого космического корабля, сочетаются / согласовываются с доступными данными радара и телеметрии. Результатом, конечно же, является таблица векторов состояний.

Соприкасающиеся элементы получены из векторов состояния только для удобства. Они не предназначены для использования в качестве основных источников данных об орбите, и каждый набор точек предназначен только для определения положения космического корабля в данный момент времени .

В течение короткого периода времени в ожидании либо дополнительных данных, либо решения о предстоящем маневре, альтернативный алгоритм для оскулирующих элементов, по-видимому, был применен к короткому промежутку дат, что привело к более плавному отображению некоторых графиков. .

Это не означает, что они обязательно менее точны в определении положения, потому что элементы предназначены для использования только для определения положения вблизи момента, для которого они рассчитаны. Начнем с того, что истинная орбита космического корабля на самом деле вовсе не коническая, особенно если учесть сложное гравитационное поле системы Земля-Луна.

После обновленного расчета внешний вид вернулся к ожидаемому, хотя это все еще в равной степени справедливо.

Что касается самих периодических вариаций: месячная и двухмесячная периодичность являются проявлением реальности обращения вокруг неравномерного гравитационного поля Луны под действием гравитационного воздействия Земли. Поскольку орбита Луны скорее эллиптическая, чем круговая, этот сильно возмущающий эффект является периодическим. Эти эффекты особенно заметны на этих графиках, потому что они не являются простыми кеплеровскими орбитами с самого начала, и попытка выразить их как таковые подчеркнет разницу.

Ниже приведено краткое сравнение полученных соприкасающихся элементов до и после регулярного запланированного обновления. Все хорошо.

введите описание изображения здесь