Рассмотрим математическое описание марковского процесса с дискретными состояниями и непрерывным временем на примере графа, изображенного на рисунке 1. Будем полагать, что все переходы системы из состояния Si в Sj происходят под воздействием простейших потоков событий с интенсивностями??ij (i, j=0, 1, 2, 3); так, переход системы из состояния S0 в S1 будет происходить под воздействием потока отказов первого узла, а обратный переход из состояния S1 в S0 -- под воздействием потока "окончаний ремонтов" первого узла и т.п.

Граф состояний системы с проставленными у стрелок интенсивностями будем называть размеченным. Рассматриваемая система S имеет четыре возможных состояния: S0, S1, S2, S3.

Вероятностью i-го состояния называется вероятность pi(f) того, что в момент t система будет находиться в состоянии Si. Очевидно, что для любого момента t сумма вероятностей всех состояний равна единице:

Рассмотрим систему в момент t и, задав малый промежуток?t, найдем вероятность p0(t+?t) того, что система в момент t+?t будет находиться в состоянии S0. Это достигается разными способами.

Система в момент t с вероятностью p0(t) находилась в состоянии S0, а за время?t не вышла из него.

Вывести систему из этого состояния можно суммарным простейшим потоком с интенсивностью (л01+л02), т.е. в соответствии с формулой, с вероятностью, приближенно равной (л01+л02)?t. А вероятность того, что система не выйдет из состояния S0, равна . Вероятность того, что система будет находиться в состоянии S0 по первому способу, равна по теореме умножения вероятностей:

Система в момент t с вероятностями p1(t) (или p2(t)) находилась в состоянии S1 или S2 и за время?t перешла в состояние S0.

Потоком интенсивностью л10 система перейдет в состояние S0 с вероятностью, приближенно равной??10?t (или??20?t). Вероятность того, что система будет находиться в состоянии S0 по этому способу, равна p1(t)??10?t. Применяя теорему сложения вероятностей, получим

Переходя к пределу при?t>0 (приближенные равенства, связанные с применением формулы, перейдут в точные), получим в левой части уравнения производную (обозначим ее для простоты):

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

Рассуждая аналогично для других состояний системы S, можно получить систему дифференциальных уравнений Колмогорова для вероятностей состояний:

Сформулируем правило составления уравнений Колмогорова. В левой части каждого из них стоит производная вероятности i-го состояния. В правой части -- сумма произведений вероятностей всех состояний (из которых идут стрелки в данное состояние) на интенсивности соответствующих потоков событий, минус суммарная интенсивность всех потоков, выводящих систему из данного состояния, умноженная на вероятность данного (i-го состояния).

В системе (14) независимых уравнений на единицу меньше общего числа уравнений. Поэтому для решения системы необходимо добавить уравнение.

Нужно задать начальные условия. Так, например, систему уравнений (14) естественно решать при условии, что в начальный момент оба узла исправны и система находилась в состоянии S0, т.е. при начальных условиях p0(0)=1, p1(0)=p2(0)=p3(0)=0.

Уравнения Колмогорова дают возможность найти все вероятности состояний как функции времени. Особый интерес представляют вероятности системы рi(t) в предельном стационарном режиме, т.е. при t>?, которые называются предельными (или финальными) вероятностями состояний.

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

Предельная вероятность состояния Si имеет четкий смысл: она показывает среднее относительное время пребывания системы в этом состоянии. Например, если предельная вероятность состояния S0, т.е. p0=0.5, то это означает, что в среднем половину времени система находится в состоянии S0.

Так как предельные вероятности постоянны, то, заменяя в уравнениях Колмогорова их производные нулевыми значениями, получим систему линейных алгебраических уравнений, описывающих стационарный режим. Для системы S с графом состояний, изображенном на рисунке 1, такая система уравнений имеет вид:

Систему (15) можно составить непосредственно по размеченному графу состояний, если руководствоваться правилом, согласно которому слева в уравнениях стоит предельная вероятность данного состояния рi, умноженная на суммарную интенсивность всех потоков, ведущих из данного состояния, а справа -- сумма произведений интенсивностей всех потоков, входящих в i-e состояние, на вероятности тех состояний, из которых эти потоки исходят.

