3.3.1. Моделирование случайных последовательностей с ПРВ семейства распределений Пирсона

Каждому конкретному распределению из семейства распределений Пирсона, которое описывается с помощью ДУ

                                              (3.10)

соответствует свой набор коэффициентов . Значит, задавая конкретное распределение, мы тем самым задаем и коэффициенты.

В качестве примера рассмотрим гамма-распределение, для которого ПРВ запишется в виде [44]

,                                 (3.11)

где α и β - параметры распределения, Г() - гамма-функция [7]. Вычисляя производную в (3.11), получим

                                                       (3.12)

Сравнивая (3.11) и (3.10), найдем   .

Рекуррентное разностное уравнение запишется в следующем виде [44]

                (3.13)

Учитывая слагаемые более высокого порядка малости, вместо (3.13) получим

.     (3.14)

Для конкретных значений  и  величина  может быть выражена через них, т.е.  и уравнение (3.13) будет содержать только один параметр .

Рис. 3.2.  График теоретической ПРВ (сплошная линия) и ПРВ, полученной с помощью алгоритма (3.14) при =0.001

Рис. 3.3.  График теоретической ПРВ (сплошная линия) и ПРВ, полученной с помощью алгоритма (3.14) при =0.003

Для экспериментальной проверки степени совпадения ПРВ последовательности (3.14) с гамма-распределением было проведено статистическое моделирование синтезированного алгоритма. На рис. 3.2 и 3.3 представлены теоретическое распределение и распределения, полученные в результате моделирования. Весь диапазон изменения был разбит на 100 одинаковых интервалов. Экспериментальные распределения получены на основе независимой выборки объемом 100000. Проведена проверка гипотез о принадлежности экспериментальных распределений гамма-распределению с помощью критерия . Для заданного числа интервалов (100) и выбранного уровня значимости () критическое значение =140. При значениях <0.003 все статистики  оказались меньше критического значения. Значит гипотеза о том, что алгоритм (3.14) дают случайные последовательности с гамма–распределением принимаются. При >0.003 некоторые реализации (3.14) приводили к неустойчивости. Для устранения этого явления, по всей видимости, необходимо использовать приближения более высокого порядка. Динамическая система, описываемая разностным уравнением (3.14), является нелинейной. Поэтому провести более подробный анализ устойчивости этих алгоритмов не представляется возможным.