ECEF Position to LLA

Вычислите геодезическую широту, долготу и высоту выше планетарного эллипсоида от положения Сосредоточенного землей зафиксированного землей (ECEF)

  • Библиотека:
  • Aerospace Blockset / Утилиты / Преобразования Осей

  • ECEF Position to LLA block

Описание

Положение ECEF с блоком LLA преобразует вектор 3 на 1 из положения ECEF (p¯) в геодезическую широту (μ¯), долгота (ι¯), и высота (h¯) выше планетарного эллипсоида. Для получения дополнительной информации о положении ECEF см. Алгоритмы.

Ограничения

  • Эта реализация генерирует геодезическую широту, которая находится между ±90 градусами и долготой, которая находится между ±180 градусами. Планета принята, чтобы быть эллипсоидальной. Путем установки выравнивания на 0, вы моделируете сферическую планету.

  • Реализация системы координат ECEF принимает, что ее источник находится в центре планеты, x - ось пересекает начало (Гринвич) меридиан и экватор, z - ось является средней осью вращения планеты (положительный на север), и y - ось завершает предназначенную для правой руки систему.

Порты

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

развернуть все

Положение в ECEF структурирует в виде вектора 3 на 1.

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

Вывод

развернуть все

Геодезическая широта и долгота, возвращенная как 2 1 вектор, в градусах.

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

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

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

Параметры

развернуть все

Устройства вывода в виде:

Модули

Положение

Экваториальный радиус

Высота

Metric (MKS)

Метры

Метры

Метры

English

Футы

Футы

Футы

Зависимости

Чтобы включить этот параметр, установите модель Planet на Earth (WGS84).

Программируемое использование

Параметры блоков: units
Ввод: символьный вектор
Значения: 'Metric (MKS)' | 'English'
Значение по умолчанию: 'Metric (MKS)'

Модель Planet, чтобы использовать, Custom или Earth (WGS84).

Программируемое использование

Параметры блоков: ptype
Ввод: символьный вектор
Значения: 'Earth (WGS84)' | 'Custom'
Значение по умолчанию: 'Earth (WGS84)'

Выравнивание планеты в виде двойного скаляра.

Зависимости

Чтобы включить этот параметр, установите модель Planet на Custom.

Программируемое использование

Параметры блоков: F
Ввод: символьный вектор
Значения: двойной скаляр
Значение по умолчанию: '1/298.257223563'

Радиус планеты в ее экваторе в виде двойного скаляра, в тех же модулях как желаемые модули для положения ECEF.

Зависимости

Чтобы включить этот параметр, установите модель Planet на Custom.

Программируемое использование

Параметры блоков: R
Ввод: символьный вектор
Значения: двойной скаляр
Значение по умолчанию: '6378137'

Алгоритмы

Положение ECEF задано как:

p¯=[p¯xp¯yp¯z].

Долгота вычисляется от положения ECEF

ι=atan(pypx).

Геодезическая широта (μ¯) вычисляется от положения ECEF с помощью метода Боуринга, который обычно сходится после двух или трех итераций. Метод начинается с исходного предположения для геодезической широты (μ¯) и уменьшаемая широта (β¯). Исходное предположение принимает форму:

β¯=atan(pz(1f)s)μ¯=atan(pz+e2(1f)(1e2)R(sinβ)3se2R(cosβ)3)

где R является экваториальным радиусом, f является выравниванием планеты, e 2 = 1− (1−f) 2, квадрат первого эксцентриситета, и:

s=px2+py2.

После того, как исходные предположения вычисляются, уменьшаемая широта (β¯) повторно вычисляется с помощью

β=atan((1f)sinμcosμ)

и геодезическая широта (μ¯) переоценен. Этот последний шаг повторяется до μ¯ сходится.

Высота (h¯) выше планетарного эллипсоида вычисляется с

h=scosμ+(pz+e2Nsinμ)sinμN,

где радиус искривления в вертикальном начале (N¯) дают

N=R1e2(sinμ)2.

Ссылки

[1] Стивенс, B. L. и Ф. Л. Льюис. Управление самолетом и симуляция, Хобокен, NJ: John Wiley & Sons, 1992.

[2] Zipfel, Питер Х., моделирование и симуляция космической динамики аппарата. Второй выпуск. Рестон, ВА: образовательный ряд AIAA, 2000.

[3] Методические рекомендации для атмосферного и систем координат транспортного средства космического полета, R-004-1992, ANSI/AIAA, февраль 1992.

Расширенные возможности

Генерация кода C/C++
Генерация кода C и C++ с помощью Simulink® Coder™.

Представлено до R2006a