Наближене формування розподілів неперервних випадкових величин
У такому разі зручно використовувати наближений спосіб перетворення РВП [0, 1] на випадкові числа з іншим розподілом. Розглянемо сутність цього способу.
Залежність щільності розподілу від можливих значень випадкової величини х зобразимо графічно на відрізку, де х змінюється від a до b Якщо межі інтервалу змінювання випадкової величини нескінченні, то початковий розподіл зрізуємо із заданою точністю. Зокрема, шукана графічна залежність може бути й експериментальною.
Розіб’ємо відрізок [a, b] на n частин (рис. 8.1) таких, що
= = . . . = =
де — координата точки розбиття.
Рис. 8.1. Графік апроксимації функції щільності розподілу
З урахуванням цього ймовірність випадкової події, яка полягає в тому, що випадкова величина X потрапить в один з інтервалів, подається у вигляді
випадкова точка може потрапити на будь-який відрізок з однаковою ймовірністю.
Функцію щільності f (x) апроксимуємо східчастою функцією так, щоб значення f (x) у кожному інтервалі були сталою величиною.
Координату випадкової точки М, яка потрапила на і-й інтервал, можна подати у вигляді
(8.10)
де r — відстань точки М від лівого кінця інтервалу.
З огляду на лінійну апроксимацію величина rє рівномірно розподіленою випадковою величиною на відрізку . Цю величину можна дістати за допомогою формули (8.4)
.
Номер інтервалу і, в якому міститься випадкова точка М, можна визначити, скориставшись схемою випробувань за «жеребкуванням» для рівноймовірних подій, що утворюють повну групу.
Тому алгоритм пошуку наступного випадкового числа, яке має розподілf (x) (після виконання попередніх розрахунків щодо апроксимації початкового розподілу), може бути такий.
1. Генеруємо РВП [0, 1].
2. За допомогою знаходимо індекс
(8.11)
3. За допомогою знаходимо
.
4. Відшукуємо наступне випадкове число
. (8.12)
Отже, щоб утворити випадкове число х, необхідно витратити два числа РВП [0, 1].
Під час вивчення цієї теми особливу увагу слід приділити генеруванню нормально розподілених випадкових чисел, оскільки нормальний (гауссів) розподіл — один з найважливіших і найчастіше застосовуваних видів неперервних розподілів.
Щільність розподілу ймовірностей для нормального розподілу