Реализуйте центробежный эффект планетарной гравитации
[
gx
gy
gz
]
= gravitycentrifugal(planet_coordinates
)
[gx
gy
gz
]
= gravitycentrifugal(planet_coordinates
, model
)
[gx
gy
gz
]
= gravitycentrifugal(planet_coordinates
,
'Custom', rotational_rate
)
[
реализует математическое представление центробежного эффекта для планетарной гравитации на основе скорости вращения планеты. Эта функция вычисляет массивы Ngx
gy
gz
]
= gravitycentrifugal(planet_coordinates
) значения гравитации по оси 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
.
|
|
|
Планетарная модель. По умолчанию это
|
|
Скалярное значение, которое задает скорость вращения планеты в радианах в секунду. Задайте этот параметр только, если |
|
Массив |
|
Массив |
|
Массив |
Вычислите центробежный эффект гравитации Земли по оси 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 )