Рассмотрим математическое описание марковского процесса с дискретными состояниями и непрерывным временем на примере случайного процесса из задачи 1, граф которого изображен на рис. 1. Будем полагать, что все переходы системы из состояния S i в S j происходят под воздействием простейших потоков событий с интенсивностями l ij (i, j=0,1,2,3); так, переход системы из состояния S 0 в S 1 будет происходить под воздействием потока отказов первого узла, а обратный переход из состояния S 1 в S 0 - под воздействием потока «окончаний ремонтов» первого узла и т.п.
Граф состояний системы с проставленными у стрелок интенсивностями будем называть размеченным (см. рис. 1). Рассматриваемая система S имеет четыре возможных состояния: S 0 , S 1 , S 2 , S 3 .
Вероятностью i-го состояния называется вероятность p i (t) того, что в момент t система будет находиться в состоянии S i . Очевидно, что для любого момента t сумма вероятностей всех состояний равна единице:
. (8)
Рассмотрим систему в момент t и, задав малый промежуток Dt , найдем вероятность p 0 (t+Dt) того, что система в момент t+ Dt будет находиться в состоянии S 0 . Это достигается разными способами.
1. Система в момент t с вероятностью p 0 (t) находилась в состоянии S 0 , а за время Dt не вышла из него.
Вывести систему из этого состояния (см. граф на рис. 1) можно суммарным простейшим потоком с интенсивностью (l 01 +l 02), т.е. в соответствии с (15.7), с вероятностью, приближенно равной (l 01 +l 02)Dt . А вероятность того, что система не выйдет из состояния S 0 , равна . Вероятность того, что система будет находиться в состоянии S 0 , по первому способу (т.е. того, что находилась в состоянии S 0 и не выйдет из него за время Dt ), равна по теореме умножения вероятностей:
p 0 (t)· .
2. Система в момент t с вероятностями р 1 (t) (или p 2 (t)) находилась в состоянии S 1 или S 2 и за время Dt перешла в состояние S 0 .
Потоком интенсивностью l 10 (или l 20 - см. рис. 1) система перейдет в состояние S 0 с вероятностью, приближенно равной l 10 Dt (или l 20 Dt ). Вероятность того, что система будет находиться в состоянии S 0 по этому способу, равна р 1 (t)× l 10 Dt (или р 2 (t) × l 20 Dt ).
Применяя теорему сложения вероятностей, получим
p 0 (t+Δt)=p 1 ·λ 10 ·Δt+p 2 (t)·λ 20 ·Δt+p 0 (t),
откуда
,
Переходя к пределу при Dt ®0 (приближенные равенства, связанные с применением формулы (7), перейдут в точные), получим в левой части уравнения производную p’ 0 (t ) (обозначим ее для простоты p’ 0):
p′ 0 = λ 10 ·p 1 +λ 20 ·p 2 +(λ 10 +λ 20)·p 0 ,
Получили дифференциальное уравнение первого порядка, т.е. уравнение, содержащее как саму неизвестную функцию, так и ее производную первого порядка.
Рассуждая аналогично для других состояний системы S, можно получить систему дифференциальных уравнений Колмогорова для вероятностей состояний:
(9)
Сформулируем правило составления уравнений Колмогорова. В левой части каждого из них стоит производная вероятности i-го состояния. В правой части - сумма произведений вероятностей всех состояний (из которых идут стрелки в данное состояние) на интенсивности соответствующих потоков событий, минус суммарная интенсивность всех потоков, выводящих систему из данного состояния, умноженная на вероятность данного (i-го состояния).
В системе (9) независимых уравнений на единицу меньше общего числа уравнений. Поэтому для решения системы необходимо добавить уравнение (8).
Особенность решения дифференциальных уравнений вообще состоит в том, что требуется задать так называемые начальные условия, т.е. в данном, случае вероятности состояний системы в начальный момент t = 0. Так, например, систему уравнений (9) естественно решать при условии, что в начальный момент оба узла исправны и система находилась в состоянии S 0 , т.е. при начальных условиях p 0 (0)=1, p 1 (0)=p 2 (0)=p 3 (0)=0.
Уравнения Колмогорова дают возможность найти все вероятности состояний как функции времени . Особый интерес представляют вероятности системы p i (t ) в предельном стационарном режиме, т.е. при t→∞, которые называются предельными (или финальными) вероятностями состояний.
В теории случайных процессов доказывается, что если число состояний системы конечно и из каждого из них можно (за конечное число шагов) перейти в любое другое состояние, то предельные вероятности существуют.
Предельная вероятность состояния S i имеет четкий смысл: она показывает среднее относительное время пребывания системы в этом состоянии. Например, если предельная вероятность состояния S 0 , т.е. p 0 =0,5, то это означает, что в среднем половину времени система находится в состоянии S 0 .
Так как предельные вероятности постоянны, то, заменяя в уравнениях Колмогорова их производные нулевыми значениями, получим систему линейных алгебраических уравнений, описывающих стационарный режим. Для системы S с графом состояний, изображенном на рис. 1, такая система уравнений имеет вид:
(10)
Систему (10) можно составить непосредственно по размеченному графу состояний, если руководствоваться правилом, согласно которому слева в уравнениях стоит предельная вероятность данного состояния p i , умноженная на суммарную интенсивность всех потоков, ведущих из данного состояния, а справа - сумма произведений интенсивностей всех потоков, входящих в i-е состояние, на вероятности тех состояний, из которых эти потоки исходят.

