exponenta event banner

Модель сферической гармонической гравитации

Реализовать сферическое гармоническое представление планетарной гравитации

  • Библиотека:
  • Аэрокосмический блок/Среда/Гравитация

  • Spherical Harmonic Gravity Model block

Описание

Блок модели сферической гармонической гравитации реализует математическое представление сферической гармонической планетарной гравитации на основе планетарного гравитационного потенциала. Он обеспечивает удобный способ описания гравитационного поля планеты вне её поверхности в сферическом гармоническом расширении.

Сферические гармоники можно использовать для изменения величины и направления сферической гравитации (-GM/r2). Наиболее значимым или наибольшим сферическим гармоническим членом является зональная гармоника второй степени, J2, на которую приходится косость планеты.

Используйте этот блок, если требуется получить более точные значения силы тяжести, чем сферические модели силы тяжести. Например, применение в условиях, не связанных с полетом в атмосфере, может потребовать более высокой точности.

Ограничения

  • Блок исключает центробежные эффекты планетарного вращения и эффекты прецессирующей системы отсчета.

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

Порты

Вход

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

Координаты фиксированного кадра от центра планеты, заданные в виде матрицы N-by-3, в выбранных единицах. Каждая строка матрицы является отдельной позицией для вычисления. Ось Z положительна к Северному полюсу. Если модель центрального тела имеет значение EGM2008 или EGM96, эта матрица содержит координаты земного центра (ECEF).

При вводе большой матрицы фиксированного кадра и значения высокой степени может возникнуть ошибка нехватки памяти. Дополнительные сведения об устранении ошибок нехватки памяти в среде MATLAB ® см. в разделе Устранение ошибок нехватки памяти.

При вводе большой матрицы фиксированного кадра можно получить ограничение максимального размера матрицы. Чтобы определить самую большую матрицу или массив, которые можно создать в среде MATLAB для вашей платформы, см. раздел Производительность и память.

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

Продукция

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

Массив значений силы тяжести по осям X, Y и Z координат фиксированного кадра в выбранных единицах длины в квадрате. Каждая строка матрицы возвращает вычисленный вектор гравитации для соответствующей строки входной матрицы.

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

Параметры

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

Единицы ввода и вывода, указанные как:

Единицы

Вход

Продукция

Metric (MKS)

Метры (м)

Meters/sec2 (м/с2)

English

Футы (фут)

Feet/sec2 (фут/с2)

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

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

Поведение входа вне диапазона, указанное как:

СтоимостьОписание

None

Никаких действий.

Warning

При появлении предупреждения в программе Diagnostic Viewer моделирование модели продолжается.

Error

MATLAB возвращает исключение, моделирование модели останавливается.

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

Параметр блока: action
Текст: символьный вектор
Значения: 'None' | 'Warning' | 'Error'
По умолчанию: 'Warning'

Планетарная модель, указанная как:

Модель центрального кузоваПримечания
EGM2008

Земля - новейшая земная сферическая гармоническая гравитационная модель от Национального агентства геопространственной разведки (NGA). Этот блок предоставляет WGS-84 версию этой гравитационной модели. Можно использовать EGM96 планетарную модель, если нужно использовать старый стандарт для Земли.

EGM96 Земля
LP100K

Луна - лучше всего для определения лунной орбиты на основе вычислительного времени, необходимого для вычисления орбит. Эта модель планеты была создана примерно в том же году, что и LP165P с аналогичными данными.

LP165P

Луна - лучше всего подходит для увеличения точности орбиты лунного полета. Эта модель планеты была создана примерно в том же году, что и LP100K с аналогичными данными.

GMM2B

Марс

Custom

Позволяет указать собственную планетарную модель. Этот параметр включает параметр MAT-файла центрального тела.

EIGENGL04C

