gravitycentrifugal

Реализовать центробежный эффект планетарной гравитации

Синтаксис

[gx gy gz] = gravitycentrifugal(planet_coordinates) [gx gy gz] = gravitycentrifugal(planet_coordinates, model ) [gx gy gz] = gravitycentrifugal(planet_coordinates, 'Custom', rotational_rate)

Описание

[gx gy gz] = gravitycentrifugal(planet_coordinates) реализует математическое представление центробежного эффекта для планетарной гравитации на основе скорости планетарного вращения. Эта функция вычисляет массивы N значения силы тяжести по осям X, Y и Z в координатах Планетоцентрированная планета (Planet-Centered Planet) - фиксированные координаты для планеты. Эти вычисления выполняются с использованием planet_coordinatesодин M-на-3 массив координат Planet-Centered Planet-Fixed. Центробежная сила используется во вращающихся или неинерциальных системах координат. На экваторе планеты наибольшие значения гравитационного центробежного эффекта.

[gx gy gz] = gravitycentrifugal(planet_coordinates, model ) реализует математическое представление центробежного эффекта на основе планетарного гравитационного потенциала для планетарной модели, model.

[gx gy gz] = gravitycentrifugal(planet_coordinates, 'Custom', rotational_rate) реализует математическое представление центробежного эффекта на основе планетарного гравитационного потенциала с использованием пользовательской скорости вращения, rotational_rate.

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

planet_coordinates

M-на-3 массив Planet-Centered Planet - фиксированные координаты в метрах. z-аксис положителен к Северному полюсу. Если model является 'Earth'Координаты планеты - координаты ECEF.

model

Планетарная модель. По умолчанию: 'Earth'. Укажите один из следующих параметров:

'Mercury'
'Venus'
'Earth'
'Moon'
'Mars'
'Jupiter'
'Saturn'
'Uranus'
'Neptune'
'Custom'

'Custom' требуется указать собственную планетарную модель с помощью rotational_rate параметр.

rotational_rate

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

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

`gx`	Массив `M` значения гравитации по оси X Планеты-Центрированной Планеты - Фиксированные координаты в метрах в секунду в квадрате (^м/с2).
`gy`	Массив `M` значения гравитации по оси Y Планеты-Центрированной Планеты - Фиксированные координаты в метрах в секунду в квадрате (^м/с2).
`gz`	Массив `M` значения гравитации по оси Z Планеты-Центрированной Планеты - Фиксированные координаты в метрах в секунду в квадрате (^м/с2).

Примеры

Рассчитайте центробежный эффект гравитации Земли по оси X на экваторе на поверхности Земли:

gx = gravitycentrifugal( [-6378.1363e3 0 0] )

Рассчитать центробежный эффект гравитации Марса при 15000 м над экватором и 11000 м над Северным полюсом:

p  = [2412.648e3 -2412.648e3 0; 0 0 3376.2e3]
[gx, gy, gz] = gravitycentrifugal( p, 'Mars' )

Рассчитайте прецессирующий центробежный эффект гравитации для Земли на высоте 15000 м над экватором и 11000 м над Северным полюсом. В этом примере используется пользовательская планетарная модель в юлианскую дату 2451545:

p       = [2412.648e3 -2412.648e3 0; 0 0 3376e3]
% Set julian date to January 1, 2000 at noon GMT
JD      = 2451545
% Calculate precession rate in right ascension in meters
pres_RA = 7.086e-12 + 4.3e-15*(JD - 2451545)/36525
% Calculate the rotational rate in a precessing reference
% frame
Omega   = 7.2921151467e-5 + pres_RA
[gx, gy, gz] = gravitycentrifugal( p, 'custom', Omega )

См. также

gravitysphericalharmonic | gravitywgs84 | gravityzonal

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

Документация