- •А.И. Ходанович математическое моделирование на компьютере
- •История математического моделирования на компьютере
- •Задачи аналитического и численного моделирования
- •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.14 Физический опыт и модель измерений
Проделайте следующий физический эксперимент. Установите на полу вертикально палку (можно линейку) и, плавно, слегка наклонив, без толчка отняв руку, отпустите ее. Палка упадет. Измерьте время падения, используя системный таймер компьютера, и сравните результат с теорией.
Рассмотрим физическую модель задачи, т.е. идеализируем объект рассмотрения: вместо палки рассмотрим стержень, не имеющий толщины, установленный вертикально на горизонтальной поверхности; стержень находится в поле тяжести, точка опоры фиксирована, а среда не оказывает сопротивления движению (рис. 15).
Рис.15. Геометрическая иллюстрация физической модели
Из
закона сохранения энергии при вращении
стержня
для
угловой скорости стержня имеем
.
Учитывая
,
по теореме о производной обратной
функции приходим к дифференциальному
уравнению
.
Учитывая начальные условия, для времени
падения имеем
[29].
Заметим, что при нулевой начальной скорости для определения времени падения необходимо численное интегрирование
,
где момент инерции
.
Слайд 13. Численное интегрирование в среде Maple
2.15 Водяные часы
Требуется найти форму водяных часов, при которой уровень воды убывает в чаше с постоянной скоростью. Изобразить часы средствами компьютерной графики.
Известно, что древние водяные часы представляли собой чашу, из которой через небольшое отверстие на дне вытекала вода. Такие часы использовались в греческих и римских судах для хронометрирования речей адвокатов. Приводимая задача из гидродинамики иллюстрирует достаточно распространенный дифференциальный метод, а также связь между физическим содержанием и геометрией [3].
Предварительно остановимся на некоторых общих теоретических рассуждениях. Итак, рассмотрим сосуд (рис. 16), площадь горизонтального сечения которого является произвольной функцией расстояния сечения от дна сосуда.
Пусть высота уровня жидкости в сосуде в начальный момент времени t = 0 равна h метров. Пусть, далее, площадь сечения на высоте х равна S(x), а площадь отверстия на дне сосуда есть s.
И
звестно,
что скорость истечения жидкости v
в тот момент, когда высота ее уровня
равна х,
определяется равенством
,
k-
коэффициент скорости истечения жидкости
из отверстия, зависящий от формы
отверстия.
Рис.16. Дифференциальный метод в гидродинамике
На физически малом промежутке времени dt истечение жидкости можно считать равномерным, а потому за время dt вытечет столбик жидкости, высота которого vdt и площадь сечения s, что в свою очередь вызовет понижение уровня жидкости в сосуде на –dx.
В
соответствии с уравнением неразрывности
приходим к дифференциальному уравнению
,
которое можно переписать в виде
.
Задача
решается с помощью выведенного выше
уравнения, которое перепишем в виде
.
Именно, учитывая, что чашу можно
рассматривать как поверхность вращения,
в соответствии с обозначениями рис. 16
из последнего уравнения получаем, что
,
где
-
проекция скорости свободной поверхности
жидкости на ось х, которая по условию
задачи есть величина постоянная. Возведя
обе части уравнения в квадрат, приходим
к уравнению
,
где
.
Последнее означает, что форма поверхности
водяных часов получается вращением
кривой
вокруг оси х.
