Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Антонов. Системный анализ. Учебник для вузов.doc
Скачиваний:
449
Добавлен:
19.02.2016
Размер:
3.82 Mб
Скачать

^2/(0.0

Э/(0,О

<Я(0,

Э03

Эв

, 9 - вектор параметров функции F( 9, t), і-I, п.

где а, = -M

Метод максимального правдоподобия широко используется при оце­нивании параметров сложных систем. Этот метод служит основой про­цедур проверки статистических гипотез и доверительного интерваль­ного оценивания. Оценки характеристик, получаемые методом макси­мального правдоподобия, обладают рядом важных свойств, таких как несмещенность, асимптотическая эффективность, состоятельность [35].

Рассмотрим применение метода максимального правдоподобия для решения задач оценивания параметров законов распределения. Вид закона распределения исследуемой случайной величины будем предпо­лагать известным.

Определим функцию правдоподобия параметра 9 как неотрицатель­ную вещественную функцию L(9, t), заданную на множестве QxT, про­порциональную функции плотности распределения:

WO = Ylf(BJi), і = й, (7.1)

где 0 - область определения вектора параметров Q; T- область опре­деления наблюдаемой случайной величины t, по результатам наблюде­ния .за которой про-изводится оценивание параметров 9; Tj - реализа­ция случайной величины V, 9 - в общем случае вектор параметров за­кона распределения.

Оценкой максимального правдоподобия для заданной функции прав­доподобия L(9, t) является функция Q(T), удовлетворяющая соотноше­нию [36]:

£{0(7-),7} = sup L(Q,T).

Єє©

Для нахождения оценки 9 решают уравнение

= 0. (7.2)

эе

Поскольку In L(9, t) при фиксированных TvT2, ...,Tn достигает макси­мума при том же значении 9, что и L(9, t), значения 9 можно опреде­лять, решая уравнение

Э1пЦЭ,7-)_0

Э0

Для величины In L(9, 0 используют обозначение /(9, /)•

7.2. Использование метода максимального правдоподобия для оценивания параметров законов распределения

где Z71(Z), F2(t), H(f) - некоторые функции, удовлетворяющие следую­щим условиям:

  • F1(I), F2(I) интегрируемы на (-«>, °°) и интеграл J H(t)f(Q,t)dt < P,

причем р не зависит от 9;

  • для каждого 9 из 0 интеграл

конечен и положителен.

Используя метод максимального правдоподобия при наличии боль­ших объемов выборки, можно произвести оценивание дисперсий век­тора оценок 9. Для этого составляется информационная матрица Фи­шера

Одним из важных достоинств метода максимального правдоподо­бия является то, что он позволяет получить асимптотически нормаль­ные и эффективные оценки параметров функции распределения случай­ной величины t. Для этого необходимо, чтобы выполнялся ряд условий, называемых условиями регулярности [32].

Для каждого 9, принадлежащего некоторому невырожденному ин­тервалу 0, существуют производные

Э In /(0,0. Э2 In/(0,0. Э3 In/(0,о.

Э0 : Э02 ’ Э03

faU

а\г

aIn

Il

а2\

а22

.. а2п

.«»1

“„2 '

.. ат

для каждого 9 из 0 имеем

Э2 InL

159,30, J

Э3 In /(0,0

<^0);

Э02

Дисперсии и ковариации оценок 9.и 9 определяются из ковариаци­онной матрицы V, которая является обратной матрице I:

ние

/=1

:0.

ЭХ.

А.

A2

Cf

А.

A2

•• D2n

А,

A2

Dm

V= I1=



D(Qi) при i = j;

где Dy =

Cov(QitQj) при ІФ j.

Приведем формулы для определения дисперсии параметров зако­нов распределения для двухпараметрических плотностей. Пусть а и (3 - оцениваемые параметры плотности / (а, (3, t). Пусть определены элементы информационной матрицы

  1. =

Экспоненциальное распределение

Рассмотрим имеющее важное прикладное значение экспоненциаль­ное распределение. Например, данное распределение широко исполь­зуется в теории надежности для описания случайной величины наработки до отказа. Плотность экспоненциального распределения имеет вид f3(X, t); функция распределения F (А,, 0 = 1- exp(-Ai). Параметр X назы­вается интенсивностью отказов. Запишем функцию правдоподобия

п Ґ п ^

L3 (А, 0 = ]^[ Аехр(-XTj) = X" ехр -X^T1 .

Li Iv .-=1 ,

Для определения оценки параметра X необходимо решить уравне-

nln А-А^7]

dl(X,t)

ЭХ

Э/(а,р,0

(7.3)

(7-4)

(7.5)

(7-6)

где

Э^(аДО__Э2/(ос,р,0

> a2l aU і <Я22 =

Эа212 ЭаЭр Дискриминант матрицы равен

Dis = аиа22 -Ctf1.

Тогда дисперсию параметров а и P определим по формулам

эр2

D = _f3L. Dis ’

D =--

Dis

Ковариация будет вычисляться следующим образом:

cov(a,P) = —^-.

Dis

Перейдем к рассмотрению примеров применения метода макси­мального правдоподобия для оценивания параметров некоторых зако­нов распределения, имеющих важное значение в задачах системного анализа.

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

1Lt,

Дисперсия оценки параметра X характеризует точность этого парамет­ра и равняется

>2

ДА.] = —■

п

Для интенсивности отказов, зная оценку параметра и ее дисперсию, можно определить доверительный интервал с заданной доверительной вероятностью. Если обозначить верхнюю и нижнюю оценки интенсив­ности отказов через Afl и А,н, то можно определить

Ab = A+-pfp, Ah = X--j=t?, л/и

где - табулированная величина, которая зависит от уровня доверитель­ной вероятности P и определяется обратным интерполированием рас­пределения Стьюдента [37, табл. 6.2].

Определение среднего времени между реализациями событий про­изводится по формуле

интервальные оценки равны

Tu=J-, Tb = J-.

К

Вероятность безотказной работы определяется следующим обра­зом:

P3 (г Д) = ехр (-Xf), соответственно интервальные оценки вычисляются так:

pu(t,X) = e\ p(-V); Pb (їД) = ехр (-V)-

Рассмотрим далее пример оценивания параметров нормального закона распределения.

Нормальное распределение

ВВЕДЕНИЕ 5