ecef2eci

Векторы положения и скорости в центрированном на Земле инерционном среднем-экваторе среднем-равноденствии

Описание

пример

[r_eci] = ecef2eci(utc,r_ecef) вычисляет вектор положения в Координате в геоцентрической инерциальной системе координат (ECI) системе для заданного вектора положения в Earth-centered Earth-fixed (ECEF) системе координат в определенном скоординированном универсальном времени (UTC). Для получения дополнительной информации о центрированной Землей системе фиксированных координат смотрите Алгоритмы.

[r_eci,v_eci] = ecef2eci(___,v_ecef) вычисляет векторы положения и скорости для векторов заданного положения и скорости.

[r_eci,v_eci,a_eci] = ecef2eci(___,a_ecef) вычисляет положение, скорость, векторы ускорения для заданных векторов положения, скорости и ускорения.

[r_eci,v_eci,a_eci] = ecef2eci(___,Name,Value) вычисляет положение, скорость и векторы ускорения с более высокой точностью, используя параметры ориентации Земли.

Примеры

свернуть все

Преобразуйте положение и скорость ECEF в ECI в 12:00 4 января 2019 года.

r_ecef = [-5762640 -1682738 3156028];
v_ecef = [3832 -4024 4837];
utc = [2019 1 4 12 0 0];
[r_eci, v_eci] = ecef2eci(utc, r_ecef, v_ecef);
r_eci =
   1.0e+06 *
   -2.9818
    5.2070
    3.1616

v_eci =
   1.0e+03 *
   -3.3837
   -4.8870
    4.8430

Преобразуйте позицию ECEF в ECI в 12:00 4 января 2019 года, включая эффекты полярного движения.

r_ecef = [-5762640 -1682738 3156028];
utc = [2019 1 4 12 0 0];
mjd = mjuliandate(utc);
pm = polarMotion(mjd, 'action', 'none')*180/pi;
r_eci = ecef2eci(utc, r_ecef, 'pm', pm);
r_eci = 
   1.0e+06 * 
   -2.9818 
   5.2070 
   3.1616

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

свернуть все

UTC в году, месяце, дне, часе, минутах и секундах порядка, заданный как массив 1 на 6 значений UTC:

Значение времениВойти
ГодДвойное значение, которое является целым числом, большим 1, таким как 2013.
МесяцДвойное значение, которое является целым числом, больше 0, в области значений 1 на 12.
ДеньДвойное значение, которое является целым числом, больше 0, в области значений 1 на 31.
ЧасДвойное значение, которое является целым числом, больше 0, в области значений 1 на 24.
Минута и секундаДвойное значение, которое является целым числом, больше 0, в области значений 1 на 60.

Пример: [2000 1 12 4 52 12.4]

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

Массив компонентов положения ECEF, заданный как массив 3 на 1.

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

Компоненты скорости ECEF, заданные как массив 3 на 1.

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

Компоненты ускорения ECEF, заданные как массив 3 на 1.

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

Аргументы в виде пар имя-значение

Задайте необязательные разделенные разделенными запятой парами Name,Value аргументы. Name - имя аргумента и Value - соответствующее значение. Name должны находиться внутри кавычек. Можно задать несколько аргументов в виде пар имен и значений в любом порядке Name1,Value1,...,NameN,ValueN.

Пример: 'dUT1',0.234

Различие между международным атомным временем (TAI) и UTC, заданная в виде скаляра, в секундах.

Пример: 32

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

Различие между UTC и универсальным временем (UT1), заданная в виде скаляра, в секундах.

Пример: 0.234

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

Полярные смещения из-за движения земной коры по оси x - и y - в степенях.

Совет

Чтобы вычислить перемещение, используйте polarMotion функция.

Пример: pm = polarMotion(mjd, 'action', 'none')*180/pi;

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

Корректировка положения Небесного Промежуточного Полюса (CIP), в степенях, задается как разделенная разделенными запятой парами, состоящая из dCIP и массив M -by-2. Это расположение (dDeltaX, dDeltaY) находится вдоль осей x - и y -. По умолчанию эта функция принимает массив нулей 1 на 2.

Исторические значения см. на веб-сайте International Earth Rotation and Reference Systems Service (https://www.iers.org) и перейдите на страницу Earth orientation data Data/Products.

  • M массив -by-2

    Задайте M -by-2 массив значений корректировки местоположения, где M количество матриц косинуса направления или преобразования, которые будут преобразованы. Каждая строка соответствует одному набору dDeltaX и dDeltaY значений.

Пример: [-0.2530e-6 -0.0188e-6]

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

Избыточная продолжительность дня (различие между астрономически определенной длительностью дня и 86400 секунд СИ), заданная в виде скаляра, в секундах.

Пример: 32

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

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

свернуть все

Компоненты положения ECI, заданные как массив 3 на 1.

Компоненты скорости ECI, заданные как массив 3 на 1.

Компоненты ускорения ECI, заданные как массив 3 на 1.

Алгоритмы

ecef2eci функция использует эти ориентированные на Землю системы координат:

  • Earth Centered Inertial Frame (ECI) - Используемая инерционная система координат является Международной небесной системой отсчета (ICRF). Этот кадр можно рассматривать как равный системе координат ECI, реализованной в J2000 (1 января 2000 г. 12:00:00 TT). Для получения дополнительной информации см. раздел «Координаты ECI».

  • Earth-centred Earth-fixed Frame (ECEF) - используемая фиксированная система координат является Международной наземной системой отсчета (ITRF). Эта опорная система координат реализуется сокращением IAU2000/2006 от системы координат ICRF. Для получения дополнительной информации см. раздел «Координаты ECEF».

Ссылки

[1] Валладо, Д. А. Основы астродинамики и применения. альг. 4. Нью-Йорк: McGraw-Hill, 1997.

[2] Gottlieb, R. G., «Fast Gravity, Gravity Partials, Normalized Gravity, Gradient Torque and Magnetic Field: Derivation, Code and Data», Технический отчет НАСА.

[3] Konopliv, A. S., S. W. Asmar, E. Carranza, W. L. Sjogen, D. N. Yuan., «Недавние гравитационные модели в результате миссии лунного изыскателя, Icarus», Vol. 150, no. 1,

[4] Lemoine, F. G., D. E. Smith, D.D. Роулендс, М. Т. Цубер, Г. А. Нойман, и Д. С. Чинн, «Улучшенное решение гравитационного поля Марса (GMM-2B) от Mars Global Surveyor», Journal Of Geophysical Research, Vol. 106, No. E10, стр. 23359-23376, 25 октября 2001 года.

[5] Seidelmann, P.K., Archinal, B.A., A 'hearn, M.F. et al. Доклад Рабочей группы МАС/МАГ по картографическим координатам и элементам ротации: 2006 год. Небесный Mech Dyn Astr 98, 155-180 (2007).

Введенный в R2019a
Для просмотра документации необходимо авторизоваться на сайте