Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

УчебПособие (Теория надежности)2011

.pdf
Скачиваний:
105
Добавлен:
07.03.2015
Размер:
2 Мб
Скачать
☆

 

 

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