gol_met_1-6
.pdf2.Задание.
3.Основные формулы и положения.
4.Расчет указанной вероятности (численный).
5.Выводы по работе.
ЗАДАНИЯ ПО РАБОТЕ
Из 1000 обследованных подшипников передней подвески автомобилей 900 подшипников выработали ресурс в исправном состоянии и 100 – в неисправном.
Все подшипники были обследованы по следующим признакам:
–общий уровень вибрации;
–температура;
–загрязнение смазки.
У70% исправных подшипников общий уровень вибрации лежал в диапазоне от 0,25 до 0,5 g, у 20% исправных подшипников – от 0,5 до 0,75 g и у 10% – >0,75g.
У80% исправных подшипников температура лежала в диапазоне 50–70 град, у 10% – в диапазоне 70–90 град. И у 10% – >90 град.
У90% исправных подшипников загрязнение смазки было в пределах нормы.
У80% неисправных подшипников наблюдалась вибрация >0,75 g, у
15% неисправных подшипников вибрация в диапазоне 0,5–0,75g.
У85% неисправных подшипников температура была >90 град, у 8% неисправных подшипников – в диапазоне 70–90 град.
У70% неисправных подшипников загрязнение смазки было выше нормы.
Рассчитать:
1.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,25–0,5g, температуры – 50–70 град, загрязнения смазки в пределах нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,25–0,5g, температура 50–70 град, загрязнение смазки в пределах нормы.
2.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,5–0,75g, температуры – 50–70 град, загрязнения смазки в пределах нормы.
11
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,5–0,75g, температура – 50–70 град, загрязнение смазки в пределах нормы.
3.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне >0,75g, температуры – 50–70 град, загрязнения смазки в пределах нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация >0,75g , температура – 50–70 град, загрязнение смазки в пределах нормы.
4.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,25–0,5g, температуры – 70–90 град, загрязнения смазки в пределах нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,25–0,5g, температура – 70–90 град, загрязнение смазки в пределах нормы.
5.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,5–0,75g, температуры – 70–90 град, загрязнения смазки в пределах нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,5–0,75g, температура – 70–90 град, загрязнение смазки в пределах нормы.
6.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне >0,75g, температуры – 70–90 град, загрязнения смазки в пределах нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация – >0,75g, температура – 70–90 град, загрязнение смазки в пределах нормы.
12
7.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,25–0,5g, температуры – >90 град, загрязнения смазки в пределах нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,25–0,5g, температура – >90 град, загрязнение смазки в пределах нормы.
8.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,5–0,75g, температуры – >90 град, загрязнения смазки в пределах нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,5–0,75g, температура – >90 град., загрязнение смазки в пределах нормы.
9.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне >0,75g, температуры – >90 град, загрязнения смазки
впределах нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация >0,75g, температура – >90 град., загрязнение смазки в пределах нормы.
10.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,25–0,5g, температуры – 50–70 град, загрязнения смазки в выше нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,25–0,5g, температура – 50–70 град, загрязнение смазки выше нормы.
11.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,5–0,75g, температуры – 50–70 град, загрязнения смазки в выше нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было
13
исправное состояние и наблюдались: вибрация 0,5–0,75g, температура – 50–70 град, загрязнение смазки выше нормы.
12.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне >0,75g, температуры – 50–70 град., загрязнения смазки в выше нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация >0,75g, температура – 50–70 град, загрязнение смазки выше нормы.
13.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,25–0,5g, температуры – 70–90 град, загрязнения смазки в выше нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,25–0,5g, температура – 70–90 град, загрязнение смазки выше нормы.
14.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,5–0,75g, температуры – 70–90 град., загрязнения смазки в выше нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,5–0,75g, температура – 70–90 град, загрязнение смазки выше нормы.
15.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне – >0,75g, температуры – 70–90 град, загрязнения смазки в выше нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация – >0,75g, температура – 70–90 град, загрязнение смазки выше нормы.
16.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,25–0,5g, температуры – 70–90 град, загрязнения смазки выше нормы.
14
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,25–0,5g, температура – >90 град, загрязнение смазки выше нормы.
17.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне 0,5–0,75g, температуры – 70–90 град, загрязнения смазки в выше нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,5–0,75g, температура – >90 град, загрязнение смазки выше нормы
18.Вероятность исправного состояния подшипника при наблюдении вибрации в диапазоне – >0,75g, температуры – >90 град, загрязнения смазки выше нормы.
Уточнить априорные вероятности появления исправного и неисправного состояний, а также условные вероятности признаков, если в результате обследования 1001 подшипника установлено, что у него было исправное состояние и наблюдались: вибрация 0,25–0,5g, температура – 50–70 град, загрязнение смазки выше нормы.
Практическая работа № 2
МЕТОД МИНИМАЛЬНОГО РИСКА
Цель работы: изучение метода минимального риска для диагностики технического состояния исследуемых систем и объектов.
МЕТОДИЧЕСКИЕ УКАЗАНИЯ
Вероятность принятия ошибочного решения слагается из вероятностей ложной тревоги и пропуска дефекта. Если приписать «цены» этим ошибкам, то получим выражение для среднего риска
∞ |
x0 |
|
R = C21P1 ∫ f (x / D1)dx + C12 P2 |
∫ f (x / D2 )dx. |
(1) |
x0 |
−∞ |
|
Разумеется, цена ошибки имеет условное значение, но она должна учесть предполагаемые последствия ложной тревоги и пропуска дефекта.
15
В задачах надежности стоимость пропуска дефекта обычно существенно больше стоимости ложной тревоги (C12 >> C21). Иногда вводится цена правильных решений С11 и С22 , которая для сравнения со стоимостью потерь (ошибок) принимается отрицательной. В общем случае средний риск (ожидаемая величина потери) выражается равенством
x0 |
∞ |
|
|
|
R = C11P1 ∫ f (x / D1 )dx + C21P1 |
∫ f (x / D1 )dx + |
|
||
−∞ |
x0 |
|
|
|
x0 |
∞ |
|
(2) |
|
+ C12 P2 ∫ f (x / D2 )dx + C22 P2 ∫ |
f (x / D2 )dx. |
|||
|
||||
−∞ |
x0 |
|
|
Величина x, предъявляемая для распознавания, является случайной и потому равенства (1) и (2) представляют собой среднее значение (математическое ожидание) риска.
f(x / D1) |
|
|
f(x / D1) |
|
|||
|
|
|
|
||||
|
|
|
|
|
|
f(x / D1) |
|
x0 max 0 |
|
|
x0 min |
|
|
x |
|
x1 |
|||||||
x |
|||||||
|
|
|
|
2 |
|
Рис. 1. Точки экстремума среднего риска ошибочных решений
Найдем граничное значение x0 из условия минимума среднего риска. Дифференцируя (2) по x0 и приравнивая производную нулю, получим сначала условие экстремума
|
dR |
=C P f (x |
/ D ) −C P f |
(x |
/ D ) +C |
P f (x |
/ D |
) −C P f (x |
/ D |
) =0 (3) |
||||||||
|
|
|||||||||||||||||
11 1 |
0 |
1 |
21 1 |
0 |
|
1 |
12 2 |
0 |
2 |
22 2 |
0 |
2 |
|
|||||
|
dx0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
f (x0 / D1) |
= |
(C12 −C22 )P2 |
. |
|
|
|
(4) |
||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
f (x |
/ D |
|
) |
|
(C |
−C |
)P |
|
|
|
|||
|
|
|
|
|
0 |
2 |
|
|
21 |
11 |
|
1 |
|
|
|
|
|
Это условие часто определяет два значения x0, из которых одно соответствует минимуму, второе – максимуму риска (рис. 1). Соотношение (4) является необходимым, но недостаточным условием миниму-
16
ма. Для существования минимума R в точке x = x0 вторая производная
должна быть положительной d 2R > 0 , что приводит к следующему ус- dx02
ловию относительно производных плотностей распределений:
f '(x0 / D1) |
< |
(C12 |
−C22 )P2 |
. |
(5) |
|
|
|
|
||||
f '(x0 / D2 ) (C21 −C11)P1 |
||||||
|
Если распределения f(x, D1) и f(x, D2) являются, как обычно, одномодальными (т. е. содержат не более одной точки максимума), то при
x1 < x0 < |
|
2 |
(6) |
x |
условие (5) выполняется. Действительно, в правой части равенства стоит положительная величина, а при x > x1 производная f '(x / D1) , тогда как
при x < x2 значение f '(x / D2 ) .
В дальнейшем под x0 будем понимать граничное значение диагностического параметра, обеспечивающее по правилу (5) минимум среднего риска. Будем также считать распределения f(x / D1) и f(x / D2) одномодальными («одногорбыми»).
Из условия (4) следует, что решение об отнесении объекта x к состоянию D1 или D2 можно связать с величиной отношения правдоподобия. Напомним, что отношение плотностей вероятностей распределения x при двух состояниях называется отношением правдоподобия.
По методу минимального риска принимается следующее решение о состоянии объекта, имеющего данное значение параметра x:
x D1, |
если |
|
f (x0 / D1) |
> |
(C12 −C22 )P2 |
; |
(7) |
|||
|
f (x0 / D2 ) |
|
|
|||||||
|
|
|
|
|
(C21 −C11)P1 |
|
||||
x D2, |
если |
f (x0 / D1) |
< |
(C12 −C22 )P2 |
. |
(8) |
||||
|
|
|||||||||
|
|
f (x0 / D2 ) |
|
(C21 −C11)P1 |
|
Эти условия вытекают из соотношений (5) и (4).
Условие (7) соответствует x < x0, условие (8) x > x0. Величина
λ = (C12 −C22 )P2 представляет собой пороговое значение для отноше-
(C21 −C11)P1
17
ния правдоподобия. Напомним, что диагноз D1 соответствует исправному состоянию, D2 – дефектному состоянию объекта; C21 – цена ложной тревоги; C12 – цена пропуска цели (первый индекс – принятое состояние, второй – действительное); C11 < 0,C22 – цены правильных решений (условные выигрыши). В большинстве практических задач условные выигрыши (поощрения) для правильных решений не вводятся и тогда
λ = C12 P2 / C21P1. |
(9) |
Часто оказывается удобным рассматривать не отношение правдоподобия, а логарифм этого отношения. Это не изменяет результата, так как логарифмическая функция возрастает монотонно вместе со своим аргументом. Расчет для нормального и некоторых других распределений при использовании логарифма отношения правдоподобия оказывается несколько проще.
Рассмотрим случай, когда параметр x имеет нормальное распределение при исправном D1 и неисправном D2 состояниях. Рассеяние параметра (величина среднеквадратичного отклонения) принимается одинаковым.
В рассматриваемом случае плотности распределений
|
|
|
|
|
|
|
( x− |
|
|
)2 |
|
|
||
|
|
|
|
|
− |
x |
|
|
||||||
|
|
|
1 |
|
|
1 |
|
|
||||||
f (x / D1) = |
|
|
|
e |
|
|
2σ2 |
; |
||||||
σ |
2π |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
( x− |
|
)2 |
|
|||
|
|
|
|
|
|
− |
x |
|||||||
|
|
|
1 |
|
|
2 |
|
|
||||||
f (x / D2 ) = |
|
|
e |
|
|
2σ2 |
. |
|||||||
σ |
2π |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Внося эти соотношения в равенство (4), получаем после логарифмирования
ln |
f (x0 / D1) |
= − |
1 |
[2x |
|
|
σ2 |
||||
|
f (x |
/ D ) |
|
0 |
|
0 |
2 |
2 |
|
Из этого уравнения
|
= |
1 |
|
|
|
+ |
|
) − |
x |
|
|
x |
|
x |
|||
|
||||||||
0 |
|
2 |
|
1 |
2 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
(C12 −C22 )P2 . |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
(x2 − x1) + x1 |
− x2 |
] = ln |
|
|
|
|
|
|||||||||||||||||
(C21 −C11)P1 |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
s2 |
P |
|
|
|
(C |
|
−C |
22 |
) |
|||||||||||||
|
|
|
|
|
|
|
(ln |
2 |
|
+ln |
|
12 |
|
|
|
. |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
− |
|
|
P |
|
|
|
(C |
|
−C |
) |
||||||||||||
x |
x |
|
|
21 |
||||||||||||||||||||
2 |
1 |
|
1 |
|
|
|
|
|
|
|
|
11 |
|
|
18
При x < x0, x D1; при x > x0 x D2. П р и м е р
Диагностика состояния трансмиссии газотурбинного двигателя осуществляется по содержанию железа в масле. Для исправного состоя-
ния среднее значение составляет x1 = 5 (5 г железа на 1 т масла) и среднеквадратичное отклонение σ1 = 2. При наличии дефекта подшипников и других деталей (неисправное состояние) эти значения равны
x2 = 12, σ2 = 3 . Распределения предполагаются нормальными. Требуется определить предельное содержание железа в масле, выше
которого двигатель подлежит снятию с эксплуатации и разборке (во избежание опасных последствий). По статистическим данным неисправное состояние трансмиссий наблюдается у 10% двигателей.
Примем, что отношение стоимостей пропуска цели и ложной трево-
ги C12 = 20 , и откажемся от «вознаграждения» правильных решений
C21
(C11 = C22 = 0). Из условия (4) получаем
f (x0 / D1) |
= 20 |
0,1 |
= 2, 22 . |
f (x0 / D2 ) |
|
||
0,9 |
|
Плотности распределения
|
|
|
|
|
1 |
|
|
− |
( x−5)2 |
|
|
||
f (x0 |
/ D1) = |
|
|
|
22 |
2 |
|
|
|||||
|
|
e |
|
|
; |
||||||||
2 |
|
|
|
||||||||||
|
|
|
|
|
π |
|
|
|
|
|
|||
|
|
|
|
1 |
|
|
e− |
( x−12)2 |
|||||
f (x |
/ D |
) = |
|
|
232 |
. |
|||||||
|
|
|
|||||||||||
0 |
|
2 |
3 2π |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
Внося эти значения в предыдущее равенство, получаем после логарифмирования
|
(x −5)2 |
|
(x −12)2 |
|
2 |
2, 22 |
|
− |
0 |
+ |
0 |
= ln |
|
|
. |
8 |
18 |
|
|
||||
|
|
|
|
3 |
Это уравнение имеет положительный корень x0 = 7,456.
Практическая часть
1. Изучить методические указания и получить задание.
19
2.Рассчитать предельное значение диагностического параметра, выше которого исследуемый объект подлежит снятию с эксплуатации по методу минимального риска.
3.Оформить отчет о практической работе.
4.Защитить отчет о практической работе при собеседовании с преподавателем.
Отчет должен содержать:
1.Цель работы.
2.Задание.
3.Основные формулы и положения.
4.Расчет указанного предельного значения диагностического параметра (численный).
5.Выводы по работе.
Практическая работа № 3
МЕТОД МИНИМАЛЬНОГО ЧИСЛА ОШИБОЧНЫХ РЕШЕНИЙ
Цель работы: изучение метода минимального числа ошибочных решений для диагностики технического состояния исследуемых систем и объектов.
МЕТОДИЧЕСКИЕ УКАЗАНИЯ
Вероятность ошибочного решения определяется как
∞ |
x0 |
|
Pош = P1 ∫ f (x / D1)dx + P2 |
∫ f (x / D2)dx . |
(1) |
x0 |
−∞ |
|
Из условия экстремума этой вероятности получаем
dPош |
= −P f (x |
/ D ) + P f (x / D ) = 0. |
(2) |
||||
|
|||||||
1 |
0 |
1 |
2 |
0 |
2 |
|
|
dx0 |
|
|
|
|
|
|
Условие минимума дает
d |
2 P |
= −P f '(x |
/ D ) + P f '(x |
/ D ) > 0 |
|
|||
|
ош |
(3) |
||||||
|
dx2 |
|||||||
|
1 |
0 |
1 |
2 |
0 |
2 |
|
|
|
0 |
|
|
|
|
|
|
|
20