Задача 2. Найти предельные вероятности для системы S задачи 1, граф состояний которой приведен на рис. 1, при l 01 =1, l 02 =2, l 10 =2, l 13 =2, l 20 =3, l 23 =1, l 31 =3, l 32 =2.
Решение . Система алгебраических уравнений , описывающих стационарный режим для данной системы, имеет вид (10) или
3p 0 =2p 1 +3p 2 (11)
4p 1 =p 0 +3p 3
4p 2 =2p 0 +2p 3
p 0 +p 1 +p 2 +p 3 =1
(Здесь мы вместо одного "лишнего" уравнения системы (10) записали нормировочное условие (8)).
Решив систему (11), получим p 0 =0,40, p 1 =0,20, p 2 =0,27, p 3 =0,13, т.е. в предельном, стационарном режиме система S в среднем 40% времени будет находиться в состоянии S 0 (оба узла исправны), 20% - в состоянии S 1 (первый узел ремонтируется, второй работает), 27% - в состоянии S 2 (второй узел ремонтируется, первый работает) и 13% времени - в состоянии S 3 (оба узла ремонтируются).

Задача 3. Найти средний чистый доход от эксплуатации в стационарном режиме системы S в условиях задач 1 и 2, если известно, что в единицу времени исправная работа первого и второго узлов приносит доход соответственно в 10 и 6 ден.ед., а их ремонт требует затрат соответственно в 4 и 2 ден.ед. Оценить экономическую эффективность СМО имеющейся возможности уменьшения вдвое среднего времени ремонта каждого из двух узлов, если при этом придется вдвое увеличить затраты на ремонт каждого узла (в единицу времени).
Решение. Из задачи 2 следует, что в среднем первый узел исправно работает долю времени, равную p 0 +p 3 =0,40+0,27=0,67, а второй узел - p 0 +p 1 =0,40+0,20=0,60. В то же время первый узел находится в ремонте в среднем долю времени, равную p 1 +p 3 =0,20+0,13=0,33, а второй узел – p 2 +p 3 =0,27+0,13=0,40. Поэтому средний чистый доход в единицу времени от эксплуатации системы, т.е. разность между доходами и затратами, равен
Д=0,67 ×10+0,60×6-0,33 ×4-0,40×2=8,18 ден.ед.
Уменьшение вдвое среднего времени ремонта каждого из узлов в соответствии с (6) будет означать увеличение вдвое интенсивностей потока "окончаний ремонтов" каждого узла, т.е. теперь l 10 =4, l 20 =6, l 31 =6, l 32 =4 и система линейных алгебраических уравнений (10), описывающая стационарный режим системы вместе с нормировочным условием (8) примет вид:
3p 0 =4p 1 +6p 2
6p 1 =p 0 +6p 3
7p 2 =2p 0 +4p 3
p 0 +p 1 +p 2 +p 3 =1
Решив систему, получим p 0 =0,60, p 1 =0,15, p 2 =0,20, p 3 =0,05.
Учитывая, что p 0 +p 2 =0,60+0,20=0,80, p 0 +p 1 =0,60+0,15=0,75, p 1 +p 3 =0,15+0,05=0,20, p 2 +p 3 =0,20+0,05=0,25, а затраты на ремонт первого и второго узла составляют теперь соответственно 8 и 4 ден. ед., вычислим средний чистый доход в единицу времени:Д 1 =0,80 ×10+0,75×6-0,20 ×8-0,25×4=9,9 ден.ед.
Так как Д 1 больше Д (примерно на 20%), то экономическая целесообразность ускорения ремонтов узлов очевидна.

Пример . Техническое устройство может находиться в одном из трех состояний S 0 , S 1 , S 2 . Интенсивность потоков, переводящих устройство из состояния, заданы в таблице.

Необходимо построить размеченный граф состояний, записать систему уравнений Колмогорова, найти финальные вероятности и сделать анализ полученных решений.
Размеченный граф состояний имеет вид.

p 0 (t) + p 1 (t) + p 2 (t) = 1

p 0 (t) + p 1 (t) + p 2 (t) = 1

2p 0 -3p 1 = 0
2p 0 +2p 1 -3p 2 =0
p 0 + p 1 + p 2 = 1
Решим СЛАУ с помощью метода Гаусса.
Вывод: При достаточно большом времени работы техническое устройство с вероятностью p 0 = 0.36 будет находиться в состоянии S 0 , с вероятностью p 1 = 0.24 в состоянии S 1 и с вероятностью p 2 = 0.4 в состоянии S 2 .

