|
Текст подпрограммы и версий ( Фортран ) gsb1r.zip |
Тексты тестовых примеров ( Фортран ) tgsb1r.zip |
|
Текст подпрограммы и версий ( Си ) gsb1r_c.zip |
Тексты тестовых примеров ( Си ) tgsb1r_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) gsb1r_p.zip |
Тексты тестовых примеров ( Паскаль ) tgsb1r_p.zip |
Генерация массива псевдослучайных чисел, имеющих бета-распределение.
Генерируются псевдослучайные числа, с плотностью распределения
f(x) = x(P-1)*(1-x)(Q-1)*(Г(P+Q)/Г(P)*(Q)) , 0≤x≤1 , P>0, Q>0.
Здесь:
P и Q - параметры бета-распределения,
Г - гамма функция.
| 1. |
Johnk,M.D., "Erzung von betaverteiler und gammaverteiler zufallzahlen", Metrica, 8(2), 1964. |
| 2. |
Phillips, Don T., and Beightler, Charles S., "Procedure for generating gamma variates with non-integer parameter sets", Journal of Statistical Computation and Simulation, 1, 1972, 203. |
SUBROUTINE GSB1R (ISEED, P, Q, N, R)
Параметры
| ISEED - | целая переменная, значение которой перед обращением к подпрограмме может быть любым целым числом в пределах [1,2147483646]; по окончании работы ей присваивается новое значение, котоpое может быть использовано при последующем вхождении в подпрограмму; |
| P - | заданный параметр бета - распределения (тип: вещественный); |
| Q - | заданный параметр бета - распределения (тип: вещественный); |
| N - | заданное количество генерируемых псевдослучайных чисел (тип: целый); |
| R - | вещественный массив длины N, содержащий вычисленные псевдослучайные числа. |
Версии: нет
Вызываемые подпрограммы
| GSU1R - | генерация массива псевдослучайных чисел, pавномеpно распределенных в интервале (0,1); |
Замечания по использованию
| Время работы подпрограммы pастет с pостом P и Q. |
DIMENSION R(3)
ISEED = 1432518450
P = 0.5
Q = 0.8
N = 3
CALL GSB1R (ISEED, P, Q, N, R)
Результаты:
R(1) = 0.364229459422
R(2) = 0.959182891531
R(3) = 0.132878061951