Геопотенциальная модель

редактировать

В геофизике модель геопотенциала представляет собой теоретический анализ измерения и расчета эффекты Земли гравитационного поля.

Содержание
  • 1 Закон Ньютона
    • 1.1 Случай однородной сферы
  • 2 Отклонения гравитационного поля Земли от гравитационного поля однородная сфера
  • 3 Рекурсивные алгоритмы, используемые для численного распространения орбит космических аппаратов
    • 3.1 Доступные модели
  • 4 Сферические гармоники
  • 5 Ссылки
  • 6 Внешние ссылки
Закон Ньютона
Диаграмма двух массы, притягивающие друг друга

Закон всемирного тяготения Ньютона утверждает, что гравитационная сила F, действующая между двумя точечными массами m1и m 2 с центром масс разделение r определяется выражением

F = - G m 1 m 2 r 2 r ^ {\ displaystyle \ mathbf {F} = -G {\ frac {m_ {1} m_ {2}} {r ^ {2 }}} \ mathbf {\ hat {r}}}\mathbf {F} =-G{\frac {m_{1}m_{2}}{r^{2}}}\mathbf {\hat {r}}

где G - гравитационная постоянная и r̂ - радиальный единичный вектор. Для объекта с непрерывным распределением массы каждый элемент массы dm можно рассматривать как точечную массу, поэтому интеграл объема по протяженности объекта дает:

F ¯ = - G ∫ V ρ r 2 r ^ dxdydz {\ displaystyle \ mathbf {\ bar {F}} = -G \ int \ limits _ {V} {\ frac {\ rho} {r ^ {2}}} \ mathbf {\ hat {r} } \, dx \, dy \, dz}\mathbf {\bar {F}} =-G\int \limits _{V}{\frac {\rho }{r^{2}}}\mathbf {\hat {r}} \,dx\,dy\,dz

(1)

с соответствующим гравитационным потенциалом

u = - G ∫ V ρ rdxdydz {\ displaystyle u \ = -G \ int \ limits _ {V} {\ frac {\ rho} {r}} \, dx \, dy \, dz}{\displaystyle u\ =-G\int \limits _{V}{\frac {\rho }{r}}\,dx\,dy\,dz}

(2)

где ρ = ρ (x, y, z) - массовая плотность в элементе объема и направление от элемента объема к точечной массе.

Случай однородной сферы

В частном случае сферы со сферически симметричной плотностью массы, тогда ρ = ρ (s), т.е. плотность зависит только от радиального расстояния

s знак равно х 2 + у 2 + г 2. {\ displaystyle s = {\ sqrt {x ^ {2} + y ^ {2} + z ^ {2}}} \,.}s={\sqrt {x^{2}+y^{2}+z^{2}}}\,.

Эти интегралы можно вычислить аналитически. Это теорема оболочек, согласно которой в данном случае:

F ¯ = - GMR 2 r ^ {\ displaystyle {\ bar {F}} = - {\ frac {GM} {R ^ { 2}}} \ {\ hat {r}}}{\bar {F}}=-{\frac {GM}{R^{2}}}\ {\hat {r}}

(3)

с соответствующим потенциалом

u = - GM r {\ displaystyle u = - {\ frac {GM} {r }}}u=-{\frac {GM}{r}}

(4)

где M = ∫ V ρ (s) dxdydz - полная масса сферы.

Отклонения гравитационного поля Земли от поля однородной сферы

На самом деле Земля не совсем сферическая, в основном из-за ее вращения вокруг полярной оси, что делает ее форму слегка сжатой. Если бы эта форма была точно известна вместе с точной плотностью массы ρ = ρ (x, y, z), интегралы (1) и (2) можно было бы оценить численными методами, чтобы найти более точную модель для Земли. гравитационное поле. Однако на самом деле ситуация обратная. Наблюдая за орбитами космического корабля и Луны, гравитационное поле Земли может быть определено довольно точно, и наилучшая оценка массы Земли получается путем деления произведения GM, определенного в результате анализа орбиты космического корабля, на значение G, определенное на более низкую относительную величину. точность при использовании других физических методов.

Из определяющих уравнений (1) и (2) ясно (взяв частные производные от подынтегрального выражения), что вне тела в пустом пространстве следующие дифференциальные уравнения действительны для поля, вызванного по телу:

∂ F x ∂ x + ∂ F y ∂ y + ∂ F z ∂ z = 0 {\ displaystyle {\ frac {\ partial F_ {x}} {\ partial x}} + {\ frac {\ partial F_ {y}} {\ partial y}} + {\ frac {\ partial F_ {z}} {\ partial z}} = 0}{\frac {\partial F_{x}}{\partial x}}+{\frac {\partial F_{y}}{\partial y}}+{\frac {\partial F_{z}}{\partial z}}=0

(5)

∂ 2 u ∂ x 2 + ∂ 2 u ∂ Y 2 + ∂ 2 u ∂ Z 2 знак равно 0 {\ displaystyle {\ frac {\ partial ^ {2} u} {\ partial x ^ {2}}} + {\ frac {\ partial ^ { 2} u} {\ partial y ^ {2}}} + {\ frac {\ partial ^ {2} u} {\ partial z ^ {2}}} = 0}{\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}+{\frac {\partial ^{2}u}{\partial z^{2}}}=0

(6)

Функции формы ϕ = R (r) Θ (θ) Φ (φ) {\ Displaystyle \ phi \ = \ R (r) \ \ Theta (\ theta) \ \ Phi (\ varphi)}\phi \ =\ R(r)\ \Theta (\theta)\ \Phi (\varphi), где (r, θ, φ) - сферические координаты, которые удовлетворяют уравнению в частных производных (6) (уравнение Лапласа ), называются сферической гармоникой. функции.

Они принимают вид:

g (x, y, z) = 1 rn + 1 P nm (sin ⁡ θ) cos ⁡ m φ, 0 ≤ m ≤ n, n = 0, 1, 2,… h (x, y, z) = 1 rn + 1 P nm (sin ⁡ θ) sin ⁡ m φ, 1 ≤ m ≤ N, n = 1, 2,… {\ displaystyle {\ begin {align} g (x, y, z) = {\ frac {1} {r ^ {n + 1}}} P_ {n} ^ {m} (\ sin \ theta) \ cos m \ varphi \,, \ quad 0 \ leq m \ leq n \,, \ quad n = 0,1,2, \ dots \\ h (x, y, z) = {\ frac {1} {r ^ {n + 1}}} P_ {n} ^ {m} (\ sin \ theta) \ sin m \ varphi \,, \ quad 1 \ leq m \ leq n \,, \ quad n = 1,2, \ dots \ end {align}}}{\begin{aligned}g(x,y,z)={\frac {1}{r^{n+1}}}P_{n}^{m}(\sin \theta)\cos m\varphi \,,\quad 0\leq m\leq n\,,\quad n=0,1,2,\dots \\h(x,y,z)={\frac {1}{r^{n+1}}}P_{n}^{m}(\sin \theta)\sin m\varphi \,,\quad 1\leq m\leq n\,,\quad n=1,2,\dots \end{aligned}}

(7)

где сферические координаты (r, θ, φ), указанные здесь в декартовых (x, y, z) для справок:

x = r cos ⁡ θ cos ⁡ φ y = r cos ⁡ θ sin ⁡ φ z = r sin ⁡ θ, {\ Displaystyle {\ begin {выровнено} x = r \ cos \ theta \ cos \ varphi \\ y = r \ cos \ theta \ sin \ varphi \\ z = r \ sin \ theta \,, \ end {выравнивается} }}{\begin{aligned}x=r\cos \theta \cos \varphi \\y=r\cos \theta \sin \varphi \\z=r\sin \theta \,,\end{aligned}}

(8)

также P n - это многочлены Лежандра, а P n для 1 ≤ m ≤ n - связанные функции Лежандра.

