|
Раздел "Математика\Statistics Toolbox"
Список функций Statistics Toolbox
\ \ Генерация псевдослучайных чисел по заданному закону распределения
|
Функция генерации псевдослучайных чисел распределенных по закону хи-квадрат () |
Синтаксис
R = chi2rnd(V)
R = chi2rnd(V,m)
R = chi2rnd(V,m,n)
Описание
R = chi2rnd(V) функция предназначена для генерации псевдослучайного числа распределенного по закону хи-квадрат для каждого значения параметра V.
R = chi2rnd(V,m) позволяет получить вектор псевдослучайных чисел на m элементов распределенных по закону хи-квадрат для параметра V, где m - вектор размерностью 12 определяющий размерность матрицы R.
R = chi2rnd(V,m,n) позволяет получить матрицу псевдослучайных чисел с размерностью m-n элементов распределенных по закону хи-квадрат для параметра V.
Примеры использования функции генерации псевдослучайных чисел
Генерация одного числа соответствующего определенному значению параметра V.
>> V=10;
>> R = chi2rnd(V)
R =
4.4268
>> V=[10 20 30 40 50];
>> R = chi2rnd(V)
R =
13.9829 23.2734 26.5743 29.3867 57.3204
Генерация вектора с размерностью 15 элементов псевдослучайных чисел.
>> V=10;
>> m=[1 5];
>> R = chi2rnd(V,m)
R =
10.2302 15.6096 9.4732 3.3293 9.2048
Другой вариант генерации вектора с размерностью 15.
>> V=10;
>> m=1; n=5;
>> R = chi2rnd(V,m,n)
R =
3.6552 11.6533 13.8338 9.8985 10.0190
Генерация матрицы с размерностью 55 элементов.
>> V=10;
>> m=[5 5];
>> R = chi2rnd(V,m)
R =
16.7829 9.0128 6.9741 14.9470 5.1372
11.1580 7.5976 6.8432 10.2689 11.8798
10.9589 6.5304 14.8066 11.7760 14.8122
8.5577 10.9080 9.5915 10.6171 9.9286
5.8218 9.3833 13.1405 5.1395 6.5166
Другой вариант генерации матрицы с размерностью 55.
>> V=10;
>> m=5; n=5;
>> R = chi2rnd(V,m,n)
R =
8.0671 5.7285 5.6716 6.3146 4.3965
15.8682 5.2127 11.7063 6.7905 8.6128
7.5513 29.6239 6.7613 5.0805 13.1623
4.6136 7.9643 6.3317 7.4396 9.4257
11.1463 12.4863 9.4614 6.4557 13.1582
Оценка качества генерации псевдослучайных чисел в виде наложения гистограммы и функции случайной величины, приведенной к абсолютному масштабу.
>> V=10;
>> N=9;
>> R = chi2rnd(V,[1 100]);
>> hist(R, N)
>> grid on
>> X=min(R):(max(R)-min(R))/100:max(R);
>> f= chi2pdf(X, V);
>> ff=f*100*((max(R)-min(R))/N);
>> hold on
>> plot(X,ff,'r')
>> hold off
\ \
|