Реализуйте центробежный эффект планетарной силы тяжести
[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, M- 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] )
Вычислите центробежный эффект силы тяжести Марса на уровне 15 000 м по экватору и 11 000 м по Северному полюсу:
p = [2412.648e3 -2412.648e3 0; 0 0 3376.2e3] [gx, gy, gz] = gravitycentrifugal( p, 'Mars' )
Вычислите прецессирующий центробежный эффект силы тяжести для Земли на уровне 15 000 м по экватору и 11 000 м по Северному полюсу. Этот пример использует пользовательскую планетарную модель в дату Джулиана 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 )