9.4.2. Генераторы случайных чисел

С функцией-генератором случайных чисел Random связано много условностей. Она всегда возвращает равномерно распределенное случайное число. Диапазон или, вернее, интервал распределения таких чисел задается видом вызова функции. Если Random вызывается без аргумента, то она вернет вещественное случайное число r в диапазоне 0 <= r < 1. Но если вызов содержал аргумент N типа Word — Random( N ), то вернется случайная величина (целая!) из диапазона 0 <= r < N. При N=0 функция вырождается и возвращает константу 0.

У функции-генератора случайных чисел Random есть одна неприятная особенность: при последовательных запусках программы выдавать одинаковые случайные последовательности. Другими словами, пока программа работает, любое обращение к Random даст вполне «случайное» значение (период повторяемости, конечно, есть, но очень большой). Но если запустить программу второй раз, то ее вычисления в точности повторят предыдущий запуск: «случайные» значения будут теми же, что и в прошлый раз. А это значит, что многократный запуск программ, моделирующих случайные процессы наберет статистику, одинаковую с единичным прогоном. Во избежание этого рекомендуем включать в начало программы вызов процедуры Randomize. Эта процедура записывает случайное число (взятое со встроенного таймера) в так называемую «затравку» случай-

- 180 -

ной последовательности, сбивая тем самым ее на новые числа. Того же эффекта можно достичь, если записывать произвольное (но каждый раз разное!) значение в системную переменную RandSeed типа LongInt.

Функции Random генерируют случайную равномерно распределенную на интервале последовательность чисел. На ее основе можно сконструировать последовательности с другими законами распределения. Наиболее часто требуется нормальное распределение. Один из вариантов генератора показан на рис. 9.7.

| { Функция Gauss возвращает случайное вещественное число }

| {распределенное по нормальному закону с математическим }

| {ожиданием mo и средне-квадратическим отклонением Sigma}

| FUNCTION Gauss( mo, Sigma : Real ) : Real;

| VAR

| a, b, r, sq : Real;

| BEGIN

| repeat

| a := 2 * Random - 1;

| b := 2 * Random - 1;

| r := Sqr(a) + Sqr(b)

| until ( r < 1 );

| sq := Sqrt( -2 * Ln(r)/r);

| Gauss := mo + Sigma*a*sq

| END;

Рис. 9.7

Функцию Gauss несложно модифицировать, чтобы она возвращала целые значения.