- •А.И. Ходанович математическое моделирование на компьютере
- •История математического моделирования на компьютере
- •Задачи аналитического и численного моделирования
- •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
- •Литература
- •Ходанович
3.8 Задача с соударениями
Рассматриваются упругие соударения шарика с равномерно движущейся стенкой. Шарик упруго отскакивает между двумя стенками. В случае неподвижных стенок скорость при отскоке будет одинаковая. Если одна из стенок движется, то скорость шарика после отскока от нее будет изменяться, в зависимости от направления движения стенки. При сближении стенок расстояние между ними уменьшается, скорость отскока- увеличивается, а произведение I=Lv – практически cохраняется (с точностью до 0,5%), т.е. I = inv.
Слайд 20. Иллюстрация точности адиабатического инварианта
Указание. В программу включить цикл по числу соударений i=1...N для следующей системы уравнений:
Знак в уравнениях зависит от направления движения стенки, причем, U<<V. Параметры задачи: N=1000, Lo =10 м, Vo =5 м/с, U=0.01 м/с .
3.9 Хаотические колебания при соударениях
Рассматривается нелинейная математическая модель, в которой шарик подскакивает на колеблющейся опоре (модель ускорения частиц в электромагнитном поле). Сделав некоторое предположение о потерях энергии при каждом соударении, можно получить следующие разностные уравнения (двумерное отображение):
Рис.25. Модель Ферми ускорения космических лучей
Здесь ф- фаза опоры в момент соударения, v- безразмерная скорость после соударения, к- потери энергии при соударении, величина А пропорциональна амплитуде колеблющейся опоры [45].
Случай консервативного хаоса исследован Лихтенбергом и Либер-маном как модель ускорения электронов в электромагнитных полях. По решениям системы уравнений временной контур координаты гравитирующей массы восстанавливается по законам кинематики. В эксперименте необходимо варьировать начальные условия.
Слайд 21. Гистограмма измерений в динамике хаотических колебаний
3.10 Представление колебаний рядом Фурье
Рассматриваются колебания шарика при упругих отскоках от неподвижной опоры (слайд 22).
Слайд 22. Разложение периодической функции в ряд Фурье
Указание. Определите период колебаний по формулам кинематики и постройте график координаты от времени, используя разложение функции в ряд Фурье. Данные ввести самостоятельно.
3.11 Два тела на пружине
Решение в общем виде допускает задача о движении системы, состоящей из двух взаимодействующих частиц. Рассмотрим наиболее простую систему двух тел (материальных точек) соединенных легкой пружиной.
В качестве предварительного шага к решению этой задачи покажем, каким образом она может быть упрощена путем рассмотрения движения системы частиц относительно центра масс с учетом обобщенного принципа инерции. Другими словами задачу двух тел можно свести к задаче движения одного тела приведенной массы в потенциальном поле.
Потенциальная энергия взаимодействия двух частиц зависит лишь от расстояния между ними. Поэтому функция Лагранжа системы имеет вид .
Введем расстояние между точками и поместим начало координат в центре масс, что дает . Из двух последних равенств находим: ; . Нетрудно видеть, что функция Лагранжа , где , - приведенная масса.
Функция L формально совпадает с функцией Лагранжа одной материальной точки с приведенной массой m, движущейся в потенциальном поле, в нашем случае, , симметричном относительно неподвижного начала координат. Поэтому собственная частота колебаний материальных точек определяется по известной формуле пружинного маятника .
Рис.26. Колебания в формализме Лагранжа
Неизменность законов механики по отношению к линейным преобразованиям координат и времени в ИСО предполагает неизменность вида функций от координат, скоростей и ускорений, фигурирующих в уравнениях, т.е. инвариантность к заданному преобразованию, например, к сдвигу координаты (однородность пространства). Такие уравнения называются ковариантными относительно группы преобразований. Примером может служить система динамических уравнений в рассматриваемой задаче двух тел на основе дифференциальных законов Ньютона: , где L- длина ненагруженной пружины. Результат численного интегрирования системы уравнений для равных масс на компьютере приведен на рис. 26. Заметим, что при более сложных преобразованиях исходной системы отсчета, например, при переходе к полярным координатам ковариантность уравнений движения в форме Ньютона утрачивается.
Указание. Провести вычислительный эксперимент при разных параметрах и начальных условиях. Проверить теорему о движении центра масс.