
- •А.И. Ходанович математическое моделирование на компьютере
- •История математического моделирования на компьютере
- •Задачи аналитического и численного моделирования
- •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
- •Литература
- •Ходанович
1.4 Сопротивление бесконечной электрической цепи
В силу симметрии добавление или исключение звена в бесконечной цепочке сопротивлений не меняет общего сопротивления цепи:
.
Решая квадратное уравнение
для
,
получаем
(рис.1) [9].
Рис.1.
Физическая симметрия в задаче о
сопротивлении
электрической цепи
1.5 Сопротивление конечной электрической цепи
Сопротивление
цепочки из n
звеньев выражается рекуррентной формулой
,
R1=2R.
(cлайд
4)
Слайд 4. Сопротивление конечной электрической цепи
1.6 «Золотое сечение» в поле тяготения
Известно, что ускорение силы тяжести при удалении от земной поверхности описывается следующей формулой (докажите!):
,
где h- высота над поверхностью Земли, R- ее радиус. При опускании тела в глубь Земли характер зависимости g от расстояния меняется:
.
Зададимся
вопросом, когда
.
Ясно, что одним из решений будет
[26]. Второе решение (которое мы можем
найти на компьютере, написав строчку
Maple-
программы):
(слайд
5).
Слайд 5. Графический метод в физике
Указание. Решить уравнения графическим методом [25]:
2. Динамическое моделирование и вычислительный эксперимент
В
соответствии со вторым законом Ньютона
дифференциальное уравнение движения:
.
Сумма проекций сил в правой части может
быть постоянной, а может зависеть от
координаты, скорости, времени, параметров.
Это дифференциальное уравнение,
которое аналитически легко решается
только в простейших случаях. Численный
приближенный метод основан на замене
дифференциальных уравнений
конечно-раз-ностными алгебраическими
уравнениями с заданной точностью.
Наиболее распространенными
до сих пор являются методы Рунге-Кутты
реализованные в стандартных процедурах
современных математических пакетов.
В 1885 г. К. Рунге высказал основную идею метода численного ре-шения систем дифференциальных уравнений, которую в 1901 г. развил и усовершенствовал В. Кутта. На этой основе была разработана серия вычислительных методов различной степени точности, ориентированных на случай дифференциальных систем [47].
Следует обратить внимание на то, что, в зависимости от содержания, задачи могут решаться в разных системах единиц, а также в безразмерных величинах.
Равенство тяжелой и инертной массы в свободном падении
Масса
является физической величиной,
характеризующая взаимодействие тел,
определяет инертные и гравитационные
свойства. Под инертностью материальных
тел понимается их способность сохранять
состояние покоя или равномерного и
прямолинейного движения. Инерт-ная
масса является мерой этой способности,
т.е. мерой инертности тела. Она определяется
как коэффициент пропорциональности
между силой и ускорением в уравнении
движения классической механики:
.
Инертная масса является также коэффициентом
пропорциональности между импульсом
(количеством движения) тела и его
скоростью:
.
Под
гравитационным свойством материальных
тел понимается их способность
создавать в окружающем пространстве
гравитационное поле
и испытывать воздействие полей тяготения,
создаваемых другими материальными
телами. В этом отношении они аналогичны
электрически заряженным
телам, создающим вокруг себя электрическое
поле и испытывающим
воздействие электрических полей,
создаваемых другими заряженными
телами. Гравитационная масса выступает
как гравитационный заряд.
Согласно закону всемирного тяготения
сила гравитационного взаимодействия
точечных тел пропорциональна произведению
их гравитационных масс m1(g)
и m2(g)
и
обратно пропорциональна квадрату
расстояния
между ними
.
Инерция тел и их способность возбуждать в окружающем пространстве гравитационные поля не могут априори (до опыта) рассматриваться как взаимосвязанные и тем более тождественные свойства тел. Поэтому возникает вопрос, какова взаимосвязь между инертной и гравитационной массами тел? Равны эти массы или пропорциональны друг к другу- эти утверждения не из каких теоретических положений механики не вытекают. Этот вопрос может быть решен только опытным путем.
Слайд 6. Визуальное моделирование
свободного падения с графическими примитивами
Слайд 7. Pascal- программа измерения ускорения свободного падения
Рассмотрим
свободное падение тел в однородном поле
тяготения Земли. Так как тела падают
под действием силы тяжести, то по второму
закону Ньютона m(i)a=m(g)g.
Тогда ускорение, с которым тела падают
на
Землю:
.
Опытным путем можно установить закон, согласно которому все тела свободно падают на Землю с одинаковым ускорением а, равным g=9,8 м/с2, из чего следует, что инертная и гравитационная массы равны. Кроме того, используя компьютер можно измерить время и ускорение свободного падения, а также проиллюстрировать исторический закон нечетных чисел Галилея.