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

Уравнения Колмогорова.
Предельные вероятности состояний

Рассмотрим математическое описание марковского процесса с дискретными состояниями и непрерывным временем* на примере случайного процесса из примера 1, граф которого изображен на рис. 1. Будем полагать, что все переходы системы из состояния S_i в S_j происходят под воздействием простейших потоков событий с интенсивностями lambda_{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 сумма вероятностей всех состояний равна единице:

sum_{i=0}^{3}p_i(t)=p_0(y)+p_1(t)+p_2(t)+p_3(t)=1.

(8)

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

1. Система в момент t с вероятностью p_0(t) находилась в состоянии S_0, а за время Delta t не вышла из него.

Вывести систему из этого состояния (см. граф на рис. 1) можно суммарным простейшим потоком с интенсивностью (lambda_{01}+lambda_{02}), т.е. в соответствии с формулой (7), с вероятностью, приближенно равной (lambda_{01}+lambda_{02})Delta t. А вероятность того, что система не выйдет из состояния S_0, равна [1-(lambda_{01}+lambda_{02})Delta t]. Вероятность того, что система будет находиться в состоянии S_0 по первому способу (т.е. того, что находилась в состоянии S_0 и не выйдет из него за время Delta t), равна по теореме умножения вероятностей:

p_0(t)[1-(lambda_{01}+lambda_{02})Delta t].

2. Система в момент t с вероятностями p_1(t) (или p_2(t)) находилась в состоянии S_1 или S_2 и за время Delta t перешла в состояние S_0.

Потоком интенсивностью lambda_{10} (или lambda_{20} — с- рис. 1) система перейдет в состояние S_0 с вероятностью, приближенно равной lambda_{10}Delta t (или lambda_{20}Delta t). Вероятность того, что система будет находиться в состоянии S_0 по этому способу, равна p_1(t)lambda_{10}Delta t (или p_2(t)lambda_{20}Delta t).

Применяя теорему сложения вероятностей, получим

p_0(t+Delta t)= p_1(t)lambda_{10}Delta t+ p_2(t)lambda_{20}Delta t+p_0(t)[1- (lambda_{01}+lambda_{02})Delta t],

откуда

frac{p_0(t+Delta t)-p_0(t)}{Delta t}= p_1(t)lambda_{10}+ p_2(t)lambda_{20}-(lambda_{01}+lambda_{02})p_0(t).

Переходя к пределу при Delta tto0 (приближенные равенства, связанные с применением формулы (7), перейдут в точные), получим в левой части уравнения производную p'_0(t) (обозначим ее для простоты p'_0):

p'_0=lambda_{10}p_1+lambda_{20}p_2-(lambda_{01}+lambda_{02})p_0.

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

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

begin{cases} p'_0=lambda_{10}p_1+lambda_{20}p_2-(lambda_{01}+lambda_{02})p_0,\[2pt] p'_1=lambda_{01}p_0+lambda_{31}p_3-(lambda_{10}+lambda_{13})p_1,\[2pt] p'_2=lambda_{02}p_0+lambda_{32}p_3-(lambda_{20}+lambda_{23})p_2,\[2pt] p'_3=lambda_{13}p_1+lambda_{23}p_2-(lambda_{31}+lambda_{32})p_3. end{cases}

(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) в предельном стационарном режиме, т.е. при ttoinfty, которые называются предельными (или финальными) вероятностями состояний.

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

Предельная вероятность состояния S_i имеет четкий смысл: она показывает среднее относительное время пребывания системы в этом состоянии. Например, если предельная вероятность состояния S_0, т.е. p_0=0,!5, то это означает, что в среднем половину времени система находится в состоянии S_0.

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

begin{cases}(lambda_{01}+lambda_{02})p_0= lambda_{10}p_1+lambda_{20}p_2,\[2pt] (lambda_{10}+lambda_{13})p_1= lambda_{01}p_0+lambda_{31}p_3,\[2pt] (lambda_{20}+lambda_{23})p_2= lambda_{02}p_0+lambda_{32}p_3,\[2pt] (lambda_{31}+lambda_{32})p_3= lambda_{13}p_1+lambda_{23}p_2. end{cases}

(10)

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


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

lambda_{01}=1,quad lambda_{02}=2,quad lambda_{10}=2,quad lambda_{13}=2,quad lambda_{20}=3,quad lambda_{23}=1,quad lambda_{31}=3,quad lambda_{32}=2.

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

begin{cases}3p_0=2p_1+3p_2,\ 4p_1=p_0+3p_3,\ 4p_2=2p_0+2p_3,\ p_0+p_1+p_2+p_3=1.end{cases}

(11)

(Здесь мы вместо одного “лишнего” уравнения системы (10) записали нормировочное условие (8)).

Решив систему (11), получим p_0=0,!4,~p_1=0,,2,~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,!4+0,!27=0,!67, а второй узел — p_0+p_1=0,!4+0,!2=0,!6. В то же время первый узел находится в ремонте в среднем долю времени, равную p_1+p_3=0,!2+0,!13=0,!33, а второй узел — p_2+p_3=0,!27+0,!13=0,!4. Поэтому средний чистый доход overline{D} в единицу времени от эксплуатации системы, т.е. разность между доходами и затратами, равен

overline{D}= 0,!67cdot10+0,!6cdot6-0,!33cdot4-0,!4cdot2= 8,!18 ден. ед.

Уменьшение вдвое среднего времени ремонта каждого из узлов в соответствии с (6) будет означать увеличение вдвое интенсивностей потока “окончаний ремонтов” каждого узла, т.е. теперь lambda_{10}=4, lambda_{20}=6, lambda_{31}=6, lambda_{32}=4 и система линейных алгебраических уравнений (10), описывающая стационарный режим системы S, вместе с нормировочным условием (8) примет вид:

begin{cases}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.end{cases}

Решив систему, получим p_0=0,!6,~p_1=0,!15,~p_2=0,!2,~p_3=0,!05.

Учитывая, что p_0+p_2=0,!8,~p_0+p_1=0,!75,~p_1+p_3=0,!2,~p_2+p_3=0,!25, а затраты на ремонт первого и второго узла составляют теперь соответственно 8 и 4 ден.ед., вычислим средний чистый доход overline{D}_1 в единицу времени:

overline{D}_1=0,!8cdot10+0,!75cdot6-0,!2cdot8-0,!25cdot4=9,!9 ден.ед.

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


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

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

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

Граф состояний процесса гибели и размножения

Рассмотрим упорядоченное множество состояний системы S_0,S_1,S_2,ldots,S_k. Переходы могут осуществляться из любого состояния только в состояния с соседними номерами, т.е. из состояния S_k возможны переходы только либо в состояние S_{k-1}, либо в состояние S_{k+1}.

Предположим, что все потоки событий, переводящие систему по стрелкам графа, простейшие с соответствующими интенсивностями lambda_{k,k+1} или lambda_{k+1,k}.

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

В соответствии с правилом составления таких уравнений (см. 13) получим: для состояния S_0

lambda_{01}cdot p_0=lambda_{10}cdot p_1,

(12)

для состояния S_1 имеем (lambda_{12}+lambda_{10})p_1=lambda_{01}p_0+ lambda_{21}p_2, которое с учетом (12) приводится к виду

lambda_{12}cdot p_1=lambda_{21}cdot p_2.

(13)

Аналогично, записывая уравнения для предельных вероятностей других состояний, можно получить следующую систему уравнений:

begin{cases}lambda_{01}cdot p_0=lambda_{10}cdot p_1,\ lambda_{12}cdot p_1=lambda_{21}cdot p_2,\ cdotscdotscdotscdotscdots\ lambda_{k-1,k}cdot p_{k-1}=lambda_{k,k+1}cdot p_k,\ cdotscdotscdotscdotscdots\ lambda_{n-1,n}cdot p_{n-1}=lambda_{n,n+1}cdot p_n.end{cases}

(14)

к которой добавляется нормировочное условие

p_0+p_1+p_2+ldots+p_n=1.

(15)

При анализе численности популяций считают, что состояние S_k соответствует численности популяции, равной k, и переход системы из состояния S_k в состояние S_{k+1} происходит при рождении одного члена популяции, а переход в состояние S_{k-1} – при гибели одного члена популяции.

Решая систему (14), (15), можно получить

