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

• Заусаев, Артем Анатольевич
• кандидат физико-математических наук
• 2005

Оглавление диссертации кандидат физико-математических наук Заусаев, Артем Анатольевич

ВВЕДЕНИЕ

ГЛАВА 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 гг. в электронном и печатном вариантах.

ЗАКЛЮЧЕНИЕ

Список литературы диссертационного исследования кандидат физико-математических наук Заусаев, Артем Анатольевич, 2005 год

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.