Пример .
Техническое устройство может находиться в одном из трех состояний S 0 , S 1 , S 2 . Интенсивность потоков, которые переводят устройства из одного состояния во второе, известны λ 01 =2, λ 10 =4, λ 21 =2, λ 12 =3, λ 20 =4.
Необходимо построить размеченный граф состояний, записать систему уравнений Колмогорова, найти финальные вероятности и сделать анализ полученных решений.
Размеченный граф состояний имеет вид.

По графу запишем систему уравнений Колмогорова в общем виде:

Вместо интенсивности потоков λ ij запишем их конкретные значения и получим искомую систему:

Чтобы найти финальные вероятности состояний, в уравнениях Колмогорова отбросим первое уравнения, а по остальным составим систему алгебраических уравнений:
2p 0 -7p 1 +2p 2 =0
3p 1 -6p 2 =0
p 0 +p 1 +p 2 =1
Делим первое уравнение на 2, а второе на 3 и получим систему
p 0 -7p 1 +2p 2 =0
3p 1 -6p 2 =0
p 0 +p 1 +p 2 =1
Из третьего уравнения вычитаем первое
p 0 -3.5p 1 +p 2 =0
p 1 -2p 2 =0
4.5p 1 =1
Отсюда получим p 1 =0,22, p 2 =0,11 и p 0 =0,67.
Вывод: При достаточно большом времени работы техническое устройство с вероятностью p 0 = 0,67 будет находиться в состоянии S 0 , с вероятностью p 1 = 0,22 в состоянии S 1 и с вероятностью p 2 = 0,11 в состоянии S 2 .

Процесс гибели и размножения

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

Рис. 4
Рассмотрим упорядоченное множество состояний системы S 0 , S 1 , S 2 , …, S k . Переходы могут осуществляться из любого состояния только в состояния с соседними номерами, т.е. из состояния S k возможны переходы только либо в состояние S k-1 , либо в состояние S k+1 . (При анализе численности популяций считают, что состояние S k соответствует численности популяции, равной k, и переход системы из состояния S k в состояние S k+1 происходит при рождении одного члена популяции, а переход в состояние S k-1 , - при гибели одного члена популяции).
Предположим, что все потоки событий, переводящие систему по стрелкам графа, простейшие с соответствующими интенсивностями l k, k+1 или l k+1, k .
По графу, представленному на рис. 4, составим и решим алгебраические уравнения для предельных вероятностей состояний (их существование вытекает из возможности перехода из каждого состояния в каждое другое и конечности числа состояний).
В соответствии с правилом составления таких уравнений (см. 13) получим: для состояния S 0
λ 01 p 0 = λ 10 p 1 (12)
для состояния S 1 – (l 12 +l 10)p 1 =l 01 p 0 +l 21 p 2 , которое с учетом (12) приводится к виду
λ 12 p 1 = λ 21 p 2 (13)
Аналогично, записывая уравнения для предельных вероятностей других состояний, можно получить следующую систему уравнений:
(14)
к которой добавляется нормировочное условие
p 0 +p 1 +p 2 +...+p n =1 (15)
Решая систему (14), (15), можно получить (16)
(17)
Легко заметить, что в формулах (17) для p 1 , p 2 , …, p n коэффициенты при p 0 есть слагаемые, стоящие после единицы в формуле (16). Числители этих коэффициентов представляют произведение всех интенсивностей, стоящих у стрелок, ведущих слева направо до данного состояния S k (k=1, 2, …, n), а знаменатели - произведение всех интенсивностей, стоящих у стрелок, ведущих справа налево до состояния S k .

Задача 4. Процесс гибели и размножения представлен графом (рис. 5). Найти предельные вероятности состояний.

Рис. 5

Решение. По формуле (16) найдем

по (17) – т.е. в установившемся, стационарном режиме в среднем 70,6% времени система будет находиться в состоянии S 0 , 17,6% - в состоянии S 1 и 11,8% - в состоянии S 2 .

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

Решение:

Система может находиться в четырех состояниях, так как у каждого аппарата по продаже билетов есть два состояния (быть занятым или свободным). Пусть S 0 - оба аппарата заняты; S 1 - 1-ый занят, 2-ой свободен; S 2 - 1-ый свободен, 2-ой занят; S 3 - оба аппарата свободны. Построим граф состояний, отметив на нем все возможные состояния кругами, а возможные переходы из состояния в состояние обозначим стрелками. Получаем, что переход из S 0 в S 3 возможен либо через S 1 , либо через S 2 , либо напрямик, как показано на рисунке 4.

