gravitycentrifugal

Реализуйте центробежный эффект планетарной силы тяжести

Описание

Планетарный гравитационный потенциал на основе планетарного уровня вращения

пример

[gx gy gz] = gravitycentrifugal(planet_coordinates) реализует математическое представление центробежного эффекта для планетарной силы тяжести на основе планетарного уровня вращения. Эта функция вычисляет массивы N значения силы тяжести в x - ось, y - ось, и z - ось Сосредоточенных Планетой Зафиксированных Планетой координат для планеты. Функция выполняет эти вычисления с помощью planet_coordinates. Вы используете центробежную силу во вращении или неинерционных системах координат. Сила тяжести центробежные значения эффекта является самой большой на экватор планеты.

Планетарный гравитационный потенциал на основе заданной планетарной модели

пример

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

Планетарный гравитационный потенциал на основе пользовательского вращательного уровня

пример

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

Примеры

свернуть все

Вычислите центробежный эффект Наземной силы тяжести в оси X на экватор на поверхности Земли.

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

Вычислите центробежный эффект силы тяжести Марса на уровне 15 000 м по экватору и 11 000 м по Северному полюсу.

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

Вычислите прецессирующий центробежный эффект силы тяжести для Земли на уровне 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)
gx =
    0.0128
         0

gy =
   -0.0128
         0

gz =
     0
     0

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

свернуть все

Сосредоточенные планетой Зафиксированные Планетой координаты в виде M-by-3 массив в метрах. z - ось положительна к Северному полюсу. Если model 'Earth', координаты планеты являются координатами ECEF.

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

Планетарная модель в виде:

  • 'Mercury'

  • 'Venus'

  • 'Earth'

  • 'Moon'

  • 'Mars'

  • 'Jupiter'

  • 'Saturn'

  • 'Uranus'

  • 'Neptune'

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

Пользовательская планетарная модель в виде аргумента значения имени, где значение задает планетарный вращательный уровень в радианах в секунду.

Пример: 'Custom',Omega

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

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

свернуть все

Значения силы тяжести в x - ось Сосредоточенных Планетой Зафиксированных Планетой координат, возвращенных как массив значений силы тяжести M в метрах в секунду, придали квадратную форму (m/s2).

Значения силы тяжести в y - ось Сосредоточенных Планетой Зафиксированных Планетой координат, возвращенных как массив значений силы тяжести M в метрах в секунду, придали квадратную форму (m/s2).

Значения силы тяжести в z - ось Сосредоточенных Планетой Зафиксированных Планетой координат, возвращенных как массив значений силы тяжести M в метрах в секунду, придали квадратную форму (m/s2).

Смотрите также

|

Введен в R2010a