1 Определение и прогнозирование орбит низких спутников. Д.т.н. проф. А.И. Назаренко (НТЦ "Космонит"). 1. Введение 2. Усовершенствование методики определения начальных условий и прогноза движения 3. Определение и прогнозирование орбиты КА «Чибис» 4.Литература

2 1. Введение Во-первых, все модели атмосферы построены по измерительной информации, накопленной на большом предшествующем интервале времени. Поэтому по содержанию они являются моделями климата. Очевидно, что погода (и в нижней, и в верхней атмосфере) не совпадает с данными климата. Во-вторых, аэродинамические характеристики большинства спутников являются очень разными по величине и меняются во времени непредвиденным образом. В течение нескольких последних десятилетий точность прогнозирования движения низких спутников остается неизменной. А именно, СКО погрешностей прогноза составляют 10% от величины атмосферных возмущений. Этот факт объясняется двумя причинами.

3 МНК разработан 200 лет назад, когда искусственных спутников еще не было. Характерной особенностью движения околоземных ИСЗ является существенное влияние возмущающих факторов, оценка которых не поддается математическому описанию с необходимой точностью. Типичный пример такого рода возмущений – торможение в атмосфере. При использовании МНК влияние возмущающих факторов проявляется в необходимости выбора оптимального, т.н. мерного интервала. Погрешностей определения и прогнозирования орбит низких спутников обусловлены непредсказуемы- ми вариациями торможения на интервале обработки измерений и при прогнозе, а также невозможностью корректного учета этих вариаций в методе наименьших квадратов.

4 2. Усовершенствование методики определения начальных условий и прогноза движения Основы усовершенствованной методики опубликованы автором почти 40 лет назад. В 70-х годах эта методика была реализована в Российском центре космического пространства для определения и прогнозирования орбит низких спутников. В дальнейшем эта методика была усовершенствована. Характерной особенностью развитой методики является учет статистических характеристик атмосферных возмущений на интервале обработки измерений и при прогнозировании движения. Результаты изучения атмосферных возмущений опубликованы в ряде статей. Расчет статистических характеристик атмосферных возмущений при прогнозировании движения реализован в виде специального программного модуля. Принимается, что автокорреляционная функция атмосферных возмущений имеет вид

5 Вычисление матриц взаимной корреляции погрешностей прогнозирования вектора состояния на моменты времени (ti и tl) выполняется по формуле Задача оценки вектора состояния x (n 1) по измерениям Z (k 1) рассматривается в классической постановке. Учитывается возможность существования неких мешающих параметров q (m 1). В соответствии с ММП искомая оценка выражается следующим образом Значение мешающих параметров (шумов) вычисляется после построения оценки x на основе остаточных невязок с использованием соотношения вида где F - некоторая матрица.

6 Рассмотрены три подхода к оценке вектора состояния, отличающихся способом учета мешающих параметров (например, торможения в атмосфере): 1. Без учета мешающих параметров. В процессе оценки вектора состояния влияние мешающих параметров не учитывается. 2. Параметризация. Вектор мешающих параметров вводится в состав расширенного вектора состояния и затем применяется МНК. 3. Без параметризации (оптимальная фильтрация измерений, ОФМ). Установлены сравнительные соотношения между погрешностями оценок вектора состояния при использовании различных методов уточнения начальных условий по изменениям. Результаты анализа представлены на рисунке.

7 На основе метода ОФИ была разработана компьютерная программа для определения НУ по исходным TLE и прогнозирования движения низких спутников. Вектор состояния включает проекции радиуса-вектора (R) на оси инерциальной геоцентрической системы координат, соответствующие компоненты вектора скорости (V) в заданный момент времени t и оценку баллистического коэффициента Sb.

8 3. Определение и прогнозирование орбиты КА «Чибис» Данные о солнечной и геомагнитной активности в 1014 г

9 Оценки баллистического коэффициента изменяются, в основном, в соответствии с вариациями индексов солнечной и геомагнитной активности. Для метода ОФИ и МНК средние значения Sb соответственно составили и м 2/ кг. Для метода ОФИ и МНК СКО отклонений оценок Sb от среднего соответственно составили и м 2/ кг (20% и 36%).

