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

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

Синтаксис

W = wishrnd(SIGMA,df)
W = wishrnd(SIGMA,df,D)
[W,D] = wishrnd(SIGMA,df)

Описание

W = wishrnd(SIGMA,df) служит для генерации матицы случайных чисел распределенных по распределению Уишарта для ковариационной матрицы SIGMA с числом степеней свободы df.

W = wishrnd(SIGMA,df,D) в отличии от предыдущего варианта синтаксиса задается коэффициент Холецкого D ковариационной матрицы SIGMA. В случае нескольких вызовов функции wishrnd с одинаковыми значениями ковариационной матрицы для увеличения эффективности алгоритма генерации целесообразно явно задавать коэффициент Холецкого D.

[W,D] = wishrnd(SIGMA,df) в отличии от первого варианта вызова функции кроме матрицы псевдослучайных чисел W возвращается значение коэффициента Холецкого DI, который может быть использован при последующих вызовах функции wishrnd.

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

Генерация псевдослучайных чисел по распределению Уишарта.

>> SIGMA=pascal(3)
SIGMA =
     1     1     1
     1     2     3
     1     3     6
>> df=5
df =
     5
>> W=wishrnd(SIGMA,df)
W =
    4.3741    1.7676   -1.2564
    1.7676    2.1324    3.4334
   -1.2564    3.4334   16.1057
   
Генерация псевдослучайных чисел по распределению Уишарта и расчет коэффициента Холецкого DI.

>> SIGMA=pascal(3)
SIGMA =
     1     1     1
     1     2     3
     1     3     6
>> df=5
df =
     5
>> [W,DI]=wishrnd(SIGMA,df)
W =
    3.4728    2.8301    2.1133
    2.8301    5.5233    8.7450
    2.1133    8.7450   22.9138
DI =
     1     1     1
     0     1     2
     0     0     1
         
Расчет псевдослучайных чисел согласно распределению Уишарта и использование 
коэффициента Холецкого DI при повторной генерации.

>> SIGMA=pascal(5)
SIGMA =
     1     1     1     1     1
     1     2     3     4     5
     1     3     6    10    15
     1     4    10    20    35
     1     5    15    35    70
>> df=10
df =
    10
>> [W,DI]=wishrnd(SIGMA,df)
W =
    10.7636    8.9244    4.7621    -6.4387        -32.4934
    8.9244     13.6604   14.7491   10.0681      -4.8925
    4.7621     14.7491   24.6901   34.7112      42.2157
   -6.4387     10.0681   34.7112   74.8364     135.2491
   -32.4934   4.8925     42.2157   135.2491   304.2332
DI =
     1     1     1     1     1
     0     1     2     3     4
     0     0     1     3     6
     0     0     0     1     4
     0     0     0     0     1
>> W=wishrnd(SIGMA,df,DI)
W =
   10.1921   15.0520   15.6658    9.8037         -2.7973
   15.0520   28.1868   36.0149    32.3553       18.0655
   15.6658   36.0149   56.1620    66.3571       67.0877
    9.8037    32.3553   66.3571    106.9305     157.7394
   -2.7973   18.0655    67.0877    157.7394     311.6346
 

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