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

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

Синтаксис

[M,V] = weibstat(A,B)

Описание

[M,V] = weibstat(A,B) функция служит для расчета математического ожидания и дисперсии распределения Вейбулла с заданными параметрами А и В. Размерность векторов и матриц А и В должна быть одинаковой. Скалярный входной параметр увеличивается до матрицы постоянных значений с размерностью второго параметра. Размерность матриц математического ожидания M и дисперсии V равна размерности входных параметров.

Математическое ожидание M распределения Вейбулла с заданными параметрами А и В определяется по формуле

Дисперсия V распределения Вейбулла с заданными параметрами А и В определяется по формуле

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

Расчет математического ожидания и дисперсии распределения Вейбулла для пары параметров А и В.

>> A=1
A =
     1
>> B=2
B =
     2
>> [M,V] = weibstat(A,B)
M =
    0.8862
V =
    0.2146
        
Расчет математического ожидания и дисперсии распределения Вейбулла для матриц параметров А и В.

>> A=[1 2 3; 4 5 6]
A =
     1     2     3
     4     5     6
>> B=[2 3 4; 5 6 7]
B =
     2     3     4
     5     6     7
>> [M,V] = weibstat(A,B)
M =
    0.8862    0.7088    0.6887
    0.6958    0.7094    0.7242
V =
    0.2146    0.0664    0.0373
    0.0254    0.0189    0.0148
        
Расположение математического ожидания M и значений (M-2); (M-); (M+); (M+2); (M+3) на графике 
функции плотности вероятности,  - среднее квадратическое отклонение.

>> A=1;
>> B=2;
>> X=0:0.01:3;
>> Y= weibpdf(X,A,B);
>> plot(X,Y)
>> grid on
>>  [M,V] = weibstat(A,B)
M =
    0.8862
V =
    0.2146
>> sigma=sqrt(V)
sigma =
    0.4633
>> H=line ([M M], [0 weibpdf(M,A,B)+0.05]);
>> set(H,'LineWidth',3, 'Color','k')
>> H1=line ([M-sigma M-sigma], [0 weibpdf(M-sigma,A,B)+0.05]);
>> set(H1,'Color','g')
>> H2=line ([M+sigma M+sigma], [0 weibpdf(M+sigma,A,B)+0.05]);
>> set(H2,'Color','g')
>> H3=line ([M-2*sigma M-2*sigma], [0 weibpdf(M-2*sigma,A,B)+0.05]);
>> set(H3,'Color','m')
>> H4=line ([M+2*sigma M+2*sigma], [0 weibpdf(M+2*sigma,A,B)+0.05]);
>> set(H4,'Color','m')
>> H5=line ([M+3*sigma M+3*sigma], [0 weibpdf(M+3*sigma,A,B)+0.05]);
>> set(H5,'Color','c')

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