Разработка и применение алгоритмов численно-аналитического метода вычисления положений искусственных спутников Земли тема автореферата и диссертации по астрономии, 01.03.01 ВАК РФ

Чазов, Вадим Викторович АВТОР
доктора физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2012 ГОД ЗАЩИТЫ
   
01.03.01 КОД ВАК РФ
Диссертация по астрономии на тему «Разработка и применение алгоритмов численно-аналитического метода вычисления положений искусственных спутников Земли»
 
Автореферат диссертации на тему "Разработка и применение алгоритмов численно-аналитического метода вычисления положений искусственных спутников Земли"

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ имени М.В. ЛОМОНОСОВА

ГОСУДАРСТВЕННЫЙ АСТРОНОМИЧЕСКИЙ ИНСТИТУТ имени П. К. ШТЕРНБЕРГА

ЧАЗОВ Вадим Викторович

РАЗРАБОТКА И ПРИМЕНЕНИЕ АЛГОРИТМОВ ЧИСЛЕННО-АНАЛИТИЧЕСКОГО МЕТОДА

ВЫЧИСЛЕНИЯ ПОЛОЖЕНИЙ ИСКУССТВЕННЫХ СПУТНИКОВ ЗЕМЛИ

Специальность 01.03.01. — Астрометрия и небесная механика

Автореферат

диссертации на соискание учёной степени доктора физико-математических наук

1 4 ПАР 2013

Москва - 2012

005050708

005050708

Работа выполнена в отделе астрометрии и службы времени Государственного астрономического института имени П.К. Штернберга Московского государственного университета имени М.В. Ломоносова. Официальные оппоненты:

Защита состоится 4 апреля 2013 года в 14 часов на заседании Диссертационного совета Д501.001.86, созданного на базе Московского государственного университета имени М.В. Ломоносова.

Адрес: 119991, Москва, Университетский проспект, дом 13, ГАИШ МГУ.

С диссертацией можно ознакомиться в Научной библиотеке МГУ по адресу: Москва, Ломоносовский проспект, дом 27, Фундаментальная библиотека.

Автореферат разослан 4 марта 2013 года. Учёный секретарь

доктор физико-математических наук, профессор Московского авиационного института (Национального исследовательского университета)

МАРКОВ Юрий Георгиевич; доктор технических наук, профессор, заведующий отделом Института проблем информатики РАН

СИНИЦЫН Игорь Николаевич; доктор технических наук, профессор Московского государственного университета геодезии и картографии

ЯШКИН Станислав Николаевич.

Ведущая организация

Институт астрономии Российской академии наук.

Диссертационного совета

доктор физико-математических наук

С.О. Алексеев.

Общая характеристика работы

Актуальность темы

В диссертации представлена новая численно-аналитическая теория движения искусственных спутников Земли (ИСЗ). Теория была применена для определения параметров движения космических объектов и получения достоверных оценок геодинамических параметров на основе наблюдений.

Численно-аналитический подход позволяет использовать преимущества как аналитических методов, так и метода численного интегрирования уравнений движения.

Классическая форма аналитического способа вычисления положений небесных объектов заключается в следующей процедуре: на вековые изменения параметров орбиты накладываются долгопернодические неравенства и короткопериодические возмущения. Преимуществом аналитического метода является возможность увеличения скорости расчёта положений объекта при условии существенных ограничений на точность прогноза.

Суть метода численного интегрирования: на коротких интервалах времени параметры движения объекта аппроксимируют полиномами. Коэффициенты полиномов определяют путём вычисления значений правых частей дифференциальных уравнений и разностей этих значений в специальных точках внутри короткого интервала - одного шага интегрирования. Преимуществом численных методов является высокая точность вычислений. Малый шаг интегрирования, обеспечивая хорошую точность, требует значительных расходов вычислительного времени.

Численно-аналитический подход объединяет оба метода: часть неравенств, имеющих "короткий" период, определяется аналитически, а долгопернодические, резонансные и вековые слагаемые возмущающей функции составляют эволюционный гамильтониан - основу численного интегрирования "осреднён-ных" уравнений движения. Шаг интегрирования "осреднённых" уравнений может быть выбран достаточно большим.

Запуск Первого искусственного спутника Земли вызвал интерес учёных к задаче определения орбит объектов, возмущаемых аномалиями геопотенциала, притяжением Луны и Солнца, потенциалом, обусловленным приливами упругой Земли, сопротивлением атмосферы и световым давлением.

Профессором Е.П.Аксёновым была построена аналитическая теория движения ИСЗ на основе решения обобщённой задачи двух неподвижных центров. Результаты представлены в монографии [1).

Было открыто и получило признание обобщение одного из методов теории возмущений - метода канонических преобразований [5].

Тогда же был успешно применён численно-аналитический метод расчёта спутниковых орбит. Важные результаты в этом направлении получены профессором М.Л.Лидовым и его учениками, некоторые итоги подведены в обзорном докладе [3], там же намечены перспективы исследований.

Тем не менее, в наши дни все центры обработки высокоточных наблюдений ИСЗ проводят расчёты с помощью программ, в которых прогноз положений объектов выполняется алгоритмом численного интегрирования дифференциальных уравнений движения [9].

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

В качестве исходных материалов предлагаемого исследования были использованы:

рекомендации Международного астрономического союза, стандартные соглашения Международной службы вращения Земли [6], ■ лазерные и позиционные наблюдения космических объектов [7], база данных о точных положениях навигационных спутников [8], база данных об элементах орбит объектов [10].

Цель работы

Целью исследований является решение следующих задач:

вывод формул для вычисления параметров промежуточной орбиты, основанной на решении обобщённой задачи двух неподвижных центров, с точностью, ограниченной только возможностями компьютера;

разработка метода преобразования возмущающей функции на основе параметров промежуточной орбиты;

дифференцирование и интегрирование слагаемых возмущающей функции;

составление и численное интегрирование "осреднённых" уравнений движения и учёт короткопериодических неравенств с использованием "несингулярных" элементов орбиты;

обработка высокоточных лазерных измерений топоцентрических дальностей до спутников Лагеос и Лагеос-2 на длительных интервалах времени и оценка значений геодинамических параметров;

разработка методики предсказания ситуаций сближения ИСЗ;

фильтрация позиционных наблюдений и оценка значения отношения средней площади поверхности к массе спутника.

Стандартные соглашения Международной службы вращения Земли содержат рекомендации по обработке наблюдений искусственных спутников Земли с помощью метода численного интегрирования. Разработка методики аналитического решения задачи с учётом всех рекомендаций астрономических организаций также является целью исследования.

В тексте диссертации термин алгоритм объединяет несколько понятий: это и связанные между собой формулы и соотношения, и последовательность действий, и'реализация процедуры вычислений на компьютере. Совокупность алгоритмов, предназначенная для решения поставленной задачи, в тексте называется "программным обеспечением" или "пакетом программ".

Научная новизна

Научная новизна работы заключается:

в совокупности формул и соотношений для вычислений в рамках обобщённой задачи двух неподвижных центров с максимально возможной методической точностью;

в методе преобразования возмущающей функции задачи, заключающейся в построении всё более сложных конструкций на основе простых начальных соотношений;

в способе вычислений, позволяющим учитывать возмущающее действие факторов различного происхождения одним набором формул;

в совокупности алгоритмов построения моделей движения космических объектов с помощью численно-аналитического метода;

в методике расчёта ситуаций опасных сближений объектов.

Достоверность результатов

Достоверность полученных результатов подтверждается примерами обработки наблюдений различных объектов., сравнением с данными, предоставляемыми Международной службой вращения Земли, и сопоставлением с параметрами движения ИСЗ, публикуемыми в Интернете.

Практическое значение

Практическая ценность диссертации определяется тем. что:

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

алгоритмы применяются для вычисления целеуказаний на Звенигородской научной базе ИНАСАН и филиале ИНАСАН на пике Терскол;

алгоритмы использовались для планирования космических экспериментов и предварительной редукции результатов наблюдений спутника Метеор-ЗМ сотрудниками Центральной аэрологической обсерватории ГМЦ;

было подготовлено учебное пособие "Модель движения ИСЗ", являющееся реконструкцией и расширением монографии профессора Е.П.Аксёнова [1].

Положения, выносимые на защиту На защиту выносятся следующие положения.

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

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

3. Впервые разработан метод аналитического интегрирования слагаемых возмущающей функции, зависящих от позиционных и угловых переменных промежуточной орбиты. Представлен способ вычисления правых частей при численном интегрировании "осреднённых" уравнений движения в "несингулярных" элементах орбиты. Методы используют алгоритмы точного вычисления как параметров промежуточной орбиты, так и частных производных от любых параметров промежуточной орбиты по позиционным и угловым переменным, а также по каноническим переменным действие-угол.

4. С помощью тензорного преобразования получены приближённые формулы связи между координатами небесных тел Солнечной системы, полученными при использовании "изотропных" координатных условий с одной стороны, и "гармонических" координатных условий с другой стороны. Общий вывод состоит в следующей рекомендации: в прикладных задачах достаточно записать релятивистские уравнения движения искусственного спутника Земли на основе "гармонических" координатных условий, а при вычислении возмущающих сил использовать современные численные эфемериды Солнца, Луны и планет, полученные в метрике с "изотропными" координатными условиями.

