УчебПособие (Теория надежности)2011
.pdf
|
|
1 |
|
|
t i |
. |
|||
|
||||
n |
||||
|
|
qi j |
(5.11) |
|
|
|
|
||
|
|
j i |
|
|
|
|
j 1 |
|
|
При помощи уравнений (5.4), (5.10), (5.11) и условия нормировки можно вычислить все характеристики состояний, исходя из интенсивности переходов, определяющих данную систему.
В заключение данного раздела приведем алгоритм определения вероятностей нахождения объекта в каждом из возможных состояний
в любой момент времени ( Pi (t)
каждого состояния распределены по 1.Составляется диаграмму
временной диаграмме параметрами qi j .
) (при этом продолжительности
экспоненциальному закону): пространства состояний по состояний с неизвестными
2.Из временной диаграммы состояний определяются параметры
ti , i и i j .
3.Определяются асимптотические значения вероятностей
нахождения объекта в каждом из состояний Pi () ti i . 4. Определяются интенсивности переходов изодного состояния
|
в другое q |
|
i j |
. |
|
|
||
|
|
|
|
|||||
|
|
i j |
|
Pi |
( ) |
|
|
|
|
|
|
|
|
|
|||
5. |
Составляются уравнения Чемпена-Колмогорова и условия |
|||||||
|
нормировки. |
|
|
|
|
|
|
|
6. |
Из |
полученной |
системы |
уравнений |
находятся |
|||
вероятности Pi (t) .
5.1.4. Объединение состояний
Во многих приложениях определение асимптотических (стационарных) значений вероятностей состояний с помощью моделей пространств состояний можно упростить, если объединить некоторые множества состояний для того, чтобы они образовали единичное совокупное состояние.
После объединения групп состояний возникает новый процесс с некоторыми новыми состояниями (совокупными состояниями) и с новыми переходами (между совокупными состояниями). В
71
большинстве случаев новый процесс уже не будет Марковским, потому что продолжительности пребывания системы в совокупных состояниях, вообще говоря, распределены не по экспоненциальному закону.
Новый процесс останется Марковским только при удовлетворении условия объединяемости – условия, при котором интенсивности перехода к любому другому состоянию или к группе объединенных состояний одинакова для всех состояний, входящих в рассматриваемую группу. При выполнении этого условия можно рассчитать вероятности нахождения системы в любом состоянии не только для стационарных, но и для нестационарных процессов.
Поясним сказанное выше на примере объединения двух состояний системы в одно (рис. 5.4).
Процесс останется Марковским, то есть интенсивности переходов q6 1, q6 2, и q6 3 останутся независимыми от времени при
выполнении следующих условий: q4 1 = q5 1, q4 2 = q5 2, и q4 3 = = q5 3. При этом неважно, равны или нет интенсивности q4 1, q4 2 и
q4 3.
Выведем формулы для вычисления вероятностей и частоты возникновения совокупных состояний. Рассмотрим схему (рис. 5.5), в которой несколько состояний j объединены в одно совокупное состояние J.
Вероятность состояния J, обозначаемая PJ, равна сумме всех вероятностей Pj:
PJ |
Pj . |
(5.12) |
j J |
|
|
|
|
Вероятности Pj можно складывать, так как события, заключающиеся в том, что система находится в одном из состояний j, взаимоисключающие. Частота возникновения состояний J, обозначаемая J, равна сумме всех частот перехода из состояния j в состояние i, находящееся вне J:
J |
j,i |
Pj q j i |
Pj q j i , |
(5.13) |
i J j J |
i J j J |
j J |
i J |
|
где при втором преобразовании была выполнена подстановка по формуле (5.8).
Чтобы найти решение, используя пространство состояний, получившееся после объединения состояний j, необходимо знать обозначенные на рис. 5.5 интенсивности переходов qi J и qJ i. Их можно вычислить, исходя из того, что частота переходов от i к J должна быть такой же, как частота переходов от i ко всем состояниям j
72
до их объединения, и аналогично для переходов от J к i. Согласно (5.8) эти требования можно представить как
P q |
i J |
P |
q |
i j |
; |
P q |
J i |
|
|
P q |
j i |
, |
(5.14) |
|
i |
|
i |
|
|
j |
|
j J |
j |
|
|||||
|
|
|
j J |
|
|
|
|
|
|
|
|
|
|
|
|
откуда |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
qi J |
qi j |
|
|
|
|
|
(5.15) |
|||
|
|
|
|
|
|
j J |
|
|
|
|
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Pj q j i |
|
|
|
|
|
|
|
|
|
|
q J i |
j J |
|
|
. |
|
|
|
|
(5.16) |
||
|
|
|
|
Pj |
|
|
|
|
||||||
|
|
|
|
|
|
|
j J |
|
|
|
|
|
|
|
Если удовлетворены условия объединяемости, то есть если qj i одинаковы для всех j, то уравнение (5.16) упрощается:
|
qJ i |
q j i для |
j J . |
(5.17) |
|
q4 1 |
1 |
|
|
4 |
q5 1 |
|
|
q6 1 1 |
|
q4 2 |
|
q6 2 |
|
|
|
2 |
6 |
|
|
q5 2 |
|
2 |
|
5 |
|
|
||
|
q |
q4 3 |
|
q |
|
|
|
6 3 |
|
|
5 3 |
3 |
|
|
|
|
|
|
|
|
|
|
|
3 |
Рис. 5.4. Объединение состояний 4 и 5 в состояние 6 при условии, что процесс должен остаться Марковским
j |
qi j |
|
qJ i |
i |
|
||
|
q j i |
J |
i |
|
|
|
qi J |
J
Рис. 5.5. Объединение состояний j в совокупное состояние J
73
В заключение найдем интенсивность переходов между двумя совокупными состояниями I и J, каждое из которых состоит из нескольких исходных состояний (естественно непересекающихся). Схематически такая система изображена на рис. 5.6; согласно уравнениям (5.15) и (5.16) эти интенсивности равны:
|
Pi |
qi j |
|
|
|
Pj q j i |
|
||||||
qI J |
i I |
j J |
|
, |
qJ I |
j J |
i I |
|
. |
(5.18) |
|||
|
Pi |
|
|
|
Pj |
|
|||||||
|
i I |
|
|
|
|
|
|
|
j J |
|
|
|
|
Если удовлетворены условия объединяемости, то эти уравнения |
|||||||||||||
приводятся к виду |
|
|
|
|
|
|
|
|
|
|
|
|
|
qJ I q j i |
для |
j J ; qI J qi j |
для i I , |
(5.19) |
|||||||||
i I |
|
|
|
|
|
|
j J |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||
поскольку внутри I все qi j |
не зависят от i, а внутри J все qj i не |
||||||||||||
зависят от j. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
qI |
|
|
|
|
|
|
|
|
|
|
|
|
|
J |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
i |
|
|
|
|
qJ |
|
|
|
j |
|
|
|
|
|
|
|
|
|
I |
|
|
|
|
||
|
|
|
|
I |
|
|
|
|
|
J |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Рис. 5.6. Переходы между двумя совокупными состояниями
Пример 1. Рассмотрим пространство состояний дублированной системы с неявным резервированием при ограниченном аварийном ремонте (возможен только их последовательный ремонт). Предположим, что оба элемента имеют одинаковую интенсивность отказов и восстановления . Пространство состояний тогда приобретает вид, показанный на рис. 5.7, а. После объединения двух состояний, характеризуемых отказом одного элемента в состоянии А, и двух состояний, характеризуемых отказом двух элементов в состоянии В, получается диаграмма пространства состояний, изображенная на рис. 5.7, б.
Легко убедиться в том, что при этом объединении состояний условия объединяемости удовлетворяются.
Для состояния А:
q1 0 q2 0 и q1 В q2 В – условие объединяемости выполняется.
74
Для состояния B:
q3 A q4 A – условие объединяемости выполняется.
Вычислим интенсивности переходов в новой диаграмме. Согласно уравнению (5.15)
q0 A 2 .
Согласно (5.16) или (5.17)
|
|
(P |
P ) |
qA 0 |
|
1 |
2 |
P |
P |
||
|
|
1 |
2 |
И наконец, согласно (5.18) или (5.19)
|
(P P ) |
|
|
qA B |
1 |
2 |
и qB A |
P |
P |
||
|
1 |
2 |
|
.
(P P ) |
|
||
3 |
4 |
. |
|
P P |
|||
|
|||
3 |
4 |
|
|
a отказал
|
|
|
|
|
|
|
|
|
|
|
a отказал, затем b |
||||
|
|
|
|
|
1 |
|
|
|
|
|
3 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
отказал |
|||
Нет |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
отказов |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
2 |
|
|
4 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
b отказал, затем а |
||||||
|
|
|
|
b отказал |
отказал |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
а) |
|
|
|
|
|
||
Нет отказов |
|
Отказ одного из |
|
|
Отказ двух |
||||||||||
|
|
элементов |
|
|
элементов |
||||||||||
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
A |
|
|
|
|
B |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
б) |
|
|
|
|
|
||
Рис. 5.7. Объединение состояний (1 2 A , |
3 4 B ) дублированной системы с |
||||||||||||||
неявным резервированием при условии возможности лишь последовательного их ремонта: a – исходная диаграмма пространства состояний; б – конечная диаграмма пространства состояний
Модель, изображенная на рис. 5.7, б, очевидно, значительно проще использовать для расчета, чем модель 5.7, а. Однако некоторая информация, имеющаяся в исходной диаграмме, маскируется в
75
упрощенной диаграмме, например тот факт, что одни ремонты элементов могут выполняться лишь последовательно, становится неразличимым.
Иногда может оказаться целесообразным объединить состояния, если даже условия объединяемости и не выполняются. В этих случаях уравнения, получаемые из диаграммы после объединения состояний, не могут быть использованы для определения неустановившихся значений вероятностей состояний. Несмотря на это, асимптотические значения вероятности состояний и частоты возникновения состояний вычисляются верно.
Пример 2. Рассмотрим пространство состояний, изображенное на рис. 5.7, а. Если объединить состояния 3 и 4, но не объединять состояния 1 и 2, то получится схема, изображенная на рис. 5.8.
|
|
|
a отказал |
|
|
|
1 |
|
|
|
|
|
|
|
Нет отказов |
|
|
2 |
Оба элемента |
|
|
|
|
отказали |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
4 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
b отказал |
|
Рис. 5.8. Объединение состояний, отвечающий отказу двух элементов (рис. 5.7, а)
При таком объединении условия объединяемости не выполняются:
(q3 2 ) (q4 2 0) и (q3 1 0) (q4 1 ) .
Но, несмотря на это, асимптотические значения вероятности состояний вычисляются верно. Интенсивности переходов для диаграммы, изображенной на рис. 5.9, находятся следующим образом.
76
Согласно выражению (5.15)
q1 B q2 B .
Согласно (5.16) |
|
|
|
|
|
|
P P 0 |
|
|
|
|
qB 1 |
4 |
3 |
|
|
. |
(P P ) |
2 |
||||
|
3 |
4 |
|
|
|
Причем вероятности сокращаются из соображений симметрии, так как P3=P4. Аналогично
|
P 0 P |
|
|
|
|
qB 2 |
4 |
3 |
|
|
. |
(P P ) |
2 |
||||
|
3 |
4 |
|
|
|
Пример 3. Временная диаграмма состояний системы, которая может находиться в пяти состояниях, представлена на рис. 5.9.
S
S5
S4
S3
S2 |
|
|
|
|
|
S1 |
t1 |
t2 |
t3 |
t4 t5 |
8760 t (ч) |
|
|||||
|
|
Рис. 5.9. Временная диаграмма состояний системы |
|
||
Известно следующее:
-продолжительность пребывания системы в каждом из состояний подчинена экспоненциальному закону;
-средняя продолжительность пребывания системы в каждом из состояний:
t1 500 ч; t2 |
100 ч; t3 50 ч; t4 200 ч; t5 377,5 ч; |
- S1 S2 S6 |
– рабочее состояние; S3 S4 S7 – состояние |
аварийного ремонта; S5 – состояние планово-предупредительного ремонта.
Задание:
1. Составить диаграмму пространства состояний по данной временной диаграмме состояний системы.
Определить:
- частоту появления каждого из состояний ( 1, 2, 3, 4, 5);
77
-вероятность нахождения системы в каждом из пяти состояний
встационарном режиме работы (P1( ), P2( ), P3( ), P4( ), P5( ));
-интенсивности переходов из состояния в состояние (q1 2, q1 3,
q3 2, q4 5 и т. д. – год-1).
2. Составить диаграмму пространства состояний для данной
системы, которая имеет три состояния: S5, S6, S7. Определить с помощью этой диаграммы вероятности нахождения системы в работоспособном состоянии, состоянии аварийного ремонта и в состоянии планово-предупредительного ремонта для установившегося (стационарного) процесса.
Решение:
Составим диаграмму пространства состояний по данной временной диаграмме состояний системы (рис 5.10).
|
|
1 (q3 2 q3 4) t |
|
q1 3 t |
S3 |
1 (q1 3 q1 2 q1 5) t |
|
|
|
q3 2 t |
q3 4 t |
|
|
|
|
S1 |
|
|
|
1 q4 2 t |
|
q2 3 t |
|
|
q2 1 t |
S4 |
q |
t |
|
1 2 |
|
|
|
|
q4 2 t |
q1 5 t
S2
q5 2 t
1 (q2 1 q2 3 q2 5) t |
1 q5 2 t |
|
q2 5 t S5
Рис. . 5.10. Диаграмма пространства состояний системы
78
Определим частоту появления каждого из состояний.
n
i i j ; j 1
1 1 2 1 3 1 5 4 3 2 9 ;2 2 1 2 3 2 5 9 6 2 17 ;3 3 2 3 4 6 3 9 ;4 4 2 3 ;5 5 2 4 .
Определим вероятности нахождения системы в каждом из пяти состояний в стационарном режиме работы.
P ( ) ti . |
|
|
|
|
|
|
|
|
|
|
|
|||
i |
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
P ( ) t1 |
9 500 |
0,514 ; |
||||||||||||
|
|
|
|
|
|
|||||||||
1 |
1 |
|
8760 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
||||||
P ( ) |
|
t2 |
|
17 100 |
0,194 ; |
|||||||||
2 |
|
|
|
|
||||||||||
2 |
|
|
8760 |
|
|
|
||||||||
|
|
|
|
|
|
|
||||||||
P ( ) |
|
t3 |
|
9 50 |
0,051 ; |
|||||||||
|
|
|
|
|||||||||||
3 |
3 |
|
8760 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
||||||
P ( ) |
|
t4 |
|
3 200 |
0,068 ; |
|||||||||
4 |
|
|||||||||||||
4 |
|
|
8760 |
|
|
|
||||||||
|
|
|
|
|
|
|
||||||||
P ( ) |
|
t5 |
|
4 377,5 |
0.172 . |
|||||||||
|
|
|||||||||||||
5 |
5 |
|
8760 |
|
|
|||||||||
|
|
|
|
|
|
|||||||||
Определим интенсивности переходов из состояния в состояние.
qi j ( ) i j .
P ( ) i
q |
( ) 1 2 |
|
4 |
|
|
7, 782 ; |
||||
|
|
|
||||||||
1 2 |
|
P ( ) |
|
|
|
0,514 |
|
|
||
|
|
1 |
|
|
|
|
|
|
|
|
q |
( ) 2 1 |
|
|
9 |
|
|
46,392 ; |
|||
|
|
|
|
|||||||
2 1 |
|
P2 ( ) |
|
|
|
0,194 |
|
|
||
|
|
|
|
|
|
|
||||
q |
( ) 1 3 |
|
|
3 |
|
5,832 ; |
||||
|
|
|
|
|||||||
1 3 |
|
P ( ) |
0,514 |
|
|
|
||||
|
|
1 |
|
|
|
|
|
|
|
|
q |
|
( ) 2 3 |
|
6 |
|
|
30,928 ; |
|||
|
|
|
|
|
||||||
2 3 |
P2 ( ) |
0,194 |
|
|||||||
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
79 |
|
q |
|
|
|
( ) 1 5 |
|
|
|
|
|
2 |
|
3,891; |
|||||||||
|
|
|
|
|
|
|
|
||||||||||||||
1 5 |
|
|
|
P ( ) |
0,514 |
|
|
|
|
||||||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
q |
|
|
|
|
( ) 2 5 |
|
2 |
|
|
|
10,309 ; |
||||||||||
|
|
|
|
|
|
|
|||||||||||||||
2 5 |
|
|
P2 ( ) |
|
|
|
|
|
0,194 |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
q |
|
|
|
|
( ) 3 2 |
|
|
|
|
6 |
|
|
|
117, 647 ; |
|||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
3 2 |
|
|
P ( ) |
|
|
|
|
|
0,051 |
|
|
||||||||||
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
q |
|
|
|
|
( ) 4 2 |
|
3 |
|
|
|
44,118 ; |
||||||||||
|
|
|
|
|
|
|
|
||||||||||||||
4 2 |
|
|
|
P4 ( ) |
|
|
|
|
|
0,068 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
q |
|
|
|
|
( ) 3 4 |
|
|
|
3 |
|
|
|
58,824 ; |
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
3 4 |
|
|
P ( ) |
|
|
|
|
|
0,051 |
|
|
||||||||||
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
q |
|
|
|
|
( ) 5 2 |
|
|
4 |
|
|
|
23, 256 . |
|||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
5 2 |
|
|
P ( ) |
|
|
|
|
|
0,172 |
|
|
||||||||||
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
Проведем проверку правильности определения значений |
|||||||||||||||||||||
интенсивностей перехода системы из i-го состояния. |
|||||||||||||||||||||
q |
|
8760 |
|
|
8760 |
|
17,52 год-1; |
||||||||||||||
|
|
|
|
||||||||||||||||||
1 |
|
|
|
|
t1 |
500 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
q1 q1 2 q1 3 q1 5 |
7,782 5,837 3,891 17,51 год-1. |
||||||||||||||||||||
q |
|
8760 |
|
8760 |
|
87,6 год-1; |
|||||||||||||||
|
|
||||||||||||||||||||
2 |
|
|
|
|
t2 |
100 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
q2 q2 1 q2 3 q2 5 46,392 30,928 10,309 87,629 год-1.
q |
8760 |
|
8760 |
175, 2 год-1; |
|
|
|||
3 |
t3 |
50 |
|
|
|
|
|||
q3 q3 2 q3 4 117,647 58,824 176, 471 год-1.
q4 8760 8760 43,8 год-1; t4 200
q4 q4 2 44,118 год-1.
q |
8760 |
|
|
8760 |
23, 205 год-1; |
|
|
||||
5 |
t5 |
|
377,5 |
|
|
|
|
|
|||
q5 q5 2 |
23, 256 год-1. |
||||
80
