- •А.И. Ходанович математическое моделирование на компьютере
- •История математического моделирования на компьютере
- •Задачи аналитического и численного моделирования
- •1. Симметрия и модели с числовыми последовательностями
- •1.1 Формула пути равнопеременного движения
- •1.3 «Золотое сечение» и числа Фибоначчи
- •1.4 Сопротивление бесконечной электрической цепи
- •1.5 Сопротивление конечной электрической цепи
- •1.6 «Золотое сечение» в поле тяготения
- •Равенство тяжелой и инертной массы в свободном падении
- •2.2 Реальная динамика тела в среде
- •2.3 Цепная линия и метод кратных скобок Пуассона
- •2.4 Баллистическая кривая
- •2.5 Модели космической динамики
- •2.6 Посадка «Гюйгенса» на Титан
- •Задача Кеплера
- •Слайд 10. Компьютерно-графическое моделирование в небесной механике
- •2.8 Орбита в пространстве скоростей
- •2.9 Конические сечения в поле тяготения
- •2.10 Прецессия орбиты при малом возмущении
- •2.11 Задача трех тел в небесной механике
- •2.12 Задача о брахистохроне
- •Слайд 11. Иллюстрация экстремальной траектории в поле тяготения
- •2.13 Управление таймером в режиме реального времени
- •2.14 Физический опыт и модель измерений
- •2.15 Водяные часы
- •2.16 Фракталы в комплексной плоскости
- •Слайд 14. Комплексные отображения в нелинейной динамике
- •3. Колебания и адиабатические инварианты
- •3.1 Гармонический осциллятор
- •3.2 Задача «Бездонный колодец»
- •3.3 Фигуры Лиссажу
- •3.4 Параметрические колебания
- •3.5 Маятник Капицы
- •3.6 Хаотическое поведение маятника и отображение Пуанкаре
- •3.8 Задача с соударениями
- •3.9 Хаотические колебания при соударениях
- •3.10 Представление колебаний рядом Фурье
- •3.11 Два тела на пружине
- •3.12 Цепочка Ферми-Паста-Улама
- •3.13 Линейная цепочка связанных осцилляторов
- •3.14 Колебания мембраны в интерактивной графике
- •3.15 Вращение частицы в трехмерной графике
- •4. Вероятностно-статистическая линия математического моделирования
- •4.1 Нормальное распределение и статистический критерий
- •4.2 Распределение расстояний между молекулами газа
- •4.3 Радиоактивный распад изотопа
- •4.4 Флуктуации и биномиальный случай
- •4.5 Статистическое распределение с бесконечной дисперсией
- •4.6 Закон больших чисел и центральная предельная теорема
- •4.7 Статистические усреднения в эксперименте
- •4.8 Диффузия в модели случайных блужданий
- •4.9 Иерархия временных масштабов в броуновской динамике
- •4.10 Метод Монте-Карло с генератором случайных чисел
- •4.11 Алгоритмы числа
- •4.12 Модели плоских фигур
- •4.13 Площадь круга с малым возмущением
- •4.14 «Игла Бюффона»
- •4.15 Выделение сигнала на фоне «шумов»
- •4.16 Регрессионная модель эксперимента и метод наименьших квадратов
- •4.17 Модели линейной регрессии в планировании эксперимента
- •5. Дидактические игры математического моделирования
- •5.1 «Телепат»
- •5.2 Игрушка «Раскидай»
- •5.3 Игрушка «Ванька-встанька»
- •5.4 Игра Баше
- •5.5 «Перевернутая» игра Баше
- •5.6 Игра в 15
- •Литература
- •Ходанович
2.9 Конические сечения в поле тяготения
Исследуйте «открытые гиперболические орбиты». Определите радиус, координаты центра и центральный угол кругового годографа вектора скорости, а также зависимость угла между асимптотами гиперболической траектории от начальной скорости.
Рис.11. Параболическая траектория в поле тяготения
Указание. Общие свойства конических сечений приведены в [33]. Исследовать траектории в поле тяготения при различных начальных условиях. Построить графики конических сечений.
2.10 Прецессия орбиты при малом возмущении
Вызванное искажение гравитационного поля Земли в результате суточного вращения приводит к прецессии орбиты спутника, т.е. к срав-нительно медленному вращению эллиптической орбиты с фокусом в центре Земли вокруг этого фокуса. Из законов сохранения энергии и момента импульса следует, что максимальное и минимальное расстояния спутника от Земли не изменяется от витка к витку при прецессии его орбиты.
Искажение центрального поля тяготения, вызванное вытянутостью вдоль оси, описывается малым дополнительным членом обратно пропорциональным четвертой степени расстояния. Отклонения от закона обратных квадратов для силы притяжения могут быть учтены малой поправкой d в уравнении: F~1 / r2+d.
У Земли реальная сплюснутость мала, поэтому большая ось эллиптической орбиты за один поворот спутника поворачивается на очень малый угол, однако этот эффект накапливается за большое число оборотов по орбите.
Рис.12. Прецессия орбиты при малом гравитационном возмущении
Указание. Для большей наглядности начальные условия нужно задавать так, чтобы орбита спутника была сильно вытянутой (начальная скорость должна заметно отличаться от круговой).
2.11 Задача трех тел в небесной механике
При изучении астрофизики необходимо отметить, что существующая иерархия масс небесных тел позволяет в ряде случаев рассматривать квазизамкнутые системы, состоящие из части взаимодействующих (гравитирующих) тел, обладающих наибольшими массами, поскольку на них наименее массивные объекты оказывают незначительные воздействия. А лишь после того, как поступательные движения массивных тел (например, больших планет Солнечной системы) определены, удается исследовать отдельную задачу о движении менее массивных тел (астероидов, комет) уже при условии, что силы, управляющие их движениями, будут известны.
В связи с этим фундаментальную основу небесной механики составляет задача N гравитирующих материальных точек (N тел), состоящая в определении по начальным положениям и скоростям характера движения в трехмерном евклидовом пространстве изолированной системы, состоящей из N материальных точек с фиксированными массами, притягивающих друг друга по закону всемирного тяготения [15].
Но необычайная сложность поиска общего решения уже задачи 3-х тел, математическая модель которой представляет собой систему нелинейных дифференциальных уравнений, привела к тому, что именно эта знаменитая “ньютонова проблема трех тел” на протяжении трех веков играла определяющую роль в развитии математических методов и в возникновении новых плодотворных направлений в математике, небесной механике, вычислительной физике.
Особенность многомерных канонических систем с числом степеней свободы больше двух связана и с явлением динамического хаоса (стохастичностью), когда траектории системы могут представлять собой реализацию случайных временных процессов, несмотря на то, что в системе непосредственно отсутствуют какие-либо внешние случайные силы. Главная черта хаотических систем состоит в том, что малое возмущение начальных условий приводит за конечное время к непредсказуемости результирующего движения. Динамический хаос обусловлен существованием специфической локальной неустойчивости относительно малых возмущений орбит системы, он присущ только нелинейным системам и непосредственно связан с их неинтегрируемостью.
Рис.13. Хаотическая траектория в задаче трех тел
Орбиты можно разделить на устойчивые и неустойчивые. Устойчивыми орбитами могут быть открытые петли, охватывающие обе звезды, орбиты в виде восьмерок или кеплероподобные орбиты вокруг только одной звезды. В случае неустойчивых (хаотических) орбит возможно падение на одну из звезд [46].
Указание. В динамических уравнениях необходимо учесть все силы, действующие в системе многих тел и модифицировать программу вычислительного эксперимента (слайд 10) .