5. Выражения для учёта возмущений от приливов упругой Земли и океанических приливов, записанные во вращающейся системе отсчёта и рекомендованные Международной службой вращения Земли, преобразованы к системе отсчёта, связанной с истинным экватором даты.

6. В результате анализа массива высокоточных лазерных измерений то-поцентрнческих дальностей до ИСЗ Лагеос и Лагеос-2 получены оценки параметров вращения Земли и скоростей смещения измерительных пунктов в Земной системе отсчёта. Для некоторых объектов на основе наблюдений определены оценки отношения средней площади поверхности к массе спутника.

7. Разработана методика предсказания ситуаций сближения объектов искусственного происхождения в околоземном пространстве. Методика была применеЕ1а на этапе проектирования параметров орбиты метеорологического "стационарного" спутника.

Публикации в журналах и сборниках трудов

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

1. Герасимов И.А., Чазов В.В., Рыхлова Л.В., Тагаева Д.А. Построение теории движения тел Солнечной системы, основанной на универсальном методе вычисления возмущающей функции. //Астрономический вестник. 2000. Т.34. Номер 6. С.559-566.

2. Герасимов И.А., Чазов В.В., Тагаева Д.А. Применение универсального метода вычисления возмущающей функции в численно-аналитической теории движения малых планет. //Вестник Московского университета. Серия 3. Физика. Астрономия. 2000. Номер 3. С.55-57.

3. Бахтигараев Н.С., Чазов В.В. Информационное обеспечение космических экспериментов на основе численно-аналитической теории движения искусственных спутников Земли. //Космические исследования. 2005. Т.43. Номер 5. С.386-389.

4. Чазов В.В., Герасимов И.А.; Соловьёва О.Д. Изотропные и гармонические координатные условия в пространстве-времени Солнечной системы. //Вестник Московского университета. Серия 3. Физика. Астрономия. 200G. Номер 2. C.6G-68.

5. Клишин А.Ф.. Чазов В.В., Бахтигараев Н.С., Костюк Н.Д. Об оценке

' уровня техногенной опасности в зоне размещения КА "Электро-Л".

//Вопросы радиоэлектроники. Серия "Радиолокационная техника". 2007. Выпуск 2. С.40-46.

6. Чазов В.В.. Бахтигараев Н.С., Косткж Н.Д. Наблюдения спутника "Молния 3-39" в Звенигородской обсерватории ИНАСАН и определение времени падения. //Вестник СибГАУ. 2011. Выпуск 6(39). С.183-185.

7. Бахтигараев Н.С., Лёвкина П.А., Сергеев A.B., Чазов В.В. Наблюдения неизвестного фрагмента космического мусора в Терскольской обсерватории. //Вестник СибГАУ. 2011. Выпуск 6(39). С.186-189.

8. Гаипова А.Н.. Чазов В.В. Комплекс программ Лепта. //Измерительная техника. 1991. Т.6. С.30-30.

9. Герасимов И.А., Чазов В.В. Переменные действие-угол в обобщённой задаче двух неподвижных центров. //Труды ГАИШ. 1988. Т.59. С. 4652.

10. Чазов В.В. Основные алгоритмы численно-аналитической теории движения искусственных спутников Земли. //Труды ГАИШ. 2000. Т.68. С. 5-20.

11. Чазов В.В. Создание численно-аналитической теории движения небесных тел. /Труды конференции "Околоземная астрономия - 2003'', Терскол. 2003. С. 171-175.

12. Бахтигараев Н.С., Чазов В.В. Компьютерное моделирование условий наблюдений небесных тел. //Кинематика и физика небесных тел. 2003. Номер 4. С. 105-107.

13. Бахтигараев Н.С., Чазов В.В. Моделирование движения космических аппаратов с учётом рекомендаций Международного астрономического союза. /Труды конференции "Околоземная астрономия - 2005", Казань. 2005. С.281-285.

Апробация результатов

Основные результаты, полученные в диссертации, докладывались автором и соавторами на семинарах и конференциях:

• Координационный совет по небесной механике и Координационный совет по астрометрии ГАИШ МГУ;

• семинар ИНАСАН "Проблемы происхождения и эволюции кометно-астероидного вещества в Солнечной системе и астероидная опасность" ;

• конференция "Астрометрия, геодинамика и небесная механика на пороге XXI века", Санкт-Петербург, 2000. год ; тезисы доклада "О новой теории движения тел Солнечной системы, основанной на универсальном методе вычисления пертурбационной функции" опубликованы в сборнике на стр.249 (совместно с И.А.Герасимовым и Д.А.Тагаевой);

• конференция "АСТРОЭКО-2002: Состояние и перспективы международных исследований по наблюдательной астрономии, экологии и экстремальной физиологии в Приэльбрусье", Терскол, 2002 год;

• конференция "Околоземная астрономия - 2003", Терскол, 2003 год;

• Всероссийская астрономическая конференция "Горизонты Вселенной", Москва, 2004 год; тезисы доклада "104 года стандарта вычислений астрометрии" опубликованы в сборнике "Труды ГАИШ", 2004, том 75, стр.216 (совместно с И.А.Герасимовым);

• конференция "Околоземная астрономия - 2005", Казань, 2005 год.

Структура диссертации

В изложении на бумаге вся завершённая работа занимает три главы и 210 страниц с предисловием, послесловием и приложением. Текст содержит 30 таблиц, 34 рисунка и список литературы из 208 наименований.

Содержание работы

В первом разделе диссертации выполнена постановка задачи.

Уделено внимание рекомендациям Международного астрономического союза об использовании числовых значений астрономических постоянных и о системах отсчёта пространства-времени и рекомендациям Международной службы вращения Земли по составу ускорений, возмущающих движение спутника. Учёт упомянутых рекомендаций необходим для того, чтобы сделать результаты обработки различных экспериментов сопоставимыми.

Подчёркивается необходимость использования неинерциалъной системы отсчёта, связанной с истинным экватором даты.

Величина угла наклонения спутника, вычисляемая относительно неподвижного экватора геоцентрической небесной опорной системы отсчёта, имеет вековую составляющую. Этот эффект, являющийся следствием прецессии оси вращения Земли, не влияет на алгоритмы численного интегрирования, но создаёт дополнительные трудности для аналитических исследований.

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

Выбранная неинерциальная система вращается относительно небесной опорной системы отсчёта с угловой скоростью . С точностью до первых степеней малых величин компоненты вектора Г2(£) равны

(1)

При этих ограничениях уравнения движения имеют вид

¿Г _ Ж £/р _дК Л ~ +др' <И ~~ дг*

где г - вектор в подвижной системе отсчёта, р = г+[Пхг] - обобщённый импульс, К - гамильтониан задачи,

К=^|р,|2-И^-Й.[гхр]-Д. (3)

№ - силовая функция обобщённой задачи двух неподвижных центров, Я -возмущающая функция второго порядка малости относительно сжатия.

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

(4)

где К = — а\ + К2 - гамильтониан системы, К2 = й-[г х p]+R - возмущающий гамильтониан второго порядка малости относительно сжатия. В случае отсутствия возмущений К2 = 0, и система имеет решение

Щ = и (i0), l{t) = i(t-t0) + lQ(t0),

G(t) = G0(t0), g(t) = g(t-to)+go(to), (5)

H{t) = Ho(to), h{t) = h(t-t0) + h0(t0), причём скорости изменения угловых переменных

; п0 . пйи n0fi

являются постоянными величинами.

Отметим следующий факт: необходимо знать только частные производные от гамильтониана /sT по каноническим элементам, представлять же гамильтониан в виде явной функции переменных действия L, G, Н нет никакой необходимости.

dL dK + di ■ dl Ж

dt dt ~ dL '

dO dK dg dy _ dK

~dt dt ~ 3G '

dH 8K + dh' dh dK

dt dt ~~ дН '

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

Выполним переход от оскулирующих элементов Ь, О, Н, I, д, И к новым переменным V, С, Н', V, д', Ы. Функцию преобразования и вековой гамильтониан обозначим, соответственно, 52 и . Эти функции, как и возмущающий гамильтониан К2, имеют второй порядок малости относительно сжатия.

В новом гамильтониане Щ (V, С, Н', I', д', И') оставим те слагаемые возмущающей функции К2, период изменения которых превышает несколько суток. Функцию преобразования 5'а найдём по формуле

Яг = ¡(.К2-К1)йг. (7)

С точностью до второго порядка относительно малого параметра, сжатия Земли, выражения 52, К2. Кд = -ах, К2 можно считать функциями как оскулирующих элементов Ь, С, Я, I, д, Л, так и новых переменных и, С", Я', V, д', Л'.

Новые переменные могут быть названы "средними" или "сглаженными" элементами орбиты. В изменениях "сглаженных" параметров отсутствуют вариации с малым периодом. Между двумя наборами переменных с точностью до второго порядка малости существует простая связь:

Ь - Ь + 31' - дЬ"

Уравнения движения, записанные в новых переменных, назовём "осред-нёнными" уравнениями. Они сохраняют каноническую форму и имеют вид:

¿V _ дК\ Ж _ дох _ сЩ ~ л " а и эй'

