Курсовая работа: Моделирование систем массового обслуживания
Для описания процессов с непрерывным временем используется модель в виде так называемой марковской цепи с дискретными состояниями системы, или непрерывной марковской цепью.
Глава II . Уравнения описывающие системы массового обслуживания
2.1 Уравнения Колмогорова
Рассмотрим математическое описание марковского случайного процесса с дискретными состояниями системы So , Sl , S2 (см. рис. 6.2.1) и непрерывным временем. Полагаем, что все переходы системы массового обслуживания из состояния Si в состояние Sjпроисходят под воздействием простейших потоков событий с интенсивностями λij , а обратный переход под воздействием другого потока λij ,. Введем обозначение pi как вероятность того, что в момент времени t система находится в состоянии Si . Для любого момента времени tсправедливо записать нормировочное условие—сумма вероятностей всех состояний равна 1:
2
Σpi (t)=p0 (t)+ p1 (t)+ p2 (t)=1
i =0
Проведем анализ системы в момент времени t, задав малое приращение времени Δt, и найдем вероятность р1 (t+ Δt) того, что система в момент времени (t+ Δt) будет находиться в состоянии S1 которое достигается разными вариантами:
а) система в момент t с вероятностью p1 (t) находилась в состоянии S1 и за малое приращение времени Δt так и не перешла в другое соседнее состояние — ни в S0 , ни bS2 . Вывести систему из состояния S1 можно суммарным простейшим потоком c интенсивностью (λ10 +λ12 ), поскольку суперпозиция простейших потоков также является простейшим потоком. На этом основании вероятность выхода из состояния S1 за малый промежуток времени Δtприближенно равна (λ10 +λ12 )* Δt. Тогда вероятность невыхода из этого состояния равна [1 -(λ10 +λ12 )* Δt].Bсоответствии с этим вероятность того, что система останется в состоянии Siна основании теоремы умножения вероятностей, равна:
p1 (t) [1 -(λ10 +λ12 )* Δt];
б)система находилась в соседнем состоянии So и за малое время Δt перешла в состояние So Переход системы происходит под воздействием потока λ01 с вероятностью, приближенно равной λ01 Δt
Вероятность того, что система будет находиться в состоянии S1 , в этом варианте равна po (t)λ01 Δt;
в) система находилась в состоянии S2 и за время Δt перешла в состояние S1 под воздействием потока интенсивностью λ21 с вероятностью, приближенно равной λ21 Δt. Вероятность того, что система будет находиться в состоянии S1 , равна p2 (t) λ21 Δt.
Применяя теорему сложения вероятностей для этих вариантов, получим выражение:
p2 (t+Δt)= p1 (t) [1 -(λ10 +λ12 )* Δt]+ po (t)λ01 Δt+p2 (t) λ21 Δt ,
которое можно записать иначе:
p2 (t+Δt)-p1 (t)/ Δt= po (t)λ01 + p2 (t) λ21 - p1 (t) (λ10 +λ12 ) .
Переходя к пределу при Δt-> 0, приближенные равенства перейдут в точные, и тогда получим производную первого порядка
dp2 /dt= p0 λ01 +p2 λ21 -p1 (λ10 +λ12 ) ,
что является дифференциальным уравнением.
Проводя рассуждения аналогичным образом для всех других состояний системы, получим систему дифференциальных уравнений, которые называются уравнениями А.Н. Колмогорова:
dp0 /dt= p1 λ10 ,
dp1 /dt= p0 λ01 +p2 λ21 -p1 (λ10 +λ12 ) ,
dp2 /dt= p1 λ12 +p2 λ21 .
Для составления уравнений Колмогорова существуют общие правила.
Уравнения Колмогорова позволяют вычислить все вероятности состояний СМО Si в функции времени pi (t). В теории случайных процессов показано, что если число состояний системы конечно, а из каждого из них можно перейти в любое другое состояние, то существуют предельные (финальные) вероятности состояний, которые показывают на среднюю относительную величину времени пребывания системы, в этом состоянии. Если предельная вероятность состояния S0 – равна p0 = 0,2, то, следовательно, в среднем 20% времени, или 1/5 рабочего времени, система находится в состоянии So . Например, при отсутствии заявок на обслуживание к = 0, р0 = 0,2,; следовательно, в среднем 2 ч в день система находится в состоянии So и простаивает, если продолжительность рабочего дня составляет 10 ч.
Поскольку предельные вероятности системы постоянны, то заменив в уравнениях Колмогорова соответствующие производные нулевыми значениями, получим систему линейных алгебраических уравнений, описывающих стационарный режим СМО. Такую систему уравнений составляют по размеченному графу состояний СМО по следующим правилам: слева от знака равенств?