Английская Википедия:Gauss's method

Шаблон:Short description In orbital mechanics (a subfield of celestial mechanics), Gauss's method is used for preliminary orbit determination from at least three observations (more observations increases the accuracy of the determined orbit) of the orbiting body of interest at three different times. The required information are the times of observations, the position vectors of the observation points (in Equatorial Coordinate System), the direction cosine vector of the orbiting body from the observation points (from Topocentric Equatorial Coordinate System) and general physical data.

Carl Friedrich Gauss developed important mathematical techniques (summed up in Gauss's methods) which were specifically used to determine the orbit of Ceres. The method shown following is the orbit determination of an orbiting body about the focal body where the observations were taken from, whereas the method for determining Ceres' orbit requires a bit more effort because the observations were taken from Earth while Ceres orbits the Sun.

Observer position vector

The observer position vector (in Equatorial coordinate system) of the observation points can be determined from the latitude and local sidereal time (from Topocentric coordinate system) at the surface of the focal body of the orbiting body (for example, the Earth) via either: Шаблон:Multiple image <math display="block">\begin{align} \mathbf{R_n} &= \left [ {R_e \over \sqrt{1-(2f-f^2)\sin^2\phi_n}}+H_n \right ] \cos\phi_n(\cos\theta_n \ \mathbf{\hat{I}}+\sin\theta_n \ \mathbf{\hat{J}})+\left [ {R_e(1-f)^2 \over \sqrt{1-(2f-f^2)\sin^2\phi_n}}+H_n \right ] \sin\phi_n \ \mathbf{\hat{K}} \\

            &= \left [ {R_e \over \sqrt{1-e^2\sin^2\phi_n}}+H_n \right ] \cos\phi_n(\cos\theta_n \ \mathbf{\hat{I}}+\sin\theta_n \ \mathbf{\hat{J}})+\left [ {R_e(1-e^2) \over \sqrt{1-e^2\sin^2\phi_n}}+H_n \right ] \sin\phi_n \ \mathbf{\hat{K}}

\end{align}</math> or <math display="block">\mathbf{R_n} = r_e\cos\phi'_n\cos\theta_n \ \mathbf{\hat{I}}+r_e\cos\phi'_n\sin\theta_n \ \mathbf{\hat{J}}+r_e\sin\phi'_n \ \mathbf{\hat{K}}</math> where,

<math>\mathbf{R_n}</math> is the respective observer position vector (in Equatorial Coordinate System)
<math>R_e</math> is the equatorial radius of the central body (e.g., 6,378 km for Earth)
<math>r_e</math> is the geocentric distance
<math>f</math> is the oblateness (or flattening) of the central body (e.g., 0.003353 for Earth)
<math>e</math> is the eccentricity of the central body (e.g., 0.081819 for Earth)
<math>\phi_n</math> is the geodetic latitude (the angle between the normal line of horizontal plane and the equatorial plane)
<math>\phi'_n</math> is the geocentric latitude (the angle between the radius and the equatorial plane)
<math>H_n</math> is the geodetic altitude
<math>\theta_n</math> is the local sidereal time of observation site

Orbiting body direction cosine vector

Файл:Ra and dec on celestial sphere.png

Right ascension (blue) and declination (green) as seen from outside the celestial sphere

The orbiting body direction cosine vector can be determined from the right ascension and declination (from Topocentric Equatorial Coordinate System) of the orbiting body from the observation points via:

<math display="block">\mathbf{\hat\boldsymbol{\rho}_n} = \cos\delta_n\cos\alpha_n \ \mathbf{\hat{I}}+\cos\delta_n\sin\alpha_n \ \mathbf{\hat{J}}+\sin\delta_n \ \mathbf{\hat{K}}</math> where,

<math>\mathbf{\hat\boldsymbol{\rho}_n}</math> is the respective unit vector in the direction of the position vector <math>\rho</math> (from observation point to orbiting body in Topocentric Equatorial Coordinate System)
<math>\delta_n</math> is the respective declination
<math>\alpha_n</math> is the respective right ascension

Algorithm

The initial derivation begins with vector addition to determine the orbiting body's position vector. Then based on the conservation of angular momentum and Keplerian orbit principles (which states that an orbit lies in a two dimensional plane in three dimensional space), a linear combination of said position vectors is established. Also, the relation between a body's position and velocity vector by Lagrange coefficients is used which results in the use of said coefficients. Then with vector manipulation and algebra, the following equations were derived. For detailed derivation, refer to Curtis.^[1]

NOTE: Gauss's method is a preliminary orbit determination, with emphasis on preliminary. The approximation of the Lagrange coefficients and the limitations of the required observation conditions (i.e., insignificant curvature in the arc between observations, refer to Gronchi^[2] for more details) causes inaccuracies. Gauss's method can be improved, however, by increasing the accuracy of sub-components, such as solving Kepler's equation. Another way to increase the accuracy is through more observations.

Step 1

Calculate time intervals, subtract the times between observations: <math display="block">\begin{align} \tau_1 &= t_1 - t_2 \\ \tau_3 &= t_3 - t_2 \\ \tau &= t_3 - t_1 \end{align}</math> where

<math>\tau_n</math> is the time interval
<math>t_n</math> is the respective observation time

Step 2

Файл:Cross product vector.svg

The cross product in respect to a right-handed coordinate system

Calculate cross products, take the cross products of the observational unit direction (order matters): <math display="block">\begin{align} \mathbf{p_1} &= \mathbf{\hat\boldsymbol{\rho}_2} \times \mathbf{\hat\boldsymbol{\rho}_3} \\ \mathbf{p_2} &= \mathbf{\hat\boldsymbol{\rho}_1} \times \mathbf{\hat\boldsymbol{\rho}_3} \\ \mathbf{p_3} &= \mathbf{\hat\boldsymbol{\rho}_1} \times \mathbf{\hat\boldsymbol{\rho}_2} \end{align}</math> where

<math>\mathbf{a} \times \mathbf{b}</math> is the cross product of vectors <math>\mathbf{a} \text{ and } \mathbf{b}</math>
<math>\mathbf{p_n}</math> is the respective cross product vector
<math>\mathbf{\hat\boldsymbol{\rho}_n}</math> is the respective unit vector

Step 3

Файл:Parallelepiped volume.svg

</math>, describes the volume.

Calculate common scalar quantity (scalar triple product), take the dot product of the first observational unit vector with the cross product of the second and third observational unit vector:

<math display="block">D_0 = \mathbf{\hat\boldsymbol{\rho}_1} \cdot \mathbf{p_1} = \mathbf{\hat\boldsymbol{\rho}_1} \cdot (\mathbf{\hat\boldsymbol{\rho}_2} \times \mathbf{\hat\boldsymbol{\rho}_3})</math> where

<math>\mathbf{a} \cdot \mathbf{b}</math> is the dot product of vectors <math>\mathbf{a}</math> and <math>\mathbf{b}</math>
<math>D_0</math> is the common scalar triple product
<math>\mathbf{p_n}</math> is the respective cross product vector
<math>\mathbf{\hat\boldsymbol{\rho}_n}</math> is the respective unit vector

Step 4

Calculate nine scalar quantities (similar to step 3): <math display="block">\begin{align} D_{11} &= \mathbf{R_1} \cdot \mathbf{p_1} & D_{12} &= \mathbf{R_1} \cdot \mathbf{p_2} & D_{13} &= \mathbf{R_1} \cdot \mathbf{p_3} \\ D_{21} &= \mathbf{R_2} \cdot \mathbf{p_1} & D_{22} &= \mathbf{R_2} \cdot \mathbf{p_2} & D_{23} &= \mathbf{R_2} \cdot \mathbf{p_3} \\ D_{31} &= \mathbf{R_3} \cdot \mathbf{p_1} & D_{32} &= \mathbf{R_3} \cdot \mathbf{p_2} & D_{33} &= \mathbf{R_3} \cdot \mathbf{p_3} \end{align}</math> where

<math>D_{mn}</math> is the respective scalar quantities
<math>\mathbf{R_m}</math> is the respective observer position vector
<math>\mathbf{p_n}</math> is the respective cross product vector

Step 5

Calculate scalar position coefficients:

<math display="block">\begin{align} A &= \frac{1}{D_0} \left ( -D_{12} \frac{\tau_3}{\tau}+D_{22}+D_{32} \frac{\tau_1}{\tau} \right ) \\ B &= \frac{1}{6D_0} \left [ D_{12} \left ( \tau_3^2-\tau^2 \right ) \frac{\tau_3}{\tau}+D_{32} \left (\tau^2-\tau_1^2 \right ) \frac{\tau_1}{\tau} \right ] \\ E &= \mathbf{R_2} \cdot \mathbf{\hat\boldsymbol{\rho}_2} \end{align}</math> where

<math>A</math>, <math> B </math>, and <math> E</math> are scalar position coefficients
<math>D_0</math> is the common scalar quantity
<math>D_{mn}</math> is the respective scalar quantities
<math>\tau_n</math> is the time interval
<math>R_n</math> is the respective observer position vector
<math>\mathbf{\hat\boldsymbol{\rho}_n}</math> is the respective unit vector

Step 6

Calculate the squared scalar distance of the second observation, by taking the dot product of the position vector of the second observation: <math display="block">{R_2}^2 = \mathbf{R_2} \cdot \mathbf{R_2}</math> where

<math>{R_2}^2</math> is the squared distance of the second observation
<math>\mathbf{R_2}</math> is the position vector of the second observation

Step 7

Calculate the coefficients of the scalar distance polynomial for the second observation of the orbiting body: <math display="block">\begin{align} a &= -\left ( A^2+2AE+{R_2}^2 \right ) \\ b &= -2\mu B(A+E) \\ c &= -\mu^2 B^2 \end{align}</math> where

<math>a \text{, } b \text{, and } c</math> are coefficients of the scalar distance polynomial for the second observation of the orbiting body
<math>A \text{, } B \text{, and } E</math> are scalar position coefficients
<math>\mu</math> is the gravitational parameter of the focal body of the orbiting body

Step 8

Find the root of the scalar distance polynomial for the second observation of the orbiting body: <math display="block"> {r_2}^8+a{r_2}^6+b{r_2}^3+c = 0 </math> where

<math>r_2</math> is the scalar distance for the second observation of the orbiting body (it and its vector, r₂, are in the Equatorial Coordinate System)
<math>a \text{, } b \text{, and } c</math> are coefficients as previously stated

Various methods can be used to find the root, a suggested method is the Newton–Raphson method. The root must be physically possible (i.e., not negative nor complex) and if multiple roots are suitable, each must be evaluated and compared to any available data to confirm their validity.

Step 9

Calculate the slant range, the distance from the observer point to the orbiting body at their respective time: <math display="block">\begin{align} \rho_1 &= \frac{1}{D_0} \left [ \frac{6 \left ( D_{31} \dfrac{\tau_1}{\tau_3}+D_{21} \dfrac{\tau}{\tau_3} \right ) {r_2}^3+\mu D_{31} \left ( \tau^2-{\tau_1}^2 \right ) \dfrac{\tau_1}{\tau_3}}{6{r_2}^3+\mu \left ( \tau^2-{\tau_3}^2 \right ) } - D_{11}\right ] \\ \rho_2 &= A+ \frac{\mu B}{{r_2}^3} \\ \rho_3 &= \frac{1}{D_0} \left [ \frac{6 \left ( D_{13} \dfrac{\tau_3}{\tau_1}-D_{23} \dfrac{\tau}{\tau_1} \right ) {r_2}^3+\mu D_{13} \left ( \tau^2-{\tau_3}^2 \right ) \dfrac{\tau_3}{\tau_1}}{6{r_2}^3+\mu \left ( \tau^2-{\tau_1}^2 \right ) } - D_{33}\right ] \end{align}</math> where

<math>\rho_n</math> is the respective slant range (it and its vector, <math>\mathbf{\rho_n}</math>, are in the Topocentric Equatorial Coordinate System)
<math>D_0</math> is the common scalar quantity
<math>D_{mn}</math> is the respective scalar quantities
<math>\tau_{(n)}</math> is the time interval.
<math>r_2</math> is the scalar distance for the second observation of the orbiting body
<math>\mu</math> is the gravitational parameter of the focal body of the orbiting body

Step 10

Calculate the orbiting body position vectors, by adding the observer position vector to the slant direction vector (which is the slant distance multiplied by the slant direction vector):

<math display="block">\begin{align} \mathbf{r_1} &= \mathbf{R_1}+\rho_1\mathbf{\hat\boldsymbol{\rho}_1} \\[1.7ex] \mathbf{r_2} &= \mathbf{R_2}+\rho_2\mathbf{\hat\boldsymbol{\rho}_2} \\[1.7ex] \mathbf{r_3} &= \mathbf{R_3}+\rho_3\mathbf{\hat\boldsymbol{\rho}_3} \end{align}</math> where

<math>\mathbf{r_n}</math> is the respective orbiting body position vector (in Equatorial Coordinate System)
<math>\mathbf{R_n}</math> is the respective observer position vector
<math>\rho_n</math> is the respective slant range
<math>\mathbf{\hat\boldsymbol{\rho}_n}</math> is the respective unit vector

Step 11

Calculate the Lagrange coefficients: <math display="block">\begin{align} f_1 &\approx 1-\frac{1}{2}\frac{\mu}{{r_2}^3}{\tau_1}^2 \\ f_3 &\approx 1-\frac{1}{2}\frac{\mu}{{r_2}^3}{\tau_3}^2 \\ g_1 &\approx \tau_1-\frac{1}{6}\frac{\mu}{{r_2}^3}{\tau_1}^3 \\ g_3 &\approx \tau_3-\frac{1}{6}\frac{\mu}{{r_2}^3}{\tau_3}^3 \end{align}</math> where,

<math>f_1</math>, <math>f_3</math>, <math>g_1</math> and <math>g_3</math> are the Lagrange coefficients (these are just the first two terms of the series expression based on the assumption of small time interval)
<math>\mu</math> is the gravitational parameter of the focal body of the orbiting body
<math>r_2</math> is the scalar distance for the second observation of the orbiting body
<math>\tau_{(n)}</math> is the time interval

Step 12

Calculate the velocity vector for the second observation of the orbiting body:

<math display="block">\mathbf{v_2} = \frac{1}{f_1 g_3 - f_3 g_1} \left(-f_3\mathbf{r_1}+f_1\mathbf{r_3}\right)</math> where

<math>\mathbf{v_2}</math> is the velocity vector for the second observation of the orbiting body (in Equatorial Coordinate System)
<math>f_1</math>, <math>f_3</math>, <math>g_1</math> and <math>g_3</math> are the Lagrange coefficients
<math>\mathbf{r_n}</math> is the respective orbiting body position vector

Step 13

The orbital state vectors have now been found, the position (r₂) and velocity (v₂) vector for the second observation of the orbiting body. With these two vectors, the orbital elements can be found and the orbit determined.

References

Шаблон:Reflist

Der, Gim J.. "New Angles-only Algorithms for Initial Orbit Determination." Advanced Maui Optical and Space Surveillance Technologies Conference. (2012). Print.

Шаблон:Carl Friedrich Gauss

↑ Curtis, Howard D. Orbital Mechanics for Engineering Students. Oxford: Elsevier Butterworth-Heinemann, 2005. Print.
↑ Gronchi, Giovanni F.. "Classical and modern orbit determination for asteroids." Proceedings of the International Astronomical Union2004.IAUC196 (2004): 1-11. Print.

[1] Curtis, Howard D. Orbital Mechanics for Engineering Students. Oxford: Elsevier Butterworth-Heinemann, 2005. Print.

[2] Gronchi, Giovanni F.. "Classical and modern orbit determination for asteroids." Proceedings of the International Astronomical Union2004.IAUC196 (2004): 1-11. Print.

[1]

[2]

Партнерские ресурсы
Криптовалюты	Обмен криптовалют - www.bestchange.ru Криптовалютная биржа CoinEx Криптовалютная биржа Binance HIVE OS - операционная система для майнинга e4pool - Мультивалютный пул для майнинга.
Магазины	AliExpress — глобальная виртуальная (в Интернете) торговая площадка, предоставляющая возможность покупать товары производителей из КНР; computeruniverse.net - Интернет-магазин компьютеров(Промо код 5 Евро на первую покупку:FWWC3ZKQ);
Хостинг	DigitalOcean - американский провайдер облачных инфраструктур, с главным офисом в Нью-Йорке и с центрами обработки данных по всему миру;
Разное	Викиум - Онлайн-тренажер для мозга Like Центр - Центр поддержки и развития предпринимательства. Gamersbay - лучший магазин по бустингу для World of Warcraft. Ноотропы OmniMind N°1 - Усиливает мозговую активность. Повышает мотивацию. Улучшает память. Санкт-Петербургская школа телевидения - это федеральная сеть образовательных центров, которая имеет филиалы в 37 городах России. Lingualeo.com — интерактивный онлайн-сервис для изучения и практики английского языка в увлекательной игровой форме. Junyschool (Джунискул) – международная школа программирования и дизайна для детей и подростков от 5 до 17 лет, где ученики осваивают компьютерную грамотность, развивают алгоритмическое и креативное мышление, изучают основы программирования и компьютерной графики, создают собственные проекты: игры, сайты, программы, приложения, анимации, 3D-модели, монтируют видео. Умназия - Интерактивные онлайн-курсы и тренажеры для развития мышления детей 6-13 лет SkillBox - это один из лидеров российского рынка онлайн-образования. Среди партнеров Skillbox ведущий разработчик сервисного дизайна AIC, медиа-компания Yoola, первое и самое крупное русскоязычное аналитическое агентство Tagline, онлайн-школа дизайна и иллюстрации Bang! Bang! Education, оператор PR-рынка PACO, студия рисования Draw&Go, агентство performance-маркетинга Ingate, scrum-студия Sibirix, имидж-лаборатория Персона. «Нетология» — это университет по подготовке и дополнительному обучению специалистов в области интернет-маркетинга, управления проектами и продуктами, дизайна, Data Science и разработки. В рамках Нетологии студенты получают ценные теоретические знания от лучших экспертов Рунета, выполняют практические задания на отработку полученных навыков, общаются с экспертами и единомышленниками. Познакомиться со всеми продуктами подробнее можно на сайте https://netology.ru, линейка курсов и профессий постоянно обновляется. StudyBay Brazil – это онлайн биржа для португалоговорящих студентов и авторов! Студент получает уникальную работу любого уровня сложности и больше свободного времени, в то время как у автора появляется дополнительный заработок и бесценный опыт. Автор24 — самая большая в России площадка по написанию учебных работ: контрольные и курсовые работы, дипломы, рефераты, решение задач, отчеты по практике, а так же любой другой вид работы. Сервис сотрудничает с более 70 000 авторов. Более 1 000 000 работ уже выполнено. StudyBay – это онлайн биржа для англоязычных студентов и авторов! Студент получает уникальную работу любого уровня сложности и больше свободного времени, в то время как у автора появляется дополнительный заработок и бесценный опыт.