Первые сферические гармоники с n = 0,1,2,3 представлены в таблице ниже.

nСферические гармоники
01 r {\ displaystyle {\ frac {1} {r}}}{\frac {1}{r}}
11 r 2 P 1 0 (sin ⁡ θ) = 1 r 2 sin ⁡ θ {\ displaystyle {\ frac {1} {r ^ {2}}} P_ {1} ^ {0} (\ sin \ theta) = {\ frac {1} {r ^ {2}}} \ sin \ theta}{\frac {1}{r^{2}}}P_{1}^{0}(\sin \theta)={\frac {1}{r^{2}}}\sin \theta
1 р 2 п 1 1 (грех ⁡ θ) соз ⁡ φ знак равно 1 р 2 соз ⁡ θ соз ⁡ φ {\ displaystyle {\ frac {1} {r ^ {2}}} P_ {1} ^ {1} (\ sin \ theta) \ cos \ varphi = {\ frac {1} {r ^ {2}}} \ cos \ theta \ cos \ varphi}{\frac {1}{r^{2}}}P_{1}^{1}(\sin \theta)\cos \varphi ={\frac {1}{r^{2}}}\cos \theta \cos \varphi
1 r 2 P 1 1 (sin ⁡ θ) sin ⁡ φ = 1 р 2 соз ⁡ θ грех ⁡ φ {\ displaystyle {\ frac {1} {r ^ {2}}} P_ {1} ^ {1} (\ sin \ theta) \ sin \ varphi = {\ frac {1 } {r ^ {2}}} \ cos \ theta \ sin \ varphi}{\frac {1}{r^{2}}}P_{1}^{1}(\sin \theta)\sin \varphi ={\frac {1}{r^{2}}}\cos \theta \sin \varphi
21 r 3 P 2 0 (sin ⁡ θ) = 1 r 3 1 2 (3 sin 2 ⁡ θ - 1) {\ displaystyle {\ frac {1} {r ^ {3}}} P_ {2} ^ {0} (\ sin \ theta) = {\ frac {1} {r ^ {3}}} {\ frac {1} { 2}} (3 \ sin ^ {2} \ theta -1)}{\ frac {1} {r ^ {3}}} P_ {2} ^ {0} (\ sin \ theta) = {\ frac {1} {r ^ {3}}} {\ frac {1} {2}} (3 \ sin ^ {2} \ theta -1)
1 r 3 P 2 1 (sin ⁡ θ) cos ⁡ φ = 1 r 3 3 sin ⁡ θ cos ⁡ θ cos ⁡ φ {\ displaystyle {\ frac {1} {r ^ {3}}} P_ {2} ^ {1} (\ sin \ theta) \ cos \ varphi = {\ frac {1} {r ^ {3}}} 3 \ sin \ theta \ cos \ theta \ \ cos \ varphi}{\frac {1}{r^{3}}}P_{2}^{1}(\sin \theta)\cos \varphi ={\frac {1}{r^{3}}}3\sin \theta \cos \theta \ \cos \varphi
1 r 3 P 2 1 (sin ⁡ θ) sin ⁡ φ = 1 r 3 3 sin ⁡ θ cos ⁡ θ sin ⁡ φ {\ di splaystyle {\ frac {1} {r ^ {3}}} P_ {2} ^ {1} (\ sin \ theta) \ sin \ varphi = {\ frac {1} {r ^ {3}}} 3 \ грех \ тета \ соз \ тета \ грех \ varphi}{\frac {1}{r^{3}}}P_{2}^{1}(\sin \theta)\sin \varphi ={\frac {1}{r^{3}}}3\sin \theta \cos \theta \sin \varphi
1 r 3 P 2 2 (грех ⁡ θ) соз ⁡ 2 φ = 1 r 3 3 cos 2 ⁡ θ cos ⁡ 2 φ {\ displaystyle {\ frac { 1} {r ^ {3}}} P_ {2} ^ {2} (\ sin \ theta) \ cos 2 \ varphi = {\ frac {1} {r ^ {3}}} 3 \ cos ^ {2 } \ theta \ \ cos 2 \ varphi}{\frac {1}{r^{3}}}P_{2}^{2}(\sin \theta)\cos 2\varphi ={\frac {1}{r^{3}}}3\cos ^{2}\theta \ \cos 2\varphi
1 r 3 P 2 2 (грех ⁡ θ) грех ⁡ 2 φ = 1 r 3 3 cos 2 ⁡ θ sin ⁡ 2 φ {\ displaystyle {\ frac {1} {r ^ {3}}} P_ {2} ^ {2} (\ sin \ theta) \ sin 2 \ varphi = {\ frac {1} {r ^ {3}}} 3 \ cos ^ {2} \ тета \ грех 2 \ varphi}{\frac {1}{r^{3}}}P_{2}^{2}(\sin \theta)\sin 2\varphi ={\frac {1}{r^{3}}}3\cos ^{2}\theta \sin 2\varphi
31 р 4 п 3 0 (грех ⁡ θ) = 1 р 4 1 2 грех ⁡ θ (5 грех 2 ⁡ θ - 3) {\ displaystyle {\ frac {1} {r ^ {4}}} P_ {3} ^ {0} (\ sin \ theta) = {\ frac {1} {r ^ {4}}} {\ frac {1} {2}} \ sin \ theta \ (5 \ sin ^ {2} \ theta -3)}{\frac {1}{r^{4}}}P_{3}^{0}(\sin \theta)={\frac {1}{r^{4}}}{\frac {1}{2}}\sin \theta \ (5\sin ^{2}\theta -3)
1 r 4 P 3 1 (sin ⁡ θ) cos ⁡ φ = 1 r 4 3 2 (5 sin 2 ⁡ θ - 1) cos ⁡ θ cos ⁡ φ {\ displaystyle {\ frac {1} {r ^ {4}}} P_ {3} ^ {1} (\ sin \ theta) \ cos \ varphi = {\ frac {1} {r ^ {4} }} {\ frac {3} {2}} \ (5 \ \ sin ^ {2} \ theta -1) \ cos \ theta \ cos \ varphi}{\frac {1}{r^{4}}}P_{3}^{1}(\sin \theta)\cos \varphi ={\frac {1}{r^{4}}}{\frac {3}{2}}\ (5\ \sin ^{2}\theta -1)\cos \theta \cos \varphi
1 r 4 P 3 1 (sin ⁡ θ) sin ⁡ φ = 1 r 4 3 2 (5 sin 2 ⁡ θ - 1) соз ⁡ θ грех ⁡ φ {\ displaystyle {\ frac {1} {r ^ {4}}} P_ {3} ^ {1} (\ sin \ theta) \ sin \ varphi = {\ frac {1} {r ^ {4}}} {\ frac {3} {2}} \ (5 \ \ sin ^ {2} \ theta -1) \ cos \ theta \ sin \ varphi}{\frac {1}{r^{4}}}P_{3}^{1}(\sin \theta)\sin \varphi ={\frac {1}{r^{4}}}{\frac {3}{2}}\ (5\ \sin ^{2}\theta -1)\cos \theta \sin \varphi
1 r 4 п 3 2 (грех ⁡ θ) соз ⁡ 2 φ = 1 р 4 15 грех ⁡ θ соз 2 ⁡ θ соз ⁡ 2 φ {\ displaystyle {\ frac {1} {r ^ {4}}} P_ {3} ^ {2} (\ sin \ theta) \ cos 2 \ varphi = {\ frac {1} {r ^ {4}}} 15 \ sin \ theta \ cos ^ {2} \ theta \ cos 2 \ varphi}{\frac {1}{r^{4}}}P_{3}^{2}(\sin \theta)\cos 2\varphi ={\frac {1}{r^{4}}}15\sin \theta \cos ^{2}\theta \cos 2\varphi
1 р 4 п 3 2 (грех ⁡ θ) грех ⁡ 2 φ = 1 r 4 15 грех ⁡ θ соз 2 ⁡ θ грех ⁡ 2 φ {\ displaystyle {\ frac {1} {r ^ {4}}} P_ {3} ^ {2} (\ sin \ theta) \ sin 2 \ varphi = {\ frac {1} {r ^ {4}}} 15 \ sin \ theta \ cos ^ {2} \ theta \ sin 2 \ varphi}{\ frac {1 }{r^{4}}}P_{3}^{2}(\sin \theta)\sin 2\varphi ={\frac {1}{r^{4}}}15\sin \theta \cos ^{2}\theta \sin 2\varphi
1 r 4 P 3 3 (грех ⁡ θ) cos ⁡ 3 φ = 1 r 4 15 cos 3 ⁡ θ cos ⁡ 3 φ {\ displaystyle {\ frac {1} {r ^ {4}} } P_ {3} ^ {3} (\ sin \ theta) \ cos 3 \ varphi = {\ frac {1} {r ^ {4}}} 15 \ cos ^ {3} \ theta \ cos 3 \ varphi}{\frac {1}{r^{4}}}P_{3}^{3}(\sin \theta)\cos 3\varphi ={\frac {1}{r^{4}}}15\cos ^{3}\theta \cos 3\varphi
1 р 4 P 3 3 (грех ⁡ θ) грех ⁡ 3 φ = 1 r 4 15 cos 3 ⁡ θ sin ⁡ 3 φ {\ displaystyle {\ frac {1} {r ^ {4}}} P_ { 3} ^ {3} (\ sin \ theta) \ sin 3 \ varphi = {\ frac {1} {r ^ {4}}} 15 \ cos ^ {3} \ theta \ sin 3 \ varphi}{\frac {1}{r^{4}}}P_{3}^{3}(\sin \theta)\sin 3\varphi ={\frac {1}{r^{4}}}15\cos ^{3}\theta \sin 3\varphi