_ I ' (9)

<Й дд> ' Л ас ЭС '

¿Н' _ (Щ М = даг дЩ

~ШГ ~~ + ан ' м +эн' он'' - 12-

Материал второго раздела составляют основные алгоритмы модели движения ИСЗ в численно-аналитической форме.

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

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

Для описания движения искусственных спутников приняты во внимание возмущения, обусловленные гравитационным полем Земли, притяжением Луны и Солнца, приливными деформациями упругой Земли, океаническими приливами и давлением солнечной радиации.

Показано, что возмущающие функции, обусловленные различными эффектами, зависят от пяти "начальных'' функций координат объекта:

г г0 г х у Го ' г ' г ' г ' г

Далее выражения для "начальных" функций записаны при помощи соотношений промежуточной орбиты, основанной на решении обобщённой задачи двух неподвижных центров.

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

Пусть

а, е, ё, в , а - параметры промежуточной орбиты спутника, ф, Е, П - угловые переменные промежуточной орбиты, £ф - гринвичское истинное звёздное время, 1м ■. Рм ■> В, Пм ~ фундаментальные аргументы.

Пусть известны также приближённые численные значения позиционных параметров орбиты объекта ао, ео, ¿о и знак параметра аз, зависящий от величины угла наклонения орбиты.

С помощью значений а0, е0, ¿о образуем одномерные массивы, содержащие численные значения как параметров промежуточной орбиты, так и частных производных высших порядков по параметрам а, е, S.

Определим элементарное слагаемое как структуру вида

■д = к\ф + к2ф 4- faE 4- /с4П + faS@ + falM + k7l's + ksFM + faD + к10Пм, (11) где А(о0, ео, 50) - одномерный массив действительных чисел, полученный на основе начальных значений ао, ео ■. ¿о 1

показатель степени j\ может принимать положительные и отрицательные

. (min) - ■ ^ ■ (max)

целые значения и нуль и удовлетворяет условию ]г ^ j\ ^ зг , показатели степени j2, h Л ■. h , Л могут принимать только положитель-

п ^ ■ ^ -{max)

ные целые значения и нуль и удовлетворяют условиям U ^ jk ^ jk , коэффициенты fa, fa, fa, fa, fa, fa, fa, fa, fa, fa0 принимают любые положительные и отрицательные целые значения и нуль.

Отличительной чертой алгоритма является способ конструирования возмущающей функции на основе "начальных" функций координат. В предлагаемом алгоритме формулы промежуточного движения учитываются с максимально возможной точностью.

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

Выполнено преобразование осреднённых уравнений движения в систему "несингулярных" элементов орбиты:

a', l' + g' + h', e'cos(g' + h'), e'sin{g' + h'), 6'cosh!, 5'smh'. (12) Штрихи означают, что все величины являются "средними".

Преобразуем "осреднёаные" уравнения к виду

да' дЩ да!_дЩ да' дЮ2 + dL' QV да\ да\ да\

+ ~Ш7 + dG'+ dip dig _ dtq _ эк}

dL' ~ ~дО ~ Тш7' de'

