- •А.И. Ходанович математическое моделирование на компьютере
- •История математического моделирования на компьютере
- •Задачи аналитического и численного моделирования
- •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
- •Литература
- •Ходанович
Задача Кеплера
|
В результате длительной обработки многолетних наблюдений датского астронома Тихо Браге (1546-1601) немецкий ученый Кеплер эмпирически установил три закона планетных движений. Эти законы формулируются следующим образом:
|
1) каждая планета движется по эллипсу, в одном из фокусов которого находится Солнце;
2) радиус-вектор планеты в равные времена описывает равные площади;
3) квадраты времен обращений планет относятся как кубы больших осей эллиптических орбит, по которым они движутся вокруг Солнца.
Первые два закона были опубликованы Кеплером в 1609 г., последний- в 1619 г. Законы Кеплера естественным путем привели Ньютона к открытию закона всемирного тяготения.
Из первого закона Кеплера следует, что траектория планеты- плоская кривая. С учетом этого обстоятельства, из второго закона Кеплера следует, что сила, заставляющая планету двигаться по замкнутым орбитам, направлена к Солнцу.
При рассмотрении планетных движений можно не учитывать движение Солнца, считая его массу бесконечно большой по сравнению с массой планеты. При формальном учете массы солнца динамика такова, как если бы Солнце оставалось неподвижным, но гравитационная постоянная увеличилась в (1 + т/М) раз. Поэтому для относительного движения первый и второй законы Кеплера остаются справедливыми. Зато третий закон должен быть уточнен. Для этого достаточно гравитационную постоянную G заменить на G(1 +m/M). Это приводит к соотношению . То обстоятельство, что для планет Солнечной системы закон выполняется с большой точностью, связано с тем, что масса планеты очень мала по сравнению с массой Солнца.
Динамические уравнения движения массы в поле тяготения имеют вид:
После несложных преобразований система уравнений приводится к виду в декартовых координатах:
Рис.10. К задаче Кеплера
Указание. Исследовать движения планет Солнечной системы и законы Кеплера в задаче двух тел. Построить траектории.
Не умаляя общности решения, можно задать следующие начальные условия (t=0): x=a, y=0, vx=0, vy=v(0). Начальную скорость удобно задавать в единицах скорости, соответствующей круговой орбите. Начальная скорость больше круговой (из интервала от 1 до 1.4) и начальная точка траектории является перигеем.
Слайд 10. Компьютерно-графическое моделирование в небесной механике
2.8 Орбита в пространстве скоростей
Линия, соединяющая концы всех векторов скорости, проведенных из одной точки, называется годографом вектора скорости (траектория в пространстве скоростей). Кеплерово движение по замкнутым траекториям обладает интересным свойством: для любых эллиптических траекторий годограф вектора скорости всегда представляет собой окружность с радиусом u=GmM/L (L- величина углового момента) [16].
Для круговой орбиты центр этой окружности совпадает с общим началом всех векторов скоростей. Для эллиптической орбиты центр окружности смещен из общего начала и расположен на векторе скорости в перигее орбиты. Данное свойство кеплерова движения позволяет находить значение скорости в любой точке эллиптической траектории, если известны скорости в перигее и апогее.
Указание. Постройте орбиту в пространстве скоростей (Vx,Vy). Убедитесь в том, что в каждый момент времени радиус-вектор перпендикулярен скорости u. Объясните этот факт.