Генерування випадковості чисел

6*. , де

Фортран-програму можна викликати оператором

CALL UNIFRM (A,B,X)

Значення А і В задаються на вході

SUBROUTINE UNIFRM (A,B,X)

1 CALL RANDUM (IX, IY, RN)

2 X = A+(B-A)*RN

3 RETURN

4 END

Методи генерування псевдовипадкових чисел, розподілених по заданому закону.

Найбільш поширеним є методи:

1)зворотної функції;

2)табличний;

3)метод, що базується на функціональних особливостях генерованих розподілів.

1)Метод полягає в використанні наступної теореми:

Якщо випадкова величина має густину розподілу f(z), розподіл випадкової величини

3. є рівномірним в інтервалі 0…1 (а – нижня границя діапазону зміни випадкової величини z).

Для того щоб отримати випадкове число zj з щільністю (густиною) f(z), необхідно знайти рішення рівняння

4.

Наприклад, для експоненційного закону

4*. - функція щільності для експоненційного закону.

Використовуємо формулу (див. ст. 53, рівн. 3

від лівої до правої частини візьмемо ln.

Звідси -

5.