p_0={left(1+ frac{lambda_{01}}{lambda_{10}}+ frac{lambda_{12}lambda_{01}}{lambda_{21}lambda_{10}}+ldots+ frac{lambda_{n-1,n}cdotslambda_{12}lambda_{01}}{lambda{n,n-1}cdotslambda_{21}lambda_{10}}right)!}^{-1},

(16)

p_1=frac{lambda_{01}}{lambda_{10}}p_0,quad p_2=frac{lambda_{12}lambda_{01}}{lambda_{21}lambda_{10}}p_0,quad ,ldots,quad p_n=frac{lambda_{n-1,n}cdotslambda_{12}lambda_{01}}{lambda{n,n-1}cdotslambda_{21}lambda_{10}}p_0.

(17)

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

Легко заметить, что в формулах (17) для p_1,p_2,ldots,p_n коэффициенты при p_0 есть слагаемые, стоящие после единицы в формуле (16). Числители этих коэффициентов представляют произведение всех интенсивностей, стоящих у стрелок, ведущих слева направо до данного состояния S_k,(k=1,2,ldots,n), а знаменатели — произведение всех интенсивностей, стоящих у стрелок, ведущих справа налево до состояния S_k.


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

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

p_0={left(1+frac{1}{4}+frac{2cdot1}{3cdot4}right)!}^{-1}= 0,!706, по (17) p_1=frac{1}{4}cdot0,!706= 0,!176,~ p_2=frac{2cdot1}{3cdot4}cdot 0,!706=0,!118,

т.е. в установившемся, стационарном режиме в среднем 70,6% времени система будет находиться в состоянии S_0, 17,6% — в состоянии S_1 и 11,8% — в состоянии S_2.

Математический форум (помощь с решением задач, обсуждение вопросов по математике).

Кнопка "Поделиться"

Если заметили ошибку, опечатку или есть предложения, напишите в комментариях.

Уравнения Колмогорова. Предельные вероятности состояний

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

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

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

.                                                                   (8)

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

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

Вывести систему из этого состояния (см. граф на рис. 1) можно суммарным простейшим потоком с интенсивностью (l01+l02), т.е. в соответствии с (15.7), с вероятностью, приближенно равной (l01+l02)Dt. А вероятность того, что система не выйдет из состояния S0, равна [1-(l01+l02)Dt]. Вероятность того, что система будет находиться в состоянии S0, по первому способу (т.е. того, что находилась в состоянии S0 и не выйдет из него за время Dt), равна по теореме умножения вероятностей:

p0(t)·[1-(λ0102)*Δt].

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

Потоком интенсивностью l10 (или l 20см. рис. 1) система перейдет в состояние S0 с вероятностью, приближенно равной l10Dt (или l20Dt). Вероятность того, что система будет находиться в состоянии S0 по этому способу, равна р1(t)×l10Dt (или р2(t)×l20Dt).

Применяя теорему сложения вероятностей, получим

p0(t+Δt)=p1·λ10·Δt+p2(t)·λ20·Δt+p0(t)[1-(λ0102)·Δt],

откуда

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

p′0 = λ10·p120·p2+(λ1020)·p0,

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

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

(9)

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

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

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

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

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

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

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

(10)

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

Задача 2. Найти предельные вероятности для системы S задачи 1, граф состояний которой приведен на рис. 1, при l01=1, l02=2, l10=2, l13=2, l20=3, l23=1, l31=3, l32=2.

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

3p0=2p1+3p2  (11)

4p1=p0+3p3

4p2=2p0+2p3

p0+p1+p2+p3=1

(Здесь мы вместо одного “лишнего” уравнения системы (10) записали нормировочное условие (8)).

Решив систему (11), получим p0=0,40, p1=0,20, p2=0,27, p3=0,13, т.е. в предельном, стационарном режиме система S в среднем 40% времени будет находиться в состоянии S0 (оба узла исправны), 20% — в состоянии S1 (первый узел ремонтируется, второй работает), 27% — в состоянии S2 (второй узел ремонтируется, первый работает) и 13% времени — в состоянии S3 (оба узла ремонтируются).

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

Решение. Из задачи 2 следует, что в среднем первый узел исправно работает долю времени, равную p0+p3=0,40+0,27=0,67, а второй узел — p0+p1=0,40+0,20=0,60. В то же время первый узел находится в ремонте в среднем долю времени, равную p1+p3=0,20+0,13=0,33, а второй узел – p2+p3=0,27+0,13=0,40. Поэтому средний чистый доход в единицу времени от эксплуатации системы, т.е. разность между доходами и затратами, равен 

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

Уменьшение вдвое среднего времени ремонта каждого из узлов в соответствии с (6) будет означать увеличение вдвое интенсивностей потока “окончаний ремонтов” каждого узла, т.е. теперь l10=4, l 20=6, l31 =6, l32=4 и система линейных алгебраических уравнений (10), описывающая стационарный режим системы вместе с нормировочным условием (8) примет вид:

3p0=4p1+6p2

6p1=p0+6p3

7p2=2p0+4p3

p0+p1+p2+p3=1

Решив систему, получим p0=0,60, p1=0,15, p2=0,20, p3=0,05.

Учитывая, что p0+p2=0,60+0,20=0,80, p0+p1=0,60+0,15=0,75, p1+p3=0,15+0,05=0,20,
p2+p3=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%), то экономическая целесообразность ускорения ремонтов узлов очевидна.

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

Задача Интенсивности потоков
λ01 λ02 λ10 λ12 λ20 λ21
78 2 2 1 2 3 0

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

Размеченный граф состояний имеет вид.

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





p0(t) + p1(t) + p2(t) = 1

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







p0(t) + p1(t) + p2(t) = 1

Чтобы найти финальные вероятности состояний, в уравнениях Колмогорова отбросим первое уравнения, а по остальным составим систему алгебраических уравнений:

2p0-3p1 = 0

2p0+2p1-3p2=0

p0 + p1 + p2 = 1

Решим СЛАУ с помощью метода Гаусса.

Вывод: При достаточно большом времени работы техническое устройство с вероятностью p0 = 0.36 будет находиться в состоянии S0, с вероятностью p1 = 0.24 в состоянии S1 и с вероятностью p2 = 0.4 в состоянии S2.

Пример.

Техническое устройство может находиться в одном из трех состояний S0, S1, S2. Интенсивность потоков, которые переводят устройства из одного состояния во второе, известны λ01=2, λ10=4, λ21=2, λ12=3, λ20=4.

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

Размеченный граф состояний имеет вид.

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

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



Чтобы найти финальные вероятности состояний, в уравнениях Колмогорова отбросим первое уравнения, а по остальным составим систему алгебраических уравнений:

2p0-7p1+2p2=0

3p1-6p2=0

p0+p1+p2=1

Делим первое уравнение на 2, а второе на 3 и получим систему

p0-7p1+2p2=0

3p1-6p2=0

p0+p1+p2=1

Из третьего уравнения вычитаем первое

p0-3.5p1+p2=0

p1-2p2=0

4.5p1=1

Отсюда получим p1=0,22, p2=0,11 и p0=0,67.

Вывод: При достаточно большом времени работы техническое устройство с вероятностью p0 = 0,67 будет находиться в состоянии S0, с вероятностью p1 = 0,22 в состоянии S1 и с вероятностью p2 = 0,11 в состоянии S2.

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

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

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

Рис. 4

Рассмотрим упорядоченное множество состояний системы S0, S1, S2, …, Sk. Переходы могут осуществляться из любого состояния только в состояния с соседними номерами, т.е. из состояния Sk возможны переходы только либо в состояние Sk-1, либо в состояние Sk+1. (При анализе численности популяций считают, что состояние Sk соответствует численности популяции, равной k, и переход системы из состояния Sk в состояние Sk+1 происходит при рождении одного члена популяции, а переход в состояние Sk-1, — при гибели одного члена популяции).

Предположим, что все потоки событий, переводящие систему по стрелкам графа, простейшие с соответствующими интенсивностями lk, k+1 или lk+1, k.

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

В соответствии с правилом составления таких уравнений (см. 13) получим: для состояния S0

λ01p0 = λ10p1  (12)

для состояния S1 – (l12+l10)p1=l01 p0+l21p2, которое с учетом (12) приводится к виду

λ12p1 = λ21p2 (13)

Аналогично, записывая уравнения для предельных вероятностей других состояний, можно получить следующую систему уравнений:

(14)

к которой добавляется нормировочное условие

p0+p1+p2+…+pn=1 (15)

Решая систему (14), (15), можно получить

 (16)

(17)

Легко заметить, что в формулах (17) для p1, p2, …, pn коэффициенты при p0 есть слагаемые, стоящие после единицы в формуле (16). Числители этих коэффициентов представляют произведение всех интенсивностей, стоящих у стрелок, ведущих слева направо до данного состояния Sk (k=1, 2, …, n), а знаменатели — произведение всех интенсивностей, стоящих у стрелок, ведущих справа налево до состояния Sk.

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



Рис. 5

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



по (17) – т.е. в установившемся, стационарном режиме в среднем 70,6% времени система будет находиться в состоянии S0, 17,6% — в состоянии S1 и 11,8% — в состоянии S2.

При
анализе случайных процессов с
дискретными состояниями используют
граф
состояний
,
в котором каждое состояние S0,
S1,
…,
Sn
исследуемой системы изображают в
виде прямоугольника, а переход из
состояния Si
в состояние Sj
под воздействием простейших потоков
событий показывают дугой (стрелкой).
Если около дуг записаны интенсивности
потоков ij,
то граф состояния называют размеченным.

Пусть
двухканальная СМО может находиться
в одном из трех состояний: S0
– оба канала свободны; S1
– один канал занят обслуживанием
(любой из двух), а другой свободен;
S2
– оба канала заняты обслуживанием.
Переходы системы из состояния S0
в состояние S1
и из S1
в S2
происходят под воздействием простейших
потоков заявок с интенсивностями
01
и 12
соответственно, а из S1
в S0
и из S2
в S1
– под воздействием потоков обслуживания
с интенсивностями 01
и 12.
Тогда размеченный граф состояний
СМО имеет вид, изображенный на рис.27.1.

Рис.27.1. Размеченный
граф состояний

Пребывание
СМО в том или ином состоянии носит
вероятностный характер. Обозначим
pi(t)
вероятность, что в момент времени t
СМО находится в состоянии Si.
Поскольку нахождение СМО в состояниях
S0,
S1
и S2
образует полную группу событий, то


.

Зададим
малое приращение времени t
и найдем вероятность того, что СМО
в момент времени t
+ t
будет находиться в состоянии S1.
Это достигается разными вариантами.

а)
В момент t
СМО находилась в состоянии S1.
Она могла быть выведена из этого
состояния простейшим потоком с
суммарной интенсивностью 10
+ 12.
Как показано выше, вероятность этого
приближенно равна

(10
+ 12)t.

Тогда
вероятность невыхода из S1
приближенно равна

1
– (10
+ 12)t.

Согласно
теореме умножения вероятностей,
вероятность того, что СМО останется
в состоянии S1
в момент t
+
t
приближенно равна:

.

б)
В момент t
СМО находилась в состоянии S0.
Переход СМО из S0
в S1
происходит под воздействием потока
с интенсивностью 01
с вероятностью приближенно равной
01t.
Тогда вероятность, что при этом СМО
будет находиться в состоянии S1
в момент t
+ t
в этом варианте приближенно равна
p0(t)01t.

в)
В момент t
СМО находилась в состоянии S2.
Переход СМО из S2
в S1
происходит под воздействием потока
с интенсивностью 21
с вероятностью приближенно равной
21t.
Тогда вероятность, что при этом СМО
будет находиться в состоянии S1
в момент t
+ t
в этом варианте приближенно равна
p2(t)21t.

Применяя
теорему сложения вероятностей для
всех этих вариантов, можно записать:

.

Раскрыв
скобки, перенеся p1(t)
в правую часть приближенного равенства
и поделив уравнение на t,
получим:

.

При
переходе к пределу при t

0 приближенные равенства перейдут в
точные и получится дифференциальное
уравнение:

.

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

В
общем случае система дифференциальных
уравнений Колмогорова может быть
построена по следующему алгоритму:

  • В
    левой части уравнения стоит производная
    по времени вероятности i-го
    состояния.

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

Такое
дифференциальное уравнение записывается
для каждого состояния. При задании
начальных условий можно решить систему
дифференциальных уравнений Колмогорова
и найти изменение во времени
вероятностей состояний pi(t).

Можно
показать, что при t


эти вероятности независимо от
начальных условий стремятся к некоторым
предельным значениям, которые
называются предельными
(или финальными) вероятностями
СМО. Смысл предельных вероятностей
состоит в том, что они показывают
среднее относительное время нахождения
СМО в данном состоянии.

Так
как предельные вероятности не
изменяются во времени, то, заменяя в
уравнениях Колмогорова производные
нулями, получим систему алгебраических
уравнений для определения предельных
вероятностей.

Можно показать,
что ранг матрицы такой системы на
единицу меньше числа уравнений,
поэтому при решении следует отбросить
одно из уравнений и добавить так
называемое условие нормировки

.

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

Пример.
Расчетный узел в отделе игрушек
детского универмага состоит из 2
кассовых аппаратов. Размеченный граф
состояний такой СМО имеет вид:

S0
– обе кассы свободны;
S1
– одна касса свободна, а другая
занята;
S2
– обе кассы заняты. Найти относительное
время пребывания СМО в каждом
состоянии, если
01
= 60;
12
= 48;
10
=
21
= 36 (покупателей/час).

Решение.
Относительное время пребывания СМО
в каждом состоянии равно предельным
вероятностям состояний. Составляем
систему уравнений:

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

откуда

Решая, находим

.

Следовательно,
20,4 % рабочего времени обе кассы
свободны, 34,1 % времени занята одна
касса и 45,5 % времени заняты обе кассы.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]

  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #

Уравнения Колмогорова.
Предельные вероятности состояний

Рассмотрим математическое описание марковского процесса с дискретными состояниями и непрерывным временем* на примере случайного процесса из примера 1, граф которого изображен на рис. 1. Будем полагать, что все переходы системы из состояния в происходят под воздействием простейших потоков событий с интенсивностями ; так, переход системы из состояния в будет происходить под воздействием потока отказов первого узла, а обратный переход из состояния в — под воздействием потока “окончаний ремонтов” первого узла и т.п.

Граф состояний системы с проставленными у стрелок интенсивностями будем называть размеченным (см. рис. 1). Рассматриваемая система имеет четыре возможных состояния: .

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

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

1. Система в момент с вероятностью находилась в состоянии , а за время не вышла из него.

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

2. Система в момент с вероятностями (или ) находилась в состоянии или и за время перешла в состояние .

Потоком интенсивностью (или — с- рис. 1) система перейдет в состояние с вероятностью, приближенно равной (или ). Вероятность того, что система будет находиться в состоянии по этому способу, равна (или ).

Применяя теорему сложения вероятностей, получим

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

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

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

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

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

Особенность решения дифференциальных уравнений вообще состоит в том, что требуется задать так называемые начальные условия, т.е. в данном случае вероятности состояний системы в начальный момент . Так, например, систему уравнений (9) естественно решать при условии, что в начальный момент оба узла исправны и система находилась в состоянии , т.е. при начальных условиях .

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

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

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

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

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

Пример 2. Найти предельные вероятности для системы из примера 1, граф состояний которой приведен на рис. 1, при

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

(Здесь мы вместо одного “лишнего” уравнения системы (10) записали нормировочное условие (8)).

Решив систему (11), получим , т.е. в предельном, стационарном режиме система в среднем 40% времени будет находиться в состоянии (оба узла исправны), 20% — в состоянии (первый узел ремонтируется, второй работает), 27% — в состоянии (второй узел ремонтируется, первый работает) и 13% времени — в состоянии (оба узла ремонтируются)

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

Решение. Из примера 2 следует, что в среднем первый узел исправно работает долю времени, равную , а второй узел — . В то же время первый узел находится в ремонте в среднем долю времени, равную , а второй узел — . Поэтому средний чистый доход в единицу времени от эксплуатации системы, т.е. разность между доходами и затратами, равен

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

Решив систему, получим .

Учитывая, что , а затраты на ремонт первого и второго узла составляют теперь соответственно 8 и 4 ден.ед., вычислим средний чистый доход в единицу времени:

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

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

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

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

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

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

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

В соответствии с правилом составления таких уравнений (см. 13) получим: для состояния

для состояния имеем , которое с учетом (12) приводится к виду