Рисунок 4 - Граф состояний аппаратов по продаже билетов

Найти предельные вероятности для системы S, граф которой изображен на рисунке.

Решение:

В теории случайных процессов доказывается, что если число состояний системы конечно и из каждого из них можно (за конечное число шагов) перейти в любое другое состояние, то предельные вероятности существуют. Их можно найти из уравнений Колмогорова, составив систему по данному размеченному графу состояний, по следующему правилу:

Слева в уравнении стоит предельная вероятность данного состояния p i , умноженная на суммарную интенсивность всех потоков, ведущих из данного состояния, а справа - сумма произведений интенсивностей всех потоков, входящих в данное состояние, на вероятности тех состояний, из которых эти состояния выходят.

Кроме этого надо учитывать, что сумма всех вероятностей данной конечной системы равна единице. Составим уравнения для состояний S 1 и S 2 (уравнение для состояния S 0 - «лишнее»):

Ответ: Система примерно 66,67% времени пребывает в состоянии S 0 , 25% - в состоянии S 1 и 8,33% времени находится в состоянии S 2 .

Найти валовой выпуск для сбалансированной многоотраслевой экономики в модели Леонтьева, если дана матрица прямых затрат А и вектор конечного потребления У:

Решение:

Для сбалансированной многоотраслевой экономики выполняется следующее соотношение:

Выразим валовой выпуск через конечное потребление и матрицу затрат:

Находим матрицу, обратную к (Е - А):

Найдем валовой выпуск:

Ответ: Валовой выпуск равен (811,3; 660,4).

*При решении задач использовался

7 ГЛАВА

Элементы теории случайных процессов и теории массового обслуживания

В главе рассматриваются:

определение случайного процесса и его характеристики, понятие марковского случайного процесса;

основные понятия теории массового обслуживания;

потоки событий;

уравнение Колмогорова;

СМО с отказами;

метод Монте-Карло.

Типовые задачи

Пример 7.1

Случайный процесс определяется формулой X (t ) = Xcoswt , где X – случайная величина. Найти основные характеристики этого процесса, если М(Х) = а, D (X ) = а 2 .

Решение

На основании свойств математического ожидания и дисперсии имеем:

Корреляционную функцию найдем по формуле (7.1):

Нормированную корреляционную функцию найдем по формуле (7.2):

Пример 7.2

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

Решение

Возможные состояния системы: S 0 – оба узла исправны; S 1 – первый узел ремонтируется, второй исправен; S 2 – второй узел ремонтируется, первый исправен; S 3 - оба узла ремонтируются. Граф системы приведен на рис. 7.4.

Стрелка, направленная, например, из S 0 в S 1 , означает переход системы в момент отказа первого узла, из S 1 в S 0 – переход в момент окончания ремонта этого узла.

На графе отсутствуют стрелки из S 0 в S 3 и из S 1 в S 2 . Это объясняется тем, что выходы узлов из строя предполагаются независимыми друг от друга и, например, вероятностью одновременного выхода из строя двух узлов (переход из S 0 в S 3 ) или одновременного окончания ремонтов двух узлов (переход из S 3 в S 0 ) можно пренебречь.

Пример 7.3

На автоматическую телефонную станцию поступает простейший поток вызовов с интенсивностью альфа = 1,2 вызовов в минуту. Найти вероятность того, что за две минуты:

а) не придет ни одного вызова;

б) придет ровно один вызов;

в) придет хотя бы один вызов.

Решение

а) Случайная величина X – число вызовов за две минуты – распределена по закону Пуассона с параметром λτ = 1,2*2 = 2,4. Вероятность того, что вызовов не будет (m = 0), по формуле (7.5):

б) Вероятность одного вызова (m = 1):

в) Вероятность хотя бы одного вызова:

Пример 7.4

Найти предельные вероятности для системы S из примера 7.2, граф состояний которой приведен на рис. 7.4, при λ 01 = 1, λ 02 = 2, λ 10 = 2, λ 13 = 2, λ 20 = 3, λ 23 = 1, λ 31 = 3, λ 32 = 2.

Решение

Система алгебраических уравнений, описывающих стационарный режим для данной системы, имеет вид (7.14) или

(7.15)

(Здесь вместо одного «лишнего»уравнения системы (7.14) записали нормировочное условие(7.12).)

Решив систему (7.15), получим p 0 = 0,40, p 1 = 0,20, p 2 = 0,27, p 3 = 0,13, т.е. в предельном стационарном режиме системе S в среднем 40% времени будет находиться в состоянии S 0 (оба узла исправны), 20% - в состоянии S 1 (первый узел ремонтируется, второй работает), 27% - в состоянии S 2 (второй узел ремонтируется, первый работает) и 13% времени – в состоянии S 3 (оба узла ремонтируются).

