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

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

Синтаксис

[M,V] = ncx2stat(NU,DELTA)

Описание

[M,V] = ncx2stat(NU,DELTA) функция служит для расчета математического ожидания и дисперсии смещенного хи-квадрат распределения с заданными параметрами числа степеней свободы NU и параметра смещения DELTA. Размерность векторов и матриц NU, DELTA должна быть одинаковой. Скалярный входной параметр увеличивается до матрицы постоянных значений с размерностью второго параметра. Размерность матриц математического ожидания M и дисперсии V равна размерности входных параметров.

Математическое ожидание M смещенного хи-квадрат распределения с заданными параметрами NU и DELTA определяется по формуле

M=DELTA+NU

Дисперсия V смещенного хи-квадрат распределения с заданными параметрами NU и DELTA определяется по формуле

V=2*(NU+2*DELTA)²

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

Расчет математического ожидания и дисперсии смещенного хи-квадрат распределения для сочетания значений
параметров NU, DELTA.

>> NU=5
NU =
     5
>> DELTA =2
DELTA =
     2
>>[M,V] = ncx2stat (NU,DELTA)
M =
    7
V =
    18
        
Расчет математического ожидания и дисперсии смещенного хи-квадрат распределения для матриц параметров NU, DELTA.

>> NU=[10 20; 30 40]
NU =
    10    20
    30    40
>> DELTA=[1 2; 3 4]
DELTA =
     3     5
     7     9
>>[M,V] = ncx2stat (NU,DELTA)
M =
    11    22
    33    44
V =
    24    48
    72    96
        
Расположение математического ожидания M и значений (M-3); (M-2); (M-); (M+); (M+2); (M+3) на графике функции 
плотности вероятности,  - среднее квадратическое отклонение.

>> DELTA=5;
>> NU=10;
>> X=1:0.01:15;
>> Y= ncx2pdf(X,NU,DELTA);
>> plot(X,Y)
>> grid on
>>[M,V] = ncx2stat(NU,DELTA)
M =
   15
V =
   40
>> sigma=sqrt(V)
sigma =
   6.3246
>> H1=line ([M-3*sigma M-3*sigma], [0 ncx2pdf((M-3*sigma),NU,DELTA)+0.02]);
>> set(H1,'Color','m')
>> H2=line ([M-2*sigma M-2*sigma], [0 ncx2pdf((M-2*sigma),NU,DELTA)+ 0.02]);
>> set(H2,'Color','g')
>> H3=line ([M-sigma M-sigma], [0 ncx2pdf((M-sigma),NU,DELTA)+0.02]);
>> set(H3,'Color','c')
>> H4=line ([M M], [0 ncx2pdf((M),NU,DELTA)+0.02]);
>> set(H4,'LineWidth',3, 'Color','k')
>> H5=line ([M+sigma M+sigma], [0 ncx2pdf((M+sigma),NU,DELTA)+0.02]);
>> set(H5,'Color','c')
>> H6=line ([M+2*sigma M+2*sigma], [0 ncx2pdf((M+2*sigma),NU,DELTA)+0.02]);
>> set(H6,'Color','g')
>> H7=line ([M+3*sigma M+3*sigma], [0 ncx2pdf((M+3*sigma),NU,DELTA)+0.02]);
>> set(H7,'Color','m')

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