
- •1.Компоненты компьютерного эксперимента.
- •2.Интегральный закон распределения вероятностей.
- •3.Дифференциальный закон распределения вероятностей.
- •4.Сложение случайных чисел и центральная предельная теорема.
- •5. Розыгрыш дискретной случайной величины
- •6. Метод обратной функции
- •7.Метод Неймана
- •8. Метод суперпозиции
- •9. Метод Батлера
- •12. Метод Монте-Карло.
- •13.Общая схема ммк.
- •15.Эффективность и погрешность ммк.
- •16.Задача переноса частиц в среде
- •17. Типы взаимодействий элементарных частиц
- •19.Электрон-фотонный ливень
- •20 .Схема использования ммк в задаче переноса.
- •18. Превращения взаимодействующих частиц в среде за счет электромагнитных взаимодействий.
- •23. Моделирование комптоновского рассеяния.
- •22.Моделирование типа дискретного взаимодействия.
- •10. Модифицированный метод суперпозиции
- •25. Системы отсчета и моделирование характеристик процессов столкновений
- •11. Моделирование специальных распределений
- •21. Моделирование пробега электронов и позитронов.
- •24.Сечение процессов взаимодействия элементарных частиц.
- •Вопросник
- •Компоненты компьютерного эксперимента.
- •Интегральный закон распределения вероятностей.
- •Дифференциальный закон распределения вероятностей.
4.Сложение случайных чисел и центральная предельная теорема.
Пусть
имеется N
одинаковых независимых случайных
величин
,
так что распределения вероятностей
этих величин совпадают (они могут быть
как дискретные, так и непрерывные).
Значит их математические ожидания и
дисперсии совпадают:
M
= M
=
… =M
=
m
D = D =… = D = b2
Введем случайную величину pN, которая является суммой случайных величин
pN=
Так как для произвольных независимых случайных величин ξ и ν выполняются соотношения
M(ξ+ ν)=M ξ+M ν и D(ξ+ ν)=D ξ+D ν
то
M pN=M( )=Nm
D pN=D( )=N b2
Теорема: При больших N (т.е. при N→∞) случайная величина pN описывается нормальным распределением с математическим ожиданием равным Nm и дисперсией D= N b2
т.е. плотность распределения
где
a= Nm,
Следовательно, вероятность найти случайную величину pN в пределах от α до β:
5. Розыгрыш дискретной случайной величины
Дискретная
случайная величина χ, как правило,
задается законом распределения в виде
таблицы
,
где x1,
x2,
. . . , xn
- значения случайной величины χ, а p1,
p2,
. . . , pn
– вероятности появления этих значений.
Другими словами Prob(χ = xi)
= pi.
Наша задача состоит в воспроизведении
значений x1,
x2,
. . . , xn
случайной величины χ, с вероятностями
появления равными p1,
p2,
. . . , pn
соответсвенно. Алгоритм розыгрыша
дискретной случайной величины
1. Необходимо разбить интервал (0,1) на n - интервалов, длины которых равны p1, p2, ..., pn соответственно. 2. Разыгрываем значение стандартной случайной величины γ. 3. Проверяем условие попадания значения γ в i-тый интервал длинной pi: p1 + p2 + . . . + pi−1 < γ < p1 + p2 + . . . +pi. Если неравенство верно, то приписываем случайной величине χ значение, соответствующее данному интервалу т.е. xi. 4. При необходимости повторяем шаги 1 − 3 несколько раз с помощью
новых
значений стандартной величины γ для
получения других значений дискретной
случайной величины с законом распределения
вероятностей.На рисунке представлено
графическая иллюстрация алгоритма
розыгрыша дискретной величины. Для
экономии времени выгодно расположить
значения дискретной случайной величины
xi,
так чтобы выполнялось соотношение p1
< p2
< ... < pn.
6. Метод обратной функции
Метод обратной функции основан на теореме:
Если
случайная величина ξ имеет плотность
распределения вероятности p(x) > 0 на
интервале (a, b) (a < x < b), то случайная
величина ξ удовлетворяющая уравнению
F(ξ) = γ (1), где F (x) функция распределения
F(x) =
,
а γ стандартная равномерно распределенная
величина наинтервале (0, 1) величина
имеет плотность распределения p(x).
Как следует из этой теоремы, с помощью генератора случайных чисел γ уравнениеγ = F(x)
=
(2)позволяет
воспроизвести значения случайной
величины ξ
с плотностью распределенияp(x)
на интервале (a,
b)
путем обращения(инверсии) уравнения
(2):x
=
(γ)
(3).
Такой метод получения называют методом обратной функции(или методом инверсии).