Пример 7.5

Найти средний чистый доход от эксплуатации в стационарном режиме системы S в условиях примеров 7.2 и 7.4, если известно, что в единицу времени исправная работа первого и второго узлов приносит доход соответственно в 10 и 6 ден. ед., а их ремонт требует затрат соответственно в 4 и 2 ден. ед. Оценить экономическую эффективность имеющейся возможности уменьшения вдвое среднего времени ремонта каждого из двух узлов, если при этом придется вдвое увеличить затраты на ремонт каждого узла (в единицу времени).

Решение

Из примера 7.4 следует, что в среднем первый узел исправно работает долю времени, равную p 0 + р 2 = 0,40 + 0,27 = 0,67, а второй узел – p 0 + р 1 = 0,40 + 0,20 = 0,60. в то же время первый узел находится в ремонте в среднем долю времени, равную р 1 + р 3 = 0,20 + 0,13 = 0,33, а второй узел – р 2 +р 3 = =0,27+0,13=0,40. Поэтому средний чистый доход в единицу времени от эксплуатации системы, т.е. разность между доходами и затратами, равен

D = 0,67*10+0,60*6-0,33*4-0,40*2 = 8,18 ден. ед.

Уменьшение вдвое среднего времени ремонта каждого из узлов в соответствии с (7.10) будет означать увеличение вдвое интенсивностей потока «окончаний ремонтов» каждого узла, т.е. теперь λ 10 = 4, λ 20 = 6, λ 31 = 6, λ 32 = 4 и система линейных алгебраических уравнений (7.14), описывающая стационарный режим системы S , вместе с нормировочным условием (7.12) примет вид:

Решив систему, получим р 0 = 0,60, р 1 = 0,15, р 2 = 0,20, р 3 = 0,05.

Учитывая, что р 0 + р 2 = 0,60 + 0,20 = 0,80, р 0 + р 1 = 0,60 + 0,15 = 0,75, р 1 + р 3 = 0,15 + 0,05 = 0,20, р 2 + р 3 = 0,20 + 0,05 = 0,25, а затраты на ремонт первого и второго узлов составляют теперь соответственно 8 и 4 ден. ед., вычислим средний чистый доход в единицу времени:

D 1 = 0,80*10+0,75*6-0,20*8-0,25*4 = 9,9 ден. ед.

Так как D 1 больше D (примерно на 20%), то экономическая целесообразность ускорения ремонтов узлов очевидна.

Пример 7.6

Процесс гибели и размножения представлен графом (рис.7.8). Найти предельные вероятности состояний.

Решение

По формуле (7.20) найдем

,

,
,

т.е. в установившемся стационарном режиме в среднем 70,6% времени система будет находится в состоянии S 0 , 17,6% - в состоянии S 1 и 11,8% - в состоянии S 2 .

Пример 7.7

Известно, что заявки на телефонные переговоры в телевизионном ателье поступают с интенсивностью Я, равной 90 заявок в час, а средняя продолжительность разговора по телефону t o 6 – 2 мин. Определить показатели эффективности работы СМО (телефонной связи) при наличии одного телефонного номера.

Решения

Имеем λ = 90 (1/ч), t об = 2 мин. Интенсивность потока обслуживаний μ = 1/ t об = 1/2 = 0,5 (1/мин) = 30 (1/ч). По (7.24) относительная пропускная способность СМО Q = 30/(90+30) = 0,25, т.е. в среднем только 25% поступающих заявок осуществят переговоры по телефону. Соответственно вероятность отказа в обслуживании составит P ОТК =0,75 (см. (7.25)). Абсолютная пропускная способность СМО по (7.26) A=90*0,25=22,5, т.е. в среднем в час будут обслужены 22,5 заявки на переговоры. Очевидно, что при наличии только одного телефонного номера

СМО будет плохо справляться с потоком заявок.

Пример 7.8

В условиях примера 7.7 определить оптимальное число телефонных номеров в телевизионном ателье, если условием оптимальности считать удовлетворение из каждых 100 заявок на переговоры в среднем не менее 90 заявок.

Решение

Интенсивность нагрузки канала по формуле (7.28) р =90/30=3, т.е. за время среднего (по продолжительности) телефонного разговора tоб = 2 мин поступает в среднем 3 заявки на переговоры.

Будем постепенно увеличивать число каналов (телефонных номеров) п = 2, 3, 4,... и определим по формулам (7.29), (7.32), (7.33) для получаемой n-канальной СМО характеристики обслуживания. Например, при п = 2 p o =(1 + 3 + 3 2 /2!) -1 = 0,118 ≈ 0,12;

