О чувствительности характеристик надежности систем к виду функций распределения времени безотказной работы и восстановления их элементов

УДК 519.718: 658.652

О чувствительности характеристик надежности систем к виду функций распределения времени безотказной работы и восстановления их элементов

В. В. Рыков*, Ань Нгиа Ча^

* РГУ нефти и газа имени И.М. Губкина Ленинский просп., д. 65, корп. 1, Москва, Россия, 119991 ^ Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198

В статье рассматривается проблема чувствительности характеристик надёжности систем (M2/GI/1} и (GI2/M/1) к виду функций распределения (ф.р.) времени безотказной работы (в.б.р.) и времени восстановления их элементов при ограничениях на доступность восстановления. Для этих систем представлены дифференциальные уравнения в частных производных для нестационарных и обыкновенные дифференциальные уравнения для стационарных вероятностей микросостояний состояний. Для стационарных вероятностей микро- и макросостояний и получены аналитические выражения их зависимости от вида распределений в.б.р. и времени восстановления, которые явно зависят от вида ф.р. не показательно распределённых исходных характеристик систем через их производящие функции (преобразования Лапласа—Стилтьеса соответствующих ф.р.) в точках, равных интенсивности показательно распределённой характеристики.

С помощью специально разработанного программного средства в среде MATLAB проведено численное исследование чувствительности вероятности отказа системы от вида функций распределения времени безотказной работы и восстановления её элементов и сравнение полученных результатов с соответствующими характеристиками для простейшей марковской модели (М2/М/1). Проведённое исследование показало, что эта чувствительность незначительна и становится исчезающе малой при «быстром» восстановлении. В частности, в результате численного анализа с помощью указанного программного средства показано, что когда в качестве общего распределения GI используются Гамма-распределение (Г) или распределение Вейбулла-Гнеденко (W — G), вероятности отказа систем (M2/GI/1} и (GI2/М/1) быстро сходятся к нулю с ростом скорости восстановления.

Ключевые слова: надёжность систем, вероятности отказа, чувствительность к виду функций распределения, микро- и макросостояния, стационарные и не стационарные вероятности.

1. Введение

Одним из первых результатов о нечувствительности характеристик систем массового обслуживания к виду ф.р. времени обслуживания была теорема Б. А. Севастьянова [1] о справедливости формул Эрланга для произвольных законов распределения длительностей разговоров. И. Н. Коваленко [2] показал, что при показательных распределениях времени безотказной работы (в.б.р.) элементов резервированной системы необходимым и достаточным условием нечувствительности стационарных характеристики её надёжности (стационарных вероятностей состояний) к виду ф.р. времени восстановления является возможность немедленного начала восстановления отказавшего элемента (т.е. достаточное количество восстанавливающих устройств). В [3] показано, что это условие достаточно для нечувствительности стационарных характеристик надёжности систем к виду как ф.р. в.б.р., так и восстановления. Однако при наличии ограничений на возможность восстановления (в случае ограниченного числа восстанавливающих устройств) этот факт не верен. Теоретическое исследование нечувствительности стационарных характеристик надёжности к виду ф.р. в.б.р. и восстановления в случае ограниченного числа восстанавливающих устройств — не простая задача.

В настоящей статье исследуется влияние вида ф.р. в.б.р. и времени восстановления на стационарные характеристики надёжности систем при ограниченном

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

Не смотря на явную зависимость стационарных вероятностей состояний системы от вида ф.р. в.б.р. и времени восстановления численный анализ показал, что при «быстром» восстановлении эта зависимость становится исчезающе малой.

В следующем разделе вводятся обозначения и приводятся характеристики надёжности для простейшей марковской модели холодного дублирования с одним восстанавливающим устройством. Затем в двух последующих разделах исследуются модели холодного резервирования в случае, когда только одно из распределений в.б.р. или времени восстановления отлично от показательного. Наконец в последнем разделе проводится соответствующие численные исследования и приводится сравнительный анализ с марковским случаем.

