ГОСТ Р ИСО 28640- 2012
6.7.2.2 Алгоритм для с = к (к - целое число)
Используя независимые равномерные случайные числа Ut, U>, ..., {/iT применяют
формулу
Y = a -b \n \(\-U l)(\-U2)... (1-Ut)}.
П рим ечание - Этим методом для а = 0 н b = 2 может быть получено распределение у’ с четным
числом степеней свободы.
6.7.2.3 Алгоритм для с = к+ 1/2 (А’- целое число)
Используя стандартное нормальное случайное число Z(Jи равномерные случайные числа
U |. U
2
, ..., Uк, применяют формулу
J’=a +ft[z‘/
2
-ln{(|-£/IXl-t/!).-(l-t/.)!].
В случае, когда к = 0, член с логарифмом исчезает.
П рим ечание
-
Тем же методом при
а=0
и
Ь= 2
может быть получено распределение
у~с
нечетным
числом степеней свободы.
6.7.2.4 Приближенный метод для с > 1/3
а) Задают г =с - 1/3. s =y[r . t =г - Нп(г), р - 1/(3-^) и q - -3 ■Jr.
b) Геиернруют стандартное нормальное случайное число Z.
c) Если Z < с/, то переходят к Ь).
d) Вычисляют Y= {pZ. +,v)\ V= Z‘/2 и генерируют U.
e) Если (Y -r)‘/Y - V<U, выполняют Y:=о +bY(конец).
0Вычисляют W= Y - rln(У )-/ - V.
g) Если W <, U. то выполняют Y :=a +bY (конец).
h) Если IV> —In(1.0 - U). то переходят к Ь).
П рим ечание - Метод основан на преобразовании Унлсона-Хилферти, приводящем / -распределение к
приближенному стандартному нормальному распределению. Точность такого приближения зависит от значения
параметра с. Идея преобразования состоит в следующем: абсолютная разность между процентной точкой
приближенного и точного распределений всегда меньше 0.2.
6.7.2.5 Точный метод генерации Ченга для с > 1/2
а)Задают q = с - In 4 и г = с + *Jlc- 1.
2
b)Генерируют стандартные равномерные случайные числа Ui и U .
’
с)
Вычисляют V =dn
\-и,
, W= г ехр((/,), Z = и !и , ,R =q +rV - IV.
d) Если R >4.5Z - (1 + 1п4.5), то вычисляют Y=а+ b\V (выход).
e) Если R >InZ. то вычисляют Y- а + MV(выход).
0Геиернруют стандартные равномерные случайные числа (/, и (Л и вычисляют
II