Аналогично, записывая уравнения для предельных вероятностей других состояний, можно получить следующую систему уравнений:

к которой добавляется нормировочное условие

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

Решая систему (14), (15), можно получить

Легко заметить, что в формулах (17) для коэффициенты при есть слагаемые, стоящие после единицы в формуле (16). Числители этих коэффициентов представляют произведение всех интенсивностей, стоящих у стрелок, ведущих слева направо до данного состояния , а знаменатели — произведение всех интенсивностей, стоящих у стрелок, ведущих справа налево до состояния .

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

Потоки событий. Уравнения Колмогорова

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

Основные понятия и классификация систем массового

Обслуживания

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

Всякая СМО включает четыре элемента: входящий поток требований, их очередь на обслуживание, обслуживающее устройство, выходящий поток обслуженных требований.

Требованием (клиентом, заявкой) в СМО называется каж­дый отдельный запрос на выполнение какой-либо работы. Обслуживание — это выполнение работы по удовлетворе­нию поступившего требования. Объект, выполняющий обслужи­вание требований, называется обслуживающим устройством (прибором) или каналом обслуживания.

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

СМО классифицируются по разным признакам.

1. По числу каналов обслуживания СМО делятся на одноканальные и многоканальные.

2. В зависимости от условий ожидания требованием нача­ла обслуживания различают СМО с потерями (отказами) и СМО с ожиданием.

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

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

Для СМО с ожиданием основными характеристиками яв­ляются математические ожидания длины очереди и времени ожидания. Примером системы с ожиданием может служить процесс восстановления телевизоров в ремонтной мастерской.

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

В качестве характеристик эффективности может приме­няться вероятность отказа как в системах с потерями (или харак­теристики времени ожидания) и в системах с ожиданием.

3. По дисциплине обслуживания СМО делятся на систе­мы с приоритетом в обслуживании и на системы без приоритета в обслуживании.

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

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

5. По месту нахождения источника требований СМО де­лятся на разомкнутые (когда источник требования находится вне системы) и замкнутые (когда источник находится в самой сис­теме).

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

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

Методы и модели исследования СМО можно условно раз­бить на аналитические и статистические (имитационного моде­лирования процессов массового обслуживания).

Аналитические методы позволяют получить характеристи­ки системы как некоторые функции от параметров ее функционирования. Благодаря этому появляется возможность проводить качественный анализ влияния отдельных факторов на эффектив­ность работы СМО.

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

Потоки событий. Уравнения Колмогорова

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

Для простейшего потока частота поступления требований в систему подчиняется закону Пуассона, то есть вероятность поступления за время t, ровно kтребований задается форму­лой:

Pk= e – λt

где λ – параметр потока (интенсивность – среднее число требований, поступивших в систему за единицу времени).

Простейший поток обладает тремя основными свойствами: ординарностью, стационарностью и отсутствием последействия.

Ординарность потока означает практическую невозмож­ность одновременного поступления двух и более требований. Например, достаточно малой является вероятность того, что из группы станков, обслуживаемых бригадой ремонтников, одно­временно выйдут из строя несколько станков.

Стационарным называется поток, вероятные характеристики которого не зависят от времени. Например, математи­ческое ожидание числа требований, поступающих в систему в единицу времени, является величиной постоянной и называется интенсивностью потока. Таким образом, вероятность поступления в систему определен­ного количества требований в течение заданного промежутка времени t зависит от его величины и не зависит от начала его отсчета на оси времени.

Отсутствие последействия означает, что число требований, поступивших в систему до момента t, не определяет того, сколько требований поступит в систему за время t + t.

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

Вероятность того, что на временном интервале t = τ не поступит ни одного требования p0 определяется как:

Тогда вероятность того, что на этом же временном интервале появится хотя бы одно требование определяется соотношением

p(t) = 1p0 = 1e – λτ .

Вероятность попадания на элементарный временной интервал, т.е. когда τ = Δt, хотя бы одного события (требования) потока, можно определить, заменив функцию e – λτ двумя первыми числами её разложения в ряд Маклорена по степеням Δt. Действительно:

p(Δt) = 1e – λΔt ≈ λ Δt.

Важной характеристикой СМО является время обслужива­ния требований в системе. Время обслуживания является, как правило, случайной величиной и, следовательно, может быть описано законом распределения. Наибольшее распространение в теории и, особенно в практических приложениях, получил экс­поненциальный закон. Для этого закона функция распределения вероятностей имеет вид:

р (tt ,

т.е. вероятность того, что время обслуживания не превосходит некоторой величины t, где μ – интенсивность потока обслуживания тре­бований в системе (среднее число требований, обслуженных в единицу времени).

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

p(Δt) = 1e -Δ t ≈ μλ Δt.

При анализе случайных процессов с дискретными состояниями пользуются графиком состояний, где прямоугольником изображаются состояния, а переходы из состояния в состояние – стрелками. У стрелок обычно проставляются значения интенсивностей λijij) перехода системы из состояния Si в состояние Sj, которые происходят под воздействием простейших потоков событий.

Рассмотрим систему, которая может находиться в двух состояниях: S0– система исправна; S1 – система находится в состоянии отказа и ремонтируется (рис. 8.1).

Рис. 8.1. Граф системы

Будем характеризовать состояние системы (S0,S1) вероятностями состояния р0 (t) и р1 (t). Очевидно, что

Найдем вероятность того, что система в момент (t+Δt) будет находиться в состоянии S0.

Это возможно, во-первых, в том случае, если система в момент времени t с вероятностью р0(t)находилась в состоянии S0 и за время Δtиз него не вышла. Вероятность выхода системы за время Δt из состояния S0 в состояние S1 определяется как λ· Δt. Противоположная вероятность (что система не выйдет из S0) определяется как (1 – λΔt). Вероятность того, что система, находившаяся в состоянии S0 с вероятностью р0(t), за время Δt не выйдет из него, равна по теореме умножения вероятностей

Во-вторых, система в момент времени t находилась с вероятностью р1(t)в состоянии S1 и за интервал времени Δt с вероятностью μ·Δt перешла в состояние S0, т.е.

Вероятность р0(t + Δt) нахождения системы в состоянии S0 момент времени (t + Δt)по какому-либо из двух рассмотренных способов равна сумме рассмотренных вероятностей

Преобразуем соотношение к виду

Переходя к пределу при Δt→0, получим

По аналогии составляется уравнение, описывающее вероятность того, что система в момент (t + Δt)будет находиться в состоянии S1, но проще это найти из условия нормирования

С учетом этого условия система уравнений для двух состояний графа имеет вид

Задав начальные условия, можно решить систему уравнений и найти систему функций времени рi(t), где i – номер состояния.

Для достаточно большого значения tраспределение вероятностей стабилизируется и практически не зависит от времени, т.е.

Тогда система уравнений упрощается (стационарный режим)

Откуда вероятности состояний установившегося процесса определяются как

В частности, если μ=2; λ=1, то Р0=0,67; Р1=0,33. Таким образом, в среднем система будет находиться в рабочем состоянии 67%, а в состоянии ремонта 33% времени.

В общем случае система уравнений Колмогорова может быть составлена по следующему алгоритму.

1. В левой части каждого уравнения стоит производная вероятности i-го состояния.

2. В правой части каждого уравнения стоит:

2.1. Сумма произведений вероятностей всех состояний, из которых идут стрелки в i-е состояние, на интенсивности соответствующих потоков событий.

2.2. Минус суммарная интенсивность всех потоков, выводящих систему из данного состояния, умноженная на вероятность i-го состояния.

Система Колмогорова состоит из независимых уравнений и условия нормирования.

СМО с ожиданием

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

Задачи с очередями являются типичными в производствен­ных условиях, например при организации наладочных и ремон­тных работ, при многостаночном обслуживании и т.д.

Постановка задачи в общем виде выглядит следующим образом.

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

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

Замкнутая СМО с ожиданием

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

Такую систему можно классифицировать как многоканальную СМО (n – каналов) и ограниченной данной очереди l, причем n +l = m.

Граф состояний такой системы изображен на рис. 8.2.

Рис. 8.2. Граф состояний многоканальной СМО с ограниченной очередью

Состояния данной системы означают:

S0 – отсутствие требований в системе;

S1 – одно требование обслуживается, очереди нет;

S2 – два требования обслуживаются, очереди нет;