Q = 1-(3 2 /2!)*0,118 ≈ 0,471; А = 90*0,471 = 42,4.

Значение характеристик СМО сведем в табл. 7.1.

Таблица 7.1

Характеристика обслуживания	Обозначение	Число каналов (телефонных номеров)
Относительная пропускная способность
Абсолютная пропускная способность

По условию оптимальности Q ≥ 0,9, следовательно, в телевизионном ателье необходимо установить 5 телефонных номеров (в этом случае Q = 0,90 – см. табл. 7.1). При этом в час будут обслуживаться в среднем 80 заявок (А = 80,1), а среднее число занятых телефонных номеров (каналов) по формуле (7.34) = 80, 1/30 = 2,67.

Пример 7.9

В вычислительный центр коллективного пользования с тремя ЭВМ поступают заказы от предприятий на вычислительные работы. Если работают все три ЭВМ, то вновь поступающий заказ не принимается, и предприятие вынуждено обратиться в другой вычислительный центр. Среднее время работы с одним заказом составляет 3 ч. Интенсивность потока заявок 0,25 1/ч. Найти предельные вероятности состояний и показатели эффективности работы вычислительного центра.

Решение

По условию п=3, λ =0,25 (1/ч), t o 6 =3 (ч). Интенсивность потока обслуживаний μ=1/ t o 6 =1/3=0,33. Интенсивность нагрузки ЭВМ по формуле (7.28) р =0,25/0,33=0,75. Найдем предельные вероятности состояний:

по формуле (7.29):

p 0 =(1+0,75+0,75 2 /2!+0,75 3 /3!) -1 =0,476;

по формуле (7.30):

p 1 = 0,75*0,476 = 0,357; p 2 = (0,75 2 /2!)*0,476 = 0,134; p 3 = (0,75 3 /3!)*0,476 = 0,033,

т.е. в стационарном режиме работы вычислительного центра в среднем 47,6% времени нет ни одной заявки, 35,7% – имеется одна заявка (занята одна ЭВМ), 13,4% – две заявки (две ЭВМ), 3,3% времени – три заявки (заняты три ЭВМ).

Вероятность отказа (когда заняты все три ЭВМ), таким образом, Р 0ТК = р 3 = 0,033.

Согласно формуле (7.32) относительная пропускная способность центра Q = 1 – 0,033 = 0,967, т.е. в среднем из каждых 100 заявок вычислительный центр обслуживает 96,7 заявок.

По формуле (7.33) абсолютная пропускная способность центра А = 0,250,967 = 0,242, т.е. в один час в среднем обслуживается 0,242 заявки.

Согласно формуле (7.34) среднее число занятых ЭВМ к = = 0,242/0,33 = 0,725, т.е. каждая из трех ЭВМ будет занята обслуживанием заявок в среднем лишь на 72,5/3 = 24,2%.

При оценке эффективности работы вычислительного центра необходимо сопоставить доходы от выполнения заявок с потерями от простоя дорогостоящих ЭВМ (с одной стороны, здесь высокая пропускная способность СМО, а с другой – значительный простой каналов обслуживания) и выбрать компромиссное решение.

Случайный процесс определяется формулой X(t)= Хе(- t ) (t > 0), где X – случайная величина, распределенная по нормальному закону с параметрами а и а 2 . Найти математическое ожидание, дисперсию, корреляционную и нормированную корреляционную функции случайного процесса.

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

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

Среднее число заказов на такси, поступающих на диспетчерский пункт в одну минуту, равно 3. Найти вероятность того, что за две минуты поступит:

а) 4 вызова;

б) хотя бы один;

в) ни одного вызова.

(Поток заявок простейший.)

Найти предельные вероятности для систем S , граф которых изображен на рис. 7.11 и 7.12.

Рассматривается круглосуточная работа пункта проведения профилактического осмотра автомашин с одним каналом (одной группой проведения осмотра). На осмотр и выявление дефектов каждой машины затрачивается в среднем 0,5 ч. На осмотр поступает в среднем 36 машин в сутки. Потоки заявок и обслуживаний – простейшие. Если машина, прибывшая в пункт осмотра, не застает ни одного канала свободным, она покидает пункт осмотра необслуженной. Определить предельные вероятности состояний и характеристики обслуживания профилактического пункта осмотра.

Решить задачу 7.15 для случая п = 4 канала (групп проведения осмотра). Найти минимальное число каналов, при котором относительная пропускная способность пункта осмотра будет не менее 0,9.

