Я задаю этот вопрос в связи с примером, показанным в книге Говарда Кертиса «Орбитальная механика для студентов-инженеров» (третье издание), стр. 659, пример 12.1.
Пример демонстрирует, как рассчитать затухание орбиты с учетом кеплеровских элементов и баллистического коэффициента. Процесс происходит примерно так:
Преобразование элементов Кеплера ==> [R, V] Вектор состояния в перифокальной системе отсчета ==> [R, V] вектор состояния в геоцентрической экваториальной системе отсчета ==> численное решение уравнения движения:
где p — возмущающее ускорение сопротивления:
После численного решения ОДУ строится график зависимости высоты от времени.
Мой вопрос: зачем нужно преобразовывать кеплеровские элементы в вектор состояния в геоцентрической экваториальной системе отсчета. Я понимаю, что представление вектора состояния необходимо для работы с уравнениями движения, но почему нельзя просто решить уравнения движения с векторами состояния в перифокальной системе отсчета? Есть ли особая причина/преимущество в этом?
Я не специалист, но вот мои догадки: геоцентрическая экваториальная система отсчета облегчает выражение , скорость космического корабля относительно атмосферы, с точки зрения и :
Ну, а если под перифокальной системой понимать ту, что для начальной орбиты, с теми же направлениями осей (я не уверен, что именно так употребляется этот термин), то формула для та же. Единственная разница в том, что больше не указывает в направлении оси, вам нужно сначала вычислить его направление. Но, конечно, вы можете решить уравнения и в этой системе отсчета. Заметим, однако, что, поскольку и сопротивление не должно лежать в плоскости орбиты, плоскость может меняться со временем, что снижает полезность перифокальной рамки (она не будет перифокальной для более поздних орбит).
А если вы хотите, чтобы система координат, оси которой меняют направление, отслеживала ориентацию орбиты, уравнения усложнялись.
ооо
ооо