10 Сравнение оценок Sb метода ОФИ и оценок Bstar Хотя оценки Sb и Bstar отличаются на порядок (по-видимому, они имеют разную размерность), их вариации является похожими: моменты времени максимумов и минимумов практически совпадают. Нормированные СКО от среднего равны: 20% для Sb и 30% для Bstar. Таким образом, из приведенных выше материалов сравнения можно сделать вывод, что применение метода ОФИ обеспечивает получение более корректных оценок баллистического коэффициента по сравнению с данными TLE и результатами применения МНК.

11

12 Определение времени существования КА Чибис Из этих данных видно, что оценки времени существования находятся в противофазе с колебаниями индекса F10.7. Существенное приближение расчетного времени падения по НУ за март и апрель объясняется непредвиденным увеличением солнечной активности в марте.

13 Заключение История разработки изложенного выше метода связана с именем П.Е. Эльясберга. Более 40 лет назад основные положения метода впервые были доложены и обсуждены на его семинарах в ИКИ и МГУ. В 1973 г П.Е опубликовал нашу большую статью на эту тему в сборнике «Прикладные задачи космической баллистики». Последующее развитие метода применительно к учету атмосферных возмущений, его испытания по модельной и реальной информации подтвердили реальную возможность существенного повышения точности определения и прогнозирования орбит низких спутников. Содержание рассматриваемого метода, учитывающего взаимосвязь результатов уточнения и прогнозирования параметров орбиты, позволяет расширить область его применения, а именно, использовать для прогнозирования движения более высоких спутников.

14 Литература 1. А.И. Назаренко, Л.Г. Маркова. Методы определения и прогнозирования орбит ИСЗ при наличии погрешностей в математическом описании движения. Сб. Прикладные задачи космической баллистики. Изд-во «НАУКА» Москва 1973, с.с A.I. Nazarenko. Determination and Prediction of Satellite Motion at the End of the Lifetime // International Workshop on Salyut-7/Kosmos-1686 Reentry, ESOC, Darmstadt (G), 9 April A.I. Nazarenko. Determination and Prediction of Orbits with Due Account of Disturbances as a «Color» Noise // AAS/AIAA Space Flight Mechanics Meeting. Monterey, CA, February AAS A.I. Nazarenko. Accuracy of Determination and Prediction Orbits in LEO // Estimation Errors Depending on Accuracy and Amount of Measurements, Seventh US/Russian Space Surveillance Workshop, Monterey, October-November A.I. Nazarenko, Accuracy of orbit determination and prediction for SOs in LEO. Dependence of estimate errors from accuracy and number of measurements. 26-th IADC, Moscow, (2008) 6. A.I. Nazarenko, Increasing the accuracy of orbit forecasting on the basis of improvement of statistical methods for processing measurements. Fifth European Conference on Space Debris, ESA SP-672, V.S. Yurasov, A.I. Nazarenko, P.J. Cefola and K.T. Alfriend, "Results and Issues of Atmospheric Density Correction, "Journal of Astronautical Society, Vol. 52, No. 3, July-September А.И. Назаренко. Применение метода оптимальной фильтрации измерений для уточнения и прогнозирования орбит космических аппаратов. Вестник, научно-технический журнал ФГУП НПО им. С.А. Лавочкина. 2, А.И. Назаренко. Погрешности прогнозирования движения спутников в гравитационном поле Земли, Москва, Институт космических исследований РАН, А.И. Назаренко, А. Г. Клименко. Патент на изобретение , Российская Федерация. Способ определения и прогнозирования движения космического аппарата на низких орбитах, подверженного влиянию торможения в атмосфере. ОАО «Российские космические системы». Заявка , приоритет изобретения Hoots, F.R. and R.L. Roehrich. Models for Propagation of NORAD Element Sets. Spacetrack Report No. 3, Aerospace Defense Command, United States Air Force, December 1980.