|
Раздел "Математика\Statistics Toolbox"
Список функций Statistics Toolbox
\ \ Функции планирования эксперимента
|
|
Функция генерации псевдослучайных чисел по многомерному нормальному распределению |
Синтаксис
R = mvnrnd(MU,SIGMA)
R = mvnrnd(MU,SIGMA,cases)
Описание R = mvnrnd(MU,SIGMA) генерирует матрицу n-d псевдослучайных чисел распределенных по многомерному нормальному закону с параметрами математического ожидания MU и ковариации SIGMA. Размерность матрицы MU n-d. Функция mvnrnd генерирует каждый ряд R на основе соответствующего ряда значений MU. Матрица SIGMA должна быть квадратной и положительно определенной. SIGMA может быть задана матрицей с размерностью d-d или 3-х мерным массивом с размерностью d-d-n. Если SIGMA задана как 3-х мерный массив, то каждый ряд R генерируется с использованием страницы массива SIGMA, т.е. R(i,:) генерируется с использованием MU(i,:) и SIGMA(:,:,i). Если MU задан как вектор с размерностью 1-d, то функция mvnrnd формирует матрицу размерности соответствующей размерности SIGMA.
r = mvnrnd(MU,SIGMA,cases) генерирует матрицу псевдослучайных чисел с размерностью cases-d распределенных по многомерному нормальному распределению для вектора средних MU с размерностью 1-d и матрицы SIGMA с размерностью d-d.
Примеры использования функции генерации псевдослучайных чисел
Генерация псевдослучайных чисел по многомерному нормальному распределению.
>> MU=[0 1 2];
>> SIGMA=pascal(3)
SIGMA =
1 1 1
1 2 3
1 3 6
>> R = mvnrnd(MU,SIGMA)
R =
-0.8436 0.6541 3.6404
Генерация псевдослучайных чисел по многомерному нормальному распределению при условии, что SIGMA задана как 3-х мерный массив.
>> MU=[0 0 0; 1 1 1; 2 2 2];
>> SIGMA = cat(3, pascal(3), pascal(3) , pascal(3))
SIGMA(:,:,1) =
1 1 1
1 2 3
1 3 6
SIGMA(:,:,2) =
1 1 1
1 2 3
1 3 6
SIGMA(:,:,3) =
1 1 1
1 2 3
1 3 6
>> R = mvnrnd(MU,SIGMA)
R =
-0.3304 -0.8303 -1.3661
0.8252 -0.1320 0.2033
2.4409 3.7219 4.5051
Генерация двумерного нормального распределения и его графическое представление.
>> mu = [2 3];
>> sigma = [1 1.5; 1.5 3];
>> r = mvnrnd(mu,sigma,100);
>> plot(r(:,1),r(:,2),'+')
>> grid on
\ \
|