Рассмотри систему надёжности из п элементов, подверженных отказам (источников заявок), и т ремонтных единиц. Обобщая символику Кендалла, такие системы будем обозначать через {GIn/GI/m(RE)), где угловые скобки { ) означают, что рассматривается замкнутая система обслуживания, с ограниченным числом источников заявок. Символ GI (General Independent) используется для указания рекуррентного потока отказов на первом месте или рекуррентного механизма восстановления элементов на втором и заменяется символом М для показательных распределений в.б.р. или времени восстановления или специальными символами для иных распределений, например Ek(X) для распределения Эрлан-га с параметрами (к, Л). Наконец, последний символ указывает число ремонтных единиц (обслуживающих приборов). Таким образом, дублированная система с одним ремонтным прибором и произвольными законами распределения в.б.р. и восстановления обозначается как {GI2/GI/I). При этом необходимо дополнительно оговаривать вид резервирования (горячее, холодное или тёплое).

Будем обозначать ф.р. в.б.р. и восстановления элементов системы через А(х) и В(х) соответственно. Предполагается, что они абсолютно непрерывны, а их производные, плотности распределения (п.р.), обозначаются через а(х) и Ь(х). Соответствующие средние значения и интенсивности отказов и восстановлений обозначаются

Наконец для производящей функции моментов (п.ф.м.) в.б.р. и восстановления (преобразований Лапласа их плотностей а(х) и Ь(х))) используются обозначения

2. Постановка задачи и обозначения

a

!

е 3Xa

!

е sxb(x)dx.

Для получения достаточно прозрачных численных результатов остановимся на системах холодного дублирования, когда п = 2 с одним восстанавливающим устройством т = 1, и, чтобы исследовать чувствительность модели по виду ф.р.

в.б.р. и времени восстановления для различных моделей {М2/СИ/1) и (СИ2/М /1) и сравнить их с соответствующими характеристиками для простейшей марковской модели (М2/М/1), напомним сначала стационарные характеристики для этой модели.

3. Марковская модель системы с ограниченным числом

ремонтных устройств

Для простейшей марковской модели (M2/M/I) случайный процесс

N (t) = число отказавших элементов в момент времени t

является процессом гибели и размножения с множеством состояний Е = {0, 1, 2}. Для вероятностей его состояний ^í(t), ni(t) = (t) = г} справедливы дифференциальные уравнения Колмогорова, в которых верхней точкой обозначаются производные по времени

-к o(t) = —a-Ko(t) + ¡3^1 (t),

кi

-к2 (t) = аж\ (t) - ¡3-K2 (t).

Решение этой системы следует искать с начальным условием, определяемым начальным состоянием процесса.

Стационарные вероятности состояний = lim^,^ ^í(t), существование которых следует из конечности множества состояний процесса и его неразложимости, удовлетворяют соответствующей системе уравнений глобального баланса

0 = —ал0 +

0 = amo -

0 = аж1 —

а

с условием нормировки wo + + ^2 = 1. С использованием обозначения р = —

и

решение этих уравнений представимо в виде

1 Р Р2 fos

по = -7-.-¡—2, п1 = -¡—2 , п2 = т—:-¡—2.

1+ р + р2 1+ р + р2 1+ р + р2

Заметим, что вероятность отказа системы совпадает с величиной Л2.

Рассмотрим теперь системы, в которых только одно из распределений (в.б.р. или восстановления) отличны от показательного. Для исследования этих систем используется метод дополнительных переменных [4]. Ограничимся для простоты случаем с двумя приборами и одним восстанавливающим устройством.

4. Система холодного дублирования [М2/01 /1)

Рассмотрим систему {М2/СИ/1) холодного дублирования с одним ремонтным устройством. В качестве дополнительной переменной используем время, затраченное на обслуживание отказавшего элемента, и обозначим через жо(р) вероятность того, что в момент времени £ система находится в состоянии % = 0, и через х) плотность распределения (п.р.) вероятностей того, что в момент времени £ система находится в состоянии % (г = 1,2) и время, затраченное на обслуживание

отказавшего элемента, находится в интервале (х,х + ёж]

М*) = (¿) = 0},

ж^; х)ёх = (г) = ц х < X(г) < ж + ёж} (г = 1,2).

Теорема 1. Функции ж 1(1; х) удовлетворяют системе дифференциальных уравнений в частных производных

йкр

<и

ъ

= —a'к0(t) + J (Ь; и)/3(и)ёи,

джх {Ъ; х) джх {Ъ; х)

—Ъъ— + дх = —(а +13 (х))^(р;х),

X) X)

т

с граничным условием

+

дх

—¡3

г

(I; 0) = ап0(1) + J и)@(и)ёи.

Доказательство. С помощью формулы полной вероятности путём сравнения вероятностей состояний процесса в моменты времени £ и £ + А легко выводится система уравнений в конечных разностях

