Рассмотрим четырехканальную систему массового обслуживания (n = 3) с максимальной длиной очереди равной трем (m = 2). В СМО поступает простейший поток заявок со средней интенсивностью л = 4.0 и показательным законом распределения времени между поступлением заявок. Поток обслуживаемых в системе заявок является простейшим со средней интенсивностью м = 1.0 и показательным законом распределения временем обслуживания.
Данная система имеет 9 состояний, обозначим их:
S0 - все каналы пусты, очередь пуста;
S1 - 1 канал занят, очередь пуста;
S2 - 2 канала заняты, очередь пуста;
S3 - 3 канала заняты, очередь пуста;
S4 - 3 канала заняты, в очереди 1 заявка;
S5 - 3 канала заняты, в очереди 2 заявки;
Вероятности прихода системы в состояния S0, S1, S2, …, S5 соответственно равны Р0, Р1, Р2, …, Р5.
Граф состояний системы массового обслуживания представляет собой схему гибели и размножения. Все состояния системы можно представить в виде цепочки, в которой каждое из состояний связано с предыдущим и последующим.
Рис. 3
Для построенного графа запишем уравнения Колмогорова:
Чтобы решить данную систему зададим начальные условия:
Систему уравнений Колмогорова (систему дифференциальных уравнений) решим численным методом Эйлера с помощью программного пакета Maple 8 (см. Приложение 1).
Метод Эйлера
где- в нашем случае, это правые части уравнений Колмогорова, n=7.
(1)
Выберем шаг по времени . Предположим , где Т - это время, за которое система выходит на стационарный режим. Отсюда получаем число шагов
.
Последовательно N раз вычисляя по формуле (1) получим зависимости вероятностей состояний системы от времени, приведенной на рис.4. Очевидно, что уже при система выходит на стационарный режим. Значения вероятностей СМО при равны:
Зависимости вероятностей состояний системы от времени
Рис. 4
2.2 Финальные вероятности системы
При достаточно большом времени протекания процессов в системе () могут устанавливаться вероятности состояний, не зависящие от времени, которые называются финальными вероятностями, т.е. в системе устанавливается стационарный режим. Если число состояний системы конечно, и из каждого из них за конечное число шагов можно перейти в любое другое состояние, то финальные вероятности существуют, т.е.
Т.к. в стационарном состоянии производные по времени равны 0, то уравнения для финальных вероятностей получаются из уравнений Колмогорова путем приравнивания правых частей 0. Запишем уравнения для финальных вероятностей для нашей СМО.
Решим данную систему линейных уравнений с помощью программного пакета Maple 8 (см. Приложение 1).
Получим финальные вероятности системы:
Сравнение вероятностей, полученных из системы уравнений Колмогорова при , с финальными вероятностями показывает, что ошибки
равны:
Т.е. достаточно малы. При увеличении интервала времени до , погрешности должны были стать еще меньше.
Это подтверждает правильность полученных результатов.
2.3 Расчет показатели эффективности системы по финальным вероятностям
Найдем показатели эффективности системы массового обслуживания. Наиболее важными являются следующие показатели:
1) Вероятность отказа в обслуживании заявки, т.е. вероятность того, что заявка покидает систему не обслуженной. В нашем случае заявке отказывается в обслуживании, если все 4 канала заняты, и очередь максимально заполнена (т.е. 3 человека в очереди), это соответствует состоянию системы S7. Т.к. вероятность прихода системы в состояние S7 равна Р7 , то
2) Относительная пропускная способность - это средняя доля поступивших заявок, обслуживаемых системой.
3) Абсолютная пропускная способность - это среднее число заявок, обслуживаемых в единицу времени.
4) Длина очереди, т.е. среднее число заявок в очереди. Длина очереди равна сумме произведений числа человек в очереди на вероятность соответствующего состояния.
5) Среднее время пребывания заявки в очереди определяется формулой Литтла
6) Среднее число занятых каналов определяется следующим образом:
Глава 3. Имитационное моделирование СМО
3.1 Алгоритм метода имитационного моделирования СМО (пошаговый подход)
Рассмотрим четырехканальную систему массового обслуживания (n = 3) с максимальной длиной очереди равной четырем (m = 2). В СМО поступает простейший поток заявок со средней интенсивностью л = 4 и показательным законом распределения времени между поступлением заявок. Поток обслуживаемых в системе заявок является простейшим со средней интенсивностью м = 1 и показательным законом распределения временем обслуживания.
Для имитации СМО воспользуемся одним из методов статистического моделирования - имитационным моделированием. Будем использовать пошаговый подход. Суть этого подхода в том, что состояния системы рассматриваются в последующие моменты времени, шаг между которыми является достаточно малым, чтобы за его время произошло не более одного события.
Выберем шаг по времени (). Он должен быть много меньше среднего времени поступления заявки () и среднего времени ее обслуживания (), т.е.
, где (3.1.1)
Исходя из условия (3.1.1) определим шаг по времени .
Время поступления заявки в СМО и время ее обслуживания являются случайными величинами. Поэтому, при имитационном моделировании СМО их вычисление производится с помощью случайных чисел.
Рассмотрим поступление заявки в СМО. Вероятность того, что на интервале в СМО поступит заявка, равна:
Сгенерируем случайное число , и, если , то будем считать, что заявка на данном шаге в систему поступила, если , то не поступила.
В программе это осуществляет функция . Интервал времени примем постоянным и равным 0,001, тогда отношение будет равно 1000 . Если заявка поступила, то она принимает значение «истина», в противном случае значение «ложь».
function sob: boolean;
var r:real;
begin
r := Random(1000)/1000;
if r <= (i_deltaT*r_Lamda) then
sob:= true
else
sob:= false;
end;
Рассмотрим теперь обслуживание заявки в СМО. Время обслуживания заявки в системе определяется выражением
,
где - случайное число. В программе время обслуживания определяется с помощью функции .
function tok: real;
t_ob: real;
r:= Random(1000)/1000;
t_ob:= -1/r_Mu*ln(1-r);
tok:= t_ob;
Алгоритм метода имитационного моделирования можно сформулировать следующим образом. Время работы СМО (Т) разбивается на шаги по времени dt, на каждом из них выполняется ряд действий. Вначале определяются состояния системы (занятость каналов, длина очереди), затем, с помощью функции , определяется, поступила ли на данном шаге заявка или нет.
Если поступила, и, при этом имеются свободные каналы, то с помощью функции генерируем время обработки заявки и ставим ее на обслуживание. Если все каналы заняты, а длина очереди меньше 4, то помещаем заявку в очередь, если же длина очереди равна 4, то заявке будет отказано в обслуживании.
В случае, когда на данном шаге заявка не поступала, а канал обслуживания освободился, проверяем, есть ли очередь. Если есть, то из очереди заявку ставим на обслуживание в свободный канал. После проделанных операций время обслуживания для занятых каналов уменьшаем на величину шага dt. По истечении времени Т, т.е., после моделирования работы СМО, вычисляются показатели эффективности работы системы и результаты выводятся на экран.
3.2 Блок-схема программы
Блок-схема программы, реализующей описанный алгоритм, приведена на рис.5. Распишем некоторые блоки более подробно.
Блок 1. Задание начальных значений параметров. Пользователем задаются значения: Time - время работы системы,
r_la - интенсивность потока поступления заявок (),
r_mu - интенсивность потока обслуживания заявок ().
Значения, задаваемые программой:
i_deltaT:=0.001 - шаг по времени;
i_post:=0;
i_otk:= 0;
i_obsl:=0;
i_och:=0;
i_tobs:=0;
t:=0 - параметр, следящий, сколько времени проработала система;
i_obsl:=0 - число обслуженных заявок;
i_tobs:=0 - общее время обслуживания заявок в системе;
i_post:=0 - число поступивших в СМО заявок;
i_otk:=0 - число заявок, которым было отказано в обслуживании;
sost:='s1' - состояние системы (изначально устанавливаем СМО в состояние s1, т.е. все каналы свободны);
for i :=1 to 3 do
t_och[i] := 0;- обнуляем время пребывания СМО в состояниях с длиной очереди 1, 2.
t_okonch[i] := 0; - время окончания обслуживания заявки во всех 3 каналах устанавливаем в 0, т.к. каналы пусты;
Рис.5.
Блок 3. Задание состояний системы. Выделим у данной 3-х канальной системы 9 различных состояний: s1, s2,.. s9. СМО находится в состоянии s1, когда система свободна; s2..s7 - хотя бы один канал свободен; в состоянии s8, когда все каналы заняты, и есть место в очереди; в состоянии s9 - все каналы заняты, и очередь достигла максимальной длины (och=2).
Задаем состояния системы:
_sost:=1;
for i:=0 to 3 do
if t_okonch[i+1]>0 then i_temp:= 1
else i_temp:= 0;
_sost:=_sost + i_temp* round(Exp(i*ln(2)));
if _sost<7 then begin Sost:= _sost; exit; end;
if (_sost= 7) and (i_och< 2) then Sost:= _sost
Страницы: 1, 2, 3, 4, 5, 6
