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

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

Синтаксис

W=iwishrnd(SIGMA,df)
W=iwishrnd(SIGMA,df,DI)
[W,DI]=iwishrnd(SIGMA,df)

Описание

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

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

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

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

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

>> SIGMA=pascal(3)
SIGMA =
     1     1     1
     1     2     3
     1     3     6
>> df=5
df =
     5
>> W=iwishrnd(A,df)
W =
    0.6011    0.1055    0.0229
    0.1055    2.1282    0.8601
    0.0229    0.8601    0.3753

Генерация псевдослучайных чисел по обратному распределению Уишарта и расчет коэффициента Холецкого DI.

>> SIGMA=pascal(3)
SIGMA =
     1     1     1
     1     2     3
     1     3     6
>> df=5
df =
     5
>> [W,DI]=iwishrnd(SIGMA,df)
W =
    2.4815    2.2623    0.5948
    2.2623    3.6664    1.5391
    0.5948    1.5391    0.7814
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]=iwishrnd(SIGMA,df)
W =
    1.2555    2.9061    3.6567    1.9071    0.3948
    2.9061    7.5631    9.6383    5.1410    1.0714
    3.6567    9.6383   13.0949    7.3482    1.5956
    1.9071    5.1410    7.3482    4.3615    0.9920
    0.3948    1.0714    1.5956    0.9920    0.2357
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=iwishrnd(SIGMA,df,DI)
W =
    2.6936    1.6496   -0.7085      -1.4002    -0.4610
    1.6496    7.8608    10.2487    6.1963      1.4143
   -0.7085   10.2487   17.2284    11.7529    2.8645
   -1.4002    6.1963    11.7529    8.3968      2.0935
   -0.4610    1.4143    2.8645      2.0935      0.5297

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