|
|
|
11 |
|
|
Р( Ди / k1 |
k2 ) |
|
0,7 0,005 |
0,01346 . |
|
|
|
|
|||
|
0,005 0,3 0,855 |
||||
|
0,7 |
|
|||
3 Диагностическая матрица в методе Байеса (метод Биргера)
Для определения вероятностей диагнозов по методу Байеса необходимо составить диагностическую матрицу (табл. 1), которая формируется на основе предварительного статистического материала. В этой таблице содержатся вероятности разрядов признаков при различных диагнозах.
Если признаки двухразрядные (простые признаки “да-нет”), то в таблице достаточно указать вероятность появления признака P(kj/Di). Bepоятность отсутствия признака P(k j / Di ) 1 P(k j / Di ) .
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 1 |
||
|
|
|
|
Пример диагностической матрицы в методе Байеса |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Признак kj |
|
|
|
|
|
||
|
|
|
|
|
|
|
k1 |
|
|
|
|
k2 |
|
|
k3 |
|
|
|
|
|
|
Диагноз |
|
|
|
|
|
|
|
|
|
|
|
|
|
P(Di ) |
|
|
|
|
Di |
) |
|
) |
|
) |
) |
|
) |
|
) |
) |
) |
) |
|
|
|
|
|
|
i |
|
i |
|
i |
i |
|
i |
|
i |
i |
i |
i |
|
|
|
|
|
|
/D |
|
/D |
|
/D |
/D |
|
/D |
|
/D |
/D |
/D |
/D |
|
|
|
|
|
|
11 |
|
12 |
|
13 |
21 |
|
22 |
|
23 |
24 |
31 |
32 |
|
|
|
|
|
|
P(k |
|
P(k |
|
P(k |
P(k |
|
P(k |
|
P(k |
P(k |
P(k |
P(k |
|
|
|
|
|
D1 |
0,8 |
|
0,2 |
|
0 |
0,1 |
|
0,1 |
|
0,6 |
0,2 |
0,2 |
0,8 |
0,3 |
|
|
|
|
D2 |
0,1 |
|
0,7 |
|
0,2 |
0 |
|
0 |
|
0,3 |
0,7 |
0,1 |
0,9 |
0,1 |
|
|
|
|
… |
… |
… |
… |
… |
… |
… |
… |
… |
… |
… |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
Однако более удобно использовать единообразную форму полагая, |
|||||||||||||||
например, |
для |
|
двухразрядного |
|
признака |
P(k j / Di ) P(k j1 / Di ) , |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
P(k |
j / Di ) P(k j 2 / Di ) . |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
m j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отметим, что P(k js / Di ) 1, где mj |
– число разрядов признака kj. Сумма |
||||||||||||||
s 1
вероятностей всех возможных реализаций признака равна единице.
В диагностическую матрицу включены априорные вероятности диагнозов. Процесс диагностирования в методе Байеса состоит в формировании диагностической матрицы и ее анализе. Важно предусмотреть возможность
12
уточнения таблицы в процессе диагностики. Для этого в памяти ЭВМ следует хранить не только значения P(kjs /Di ), но и следующие величины: N – общее число объектов, использованных для составления диагностической матрицы; Ni
– число объектов с диагнозом Di; Nij – число объектов с диагнозом Di, обследованных по признаку kj. Если поступает новый объект с диагнозом Dμ, то проводится корректировка прежних априорных вероятностей диагнозов следующим образом:
|
Ni |
|
|
N |
|
|
|
|
|
||||
|
|
|
P(Di ) |
|
|
; i 1, 2, ..., n; i ; |
|
||||||
|
N 1 |
N 1 |
|
||||||||||
P(Di |
) |
N |
1 |
|
|
N |
1 |
|
(20) |
||||
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
P(D ) N 1 N 1; i . |
|
||||||||
|
N |
|
|||||||||||
Далее вводятся поправки к вероятностям признаков. Пусть у диагностируемого объекта с диагнозом Dμ выявлен разряд r признака kj. Тогда для дальнейшей диагностики принимаются новые значения вероятности интервалов признака kj при диагнозе Dμ:
|
|
|
N |
|
|
r; |
|
|||||
P(k js |
/ D ) |
|
|
j |
|
; s |
|
|||||
N |
|
1 |
|
|||||||||
|
|
|
j |
|
|
|
|
|
|
|||
P(k js / D ) |
|
|
|
|
|
|
|
|
|
|
(21) |
|
|
|
N j |
|
|
|
1 |
|
|
||||
P(k |
|
/ D ) |
|
|
|
|
; s r. |
|||||
js |
|
|
|
|
|
|
1 |
|||||
|
|
N |
j |
1 N |
j |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
||
Условные вероятности признаков при других диагнозах корректировки не требуют.
4 Решающее правило Решающее правило – правило, в соответствии с которым принимается
решение о диагнозе. В методе Байеса объект с комплексом признаков K* относится к диагнозу с наибольшей (апостериорной) вероятностью
K * D , если P(D / K * ) P(D |
/ K * ) ( j 1,2, ..., n; i j) . |
(22) |
||
i |
i |
j |
|
|
Символ Є, применяемый в функциональном анализе, означает принадлежность множеству. Условие (22) указывает, что объект, обладающий данной реализацией комплекса признаков K* или, короче, реализация K*
13
принадлежит диагнозу (состоянию) Di. Правило (22) обычно уточняется введением порогового значения для вероятности диагноза:
P(D / K * ) P , |
(23) |
|
i |
i |
|
где Рi – заранее выбранный уровень распознавания для диагноза Di. При этом вероятность ближайшего конкурирующего диагноза не выше 1 Рi. Обычно принимается Pi 0,9 . При условии
P(D / K * ) P |
(22) |
|
i |
i |
|
решение о диагнозе не принимается (отказ от распознавания).
Метод Байеса имеет недостатки, заключающиеся в том, что он имеет повышенные погрешности при распознавании редких диагнозов. Для их уменьшения в практических расчетах целесообразно провести диагностику, положив диагнозы равновероятными
P(D ) 1 |
n |
. |
(24) |
i |
|
|
|
|
|
|
Тогда наибольшим значением апостериорной вероятности будет обладать диагноз Di , для которого P(K*/Di) максимальна:
K * D , если P(K * / D ) P(K * / D |
) (j = 1,2 ,…, n; i≠j). |
(25) |
||
i |
i |
j |
|
|
Иными словами, устанавливается диагноз Di, если данная совокупность признаков чаще встречается при диагнозе Di, чем при других диагнозах. Такое решающее правило соответствует методу максимального правдоподобия. Из предыдущего вытекает, что этот метод является частным случаем метода Байеса при одинаковых априорных вероятностях диагнозов. В методе максимального правдоподобия “частные” и “редкие” диагнозы равноправны. Для надежности распознавания условие (24) должно быть дополнено пороговым значением
P(K * / D ) P , |
(26) |
i i |
|
где Pi – заранее выбранный уровень распознавания для диагноза Di.
14
Задача 3. Метод Биргера в определении состояния объекта
Пусть при наблюдении за газотурбинным двигателем проверяются два признака: k1 повышение температуры газа за турбиной более чем на 50 0С и k2увеличение времени выхода на максимальную частоту вращения более чем на 5 с. Предположим, что для данного типа двигателей появление этих признаков связано либо с неисправностью топливного регулятора (состояние D1), либо с увеличением радиального зазора в турбине (состояние D2).
При нормальном состоянии двигателя (состояние D3) признак k1 не наблюдается, a признак k2 наблюдается в 5 % случаев. На основании статистических данных известно, что 80 % двигателей вырабатывают ресурс в нормальном состоянии, 5 % двигателей имеют состояние D1 и 15 % – состояние D2. Известно также, что признак k1 встречается при состоянии D1 в 20 %, а при состоянии D2 в 40 % случаев; признак k2 при состоянии D1 встречается в 30 %, а при состоянии D2 – в 50 % случаев.
Решение. Сведем эти данные в диагностическую таблицу 2. Найдем сначала вероятности состояний двигателя, когда обнаружены оба признака k1 и k2. Для этого, считая признаки независимыми, применим формулу (13).
Вероятность диагноза D1 при условии, что проявились оба признака k1 и k2
P(D1 |
/ k1 k2 ) |
|
|
0,05 0,2 0,3 |
|
|
0,09 . |
||
|
|
|
|
|
|||||
0,05 0,2 0,3 0,15 0,4 0,5 |
|
|
|||||||
|
|
|
0,8 0 0,05 |
||||||
Аналогично получим |
|
|
|
|
|
|
|||
|
|
P(D2 / k1 |
k2 ) 0,91; P(D3 / k1 |
k2 ) 0. |
|||||
|
|
|
|
|
|
|
|
Таблица 2 |
|
Вероятности признаков и априорные вероятности состояний |
|||||||||
|
|
|
|
|
|
|
|
||
Di |
|
P(k1/Di) |
|
P(k2/Di) |
|
P(Di) |
|
||
|
|
|
|
|
|
|
|
|
|
D1 |
|
0,2 |
|
|
0,3 |
|
0,05 |
|
|
|
|
|
|
|
|
|
|
|
|
D2 |
|
0,4 |
|
|
0,5 |
|
0,15 |
|
|
|
|
|
|
|
|
|
|
|
|
D3 |
|
0,0 |
|
|
0,05 |
|
0,80 |
|
|
|
|
|
|
|
|||||
Определим вероятность состояний двигателя, если обследование |
|||||||||
показало, что |
повышение температуры не наблюдается (признак k1 |
||||||||
15
отсутствует), но увеличивается время выхода на максимальную частоту вращения (признак k2 наблюдается). Отсутствие признака k1 есть признак
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
наличия противоположного события k1 : |
P(k1 / Di ) 1 P(k1 / Di ) . |
|||||||||||||||
Для расчета применяют также формулу (13), но значение P(k1/Di) в |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
диагностической таблице заменяют на P(k1 / Di ) . В этом случае |
||||||||||||||||
|
|
|
|
|
|
0,05 0,8 0,3 |
|
|
||||||||
P(D1 / k1 k2 ) |
|
|
|
|
0,12 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
0,05 0,8 0,3 |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
0,15 0,6 0,5 0,8 1 0,05 |
|||||||||||||
|
|
|
|
|
|
|
|
|
||||||||
и аналогично P(D2 / k1 |
k2 ) 0,46; |
|
P(D3 / k1 |
k2 ) 0,41. Вычислим |
||||||||||||
вероятности состояний в том случае, когда оба признака отсутствуют. Аналогично предыдущему получим
|
|
|
|
|
|
|
|
|
|
|
0,05 0,8 0,7 |
|
|
|
|
|||
P(D1 / k1 |
k 2) |
|
|
|
|
|
|
0,03; |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
0,05 |
0,8 0,7 0,15 0,6 0,5 0,8 1 0,95 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
P(D2 / k1 |
k2 ) 0,05; P(D3 |
/ k1 |
k2 ) 0,92. |
||||||||||||
Отметим, что вероятности состояний D1 |
и D2 |
отличны от нуля, так как |
||||||||||||||||
рассматриваемые признаки не являются для них детерминирующими. Из проведенных расчетов можно установить, что при наличии признаков k1 и k2 в двигателе с вероятностью 0,91 имеется состояние D2, т. е. увеличение радиального зазора. При отсутствии обоих признаков наиболее вероятно нормальное состояние (вероятность 0,92). При отсутствии признака k1 и наличии признака k2 вероятности состояний D2 и D3, примерно одинаковы (0,46 и 0,41) и для уточнения состояния двигателя требуется проведение дополнительных обследований.
Сводное решение данной задачи в виде электронной таблицы Excel представлено на рис. 3.
5 Метод последовательного анализа
При использовании метода Байеса для распознавания состояний D1 и D2 следует составить отношение (для независимых признаков):
P(D2 / K * ) .
P(D1 / K * )