Snn требований обслуживаются, очереди нет;

Sn+1n требований обслуживаются, одно требование стоит в очереди;

Sn+ln требований обслуживаются, l требований стоят в очереди.

Система уравнений вероятностей состояний в стационарном режиме для цепочки S0 Sn будет:

Для цепочки состояний Sn+1Sn+lсистема уравнений стационарного режима будет:

В качестве основных критериев, характеризующих качество функционирования рассматриваемой системы, выберем: 1) от­ношение средней длины очереди к наибольшему числу требова­ний, находящихся одновременно в обслуживающей системе – коэффициент простоя обслуживаемого объекта; 2) отношение среднего числа незанятых обслуживающих каналов к их общему числу – коэффициент простоя обслуживаемого канала.

Рассмотрим расчет необходимых вероятностных характе­ристик (показателей качества функционирования) замкнутой СМО.

1. Вероятность того, что в системе находится k требований при условии, когда их число не превышает числа обслуживаю­щих аппаратов п:

Pk=αkP0, (1 k n),

где αk= или αk= ;

;

– интенсивность поступления требований в си­стему от одного источника;

об – средняя продолжительность обслуживания одного требования;

т – наибольшее возможное число требований, находящих­ся в обслуживающей системе одновременно (m=n+l);

п– число обслуживающих аппаратов;

Р0 – вероятность того, что все обслуживающие аппараты свободны.

2. Вероятность того, что в системе находится k требований при условии, когда их число больше числа обслуживающих ап­паратов:

Pk= kP0, (n k m),

где k= ( об) k или αk= .

3. Вероятность того, что все обслуживающие аппараты сво­бодны, определяется из условия

4. Среднее число требований, ожидающих начала обслужи­вания (средняя длина очереди):

.

5. Коэффициент простоя требования в ожидании обслужи­вания:

a1= .

6. Вероятность того, что все обслуживающие аппараты заняты:

Pотк= .

7. Среднее число требований, находящихся в обслуживающей системе (обслуживаемых и ожидающих обслуживания):

A2= = .

8. Коэффициент полного простоя требований на обслуживании и в ожидании обслуживания:

a2= .

9. Среднее время простоя требования в очереди на обслуживание:

Tож=a1/ .

10. Среднее число свободных обслуживающих аппаратов:

A3= .

11 . Коэффициент простоя обслуживающих аппаратов:

a3 = .

12. Вероятность того, что число требований, ожидаю­щих обслуживания, больше некоторого числа В (вероятность того, что в очереди на обслуживание находится более В тре­бований):

PB = = .

Рассмотрим пример расчета характеристик замкнутой СМО.

Пример 8.1. Оптовый склад строительных материалов обслу­живает шесть предприятий-потребителей материалов. Каждый из потребителей направляет на склад автомашину за мате­риалами в среднем один раз в смену (продолжительность сме­ны 8 ч). На складе имеется один автопогрузчик, который исполь­зуется только для погрузки материалов на прибывающие авто­машины. Прибывшая на склад автомашина становится в оче­редь, если автопогрузчик занят погрузкой другой автомашины. Обработка статистических данных о продолжительности погруз­ки одной автомашины и проверка соответствующей гипотезы по­казали, что продолжительность погрузки одной автомашины подчиняется показательному закону распределения и составля­ет в среднем 48 мин (0,1 смены). Статистическое исследование потока автомашин показа­ло, что число автомашин, поступающих на склад в единицу вре­мени, подчиняется пуассоновскому закону распределения. Требуется провести расчет характеристик функционирова­ния приведенной производственной системы как СМО.

Решение. Рассчитаем основные параметры системы для условий задачи.

Вероятность того, что все обслуживающие ап­параты свободны (на складе нет автомашин) определяется как P0, λ=1, μ=0,1.

Вероятность того, что на складе одна автомашина:

P1= 0,1P0=0,6P0,

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

P2= 0,1 2 P0=0,3P0.

Рассчитывая аналогично, получим: Р3=0,12Р0; Р4=0,036Р0;Р5= 0,0072Р0; Р6= 0,0007Р0. Так как сумма вероятнос­тей нахождения системы в любом из состояний равна 1, т.е.

=1,

то P0(1 + 0,6 + 0,3 + 0,12 + 0,036 + 0,0072 + 0,0007) = 2,0639; Р0 = 1.

Отсюда находим Р0 = 0,4845.

Дальнейшие расчеты затруднений не вызывают. Напри­мер, средняя длина очереди равна

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

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

Содержание:

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

Теория случайных процессов – это раздел математической науки, который изучает закономерности случайных явлений в динамике их развития.

Определение случайного процесса и его характеристики

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

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

Количество реализаций определенного случайного процесса изображено на рис. 4.1. Пусть сечение процесса при данном является непрерывной случайной величиной. Тогда случайный процесс при данном определяется плотностью вероятности

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

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

Таких сечений бесконечно много, но для задания случайного процесса удается ограничиться сравнительно небольшим количеством сечений.

Случайный процесс имеет порядок если он полностью определяется плотностью общего распределения произвольных сечений процесса, то есть плотностью -мерной случайной величины где – сечение случайного процесса в момент времени

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

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

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

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

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

Корреляционной функцией случайного процесса называется детерминированная функция

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

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

Поэтому рассматривается также нормированная корреляционная функция случайного процесса.

Нормированной корреляционной функцией случайного процесса называется функция

Пример. Случайный процесс определяется формулой где – случайная величина. Найти основные характеристики этого процесса, если

Решение. Согласно свойствам математического ожидания и дисперсии получим:

Находим далее корреляционную функцию

а также нормированную корреляционную функцию

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

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

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

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

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

Заявки поступают в СМО конечно нерегулярно, а случайно, образуя так называемый случайный поток заявок (ссылок). Обслуживание заявок также длится в течение определенного случайного времени. Учитывая случайность потока заявок и время обслуживания, СМО загружаются неравномерно: в определенные периоды накапливается очень много заявок (они или стают в очередь, или оставляют СМО не обслуженными), в другие периоды СМО работает с малой загрузкой или простаивает.

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

Показателями эффективности СМО являются:

  • – среднее количество заявок, которые она обслуживает за единицу времени;
  • – среднее количество заявок в очереди;
  • – среднее время ожидания обслуживания;
  • – вероятность отказа в обслуживании без ожидания;
  • – вероятность того, что количество заявок в очереди превышает определенное значение и т.д.

СМО делятся на два основных класса: СМО с отказами и СМО с ожиданием (очередью).

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

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

Процесс работы СМО представляет собой случайный процесс.

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

Процесс функционирования СМО представляет собой случайный процесс с дискретными состояниями и непрерывным временем.

Математический анализ работы СМО существенно упрощается, если процесс этой работы – марковский.

Понятие марковского процесса

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

Пример. Система – счетчик в такси. Состояние системы в момент характеризуется количеством километров, пройденных автомобилем до данного момента. Пусть в момент счетчик показывает Вероятность того, что в момент счетчик будет показывать то или иное количество километров зависит от но не зависит от того, в какие моменты времени изменялись показатели счетчика до момента

Некоторые процессы можно приблизительно считать марковскими.

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

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

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

Решение. Возможные состояния системы: – оба узла исправны; – первый узел ремонтируется, а второй исправный; второй узел ремонтируется, а первый исправный; – оба узла ремонтируются.

Граф системы приведен на рис. 4.2.

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

Для математического описания марковского случайного процесса с дискретными состояниями и непрерывным временем, которое происходит в СМО, рассмотрим одно из важных понятий теории вероятностей – понятие потока событий.

Простейший поток событий

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

Поток называется простейшим, если он имеет такие свойства:

1) стационарность – вероятность того, что за некоторый промежуток времени произойдет то или иное количество событий, зависит только от длины промежутка и не зависит от начала его отсчета, то есть интенсивность потока постоянная;

2) отсутствие последействия – вероятность наступления некоторого количества событий в произвольном промежутке времени не зависит от того, какое количество событий произошло до начала этого промежутка;

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

Если поток событий простейший, то вероятность того, что за промежуток времени событие наступит раз, определяется формулой: где – интенсивность потока. Эта формула отражает все свойства простейшего потока, а следовательно, является его математической моделью.

Пример. Среднее количество заявок, поступающих на комбинат бытового обслуживания за 1 час равно 4. Найти вероятность того, что за 3 часа поступит: 1) 6 заявок; 2) менее 6 заявок; 3) не менее 6 заявок.