Земля - Поддерживает модель гравитационного поля, EIGEN-GL04C (http://icgem.gfz-potsdam.de/tom_longtime). Эта модель является обновлением до EIGEN-CG03C.

Дополнительные сведения о системе координат фиксированного кадра для центральных тел см. в разделе Алгоритмы.

При определении собственной планетарной модели параметр Градус (Degree) ограничивается максимальным значением для int16. При вводе большой степени может возникнуть ошибка нехватки памяти. Дополнительные сведения об устранении ошибок нехватки памяти в среде MATLAB см. в разделе Устранение ошибок нехватки памяти.

Зависимости

Установка для этого параметра значения Custom включает MAT-файл центрального тела.

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

Параметр блока: ptype
Текст: символьный вектор
Значения: 'EGM2008' | 'EGM96' | 'LP100K' | 'LP165P' | 'GMM2B' | 'Custom' | 'EIGENGL04C'
По умолчанию: 'EGM2008'

Степень гармонической модели, заданная как скаляр:

Модель центрального кузоваРекомендуемая степеньМаксимальная степень

EGM2008

120

2159

EGM96

70

360

LP100K

60

100

LP165P

60

165

GMM2B

60

80

EIGENGL04C

70

360

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

Параметр блока: degree
Текст: символьный вектор
Значения: scalar
По умолчанию: '120'

MAT-файл центрального тела, содержащий определения для пользовательской планетарной модели. aerogmm2b.mat файл в Aerospace Blockset™ является MAT-файлом по умолчанию для пользовательской планетарной модели.

Этот файл должен содержать:

ПеременнаяОписание
Ре

Скаляр экваториального радиуса планеты в метрах (м).

GM

Скаляр планетарного гравитационного параметра в метрах в кубах в секунду в квадрате (м3/с2)

степень

Скаляр максимальной степени.

C

(степень + 1) -by- (степень + 1) матрица, содержащая нормализованную матрицу сферических гармонических коэффициентов, С.

S

(степень + 1) -by- (степень + 1) матрица, содержащая нормализованную матрицу сферических гармонических коэффициентов, S.

При использовании большого значения градуса может возникнуть ошибка нехватки памяти. Дополнительные сведения об устранении ошибок нехватки памяти в среде MATLAB см. в разделе Устранение ошибок нехватки памяти.

Зависимости

Чтобы включить этот параметр, задайте для модели центрального тела значение Custom.

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

Параметр блока: datafile
Текст: символьный вектор
Значения: 'aerogmm2b.mat' | MAT-файл
По умолчанию: 'aerogmm2b.mat'

Алгоритмы

Блок сферической гармонической гравитации работает в системе координат фиксированного кадра для центральных тел:

  • Земля - система координат с фиксированным кадром представляет собой систему координат с фиксированным земным ядром (ECEF).

  • Луна - система координат фиксированного кадра является системой главной оси (PA), ориентация определяется планетарной эфемеридной DE403 JPL.

  • Марс - система координат фиксированного кадра определяется направлениями полюсов вращения и простых меридианов, определенных в [14].

Ссылки

[1] Готлиб, Роберт Г., «Быстрая гравитация, гравитационные части, нормализованная гравитация, гравитационный градиентный крутящий момент и магнитное поле: деривация, код и данные». NASA-CR-188243. Хьюстон, Техас: Космический центр им. Линдона Джонсона НАСА, февраль 1993 года.

[2] Вальядо, Давид. Основы астродинамики и приложений. Нью-Йорк: Макгроу-Хилл, 1997.

[3] «Всемирная геодезическая система Министерства обороны 1984, ее определение и связь с местными геодезическими системами». NIMA TR8350.2.

[4] Коноплив, А.С., У. Асмар, Э. Карранза, У.Л. Шьогрен и Д.Н. Юань. «Недавние гравитационные модели в результате миссии лунного старателя», Икар, 150, № 1 (2001): 1-18.

[5] Lemoine, F. G., Д. Э. Смит, D.D. Роулендс, М. Т. Цубер, Г. А. Нейман и Д. С. Чинн. «Усовершенствованное решение гравитационного поля Марса (GMM-2B) от Mars Global Surveyor». Журнал геофизических исследований 106, np E10 (25 октября 2001): pp 23359-23376.

[6] Кеньон С., Дж. Фактор, Н. Павлис и С. Холмс. «К следующей гравитационной модели Земли.» Общество геофизиков-разведчиков 77-е ежегодное собрание, Сан-Антонио, Техас, 23-28 сентября 2007 года.

[7] Павлис, Н.К., С.А. Холмс, С.К. Кеньон и Дж. К. Фактор, «Гравитационная модель Земли до степени 2160: EGM2008.» Представлен на Генеральной ассамблее Европейского союза наук о Земле 2008 года, Вена, Австрия, 13-18 апреля 2008 года.

[8] Грубер, Т. и А. Кёль. «Проверка области силы тяжести EGM2008 с выравниванием GPS и океанографическими исследованиями». Представлен на Международном симпозиуме МАГ по гравитации, геоиду и наблюдению Земли 2008, Ханья, Греция, 23-27 июня 2008 года.

[9] Förste, C., Flechtner et al., «Средняя модель глобального гравитационного поля из комбинации данных о полетах спутников и альтиметрии/гравметрии поверхности - EIGEN-GL04C.» Тезисы геофизических исследований 8, 03462, 2006.

[10] Хилл, К. А. «Автономная навигация на орбитах точек либрации». Докторская диссертация, Колорадский университет, Боулдер. 2007.

[11] Коломбо, Оскар Л. «Численные методы гармонического анализа на сфере». Доклады Департамента геодезических наук, Доклад № 310, Университет штата Огайо, Колумбус, ок., март 1981 года.

[12] Коломбо, Оскар Л. «Глобальное картирование гравитации двумя спутниками». Нидерландская геодезическая комиссия 7, № 3, Делфт, Нидерланды, 1984., Отчеты отдела геодезических наук. Отчет № 310. Колумб: Университет штата Огайо, март 1981 года.

[13] Джонс, Брэндон А. «Эффективные модели оценки и оценки гравитационного поля». Докторская диссертация, Колорадский университет, Боулдер. 2010.

[14] Доклад Рабочей группы МАС/МАГ по картографическим координатам и элементам ротации: 1991 год.

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

Создание кода C/C + +
Создайте код C и C++ с помощью Simulink ® Coder™

.
Представлен в R2010a