Чт Модель гравитационного потенциала Земли представляет собой сумму

u = - μ r + ∑ n = 2 N z J n P n 0 (sin ⁡ θ) rn + 1 + ∑ n = 2 N t ∑ m = 1 n P нм (грех ⁡ θ) (С нм соз ⁡ м φ + S нм грех ⁡ м φ) rn + 1 {\ displaystyle u = - {\ frac {\ mu} {r}} + \ sum _ {n = 2} ^ {N_ {z}} {\ frac {J_ {n} P_ {n} ^ {0} (\ sin \ theta)} {r ^ {n + 1}}} + \ sum _ {n = 2} ^ {N_ {t}} \ sum _ {m = 1} ^ {n} {\ frac {P_ {n} ^ {m} (\ sin \ theta) (C_ {n} ^ {m} \ cos m \ varphi + S_ {n} ^ {m} \ sin m \ varphi)} {r ^ {n + 1}}}}u=-{\frac {\mu }{r}}+\sum _{n=2}^{N_{z}}{\frac {J_{n}P_{n}^{0}(\sin \theta)}{r^{n+1}}}+\sum _{n=2}^{N_{t}}\sum _{ m=1}^{n}{\frac {P_{n}^{m}(\sin \theta)(C_{n}^{m}\cos m\varphi +S_{n}^{m}\sin m\varphi)}{r^{n+1}}}

(9)

где μ = GM {\ displaystyle \ mu = GM }\mu =GM, а координаты (8) являются относительно стандартной геодезической системы отсчета расширен в пространстве с началом в центре эллипсоида и оси в направлении полярная ось.

зональные термины относятся к терминам формы:

P n 0 (sin ⁡ θ) rn + 1 n = 0, 1, 2,… {\ displaystyle {\ frac {P_ {n} ^ {0} (\ sin \ theta)} {r ^ {n + 1}}} \ quad n = 0,1,2, \ dots}{\frac {P_{n}^{0}(\sin \theta)}{r^{n+1}}}\quad n=0,1,2,\dots

и тессеральные термины термины относятся к терминам формы:

P nm (sin ⁡ θ) cos ⁡ m φ rn + 1, 1 ≤ m ≤ nn = 1, 2,… {\ displaystyle {\ frac {P_ {n) } ^ {m} (\ sin \ theta) \ cos m \ varphi} {r ^ {n + 1}}} \,, \ quad 1 \ leq m \ leq n \ quad n = 1,2, \ dots}{\frac {P_{n}^{m}(\sin \theta)\cos m\varphi }{r^{n+1}}}\,,\quad 1\leq m\leq n\quad n=1,2,\dots
P нм (грех ⁡ θ) грех ⁡ м φ rn + 1 {\ displaystyle {\ frac {P_ {n} ^ {m} (\ sin \ theta) \ sin m \ varphi} {r ^ {n + 1}}}}{\frac {P_{n}^{m}(\sin \theta)\sin m\varphi }{r^{n+1}}}

Зональные и тессеральные термины для n = 1 опущены в (9). Коэффициенты при n = 1 с членами как m = 0, так и m = 1 соответствуют произвольно ориентированному дипольному члену в многополюсном разложении. Гравитация физически не проявляет дипольного характера, поэтому интеграл, характеризующий n = 1, должен быть равен нулю.

различным коэффициентам J n, C n, S n затем присваиваются значения, для которых наилучшее возможное совпадение между вычисленными и получены наблюдаемые орбиты космических аппаратов.

As P n (x) = −P n (−x) ненулевые коэффициенты J n для нечетного n соответствуют отсутствию симметрии «север-юг» относительно экваториальной плоскости для распределения массы Земли. Ненулевые коэффициенты C n, S n соответствуют отсутствию вращательной симметрии вокруг полярной оси для распределения массы Земли, то есть «трехосности» Земли.

Для больших значений n приведенные выше коэффициенты (которые делятся на r в (9)) принимают очень большие значения, когда, например, в качестве единиц используются километры и секунды. В литературе принято вводить произвольный «опорный радиус» R, близкий к радиусу Земли, и работать с безразмерными коэффициентами

J n ~ = - J n μ R n {\ displaystyle {\ tilde {J_ {n}) }} = - {\ frac {J_ {n}} {\ mu \ R ^ {n}}}}{\tilde {J_{n}}}=-{\frac {J_{n}}{\mu \ R^{n}}}
C нм ~ = - C нм μ R n {\ displaystyle {\ tilde {C_ {n} ^ {m}}} = - {\ frac {C_ {n} ^ {m}} {\ mu \ R ^ {n}}}}{\tilde {C_{n}^{m}}}=-{\frac {C_{n}^{m}}{\mu \ R^{n}}}
S нм ~ = - S нм μ R n {\ displaystyle {\ тильда {S_ {n} ^ {m}}} = - {\ frac {S_ {n} ^ {m}} {\ mu \ R ^ {n}}}}{\tilde {S_{n}^{m}}}=-{\frac {S_{n}^{m}}{\mu \ R^{n}}}

и записать потенциал как

u = - μ r (1 + ∑ n = 2 N z J n ~ P n 0 (sin ⁡ θ) (r R) n + ∑ n = 2 N t ∑ m = 1 n P nm (sin ⁡ θ) (C нм ~ соз ⁡ м φ + S нм ~ грех ⁡ м φ) (r R) n) {\ displaystyle u = - {\ frac {\ mu} {r}} \ left (1+ \ sum _ {n = 2} ^ {N_ {z}} {\ frac {{\ tilde {J_ {n}}} P_ {n} ^ {0} (\ sin \ theta)} {{({\ frac {r} {R }})} ^ {n}}} + \ sum _ {n = 2} ^ {N_ {t}} \ sum _ {m = 1} ^ {n} {\ frac {P_ {n} ^ {m} (\ sin \ theta) ({\ tilde {C_ {n} ^ {m}}} \ cos m \ varphi + {\ tilde {S_ {n} ^ {m}}} \ sin m \ varphi)} {{ ({\ frac {r} {R}})} ^ {n}}} \ right)}u=-{\frac {\mu }{r}}\left(1+\sum _{n=2}^{N_{z}}{\frac {{\tilde {J_{n}}}P_{n}^{0}(\sin \theta)}{{({\frac {r}{R}})}^{n}}}+\sum _{n=2}^{N_{t}}\sum _{m=1}^{n}{\frac {P_{n}^{m}(\sin \theta)({\tilde {C_{n}^{m}}}\cos m\varphi +{\tilde {S_{n}^{m}}}\sin m\varphi)}{{({\frac {r}{R} })}^{n}}}\right)

