Местоположение спутника

[Tau]

Я - блондинка, я не понимаю.
Время прохождения восходящего узла - постоянная величина на каждом витке, то есть в течение 10 с лишним часов. Как может разность времени расчета и времени прохождения узла изменяться нелинейно?

[aleksa-chizhi]

Да, вы правы, на рис. 2.5 показано расчетное время прохождения долготы восходящего узла на n-ом витке (некие ступеньки) и оно не изменяется в течении Tобр НКА ГЛОНАСС, на рис 2.1 показано время наблюдения ti, выраженное в секундах от начала суток наблюдения (пила). Если из значений графика на рис. 2.1 вычесть значения графика на рис. 2.5, то получится то что имеем на рис. 2.6.

[Tau]

Ага, я наконец-то поняла, что по оси абсцисс у вас такой масштаб: 300 сек. на единицу измерения.
Действительно есть скачки в изменении координат (рис. x_eci.jpg).
Разность tau=Ti-Tlk и должна изменяться скачками, потому что эта разность всегда должна быть в пределах от 0 до величины периода обращения спутника. Эта разность не должна быть < 0, а она у нас и в отрицательную область скачет. Это эффект неравномерной шкалы времени: при переходе через начало суток возникает отрицательная разность, если последнее прохождение восходящего узла было в предыдущую дату, а Ti измеряется в секундах текущей даты.
Это лечится введением условия: если tau < 0, то tau = tau + 86400
Изменение координат стало ровным (рис. x_eci2.jpg)

[aleksa-chizhi]

Спасибо, приятно иметь дело с профессионалом. Значит я был на правильном пути поиска бага. Я исправил программу по Вашей подсказке, ВСЕ РАБОТАЕТ ТЕПЕРЬ ПРАВИЛЬНО. Интересно почему в ИКД это не оговорено? Еще раз спасибо!!!

[Tau]

На здоровье .

Интересно почему в ИКД это не оговорено?

Вероятно, авторам казалось, что это очевидно.

[torx]

уважаемая тау, а у вас нет случайно еще примера расчета параметров движеия НКА по данным эфемерид? просто пишу программку для перерасчета а свериться в полученных данных не с чем

[Tau]

torx, смотрите прицепленный файл.

P.S. Звездное время использовано истинное.
P.P.S. В альманахе под единицей измерения "полуцикл" подразумевается коэффициент при "пи" (количество пи).

[torx]

спасибо

[t_rex]

Уважаемая Tau, Вы чуть ли не единственная, кто есть в рунете на форумах и действительно "шарит" в теме    не подскажете такой момент, набили с коллегой алгоритм расчета по альманаху согласно ИКД 5.1, для проверки использовали данные из примера того же ИКД. Результаты получаются совсем разные. механические ошибки исключили сравнением текстов программ и параллельной отладкой. значит ошибки логические.
что делали: взяли данные альманаха, проделали с ними все, что описано в ИКД, результаты другие. может нужно переводить в другую СК? или есть какие-то неочевидные для не-астрономов хитрости? если знаете - подскажите, пожалуйста.

вот кусок из примера с альманаха и наш результат (в метрах)

[Tau]

Углы, данные в альманахе в полуциклах, на ПИ умножили, чтобы в радианы перевести?

[t_rex]

Умножил с Вашего совета, получилось!!! Tau,  спасибо Вам большое человеческое, сами бы до этого дооооолго доходили 

[t_rex]

Уважаемая Tau, а не подскажете, где брать нутации для вычисления истинного звёздного времени?

[Tau]

t_rex, http://law.rufox.ru/view/standarti/1672.htm

[Веселый гном]

2) Вычисление среднего звездного гринвического времени на полночь даты MD
 

A1 = 24110.54841
 A2 = 8640184.812
 A3 = 0.093104
 A4 = 0.0000062
 T0 = (MD - 51544.5) / 36525
 Sg = A1 + A2 * T0 + A3 * T0 ^ 2 - A4 * T0 ^ 3
 Sg - звездное время в секундах времени.
 S0 = Sg *15 / 3600 ' в градусах

Пример: 14 января 2010
MD = 55210
Sg=  891204.250330178 sec
S0=  3713.35104304241 deg = 113.351043042409 deg

На 14 января 2010 у меня результат сошелся с вашим S0=  3713.35104304241
Но в примере интерфейсного документа Глонасс на дату 06.09.2001  S0 = 6.02401539573, а по вашей формуле получается 705,15528. В чем ошибка?

[Tau]

Но в примере интерфейсного документа Глонасс на дату 06.09.2001  S0 = 6.02401539573, а по вашей формуле получается 705,15528. В чем ошибка?

Приведите звездное время 705,15528 в диапазон от 0 до 360 градусов, потом переведите в радианы и сравнивайте с S0 из интерфейсного документа.

[dmtr]

Кстати, слышал что ИКД для ГЛОНАСС содержат ошибки и опечатки (в том числе и ИКД 5.1 от 2008 года).
Я набрал алгоритм, который там предлагается, но поправки у меня не сошлись.
Если кто уже пробовал проделать тоже самое, скажите: там правда опечатки есть или я просто такой криволапый?)

[Tau]

Да тут уж толпа народа набирала этот алгоритм, и ничего - всё сошлось. Явных опечаток не обнаружено.
А что вы подразумеваете под поправками?

[Веселый гном]

Да тут уж толпа народа набирала этот алгоритм

А может даже есть где готовые коды на си или делфи? Конкретно расчет координат КА используя все нутации-шмутации и поправки?

[dmtr]

Грубо параметры орбиты сходятся, то есть у меня координаты в конце концов не сходятся на 200-300 метров.
Я так считаю что не сходятся поправки к элементам орбиты: di (поправка к наклонению), da (поправка к полуоси) и другие. Именно: пункты 3 и 4 в алгоритме рассчета по альманаху.

[Веселый гном]

Рассчитал координаты Глонасов в абсолютной геоцентрической системе координат как букваре ICD_GLONASS_rus.pdf
А как теперь перевести координаты в человеческий вид: широта, долгота, высота?