п0(г + А) = п0(г)(1 — аА) + ^ Ж1 (Ь; и)/3(и)ёи + о(А),

о

Ж1 (г + А; ж + А) = Ж1 (г; х)(1 — аА)(1 — @(х)А) + о(А), п2(г + А; х + А) = п2(г; х)(1 — Р(х)А) + щ(г; х)аА + о(А) с граничным условием

жх(г + А; А)А = ж0(Ь)аА + А ! ж2(Ь; и)/3(и)ёи.

После несложных алгебраических преобразований и перехода к пределу при А ^ 0 последняя система переходит в систему дифференциальных уравнений в частных производных (4) с граничным условием (5).

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

Следствие 1. Для вероятностей состояний системы в стационарном режиме соответствующие уравнения примут вид

ак0 = ^ (и)Р(и)ёи,

о

К1

с граничным условием

с

^1

о

Решение этой системы уравнений содержится в следующей теореме. Теорема 2. Решение системы уравнений (6) имеет вид

ж0 = С1а~1 Ъ(а), Ж1 (х) = С1&~ ах(1 — В (х)), (х) = С1 (1 — е-ах) (1 — В(х)),

с некоторой постоянной С1, которая будет определена позже исходя из условия нормировки.

Доказательство. Решение второго из уравнений (6) имеет вид

Ж1 (х) = С1&~ ах(1 — В (х)), С1 = ^1(0).

Подставляя это решение в первое из уравнений и используя для преобразования Лапласа обозначение Ь(в) = / е~вхЬ(х)ёх, находим

к0 = С1а~ 1Ъ(а).

Наконец, решение однородной части третьего из этих уравнений даёт

ъ(х) = С 2(1 — В(х)), откуда методом вариации постоянной найдём

С2(х)(1 — В(х)) — С2(х)Ь(х) = 1 —Ь(х\ С2(х)(1 — В(х)) + аж1(х),

1 — В

что после сокращений даёт

С2 (х) = аСх е~ откуда интегрированием найдём

С2(х) = С1(1 — е~ ах) + С 2,

что окончательно позволяет представить ,К2(х) в виде

*2(х) = [С1(1 — е~ах) + С2](1 — В(х)).

Наконец, граничное условие (7) даёт

сю

С1 = ^1(0) = ак0 + I п2(и)/3(и)ёи = С1Ь(а) + С1(1 — Ь(а)) + С2,

0

откуда следует, что С2 = 0 и, следовательно,

П2(х) = С1 (1 — е") (1 — В(х)).

Для вероятностей макросостояний ж0, = / ^(х^х (г = 1, 2) путём интегрирования получим утверждение теоремы.

Следствие 2. Вероятности макросостояний равны

Ь

а а а

Для определения постоянной С\ воспользуемся условием нормировки

1 = Сл

а-Ч

Из последнего соотношения с использованием обозначения р = аЬ получим следующее следствие.

Следствие 3. Постоянная С\ определяется соотношением

а

С\ = .

р + о

Объединяя полученные результаты для стационарных макровероятностей состояний системы {М2/С1/1), получим выражения, представленные в следующем следствии.

Следствие 4. Вероятности макро-состояний равны

Ь

К0 = -=-, = -=-, К2 = -=-.

р + Ь

Замечание 1. Зависимость вероятностных характеристик системы от вида ф.р. времени ремонта носит достаточно сложный характер и выражается в терминах её преобразования Лапласа (производящей функции).

Замечание 2. При показательном распределении в.б.р. В(х) = 1 — элементов Ь(а) = ¡3(а + Р)-1 постоянная С\ имеет вид

= Р-Р

1 + р + р2 '

а распределение вероятностей состояний системы имеет вид

1

ж0 = С\а 1Ъ

1 + р + р2 '

= Сга-\1 - Ь(а))= 9

1 + р + р2' п2

Ж2 = С\а-1(аЬ - (1 - Ь(а)))= И

1 + р + р2

и совпадает с приведённым ранее для марковского случая (3).

5. Система холодного резервирования {С12/М/1)

Рассмотрим теперь аналогичную систему с произвольно распределённым в.б.р. элементов и показательно распределённым временем ремонта, при этом для простоты предполагается холодное резервирование, так что одновременно может работать только один элемент.

Используя в качестве дополнительной переменной X(1) время, прошедшее с момента очередного начала работы прибора, и вводя п.р. вероятностей состояний

ж^; х)ёх = Р{^(г) = Г; х < X(г) < ж + ёж} (г = 0,1), Ж2(1) = Р{^ (I) = 2},

аналогично предыдущему случаю для этих функций несложно получить систему дифференциальных уравнений в частных производных.

Теорема 3. Функции к^; х) удовлетворяют системе дифференциальных уравнений в частных производных