(10)

Главный член (после члена −μ / r) in (9) - это «член J 2 »:

u = J 2 P 2 0 (sin ⁡ θ) r 3 = J 2 1 r 3 1 2 (3 sin 2 ⁡ θ - 1) знак равно J 2 1 р 5 1 2 (3 z 2 - r 2) {\ displaystyle u = {\ frac {J_ {2} \ P_ {2} ^ {0} (\ sin \ theta)} { r ^ {3}}} = J_ {2} {\ frac {1} {r ^ {3}}} {\ frac {1} {2}} (3 \ sin ^ {2} \ theta -1) = J_ {2} {\ frac {1} {r ^ {5}}} {\ frac {1} {2}} (3z ^ {2} -r ^ {2})}u={\frac {J_{2}\ P_{2}^{0}(\sin \theta)}{r^{3}}}=J_{2}{\frac {1}{r^{3}}}{\frac {1}{2}}(3\sin ^{2}\theta -1)=J_{2}{\frac {1}{r^{5}}}{\frac {1}{2}}(3z^{2}-r^{2})

Относительно системы координат

φ ^ = - sin ⁡ φ x ^ + cos ⁡ φ y ^ θ ^ = - sin ⁡ θ (cos ⁡ φ x ^ + sin ⁡ φ y ^) + cos ⁡ θ z ^ r ^ = cos ⁡ θ ( соз ⁡ φ Икс ^ + грех ⁡ φ Y ^) + грех ⁡ θ Z ^ {\ Displaystyle {\ begin {align} {\ hat {\ varphi}} = - \ sin \ varphi {\ hat {x}} + \ cos \ varphi {\ hat {y}} \\ {\ hat {\ theta}} = - \ sin \ theta \ (\ cos \ varphi \ {\ hat {x}} + \ sin \ varphi {\ hat {y}}) + \ cos \ theta {\ hat {z}} \\ {\ hat {r}} = \ cos \ theta \ (\ cos \ varphi \ {\ hat {x}} \ + \ \ sin \ varphi \ {\ hat {y}}) + \ \ sin \ theta \ {\ hat {z}} \ end {align}}}{\begin{aligned}{\hat {\varphi }}=-\sin \varphi {\hat {x}}+\cos \varphi {\hat {y}}\\{\hat {\theta }}=-\sin \theta \ (\cos \varphi \ {\hat {x}}+\sin \varphi {\hat {y}})+\cos \theta {\hat {z}}\\{\hat {r}}=\cos \theta \ (\cos \varphi \ {\hat {x}}\ +\ \sin \varphi \ {\hat {y}})+\ \sin \theta \ {\hat {z}}\end{aligned}}

(11)

Рисунок 1: Единичные векторы. Это не правильно. Должна быть тета, а не лямбда φ ^, θ ^, r ^ {\ displaystyle {\ hat {\ varphi}} \, \ {\ hat {\ theta}} \, \ {\ hat {r} }}{\hat {\varphi }}\,\ {\hat {\theta }}\,\ {\hat {r}}

проиллюстрированные на рисунке 1 компоненты силы, вызванной «членом J 2 », равны

F θ = - 1 r ∂ u ∂ θ = - J 2 1 r 4 3 соз ⁡ θ грех ⁡ θ F r = - ∂ u ∂ r = J 2 1 r 4 3 2 (3 sin 2 ⁡ θ - 1) {\ displaystyle {\ begin {align} F _ {\ theta} = - {\ frac {1} {r}} \ {\ frac {\ partial u} {\ partial \ theta}} = - J_ {2} \ {\ frac {1} {r ^ {4}}} 3 \ \ cos \ theta \ \ sin \ theta \\ F_ {r} = - {\ frac {\ partial u} {\ partial r}} = J_ {2} \ {\ frac {1} {r ^ {4}}} {\ frac {3} {2}} \ left (3 \ sin ^ {2} \ theta \ - \ 1 \ right) \ end {align}}}{\begin{aligned}F_{\theta }=-{\frac {1}{r}}\ {\frac {\partial u}{\partial \theta }}=-J_{2}\ {\frac {1}{r^{4}}}3\ \cos \theta \ \sin \theta \\F_{r}=-{\frac {\partial u}{\partial r}}=J_{2}\ {\frac {1}{r^{4}}}{\frac {3}{2}}\ \left(3\sin ^{2}\theta \ -\ 1\right)\end{aligned}}

(12)

В прямоугольной системе координат (x, y, z) с единичными векторами (x̂ ŷ ẑ) компоненты силы:

F x = - ∂ u ∂ x = J 2 xr 7 (6 z 2 - 3 2 (x 2 + y 2)) F y = - ∂ u ∂ y = J 2 года 7 (6 z 2 - 3 2 (x 2 + y 2)) F z = - ∂ u ∂ z = J 2 zr 7 (3 z 2 - 9 2 (x 2 + y 2)) {\ displaystyle {\ begin {align} F_ {x} = - {\ frac {\ partial u} {\ partial x}} = J_ {2} \ {\ frac { x} {r ^ {7}}} \ left (6z ^ {2} - {\ frac {3} {2}} (x ^ {2} + y ^ {2}) \ right) \\ F_ {y } = - {\ frac {\ partial u} {\ partial y}} = J_ {2} \ {\ frac {y} {r ^ {7}}} \ left (6z ^ {2} - {\ frac { 3} {2}} (x ^ {2} + y ^ {2}) \ right) \\ F_ {z} = - {\ frac {\ partial u} {\ partial z}} = J_ {2} \ {\ frac {z} {r ^ {7}}} \ left (3z ^ {2} - {\ frac {9} {2}} (x ^ {2} + y ^ {2}) \ right) \ конец {выровнен}}}{\begin{aligned}F_{x}=-{\frac {\partial u}{\partial x}}=J_{2}\ {\frac {x}{r^{7}}}\left(6z^{2}-{\frac {3}{2}}(x^{2}+y^{2})\right)\\F_{y}=-{\frac {\partial u}{\partial y}}=J_{2}\ {\frac {y}{r^{7}}}\left(6z^{2}-{\frac {3}{2}}(x^{2}+y^{2})\right)\\F_{z}=-{\frac {\partial u}{\partial z}}=J_{2}\ { \frac {z}{r^{7}}}\left(3z^{2}-{\frac {9}{2}}(x^{2}+y^{2})\right)\end{aligned}}

(13)

Составляющие силы, соответствующие «члену J 3 »

u = J 3 P 3 0 (sin ⁡ θ) r 4 знак равно J 3 1 р 4 1 2 грех ⁡ θ (5 грех 2 ⁡ θ - 3) = J 3 1 r 7 1 2 z (5 z 2 - 3 r 2) {\ displaystyle u = {\ frac {J_ { 3} P_ {3} ^ {0} (\ sin \ theta)} {r ^ {4}}} = J_ {3} {\ frac {1} {r ^ {4}}} {\ frac {1} {2}} \ sin \ theta (5 \ sin ^ {2} \ theta -3) = J_ {3} {\ frac {1} {r ^ {7}}} {\ frac {1} {2}} z (5z ^ {2} -3r ^ {2})}u={\frac {J_{3}P_{3}^{0}(\sin \theta)}{r^{4}}}=J_{3}{\frac {1}{r^{4}}}{\frac {1}{2}}\sin \theta (5\sin ^{2}\theta -3)=J_{3}{\frac {1}{r^{7}}}{\frac {1}{2}}z(5z^{2}-3r^{2})

