
- •А.И. Ходанович математическое моделирование на компьютере
- •История математического моделирования на компьютере
- •Задачи аналитического и численного моделирования
- •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
- •Литература
- •Ходанович
4.2 Распределение расстояний между молекулами газа
Пространство заполняет газ. Вероятность встретить молекулу газа внутри бесконечно малого объема dv равна adv. Для любой молекулы в любой момент времени найдется какая-то молекула- ближайший сосед. Расстояние до ближайшего соседа есть случайная величина. В разные моменты времени она различна. Найти плотность распределения расстояния молекулы до ближайшей частицы [1].
Для
того чтобы ближайший сосед находился
на расстоянии, заключенном между r
и r
+ dr,
необходимо, 1) чтобы расстояние до
ближайшего соседа было не меньше r-
вероятность этого события равна
;
2) чтобы внутри сферического слоя 4πr2dr
была молекула- вероятность этого события
равна a4πr2dr.
Так как эти два события взаимно независимы,
то вероятность, что произойдет и то,
и другое, равна произведению их
вероятностей. Следовательно,
.
Дифференцируя по r
и имея в виду
,
получаем
.
Деля обе части на
и интегрируя, находим
.
Постоянная интегрирования находится
из условия нормировки, т.е.
.
Из
полученного выражения следует, что при
возрастании r
плотность вероятности расстояния
до ближайшего соседа растет от нуля,
достигает максимума при
,
а затем убывает, стремясь к нулю.
4.3 Радиоактивный распад изотопа
Вероятность
распада ядра не зависит от того, как
долго ядро существует, не распадаясь.
Какова вероятность распада изотопа в
момент времени t
? Найти
зависимость между постоянной распада
и временем полураспада Т1/2.
Найти плотность распределения времени
распада радиоактивного изотопа и среднее
время распада.
Если
-
вероятность распада в единицу времени,
то вероятность того, что ядро не
распадается за время dt,
равна 1 -
.
В промежутке времени t
промежуток времени dt
содержится t/dt
раз. Вероятность того, что в момент
времени t
изотоп распадется, равна, согласно
теореме умножения вероятностей,
произведению t/dt
множителей (вероятностей того, что
изотоп не распадется до данного момента
времени), т.е.
.
Считая
dt
бесконечно малым, получим после
предельного перехода
.
Полагая
,
приходим к закону радиоактивного распада
.
Если
вероятность того, что изотоп распадется
за время Т1/2,
равна 1/2, то Т1/2
называется временем полураспада.
Таким образом, находим
.
Вероятность
того, что распад произойдет в интервале
времени
[t
+ dt], равна
вероятности того, что изотоп не распадается
за время t,
умноженной на вероятность того, что он
распадается за следующий промежуток
времени dt.
Следовательно,
и
,
а среднее время распада изотопа
.
4.4 Флуктуации и биномиальный случай
Любое физическое измерение производится с определенной точностью. Причем точность измерения, как правило, определяется конструкцией прибора и методикой измерений. Но существуют принципиальные физические ограничения точности измерений связанные с флуктуацией физических величин [28].
Например, макроскопические параметры термодинамической системы в тепловом равновесии не являются строго постоянными, а испытывают малые беспорядочные колебания вблизи некоторых средних значений. Флуктуации- это хаотические отклонения физических величин от их средних значений, вызываемые тепловым движением образующих систему частиц [2].
Слайд 25. Асимптотика биномиального распределения вероятностей для газа с n=100
Основные закономерности флуктуаций видны на примере пространственного распределения молекул идеального газа внутри сосуда в состоянии термодинамического равновесия.
Вероятность
того, что в половине сосуда находится
m
частиц из полного их частиц n,
можно найти в предположении, что
вероятность нахождения одной молекулы
в определенной половине сосуда равна
½ и не зависит от положения других
частиц. В этом случае применима модель
независимых случайных испытаний с
биномиальным распределением вероятностей
(или схемой Бернулли)
,
где q=1-p.
Заметим,
что при больших n
пользоваться формулой Бернулли
затруднительно, поэтому удобнее
воспользоваться асимптотикой Лапласа
,
где
,
а
,
т.е. огибающая распределения переходит
в нормальное распределение Гаусса
(слайд 24).
Вероятность
того, что в n
независимых испытаниях событие появится
от m1
до m2
раз определяется интегральной теоремой
Лапласа
,
где
,
а
[17].
Флуктуация
физической величины характеризуется
стандартным отклонением распределения,
т.е. корнем из дисперсии. Причем, с
увеличением n
флуктуация растет
,
а относительная флуктуация убывает
.