Решение. Пусть событие – “поступление одной заявки”. Поток заявок простейший. Поэтому для решения задачи используем приведенную только что формулу, в которой Вычислим соответствующие вероятности:

Уравнения Колмогорова. Предельные вероятности состояний

Вероятностью состояния называется вероятность того, что в момент система будет находиться в состоянии

Очевидно, что для любого момента сумма вероятностей состояний равна 1:

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

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

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

Особенность решения дифференциальных уравнений вообще состоит в том, что нужно задавать так называемые начальные условия, в данном случае – вероятности состояний системы в начальный момент Так, систему (2) должны решать при условии, что в начальный момент оба узла исправны и система находилась в состоянии то есть при начальных условиях

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

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

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

Пример 1. Найти предельные вероятности для системы из последнего примера, граф состояний которой приведен на рис. 4.2. При

Решение. Система алгебраических уравнений, которая описывает стационарный режим для данной системы, принадлежит к виду (1):

Решая эту систему уравнений, получаем Следовательно, в предельном стационарном режиме система в среднем 40% времени находится в состоянии 20% – в состоянии 27% – в состоянии 13% – в состоянии

Пример 2. Найти прибыль от эксплуатации в стационаром режиме системы когда известно, что за единицу времени исправная работа первого и второго узлов приносит доход, который составляет соответственно 10 и 6 ус. ед., а их ремонт требует расходов, которые составляют соответственно 4 и 2 ус. ед.

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

Решение. Из примера 1 следует, что в среднем первый узел исправен в течение части времени, которая составляет а второй узел – в течение части В этом случае первый узел находится в ремонте в среднем часть времени, равной а второй – Поэтому средняя прибыль за единицу времени от эксплуатации системы (разница между доходом и расходами) будет такой:

Прибыль = (ус. ед.).

Уменьшение вдвое среднего времени ремонта каждого из узлов согласно с будет означать увеличение вдвое интенсивности потока “окончания ремонта” каждого узла. Следовательно, в этом случае и система линейных алгебраических уравнений (1) принимает вид:

Решая эту системы, получаем

Поскольку то расходы на ремонт первого и второго узла будут составлять соответственно 8 и 4 ус. ед. Отсюда получим среднюю прибыль за единицу времени:

(Прибыль) (ус. ед.)

(Прибыль) больше, чем Прибыль (приблизительно – на 2%), поэтому экономическая целесообразность сокращения срока ремонта узлов очевидна.

Лекции:

Присылайте задания в любое время дня и ночи в ➔

Официальный сайт Брильёновой Натальи Валерьевны преподавателя кафедры информатики и электроники Екатеринбургского государственного института.

Все авторские права на размещённые материалы сохранены за правообладателями этих материалов. Любое коммерческое и/или иное использование кроме предварительного ознакомления материалов сайта natalibrilenova.ru запрещено. Публикация и распространение размещённых материалов не преследует за собой коммерческой и/или любой другой выгоды.

Сайт предназначен для облегчения образовательного путешествия студентам очникам и заочникам по вопросам обучения . Наталья Брильёнова не предлагает и не оказывает товары и услуги.

[spoiler title=”источники:”]

http://lektsii.org/7-1858.html

http://natalibrilenova.ru/teoriya-sluchajnyih-protsessov-i-teoriya-massovogo-obsluzhivaniya/

[/spoiler]

Макеты страниц

Пусть имеется физическая система S с дискретными состояниями:

в которой протекает марковский случайный процесс с непрерывным временем (непрерывная цепь Маркова). Граф состояний показан на рис. 4.32.

Рис. 4.32

Предположим, что все интенсивности потоков событий, переводящих систему из состояния в состояние, постоянны:

другими словами, все потоки событий — простейшие (стационарные пуассоновские) потоки.

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

при любом t дающих в сумме единицу:

Поставим теперь следующий вопрос: что будет происходить с системой S при Будут ли функции стремиться к каким-то пределам? Эти пределы, если они существуют, называются предельными (или «финальными») вероятностями состояний.

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

Рис. 4.33

Рис. 4.34

На рис. 4.33 показан граф состояний, удовлетворяющий поставленному условию: из любого состояния система может рано или поздно перейти в любое другое. Напротив, для системы, граф состояний которой показан на рис. 4.34, условие не выполнено. Очевидно, что если начальное состояние такой системы то, например, состояние при может быть достигнуто, а если начальное состояние — не может.

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

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

Очевидно, предельные вероятности состояний, так же как и допредельные, в сумме должны давать единицу:

Таким образом, при в системе S устанавливается некоторый предельный стационарный режим: он состоит в том, что система случайным образом меняет свои состояния, но вероятность каждого из них уже не зависит от времени: каждое из состояний осуществляется с некоторой постоянной вероятностью. Каков смысл этой вероятности? Она представляет собой не что иное, как среднее относительное время пребывания системы в данном состоянии. Например, если у системы S три возможных состояния: причем их предельные вероятности равны 0,2, 0,3 и 0,5, это означает, что после перехода к установившемуся режиму система S в среднем две десятых времени будет находиться в состоянии три десятых — в состоянии и половину времени в состоянии Возникает вопрос: как вычислить предельные вероятности состояний

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

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

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

(так называемым «нормировочным условием») эти уравнения дают возможность вычислить все предельные вероятности

Пример 1. Физическая система 5 имеет возможные состояния: размеченный граф которых дан на рис. 4.35 (у каждой стрелки поставлено численное значение соответствующей интенсивности). Вычислить предельные вероятности состояний:

Решение. Пишем уравнения Колмогорова для вероятностей состояний:

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

Уравнения (7.4) — так называемые однородные уравнения (без свободного члена). Как известно из алгебры, эти уравнения определяют величины только с точностью до постоянного множителя. К счастью, у нас есть нормировочное условие:

которое, совместно о уравнениями (7.4), дает возможность найти все неизвестные вероятности.

Рис. 4.35

Рис. 4.36

Рис. 4.37

Действительно, выразим из (7.4) все неизвестные вероятности через одиу из них, например, через Из первого уравнения:

Подставляя во второе уравнение, получим:

Четвертое уравнение дает:

Подставляя все эти выражения вместо в нормировочное условие (7.5), получим

Отсюда

Таким образом, предельные вероятности состояний получены, они равиы:

Это значит, что в предельном, установившемся режиме система S будет проводить в состоянии в среднем одну двадцать четвертую часть времени, в состоянии — половину времени, в состоянии — пять двадцать четвертых и в состоянии — одну четверть времени.

Заметим, что решая эту задачу, мы совсем не пользовались одним из уравнений (7 4) — третьим Нетрудно убедиться, что оно является следствием трех остальных: складывая все четыре уравнения, мы получим тождественный нуль. С равным успехом, решая систему, мы могли бы отбросить любое из четырех уравнений (7.4).

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

Пример 2. Граф состояний системы показан на рис. 4.36. Написать алгебраические уравнения для предельных вероятностей состояний.

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

Чтобы в дальнейшем сразу же писать такие уравнения, полезно Запомнить следующее мнемоническое правило: «что втекает, то и вытекает», то есть для каждого состояния сумма членов, соответствующих входящим стрелкам, равна сумме членов, соответствующих выходящим, каждый член равен интенсивности потока событий, переводящего систему по данной стрелке, умноженной на вероятность того состояния, из которого выходит стрелка.

В дальнейшем мы во всех случаях будем пользоваться именно этим кратчайшим способом записи уравнений для предельных вероятностей.

Пример 3. Написать алгебраические уравнения для предельных вероятностей состояний системы S, граф состояний которой дан на рис. 4.37. Решить эти уравнения.

Решение. Пишем алгебраические уравнения для предельных вероятностей состояний;

Нормировочное условие!

Выразим с помощью первых двух уравнений через

Подставим их в нормировочное условие (7.9):

откуда

Далее, из (7 10) получим

1