равны

F θ = - 1 r ∂ u ∂ θ = - J 3 1 r 5 3 2 cos ⁡ θ (5 sin 2 ⁡ θ - 1) F р знак равно - ∂ U ∂ р знак равно J 3 1 r 5 2 грех ⁡ θ (5 sin 2 ⁡ θ - 3) {\ displaystyle {\ begin {align} F _ {\ theta} = - {\ frac { 1} {r}} {\ frac {\ partial u} {\ partial \ theta}} = - J_ {3} {\ frac {1} {r ^ {5}}} {\ frac {3} {2} } \ cos \ theta \ left (5 \ sin ^ {2} \ theta -1 \ right) \\ F_ {r} = - {\ frac {\ partial u} {\ partial r}} = J_ {3} {\ frac {1} {r ^ {5}}} 2 \ sin \ theta \ left (5 \ sin ^ {2} \ theta -3 \ right) \ end {align}}}{\begin{aligned}F_{\theta }=-{\frac {1}{r}}{\frac {\partial u}{\partial \theta }}=-J_{3}{\frac {1}{r^{5}}}{\frac {3}{2}}\cos \theta \left(5\sin ^{2}\theta -1\right)\\F_{r}=-{\frac {\partial u}{\partial r}}=J_{3}{\frac {1}{r^{5}}}2\sin \theta \left(5\sin ^{2}\theta -3\right)\end{aligned}}

(14)

и

F x = - ∂ u ∂ x = J 3 xzr 9 ( 10 z 2 - 15 2 (x 2 + y 2)) F y = - ∂ u ∂ y = J 3 yzr 9 (10 z 2 - 15 2 (x 2 + y 2)) F z = - ∂ u ∂ z Знак равно J 3 1 р 9 (4 z 2 (z 2 - 3 (x 2 + y 2)) + 3 2 (x 2 + y 2) 2) {\ displaystyle {\ begin {align} F_ {x} = - {\ frac {\ partial u} {\ partial x}} = J_ {3} {\ frac {xz} {r ^ {9}}} \ left (10z ^ {2} - {\ frac {15} {2 }} (x ^ {2} + y ^ {2}) \ right) \\ F_ {y} = - {\ frac {\ partial u} {\ partial y}} = J_ {3} {\ frac {yz } {r ^ {9}}} \ left (10z ^ {2} - {\ frac {15} {2}} (x ^ {2} + y ^ {2}) \ right) \\ F_ {z} = - {\ frac {\ partial u} {\ partial z}} = J_ {3} {\ frac {1} {r ^ {9}}} \ left (4z ^ {2} \ \ left (z ^ { 2} -3 (x ^ {2} + y ^ {2}) \ right) + {\ frac {3} {2}} (x ^ {2} + y ^ {2}) ^ {2} \ right) \ end {align}}}{\begin{aligned}F_{x}=-{\frac {\partial u}{\partial x}}=J_{3}{\frac {xz}{r^{9}}}\left(10z^{2}-{\frac {15}{2}}(x^{2}+y^{2})\right)\\F_{y}=-{\frac {\partial u}{\partial y}}=J_{3}{\frac {yz}{r^{9}}}\left(10z^{2}-{\frac {15}{2}}(x^{2}+y^{2})\right)\\F_{z}=-{\frac {\partial u}{\partial z}}=J_{3}{\frac {1}{r^{9}}}\left(4z^{2}\ \left(z^{2 }-3(x^{2}+y^{2})\right)+{\frac {3}{2}}(x^{2}+y^{2})^{2}\right)\end{aligned}}

(15)

Точные численные значения коэффициентов различаются (несколько) между разными моделями Земли, но для самых низких коэффициентов все они почти полностью совпадают.

Для JGM-3 значения следующие:

μ = 398600,440 км⋅с
J2= 1,75553 × 10 км⋅с
J3= −2,61913 × 10 км⋅с

. Например, на радиусе 6600 км (около 200 км над поверхностью Земли) J 3 / (J 2 r) составляет около 0,002, то есть поправка к «J 2 сила »из« члена J 3 »составляет порядка 2 промилле. Отрицательное значение J 3 подразумевает, что для точечной массы в экваториальной плоскости Земли гравитационная сила слегка наклонена к югу из-за отсутствия симметрии массового распределения Земли «север-юг».

Рекурсивные алгоритмы, используемые для численного распространения орбит космических аппаратов

Орбиты космических аппаратов вычисляются посредством численного интегрирования уравнения движения. Для этого необходимо вычислить гравитационную силу, т.е. градиент потенциала. Эффективные рекурсивные алгоритмы были разработаны для вычисления силы тяжести для любых N z {\ displaystyle N_ {z}}N_{z}и N t {\ displaystyle N_ {t }}N_{t}(максимальная степень зональных и тессеральных условий) и такие алгоритмы используются в стандартном программном обеспечении для распространения орбиты.

Доступные модели

Самыми ранними моделями Земли, обычно используемыми NASA и ESRO / ESA, были «Земля Годдарда. Модели », разработанные Центром космических полетов им. Годдарда, обозначались« GEM-1 »,« GEM-2 »,« GEM-3 »и так далее. Позже стали доступны «Совместные модели земной гравитации», обозначенные «JGM-1», «JGM-2», «JGM-3», разработанные Центром космических полетов Годдарда в сотрудничестве с университетами и частными компаниями. В более новых моделях обычно использовались термины более высокого порядка, чем в их предшественниках. EGM96 использует N z = N t = 360, что дает 130317 коэффициентов. Также доступна модель EGM2008.

Для обычного спутника Земли, требующего точности определения / предсказания орбиты в несколько метров, "JGM-3" усечено до N z = N t = 36 ( 1365 коэффициентов) обычно достаточно. Неточности моделирования сопротивления воздуха и, в меньшей степени, давления солнечной радиации будут превышать неточности, вызванные ошибками моделирования гравитации.

Безразмерные коэффициенты J n ~ = - J n μ R n {\ displaystyle {\ tilde {J_ {n}}} = - {\ frac {J_ {n}} {\ mu \ R ^ {n}}}}{\tilde {J_{n}}}=-{\frac {J_{n}}{\mu \ R^{n}}}, C нм ~ = - C нм μ R n {\ displaystyle {\ tilde {C_ {n} ^ {m}}} = - {\ frac {C_ {n} ^ { m}} {\ mu \ R ^ {n}}}}{\tilde {C_{n}^{m}}}=-{\frac {C_{n}^{m}}{\mu \ R^{n}}}, S нм ~ = - S нм μ R n {\ displaystyle {\ tilde {S_ {n} ^ {m}}} = - {\ frac {S_ {n} ^ {m}} {\ mu \ R ^ {n}}}}{\tilde {S_{n}^{m}}}=-{\frac {S_{n}^{m}}{\mu \ R^{n}}}для первых зональных и тессеральных терминов (с использованием R {\ displaystyle R}R= 6378,1363 км и μ {\ displaystyle \ mu}\mu = 398600,4415 км / с) модели JGM-3 равны

Зональные коэффициенты
n
2-0,1082635854D-02
30,2532435346D-05
40,1619331205D-05
50,2277161016D-06
6-0,5396484906D-06
70,3513684422D-06
80.2025187152D-06
Тессеральные коэффициенты
nmCS
21- 0.3504890360D-090.1635406077D-08
220.1574536043D-05-0.9038680729D-06
310.2192798802D-050.2680118938D-06
320,3090160446D-06-0,2114023978D-06
330.1005588574D-060,1972013239D-06
41-0,5087253036D-06- 0.4494599352D-06
420.7841223074D-070.1481554569D-06
430.5921574319D-07-0.1201129183D-07
44-0.3982395740D-080,6525605810D-08