+ — cos (g' + h)

(дач , dax dig ОК*Л

-esino+n ----------- /,q\

* \dG' dll' дG' dll' J1 (13)

de' . , , .... + -sm(g' + k')

w '\dG' dH' dG' dH'J'

J.dS' w v ■ u, (dcyi dKi\ + —r- COS h. - 5' sin h---- ,

dt \dlf dH'J'

+ -7— sill h + 6' cos h' —---- .

dt \dH' dH'J

Производные от параметров e' и 6' по времени t равны

= 1 /де'2 dig де^дЩ де'2 Oig\

dt 2e'\dL' dl' + dG' dg' + 'dfP~dhFj'

d£_ = 1 fos'2 dig ds^dig ds^dig\ dt 26'\dL> dl' + dG' dg' + dlr ~7)J7J '

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

Для заданных приближённых численных значений позиционных параметров орбиты'а0, е0, Jo и знака параметра a:i найдём слагаемые возмущающей функции и получим ряды для функции преобразования S2 и "осреднённого" гамильтониана .

da' dt

d(l' + g' + h') dt

d [ef cos(g' + h')\ dt

d [e'sin(ff' H-dt

d[5'cos/i'] dt

¿[J'sin/i'] dt

Пусть в начальный момент времени t0 известны численные значения "средних" элементов орбиты a'(ío); e'(í0), S'(t0), l'{to), g'{to), h'(to) ■

Вычислим числовые значения "несингулярных" переменных (12) в момент времени £о • С такими начальными условиями на заданном интервале времени выполним численное интегрирование "осреднённых" уравнений (13).

Для преобразований от системы сглаженных "несингулярных" элементов к "средним" параметрам орбиты e'(t), S'(t), l'(t); g'{t), h!(t) в произвольный момент времени £ будем использовать формулы

e'(í) = \J [e'cos(</ + h')]2 + [e' sin (г/ + h')}2 , S'(t) = yJ(S'cosh')2 + (ó'smh')\ «*Ш + НШ = , eos k'(t) = ^ , (15)

sin{9'(t)+m) = e's%;h\ =

l'(t) = (¿' + g' + h!) - (g'(t) + h'(t)), g'(t) = (g'(t) + h'(t)) - h'(t).

На основе "средних" параметров орбиты на любой момент времени определим числовые значения оскулирующих "несингулярных" элементов (8):

a, l + g + h, ecos(g + h), esin(g + h), ícos/i, ísin/i.

Спекулирующие элементы промежуточной орбиты е, ó, 1, д, h вычислим из соотношений, аналогичных формулам (15).

Далее с помощью формул, основанных на решении обобщённой задачи двух неподвижных центров, и числовых значений оскулирующих элементов промежуточной орбиты получим мгновенные значения прямоугольных координат и скоростей объекта.

В ходе построения численно-аналитической модели движения сделано несколько упрощающих предположений:

• аномальная часть геопотенциала и разложение потенциала, обусловленного притяжением Луны и Солнца, были ограничены по порядку и степени сферических функций;

• были использованы приближённые тригонометрические ряды для представления функций от координат возмущающих тел;

• при построении возмущающей функции в "буквенном" виде ограничена величина максимальной степени позиционных параметров;

• при интегрировании предполагается, что звёздное время и фундаментальные аргументы являются линейными функциями времени;

• алгоритм аналитического интегрирования не учитывает короткоперио-дические вариации позиционных параметров.

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

Далее выполнены дополнительные исследования.

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

Международный астрономический союз, принимая во внимание факт, что многие работы по теории относительности выполнены при использовании "гармонических'' координат, оказавшихся полезными для приложений, рекомендует выбор гармонических координатных условий.

С помощью тензорного преобразования

9а0(х°,= ¿СЛЛ **) (16)

были выведены формулы связи между координатами х'°, х'1, ха, х'г, удовлетворяющими "изотропным" координатным условиям, и "гармоническими" координатами х°, х1, х2 , х3 .

Расчёты показали, что отличия в численных значениях координат очень малы. Для пространственно-временной траектории Земли, например, они не

превосходят 3 метра на интервале 100 лет.

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

Выражения для учёта возмущений от приливов упругой Земли и океанических приливов, записанные во вращающейся системе отсчёта и рекомендованные Международной службой вращения Земли (геофизический подход), были преобразованы к системе отсчёта, связанной с истинным экватором даты (астрономический подход).

Для согласования численной и аналитической моделей возмущающих сил, обусловленных действием Луны и Солнца, использована методика разностно-полиномиальной коррекции.

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

Учтём сжатие Земли и запишем дополнительную часть возмущающей функции, обусловленной действием Луны на Землю:

* = £)'«(£). от

где

/ - гравитационная постоянная, т' - масса Луны,

г о - экваториальный радиус Земли,

/2 - коэффициент при второй зональной гармонике,

Р2{х) - полином Лежандра второго порядка.

г' - расстояние между Землёй и Луной,

х'. у' у г' - координаты Земли относительно Луны.

Дифференцируя потенциал и по переменным х', у', г', определим составляющие ускорения V/ центра масс Земли, обусловленного взаимодействием Луны и сжатия Земли. Для того, чтобы в выражениях для ускорения

перейти к координатам Луны относительно Земли, следует знак каждой из переменных х', у', г' поменять на противоположный:

(го12

гп

ох г

ди / Ч 1 — = /ш Л -75 ау ги

дг г-

15

У

15 2

15

У

(18)

Неинерциальность системы отсчёта, выбранной для изучения движения спутника, обусловлена не только вращением относительно небесной опорной системы отсчёта с угловой скоростью {!(<) (формула (1) на с.10), но и ускорением Ж, вызываемым взаимодействием Луны и сжатия Земли. Гамильтониан задачи К (с. 11) надо дополнить ещё одним слагаемым - скалярным произведением ускорения ]¥ на вектор положения объекта г:

1

К=2|Я

IV + (Ф ■ г) - • [г х р ] Возмущающая функция

% = - • г) = - х\Ух - у\Уу

Л.

(19)

имеет вид

= +

+

/ш\ / г Го ) \Го

/т \ / г

соэ гу 4- бш ю

(20)

где

А

("О 10

/т /т

Действия с функцией (20) выполняются с помощью тех же алгоритмов, которые были использованы для преобразования возмущающих функций, обусловленных притяжением Луны и Солнца.

В третьем разделе приводятся различные примеры применения алгоритмов для обработки длительных рядов высокоточных наблюдений ИСЗ.

Выполнена проверка рекуррентных соотношений для вычисления полиномов Лежандра высших порядков

Рп(х) = ■ X ■ Р^х) - Hzl. рп^х) (21)

71

с начальными условиями

-1 < х ^ +1, Poix) = 1,

Pt(x)=x, (22)

ад = 4+§,*'.

Возможны два варианта применения рекуррентной формулы.

1. Вычисление значения функции Рп{х) для произвольного п при заданном числовом значении аргумента х (численное интегрирование).

2. Нахождение численных значений коэффициентов полинома Рп(х) при различных степенях аргумента х (аналитический подход).

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

Первый вариант может быть проверен с помощью тождеств

Р„(-1) = (-1)", Р„(+1) = +1. (23)

Расчёты показали, что в результате использования формулы (21) с начальными условиями (22) и аргументами х = -1 и х = +1 тождества (23) справедливы с точностью до 15 знаков после запятой при всех порядках полинома Лежандра от п = 1 до п = 720.

•Второй способ проверки состоит в использовании формулы для производящей функции полиномов Лежандра:

2 оо

Расчёты показали (табл.1), что при различных численных значениях параметра а и аргумента х значение суммы в правой части выражения (24) соответствует значению функции в левой части с точностью до 12 значащих

Таблица 1: Производящая функция для полиномов Лежандра

а X функция разность V

0.96 +1.00 25.000000000000000 -0.988 ю- 12 755

0.96 +0.99 6.9337524528153640 -0.196 ю- 12 606

0.96 +0.95 3.2009219983223993 -0.090 ю- 12 546

0.96 +0.90 2.2727272727272727 +0.685 10" 12 565

0.96 +0.80 1.6103915660020771 -0.926 10" 12 547

0.96 +0.50 1.0197712705600052 -0.897 10" 12 561

0.96 +0.20 0.8064516129032258 -0.269 10- 12 502

0.96 +0.00 0.7213873210309515 +0.978 10- 12 576

0.96 -0.20 0.6585792122172903 -0.972 10" 12 471

0.96 -0.50 0.5890920370328413 -0.907 10" 12 550

0.96 -0.80 0.5377898796468977 +0.210 10- 12 525

0.96 -1.00 0.5102040816326531 +0.975 10" 12 660

цифр. Целое число Ытах в последней колонке означает наибольший порядок полинома Рп(х) в сумме (24). Значение а = 0.96 характерно для объектов с высотой полёта 300 километров над поверхностью Земли.

Для параметра а = 0.90, присущего геодезическим спутникам с высотой полёта более 700 километров, значение Мтах < 300.

Вывод: для заданных значений аргумента численные зна-

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

Второй вариант используется на предварительной стадии преобразования возмущающей функции.

Пусть а-"' - численные коэффициенты полинома Рп(х) = а^ ■ хг:

¿=0

(п) 2 П — 1 (п—1) П-1 (п—2)

а) ---а■--- а]

1 п 1 п

„(о) _ т . я(0) _ п „(1) ^ , . „(0) = а(1) = 0 а(2) = 3

2 ' 2 2

Сумма коэффициентов полинома любого порядка всегда равна единице

Е4п) = +1,

г=0

но величины коэффициентов достигают больших значений.

В табл.2 приводятся числовые значения некоторых коэффициентов:

Таблица 2: Величины коэффициентов полиномов

п г а; '

23 1 -0.1681880950927734

23 17 -18733148.9151477814000000

23 23 981501.4276027679440000

26 0 -0.1549810171127319

26 18 229350982.8986495730000000

26 26 7389761.9984761476500000

29 1 4.3339334428310394

29 21 2978545001.5916727500000000

29 29 56004648.0960045755000000

36 0 0.1320605995715596

36 26 -1150785784265.5309200000000000

36 36 6439404973.5978009900000000

45 1 5.3821880423697621

45 33 2523815359824958.2900000000000000

45 45 2950952798742.9865800000000000

50 0 -0.1122751726592170

50 36 -189971831414987581.0000000000000000

50 50 89609514959900.0547000000000000

Большое количество нулей после восемнадцатой значащей цифры каждого числа возникает как следствие ограниченности разрядной сетки компьютера.

Вывод: алгоритм (25) определения значений коэффициентов полинома Рп{х) при различных степенях аргумента х приводит к потере вычислительной точности в случае полиномов высоких порядков. Для п = 36 погрешность может оказаться в пятом знаке после запятой. При значениях п > 39 нет смысла использовать полученные коэффициент.

Важное место в процессе разработки и отладки программного обеспечения занимает методика сравнительных испытаний, направленная на определение точностных характеристик расчётов.

Одним из вариантов методики является сравнение результатов прогноза движения с результатами высокоточных измерений топоцентрических дальностей до искусственных спутников Земли, снабжённых уголковыми отражателями. Такие измерения проводятся на обсерваториях с помощью лазерных дальномеров.

Вариант заключается в применении алгоритма дифференциального улучшения орбит на интервалах времени продолжительностью несколько суток. Результатом вычислений являются оценки параметров модели и "невязки" - разности между измеренными и вычисленными значениями. Полученные "невязки" необходимо сравнить с результатами аналогичных вычислений, регулярно выполняемых в центрах обработки данных Международной службы лазерной далыюметрии. В табл.3 даны названия некоторых центров.

Таблица 3: Центры анализа

SAO DGFI HIT-U JCET МСС CSR DELFT

астрономическая обсерватория, Шанхай (Китай) Институт геодезических исследований (Германия) университет Хитоцубаши, Токио (Япония) отдел HACA и университет Мэриленд (США) Центр управления полётами (Россия) Центр космических исследований, Техас (США) Технологический университет, Дельфг, Нидерланды

Ежедневно, по мере поступления новых массивов измерений топоцентрических дальностей, выполняется улучшение начальных параметров движения объектов и ряда эмпирических параметров. Для всех измерительных пунктов на каждом прохождении вычисляются "невязки". В результате линейной аппроксимации таких "невязок" получают оценку "среднего смещения" (СС) дальностей, измеренных на отдельном прохождении спутника в поле видимости обсерватории, и оценку "среднего квадратического отклоне-

ния" (СКО) от аппроксимирующей прямой.

Таблица 4: Объекты и наблюдения

Спутник Ь'т'т е г п й> й %

9306102 Стелла 800 0.0007 98.476 14.27121 -2.9342 0.9695 3000

7501001 Старлет 900 0.0206 49.817 13.82205 3.3041 -3.9476 4000

8606101 Эйджисаи 1500 0.0011 50.010 12.44446 2.5461 -3.0752 8000

7603901 Лагеос 6000 0.0044 109.835 6.38665 -0.2137 0.3425 5000

9207002 Лагеос-2 6000 0.0137 52.650 6.47295 0.4370 -0.6316 5000

8900103 Эталон-1 19500 0.0012 65.312 2.13156 -0.0062 -0.0320 600

8903903 Эталон-2 19500 0.0011 64.778 2.13205 -0.0036 -0.0332 600

В табл.4 для нескольких спутников представлены средние значения минимальной высоты полёта /¿тгп в километрах, эксцентриситета орбиты е, угла наклонения орбиты г в градусах, среднего движения п в оборотах за сутки, скорости изменения аргумента перигея ш в градусах за сутки, скорости изменения долготы восходящего узла П в градусах за сутки и количество "нормальных мест" ЛГр, обычно получаемых на интервалах наблюдений протяжённостью один месяц.

Наблюдения двух спутников Лагеос поставляют основной материал для вывода осреднённых на коротких интервалах времени параметров вращения Земли и координат измерительных пунктов.

Совокупности измеренных дальностей до двух спутников Эталон недостаточно для решения этих задач в полном объёме. Обработка наблюдений этих объектов, обращающихся на орбитах спутников системы Глонасс, необходима для калибровки измерительной аппаратуры и отладки пакетов вычислительных программ.

Геодезические спутники Эйджисаи, Старлет, Стелла обращаются вокруг Земли на низких орбитах. Обработка наблюдений этих объектов затруднена необходимостью учёта гармоник геопотенциала большого порядка и непредсказуемыми вариациями плотности верхней атмосферы Земли.

Результаты расчётов, относящиеся к спутникам Лагеос и Лагеос-2, соби-

рает и публикует Европейский центр определения орбит (CODE), расположенный в городе Берн (Швейцария).

Для проведения сравнительных испытаний была выполнена фильтрация лазерных наблюдений с помощью пакета программ LENTA, разработанного на основе предлагаемых в данном исследовании алгоритмов.

В табл.5 содержатся результаты определения "среднего смещения" и "среднего квадратического отклонения" для измерений топоцентрической дальности до спутника Лагеос, выполненных на пунктах 7080 (Форт Дэвис,

Таблица 5: Лагеос, сравнительные испытания, пункты 7080 и 7105

Оценка "среднего смещения" в миллиметрах, пункт 7080

дата Np LENTA DGFI МСС HIT-U JCET SAO

05.10.2009 14 27 ±3 9 ± 3 —7 ± 3 19 ±3 35 ±4 22 ±3

09.10.2009 8 3± 10 -7± 14 -30 ± 11 -1 ±7 -21 ± 12 4 ± 8

10.10.2009 8 30 ±3 -19 ±3 —19 ± 4 —14 ± 3 23 ±6 5 ±3

11.10.2009 15 -12 ±9 —18 ± 11 -27 ±5 -6 ±3 -12 ± 8 -2 ±4

Оценка "среднего смещения" в миллиметрах, пункт 7105

дата Np LENTA DGFI МСС HIT-и JCET SAO

05.10.2009 7 9 ± 5 9 ± 6 5 ±3 6±4 11±8 13 ±6

06.10.2009 19 — 11 ± 9 21 ±4 -60 ±3 2 ± 2 -16±5 -10±3

07.10.2009 16 -34 ±4 — 14 ± 3 -16±3 —13 ± 4 —3 ± 5 -23 ±1

07.10.2009 9 -12 ±5 -4 ±5 1 ±3 -7 ±4 -18 ±6 -9 ±2

07.10.2009 12 58 ±3 -7 ±4 3± 1 5 ± 2 -3±8 -5 ±3

08.10.2009 6 31 ±6 —5 ± 3 15 ±2 5 ± 2 31 ± 1 4 ± 2

09.10.2009 4 -23 ± 1 -13 ±3 -3±3 -11 ±2 -26 ±4 -22 ±1

09.10.2009 13 -12 ±5 —3 ± 5 —5 ±4 -7 ±2 -29 ±7 -13±3

09.10.2009 18 40 ±4 -15 ±4 6 ±3 —3 ± 2 -2 ±5 -4±2

США) и 7105 (Вашингтон, США). В таблице оценки "СС" и "СКО" даны в миллиметрах. В колонке Np дано число "нормальных точек". Данные по пяти центрам обработки были взяты из отчёта CODE, данные в колонке LENTA получены в процессе вычислений.

Анализ позволяет сделать вывод, что пакет вычислительных программ LENTA уступает по точности представления наблюдений спутников типа Ла-

геос пакетам DGFI, МСС и HIT-U не более 20%. Точность пакетов LENTA, JCET и SAO при обработке высокоточных лазерных наблюдений Лагеос и Лагеос-2 находится примерно на одном уровне.

Результаты сравнения (табл.6) приводят к выводу, что пакет вычислительных программ LENTA уступает по точности представления наблюдений спутников типа Эталон и низкоорбитальных объектов пакетам DGFI, JCET и HIT-U. В колонке At дана продолжительность интервала измерений в минутах времени, оценки "СС" и "СКО'" даны в миллиметрах.

Таблица 6: Обсерватория Яррагади, сравнительные испытания

спутник дата UTC At Np LENTA DGFI

Эталон-1 2010/02/21 11:59 221.4 И 379.5 ± 10.3 16 ±31

Эталон-1 2010/02/26 17:47 150.4 15 -23.0 ± 5.6 -28 ±21

спутник дата UTC At NP LENTA HIT-U

Эйджисаи 2010/01/21 10:47 13.6 21 272.6 ± 159.1 43 ±5

Эйджисаи 2010/01/21 12:48 14.8 31 384.7 ±429.5 1 ± 4

Эйджисаи 2010/01/21 17:00 9.1 20 -1274.3 ±491.4 5 ± 5

Эйджисаи 2010/01/21 21:07 10.0 20 -3639.5 ± 896.3 25 ±6

Старлет 2010/03/01 07:13 2.6 6 1047.6 ± 9.5 22 ±6

Старлет 2010/03/01 09:00 6.6 15 -717.3 ±59.6 24 ±5

Старлет 2010/03/01 14:33 4.3 10 -1143.0 ±20.4 11 ±1

Старлет 2010/03/01 16:21 9.8 21 -929.5 ±83.4 —0 ± 4

Старлет 2010/03/01 18:11 6.6 15 636.5 ± 60.3 -8 ±4

Стелла 2010/01/21 06:08 7.8 17 -1762.9 ±3823.5 —2 ± 3

Стелла 2010/01/21 17:22 7.2 16 5780.8 ±987.1 —3 ± 6

спутник дата UTC At LENTA JCET

Эйджисаи 2010/01/21 10:47 13.6 21 272.6 ±159.1 -16 ±16

Эйджисаи 2010/01/21 12:48 14.8 31 384.7 ±429.5 74 ±13

Эйджисаи 2010/01/21 14:58 4.0 9 -850.0 ±5.3 195 ± 24

Эйджисаи 2010/01/21 17:00 9.1 20 -1274.3 ±491.4 154 ± 21

Эйджисаи 2010/01/21 19:12 3.9 9 -618.6 ±20.8 321 ±39

Эйджисаи 2010/01/21 21:07 10.0 20 -3639.5 ±896.3 59 ±19

Старлет 2010/03/01 09:00 6.6 15 -717.3 ±59.6 267 ± 131

Старлет 2010/03/01 14:33 4.3 10 -1143.0 ±20.4 348 ± 232

Старлет 2010/03/01 16:21 9.8 21 -929.5 ±83.4 126 ± 91

Вот уже более тридцати лет сеть наблюдательных станций, расположенных на всех континентах, проводит измерения наклонной дальности до искусственного спутника Земли Лагеос с помощью лазерной техники. С 1992 года среди объектов наблюдений появился ещё один спутник - Лагеос-2.

С помощью разработанного пакета программ была выполнена обработка высокоточных лазерных наблюдений этих объектов.

Ряды лазерных наблюдений объектов Лагеос и Лагеос-2 с 1994 по 2002 год были разделены на пятисуточные интервалы. Для обработки наблюдений с 2003 года по 2010 год был выбран интервал, составляющий трое суток.

На каждом интервале времени были определены шесть численных зна-

0.10 0.08 0.06 0.04 0.02 0.00

метры

Лагеос

годы

1995

1997

1999

2001

2003

2005

2007

2009

0.10 0.08 0.06 0.04 0.02 0.00

метры

Лагеос-2

годы

1995

1997 1999 2001 2003 2005 2007 2009

Рис. 1: Лагеос и Лагеос-2, стандартное отклонение

чений "сглаженных'"' элементов орбиты и три значения параметров вращения Земли: координаты полюса хр\ ур и вариация продолжительности суток ЬОБ. На рис.1 представлены оценки средних квадратических погрешностей одного измерения.

На рис.2 и 3 представлены разности двух рядов координат полюса хр, ур. Из значений, полученных на основе наблюдений, были вычтены значения, опубликованные в материалах Международной службы вращения Земли.

А т / Илрлля _ 11700\

секуклы

1995 1997 1999 2001 2003 2005 2007 2009 годы

Рис. 2: Лагеос и ШИЭ, разности координат полюса

0.004 0.000 -0.004

0.004 0.000 -0.004

секунды дуги

Лхр (Лагеос-2 - 1ЕШЗ)

' •—1-:—I-1-1_I_I I

-1-1-1-1_I ■• '

1995 1997 1999 2001 2003 2005 2007 2009 годы

секунды дуги Аур (Лагеос 2 - 1ЕЛЗ)

-1-1_1_I I

1995 1997 1999 2001 2003 2005 2007 2009 годы

Рис. 3: Лагеос-2 и ГЕРБ, разности координат полюса График на рис.4 даёт представление о вариациях продожительности су-

ток, которым соответствуют вариации угловой скорости вращения Земли, на интервале времени 16 лет. Результаты опубликованы в материалах Международной службы вращения Земли.

секунды

1995

1997

1999

2001 годы

Рис. 4: Данные Международной службы вращения Земли

0.003 0.002 0.001 0.000 -0.001

секунды

•V ■, ■

вариация продолжительности суток

1995

1997

1999

2001 годы

2003

2005

2007

2009

Рис. 5: Результаты обработки наблюдений Лагеос и Лагеос-2

На рис.5 представлены результаты определения вариаций продолжительности суток. Результаты были получены на основе обработки измерений то-поцентрических дальностей до спутников Лагсос и Лагеос-2.

Поправки координат наблюдательных станций позволяют получить оценки вековых изменений положений пунктов на поверхности Земли (табл.7). В

Таблица 7: Поле скоростей для пяти пунктов

Х(м/год) К(м/год) ¿■(м/год) <г*(м/год) оу(м/год) ^¿(м/год)

Европа, Рига. 1884 А = 24°3'32"74, <р = +56°56'54:79, Я = 24.1 м

-0.0212 -0.0290 -0.0184 0.0118 0.0309 0.0151 -0.0001 0.0026 0.0046 0.0015 0.0019 0.0055 0.0021 0.0023 0.0038 0.0017 0.0026 0.0081 Лагеос Лагеос-2 ШН?

Тихий океан, Мауи, 7210 А = 203°44'38"77, р = 20°42'2бЛ1, Я = 3061.2 м

-0.0133 -0.0011 -0.0154 0.0678 0.0419 0.0589 0.0254 0.0353 0.0338 0.0019 0.0020 0.0007 0.0039 0.0034 0.0005 0.0018 0.0018 0.0004 Лагеос Лагеос-2 1ТЮ'

Африка, Хартебестёк, 7501 А = 27°41'10'30, у = -25°53'22'.<Э5, Я = 1402.8 м

-0.0038 -0.0242 -0.0090 0.0356 0.0708 0.0187 0.0223 0.0192 0.0160 0.0012 0.0026 0.0012 0.0017 0.0036 0.0008 0.0016 0.0015 0.0007 Лагеос Лагеос-2 ГГИ?

Испания, Андалусия, Сан Фернандо, 7824 А = 353°47'40$5, хр = 36°27'54»2, Я = 91.0 м

0.0080 0.0085 -0.0025 0.0267 0.0597 0.0179 0.0263 0.0247 0.0163 0.0023 0.0029 0.0025 0.0014 0.0033 0.0006 0.0015 0.0018 0.0019 Лагеос Лагеос-2 1ТЯР

Австралия, Маунт Строило, 7825 А = 149°0'35'.'63, у = -35°18'58'.11, Я = 797.8 м

-0.0435 -0.0352 -0.0363 -0.0055 -0.0039 0.0003 0.0524 0.0356 0.0457 0.0023 0.0029 0.0012 0.0018 0.0026 0.0015 0.0026 0.0032 0.0010 Лагеос Лагеос-2 !ТПГ

строке 'ТГОР' (Международная земная система отсчета) приводятся данные Международной службы'вращения Земли.

В работе получены оценки вековых изменений положений 30 пунктов измерений дальности.

На основе предлагаемых алгоритмов разработана методика расчёта прохождений искусственных спутников Земли на некотором расстоянии друг от друга. Такой подход необходим для решения следующих задач:

1. предсказание опасных сближений объектов в заданном диапазоне высот над поверхностью Земли;

2. предсказание опасных сближений конкретного космического аппарата с другими спутниками;

3. вычисление взаимных расстояний между конкретным космическим аппаратом и спутниками, отобранными в список по какому-либо признаку.

В Терскольском филиале Института астрономии Российской академии наук проводятся регулярные наблюдения космического мусора.

16 октября 2009 года во время наблюдений был обнаружен неизвестный объект 18-й звёздной величины.

Объект получил временный номер 95334. Продолжительность наблюдений составила шесть вечеров. 21 октября объект вышел из зоны видимости. Было получено 695 отдельных положений объекта - топоцентрических значений прямых восхождений и склонений в системе стандартного экватора.

На основе наблюдений было выполнено улучшение элементов орбиты и получена оценка эмпирического коэффициента отражения Сг. Так как это значение определяется формулой

Сг = 4.5606 • 10"9 кг — .

тп

то, полагая кг « 1.0, можно оценить отношение средней площади поверхности А к массе объекта т.

Оскулирующие кеплеровские элементы орбиты объекта в системе истинного экватора даты приведены в табл.8. В последней строке дана оценка отношения средней площади объекта к его массе.

Средняя квадратическая погрешность одного измерения по прямому восхождению составила 1.40", а по склонению - 0.44".

Таблица 8: Параметры орбиты объекта 95334

дата 16.10.2009

ит 00Л 00т 00.000'

а 41715.137 ± 0.020 км

е 0.029916301 ± 0.00001

г 14.005113° ± 0.0002°

П 347.868781° ± 0.0002°

ш 163.477533° ± 0.003°

М 307.308140° ± 0.003°

А/т 2.2206 ± 0.005 ы2/кг

Космический аппарат Молния 3-39 был выведен на орбиту 20 сентября 1990 года.

В феврале 2009 года эксцентриситет орбиты космического аппарата Молния 3-39 достиг критического значения. Под действием силы сопротивления атмосферы большая полуось и высота апогея орбиты стали уменьшаться. Объект завершил свой полёт во второй половине даты 8 июля 2009 года после входа в плотные слои атмосферы.

В мае 2009 года на базе Института астрономии Российской академии наук в Звенигороде старшим научным сотрудником Бахтигираевым Н.С. были выполнены позиционные наблюдения космического аппарата Молния 3-39. Получено 313 топоцентрических положений объекта: прямое восхождение и склонение в небесной системе отсчёта.

Обработка результатов позиционных наблюдений позволила уточнить начальный вектор состояния объекта на момент 18 мая 2009 года (полночь). Средняя квадратическая погрешность одного измерения составила 1.4 секунд дуги. Был выполнен численный прогноз движения космического аппарата с целью оценки момента завершения полёта.

С учётом погрешности оценки отношения площади к массе был получен интервал вероятных дат завершения полёта: от первой половины 7 июля до первой половины 9 июля.

Космические аппараты глобальной навигационной спутниковой системы (ГЛОНАСС) снабжены уголковыми отражателями для измерения топоден-трических дальностей. Длинные ряды лазерных наблюдений различных объектов ГЛОНАСС содержатся в базе данных.

В табл.9 представлены результаты обработки наблюдений космического аппарата Глонасс 109 на интервалах, равных одному месяцу:

Таблица 9: Глояасс 109. фильтрация наблюдений

дата К Aftmn (м) A/Wc (м) а (м) А/т (ы2/кг)

2011/09/01 747 -20.368 13.022 5.124 0.0321 ± 8.4 10" -6

2011/10/01 684 -38.719 20.351 10.788 0.0327 ± 2.0 10" -5

2011/11/01 724 -36.902 30.871 11.435 0.0322 ± 2.4 10" -5

2011/12/01 539 -38.109 28.598 9.780 0.0325 ± 3.2 10- -5

2012/01/01 491 -10.911 15.106 5.190 0.0325 ± 1.8 ю- -5

2012/02/01 509 -14.730 11.456 4.673 0.0323 ± 9.9 ю- -6

2012/03/01 795 -22.915 21.870 6.942 0.0328 ± 1.2 10" -5

2012/04/01 514 -27.426 33.752 12.450 0.0330 ± 2.3 10- -5

На рис.6 представлены разности между положениями объекта Глонасс 118. вычисленными с помощью улучшенных параметров движения, полученных на основе фильтрации лазерных измерений, и данными за сентябрь 2010 года из базы NGS [8], которые послужили эталонной орбитой:

невязки по оси X (метры) 4 невязки по оси 2 (метры)

60.0 3.0 -54.0

-Lv

К 1П 1F

невязки по оси Y (метры)

26 30

90.0 12.5 -66.0

.....j&

Н 1

К 1П 1К Ot\ 9R невязки по положению (метры)

30

10 15 20 25

Рис. 6: Глонасс 118, сравнение с эталонной орбитой

В табл.10 даны оценки максимальной ошибки прогноза для искусственных спутников Земли, находящихся на орбитах с различной высотой полёта. Высота полёта Лтгп дана в километрах.

Для каждого объекта оценки были получены с помощью лазерных наблюдений. На основе ЛГР значений топоцентрических дальностей на интервале времени АТа в сутках выполняется улучшение параметров движения: шесть средних элементов орбиты, эмпирический коэффициент отражения и эмпирический коэффициент ускорения. Оценка средней квадратической погрешности одного измерения в метрах дана в столбце а.

На интервале времени ДТр суток, следующим сразу за интервалом АТа, на основе улучшенных параметров движения проводится прогноз движения объекта и сравнение с измеренными значениями топоцентрических расстояний. Абсолютная величина самой большой разности измеренных и вычисленных величин в метрах обозначена \Артах\.

Таблица 10: Оценка точности прогноза

спутник fomin ATa V a Д T„ |ДРгши|

Эталон-1 19500 30 1530 7.0 30 41.2

Глопасс 109 19500 30 720 11.4 30 268.1

Глонасс 115 19500 30 960 6.2 30 194.2

Lageos G000 30 7000 1.7 30 9.2

Ajisai (EGP) 1400 30 11000 5.8 30 149.0

Метеор-ЗМ 1000 30 380 47.8 30 5349.0

Метеор-ЗМ 1000 5 60 1.2 5 71.3

Starlette 900 30 6500 18.2 30 320.0

WestPac 830 30 1070 18.4 30 2446.3

Блиц 820 30 1680 13.9 30 146.4

Stella 800 30 3500 8.8 30 281.1

Ларец 680 30 3250 15.0 30 239.6

Tandem X 510 2 210 16.4 2 886.7

Grâce A 470 2 220 4.9 2 214.1

Champ 420 2 380 6.0 2 2682.7

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

• точность обработки высокоточных лазерных наблюдений спутников Jla-геос и Лагеос-2 с помощью вычислительных программ, основанных на предлагаемых в данной работе алгоритмах численно-аналитического метода расчёта положений космических объектов, находится примерно на одном уровне с оценками точности, публикуемыми основными центрами анализа данных;

• пакет вычислительных программ LENTA на интервалах времени порядка нескольких суток уступает по точности представления наблюдений спутников типа Эталон и низкоорбитальных объектов Эйджисаи, Старлет и Стелла пакетам программ DGFI, JCET и HIT-U, использующим метод численного интегрирования уравнений движения;

• основное преимущество численно-аналитических алгоритмов в сравнении с методом численного интегрирования заключается в сокращении времени обработки очередной порции измерительной информации.

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

Выражаю признательность профессору Нестерову Вилену Валентиновичу за простую и чёткую постановку задачи в далёком 1972 году и за предоставление возможности её решения.

Работа выполнена при финансовой и моральной поддержке гранта РФФИ № 00-02-17558, научно-образовательного проекта Института астрономии РАН и госконтракта №02.740.11.0249 "Оптический мониторинг ближнего и дальнего космического пространства роботизированной сетью телескопов МАСТЕР" ГАИШ МГУ.

Список литературы

[1] Аксёнов Е.П. Теория движения искусственных спутников Земли. М. Наука. 1977.

[2J Вашковьяк М.А. Численно-аналитический метод исследования эволюции орбит далёких спутников планет. //Письма в астрономический журнал. 2005. Т.31. №1. С.66-75.

[3] Лидов М.Л. Полу аналитические методы расчёта движения спутников. //Труды Института теоретической астрономии АН СССР. 1978. Т. 17. С.54-61.

[4] Нестеров В.В. Стандарт основных вычислений астрономии. Основные алгоритмы спутниковой геодинамики. Лекции для студентов старших курсов. РФФИ. Москва. Изд-во Янус-К. 2001.

[5] Холшевников К.В. Асимптотические методы небесной механики. Ленинград. Изд-во Ленинградского государственного университета. 1985.

[6] IERS Conventions (2000). Dennis D. McCarthy (ed.) /IERS Technical note 29. Paris. 2000.

URL: http://tai.bipm.org/iers/conv2003/conv2003.html

[7] Pearlman M.R., Degnan J.J., Bosworth J.M.

The International Laser Ranging Service. //Advances in Space Research. V.30. No.2. P. 135-143. URL: http://ilrs.gsfc.nasa.gov/

[8] Precise GPS orbits computed at the National Geodetic Service. URL: http://www.ngs.noaa.gov/orbits/

[9] Tapley B.D., Schutz B.E., Eans R.J. Satellite laser ranging and its applications. //Celestial Mechanics. 1985. V.37. No.3. P.247-261.

[10] Two line elements.

URL: littp://celestrak.com/NORAD/elements/

Подписано в печать 15.01.2013 Формат 60x88 1/16. Объем 1.0 п.л. Тираж 100 экз. Заказ № 1286 Отпечатано в ООО «Соцветие красок» 119991 г.Москва, Ленинские горы, д.1 Главное здание МГУ, к. А-102

 
Текст научной работы диссертации и автореферата по астрономии, доктора физико-математических наук, Чазов, Вадим Викторович, Москва

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

имени М. В. ЛОМОНОСОВА

ГОСУДАРСТВЕННЫЙ АСТРОНОМИЧЕСКИЙ ИНСТИТУТ

имени П. К. ШТЕРНБЕРГА

На правах рукописи 05201350599 УДК 521.13

ЧАЗОВ Вадим Викторович

РАЗРАБОТКА И ПРИМЕНЕНИЕ АЛГОРИТМОВ ЧИСЛЕННО-АНАЛИТИЧЕСКОГО МЕТОДА

ВЫЧИСЛЕНИЯ ПОЛОЖЕНИЙ ИСКУССТВЕННЫХ СПУТНИКОВ ЗЕМЛИ

Диссертация на соискание учёной степени доктора физико-математических наук

Специальность 01.03.01. — Астрометрия и небесная механика

Москва - 2012

С

Содержание

Предисловие 7

1 Постановка задачи 17

1.1 Стандартные соглашения......................17

1.1.1 Системы отсчёта.......................17

1.1.2 Пространство-время.....................19

1.1.3 Окрестности Земли......................20

1.1.4 Солнечная система......................21

1.1.5 Прецессия и нутация.....................22

1.1.6 Стандартное преобразование................23

1.1.7 Координаты обсерваторий..................25

1.2 Исходные данные...........................28

1.2.1 Эволюция элементов орбиты................32

1.2.2 Уравнения движения.....................33

1.2.3 Геопотенциал.........................35

1.2.4 Численная модель......................36

1.2.5 Приливный потенциал....................39

1.3 "Стационарный" объект.......................42

1.4 Промежуточный потенциал.....................43

1.5 Канонические преобразования...................47

1.6 "Осреднённые" уравнения......................51

1.7 Постановка задачи..........................51

2 Основные алгоритмы 52

2.1 Промежуточная орбита.......................52

2.1.1 От вектора состояния к параметрам орбиты .......52

2.1.2 Постоянные интегрирования а\, а^ , аз .........53

2.1.3 Позиционные параметры а, е, 5, 5* ...........54

2.1.4 Алгоритм операций с полиномами.............57

2.1.5 Угловые переменные ¿о , до > ^о...............62

2.1.6 От параметров орбиты к вектору состояния........63

2.1.7 Проверка алгоритмов ....................65

2.1.8 Новые результаты......................67

2.2 Частные производные........................68

2.2.1 Переменные действия....................68

2.2.2 Частные производные по элементам Ь, С, Н......70

2.2.3 Частные производные по параметрам а, е, 6 ......72

2.2.4 Производные от угловых переменных...........75

2.2.5 Новые результаты......................77

2.3 Алгоритм преобразования......................78

2.3.1 Аномальный геопотенциал .................78

2.3.2 Притяжение светил......................80

2.3.3 Эффект приливов ......................81

2.3.4 Световое давление......................81

2.3.5 Положения светил......................83

2.3.6 Начальные функции.....................83

2.3.7 Элементарное слагаемое...................87

2.3.8 Сумма слагаемых.......................89

2.3.9 Формулы рекурсии......................93

2.3.10 Новые результаты......................94

2.4 Алгоритм интегрирования......................96

2.4.1 Дифференциальные соотношения .............97

2.4.2 Рекуррентный алгоритм...................99

2.4.3 Особые слагаемые......................100

2.4.4 Новые результаты......................101

2.5 Осреднённые уравнения.......................102

2.5.1 "Несингулярные" переменные................102

2.5.2 Частные производные....................103

2.5.3 Короткопериодические неравенства ............105

2.5.4 Порядок действий ......................106

2.5.5 Функции координат.....................107

2.5.6 Список упрощений......................108

2.5.7 Новые результаты......................108

2.6 Дополнительные алгоритмы.....................109

2.6.1 Координатные условия....................109

2.6.2 Влияние приливов......................112

2.6.3 Тень Земли..........................114

2.6.4 Согласование моделей....................115

2.6.5 Косвенное ускорение.....................116

2.6.6 Новые результаты......................118

3 Применение алгоритмов 119

3.1 "Средние" элементы.........................119

3.2 Вычислительные аспекты......................120

3.3 Реализация алгоритмов.......................124

3.4 Сравнительные испытания .....................125

3.5 Спутник Лагеос............................134

3.5.1 Параметры вращения Земли ................135

3.5.2 Эмпирические коэффициенты................137

3.5.3 Поле скоростей........................138

3.5.4 Сравнение результатов, пункт 7090 ............. 139

3.5.5 Сравнение результатов, пункт 7825 ............. 140

3.5.6 Сравнение результатов, пункт 7840 ............. 141

3.6 Ситуации сближения.........................142

3.7 "Стационарный" объект.......................145

3.8 Объект Молния............................148

3.9 Объекты Глонасс...........................152

3.10 Спутник Метеор-ЗМ.........................155

3.11 Спутник Блиц ............................156

3.12 Спутник Ларец............................157

3.13 Точность прогноза..........................158

Послесловие 160

Приложение 161

Промежуточный экватор.........................161

Двустрочные элементы..........................163

"Нормальные места"............................164

Правые части................................167

Атмосферная задержка..........................169

Простое дифференцирование.......................171

Полиномы и присоединённые функции Лежандра...........172

Приближённые формулы.........................173

Переменные Дудсона ...........................174

Простые разложения............................175

Две суммы элементарных слагаемых ..................176

Функции координат............................179

Фильтрация наблюдений .........................180

Лагеос и Лагеос-2. Состав наблюдений .................184

Лагеос и Лагеос-2. Поле скоростей....................186

Литература 189

Список иллюстраций

1 Стелла, средние элементы орбиты.................32

2 Стационарный объект........................44

3 Сравнение методов..........................67

4 Сравнение элементов.........................95

5 Вариация...............................112

6 Лагеос и Лагеос-2, стандартное отклонение............134

7 Лагеос и IERS, разности координат полюса............135

8 Лагеос-2 и IERS, разности координат полюса...........135

9 Данные Международной службы вращения Земли........136

10 Результаты обработки наблюдений Лагеос и Лагеос-2......136

11 Коэффициент эмпирического ускорения..............137

12 Эмпирический коэффициент отражения..............137

13 Лагеос, пункт 7090, "невязки"....................139

14 Лагеос, пункт 7090, "невязки", линейная аппроксимация.....139

15 Лагеос, пункт 7825, "невязки"....................140

16 Лагеос, пункт 7825, "невязки", линейная аппроксимация.....140

17 Лагеос, пункт 7840, "невязки"....................141

18 Лагеос, пункт 7840, "невязки", линейная аппроксимация.....141

19 Расчёт ситуаций сближения.....................143

20 "Электро-Л", прогноз ситуаций сближения на 2007 год......144

21 Объект 95334, первый вечер.....................146

22 Объект 95334, "невязки".......................146

23 Объект 95334, прогноз........................147

24 Молния 3-39, активный режим полёта...............148

25 Молния 3-39, коррекция орбиты..................148

26 Молния 3-39, пассивный режим полёта ..............149

27 Молния 3-39, фильтрация набюдений...............150

28 Молния 3-39, прогноз ........................151

29 Молния 3-39, данные "NORAD"...................151

30 Глонасс 118, фильтрация наблюдений...............153

31 Метеор-ЗМ, фильтрация наблюдений ...............155

32 Блиц, фильтрация наблюдений...................156

33 Ларец, фильтрация наблюдений..................157

34 Небесный промежуточный экватор.................162

Список таблиц

1 Объекты и наблюдения.......................31

2 Числовые значения коэффициентов зональных гармоник .... 65

3 Коэффициенты многочленов от двух переменных........85

4 Функция го/£ как сумма элементарных слагаемых .......89

5 Функция г] как сумма элементарных слагаемых.........90

6 Начальные функции: число слагаемых...............92

7 Функции координат: число слагаемых...............107

8 Оценка погрешности прогноза на основе "средних" элементов . .119

9 Производящая функция для полиномов Лежандра........121

10 Величины коэффициентов полиномов...............122

11 Объекты и функции.........................124

12 Центры анализа ...........................125

13 Объекты анализа...........................126

14 Лагеос, сравнительные испытания, пункты 7080 и 7105.....127

15 Лагеос, сравнительные испытания, пункт 7090 .......... 128

16 Лагеос, сравнительные испытания, пункт 7941 .......... 129

17 Лагеос, сравнительные испытания, несколько пунктов......130

18 Лагеос-2, сравнительные испытания................131

19 Эталон-1, сравнительные испытания................132

20 Обсерватория Яррагади, сравнительные испытания.......133

21 Поле скоростей для пяти пунктов .................138

22 Параметры орбиты объекта 95334 ................. 145

23 Глонасс 109, фильтрация наблюдений...............152

24 Разность нуль-пунктов шкал времени...............154

25 Оценка точности прогноза......................158

26 Функция л/1 — rj2 cos w как сумма элементарных слагаемых . . 177

27 Функция л/l — г]2 sin vu как сумма элементарных слагаемых . . 178

28 Функции координат: несколько слагаемых.............179

29 Состав наблюдений..........................184

30 Поле скоростей............................186

Предисловие

В диссертации представлена новая численно-аналитическая теория движения искусственных спутников Земли (ИСЗ). Теория была применена для определения параметров движения космических объектов и получения достоверных оценок геодинамических параметров на основе наблюдений.

Численно-аналитический подход позволяет использовать преимущества как аналитических методов, так и метода численного интегрирования уравнений движения.

Классическая форма аналитического способа вычисления положений небесных объектов заключается в следующей процедуре: на вековые изменения параметров орбиты накладываются долгопериодические неравенства и короткопериодические возмущения. Преимуществом аналитического метода является возможность увеличения скорости расчёта положений объекта при условии существенных ограничений на точность прогноза.

Суть метода численного интегрирования: на коротких интервалах времени параметры движения объекта аппроксимируют полиномами. Коэффициенты полиномов определяют путём вычисления значений правых частей дифференциальных уравнений и разностей этих значений в специальных точках внутри короткого интервала - одного шага интегрирования. Преимуществом численных методов является высокая точность вычислений. Малый шаг интегрирования, обеспечивая хорошую точность, требует значительных расходов вычислительного времени.

Численно-аналитический подход объединяет оба метода: часть неравенств, имеющих "короткий" период, определяется аналитически, а долгопериодические, резонансные и вековые слагаемые возмущающей функции составляют эволюционный гамильтониан - основу численного интегрирования "осреднён-ных" уравнений движения. Шаг интегрирования "осреднённых" уравнений значительно превышает значение шага интегрирования, задаваемого в численных методах. Конкретное число, разделяющее слагаемые "короткого" и "долгого" периодов, является свободным параметром.

Запуск Первого искусственного спутника Земли вызвал интерес учёных к задаче определения орбит объектов, возмущаемых аномалиями геопотенциала, притяжением Луны и Солнца, потенциалом, обусловленным приливами упругой Земли, сопротивлением атмосферы и световым давлением.

Профессором Е.П.Аксёновым была построена аналитическая теория движения ИСЗ на основе решения обобщённой задачи двух неподвижных центров. Результаты представлены в монографии [5].

Было открыто и получило признание обобщение одного из методов теории возмущений - метода канонических преобразований [156].

Тогда же был успешно применён численно-аналитический метод расчёта спутниковых орбит [36]. Важные результаты в этом направлении получены профессором М.Л.Лидовым и его учениками, некоторые итоги подведены в обзорном докладе [96], там же намечены перспективы исследований.

Тем не менее, в наши дни все центры обработки высокоточных наблюдений ИСЗ проводят расчёты с помощью программ, в которых прогноз положений объектов выполняется алгоритмом численного интегрирования дифференциальных уравнений движения [109, 204].

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

В качестве исходных материалов предлагаемого исследования были использованы:

рекомендации Международного астрономического союза [199], стандартные соглашения Международной службы вращения Земли [189], лазерные наблюдения искусственных спутников Земли [190], позиционные наблюдения космических объектов [178], база данных об элементах орбит более 14000 объектов [206], база данных о точных положениях навигационных спутников [197].

