moonLibration

Лунные колебания

Синтаксис

angles= moonLibration(ephemerisTime)
angles= moonLibration(ephemerisTime,ephemerisModel)
angles= moonLibration(ephemerisTime,ephemerisModel,action)
[angles,rates] = moonLibration(___)

Описание

пример

angles= moonLibration(ephemerisTime) реализует Лунные углы колебания для ephemerisTime, выраженного в дни Джулиана.

Функция использует Чебышевские коэффициенты, которые обеспечивает Лаборатория реактивного движения НАСА.

Эта функция требует, чтобы вы загрузили эфемеридные данные Add-On Explorer. Для получения дополнительной информации смотрите aeroDataPackage.

angles= moonLibration(ephemerisTime,ephemerisModel) использует коэффициенты ephemerisModel, чтобы реализовать эти значения.

angles= moonLibration(ephemerisTime,ephemerisModel,action) использование action, чтобы определить сообщение об ошибке.

пример

[angles,rates] = moonLibration(___) реализует Лунные углы колебания и уровни с помощью любой комбинации входных параметров в предыдущих синтаксисах.

Примеры

свернуть все

Реализуйте углы колебания Луны на 1 декабря 1990 с DE405. Используйте функцию juliandate, чтобы вычислить значение даты входа Julian.

angles = moonLibration(juliandate(1990,12,1))
angles =
   1.0e+03 *
    0.0001    0.0004    1.8010

Задайте ephemerides (DE421) и используйте функцию juliandate для даты (1 января 2000), чтобы вычислить и Лунные углы колебания и уровни.

[angles,rates] = moonLibration([2451544.5 0.5],'421')
angles =
   1.0e+03 *
   -0.0001    0.0004    2.5643

rates =
   -0.0001    0.0000    0.2301

Входные параметры

свернуть все

Юлианские даты, для которых вычисляются положения, задали как одно из следующего:

  • Скаляр

    Укажите, что тот назначил дату Джулиана.

  • Вектор с 2 элементами

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

  • Вектор-столбец

    Задайте вектор-столбец с элементами M, где M является количеством дат Джулиана.

  • M-by-2 матрица

    Задайте матрицу, где M является количеством дат Джулиана, и второй столбец содержит прошедшие дни (Опорная дата Джулиана / протек дневные пары).

Типы данных: double

Коэффициенты Ephemerides, заданные как один из этих ephemerides, заданных Лабораторией реактивного движения:

  • '405'

    Выпущенный в 1 998. Этот ephemerides учитывает диапазон дат Джулиана 2305424.50 (9 декабря 1599) к 2 525 008,50 (20 февраля, 2201).

    Эта функция вычисляет эти ephemerides относительно Международной Астрономической Ссылочной версии 1.0 Кадра, принятой в 1 998.

  • '421'

    Выпущенный в 2 008. Этот ephemerides учитывает диапазон дат Джулиана 2414992.5 (4 декабря 1899) к 2 469 808,5 (2 января 2050).

    Эта функция вычисляет эти ephemerides относительно Международной Астрономической Ссылочной версии 1.0 Кадра, принятой в 1 998.

  • '423'

    Выпущенный в 2 010. Этот ephemerides учитывает диапазон дат Джулиана 2378480.5 (16 декабря 1799) к 2 524 624,5 (1 февраля, 2200).

    Эта функция вычисляет эти ephemerides относительно Международной Астрономической Ссылочной версии 2.0 Кадра, принятой в 2 010.

  • '430'

    Выпущенный в 2 013. Этот ephemerides учитывает диапазон дат Джулиана 2287184.5 (21 декабря 1549) к 2 688 976,5 (25 января, 2650).

    Эта функция реализует эти ephemerides относительно Международной Астрономической Ссылочной версии 2.0 Кадра, принятой в 2 010.

  • '432t'

    Выпущенный в апреле 2014. Этот ephemerides учитывает диапазон дат Джулиана 2287184.5, (21 декабря 1549) к 2 688 976,5, (25 января, 2650).

    Эта функция реализует эти ephemerides относительно Международной Астрономической Ссылочной версии 2.0 Кадра, принятой в 2 010.

Типы данных: char

Функциональное поведение, когда входные параметры вне области значений, заданной как одно из этих значений:

ЗначениеОписание
'None'Никакое действие.
'Warning'Предупреждая в Командном окне MATLAB®, симуляция модели продолжается.
'Error'MATLAB возвращает исключение, остановки симуляции модели.

Типы данных: char

Выходные аргументы

свернуть все

Лунные углы колебания, заданные как M-by-3 числовой массив. M является количеством дат Джулиана в строках. Столбцы содержат Углы Эйлера (φ θ ψ) для Лунного отношения в радианах.

Если входные параметры включают несколько дат Джулиана или эпохи, этот массив имеет одинаковое число строк как вход ephemerisTime.

Лунное колебание угловые уровни, заданные как M-by-3 числовой массив. M является количеством дат Джулиана в строках. Столбцы содержат Лунное колебание Эйлеровы угловые уровни (ω) в радианах/день.

Если входные параметры включают несколько дат Джулиана или эпохи, этот массив имеет одинаковое число строк как вход ephemerisTime.

Ссылки

[1] Folkner, W. M. Дж. Г. Уильямс, Д. Х. Боггс, “Планетарная и лунная эфемерида DE 421”, межпланетный сетевой отчет о выполнении работ JPL 24-178, 2009.

[2] Vallado, D. A. основные принципы астродинамики и приложений, McGraw-Hill, Нью-Йорк, 1997.

Введенный в R2013a