Следовательно, согласно JGM-3, J 2 = 0,1082635854 ⋅ 10–02 ⋅ 6378,1363 2 ⋅ 398600.4415 {\ displaystyle J_ {2} \ = \ 0.1082635854 \ \ cdot \ 10 ^ {- 02} \ \ cdot \ 6378.1363 ^ {2} \ \ cdot \ 398600.4415}{\displaystyle J_{2}\ =\ 0.1082635854\ \cdot \ 10^{-02}\ \cdot \ 6378.1363^{2}\ \cdot \ 398600.4415}км / с = 1,75553 ⋅ 10 10 {\ displaystyle 1.75553 \ \ cdot \ 10 ^ {10 }}{\displaystyle 1.75553\ \cdot \ 10^{10}}км / с и J 3 = - 0,2532435346 ⋅ 10 - 05 ⋅ 6378.1363 3 ⋅ 398600.4415 {\ displaystyle J_ {3} \ = \ -0.2532435346 \ \ cdot \ 10 ^ {- 05 } \ \ cdot \ 6378.1363 ^ {3} \ \ cdot \ 398600.4415}{\displaystyle J_ {3} \ = \ -0.2532435346 \ \ cdot \ 10 ^ {- 05} \ \ cdot \ 6378.1363 ^ {3} \ \ cdot \ 398600.4415} км / с = - 2.61913 ⋅ 10 11 {\ displaystyle -2.61913 \ \ cdot \ 10 ^ {11}}{\displaystyle -2.61913\ \cdot \ 10^{11}}км / с

Сферические гармоники

Ниже приводится краткое описание сферических гармоник, используемых для моделирования гравитационного поля Земли. Сферические гармоники являются производными от поиска гармонических функций вида

ϕ = R (r) Θ (θ) Φ (φ) {\ displaystyle \ phi \ = \ R (r) \ \ Theta (\ theta) \ \ Phi (\ varphi)}\phi \ =\ R(r)\ \Theta (\theta)\ \Phi (\varphi)

(16)

где (r, θ, φ) - это сферические координаты, определенные уравнениями (8). Непосредственными вычислениями получаем, что для любой функции f

∂ 2 f ∂ x 2 + ∂ 2 f ∂ y 2 + ∂ 2 f ∂ z 2 = 1 r 2 ∂ ∂ r (r 2 ∂ f ∂ r) + 1 р 2 соз ⁡ θ ∂ ∂ θ (соз ⁡ θ ∂ f ∂ θ) + 1 r 2 соз 2 ⁡ θ ∂ 2 f ∂ φ 2 {\ displaystyle {\ frac {\ partial ^ {2} f} {\ partial x ^ {2}}} \ + \ {\ frac {\ partial ^ {2} f} {\ partial y ^ {2}}} \ + \ {\ frac {\ partial ^ {2} f} {\ partial z ^ {2}}} \ = \ {1 \ над r ^ {2}} {\ partial \ over \ partial r} \ left (r ^ {2} {\ partial f \ over \ partial r} \ right) + {1 \ over r ^ {2} \ cos \ theta} {\ partial \ over \ partial \ theta} \ left (\ cos \ theta {\ partial f \ over \ partial \ theta} \ right) + {1 \ over r ^ {2} \ cos ^ {2} \ theta} {\ partial ^ {2} f \ over \ partial \ varphi ^ {2}}}{\frac {\partial ^{2}f}{\partial x^{2}}}\ +\ {\frac {\partial ^{2}f}{\partial y^{2}}}\ +\ {\frac {\partial ^{2}f}{\partial z^{2}}}\ =\ {1 \over r^{2}}{\partial \over \partial r}\left(r^{2}{\partial f \over \partial r}\right)+{1 \over r^{2}\cos \theta }{\partial \over \partial \theta }\left(\cos \theta {\partial f \over \partial \theta }\right)+{1 \over r^{2}\cos ^{2}\theta }{\partial ^{2}f \over \partial \varphi ^{2}}

(17)

Представляем выражение (16) в (17) получаем, что

r 2 ϕ (∂ 2 ϕ ∂ x 2 + ∂ 2 ϕ ∂ y 2 + ∂ 2 ϕ ∂ z 2) = 1 R ddr (r 2 d R dr) + 1 Θ соз ⁡ θ dd θ (соз ⁡ θ d Θ d θ) + 1 Φ cos 2 ⁡ θ d 2 Φ d φ 2 {\ displaystyle {\ frac {r ^ {2}} {\ phi}} \ left ( {\ frac {\ partial ^ {2} \ phi} {\ partial x ^ {2}}} + {\ frac {\ partial ^ {2} \ phi} {\ partial y ^ {2}}} + {\ frac {\ partial ^ {2} \ phi} {\ partial z ^ {2}}} \ right) \ = {\ frac {1} {R}} {\ frac { d} {dr}} \ left (r ^ {2} {\ frac {dR} {dr}} \ right) + {\ frac {1} {\ Theta \ cos \ theta}} {\ frac {d} { d \ theta}} \ left (\ cos \ theta {\ frac {d \ Theta} {d \ theta}} \ right) + {\ frac {1} {\ Phi \ cos ^ {2} \ theta}} { \ frac {d ^ {2} \ Phi} {d \ varphi ^ {2}}}}{\frac {r^{2}}{\phi }}\left({\frac {\partial ^{2}\phi }{\partial x^{2}}}+{\frac {\partial ^{2}\phi }{\partial y^{2}}}+{\frac {\partial ^{2}\phi }{\partial z^{2}}}\right)\ ={\frac {1}{R}}{\frac {d}{dr}}\left(r^{2}{\frac {dR}{dr}}\right)+{\frac {1}{\Theta \cos \theta }}{\frac {d}{d\theta }}\left(\cos \theta {\frac {d\Theta }{d\theta }}\right)+{\frac {1}{\Phi \cos ^{2}\theta }}{\frac {d^{2}\Phi }{d\varphi ^{2}}}

(18)

В качестве члена

1 R ddr (r 2 d R dr) {\ displaystyle {\ frac {1} {R}} {\ frac {d} {dr}} \ left (r ^ {2} {\ frac {dR} {dr}} \ right)}{\frac {1}{R}}{\frac {d}{dr}}\left(r^{2}{\frac {dR}{dr}}\right)

зависит только от переменной r {\ displaystyle r}rи сумма

1 Θ cos ⁡ θ dd θ (cos ⁡ θ d Θ d θ) + 1 Φ cos 2 ⁡ θ d 2 Φ d φ 2 {\ displaystyle {\ frac {1} {\ Theta \ cos \ theta}} {\ frac {d} {d \ theta}} \ left (\ cos \ theta {\ frac {d \ Theta} {d \ theta}) } \ right) + {\ frac {1} {\ Phi \ cos ^ {2} \ theta}} {\ frac {d ^ {2} \ Phi} {d \ varphi ^ {2}}}}{\frac {1}{\Theta \cos \theta }}{\frac {d}{d\theta }}\left(\cos \theta {\frac {d\Theta }{d\theta }}\right)+{\frac {1}{\Phi \cos ^{2}\theta }}{\frac {d^{2}\Phi }{d\varphi ^{2}}}

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

1 R ddr (r 2 d R dr) = λ {\ displaystyle {\ frac {1} {R}} {\ frac {d} {dr}} \ left (r ^ {2} {\ frac {dR} {dr}} \ right) \ = \ \ lambda}{\frac {1}{R}}{\frac {d}{dr}}\left(r^{2}{\frac {dR}{dr}}\right)\ =\ \lambda

