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_coordinates, an M-by-3 массив Планета-Центрированная Планета-Фиксированная Координата. Центробежная сила используется во вращающихся или неинерционных системах координат. Значения гравитационного центробежного эффекта наибольшие на экваторе планеты.

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

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

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

planet_coordinates

M-by-3 массив Планета-Центрированная Планета-Фиксированные координаты в метрах. The 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 )
Введенный в R2010a