Координаты спутника по TLE

[Archi]

ke = √GM where G is Newton’s universal gravitational constant and M is the mass
of the Earth - если GM =    398600.4418, то корень от него = 631.3481б

[Sleepwalker]

ke(XKE) уже определена в таблице констант, ничего привносить не надо.

[Archi]

Тогда a1 = (ke/n0)^(2/3)
ke = .743669161E-1
n0 = 0.061 радианы в минуту   

a1 = 1.137492, a0 = 1.13764, a0'' = 1.13718. И потом во всех формулах с a0'' используется эта цифра = 1.13718.

А эти все - Наклонение в градусах, Долгота восходящего узла в градусах, Аргумент перицентра в градусах,Средняя аномалия в градусах - их всех в формулах брать в радианах?

[Sleepwalker]

да, все углы в радианах

[Archi]

А радианы в рубежах от 0 ... 2pi , или -2pi ... 2pi тоже годится?

[Sleepwalker]

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

[Archi]

А эпоха: если t = 4133.83, а t0 = 4126.75, то там где в формулах (t-t0) = (4133.83-4126.75,) =  7.08?

[Tau]

А эпоха: если t = 4133.83, а t0 = 4126.75, то там где в формулах (t-t0) = (4133.83-4126.75,) =  7.08?

Откуда вы взяли такую эпоху и в чем она измеряется? dt должна быть в минутах, это промежуток времени между эпохой, на которую даны TLE, и моментом расчета. В документе есть примеры результатов расчета, возьмите конкретный пример и тренируйтесь на нем, пока не совпадет с приведенными там числами.

[Archi]

t это сколько дней прошло от эпохи 1 января 2000 года.
А в примере показано только, как должно быть в конце - вектор x, y, z и вектор xdot, ydot, zdot. Что посередине не известно. Мне надо что-бы можно сравнить все величины начиная от a1 до velocity. Может где есть такая возможность?

[Sleepwalker]

Самое простое это взять готовый код и сравнить.

[Archi]

А квадратные скобки в формулах  - это не модуль (если минус знак то отбрасывать его)?

[Tau]

Нет, это просто скобки,  для лучшей читабельности формул. Модулей там вообще нигде нет.

[Archi]

Я тут посмотрел код - там константы дефинированы чуть иначе - например J3 там -2.53881E-6, а в документе -.253881E-5 (запятую перенес), и q0 и S сам считает, а не берет из файла.

[Sleepwalker]

какой код?
не вижу разницы в J3.
Что за файл, из которого берутся q0,S ?

[Archi]

C++ код от сюда: http://www.zeptomoby.com/satellites/cppOrbitTools.zip, только скомпилировать не удалось через Microsoft Visual Studio 2010.

Еще нашел кое-что - например где считают δ1, там по формуле надо cos(i0)^2 (косинус от i в квадрате), а в этом коде sqr(cos(i) - тянет корень от i.

[Sleepwalker]

корень - sqrt
sqr - возведение в квадрат

в VS2005 собралось и отработало.

[Archi]

1) Не совпадает S если перигей >98 и <156. У него там в формуле m_s4 = perigee - 78.0, а в документе формула другая. От куда этот 78?

2) A3,0 по формуле -J3 * (AE^3), а тут A3,0 = -J3/CK2*(AE^3). Почему еще на CK2 делится?

[Sleepwalker]

1) Не совпадает S если перигей >98 и <156. У него там в формуле m_s4 = perigee - 78.0, а в документе формула другая. От куда этот 78?

 в документе просто описывается сразу в безразмерных величинах, а в коде сначала в размерных, а потом идет деление на радиус Земли
78 идет отсюда
s* = a``0*(1-e0)  -  s + ae
где s  = 1.01222...,
ae = 1
в радиусах Земли (-s+ae) = -78

[Sleepwalker]

2) A3,0 по формуле -J3 * (AE^3), а тут A3,0 = -J3/CK2*(AE^3). Почему еще на CK2 делится?

При расчете c3 кусок делителя(k2==CK2) перекочевал в А3.0

[Archi]

Там где Velocity вектор умножается на радиусс земли и потом еще на (1440/86400) - как эту формулу понять, почему умножается на минуты в день и делится на секунды в день?