(19)

и

1 Θ cos ⁡ θ dd θ (cos ⁡ θ d Θ d θ) + 1 Φ cos 2 ⁡ θ d 2 Φ d φ 2 = - λ {\ displaystyle {\ frac {1} {\ Theta \ cos \ theta}} {\ frac {d} {d \ theta}} \ left (\ cos \ theta {\ frac {d \ Theta} {d \ theta}} \ right) + {\ frac {1} {\ Phi \ cos ^ {2} \ theta}} {\ frac {d ^ {2} \ Phi} {d \ varphi ^ {2}}} \ = \ - \ lambda}{\frac {1}{\Theta \cos \theta }}{\frac {d}{d\theta }}\left(\cos \theta {\frac {d\Theta }{d\theta }}\right)+{\frac {1}{\Phi \cos ^{2}\theta }}{\frac {d^{2}\Phi }{d\varphi ^{2}}}\ =\ -\lambda

(20)

для некоторой константы λ {\ displaystyle \ lambda}\lambda

From ( 20) тогда следует, что

1 Θ соз ⁡ θ dd θ (соз ⁡ θ d Θ d θ) + λ cos 2 ⁡ θ + 1 Φ d 2 Φ d φ 2 = 0 {\ displaystyle {\ frac { 1} {\ Theta}} \ \ cos \ theta \ {\ frac {d} {d \ theta}} \ left (\ cos \ theta {\ frac {d \ Theta} {d \ theta}} \ right) \ + \ lambda \ \ cos ^ {2} \ theta \ + \ {\ frac {1} {\ Phi}} {\ frac {d ^ {2} \ Phi} {d \ varphi ^ {2}}} \ = \ 0}{\frac {1}{\Theta }}\ \cos \theta \ {\frac {d}{d\theta }}\left(\cos \theta {\frac {d\Theta }{d\theta }}\right)\ +\lambda \ \cos ^{2}\theta \ +\ {\frac {1}{\Phi }}{\frac {d^{2}\Phi }{d\varphi ^{2}}}\ =\ 0

Первые два члена зависят только от переменной θ {\ displaystyle \ theta}\theta , а третий - только от переменной φ {\ displaystyle \ varphi}\varphi .

Из определения φ как сферической координаты ясно, что Φ (φ) должна быть периодической с периодом 2π, и поэтому должно быть, что

1 Φ d 2 Φ d φ 2 знак равно - м 2 {\ displaystyle {\ frac {1} {\ Phi}} {\ frac {d ^ {2} \ Phi} {d \ varphi ^ {2}}} \ = \ -m ^ { 2}}{\frac {1}{\Phi }}{\frac {d^{2}\Phi }{d\varphi ^{2}}}\ =\ -m^{2}

(21)

и

1 Θ cos ⁡ θ dd θ (cos ⁡ θ d Θ d θ) + λ cos 2 ⁡ θ = m 2 {\ displaystyle {\ frac {1} {\ Theta}} \ \ cos \ theta \ {\ frac {d} {d \ theta}} \ left (\ cos \ theta {\ frac {d \ Theta} {d \ theta}} \ right) \ + \ лямбда \ \ cos ^ {2} \ theta = \ m ^ {2}}{\frac {1}{\Theta }}\ \cos \theta \ {\frac {d}{d\theta }}\left(\cos \theta {\frac {d\Theta }{d\theta }}\right)\ +\lambda \ \cos ^{2}\theta =\ m^{2}

(22)

для некоторого целого числа m в качестве семейства решений (21), тогда равны

Φ (φ) знак равно a соз ⁡ м φ + б грех ⁡ м φ {\ displaystyle \ Phi (\ varphi) \ = a \ \ cos m \ varphi \ + \ b \ \ sin m \ varphi}\Phi (\varphi)\ =a\ \cos m\varphi \ +\ b\ \sin m\varphi

(23)

При подстановке переменной

x = sin ⁡ θ {\ displaystyle x = \ sin \ theta}x=\sin \theta

уравнение (22) принимает форму

ddx ((1 - x 2) d Θ dx) + (λ - м 2 1 - x 2) Θ знак равно 0 {\ displaystyle {\ frac {d} {dx}} \ left ((1-x ^ {2}) {\ frac {d \ Theta} {dx} } \ right) + \ left (\ lambda - {\ frac {m ^ {2}} {1-x ^ {2}}} \ right) \ Theta = 0}{\frac {d}{dx}}\left((1-x^{2}){\frac {d\Theta }{dx}}\right)+\left(\lambda -{\frac {m^{2}}{1-x^{2}}}\right)\Theta =0

(24)

Из (19) следует, что в чтобы получить решение ϕ {\ displaystyle \ phi}\phi с

R (r) = 1 rn + 1 {\ displaystyle R (r) = {\ frac {1} {r ^ {n + 1}}}}R(r)={\frac {1}{r^{n+1}}}

необходимо иметь это

λ = n (n + 1) {\ displaystyle \ lambda = n (n + 1)}\ lambda = n (n + 1)

Если P n (x) является решением дифференциального уравнения

ddx ((1 - x 2) d P ndx) + n (n + 1) P n = 0 {\ displaystyle {\ frac {d} {dx} } \ left ((1-x ^ {2}) \ {\ frac {dP_ {n}} {dx}} \ right) \ + \ n (n + 1) \ P_ {n} \ = \ 0}{\frac {d}{dx}}\left((1-x^{2})\ {\frac {dP_{n}}{dx}}\right)\ +\ n(n+1)\ P_{n}\ =\ 0

(25)

следовательно, потенциал, соответствующий m = 0

ϕ = 1 rn + 1 P n (sin ⁡ θ) {\ displaystyle \ phi = {\ frac {1} {r ^ {n + 1}}} \ P_ {n} (\ sin \ theta)}\phi ={\frac {1}{r^{n+1}}}\ P_{n}(\sin \theta)

, который является вращательно-симметричным относительно оси z, является гармонической функцией

Если P нм (x) {\ displaystyle P_ {n} ^ {m} (x)}P_{n}^{m}(x)- решение дифференциального уравнения

ddx ((1 - x 2) d P nmdx) + (n ( п + 1) - м 2 1 - Икс 2) п нм = 0 {\ displaystyle {\ frac {d} {dx}} \ left ((1-x ^ {2}) \ {\ frac {dP_ {n} ^ {m}} {dx}} \ right) \ + \ \ left (n (n + 1) - {\ frac {m ^ {2}} {1-x ^ {2}}} \ right) \ P_ {n} ^ {m} \ = \ 0}{\ frac {d} {dx}} \ left ((1-x ^ {2}) \ {\ frac {dP_ {n} ^ {m}} {dx}} \ right) \ + \ \ left (n (n + 1) - {\ frac {m ^ {2}} {1-x ^ {2}}} \ right) \ P_ {n} ^ {m} \ = \ 0

(26)

с m ≥ 1 имеет потенциал

ϕ = 1 rn + 1 P nm (sin ⁡ θ) (a cos ⁡ m φ + b sin ⁡ m φ) {\ displaystyle \ phi = {\ frac {1} {r ^ {n + 1}}} \ P_ {n} ^ {m} (\ sin \ theta) \ (a \ \ cos m \ varphi \ + \ b \ \ sin m \ varphi)}\phi ={\frac {1}{r^{n+1}}}\ P_{n}^{m}(\sin \theta)\ (a\ \cos m\varphi \ +\ b\ \sin m\varphi)

(27)

, где a и b - произвольные константы, является гармонической функцией, которая зависит от φ и поэтому не осесимметрична относительно оси z

Дифференциальное уравнение (25) имеет вид дифференциальное уравнение Лежандра, для которого полиномы Лежандра определяют

