Список функций Statistics Toolbox

В оглавление \ К следующему разделу \ К предыдущему разделу

Синтаксис

R = gamrnd(A,B)
R = gamrnd(A,B,m)
R = gamrnd(A,B,m,n)

Описание

R = gamrnd(A,B) функция предназначена для генерации псевдослучайного числа по Гамма распределению для каждой пары параметров A и B. Размерность векторов или матриц параметров A и B должна быть одинаковой. Скалярный параметр увеличивается до размера остальных входных аргументов.

R = gamrnd(A,B,m) позволяет получить вектор псевдослучайных чисел на m элементов по Гамма распределению для параметров A и B.

R = gamrnd(A,B,m,n) ) позволяет получить матрицу псевдослучайных чисел с размерностью m-n элементов по Гамма распределению для параметров A и B.

Примеры использования функции генерации псевдослучайных чисел

Генерация одного числа соответствующего заданной паре значений параметров A и B.

>> A=1;
>> B=2;
>> R = gamrnd(A,B)
R =
    0.1023
>> A=[1 2 3 4 5];
>> B=5;
>> R = gamrnd(A,B)
R =
    7.3237    7.5060    9.5972   36.0827   31.6181

Генерация вектора псевдослучайных чисел с размерностью 14.

>> A=5;
>> B=7;
>> m=[1 4];
>> R = gamrnd(A,B,m)
R =
   34.2698   47.0043   23.1796   42.3183

Второй вариант генерации вектора с размерностью 14.

>> A=1;
>> B=3;
>> m=1; n=5;
>> R = gamrnd(A,B,m,n)
R =
    8.5476    3.1250    0.6205   13.8574    5.9222

Генерация матрицы псевдослучайных чисел с размерностью 44.

>> A=2;
>> B=3;
>> m=[4 4];
>> R = gamrnd(A,B,m)
R =
   10.1321    2.9526   14.7757    3.3462
    9.8778   12.1505    2.5314    8.7751
    4.4587    6.9010    1.8665    5.5206
    4.3094    2.2406    9.5973    9.3098

Другой вариант генерации матрицы с размерностью 44.

>> A=1;
>> B=3;
>> m=4; n=4;
>> R = gamrnd(A,B,m,n)
R =
    1.8816    1.6965    2.4299    0.1323
    0.9560    2.9794    1.0934    1.9469
    3.5204    1.0583    1.4278    0.3830
    0.5284    1.8123    0.6889    5.2642

Графическая оценка качества генератора псевдослучайных чисел

>> A=1;
>> B=3;
>> N=9;
>> R = gamrnd(A,B,[1 100]);
>> hist(R, N)
>> grid on
>> X=min(R):(max(R)-min(R))/100:max(R);
>> f= gampdf(X, A,B);
>> ff=f*100*((max(R)-min(R))/N);
>> hold on
>> plot(X,ff,'r')
>> hold off

В оглавление \ К следующему разделу \ К предыдущему разделу