Вычисление аргумента перицентра ИСЗ

[GreenEkatherine]

Совершенно непонятно, что вам в результате нужно  . Проверить правильность какой-то диссертации или вычислить элементы орбиты из реальных данных?

Вычислить орбиту, конечно же. Диссертацию я упомянула как источник формул.

5) Казалось, что на таком небольшом промежутке возмущения не должны быть заметны, спасибо за пояснение.
1) Известные значения высот взяты из открытых источников, как они были получены, я не знаю.
2), 4) Поняла, буду вычислять по книге.
3) Наклонение также брала из открытых источников, о давности ничего не оговаривалось.

[Upsilon]

В задачах движения ИСЗ довольно много алгоритмов и системы координат используются разные, вот что для решения конкретной задачи удобнее - то и используется. Поэтому диссертация как источник формул - не очень удачный вариант, автор диссертации использовал те формулы, которые ему более подходят для решения его задачи. Учиться лучше по учебникам .

5) Казалось, что на таком небольшом промежутке возмущения не должны быть заметны, спасибо за пояснение.

Так там и возмущения небольшие. И опять-таки, смотря откуда вы взяли вектора положения и скорости. Для их вычисления используют разные способы вычисления (аналитические разложения, интегрирование, смешанные методы). Если бы эти данные были вычислены в предположении невозмущенного движения, изменялся бы только один параметр -  истинная аномалия. А в возмущенном движении шевелятся все параметры, это нормально.

[GreenEkatherine]

Векторы положения и скорости якобы списываются со спутника. Каким методом они там рассчитываются – хороший вопрос. Спасибо, Upsilon, вы мне очень помогли. Впервые за время общения на форумах удалось получить грамотную консультацию 
Буду дальше разбираться.

[xd]

Если прямо со спутника, значит с его подсистемы навигации, вероятно GPS.
Работайте лучше с большой полуосью орбиты, высоты - это профанация

[GreenEkatherine]

Работайте лучше с большой полуосью орбиты, высоты - это профанация

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

[xd]

У нас есть две величины, связанные между собой через аддитивную константу, взятую от балды. Абсолютные погрешности остаются теми же, что вразнос может пойти?

[GreenEkatherine]

У нас есть две величины, связанные между собой через аддитивную константу, взятую от балды.

Хорошо, буду работать с большой полуосью. На первых порах казалось, что сообщили, то и нужно использовать, а сообщили высоты.

[GreenEkatherine]

Проглядела одну важную вещь. Векторы скорости и положения, как говорится в приведенной книге, должны быть в экваториальной системе координат, а я приводила их координаты в Гринвичской системе координат.

Полученные элементы орбиты, а именно долготу восходящего узла, наклонение орбиты и аргумент широты, я подставляю в матрицу поворота для перехода из орбитальной системы координат в инерциальную геоцентрическую (назовем ее M_ин).
Далее по звездному времени я строю матрицу поворота для перехода из инерциальной СК в Гринвичскую (M_гр).

Вопрос: если вычислять элементы орбиты для матрицы поворота по векторам скорости и положения, заданным в Гринвичской СК, получится ли сразу матрица вида M_гр*M_ин или же надо честно сначала переводить векторы в инерциальную СК, а затем вычислять обе матрицы?

[xd]

Можно произведение использовать, главное порядок соблюдать (умножение матриц некоммутативно)

[Upsilon]

Полученные элементы орбиты, а именно долготу восходящего узла, наклонение орбиты и аргумент широты, я подставляю в матрицу поворота для перехода из орбитальной системы координат в инерциальную геоцентрическую (назовем ее Mин).

Что вы подразумеваете под орбитальной системой координат в данном случае??

[GreenEkatherine]

Что вы подразумеваете под орбитальной системой координат в данном случае??

Систему координат с центром, совпадающим с центром масс спутника, ось OX совпадает с вектором движения спутника, OZ направлена от центра масс Земли, а OY дополняет систему векторов до правой.

[GreenEkatherine]

Можно произведение использовать, главное порядок соблюдать (умножение матриц некоммутативно)

Вопрос в том, будет ли полученная матрица матрицей произведения, если элементы орбиты посчитаны по векторам в Гринвиче?

[Upsilon]

Систему координат с центром, совпадающим с центром масс спутника, ось OX совпадает с вектором движения спутника, OZ направлена от центра масс Земли, а OY дополняет систему векторов до правой.

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

[GreenEkatherine]

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

Полный переход я выполняю из СК фокальной плоскости в Гринвичскую СК, требуются три матрицы поворота: 1) переход фокальная СК -> орбитальная СК; 2) орбитальная СК -> инерциальная СК; 3) инерциальная СК -> Гринвичская СК.

Геоцентрические (в Гринвиче) векторы положения и скорости использовались для вычисления элементов орбиты, которые в свою очередь должны участвовать в матрице 2)

[Upsilon]

Спасибо за разъяснение. Нетривиально.

[GreenEkatherine]

Спасибо за разъяснение. Нетривиально.

Не за что =) А как было бы тривиально?

[vitalyn]

Большая полуось, эксцентриситет, наклонение, долгота восх.узла, аргумент перигея, аргумент широты
---
Углы напечатаны в градусах, большая полуось в километрах.
Сравнивайте со своим результатом.

А у меня, кстати, вот что получилось:
7027.017   0.02038951     99.892     48.451     40.428     36.341
7026.938   0.02038090     99.886     48.447     40.513     36.430
7026.886   0.02037484     99.881     48.443     40.605     36.518
7026.804   0.02036610     99.876     48.439     40.695     36.607
7026.732   0.02035829     99.870     48.435     40.786     36.695
7026.649   0.02034869     99.865     48.431     40.871     36.784
7026.575   0.02034071     99.859     48.426     40.962     36.872
7026.505   0.02033312     99.854     48.422     41.046     36.961
7026.439   0.02032552     99.848     48.418     41.133     37.049
7026.354   0.02031631     99.843     48.414     41.222     37.137
7026.277   0.02030736     99.837     48.410     41.313     37.226
7026.190   0.02029799     99.832     48.405     41.402     37.314
7026.124   0.02029077     99.826     48.401     41.487     37.403
7026.051   0.02028293     99.821     48.397     41.579     37.491
7025.995   0.02027638     99.815     48.393     41.666     37.580
7025.909   0.02026689     99.810     48.388     41.757     37.668

Близко, но всё ж таки отличие неприемлемо большое. Причём, я у себя проверил - преобразовал обратно из элементов в декартовы координаты - совпало с исходными координатами.

[GreenEkatherine]

vitalyn, следует заметить, что наклонение, долгота восходящего узла и аргумент широты в Ваших расчетах совпадают, а именно они и нужны для расчета матрицы поворота из орбитальной СК в геоцентрическую, то есть для данной задачи остальные расхождения не принципиальны. Но спасибо за указание.