Одноканальная СМО с отказами представляет собой одну телефонную линию, на вход которой поступает простейший поток вызовов с интенсивностью 0,4 вызовов/мин. Средняя продолжительность разговора 3 мин.; время разговора имеет показательное распределение. Найти предельные вероятности состояний и характеристики обслуживания СМО. Сравнить пропускную способность СМО с номинальной, которая была бы, если разговор длился в точности 3 мин., а заявки шли одна за другой регулярно, без перерывов.

Имеется двухканальная простейшая СМО с отказами. На ее вход поступает поток заявок с интенсивностью 4 заявки/ч. Среднее время обслуживания одной заявки 0,8 ч. Каждая обслуженная заявка приносит доход 4 ден. ед. Содержание каждого канала обходится 2 ден. ед./ч. Выяснить, выгодно или невыгодно в экономическом отношении увеличить число каналов до трех.

Задания по вариантам

Варианты	Главы

Что будет происходить с вероятностями состояний при .Будут лиP 1 (t), P 2 (t), … стремится к каким-то пределам? Если эти пределы существуют и не зависят от начального состояния системы, то они называются финальными вероятностями состояний. В теории случайных процессов доказывается, что если число n состояний системы конечно и из каждого из них можно (за конечное число шагов) перейти в любое другое, то финальные вероятности существуют (это условие достаточно, но не необходимо для существования финальных вероятностей).

Предположим, что это условие выполнено и финальные вероятности существуют:

Будем обозначать их теми же буквами P 1 , P 2 , …, что и сами вероятности состояний, но подразумевая под ними не функции времени, а постоянные числа. Очевидно, они тоже образуют в сумме единицу:

. (4.10)

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

Например, если система S имеет три состояния S 1 , S 2 , S 3 и их финальные вероятности равны 0,2; 0,3; 0,5, это значит, что в предельном стационарном режиме система в среднем две десятых времени проводит в состоянии S 1 , три десятых – в состоянии S 2 и половину времени – в состоянии S 3 .

Как же вычислить финальные вероятности? Если вероятности P 1 , P 2 , … постоянны, то их производные равны нулю. Значит, чтобы найти финальные вероятности, нужно все левые части в уравнениях Колмогорова положить равными нулю и решить полученную систему уже не дифференциальных, а линейных алгебраических уравнений. Даже можно сразу по графу состояний написать систему алгебраических уравнений. Если перенести отрицательный член каждого уравнения из правой части в левую, то получим сразу систему уравнений, где слева стоит финальная вероятность данного состояния P i , умноженная на суммарную интенсивность всех потоков, ведущих из данного состояния, а справа – сумма произведений интенсивностей всех потоков, входящих в i – е состояние, на вероятности тех состояний, из которых эти потоки исходят.

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

(4.11)

Эту систему 4-х уравнений с 4-мя неизвестными P 0 , P 1 , P 2 , P 3 можно решить воспользовавшись так называемым нормировочным условием:

, (4.12)

при этом одно (любое) из уравнений можно отбросить (оно вытекает как следствие из остальных).

Зададимся численными значениями интенсивностей λ 1 =1, λ 2 =2, μ 1 =2, μ 2 =3 и решим систему (4.11). Отбросим четвертое уравнение, добавив вместо него нормировочное условие (4.12). Уравнения примут вид:

(4.13)

Решая их, получим т.е. в предельном, стационарном режиме системаS в среднем 40% времени будет проводить в состоянии S 0 (оба узла исправны), 20% - в состоянии S 1 (первый узел ремонтируется, второй работает), 27% - в состоянии S 2 (второй узел ремонтируется, первый работает) и 13% - в состоянии S 3 полной негодности (оба узла ремонтируются). Знание этих предельных вероятностей может помочь оценить среднюю эффективность работы системы и загрузку ремонтных органов. Предположим, что система S в состоянии S 0 приносит в единицу времени доход 8 (условных единиц), в состоянии S 1 – доход 3, в состоянии S 2 – доход 5, а в состоянии S 3 – вообще не приносит дохода. Тогда, в предельном стационарном режиме средний доход в единицу времени будет . Теперь оценим загрузку ремонтных органов (рабочих), занятых ремонтов узлов 1 и 2. Узел 1 ремонтируется долю времени, равную. Узел 2 ремонтируется долю времени
.

Здесь уже может возникнуть вопрос об оптимизации решения. Допустим, что мы можем уменьшить среднее время ремонта того или другого узла (а может быть и того, и другого), но это нам обойдется в какую-то сумму. И необходимо оценить, а окупит ли увеличение дохода, связанное с ускорением ремонта, повышенные расходы на ремонт? (для этого надо будет решить систему 4-х уравнений с 4-мя неизвестными).