P 0 (x) = 1 P n (x) = 1 2 nn! dn (x 2–1) ndxnn ≥ 1 {\ displaystyle {\ begin {выровнено} P_ {0} (x) = 1 \\ P_ {n} (x) = {\ frac {1} {2 ^ {n} n!}} \ {\ frac {d ^ {n} (x ^ {2} -1) ^ {n}} {dx ^ {n}}} \ quad n \ geq 1 \\\ конец {выровнено}} }{\begin{aligned}P_{0}(x)=1\\P_{n}(x)={\frac {1}{2^{n}n!}}\ {\frac {d^{n}(x^{2}-1)^{n}}{dx^{n}}}\quad n\geq 1\\\end{aligned}}

(28)

- решения.

Произвольный коэффициент 1 / (2n!) Выбирается так, чтобы P n (−1) = - 1 и P n (1) = 1 для нечетных n и P n (−1) = P n (1) = 1 для четного n.

Первые шесть полиномов Лежандра:

P 0 (x) = 1 P 1 (x) = x P 2 (x) = 1 2 (3 x 2 - 1) P 3 (x) = 1 2 (5 x 3 - 3 x) P 4 (x) = 1 8 (35 x 4 - 30 x 2 + 3) P 5 (x) = 1 8 (63 x 5 - 70 x 3 + 15 x) {\ Displaystyle {\ begin {align} P_ {0} (x) = 1 \\ P_ {1} (x) = x \\ P_ {2} (x) = {\ frac {1} {2}} \ слева (3x ^ {2} -1 \ справа) \\ P_ {3} (x) = {\ frac {1} {2}} \ left (5x ^ {3} -3x \ right) \\ P_ {4 } (x) = {\ frac {1} {8}} \ left (35x ^ {4} -30x ^ {2} +3 \ right) \\ P_ {5} (x) = {\ frac {1} {8}} \ left (63x ^ {5} -70x ^ {3} + 15x \ right) \\\ end {align}}}{\begin{aligned}P_{0}(x)=1\\P_{1}(x)=x\\P_{2}(x)={\frac {1}{2}}\left(3x^{2}-1\right)\\P_{3}(x)={\frac {1}{2}}\left(5x^{3}-3x\right)\\P_{4}(x)={\frac {1}{8}}\left(35x^{4}-30x^{2}+3\right)\\P_{5}(x)={\frac {1}{8}}\left(63x^{5}-70x^{3}+15x\right)\\\end{aligned}}

(29)

Решения дифференциального уравнения (26) являются ассоциированными функциями Лежандра

P nm (x) = (1 - x 2) m 2 dm P ndxm 1 ≤ m ≤ n {\ displaystyle P_ {n} ^ {m} (x) \ = (1-x ^ {2}) ^ {\ frac {m} {2}} \ {\ frac {d ^ {m} P_ {n}} {dx ^ {m}}} \ quad 1 \ leq m \ leq n}{\ displaystyle P_ {n} ^ {m} (x) \ = (1-x ^ {2}) ^ {\ frac { m} {2}} \ {\ frac {d ^ {m} P_ {n}} {dx ^ {m}}} \ quad 1 \ leq m \ leq n}

(30)

Следовательно,

P nm (sin ⁡ θ) = cos m ⁡ θ dn P ndxn (sin ⁡ θ) {\ displaystyle P_ {n} ^ {m} (\ sin \ theta) = \ cos ^ {m} \ theta \ {\ frac {d ^ {n} P_ {n}} {dx ^ {n}}} (\ sin \ theta)}P_{n}^{m}(\sin \theta)=\cos ^{m}\theta \ {\frac {d^{n}P_{n}}{dx^{n}}}(\sin \theta)
Ссылки
  • Эль Ясберг Теория полета искусственных спутников Земли, Израильская программа научных переводов (1967)
  • Лерх, Ф.Дж., Вагнер, Калифорния, Смит, Д.Э., Сэндсон, М.Л., Браунд, Дж. Э., Ричардсон, Дж. А., "Гравитационное поле" Модели Земли (GEM1 и 2) », Отчет X55372146, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1972
  • Лерх, Ф.Дж., Вагнер, Калифорния, Патни, М.Л., Сандсон, М.Л., Браунд, Дж. Э., Ричардсон, Дж. А., Тейлор, Вашингтон, "Модели гравитационного поля GEM3 и 4", Отчет X59272476, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1972
  • Лерх, Ф.Дж., Вагнер, Калифорния, Ричардсон, Дж. А., Браунд, Дж. Э., «Модели Земли Годдарда (5 и 6)», Отчет X92174145, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1974
  • Лерх, Ф.Дж., Вагнер, Калифорния, Клоско, С.М., Белотт, Р.П., Лаубшер, Р. Рейлор, WA, «Улучшение гравитационной модели с использованием Geos3 Altimetry (GEM10A и 10B)», Весеннее ежегодное собрание Американского геофизического союза 1978 года, Майами, 1978
  • Лерх, Ф.Дж., Клоско, С.М., Ла Убшер, Р.Е., Вагнер, К.А., "Улучшение гравитационной модели с использованием Geos3 (GEM9 и 10)", Журнал геофизических исследований, том. 84, В8, с. 3897-3916, 1979
  • Lerch, F.J., Putney, B.H., Wagner, C.A., Klosko, S.M., "Модели Земли Годдарда для океанографических приложений (GEM 10B и 10C)", Морская геодезия, 5 (2), с. 145-187, 1981
  • Лерх, Ф.Дж., Клоско, С.М., Патель, Великобритания, «Уточненная гравитационная модель из Лагеоса (GEML2)», 'Технический меморандум НАСА 84986, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1983
  • Лерх, Ф.Дж., Нерем, Р.С., Патни, Б.Х., Фелсентрегер, Т.Л., Санчес, Б.В., Клоско, С.М., Патель, Великобритания, Уильямсон, Р.Г., Чинн, Д.С., Чан, Дж.К., Рахлин, К.Э., Чендлер, Н.Л., Маккарти, Дж. Дж., Маршалл, Дж. А., Латке, С.Б., Павлис, Д.В., Роббинс, Дж. У., Капур, С., Павлис, Е.К., "Геопотенциальные модели Земли по данным спутникового слежения, высотомера и наблюдений за поверхностной гравитацией: GEMT3 и GEMT3S », Технический меморандум НАСА 104555, Центр космических полетов Годдарда, Гринбелт / Мэриленд, 1992
  • Лерх, Ф.Дж., Нерем, Р.С., Патни, Б.Н., Фелсентрегер, Т.Л., Санчес, Б.В., Маршалл, Дж. Клоско, С.М., Патель, ГБ, Уильямсон, Р.Г., Чинн, Д.С., Чан, Дж. К., Рахлин, К. Э., Чендлер, Н. Л., Маккарти, Дж. Дж., Лутке, С. Б., Павлис, Н. С., Павлис, ЭЦ, "А" Геопотенциальная модель на основе данных спутникового слежения, высотомера и данных поверхностной силы тяжести: GEMT3 ", Journal of Geophysical Research, Vol. 99, № B2, стр. 2815-2839, 1994
  • Нерем, Р.С., Лерх, Ф.Дж., Маршалл, Дж. А., Павлис, Е.К., Патни, Б.Х., Тэпли, Б.Д., Эанс, Р.Дж., Райс, Дж., Уоткинс, М.М., Клоско, С.М., Чан, Дж. К., Лутке, С.Б., Патель, Г.Б., Павлис, Н.К., Уильямсон, Р.Г., Рапп, Р.Х., Бианкал, Р., Ноуэль, Ф., "Разработка модели гравитации для Topex / Посейдон: Совместные гравитационные модели 1 и 2 », Журнал геофизических исследований, Vol. 99, No. C12, p. 24421-24447, 1994a
Внешние ссылки
Последняя правка сделана 2021-05-21 03:46:58
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте