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

На правах рукописи

Заусаев Артем Анатольевич

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДВИЖЕНИЯ НЕБЕСНЫХ ТЕЛ НА ОСНОВЕ ВЫСОКОТОЧНЫХ РАЗНОСТНЫХ СХЕМ

05.13.18 - Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

САМАРА-2005

Работа выполнена на кафедре «Прикладная математика и информатика» Государственного образовательного учреждения высшего профессионального образования «Самарский государственный технический университет» (СамГТУ).

Научный руководитель:

доктор физико-математических наук, профессор Радченко Владимир Павлович.

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

доктор физико-математических наук, профессор Филатов Олег Павлович.

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

Государственное образовательное учреждение высшего профессионального образования «Казанский государственный университет» (КГУ).

Защита состоится « 2 » декабря 2005 г. в 10 — часов на заседании диссертационного совета Д 212.215.05 при Государственном образовательном учреждении высшего профессионального образования «Самарский государственный аэрокосмический университет имени академика С.П. Королева» (СГАУ) по адресу: 443086, Самара, Московское шоссе, 34.

С диссертацией можно ознакомиться в библиотеке СГАУ. Автореферат разослан « » 2005 г.

Ученый секретарь . - /'- . , , Калентьев А. А.

'///„'/ у/у,-',// -4 ;

диссертационного совета - / - / - '-л ''-у

д.т.н., профессор

^ 2&&А9

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

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

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

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

Концентрировано, результаты моделирования и накопленный информационный банк данных можно представить в виде электронного каталога ко-метных орбит. В существующих на сегодняшний день каталогах (H.A. Беляев, А. Карузи, JI. Кресак и другие авторы) данные получены на основе решения стандартной задачи п тел с учетом гравитационных сил и негравитационных эффектов, с использованием методов численного интегрирования не выше 19-го порядка, при этом динамические параметры приводятся на дискретные моменты времени с интервалами в несколько десятилетий. Поэтому учет влияния релятивистских эффектов в совокупности с повышением порядка аппроксимирующих формул численного метода позволяет не только увеличить точность и эффективность вычислений, но и разработать усовершенствованную информационную технологию при создании электронного каталога кометных орбит.

Вышеперечисленное и

РОС н БИ

| ■ !■■. т w

темы диссертации. РОС НАЦИОНАЛЬНАЯ БИБЛИОТЕКА j C.I •8

Цель работы.

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

2. Создание программного обеспечения для реализации модифицированного метода Эверхарта.

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

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

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

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

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

3. Создан комплекс нового программного обеспечения для реализации модифицированного метода Эверхарта и его применения при математическом моделировании движения небесных объектов.

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

5. На основе усовершенствованной информационной технологии создан и издан новый каталог орбитальной эволюции 164 короткопериодических комет с 1900 по 2100 гт. в электронном и печатном вариантах.

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

1. Модификация метода Эверхарта численного интегрирования дифференциальных уравнений движения небесных тел.

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

3. Модифицированная математическая модель движения небесных объектов с учетом гравитационных и релятивистских эффектов.

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

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

Связь диссертационной работы с планами научных исследований. Работа выполнялась в рамках плана НИР Самарского государственного технического университета (тема «Математическое моделирование движения небесных объектов, разработка высокоточных численных методов интегрирования уравнений движения небесных тел и их программного обеспечения») и фанта Ученого совета Самарского государственного технического университета 2005 года.

Апробация работы. Основные результаты диссертации докладывались на второй, третьей, четвертой и пятой Международных конференциях молодых ученых «Актуальные проблемы современной науки» (Самара, 2001, 2002, 2003, 2004 гг.), на двенадцатой, тринадцатой Межвузовских конференциях «Математическое моделирование и краевые задачи» (Самара, 2002, 2003), на Всероссийской научной конференции «Математическое моделирование и краевые задачи» (Самара, 2004), на Всероссийской астрономической конференции «Горизонты Вселенной» (Москва, МГУ: ГАИШ, 2004), на Международной научной конференции «Актуальные проблемы математики и механики» (Казань, 2004), на шестой Международной Петрозаводской конференции

«Вероятностные методы в дискретной математике» (Петрозаводск, 2004), на пятом Всероссийском симпозиуме по прикладной и промышленной математике (Сочи, 2004), на Международном симпозиуме «Астрономия - 2005: состояние и перспективы развития» (Москва, МГУ: ГАИШ, 2005), на второй Всероссийской научной конференции «Математическое моделирование и краевые задачи» (Самара, 2005), на Всероссийской конференции «Дифференциальные уравнения и их приложения» (Самара, 2005), на Международном форуме молодых ученых «Актуальные проблемы современной науки» (Самара, 2005), на научном семинаре «Механика и прикладная математика» Самарского государственного технического университета (рук. проф. Радчен-ко В.П., 2003-2005 гг.).

Публикации. По теме диссертации опубликовано 24 работы, список которых приведен в конце автореферата.

Личный вклад автора. Автору во всех работах, опубликованных в соавторстве, в равной степени принадлежат как постановки задач, так и результаты выполненных исследований.

Структура работы. Диссертация состоит из введения, четырех глав, общих выводов, списка литературы и двух приложений, в которых приведены выдержки из печатного варианта каталога орбитальной эволюции короткопе-риодических комет и листинги разработанных программ. Общий объем диссертации 149 страниц, включая 24 рисунка и 20 таблиц. Библиографический список включает 143 наименования.

Краткое содержание работы

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

Глава 1. Аналитический обзор и постановка задачи

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

ся схема математического моделирования движения небесных тел на основе концепции A.A. Самарского.

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

Развитие численных методов решения обыкновенных дифференциальных уравнений в XX в. тесно связано с именами Н.С. Бахвалова, H.H. Боголюбова, В.М. Волосова, Д.К. Куликова, В.Ф. Мячина, A.A. Самарского, А.Н. Тихонова, O.A. Сизовой, R. Bulirsh, J.C. Butcher, C.J. Cohen, E. Everhart, E.C. Hubbard, J. Stoer и других.

В главе описываются алгоритмы численных методов Рунге-Кутга, Адам-са-Мултона, Адамса-Бэшфорта, Ричардсона, Коуэлла, Тейлора и другие, которые применяются для решения задачи о движении небесных тел. Рассматриваются вопросы сходимости и устойчивости численного метода.

При вычислении эволюции орбит комет необходимо иметь информацию о параметрах всех небесных тел, оказывающих возмущающее действие на комету. В главе анализируется проблема разработки и создания информационных фондов координат и скоростей больших планет (Меркурий-Плутон) для использования их в численном решении различных небесно-механических задач. Данному вопросу посвящены работы H.A. Беляева, Е.И. Казимирчак-Полонской, Г.А. Красинского, Е.В. Питьевой, D. Brouwer, G.M. Clemence, W.J. Eckert, X.X. Newhall, E.M. Standish, Jr. и J.G. Williams и других.

Отмечается роль электронного каталога как концентрированного результата численного эксперимента в информационных технологиях получения, накопления и систематизации новых знаний в небесной механике. Указывается, что каталоги короткопериодических комет содержат информацию о распределении комет в Солнечной системе, дают общее представление об эволюции параметров орбит с течением времени. Приводится обзор кометных каталогов, в частности, работ H.A. Беляева, Э. Галлея, А. Карузи, Л. Кресака, Б. Марсдена, И. Хасегава и других.

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

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

В главе приводится подробная схема процесса математического моделирования эволюции короткопериодических комет на основе концепции А.А.Самарского.

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

Глава 2. Разработка математической модели и метода совместного решения дифференциальных уравнений движения Солнца, больших планет, Луны я комет на основе алгоритма Эверхарта

В пункте 2.1 рассматриваются основные математические модели движения небесных объектов.

Дифференциальные уравнения движения в стандартной ньютоновской задаче к тел в векторной форме имеют вид1:

где r{x,y,z) - гелиоцентрический радиус-вектор тела, движение которого исследуется; т - масса этого тела; ^(x^y^zj - гелиоцентрический радиус-вектор /-того возмущающего тела; т1 - масса /-того возмущающего тела; к2-Гауссова постоянная,

f2 = х2 + у2 + z2, А* = (х, - х)2 +(у, - у)2 +(z, - z)2, г2 = X/2 + у,г + Z,2.

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

(О

' Чеботарев Г А Аналитические и численные методы небесной механики M ,Л Наука, 1965 368 с

8

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

рые имеют следующий вид2:

.. 2+ I y(U<

1*1 '¡J I L 'I* L k*J rjk

г

гл Ч с

I V,

У

2(1 + Г)г • _ J_

с2 ' ' 2сг

(Г, ~rj)-r,

2

IX "i fe - rJ )] b+2r)r, - a+2Y)rj DCA - т, >+

+

с

J - ■ у

1* ГУ 14=1 Гш

где г„г„ г, - радиус-векторы координат, скоростей, ускорений в барицентрической системе координат /-того тела; - k2ntj, где к2 - Гауссова постоянная и mj - массау'-того тела; rt] = - г, | ; р и у - релятивистские параметры, р = у = \\ и, =|г,| и с - скорость света.

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

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

2

Newhall X X, Standish E.M., Williams Jr and J G DE 102 a numerically integrated ephemens of the Moon and planète spanning forty-four centuries//Astron Astrophys 1983 №125 P150-167

сближения с Юпитером, приводит к лучшему согласованию теоретических расчетов с наблюдениями, чем без учета этих эффектов.

В пункте 2.2 рассматривается метод численного интегрирования Эверхар-та и производится его модификация.

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

Алгоритм метода Эверхарта4 при решения уравнения

* = (3)

состоит в представлении правой части дифференциального уравнения в виде временного ряда

= = ^ + А1( + А212 +... + А„Г. (4)

Интегрируя (4), получим выражения для определения координат и скоростей:

х = х, +ХЛ + Р, — + А, — + ... + Л--, (5)

1 1 '2 '6 "(п + 1Х» + 2)

. . /2 (3

х = х1+Р11 + А1— + А2Т + ... + А„-(6)

2 3 п +1

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

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

Алгоритм и программа численного интегрирования методом Эверхарта ранее были разработаны до 27 порядка, однако при возрастании порядка свыше 19-го не происходило улучшения точности вычислений. Автором настоя-

Совремеяные численные методы решения обыкновенных дифференциальных уравнений Подред Дж Холла, Дж У ana М Мир 1979 312 с.

4

EverhartE Implicit single methods for integrating orbits // Celestial mechanics 1974 № 10. P 35-55. 10

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

В данной диссертационной работе в алгоритм метода Эверхарта было введено дополнительное условие: начиная с 19-го порядка значение коэффициента А„ приравнивалось к нулю, что влекло выполнение необходимого условия сходимости степенных рядов. В соответствии с этой модификацией, коэффициенты А, также принимали другие значения. Данный прием позволил разработать алгоритм и программу на языке С++, эффективность и точность которых возрастали применительно к планетной задаче при увеличении порядка до 31.

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

В пункте 2.5 дается понятие времени в задачах небесной механики, описываются параметры, характеризующие траекторию движения короткопериодических комет (элементы орбиты), часть из них обозначена на рис. 1: Т момент прохождения кометой своего перигелия (Р); М - средняя аномалия,

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

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

В завершении сделаны выводы по главе 2.

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

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

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

Рис 2 Схема формирования банков данных координат и скоростей

Полученные результаты для внутренних планет (Меркурий-Марс) сопоставлены с радиолокационными наблюдениями на интервале времени с 1961 по 2000 гг., также проведены сравнения найденных значений вблизи начального и конечного моментов интегрирования с данными одной из самых современных теорий движения больших планет - DE 405 5. На основании сопоставления сделан вывод о высокой степени точности полученных координат и скоростей, в частности, для планет максимальное расхождение в координатах на 1660 год, по сравнению с DE 405, имеет место у Меркурия и составляет 8-Ю"7 а.е., для Луны - 310"6 а.е.

Далее, путем совместного численного интегрирования уравнений движения больших планет, Луны, Солнца и комет на интервале времени с 1900 по 2100 гг. создан банк данных координат, скоростей 164 короткопериодических комет на те же моменты времени, что и для планет.

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

В пункте 3.2 исследуется эффективность метода Эвсрхарта при численном интегрировании уравнений движения короткопериодических комет.

С целью проведения тестовых расчетов для определения оптимальных параметров (шага и порядка) в модифицированном методе Эверхарта при решении уравнений движения короткопериодических комет, было отобрано 10 объектов. В их число вошли кометы, обладающие характерными особенностями: обратным и прямым движением, имеющие тесные сближения с планетами-гигантами и не имеющие таковых, кометы с малыми и большими периодами обращения и т.д. Совместное численное интегрирование уравнений движения больших планет, Луны, Солнца и кометы в барицентрической системе координат проводилось от момента последнего зарегистрированного прохождения кометы через перигелий до ближайшего к 1900 году момента времени, на который известны орбитальные параметры комет. Исследование проводилось для трех порядков метода - 19, 23 и 27. Метод считался оптимальным, если при заданном порядке и шаге получен результат, удовлетворяющий требуемой точности и обеспечивающий максимальное быстродействие работы программы.

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

5 Standish EM JPL Planetary and Lunar Ephemendes, DE405/LE405 // Jet Prop Lab Technical Report IOM 312 F-

048 1998 PI-7

ный метод Эверхарта с порядками 19, 23, 27, при этом оптимальная величина шага для порядков 19 и 23 является нестабильной, в то время как метод 27 порядка с шагом 6 дней обеспечивает требуемую точность и быстродействие, то есть является наиболее эффективным.

В пункте 3.3 проведено исследование устойчивости метода Эверхарта при вычислении орбитальной эволюции короткопериодических комет.

Сравнивая результаты тестовых вычислений на основе возмущенных и невозмущенных начальных параметров 10 короткопериодических комет, делается вывод об устойчивости используемого метода. Отмечается, что к числу комет с менее устойчивыми решениями относятся кометы с большими эксцентриситетами, а также кометы, имеющие очень тесные сближения с Юпитером (менее 0,1 а.е.).

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

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

В пунктах 3.5 - 3.7 анализируются прикладные вопросы и проводится статистический анализ распределения элементов орбит 164 короткопериодических комет. На основании проведенного анализа делается ряд выводов о распределении комет в Солнечной системе.

Отмечается, что под влиянием возмущающего действия Юпитера в 30 случаях на рассматриваемом интервале времени происходит существенное изменение орбит комет. Выполненное исследование показало, что умеренные и тесные сближения комет с Юпитером - есть закономерное явление, которое реализуется со временем для каждой кометы семейства Юпитера. Время и величина сближений короткопериодических комет с большими планетами сопоставлены с данными каталога А. Карузи, Л. Кресака, Д. Валсеччи. Несмотря на различие применяемых численных теорий движения, установлено согласование данных о времени и минимальных расстояниях сближений комет с большими планетами.

В пункте 3.8 оценивается точность проведенных расчетов эволюции короткопериодических комет. Приводится сравнение результатов исследования с данными каталога Б. Марсдена, согласованными с наблюдениями. Сопоставление элементов орбит проводилось на момент прохождения кометы через перигелий, ближайший к 1900 году. Расхождения в средней аномалии М у 115 комет из 164 не превышают 0,05 градуса, что соответствует относительной погрешности 0,014%. Наибольшие отклонения установлены у четырех объектов: 16 P/Brooks 2 (6,03%), 22 P/Kopff (2,03%), 24 P/Schaumasse (1.18%), 52 P/Harrington-Abell (3.64%). Основной причиной столь существенных расхождений могут быть столкновения комет с астероидами или другие эффекты. Расхождение в элементах орбит для остальных 160 короткопериодических комет устраняется путем стандартного учета негравитационных сил в уравнениях движения.

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

В завершении делаются выводы по главе 3.

Глава 4. Разработка программного обеспечения и его компьютерная реализация при математическом моделировании движения короткопериодических комет

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

На базе сформированных банков данных координат и скоростей больших планет, Луны, Солнца и 164 короткопериодических комет (рис. 2) разработан электронный каталог орбитальной эволюции короткопериодических комет с 1900 по 2100 гг. Каталог реализован программными средствами С++ Builder, имеет оконный интерфейс и предоставляет пользователю следующие возможности:

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

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

- построение графиков изменения орбитальных элементов.

В подготовленном и изданном печатном варианте каталога приводится информация об изменении элементов орбит 164 короткопериодических комет

на интервале с 1900 по 2100 годы на момент начала года с шагом 10 лет, указываются тесные сближения с большими планетами, а также даются графики изменения орбитальных элементов на 200-летнем интервале времени.

Основные результаты выполненных исследований состоят в следующем:

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

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

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

4. Разработана новая информационная среда, в рамках которой созданы банки данных координат и скоростей больших планет, Луны, Солнца (на интервале 1660-2200 гг.) и 164 короткопериодических комет (на интервале 19002100 гг.).

5. На основе численного эксперимента:

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

- исследована и показана устойчивость применяемого численного метода;

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

- проведено сравнение и классификация короткопериодических комет по значениям элементов орбит;

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

6. На базе составленных банков данных, разработанных алгоритмов и программ создан новый открытый информационный банк данных - каталог орбитальной эволюции короткопериодических комет с 1900 по 2100 гг. в электронном и печатном вариантах.

Основные результаты диссертации опубликованы в следующих работах:

1. Заусаев А.Ф., Заусаев A.A. Каталог орбитальной эволюции короткопериодических комет с 1900 по 2100 гг. М.: Машиностроение - 1. 2005. 346 с. (личный вклад автора - 173 страницы).

2. Заусаев А.Ф., Заусаев А А. Математическое моделирование и анализ эволюции орбит 25 короткопериодических комет и их сближений с большими планетами // Вестник СамГТУ. Серия: Физико-математические науки. Вып. 16. Самара 2002 С.57-61. (авт. - 3 е.).

3. Заусаев А.Ф., Заусаев A.A., Ольхин А.Г. Численное интегрирование уравнений движения больших планет (Меркурий-Плутон) и Луны с учетом радиолокационных наблюдений // Вестник СамГТУ. Серия: Физико-математические науки. Вып. 26. Самара. 2004 С.43-47. (авт. - 2 е.).

4. Заусаев А.Ф, Заусаев A.A., Ольхин А.Г. Оценка точности метода Эверхар-та при решении уравнений движения больших планет на интервале времени 10 000 лет // Вестник СамГТУ. Серия: Физико-математические науки. Вып. 30. Самара. 2004. С.108-113. (авт. -2 е.).

5. Заусаев A.A. Разработка каталога короткопериодических комет на интервале времени 1900-2100 годы // Вестник СамГТУ. Серия: Физико-математические науки. Вып. 34. Самара. 2005. С.40-46.

6. Заусаев А. Ф., Заусаев A.A. Исследование влияния релятивистских эффектов на эволюцию орбит короткопериодических комет // Математическое моделирование и краевые задачи. Труды двенадцатой межвузовской конференции. Часть 3. Самара. 2002. С. 55-61. (авт. - 3 е.).

7. Заусаев А.Ф., Заусаев А.А, Ольхин А.Г. Исследование эффективности метода Эверхарта численного интегрирования уравнений движения больших планет // Математическое моделирование и краевые задачи. Труды тринадцатой межвузовской конференции. Часть 3. Самара. 2003. С. 70-74. (авт. - 2 е.).

8. Заусаев A.A. Исследование вклада релятивистских эффектов в эволюцию короткопериодических комет // Математическое моделирование и краевые задачи. Труды всероссийской научной конференции. Часть 3. Самара. 2004. С. 116-119.

9. Заусаев A.A. Расчет эволюции короткопериодических комет с учетом эффектов специальной теории относительности // Актуальные проблемы современной науки. Труды 5-й Международной конференции молодых ученых. Части 1,2. Самара. 2004. С. 45-47.

10. Заусаев А.Ф., Заусаев A.A., Олъхин А.Г. Исследование эффективности метода Эверхарта на примере решения задач небесной механики // Актуальные проблемы математики и механики. Труды математического центра им. Н.И.Лобачевского. Том 25. Казань: КГУ. 2004. С. 122-123. (авт. - 1 е.).

11. Заусаев A.A. Исследование устойчивости численного метода Эверхарта при интегрировании уравнений движения короткопериодических комет // Математическое моделирование и краевые задачи. Труды второй всероссийской научной конференции. Часть 3. Самара. 2005. С. 105-107.

12. Заусаев А А. Электронный каталог орбитальной эволюции короткопериодических комет с 1900 по 2100 гг. // Актуальные проблемы современной науки. Труды 1-го Международного форума молодых ученых. Части 1,2. Самара. 2005. С.113-115.

13. Заусаев A.A. Численное интегрирование уравнений движения короткопериодических комет с учетом гравитационных и негравитационных эффектов // Актуальные проблемы современной науки. Тезисы докладов 3-й Международной конференции молодых ученых. Часть 1. Самара. 2002. С.17.

14. Заусаев A.A. Использование метода Эверхарта при интегрировании уравнений движения кометы Галлея // Актуальные проблемы современной науки. Тезисы докладов 4-й Международной конференции молодых ученых. Части 1-3. Самара. 2003. С.38.

15. Заусаев A.A. Исследование влияния релятивистских эффектов на движение короткопериодических комет // ВАК-2004: Горизонты Вселенной. Тезисы докладов Всероссийской астрономической конференции. Москва: МГУ, ГАИШ. 2004. С.226.

16. Заусаев А.Ф., Заусаев A.A., Олъхин А.Г. Применение метода Эверхарта 31 порядка для решения уравнений движения больших планет // ВАК-2004: Горизонты Вселенной. Тезисы докладов на Всероссийской астрономической конференции. Москва: МГУ, ГАИШ. 2004. С.209. (авт. - 0.3 е.).

17. Заусаев А.Ф., Заусаев A.A., Олъхин А Г. Применение метода Эверхарта 31 порядка для решения уравнений движения различных небесных объектов /7 Обозрение прикладной и промышленной математики. Том 11. Вып. 3. Москва. 2004. С.636. (авт. - 0.3 е.).

18. Заусаев A.A. Влияние релятивистских эффектов на движение короткопериодических комет // Обозрение прикладной и промышленной математики. Том 11. Вып. 3. Москва. 2004. С.637.

19. Заусаев A.A. Учет эффектов специальной теории относительности при исследовании эволюции короткопериодических комет // Обозрение прикладной и промышленной математики. Том 11. Вып. 4. Москва. 2004. С.810.

20. Заусаев А. Ф., Заусаев A.A., Ольхин А.Г. Использование метода Эверхарта для решения уравнений движения больших планет и Луны // Обозрение прикладной и промышленной математики. Том 11. Вып. 4. Москва. 2004. С.810. (авт.-0.3 е.).

21. Заусаев А.Ф., Заусаев A.A., Ольхин А.Г. Применение метода Эверхарта высокого порядка к решению задач небесной механики // Астрономия -2005: состояние и перспективы развития. Тезисы докладов международного симпозиума. Труды ГАИШ. Том LXXVIII. Москва: МГУ, ГАИШ. 2005. С. 12. (авт. -0.3 е.).

22. Заусаев A.A. Разработка каталога орбитальной эволюции короткопериодических комет // Астрономия - 2005: состояние и перспективы развития. Тезисы докладов международного симпозиума. Труды ГАИШ. Том LXXVm. Москва: МГУ, ГАИШ. 2005. С. 12.

23. Заусаев A.A., Радченко В.П. О составлении каталога орбитальной эволюции короткопериодических комет с 1900 по 2100 гг. // Дифференциальные уравнения и их приложения. Тезисы докладов Всероссийской конференции Самара: СамГУ. 2005. С. 45-46. (авт. - 1 е.).

24. Заусаев А.Ф., Заусаев A.A., Ольхин А.Г. Модифицированный метод Эверхарта для решения обыкновенных дифференциальных уравнений // Обозрение прикладной и промышленной математики. Том 12. Вып. 2. Москва. 2005. С.369. (авт.-0.3 е.).

»20867

РНБ Русский фонд

2006-4 19226

Автореферат отпечатан с разрешения диссертационного совета Д 212.215.05 (протокол № 52 от 21 октября 2005 года)

Подписано в печать 24 октября 2005 г Заказ № 420 Тираж 100 экз. Отпечатано на ризографе Самарский государственный технический университет Отдел типографии и оперативной полиграфии 443100, г. Самара, ул Молодогвардейская, 244

ВВЕДЕНИЕ

ГЛАВА 1. Аналитический обзор и постановка задачи

1.1. Основные сведения о проблеме изучения движения небесных объектов.

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

1.3. Этапы исследования движения больших планет.

1.4. Основные сведения о короткопериодических кометах.

• 1.5. Обзор кометных каталогов.

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

ГЛАВА 2. Разработка математической модели и метода совместного решения дифференциальных уравнений движения Солнца, больших планет, Луны и комет на основе алгоритма Эверхарта

2.1. Дифференциальные уравнения движения.

2.2. Выбор и обоснование метода решения дифференциальных уравнений движения. Модифицированный метод Эверхарта.

2.3. Интегрируемость и устойчивость в задачах небесной механики.

2.4. Негравитационные эффекты в эволюции короткопериодических комет.

2.5. Обоснование выбора пространственно-временных систем для описания движения небесных объектов.

2.5.1. Время в задачах небесной механики.

2.5.2. Основные параметры, описывающие траекторию движения короткопериодических комет.

2.5.3. Основные системы координат, используемые при решении задач небесной механики.

2.5.4. Вычисление элементов гелиоцентрической орбиты по положению и скорости в начальный момент.

2.5.5. Вычисление прямоугольных координат и компонент скорости по элементам орбит.

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

3.1 Создание банков данных координат и скоростей больших планет (Меркурий-Плутон), Луны, Солнца и короткопериодических комет в барицентрической системе координат.

3.2 Исследование эффективности метода Эверхарта при численном интегрировании уравнений движения короткопериодических комет.

3.3 Исследование устойчивости метода Эверхарта при вычислении орбитальной эволюции короткопериодических комет.

3.4 Учет негравитационных эффектов для короткопериодических комет

3.5 Распределение и классификация короткопериодических комет по значениям критерия Тиссерана.

3.6 Анализ эволюции орбит короткопериодических комет на интервале времени с 1900 по 2100 годы.

3.7 Короткопериодические кометы, имеющие тесные сближения с Юпитером.

3.8 Оценка точности расчетов орбитальной эволюции короткопериодических комет

ГЛАВА 4. Разработка программного обеспечения и его компьютерная реализация при математическом моделировании движения короткопериодических комет и составление каталога

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

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

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

Концентрировано, результаты моделирования и накопленный информационный банк данных можно представить в виде электронного каталога кометных орбит. В существующих на сегодняшний день каталогах (H.A. Беляев, А. Карузи, JI. Кресак и другие авторы) данные получены на основе решения стандартной задачи п тел с учетом гравитационных сил и негравитационных эффектов, с использованием методов численного интегрирования не выше 19-го порядка, при этом динамические параметры приводятся на дискретные моменты времени с интервалами в несколько десятилетий. Поэтому учет влияния релятивистских эффектов в совокупности с повышением порядка аппроксимирующих формул численного метода позволяет не только увеличить точность и эффективность вычислений, но и разработать усовершенствованную информационную технологию при создании электронного каталога кометных орбит.

Вышеперечисленное и определяет актуальность темы диссертации.

Цель работы.

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

• проксимирующей формулы.

2. Создание программного обеспечения для реализации модифицированного метода Эверхарта.

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

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

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

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

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

3. Создан комплекс нового программного обеспечения для реализации модифицированного метода Эверхарта и его применения при математическом моделировании движения небесных объектов.

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

5. На основе усовершенствованной информационной технологии создан и издан новый каталог орбитальной эволюции 164 короткопериодиче-ских комет с 1900 по 2100 гг. в электронном и печатном вариантах.

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

1. Модификация метода Эверхарта численного интегрирования дифференциальных уравнений движения небесных тел.

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

3. Модифицированная математическая модель движения небесных объектов с учетом гравитационных и релятивистских эффектов.

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

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

Связь диссертационной работы с планами научных исследований. Работа выполнялась в рамках плана НИР Самарского государственного технического университета (тема «Математическое моделирование движения небесных объектов, разработка высокоточных численных методов интегрирования уравнений движения небесных тел и их программного обеспечения») и гранта Ученого совета Самарского государственного технического университета 2005 года.

Апробация работы. Основные результаты диссертации докладывались на второй, третьей, четвертой и пятой Международных конференциях молодых ученых «Актуальные проблемы современной науки» (Самара, 2001, 2002, 2003, 2004 гг.), на двенадцатой, тринадцатой Межвузовских конференциях «Математическое моделирование и краевые задачи» (Самара, 2002,

2003), на Всероссийской научной конференции «Математическое моделирование и краевые задачи» (Самара, 2004), на Всероссийской астрономической конференции «Горизонты Вселенной» (Москва, МГУ: ГАИШ, 2004), на Международной научной конференции «Актуальные проблемы математики и механики» (Казань, 2004), на шестой Международной Петрозаводской конференции «Вероятностные методы в дискретной математике» (Петрозаводск,

2004), на пятом Всероссийском симпозиуме по прикладной и промышленной математике (Сочи, 2004), на Международном симпозиуме «Астрономия -2005: состояние и перспективы развития» (Москва, МГУ: ГАИШ, 2005), на второй Всероссийской научной конференции «Математическое моделирование и краевые задачи» (Самара, 2005), на Всероссийской конференции «Дифференциальные уравнения и их приложения» (Самара, 2005), на Международном форуме молодых ученых «Актуальные проблемы современной науки» (Самара, 2005), на научном семинаре «Механика и прикладная математика» Самарского государственного технического университета (рук. Ф проф. Радченко В.П., 2003-2005 гг.).

Публикации. По теме диссертации опубликовано 24 работы, список которых приведен в конце автореферата.

Личный вклад автора. Автору во всех работах, опубликованных в соавторстве, в равной степени принадлежат как постановки задач, так и результаты выполненных исследований.

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

Основные результаты выполненных исследований состоят в следующем:

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

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

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

4. Разработана новая информационная среда, в рамках которой созданы банки данных координат и скоростей больших планет, Луны, Солнца (на интервале 1660-2200 гг.) и 164 короткопериодических комет (на интервале 1900-2100 гг.).

5. На основе численного эксперимента:

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

- исследована и показана устойчивость применяемого численного метода;

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

- проведено сравнение и классификация короткопериодических комет по значениям элементов орбит;

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

6. На базе составленных банков данных, разработанных алгоритмов и программ создан новый открытый информационный банк данных — каталог орбитальной эволюции короткопериодических комет с 1900 по 2100 гг. в электронном и печатном вариантах.

ЗАКЛЮЧЕНИЕ

1. Абалакин В.К., Аксенов Е.П., Гребенников Е.А., Демин В.Г., Рябов Ю.А. Справочное руководство по небесной механике и астродинамике. М.: Наука, 1976. 862 с.

2. Альвен X., Аррениус Г. Эволюция Солнечной системы. М.: Мир. 1979. 512 с.

3. Арнольд В.И., Козлов В.В., Нейштадт А.И. Математические аспекты классической и небесной механики. М.: Эдиториал УРСС. 2002. 416 с.

4. Батурин А.П. Определение орбит комет, наблюдавшихся в одном и двух появлениях // Исслед. по баллистике и смеж. вопр. мех. 1999. 3 . С. 128129.

5. Бахвалов Н.С. Численные методы Т.1. М.: Наука, 1975. 632 с.

6. Белый Ю.А. Иоганн Кеплер. М.: Наука. 1971.

7. Беляев Н.А, Чурюмов К.И. Комета Галлея и ее наблюдение. М.: Наука, 1985. 272 с.

8. Беляев H.A. Эволюция орбиты кометы Даниэля 1909IV за 400 лет (16602060 гг.) // Бюлл.ИТА. 1966. Т. 10. №10. С.696-710.

9. Беляев H.A., Кресак Л., Питтих Э.М., Пушкарев А.Н. Каталог коротко-периодических комет. Братислава, 1986. 398 с.

10. Бордовицина Т.В. Алгоритмы численного моделирования движения малых тел Солнечной системы // Астрометрия, геодинамика и небесная механика на пороге XXI века. СПб. 2000. С. 208-211.

11. Бордовицина Т.В. Алгоритмы численного моделирования движения малых тел Солнечной системы // Тр. ИПА РАН. 2001. 6. С.160-169.

12. Бордовицина Т.В. Современные численные методы в задачах небесной механики. М.: Наука. 1984. 136 с.

13. Брауэр Д., Клеменс Д. Методы небесной механики. М.: Мир. 1964. 516 с.

14. Брумберг В.А. Методика определения релятивистских планетных возмущений в теориях движения больших планет // Тр. ИПА РАН. 1999. 4. С. 199-224.

15. Брумберг В.А. Релятивистская небесная механика. М.: Наука, 1972. 382 с.

16. Всехвятский С. К. Физические характеристики комет. М.: Физматгиз. 1958.

17. Всехсвятский С.К. Природа и происхождение комет и метеорного вещества. М.: Просвещение. 1967.

18. Григорян О.Ф., Медведев Ю.Д., Томанов В.П. Гипотеза дополнительного ускорения в движении кометы Харрингтона-Абеля от Юпитера // Ас-трон. вестник. 2004. Т. 38. 5. С. 452-461.

19. Дубошин Г.Н. Небесная механика. Основные задачи и методы. М.: Физматгиз. 1963. 588 с.

20. Дубяго А.Д. Движение периодической кометы Брукса 1925-1960 // Уч. зап. Казанск. унив.1956. 116. 6. С. 3-31.

21. Дубяго А.Д. О вековом ускорении движения периодических комет // Астр. журн. 1948. 25. 6. С. 361-368.

22. Дубяго А.Д. Движение периодической кометы Брукса с 1883 по 1946 г. // Уч. зап. Казанск. унив.1950. 110. 8. С. 5-44.

23. Емелъяненко В.В. Исследование долговременной эволюции кометных орбит на основе отображений // Астрометрия, геодинамика и небесная механика на пороге XXI века. СПб. 2000. С. 223-224.

24. Емелъяненко Н.Ю. Динамика орбит комет при тесном сближении с Юпитером. Анализ длительности сближений // Астрон. вестник. 2003. 37. 2. С. 174-182.

25. Емелъяненко Н.Ю. Короткопериодические кометы с высоким значением постоянной Тиссерана. III. Кинематика низкоскоростных сближений // Астрон. вестн. 2003. 37. 1. С. 72-79.

26. Емельяненко Н.Ю. Эволюция комет, имеющих тесные сближения с

27. Юпитером. I. Анализ влияния ошибок исходной системы элементов //

28. Астрон. вестник. 1992. 5. С.24-29.

29. Заусаев A.A. Влияние релятивистских эффектов на движение короткопериодических комет // Обозрение прикладной и промышленной математики. Том 11. Вып. 3. Москва. 2004. С.637.

30. Заусаев A.A. Использование метода Эверхарта при интегрировании уравнений движения кометы Галлея // Актуальные проблемы современной науки. Тезисы докладов 4-й Международной конференции молодых ученых. Части 1-3. Самара. 2003. С.38.

31. М 29. Заусаев A.A. Исследование вклада релятивистских эффектов в эволюциюкороткопериодических комет // Математическое моделирование и краевые задачи. Труды всероссийской научной конференции. Часть 3. Самара. 2004. С. 116-119.

32. Заусаев A.A. Исследование влияния релятивистских эффектов на движение короткопериодических комет // ВАК-2004: Горизонты Вселенной. Тезисы докладов на Всероссийской астрономической конференции. Москва: МГУ, ГАИШ. 2004 С.226.

33. Заусаев A.A. Исследование движения короткопериодических комет численными методами // Тезисы докладов XXVIII Самарской областной научной конференции. Часть 1. Самара. 2002. С.90.

34. Заусаев A.A. Разработка каталога короткопериодических комет на интервале времени 1900-2100 годы // Вестник СамГТУ. Серия: Физико-математические науки. Вып. 34. Самара. 2005. С.40-46.

35. Заусаев A.A. Разработка каталога орбитальной эволюции короткоперио-Щ дических комет // Астрономия 2005: состояние и перспективы развития. Тезисы докладов международного симпозиума. Труды ГАИШ. Том LXXVIII. Москва: МГУ, ГАИШ. 2005. С. 12.

36. Заусаев A.A. Расчет эволюции короткопериодических комет с учетом эффектов специальной теории относительности // Актуальные проблемы современной науки. Труды 5-й Международной конференции молодых ученых. Части 1, 2. Самара. 2004. С. 45-47.

37. Заусаев A.A. Учет эффектов специальной теории относительности при исследовании эволюции короткопериодических комет // Обозрение прикладной и промышленной математики. Том 11. Вып. 4. Москва. 2004. С.810.

38. Заусаев A.A. Численное интегрирование уравнений движения короткопериодических комет // Актуальные проблемы современной науки. Тезисы докладов 2-й Международной конференции молодых ученых. Часть 1. Самара. 2001. С.34.

39. Заусаев A.A. Электронный каталог орбитальной эволюции короткопериодических комет с 1900 по 2100 гг. // Труды 1-го Международного форума молодых ученых «Актуальные проблемы современной науки». Части 1,2. Самара. 2005. С.113-115.

40. Заусаев A.A., Радченко В.П. О составлении каталога орбитальной эволюции короткопериодических комет с 1900 по 2100 гг. // Дифференциальные уравнения и их приложения. Тезисы докладов Всероссийской конференции Самара: СамГУ. 2005. С. 45-46.

41. Заусаев А.Ф. Эволюция орбит малых тел Солнечной системы, сближающихся с Землей // Автореф. доктор, диссер. М. 1994. 20 с.

42. Заусаев А. Ф., Заусаев A.A. Исследование влияния релятивистских эффектов на эволюцию орбит короткопериодических комет // Математическое моделирование и краевые задачи. Труды двенадцатой межвузовской конференции. Часть 3. Самара. 2002. С. 55-61.

43. Заусаев А.Ф., Заусаев A.A. Каталог орбитальной эволюции короткопериодических комет с 1900 по 2100 гг. М.: Машиностроение 1. 2005. 346 с.

44. Заусаев А.Ф., Заусаев A.A. Математическое моделирование и анализ эволюции орбит 25 короткопериодических комет и их сближений с большими планетами // Вестник СамГТУ. Серия: Физико-математические науки. Вып. 16. Самара. 2002. С.57-61.

45. Заусаев А.Ф., Заусаев A.A., Олъхин А.Г. Использование метода Эверхарта для решения уравнений движения больших планет и Луны // Обозрение прикладной и промышленной математики. Том 11. Вып. 4. Москва. 2004. С.810.

46. Заусаев А. Ф., Заусаев A.A., Олъхин А.Г. Модифицированный метод Эверхарта для решения обыкновенных дифференциальных уравнений // Обозрение прикладной и промышленной математики. Том 12. Вып. 2. Москва. 2005. С.369.

47. Заусаев А.Ф., Заусаев A.A., Олъхин А.Г. Оценка точности метода Эверхарта при решении уравнений движения больших планет на интервалевремени 10 ООО лет // Вестник СамГТУ. Серия: Физико-математические науки. Вып. 30. Самара. 2004. С.108-113.

48. Заусаев А. Ф., Заусаев A.A., Ольхин А.Г. Применение метода Эверхарта 31 порядка для решения уравнений движения различных небесных объектов // Обозрение прикладной и промышленной математики. Том 11. Вып. 3. Москва. 2004. С.636.

49. Заусаев А.Ф., Заусаев A.A., Ольхин А.Г. Численное интегрирование уравнений движения больших планет (Меркурий-Плутон) и Луны с учетом радиолокационных наблюдений // Вестник СамГТУ. Серия: Физико-математические науки. Вып. 26. Самара. 2004 С.43-47.

50. Казимирчак-Полонская Е.И. Обзор исследований тесных сближений короткопериодических комет с Юпитером (1770—1960) // Тр. ИТА. 1961.7. С.19-190.

51. Казимирчак-Полонская Е.И Некоторые актуальные задачи кометной астрономии с современных позиций небесной механики // Труды ИТА АН СССР. 1967. 12. С.3-23.

52. Казимирчак-Полонская Е.И. Эволюция орбит короткопериодических комет на интервале 1660-2060 гг. и роль внешних планет в этой эволюции //Астрон.журнал. Л. 1967. № 2. С.439-460.

53. Коперник Н. О вращениях небесных сфер. Малый комментарий. Послание против Вернера. Упсальская запись. М.: Наука, 1964.

54. Красинский Г.А., Питъева Е.В., Свешников М.Л., Свешникова Е.С. Уточнение эфемерид внутренних планет и Луны по радиолокационным, лазерным и мередианным измерениям 1961-1980 гг. // Бюлл. ИТА. 1982. 3. С.145-163.

55. Куликов Д. К. Интегрирование уравнений движения небесной механики на электронных вычислительных машинах по квадратурному методу Коуэлла с автоматическим выбором шага // Бюлл. „ИТА. 1960. №10. С.770-797.

56. Левин Б.Ю. Физическая теория метеоров и метеорное вещество в солнечной системе. М.: Изд-во АН. 1956. 290 с.

57. Ляпунов A.M. Общая задача об устойчивости движения. Гостехиздат. 1950.

58. Медведев Ю.Д. Обзор исследований по динамике комет // Астрометрия, геодинамика и небесная механика на пороге XXI века. СПб. 2000. С.200-203.

59. Монтенбрук О., Пфлегер Т. Астрономия на персональном компьютере. СПб.: Питер. 2002. 320 с.

60. Мысовских НИ Лекции по методам вычислений. М.: Физматгиз.1962.

61. Мячин В.Ф., Сизова O.A. Совместное интегрирование уравнений небесной механики численным методом Тейлора-Стеффенсона // Бюлл. ИТА. 1970. №5. С.389-400.

62. Ньютон И. Математические начала натуральной философии / Пер. с лат. и комментарии А.Н. Крылова. Ред. и предисл. Л.С. Полака. М. 1989.

63. Ортега Д., Пул.У. Введение в численные методы решения дифференциальных уравнений. М.: Наука. 1986. 288 е.

64. Питъева Е.В. Новая численная теория движения планет ЕРМ98 и ее сравнение с эфемеридой DE403 Лаборатории реактивного движения США // Труды ИПА РАН. 1998. 3. С. 5-23.

65. Питъева Е.В. Современные численные теории движения Солнца, Луны и больших планет // Труды ИПА РАН. 2003. С. 112-134.

66. Платонов А.К., Власова З.Я., Степанянц В.А. Многоточечный метод интегрирования с переменным шагом для обыкновенных дифференциальных уравнений // М.: 1976. 18 с. (Препринт/ИПМ:72).

67. Птолемей К. Альмагест. М.: Наука. Физматлит. 1998.

68. Пуанкаре А. Избранные труды. Т. I. М.: Наука. 1971; Т. II. М.: Наука. 1972.

69. Пуанкаре А. Лекции по небесной механике. М.: Наука. 1965. 572 с.

70. Радзиевский В.В. Анализ движения тел солнечной системы и их наблюдение. Рига. 1986. С. 126-143.

71. Радзиевский В.В. Кинематика и физика небесных тел. 1987. Т. 3. № 1. С.66-77.

72. Савченко В.В. Использование модифицированного метода Эверхарта при прогнозировании и уточнении орбит комет // Бюлл.ИТА. 1984. 6. С. 324328.

73. Самарский A.A. Математическое моделирование и вычислительный эксперимент // Современные методы математического моделирования. Сборник лекций международной конференции Математическое моделирование. 13-16 июня 2001. Самара. 2001. С. 4-12.

74. Самарский A.A., Гулин A.B. Численные методы М.: Наука. 1989. 432 с.

75. Самарский A.A., Михайлов А.П. Математическое моделирование. М.: Физматлит. 1997.

76. Современные численные методы решения обыкновенных дифференциальных уравнений. Под ред. Дж. Холла, Дж. Уатта. М.: Мир. 1979. 312 с.

77. Соловьев A.C. К проблеме захвата комет на периодические орбиты // На-учн. тр. матем. фак. ВГПУ. Вологда: Русь. 2001. С. 150-169.

78. Субботин М.Ф. Введение в теоретическую астрономию. М.: Наука. 1968. 800 с.

79. Томанов В.П. Кометная космогония. Вологда. 1989. 96 с.

80. Томанов В.П., Кузьмичев В.В. Результаты статистического анализа ко-метных орбит//Волог. гос. пед. ун-т. Вологда. 2003. 43 с.

81. Уинтнер А. Аналитические основы небесной механики. М.: Наука. 1967. 524 с.

82. Уипл Ф.Л. Семья Солнца: Планеты и спутники Солнечной системы. М.: Мир. 1984.316 с.

83. Херик С. Астродинамика. В 3-х книгах. М.: Мир. 1976-1978.

84. Цицин Ф.А. Происхождение комет: новый взгляд на старую проблему // Вселенная и мы. М. 2001. 4. С.21-29.

85. Чазов В.В. Создание численно-аналитической теории движения небесных тел // Околоземная астрономия-2003: сборник трудов конференции. Ин-т астрон. РАН. СПб.: ВВМ 2003. С.171-175.

86. Чеботарев Г.А. Аналитические и численные методы небесной механики. М.,Л.: Наука, 1965. 368 с.

87. Черницов A.M., Батурин А.П. О точности прогнозирования движения комет // Исслед. по баллистике и смеж. вопр. мех. 1998. 2 . С. 144-148.

88. Шарлъе К. Небесная механика. М.: Наука. 1966. 628 с.

89. Штеттер Х.Д. Анализ методов дискретизации для обыкновенных дифференциальных уравнений. М.: Мир. 1978.

90. Bretagnon P. Theorie du movement de l'ensemble des planets. Solution VSOP82//Astron. Astrophys. 1982. 144.

91. Bulirsch R., Stoer J. Fehlerabschatzungen und Extrapolation mit rationalen Funktionen bei Verfahren vom Richardson-Typus // Num. Math. 1964. 6. P. 413-427.

92. Bulirsh R., Stoer J. Numerical Treatment of Ordinary Differential Equations by Extrapolation Methods //Num. Math. 1966. 8. P. 1-13.

93. Butcher J.C. On the convergence of numerical solutions to ordinary differential equations //Math. Сотр. 1966. 20. P. 1-10.

94. Butcher J.C. The effective order of Runge-Kutta methods // In: Lecture Notes in Mathematics Berlin: Springer. 1969. 109.P 133-139.

95. Carusi A., Kresak L., Perozzi E., Valsecchi G.B. Long-term Evolution of Short-period Comets. Adam Hilger. Bristol. 1985.

96. Carusi A., Kresak L., Valsecchi G.B. Electronic Atlas of Dynamical Evolutions of Short-Period Comets.

97. Cohen C.J., Hubbard E.C. An Algorithm Applicable to Numerical Integration of Orbits in Multirevolution steps // Astron. J. 1960/V. 65. P. 454-456.

98. Cowell P.H., Crommelin A.C.D. Essay on the Return of Halleys Comet // Leipzig. P. 1-60.

99. D'Alembert Opuscules mathématiques, ou mémoires sur différents Sujets de Geometrie, de Mechanique d'Optique, d'Astronomie, 6 // Paris. 1773.

100. Dahlquist G. Convergence and stability in the numerical integration of ordinal differential equations // Math. Scand. 1956. 4. P. 33-53.

101. Davis P.J. Interpolation and approximation // Blaisdell. 1963.

102. Eckert W.J., Brouwer D., Clemence G.M. Coordinates of the Five Outer Planets 1653-2060 // Astr. Pap. 1951. 12.

103. Emel'yanenko V. An explicit simplistic integrator for cometary orbits // Ce-lest. Mech. and Dyn. Astron. 2002. 84. 4. C. 331-341.

104. Encke J.F. Uber mechanische Quadratur // Berl. Astr. Jahrb. fur. 1837. P.251-287.

105. Everhart E. Implicit single methods for integrating orbits // Celestial mechanics. 1974. №.10. P.35-55.

106. Fernandes J.A., Tancredi G., Rickman H., Licandro J. The population, magnitudes and size of the Jupiter family comets // Astron. and. Astrophys. 1999. 352. l.P.327-340.

107. Gauss C.F. Exposition d'une nouvelle methode de calculer les perturbations planétaires. (Nachlass) // Werke. 1906. 7. P.439-472.

108. Gear C.W. Hybrid methods for the initial value problem in ordinary differential equations // SIAM J. Numer. Anal, 1964. 2. P. 69-86.

109. Gragg W.B. On extrapolation algorithms for ordinary initial value problems // J.SIAM Numer. Anal.1965. 2. P.384-403.

110. Grigor'yan O.F., Medvedev Yu.D., Tomanov V.P. Non-gravitational effects in Harrington-Abel comet motion due to Jupiter // Proceeding of the Conference Asteroids, Comets, Meteors. Berlin. Noordwijk. ESTEC. 2002. P. 461464.

111. Halley E. Synopsis of the Astronomy of Comets // Phil. Transactions Oxford Univ. 1705.

112. Hasegava I. Catalogue of periodic comets // Mem. College Sci. Kyoto Univ. Series Phys. 1968. 32. P.37-83.

113. Henrici P. Discrete variable methods in ordinary differential equations // Wiley. 1962.

114. Higichi A., Kokubo E., Mikai T. Cometary dinamics: migration due to gas drag and Scattering by protoplanets // Proceeding of the Conference Asteroids, Comets, Meteors. Berlin. Noordwijk. ESTEC. 2002. P.453-456.

115. Hughes D.W. The magnitude distribution and evolution of short-period comets // Mon. Notic. Roy. Astron. Soc. 2002. 336. №2. C.363-372.

116. Isaacson E., Keller H.B. Analysis of numerical methods // Wiley. 1966.

117. Lagrange J.L. Mecanique analytique. Seconde partie, section 7, ch. 2. Sur la variation des elements des orbites elliptiques produites par une force d'impulsion ou par des forces acceleratrices. Oeuvres de Lagrange // Paris. 1812.12. P. 63-101.

118. Lambert J.D. Computational methods in ordinary differential equations // Wiley. 1973.

119. Laplace P.S. Traite de mecanique celeste // Paris. 1805. 4(2). P. 216—228.

120. Le Verrier U.J. Recherches sur les cometes périodiques II A. N. 1848. 26. P.375-384.

121. Lexell A.I. Ulteriores disquisitiones de tempore periodico cometae anno 1770 observati//Acta Acad. Se. Petropol. 1778. 1. P. 317-352.

122. Liniger W., Willoughby R.A. Efficient integration methods for stiff systems of ordinary differential equations // SIAM J. Numer. Anal. 1970. 7. P. 47-66.

123. Marsden B., Williams G.V. Catalogue of Cometary Orbits 1999. 13th ed. Central Bureau for Astronomical Telegrams and Minor Planet Center. 1999.

124. Marsden B.G. Nongravitational forces on comets The first fifteen years // Dynamics of comets, eds. Carusi A. and Valsecchi. Reidel. Dordrecht. P.343-352.

125. Mazeeva O.A., Emel'yanenko V.V. Variations of the Oort cloud comet flux in the planetary region // Proceeding of the Conference Asteroids, Comets, Meteors. Berlin. Noordwijk. ESTEC. 2002. P. 445-448.

126. Newhall X.X., Standish E.M., Williams Jr. and J.G. DE 102: a numerically integrated ephemeris of the Moon and planets spanning forty-four centuries //Astron.Astrophys. 1983. № 125. P. 150-167.

127. Pitjeva E. V. EPM2002 and EPM2002C two versions of high accuracy numerical planetary ephemerides constructed for TDB and TCB time scales // Communication of IAA RAN. 2003. 155. P. 1-19.

128. Pitjeva E.V. Modern numerical ephemerides of planets and the importance of ranging observations for their creation // Celest. Mech. and Dyn. Astron. 2001.80. 3-4. P.294-271.

129. Rickman H., Valsecchi G.B. From the Oort cloud to observable short-periods comets. I. The initial stage of cometary capture // Froeschle CI. Mon. Notic Roy Astron. Soc. 2001. 325. 4. P.1303-1311.

130. Schafer W.A. Semi-analytical methods for computing the orbits of asteroids and comets // Proceeding of the Conference Asteroids, Comets, Meteors. Berlin. Noordwijk. ESTEC. 2002. P.345-349.

131. Shampine L.F., Watts H.A. Block implicit one-step methods I I Math. Comp.1969. 23. P.731-740.

132. Spijker M.N. Convergence and stability of step-by-step methods for the numerical solution of initial value problems //Num. Math. 1966. 8. P.161-177.

133. Standish E.M. DE403/LE403: Annooncement // Jet Prop Lab Technical Report. IOM314. 1995. P. 10-124.

134. Standish EM JPL Planetary and Lunar Ephemerides, DE405/LE405 // Jet Prop Lab Technical Report. IOM 312.F-048. 1998. P. 1-7.

135. Standish E.M. Orientation of the JPL Ephemerides, DE200/LE200, to the Dynamical Eguinox of J2000 // Astron. Astrophys. 114. P.297-302.

136. Stoer J. Uber zwei Algorithmen zur Interpolation mit rationalen Funktionen //Num. Math. 1961. 3. P.285-304.

137. Wang Yu., Ding Kong, Dong Hong, Zhou Hong-han. Nongravitational effect on orbital motions of comets // J. Nanjng. Norm. Univ. Sei. 2000. 23. 4. P.58-61.

138. Whipple F.L. The comet model // Aph. J. 1950. 111. P. 375-394.