{ дп0

+--^-= —а(х)ж0(ц х) + рж1(Г; х),

д1 дх

д^1

т

dK2 (¿) <И

+

дх

—(а(х) + х),

ь

—Ръ2(1) + I а(и)п1(1; и)ёи

с граничным условием

г

0) = J п0(1; и)а(и)ёи + Рк2(1).

Доказательство. Доказательство проводится аналогично теореме 1 с помощью формулы полной вероятности путём составления конечно-разностных уравнений и сравнения вероятностей состояний процесса в моменты времени £ и £ + А с последующим переходом к пределу при А ^ 0.

Следствие 5. Для вероятностей состояний системы в стационарном режиме соответствующие уравнения примут вид

тг0(х) = —а(х)к0 (х) + Р,к1(х), тг 1(х) = —(а(х) + Р )к1(х),

Ъ2 = Р

1

к1(и)а(и)ёи

с граничным условием

с

п1(0) = J п0(и)а(и)ёи + /Зп2.

Следующая теорема содержит решение этой системы уравнений. Теорема 4. Решение системы уравнений (13, 14) имеет вид

П0

с некоторой постоянной С1, которая будет определена позже, исходя из условия нормировки.

Доказательство. Решение второго из этих уравнений имеет вид т(х) = С1е-13х(1 - А(х)), С1 = яч(0). При этом из третьего уравнения сразу следует, что

Ж2 = Сф-1а(р).

Наконец, однородная часть первого уравнения даёт по(х) = С0(1 - А(х)), а использование метода вариации постоянной приводит к следующему уравнению для Со(х) :

Со(х)(1 - А(х)) - Со(х)а(х) = -а(х)Со(1 - А(х)) + рС1е-/3х(1 - А(х)),

откуда следует Со(х) = С1@е-13х или С0(х) = -С1е-/Зх + С0. Таким образом имеем

Мх) = (Со - С1е-^х)(1 - А(х)).

Для определения постоянной Со воспользуемся граничным условием (14), которое принимает вид

С1 = т(0) = 1[по(х) + т(х)]а(х)<х = !\(Со - С1е-13х) + С1е-/Зх]а(х)<х = Со

и позволяет выразить постоянную С0 в виде С0 = С1. Отсюда для л0 (х) окончательно найдём

жо(х) = С 1(1 - е-^х)(1 - А(х)).

Наконец, для вероятностей макросостояний путём интегрирования получим следующее следствие.

Следствие 6. Вероятности макросостояний определяются формулами

/ сю

1 - аф)~

к0 = I к0(х)Ах = С1 а -

о

с

т = ! т(х)<х = С1р-1(1 - а(р)),

Ж2 = Сф-1а(0).

Для определения постоянной С1 следует воспользоваться условием нормировки

1 = жо + т + П2 = Сф-1 [ар - (1 - аф)) + (1 - а(р)) + а(р)] = Сф-1[ар + а(р)].

Следствие 7. Постоянная С1 имеет вид

где р-1 = а/3.

Объединение соотношений (15, 16) приводит к окончательному результату в форме следующего следствия.

Следствие 8. Вероятности макросостояний определяются формулами

^о =

= а(р)(ар + а(Р))-1 = а(р)(р-1 + а(р))-1.

Замечание 3. Зависимость вероятностных характеристик системы от вида ф.р. времени жизни носит достаточно сложный характер и выражается в терминах её производящей функции в точке, равной интенсивности ремонта.

Замечание 4. При показательном распределении в.б.р. А(х) = 1 - е-ах элементов

а(р) = а(а + р)-1 =

и распределение вероятностей состояний системы совпадает с приведённым ранее для марковского случая (3):

ар - (1 - а(Р)) р-1 - (1 + р)-1 1

ко = ~ ~

р-1 + а(Р) р-1 + р(1 + р)-1 1 + р + р2

а(Р) Р2

_ = 1 - а

К2 =

р-1 + а(р) 1 + р + р2'

6. Численный анализ систем

Для численного анализа рассмотрим два примера, когда в качестве общих распределений СИ используются Гамма-распределение (Г) с функцией плотности