Целью исследований является решение следующих задач: вывод формул для вычисления параметров промежуточной орбиты, основанной на решении обобщённой задачи двух неподвижных центров, с точностью, ограниченной только возможностями компьютера;

разработка метода преобразования возмущающей функции на основе параметров промежуточной орбиты;

дифференцирование и интегрирование слагаемых возмущающей функции; составление и численное интегрирование "осреднённых" уравнений движения и учёт короткопериодических неравенств с использованием "несингулярных" элементов орбиты;

обработка высокоточных лазерных измерений топоцентрических дальностей до спутников Лагеос и Лагеос-2 на длительных интервалах времени и оценка значений геодинамических параметров;

разработка методики предсказания ситуаций сближения ИСЗ; фильтрация позиционных наблюдений и оценка значения отношения средней площади поверхности к массе спутника.

Стандартные соглашения Международной службы вращения Земли содержат рекомендации по обработке наблюдений искусственных спутников Земли с помощью метода численного интегрирования. Разработка методики аналитического решения задачи с учётом всех рекомендаций астрономических организаций также является целью исследования.

В тексте диссертации термин алгоритм объединяет несколько понятий: это и связанные между собой формулы и соотношения, и последовательность действий, и реализация процедуры вычислений на компьютере. Совокупность алгоритмов, предназначенная для решения поставленной задачи, в тексте называется " программным обеспечением" или " пакетом программ".

