Определение ПСТ по TLE

[art193]

   Convert_Sat_State преобразовывает вектор положения, а вот во что - сказать не могу - не силен в астрономии. На Паскале процедура выглядит следующим образом:

Procedure Convert_Sat_State(var pos,vel : vector);
  var
    i : byte;
  begin
  for i := 1 to 3 do
    begin
    pos := pos*xkmper;    {kilometers}
    vel := vel*xkmper/60; {kilometers/second}
    end; {for i}
  Magnitude(pos);
  Magnitude(vel);
  end; {Procedure Convert_Sat_State}
end.

  theta почти сошлась  0,000399226429171477

координаты получились следующие:
XPOS  -0,571060703214096 широта
YPOS  199,752753185864    долгота
ZPOS  35970,1164407167

[Sleepwalker]

Все верно.

[art193]

   Да, точно... Вот что значит 00:40 местного времени, мозг отключился...
  Спасибо огромное за помощь!!! :-).

[Sleepwalker]

Всегда пожайлуста.

[art193]

Что-то опять не соответствует... Может еще разочек проверимся???

GALAXY 18 (G-18)       

1 32951U 08024A   09052.44886927  .00000023  00000-0  10000-3 0  1337
2 32951   0.0064 130.5546 0002875 206.5844 212.9170  1.00273392  2794

XPOS  0,00551897039944845 широта
YPOS  29,767980158358       долгота
ZPOS  35795,8131931231

Хотя реально должна быть другая...

[Sleepwalker]

А где все входные параметры и промежуточные результаты?

[art193]

Пример TLE для GALAXY 18:

1 32951U 08024A   09052.44886927  .00000023  00000-0  10000-3 0  1337
2 32951   0.0064 130.5546 0002875 206.5844 212.9170  1.00273392  2794

Результаты расчета SDP4:
X     -6,51104935745225
Y     -1,15253858093781
Z     0,00063627743785033

Результаты расчета Convert_Sat_State
X     -41528,3517934937
Y     -7351,04666192978
Z     4,05826339606351

theta=-2,96639474115606

time: 27.02.2009 0:29:06

XPOS  0,00551897039944845
YPOS  25,7820699367166
ZPOS  35795,8131931231

Вот...

[Sleepwalker]

time: 27.02.2009 0:29:06

UTC, MSK, ... ?

Предлагаю для отброса подозрений на пропагатор для начала проверить на эпоху TLE:
21/02/09 10:46:22.305 2454883.9489039932

[art193]

    Все работает, но как-то интересно получается... Если расчитывать ПСТ по эпохе, указанной в ТЛЕ, тогда все нормально получается, а если на другое время - криво... Как можно рассчитать на другую дату? (в пределах недельки от эпохи в ТЛЕ на падобие функции "симуляция" в Орбитроне) А данные не те получались, потому что я в качестве времени ставил текущее, а не эпоху...

[Sleepwalker]

Значит ошибка при реализации алгоритма пропагатора. Какова скорость роста ошибки от времени прогноза? И ПСТ пока можно отложить. Сравнивайте прямоугольные координаты. В Heavensat'e на последней закладке в информации по спутнику выдаются координаты XYZ с выхода пропагатора.

[art193]

 Ежедневно на 0,986 приблизительно.

[art193]

 ...градуса

[TazmanianDevil]

Добрый день, Господа Знатоки. Случилось так, что у меня тоже возникла проблема определения подспутниковой точки ИСЗ по TLE. В моем случае это необходимо только для геостационарных ИСЗ. Я честно почитал тему, и у меня возник ряд вопросов:
1) Есть программа расчета из TLE в координаты по модели SGP4 http://www.centerforspace.com/downloads/files/sgp4/AIAA-2006-6753Code.zip . Что за данные она возвращает? Точнее, что это за координаты?
2) Как от этих координат перейти к подспутниковой точке?

[Sleepwalker]

Добрый день, Господа Знатоки. Случилось так, что у меня тоже возникла проблема определения подспутниковой точки ИСЗ по TLE. В моем случае это необходимо только для геостационарных ИСЗ. Я честно почитал тему, и у меня возник ряд вопросов:
1) Есть программа расчета из TLE в координаты по модели SGP4 http://www.centerforspace.com/downloads/files/sgp4/AIAA-2006-6753Code.zip . Что за данные она возвращает? Точнее, что это за координаты?

http://www.centerforspace.com/downloads/files/pubs/AIAA-2006-6753.pdf Coordinate System

2) Как от этих координат перейти к подспутниковой точке?

Читать тему еще раз до достижения эффекта.

[Tau]

В моем случае это необходимо только для геостационарных ИСЗ.
1) Есть программа расчета из TLE в координаты по модели SGP4

SGP - для низких ИСЗ.

[TazmanianDevil]

В моем случае это необходимо только для геостационарных ИСЗ.
1) Есть программа расчета из TLE в координаты по модели SGP4

SGP - для низких ИСЗ.

А какой моделью тогда можно пользоваться?

[art193]

В моем случае это необходимо только для геостационарных ИСЗ.
1) Есть программа расчета из TLE в координаты по модели SGP4

SGP - для низких ИСЗ.

А какой моделью тогда можно пользоваться?

SDP4

[TazmanianDevil]

В моем случае это необходимо только для геостационарных ИСЗ.
1) Есть программа расчета из TLE в координаты по модели SGP4

SGP - для низких ИСЗ.

А какой моделью тогда можно пользоваться?

SDP4

Спасибо большое

[Sleepwalker]

в существующей реализации модели с www.centerforspace.com уже нет деления на sgp4 sdp4, все свалили в одну кучу и алгоритм сам переключается в зависимости от среднего движения.

[TazmanianDevil]

Все, решил поставленную задачу. Спасибо большое господам знатокам за консультацию!