. , Х(М)С-1 „

= Мге

и функция распределения Вейбулла-Гнеденко (Ш - С) с функцией плотности

р(1) = сХЧс-1е-(Х1)°,

где с — параметр формы, Л — параметр скорости.

Подстановка этих распределений в общие формулы (10) и (17) даёт следующие результаты.

Теорема 5. 1. Для системы {М2/Г/1)

МР) = 1__^ = 1 - 13(1 + <* )С

(р + Ь(а)) (1 + ср)с + ссрс+1'

где Р = Х/с, р = аЬ = 1/Р, а = а = 1. 2. Для системы {Г2 /М/1)

* (я) й(13) К2(Р) = —

р-1 + а(р) сс + р (с + РУ

сс

где X = с, р-1 = аР = Р, а = а = 1. 3. Для системы - С/1)

СЮ

Ь(а) = ! с(рд)сгс-1е-(г+(^дг)С)сИ,

о

МР) = 1 -7-^^ = 1 - ^

(р + Ь(а)) (1 + р~Ь(а)У

где д = Г(1 + 1), Р = Х/д, р = аЬ = 1/Р, а = а = 1.

4. Для системы — С2/М/1)

С

а(Р) = У сдЧ^еГ^+^сИ,

о

Ж2(в) = ~а(Р) = ~а(Р) 2(Р) р-1 + а(р) р + а(Р),

где д = Г(1 + 1), X = д, р-1 = аР = Р, а = а = 1.

Результаты численного анализа приведены в виде таблиц и графиков. В табл. 1 и 2 приведены значения вероятностей отказа п2(Р) систем (М2/Г/1) и (Г2/М/1) при различных значениях параметров с и р.

На рис. 1-4 представлены графики функций вероятностей отказа ^(Р) систем различного вида в зависимости от скорости ремонта Р при различных значениях параметра с.

Таблица 1

Стационарные вероятности отказа к2

Ъ(Р) с = 0,1 с = 0, 5 с = 1,0 с = 55,0 с =10

р = Х/с =1 0,4403 0, 3660 0, 3333 0, 2867 0,2783

р = Х/с = 10 0, 0320 0, 0127 0, 0090 0,0057 0,0053

р = Х/с = 100 5,139910-4 1,475210-4 9, 9  10-5 5,9718 10-5 5,4778 10-5

<М2/0атта/1>

1 г..............;................г...............;................г...............1................

°0 123456789 10 |3 = X / С

Рис. 1. Функции вероятности отказа ж2

Таблица 2

Стационарные вероятности отказа ) системы

Ъ(Р ) с = 0,1 с = 0, 5 с = 1, 0 с = 5, 0 с =10

р = 1(Л = с) 0,4403 0, 3660 0,3333 0,2867 0,2783

Р = 10(А = с) 0, 0593 0,0214 0, 009 4,1135  10-4 9, 7647  10-5

Р = 100(А = с) 0, 005 7, 0485  10-4 9, 9  10-5 2,448510-9 3,855410-13

Результаты численного анализа показывают, что вероятности отказа систем (М2 /Г/1) и (Г2/М/1) мало отличаются друг от друга и с ростом параметра /3 быстро сходятся к нулю тем быстрее, чем больше параметр с. В частности, при с = 1.0 график совпадает с функцией вероятности отказа системы {М2/М/1). Таким образом, несмотря на явную зависимость вероятностей отказа систем от вида ф.р. в.б.р. и времени ремонта, она становится исчезающе малой при «быстром» восстановлении.

«Сатта^М/1 >

°0 1 2 3 4 5 ~6 7 8 " 9 10 р,

Рис. 2. Функции вероятности отказа системы

<М2ЛЛ«3/1>

.....<м2/мл> —о1=1.0 ......с2=2.0 ...........сЗ=5.0 -.....с4=10

1

°0 123456789 10 |5 = Я./Г

Рис. 3. Функции вероятности отказа п2

Аналогичным образом, когда в качестве общего распределения СИ используется распределение Вейбулла-Гнеденко, результаты численного анализа показывают, что вероятности отказа систем (М2/Ш — С/1) и (^ — С2/М/1) быстро сходится к нулю с ростом параметра /3 тем быстрее, чем больше параметр с. В

<W-G2/M/1>

.....<м2/м/1>

—С1-1.0

.............с2=2.0

..............................!............................сЗ=5 0

-.....с4=10

°0 1 23456789 10 р;

