Rambler's Top100Astronet    
  по текстам   по форуму  внутри темы 		 

args[0]=message
args[1]=DB::DB::Message=HASH(0x44f1300)
Re[2]: Вопросы преподавания астрономии
21.07.2007 21:28 | Дмитрий Доценко

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

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

Тест 1: Точечная масса + безмассовая частица вокруг нее на круговой орбите. В конце (например, >>1000 орбит) измерить параметры орбиты (6 штук) и сравнить с начальными значениями. Если поменялась энергия - меняйте схему интегрирования. Если орбита больше не сферическая - то же или что-то с координатами. Если орбита меняет свою плоскость - что-то с преобразованиями координат. И т.д.

Тест 2: Две точечные массы (Звезда + планета). В приведенных координатах должно быть то же самое. Так как отношение масс порядка 10^10 (а можно поставить и больше), то начинает играть роль конечная точность вычислений. Если все ОК, сделайте орбиту вытянутой и повторите.

Тест 3. Взять более сложную задачу и сравнить результат с признанной программой (например, Mercury: http://star.arm.ac.uk/~jec/home.html). Объяснить все несовпадения.

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

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

а) Уберите все планеты, кроме Меркурия, "выключите" ОТО и посмотрите, какой стала орбита. Должна быть замкнутой и следовать законам Кеплера.

б) "Включите" ОТО (если в программе она есть в виде коррекций к закону притяжения). Воспроизводите ли вы аналитическое решение уравнений ОТО?

Если все в порядке, пишите (какие тесты провели, какие результаты, какие отклонения от теории), продолжим дискуссию.

Теперь по вашему вопросу. Насколько мне известно, сильных периодических возмущений орбиты Меркурия нет. Их просто некому вызвать.

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



Форумы >> Вопросы преподавания астрономии
Список  /  Дерево
Заголовки  /  Аннотации  /  Текст
Астронет | Научная сеть | ГАИШ МГУ | Поиск по МГУ | О проекте | Авторам

Комментарии, вопросы? Пишите: info@astronet.ru или сюда

Rambler's Top100 Яндекс цитирования