Оглавление

  • ПРЕДИСЛОВИЕ
  • ВВЕДЕНИЕ
  • 1. ОСНОВНЫЕ ПОНЯТИЯ ИССЛЕДОВАНИЯ ОПЕРАЦИЙ
  • 1. ОПЕРАЦИЯ. ЭФФЕКТИВНОСТЬ ОПЕРАЦИИ
  • 2. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ОПЕРАЦИИ
  • 3. ОБЩАЯ ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЯ ОПЕРАЦИЙ. ДЕТЕРМИНИРОВАННЫЙ СЛУЧАЙ
  • 4. ОБЩАЯ ПОСТАНОВКА ЗАДАЧИ ИССЛЕДОВАНИЯ ОПЕРАЦИЙ. ОПТИМИЗАЦИЯ РЕШЕНИЯ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ
  • 5. ОЦЕНКА ОПЕРАЦИИ ПО НЕСКОЛЬКИМ ПОКАЗАТЕЛЯМ
  • 2. ЛИНЕЙНОЕ ПРОГРАММИРОВАНИЕ
  • 1. ЗАДАЧИ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ
  • 2. Задача о загрузке станков.
  • 3. Задача о распределении ресурсов.
  • 4. Задача о перевозках.
  • 5. Задача о производстве сложного оборудования.
  • 2. ОСНОВНАЯ ЗАДАЧА ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ
  • 3. ГЕОМЕТРИЧЕСКАЯ ИНТЕРПРЕТАЦИЯ ОСНОВНОЙ ЗАДАЧИ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ
  • 4. ЗАДАЧА ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ С ОГРАНИЧЕНИЯМИ-НЕРАВЕНСТВАМИ. ПЕРЕХОД ОТ НЕЕ К ОЗЛП И ОБРАТНО
  • 5. СИМПЛЕКС-МЕТОД РЕШЕНИЯ ЗАДАЧИ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ
  • 6. ТАБЛИЧНЫЙ АЛГОРИТМ ЗАМЕНЫ БАЗИСНЫХ ПЕРЕМЕННЫХ
  • 7. ОТЫСКАНИЕ ОПОРНОГО РЕШЕНИЯ ОСНОВНОЙ ЗАДАЧИ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ
  • 8. ОТЫСКАНИЕ ОПТИМАЛЬНОГО РЕШЕНИЯ ОСНОВНОЙ ЗАДАЧИ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ
  • 9. ТРАНСПОРТНАЯ ЗАДАЧА ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ
  • 10. НАХОЖДЕНИЕ ОПОРНОГО ПЛАНА
  • 11. УЛУЧШЕНИЕ ПЛАНА ПЕРЕВОЗОК. ЦИКЛ ПЕРЕСЧЕТА
  • 12. РЕШЕНИЕ ТРАНСПОРТНОЙ ЗАДАЧИ МЕТОДОМ ПОТЕНЦИАЛОВ
  • 13. ТРАНСПОРТНАЯ ЗАДАЧА С НЕПРАВИЛЬНЫМ БАЛАНСОМ
  • 2. ТЗ с избытком заявок
  • 14. РЕШЕНИЕ ТРАНСПОРТНОЙ ЗАДАЧИ ПО КРИТЕРИЮ ВРЕМЕНИ
  • 3. ДИНАМИЧЕСКОЕ ПРОГРАММИРОВАНИЕ
  • 1. ЗАДАЧИ ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ
  • 2. ЗАДАЧА О НАБОРЕ ВЫСОТЫ И СКОРОСТИ ЛЕТАТЕЛЬНЫМ АППАРАТОМ
  • 3. ОБЩАЯ ПОСТАНОВКА ЗАДАЧИ ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ. ИНТЕРПРЕТАЦИЯ УПРАВЛЕНИЯ В ФАЗОВОМ ПРОСТРАНСТВЕ
  • 4. ЗАДАЧИ РАСПРЕДЕЛЕНИЯ РЕСУРСОВ
  • 5. ПРИМЕР РЕШЕНИЯ ЗАДАЧИ РАСПРЕДЕЛЕНИЯ РЕСУРСОВ
  • 6. ДРУГИЕ ЗАДАЧИ РАСПРЕДЕЛЕНИЯ РЕСУРСОВ
  • 1. Распределение ресурсов по неоднородным этапам
  • 2. Задача о резервировании ресурсов
  • 3. Задача распределения ресурсов между тремя и более отраслями
  • 7. РАСПРЕДЕЛЕНИЕ РЕСУРСОВ СО ВЛОЖЕНИЕМ ДОХОДОВ В ПРОИЗВОДСТВО
  • 8. РЕШЕНИЕ ЗАДАЧИ ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ С УЧЕТОМ ПРЕДЫСТОРИИ ПРОЦЕССА
  • 9. ЗАДАЧИ ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ, НЕ СВЯЗАННЫЕ СО ВРЕМЕНЕМ
  • 10. ЗАДАЧИ ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ С МУЛЬТИПЛИКАТИВНЫМ КРИТЕРИЕМ
  • 11. БЕСКОНЕЧНОШАГОВЫЙ ПРОЦЕСС ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ
  • 4. МОДЕЛИРОВАНИЕ ОПЕРАЦИЙ ПО СХЕМЕ МАРКОВСКИХ СЛУЧАЙНЫХ ПРОЦЕССОВ
  • 1. МАРКОВСКИЙ СЛУЧАЙНЫЙ ПРОЦЕСС С ДИСКРЕТНЫМИ СОСТОЯНИЯМИ
  • 2. СЛУЧАЙНЫЕ ПРОЦЕССЫ С ДИСКРЕТНЫМ И НЕПРЕРЫВНЫМ ВРЕМЕНЕМ. МАРКОВСКАЯ ЦЕПЬ
  • 3. МАРКОВСКИЙ ПРОЦЕСС С ДИСКРЕТНЫМИ СОСТОЯНИЯМИ И НЕПРЕРЫВНЫМ ВРЕМЕНЕМ. УРАВНЕНИЯ КОЛМОГОРОВА ДЛЯ ВЕРОЯТНОСТЕЙ СОСТОЯНИЙ
  • 4. ПОТОК СОБЫТИИ. ПРОСТЕЙШИЙ ПОТОК И ЕГО СВОЙСТВА
  • 5. ПОТОКИ ПАЛЬМА. ПОТОКИ ЭРЛАНГА
  • 6. ПУАССОНОВСКИЕ ПОТОКИ СОБЫТИЙ И НЕПРЕРЫВНЫЕ МАРКОВСКИЕ ЦЕПИ
  • 7. ПРЕДЕЛЬНЫЕ ВЕРОЯТНОСТИ СОСТОЯНИЙ
  • 8. ПРОЦЕСС «ГИБЕЛИ И РАЗМНОЖЕНИЯ»
  • 9. ЦИКЛИЧЕСКИЙ ПРОЦЕСС
  • 10. ПРИБЛИЖЕННОЕ СВЕДЕНИЕ НЕМАРКОВСКИХ ПРОЦЕССОВ К МАРКОВСКИМ. МЕТОД «ПСЕВДОСОСТОЯНИЙ»
  • 5. ТЕОРИЯ МАССОВОГО ОБСЛУЖИВАНИЯ
  • 1. ЗАДАЧИ ТЕОРИИ МАССОВОГО ОБСЛУЖИВАНИЯ
  • 2. КЛАССИФИКАЦИЯ СИСТЕМ МАССОВОГО ОБСЛУЖИВАНИЯ И ИХ ОСНОВНЫЕ ХАРАКТЕРИСТИКИ
  • 3. ОДНОКАНАЛЬНАЯ СМО С ОТКАЗАМИ
  • 4. МНОГОКАНАЛЬНАЯ СМО С ОТКАЗАМИ
  • 5. ОДНОКАНАЛЬНАЯ СМО С ОЖИДАНИЕМ
  • 6. МНОГОКАНАЛЬНАЯ СМО С ОЖИДАНИЕМ
  • 7. СМО С ОГРАНИЧЕННЫМ ВРЕМЕНЕМ ОЖИДАНИЯ
  • 8. ЗАМКНУТЫЕ СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ
  • 9. СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ СО «ВЗАИМОПОМОЩЬЮ» МЕЖДУ КАНАЛАМИ
  • 1. СМО с отказами
  • 2. СМО с очередью
  • 10. СИСТЕМА МАССОВОГО ОБСЛУЖИВАНИЯ С ОШИБКАМИ
  • 11. СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ С НЕ-ПУАССОНОВСКИМИ ПОТОКАМИ СОБЫТИЙ
  • 2. Одноканальная СМО с ожиданием
  • 6. МЕТОД ДИНАМИКИ СРЕДНИХ
  • 1. ИДЕЯ МЕТОДА. ОБЛАСТЬ ПРИМЕНИМОСТИ
  • 2. УЧЕТ ЗАВИСИМОСТИ ИНТЕНСИВНОСТЕЙ ПОТОКОВ СОБЫТИЙ ОТ ЧИСЛЕННОСТЕЙ СОСТОЯНИИ. ПРИНЦИП КВАЗИРЕГУЛЯРНОСТИ
  • 3. УЧЕТ ПОПОЛНЕНИЯ ЧИСЛЕННОСТЕЙ СОСТОЯНИЙ
  • 4. МЕТОД ДИНАМИКИ СРЕДНИХ ДЛЯ СИСТЕМЫ, СОСТОЯЩЕЙ ИЗ НЕОДНОРОДНЫХ ЭЛЕМЕНТОВ
  • 5. ПРЕДЕЛЬНОЕ ПОВЕДЕНИЕ СРЕДНИХ ЧИСЛЕННОСТЕЙ СОСТОЯНИЙ
  • 6. УРАВНЕНИЯ ДИНАМИКИ БОЯ (МОДЕЛЬ А)
  • 7. УЧЕТ ПОПОЛНЕНИЯ СИЛ, УПРЕЖДАЮЩЕГО УДАРА, ТЕМПА МОБИЛИЗАЦИИ И ПРОЧИХ ФАКТОРОВ В УРАВНЕНИЯХ ДИНАМИКИ БОЯ
  • 8. МОДЕЛЬ Б. СЛУЧАЙ ОТСУТСТВИЯ ПЕРЕНОСА ОГНЯ
  • 9. МОДЕЛЬ В. УЧЕТ ДЕЯТЕЛЬНОСТИ РАЗВЕДКИ И СИСТЕМЫ УПРАВЛЕНИЯ БОЕМ
  • 10. УЧЕТ ВОССТАНОВЛЕНИЯ ЕДИНИЦ В ХОДЕ БОЕВЫХ ДЕЙСТВИЙ
  • 11. УРАВНЕНИЯ ДИНАМИКИ БОЯ ДЛЯ НЕОДНОРОДНЫХ ЕДИНИЦ. ФУНКЦИИ РАСПРЕДЕЛЕНИЯ ОГНЯ
  • 12. УРАВНЕНИЯ СМЕШАННОГО ТИПА
  • 13. НЕКОТОРЫЕ УТОЧНЕНИЯ МЕТОДА ДИНАМИКИ СРЕДНИХ
  • 7. МЕТОДЫ УЧЁТА НАДЕЖНОСТИ ТЕХНИЧЕСКИХ УСТРОЙСТВ
  • 2. НАДЕЖНОСТЬ ЭЛЕМЕНТА. ПЛОТНОСТЬ РАСПРЕДЕЛЕНИЯ ВРЕМЕНИ БЕЗОТКАЗНОЙ РАБОТЫ. СРЕДНЕЕ ВРЕМЯ БЕЗОТКАЗНОЙ РАБОТЫ
  • 3. ЭКСПОНЕНЦИАЛЬНЫЙ ЗАКОН НАДЕЖНОСТИ. ИНТЕНСИВНОСТЬ ОТКАЗОВ
  • 4. ОПРЕДЕЛЕНИЕ НАДЕЖНОСТИ СИСТЕМЫ ПО НАДЕЖНОСТИ ЕЕ ЭЛЕМЕНТОВ. НАДЕЖНОСТЬ НЕРЕЗЕРВИРОВАННОЙ СИСТЕМЫ
  • 5. НАДЕЖНОСТЬ РЕЗЕРВИРОВАННОЙ СИСТЕМЫ («ГОРЯЧИЙ РЕЗЕРВ»)
  • 6. НАДЕЖНОСТЬ РЕЗЕРВИРОВАННОЙ СИСТЕМЫ («ХОЛОДНЫЙ» И «ОБЛЕГЧЕННЫЙ» РЕЗЕРВ)
  • 7. НАДЕЖНОСТЬ СИСТЕМЫ С ВОССТАНОВЛЕНИЕМ
  • 8. УЧЕТ ЗАВИСИМОСТИ ОТКАЗОВ ПРИ ОЦЕНКЕ НАДЕЖНОСТИ ТЕХНИЧЕСКИХ УСТРОЙСТВ
  • 8. МОДЕЛИРОВАНИЕ ОПЕРАЦИЙ МЕТОДОМ СТАТИСТИЧЕСКИХ ИСПЫТАНИЙ
  • 2. ЕДИНИЧНЫЙ ЖРЕБИЙ
  • 3. РОЗЫГРЫШ ЗНАЧЕНИЯ НОРМАЛЬНО РАСПРЕДЕЛЕННОЙ СЛУЧАЙНОЙ ВЕЛИЧИНЫ
  • 4. ПОЛУЧЕНИЕ СЛУЧАЙНОГО ЧИСЛА R ОТ 0 ДО 1
  • 5. ПРИМЕРЫ МОДЕЛИРОВАНИЯ СЛУЧАЙНЫХ ПРОЦЕССОВ МЕТОДОМ МОНТЕ-КАРЛО
  • 6. ОПРЕДЕЛЕНИЕ ХАРАКТЕРИСТИК СТАЦИОНАРНОГО СЛУЧАЙНОГО ПРОЦЕССА МЕТОДОМ МОНТЕ-КАРЛО ПО ОДНОЙ РЕАЛИЗАЦИИ
  • 7. ОЦЕНКА ТОЧНОСТИ ХАРАКТЕРИСТИК, ПОЛУЧЕННЫХ МЕТОДОМ МОНТЕ-КАРЛО. НЕОБХОДИМОЕ ЧИСЛО РЕАЛИЗАЦИЙ
  • 9. ИГРОВЫЕ МЕТОДЫ ОБОСНОВАНИЯ РЕШЕНИЙ
  • 1. ЗАДАЧИ ТЕОРИИ ИГР И СТАТИСТИЧЕСКИХ РЕШЕНИЙ
  • 2. ПРЕДМЕТ ТЕОРИИ ИГР. ОСНОВНЫЕ ПОНЯТИЯ
  • 3. ПЛАТЕЖНАЯ МАТРИЦА
  • 4. НИЖНЯЯ И ВЕРХНЯЯ ЦЕНА ИГРЫ. ПРИНЦИП МИНИМАКСА
  • 5. РЕШЕНИЕ ИГРЫ В СМЕШАННЫХ СТРАТЕГИЯХ
  • 6. УПРОЩЕНИЕ ИГР
  • 7. ИГРА 2Х2
  • 8. ИГРЫ 2xn И mx2
  • 9. РЕШЕНИЕ ИГР mxn
  • 10. РЕШЕНИЕ КОНЕЧНЫХ ИГР МЕТОДОМ ИТЕРАЦИЙ
  • 11. ФИЗИЧЕСКАЯ СМЕСЬ СТРАТЕГИЙ
  • 12. ЭЛЕМЕНТЫ ТЕОРИИ СТАТИСТИЧЕСКИХ РЕШЕНИЙ
  • 13. КРИТЕРИЙ, ОСНОВАННЫЙ НА ИЗВЕСТНЫХ ВЕРОЯТНОСТЯХ УСЛОВИЙ. КРИТЕРИИ ВАЛЬДА, ГУРВИЦА, СЭВИДЖА
  • 14. ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТА В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ
  • 10. МЕТОД СЕТЕВОГО ПЛАНИРОВАНИЯ
  • 2. СЕТЕВОЙ ГРАФИК КОМПЛЕКСА РАБОТ. ВРЕМЕННОЙ СЕТЕВОЙ ГРАФИК
  • 3. ФОРМАЛЬНАЯ ЗАПИСЬ (АЛГОРИТМ) ЗАДАЧИ СЕТЕВОГО ПЛАНИРОВАНИЯ
  • 4. ОПТИМИЗАЦИЯ ПЛАНА КОМПЛЕКСА РАБОТ
  • 5. СЕТЕВОЕ «ПЛАНИРОВАНИЕ ПРИ СЛУЧАЙНЫХ ВРЕМЕНАХ ВЫПОЛНЕНИЯ РАБОТ. ПРИМЕНЕНИЕ ЭЦВМ
  • ПРИЛОЖЕНИЕ
  • ЛИТЕРАТУРА

Добавить комментарий