Рис. 4. Функции вероятности отказа системы {Ш — С2/М/1)

частности, при с = 1, 0 график совпадает с функцией вероятности отказа системы (М2 /М/1).

Из графиков видно также, что вероятности отказа систем {М2/Г/1) и (М2/Ш— С/1) или (Г2/М/1) и {Ш — С2/М/1) близки друг другу, что показывает незначительную чувствительность вероятности отказа этих систем к виду распределений в.б.р. и времени восстановления.

7. Заключение

В данной статье рассмотрена проблема чувствительности характеристик надёжности систем {М2/С1 /1) и {С12/М/1) к виду ф.р. в.б.р. и времени восстановления её элементов при ограничениях на доступность восстановления. Получены явные выражения для зависимости стационарных вероятностей состояний систем от вида распределений в.б.р. и времени восстановления. Тем не менее численные расчёты показывают, что эта зависимость незначительна и становится исчезающе малой при «быстром» восстановлении.

В частности, в результате численного анализа с помощью программного средства МЛТЬЛБ показано, что когда вместо общего распределения С1 используются Гамма-распределение и распределение Вейбулла-Гнеденко, вероятности отказа систем {М2/С1 /1) и {С12/М/1) быстро сходятся к нулю с ростом скорости ремонта /3 тем быстрее, чем больше параметр с. В частности, при с = 1.0 соответствующие характеристики совпадают с вероятностью отказа системы для

{М2 /М/1).

Литература

1. Севастьянов Б. А. Предельная теорема для марковских процессов и ее применение в телефонных системах с отказом // Теория вероятностей и ее применения. — 1957. — Т. 2, № 1.

2. Коваленко И. Н. Исследования по анализу сложных систем надежности. — Киев: Наукова думка, 1976. — 210 с.

3. Rykov V. V. Multidimensional Alternative Processes as Reliability Models. — Springer, 2013. — Pp. 147-157.

4. Кениг Д., Рыков В., Штойн Д. Теория массового обслуживания. — Москва, 1979. — 115 с.

PbiKOB B.B., HaH AHtHraa O HyBCTBHTeatHOCTH xapaKTepHCTHK .

77

UDC 519.718: 658.652

On Sensitivity of Systems Reliability Characteristics to the Shape of Their Elements Life and Repair Time Distributions

V. V. Rykov*, Anh Nghia Tra^

* Gubkin Russian State University of Oil and Gas ^ Peoples' Friendship University of Russia 6, Miklukho-Maklaya str., Moscow, Russian Federation, 117198

The paper deals with the problem of the systems (M2/GI/1} and (GI2/M/1} reliability characteristics sensitivity to the shape of their elements life and repair times distributions under restrictions on the availability of recovery. Partial differential equation for the time dependent and usual differential equations for the stationary micro-state probabilities of these systems are proposed. Explicit expressions for the micro- and macro-state stationary probabilities of these systems are given and they show their strong dependability on the shape of their elements life and repair times distributions. This dependence represents in terms of moment generation functions non-exponential distribution in the point of the exponential distribution parameters.

Special software tool based on the MATLAB computer system has been developed for the numerical analysis of the system failure probability sensitivity to the shape of its elements life and recovery distributions and its comparison with the simplest Markov system. The numerical analysis shows that this dependence becomes negligible and vanishes for "fast" recovery (with recovery rate increasing). In particular, it has been shown that the failure probabilities of the systems (M2/GI/1} and (GI2 /M/1} with Gamma and Weibull-Gnedenko distributions instead of the general ones quickly converge to zero with increasing recovery rate and coincide with the simplest Markov system (M2/M/1} for special value of the particular value of the parameter c = 1.0.

Key words and phrases: reliability systems, failure probabilities, sensitivity to the shape elements life and recovery time distributions, micro- and macro-states, stationary and non-stationary probabilities.

References

1. B. A. Sevast'yanov, Limit Theorem for Markov Processes and its Application to Telephone Systems with Rejection, Probability Theory and its Applications 2 (1) (1957) 363-388, in Russian.

2. I. N. Kovalenko, Investigations on Analysis of Complex Systems Reliability, Naukova Dumka, Kiev, 1976, in Russian.

3. V. V. Rykov, Multidimensional Alternative Processes as Reliability Models, Springer, 2013, pp. 147-157.

4. D. Koenig, V. Rykov, D. Schtoyn, Queueing Theory, Gubkin University Press, Moscow, 1979, in Russian.