Для отладки программного обеспечения и оценки точностных характеристик алгоритмов была разработана методика сравнительных испытаний. Методика основана на использовании высокоточных измерений топоцентрических дальностей до ИСЗ. Испытания выполнены для объектов с различной высотой полёта над поверхностью Земли.

Исполнение задуманного стало возможным благодаря опоре на самые основные методы классиков небесной механики и астрометрии - последовательные приближения, промежуточные орбиты, канонические преобразования:

метод последовательных приближений применялся на всех этапах, от построения моделей движения объектов в численно-аналитической форме до дифференциального улучшения параметров модели на основе наблюдений;

в качестве промежуточной была принята орбита, построенная с помощью точного решения обобщённой задачи двух неподвижных центров;

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

и, в то же время, потребовало привлечения новых идей:

методическая точность вычисления параметров промежуточной орбиты, построенной на основе решения обобщённой задачи двух неподвижных центров, и частных производных высших порядков от любых величин по позиционным параметрам и каноническим элементам промежуточной орбиты должна быть ограничена только разрядной сеткой компьютера;

алгоритм разложения возмущающей функции в тригонометрический ряд с численными коэффициентами является неудовлетворительным по причине существенных вариаций числовых значений эксцентриситета орбиты и угла наклонения орбиты космических тел и должен быть заменён на алгоритм преобразования возмущающего гамильтониана в сумму "элементарных" слагаемых, каждое из которых является функцией как позиционных, так и угловых параметров орбиты объекта, а также величин, характеризующих изменения возмущающих факторов модели движения;

выражения для учёта возмущений от приливов упругой Земли и океанических приливов, записанные во вращающейся системе отсчёта ("геофизический" подход) и рекомендованные Международной службой вращения Земли, должны быть преобразованы к системе отсчёта, связанной с истинным экватором даты ("астрономический" подход).

Достоверность полученных результатов подтверждается примерами обработки наблюдений различных объектов, сравнением с данными, предоставляемыми Международной службой вращения Земли, и сопоставлением с параметрами движения ИСЗ, публикуемыми в Интернете. Научная новизна работы заключается:

в совокупности формул и соотношений для вычислений в рамках обобщённой задачи двух неподвижных центров с максимально возможной методической точностью;

в методе преобразования возмущающей функции задачи, заключающейся в построении всё более сложных конструкций на основе простых начальных соотношений;

в способе вычислений, позволяющим учитывать возмущающее действие факторов различного происхождения одним набором формул;

в совокупности алгоритмов построения моделей движения космических объектов с помощью численно-аналитического метода;

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

ал