Одноканальная смо с ожиданием и ограниченной очередью. Одноканальная смо с ожиданием Смо с m приборами и ограниченной очередью

Системы с ожиданием при неограниченном входящем потоке

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

Расчетные формулы
Вероятность того, что все каналы свободны

Вероятность того, что занято k каналов, при условии, что общее число заявок, находящихся на обслуживании, не превосходит числа каналов,

Вероятность того, что в системе находится k заявок, в случае, когда их число больше числа каналов,

Вероятность того, что все каналы заняты,

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

Средняя длина очереди

Среднее число свободных от обслуживания каналов

Пример
Автозаправочная станция с двумя колонками обслуживает пуассоновский поток машин с интенсивностью λ=0,8 машин в минуту. Время обслуживания одной машины подчиняется показательному закону со средним значением 2 минуты. В данном районе нет другой АЗС, так что очередь перед АЗС может расти практически неограниченно. Найдите:
1) среднее число занятых колонок;
2) вероятность отсутствия очереди у АЗС;
3) вероятность того, что придется ждать начала обслуживания;
4) среднее число машин в очереди;
5) среднее время ожидания в очереди;
6) среднее время пребывания машины на АЗС;
7) среднее число машин на АЗС.
Решение . По условию задачи n=2, λ=0.8; μ=1/t обсл =0.5; ρ=λ/μ=1.6
Поскольку ρ /n =0,8<1, то очередь не растет безгранично и имеет смысл говорить о предельном стационарном режиме работы системы массового обслуживания.
Находим вероятности состояний СМО:

Среднее число занятых колонок:
N зан =n-N 0 = 2-(2·p 0 +1·p 1) = 2-2·0.1111 - 0.1778 = 1.6
Вероятность отсутствия очереди у АЗС:

Вероятность того, что придется ждать начала обслуживания равна вероятности того, что все колонки заняты:
p 0 +p 1 +p 2 = 0.1111+0.1778+0.1422 = 0.4311
Среднее число машин в очереди:

Среднее время ожидания в очереди:
Среднее время пребывания машины на АЗС:
t преб =t обсл +t ож = 2+3.5556 = 5.5556 мин.
Среднее число машин на АЗС:
N зан +L оч = 1.6+2.8444 = 4.4444
Рассмотрим одноканальную СМО с ожиданиями, в которой число каналов равно единице n = 1, интенсивность поступления заявок – λ, интенсивность обслуживания равна μ. Заявка, поступившая в тот момент времени, когда канал занят, становится в очередь и ждет обслуживания. Количество мест в очереди ограничено и равно m . Если все места в очереди заняты, то заявка покидает очередь не обслуженной. Проанализируем состояние системы:

S 0 – канал свободен;
S 1 – канал занят;
S 2 – канал занят, одна заявка в очереди;
S k – канал занят, (k–1) заявок в очереди;
S m + 1 – канал занят, в очереди m заявок.

Изобразим граф состояний такой СМО (рис. 25).

Рис. 25
По формулам Эрланга найдем вероятности событий, состоящих в том, что СМО находится в состоянии S 1 , S 2 , …, S m+1:
(28)

При этом вероятность того, что заявка, прибывшая в систему, найдет ее свободной, равна
. (29)
Отношение интенсивности поступления заявок λ к интенсивности обслуживания заявок μ есть приведенная интенсивность μ, т.е.

ρ=λ/μ
Произведем замену в формулах (28) и (29) отношения λ/&mu на ρ, тогда выражения примут вид:

(30)
Вероятность Р 0 будет вычисляться по следующей формуле:
p 0 = -1 . (31)
Выражение для вероятности P 0 есть геометрическая прогрессия, сумма которой будет равна

.
Таким образом, формулы (30) и (31) позволяют определить вероятность любого события, которое может произойти в системе, т. е. определить вероятность нахождения системы в любом состоянии.
Формула для P 0 справедлива для случая, когда ρ ≠ 1 . В случае, когда ρ = 1 , т. е. интенсивность поступления заявок равна интенсивности их обслуживания, используется другая формула для вычисления вероятности того, что система свободна:

,
где m – это количество заявок, находящихся в очереди.

Определим характеристики эффективности одноканальной СМО :

вероятность того, что очередная заявка, прибывшая в систему, получит отказ Р отк;
абсолютную пропускную способность А ,
относительную пропускную способность Q ,
число занятых каналов k ,
среднее число заявок в очереди r ,
среднее число заявок, связанных с СМО, z .

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

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

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

A=Q·λ (34)
Таким образом, по формулам (32), (33), (34) можно вычислить основные показатели эффективности для любой системы массового обслуживания. теперь выведем выражения для вычисления характеристик, присущих лишь данной СМО.
Среднее число заявок в очереди r определим как математическое ожидание дискретной случайной величины, где R – число заявок в очереди.
♦ Р 2 – это вероятность того, что в очереди на обслуживание находится одна заявка;
♦ Р 3 – вероятность того, что в очереди две заявки;
♦ Р k – вероятность того, что в очереди (k–1) заявка;
♦ Р m + 1 – вероятность того что в очереди m заявок.
Тогда среднее число заявок в очереди можно вычислить следующим образом:
r =1·P 2 +2·P 3 + ... +(k-1)·P k + ... +m·P m+1 . (35)
Подставим в формулу (35) найденные ранее значения вероятностей, вычисленные в формуле (30):
r =1·ρ 2 ·p 0 +2·ρ 3 ·p 0 + ... +(k-1)·ρ k ·p 0 + ... +m·ρ m+1 ·p 0 . (35)
Вынесем за скобку вероятность P 0 и Р 2 , тогда получим итоговую формулу для вычисления среднего числа заявок в очереди на обслуживание:
r =ρ 2 ·p 0 (1+2·ρ+ ... +(k-1)·ρ k-2 + ... +m·ρ m-1)
Выведем формулу для среднего числа заявок, связанных с СМО, z , т. е. число заявок в очереди, находящихся на обслуживании. Рассмотрим общее число заявок, связанных с СМО, z как сумму двух величин среднего числа заявок в очереди r и числа занятых каналов k :

z = r +k .
Так как канал один, то число занятых каналов k может принимать значения 0 или 1. Вероятность того, что k = 0, т.е. система свободна, соответствует вероятности Р 0 , значение которой можно найти по формуле (31). Если k = 1, т.е. канал занят обслуживанием заявки, но места в очереди еще есть, то вероятность этого события можно вычислить по формуле

.
Следовательно, z будет равно:

. (37)

Одноканальная СМО с ожиданием

Система массового обслуживания имеет один канал. Входящий поток заявок на обслуживание - простейший поток с интенсивностью l. Интенсивность потока обслуживания равна m (т. е. в среднем непрерывно занятый канал будет выдавать m. обслуженных заявок). Длительность обслуживания - случайная величина, подчиненная показательному закону распределения. Поток обслуживании является простейшим пуассоновским потоком событий. Заявка, поступившая в момент, когда канал занят, становится в очередь и ожидает обслуживания.
Предположим, что независимо от того, сколько требований поступает на вход обслуживающей системы, данная система (очередь + обслуживаемые клиенты) не может вместить более N-требований (заявок), т. е. клиенты, не попавшие в ожидание, вынуждены обслуживаться в другом месте. Наконец, источник, порождающий заявки на обслуживание, имеет неограниченную (бесконечно большую) емкость.
Граф состояний СМО в этом случае имеет вид, показанный на Рис. 3.2.

Граф состояний одноканальной СМО с ожиданием (схема гибели и размножения)
Состояния СМО имеют следующую интерпретацию:
S 0 - канал свободен
S 1 - канал занят (очереди нет);
S 2 - канал занят (одна заявка стоит в очереди);
………………………………
S n - канал занят (n - 1 заявок стоит в очереди);
……………………………
S N - канал занят (N - 1 заявок стоит в очереди).
Стационарный провес в данной системе будет описываться следующей системой алгебраических уравнений :

п - номер состояния.
Решение приведенной выше системы уравнений (3.10) для нашей модели СМО имеет вид

Следует отметить, что выполнение условия стационарности для данной СМО необязательно, поскольку число допускаемых в обслуживающую систему заявок контролируется путем введения ограничения на длину очереди (которая не может превышать N - 1), а не соотношением между интенсивностями входного потока, т. е. не отношением
l/m = p
Определим характеристики одноканальной СМО с ожиданием и ограниченной длиной очереди, равной (N - 1):

Рассмотрим пример одноканальной СМО с ожиданием.
Пример 3.2. Специализированный пост диагностики представляет собой одноканальную СМО. Число стоянок для автомобилей, ожидающих проведения диагностики, ограничено и равно 3 [(N - 1) = 3]. Если все стоянки заняты, т. е. в очереди уже находится три автомобиля, то очередной автомобиль, прибывший на диагностику, в очередь на обслуживание не становится. Поток автомобилей, прибывающих на диагностику, распределен по закону Пуассона и имеет интенсивность l = 0,85 (автомобиля в час). Время диагностики автомобиля распределено по показательному закону и в среднем равно 1,05 час.
Требуется определить вероятностные характеристики поста диагностики, работающего в стационарном режиме.
Решение
1. Параметр потока обслуживании автомобилей:

2. Приведенная интенсивность потока автомобилей определяется как отношение интенсивностей l и m, т. е.

3. Вычислим финальные вероятности системы:

P 1 =ρ·P 0 = 0.893·0.248 = 0.221
P 2 =ρ 2 ·P 0 = 0.893 2 ·0.248 = 0.198
P 3 =ρ 3 ·P 0 = 0.893 3 ·0.248 = 0.177
P 4 =ρ 4 ·P 0 = 0.893 2 ·0.248 = 0.158
4. Вероятность отказа в обслуживании автомобиля:
P отк =P 4 =ρ 4 ·P 0 ≈ 0.158
5. Относительная пропускная способность поста диагностики:
q=1-P отк = 1-0.158 = 0.842
6. Абсолютная пропускная способность поста диагностики
A=λ·q = 0.85·0.842 = 0.716 (автомобиля в час)
7. Среднее число автомобилей, находящихся на обслуживании и в очереди (т.е. в системе массового обслуживания):

8. Среднее время пребывания автомобиля в системе:
9. Средняя продолжительность пребывания заявки в очереди на обслуживание:
W q =W S -1/μ = 2.473-1/0.952 = 1.423 часа
10. Среднее число заявок в очереди (длина очереди): L q = А,(1 - P N) W q = 0,85
L q =λ(1-P N)·W q = 0.85·(1-0.158)·1.423 = 1.02
Работу рассмотренного поста диагностики можно считать удовлетворительной, так как пост диагностики не обслуживает автомобили в среднем в 15,8% случаев (Р отк = 0,158). В качестве показателей эффективности СМО с ожиданием, кроме уже известных показателей - абсолютной А и относительной Q пропускной способности, вероятности отказа P отк. , среднего числа занятых каналов (для многоканальной системы) будем рассматривать также следующие: L сист. - среднее число заявок системе; Т сист. - среднее время пребывания заявки в системе; L оч. - среднее число заявок в очереди (длина очереди); Т оч. - среднее время пребывания заявки в очереди; Р зан.. - вероятность того, что канал занят (степень загрузки канала).

Одноканальная система с неограниченной очередью

На практике часто встречаются одноканальные СМО с неограниченной очередью (например, телефон-автомат с одной будкой).
Рассмотрим задачу.
Имеется одноканальная СМО с очередью, на которую не наложены никакие ограничения (ни по длине очереди, ни по времени ожидания). Поток заявок, поступающих в СМО, имеет интенсивность λ, а поток обслуживании - интенсивность μ. Необходимо найти предельные вероятности состояний и показатели эффективности СМО.
Система может находиться в одном из состояний S 0 , S 1 , S 2 , …, S k , по числу заявок, находящихся в СМО: S 0 - канал свободен; S 1 - канал занят (обслуживает заявку), очереди нет, S 2 - канал занят, одна заявка стоит в очереди; ... S k - канал занят, (k-1) заявок стоят в очереди и т.д.
Граф состояний СМО представлен на рис. 8.

Рис. 8
Это процесс гибели и размножения, но с бесконечным числом состояний, в котором интенсивность потока заявок равна λ, а интенсивность потока обслуживании μ.
Прежде чем записать формулы предельных вероятностей, необходимо быть уверенным в их существовании, ведь в случае, когда время t→∞, очередь может неограниченно возрастать. Доказано, что если ρ<1, т.е. среднее число приходящих заявок меньше среднего числа обслуженных заявок (в единицу времени), то предельные вероятности существуют. Если ρ≥1, очередь растет до бесконечности.

Для определения предельных вероятностей состояний воспользуемся формулами (16), (17) для процесса гибели и размножении (здесь мы допускаем известную нестрогость, так как ранее эти формулы были получены для случая конечного числа состояний системы). Получим(32)
Так как предельные вероятности существуют лишь при ρ < 1, то геометрический ряд со знаменателем
ρ < 1, записанный в скобках в формуле (32), сходится к сумме, равной . Поэтому
p 0 =1-ρ, (33)
и с учетом соотношений (17)
p 1 =ρ·p 0 ; p 2 =ρ 2 ·p 0 ; ... ; p k =ρ k ·p 0 ; ...
найдем предельные вероятности других состояний
p 1 =ρ·(1-ρ); p 2 =ρ 2 ·(1-ρ); ... ; p k =ρ k ·(1-ρ); ... (34)
Предельные вероятности p 0 , p 1 , p 2 , …, p k ,… образуют убывающую геометрическую профессию со знаменателем р < 1, следовательно, вероятность р 0 - наибольшая. Это означает, что если СМО справляется с потоком заявок (при ρ < 1), то наиболее вероятным будет отсутствие заявок в системе.
Среднее число заявок в системе L сист. определим по формуле математического ожидания, которая с учетом (34) примет вид
(35)
(суммирование от 1 до ∞, так как нулевой член 0·p 0 =0).
Можно показать, что формула (35) преобразуется (при ρ < 1) к виду
(36)
Найдем среднее число заявок в очереди L оч. Очевидно, что
L оч =L сист -L об (37)
где L об. - среднее число заявок, находящихся под обслуживанием.
Среднее число заявок под обслуживанием определим по формуле математического ожидания числа заявок под обслуживанием, принимающего значения 0 (если канал свободен) либо 1 (если канал занят):
L оч =0·p 0 +1·(1-p 0)
т.е. среднее число заявок под обслуживанием равно вероятности того, что канал занят:
L оч =P зан =1-p 0 , (38)
В силу (33)
L оч =P зан ρ, (39)
Теперь по формуле (37) с учетом (36) и (39)
(40)
Доказано, что при любом характере потока заявок, при любом распределении времени обслуживания, при любой дисциплине обслуживания среднее время пребывания заявки в системе (очереди) равна среднему числу заявок в системе (в очереди), деленному на интенсивность потока заявок, т.е.
(41)
(42)
Формулы (41) и (42) называются формулами Литтла. Они вытекают из того, что в предельном, стационарном режиме среднее число заявок, прибывающих в систему, равно среднему числу заявок, покидающих ее: оба потока заявок имеют одну и ту же интенсивность λ.
На основании формул (41) и (42) с учетом (36) и (40) среднее время пребывания заявки в системе определится по формуле:
(43)
а среднее время пребывания заявки в очереди
(44)

Одноканальная СМО с ожиданием без ограничения на вместимость блока ожидания

Стационарный режим функционирования данной СМО существует при t→∞ для любого п=0,1,2,… и когда l < m.Система алгебраических уравнений, описывающих работу СМО при t®¥ для любого n = 0, 1, 2...., имеет вид
Решение данной системы уравнений имеет вид
P n =(1-ρ)·ρ n , n=0,1,2,... (3.21)
где ρ=λ/μ < 1
Характеристики одноканальной СМО с ожиданием, без ограничения на длину очереди, следующие:
среднее число находящихся в системе клиентов (заявок) на обслуживание:
средняя продолжительность пребывания клиента в системе:

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

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

Решение
1. Параметр потока обслуживания m и приведенная интенсивность потока автомобилей р определены в примере 3.2:
m = 0,952; p = 0,893.
2. Вычислим предельные вероятности системы по формулам
P 0 =1-ρ = 1-0.893 = 0.107
P 1 =(1-ρ)·ρ = (1-0.893)·0.893 = 0.096
P 2 =(1-ρ)·ρ 2 = (1-0.893) 2 ·0.893 = 0.085
P 3 =(1-ρ)·ρ 3 = (1-0.893) 3 ·0.893 = 0.076
P 4 =(1-ρ)·ρ 4 = (1-0.893) 4 ·0.893 = 0.068
P 5 =(1-ρ)·ρ 5 = (1-0.893) 5 ·0.893 = 0.061
и т.д.
Следует отметить, что Р о определяет долю времени, в течение которого пост диагностики вынужденно бездействует (простаивает). В нашем примере она составляет 10,7%, так как Р о = 0,107.
3. Среднее число автомобилей, находящихся в системе (на обслуживании и в очереди):
4. Средняя продолжительность пребывания клиента в системе:

6. Средняя продолжительность пребывания автомобиля в очереди-
7. Относительная пропускная способность системы:
т. е. каждая заявка, пришедшая в систему, будет обслужена.
8. Абсолютная пропускная способность: А = lq = 0,85·1 = 0,85
Следует отметить, что предприятие, осуществляющее диагностику автомобилей, прежде всего интересует количество клиентов, которое посетит пост диагностики при снятии ограничения на длину очереди.
Допустим, в первоначальном варианте количество мест для стоянки прибывающих автомобилей было равно трем (см. пример 3.2). Частота m возникновения ситуаций, когда прибывающий на пост диагностики автомобиль не имеет возможности присоединиться к очереди:

т = l P N

В нашем примере при N = 3 + 1 = 4 и р = 0,893,
m = l Р о р 4 = 0,85·0,248·0,8934·0,134 автомобиля в час.
При 12-часовом режиме работы поста диагностики это эквивалентно тому, что пост диагностики в среднем за смену (день) будет терять 12·0,134 = 1,6 автомобиля.
Снятие ограничения на длину очереди позволяет увеличить количество обслуженных клиентов в нашем примере в среднем на 1,6 автомобиля за смену (12 ч. работы) поста диагностики. Ясно, что решение относительно расширения площади для стоянки автомобилей, прибывающих на пост диагностики, должно основываться на оценке экономического ущерба, который обусловлен потерей клиентов при наличии всего трех мест для стоянки этих автомобилей.

Многоканальная СМО с неограниченной очередью

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

Система может находиться в одном из состояний S 0 , S 1 , S 2 ,…, S k ,…, S n ,…, - нумеруемых по числу заявок, находящихся в СМО: S 0 - в системе нет заявок (все каналы свободны); S 1 - занят один канал, остальные свободны; S 2 - заняты два канала, остальные свободны;..., S k - занято k каналов, остальные свободны;..., S n - заняты все n каналов (очереди нет); S n+1 - заняты все n каналов, в очереди одна заявка;..., S n+r - заняты все n каналов, r заявок стоит в очереди,....

Граф состояний системы показан на рис. 9. Обратим внимание на то, что в отличие от предыдущей СМО, интенсивность потока обслуживаний (переводящего систему из одного состояния в другое справа налево) не остается постоянной, а по мере увеличения числа заявок в СМО от 0 до n увеличивается от величины m до nm, так как соответственно увеличивается число каналов обслуживания. При числе заявок в СМО большем, чем n, интенсивность потока обслуживании сохраняется равной nm.

среднее число заявок в очереди
, (50)
среднее число заявок в системе
L сист =L оч +ρ, (51)
Среднее время пребывания заявки в очереди и среднее время пребывания заявки в системе, как и ранее, находятся по формулам Литтла (42) и (41).
Замечание. Для СМО с неограниченной очередью при r < 1 любая заявка, пришедшая в систему, будет обслужена, т.е. вероятность отказа P отк = 0, относительная пропускная способность Q =1, а абсолютная пропускная способность равна интенсивности входящего потока заявок, т.е. А =l.

СМО с ограниченной очередью

СМО с ограниченной очередью. СМО с ограниченной очередью отличаются от рассмотренных выше задач лишь тем, что число заявок в очереди ограничено (не может превосходить некоторого заданного т). Если новая заявка поступает в момент, когда все места в очереди заняты, она покидает СМО необслуженной, т.е. получает отказ.
Очевидно: для вычисления предельных вероятностей состояний и показателей эффективности таких СМО может быть использован тот же подход, что и выше, с той разницей, что суммировать надо не бесконечную прогрессию (как, например, мы делали при выводе формулы (33)), а конечную.
Среднее время пребывания заявки в очереди и в системе, как и ранее, определяем по формулам Литтла (44) и (43).
СМО с ограниченным временем ожидания. На практике часто встречаются СМО с так называемыми "нетерпеливыми" заявками. Такие заявки могут уйти из очереди, если время ожидания превышает некоторую величину. В частности, такого рода заявки возникают в различных технологических системах, в которых задержка с началом обслуживания может привести к потере качества продукции, в системах оперативного управления, когда срочные сообщения теряют ценность (или даже смысл), если они не поступают на обслуживание в течение определенного времени.

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

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

Имеется n-канальная СМО с неограниченной очередью. Она характеризуется следующими показателями :

Предельные вероятности:

, , . . . , , ,…, ,… (10)

Вероятность того, что заявка окажется в очереди:

(11)

(13)

Среднее время нахождения в очереди:

(15)

Среднее время нахождения заявки в очереди:

Рассмотрим пример решения задачи многоканальной СМО с ожиданием.

Задача . В магазине к кассам поступает поток покупателей с интенсивностью 81 человек в час. Средняя продолжительность обслуживания кассиром одного покупателя tобсл = 2 мин. Определить предельные вероятности состояний и характеристики обслуживания узла расчета.

По условию λ=81(чел./час)= 81/60=1,35 (чел./мин.). По формулам (1, 2):

= λ/μ= λ * tобсл = 1,35 * 2 = 2,7

<1, т.е. при n > = 2,7. Таким образом, минимальное количество кассиров n =3.

Найдем характеристики обслуживания СМО при n=3.

Вероятность того, что в кассах отсутствуют покупатели, по формуле (9):

= (1+2,7+2,7 /2!+2,7 /3!+2,7 /3!(3-2,7)) = 0,025

В среднем 2,5 % времени кассиры будут простаивать.

Вероятность того, что в кассах будет очередь, определим по формуле (11):

P = (2,7 /3!(3-2,7))0,025 = 0,735

Среднее число покупателей, находящихся в очереди рассчитывается по формуле (13):

L = (2,7 /(3*3!(1-2,7/3) ))*0,025 = 7,35 (чел.)

T =7,35/1,35 = 5,44 (мин.)

Определим среднее число покупателей в кассах по формуле (15):

L =7,35+2,7=10,05 (чел.)

Среднее время нахождения покупателей в кассах находится по формуле (16):

T =10,05/1,35=7,44 (мин)

Среднее число кассиров, занятых обслуживанием покупателей, по формуле (12) =2,7.

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

Абсолютная пропускная способность узла расчета A=1,35 (чел./мин), или 81 (чел./час), т.е. 81 покупатель в час. Анализ характеристик обслуживания свидетельствует о значительной перегрузке касс при наличии трех кассиров.

Системы массового обслуживания с ограниченной очередью

Имеется n-канальная СМО с ограниченной очередью. Число заявок в очереди ограничено числом m. Если заявка поступает в момент, когда в очереди уже m заявок, она не обслуживается. Такая СМО характеризуется следующими показателями :

Предельные вероятности:

(17)

, , . . . , , ,…, (18)

Вероятность отказа:

(19)

Относительная пропускная способность:

Абсолютная пропускная способность:

Среднее число занятых каналов:

Среднее число заявок в очереди:

(23)

Среднее число заявок в системе:

Пример оптимизации СМО

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

Задача.

Определить оптимальное количество причалов в порту с минимальными затратами, если известно, что за год было обслужено 270 судов. Разгрузка одного судна длится в среднем 12 часов. Пеня за простой судна в порту составляет 100 тыс.р./сут.. Затраты на причал 150 тыс.р./сут. Расчеты приведены в таблице.

Решение.

По условию

λ=270(судов/год)=270/360=0,75(судов/сут.),

tобсл=12ч=12/24=0,5 сут.

По формулам (1, 2):

= λ/μ= λ * tобсл = 0,75 * 0,5 = 1,5

Очередь не будет возрастать до бесконечности при условии /n <1, т.е. при n > = 1,5. Таким образом, минимальное количество причалов n =2.

Найдем характеристики обслуживания СМО порта при количестве причалов n=2.

Вероятность того, что в порту отсутствуют суда, вычислим по формуле (9):

В среднем 1,4 % времени причалы будут простаивать.

Среднее число судов, находящихся в очереди рассчитывается по формуле (13):

Среднее время ожидания в очереди вычисляется по формуле (14):

T =1,93/0,75 = 2,57 (сут.)

Определим среднее число судов в порту по формуле (15):

L =1,93+1,5=3,43 (судна)

Среднее время нахождения судов в порту находится по формуле (16):

T =3,43 /0,75 =4,57 (сут)

Среднее число занятых причалов (12) =1,5.

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

Найдем суммарную пеню за простой судов в порту в сутки. Для этого перемножим пеню за простой судна в порту и среднее число судов в очереди:

= * L .

Определим затраты по обслуживанию причалов в сутки: = *n.

Для двух причалов в сутки

Суммарные затраты составят: С= + =193+300=493(ден.ед.)

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

Рассчитаем суммарные затраты для количества причалов n = 2, 3, 4. Расчеты приведены в таблице. Как видно из таблицы, минимальные затраты достигаются при n = 3. Следовательно, для минимизации затрат необходимо 3 причала.

Таблица 1.- Расчет оптимального числа причалов

Показатель	Количество причалов

Интенсивность потока судов	0,75	0,75	0,75
Интенсивность обслуживания судов	0,5	0,5	0,5
Интенсивность нагрузки причала	1,5	1,5	1,5
Вероятность, что все причалы свободны	0,14	0,21	0,22
Среднее число судов в очереди	1,93	0,24	0,04
Среднее время пребывания судна в очереди, сут.	2,57	0,32	0,06
Среднее число судов в порту	3,43	1,74	1,54
Среднее время пребывания судна в порту, сут	4,57	2,32	2,06
Пеня за простой судна в порту, ден.ед./сут. ()	100,00	100,00	100,00
Затраты по обслуживанию причала в сутки, ден.ед./сут. ()	150,00	150,00	150,00
Суммарная пеня за простой судов в порту в сутки, ден.ед. ()	192,86	23,68	4,48
Суммарные затраты по обслуживанию причалов в сутки, ден.ед. ()	300,00	450,00	600,00
Суммарные затраты, ден.ед.(С)	492,86	473,68	604,48

Варианты заданий

Таблица 2 - Варианты заданий

Номер варианта

Задача
Номер варианта

Задача

1. В парикмахерской в зависимости от сложности стрижки, мастер выполняет работу в среднем за 30 мин. Посетители приходят в среднем через 25 мин. За каждый час работы мастер зарабатывает 300 ден.ед.. Очередь ограничена до 4 человек. Если в очереди больше 4 человек, клиент уходит, и потери за час составляют 150 ден.ед. Определить предельные вероятности состояний и характеристики обслуживания. Определить оптимальное количество мастеров.

2. Автомобили подъезжают на АЗС со средней частотой 2 автомобиля за 5 минут. Заправка автомобиля в среднем длится 3 минуты. Определить предельные вероятности состояний и характеристики обслуживания. Определить количество колонок, чтобы средняя длина очереди не превышала 3 авт.

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

4. В срочной мастерской по починке обуви в зависимости от сложности ремонта мастеру требуется в среднем 15 мин. Посетители приходят в среднем через каждые 14 мин. Определить предельные вероятности состояний и характеристики обслуживания. Определить количество мастеров, чтобы средняя длина очереди не превышала 5 заказов.

5. В справочной оператор дает справку в среднем за 4 мин. Звонки поступают каждые 3мин. Если операторы заняты, то звонок не обслуживается. Определить вероятности состояний и характеристики обслуживания справочной. Определить количество каналов, чтобы относительная пропускная способность была не меньше 0,75.

6. В зависимости от количества продуктов у покупателя кассиру в магазине требуется в среднем на один чек 2 мин. Покупатели подходят к кассе с интенсивностью 81 человек/час. Определить предельные вероятности состояний и характеристики обслуживания. Определить количество кассиров, чтобы средняя длина очереди не превышала 4 покупателей.

7. Диспетчеру в АТП в зависимости от типа автомобиля требуется в среднем на выдачу одного маршрутного листа 20 минут. Заявки на автомобили поступают в среднем через каждые 30 минут. Определить предельные вероятности состояний и характеристики обслуживания. Определить количество диспетчеров, чтобы средняя длина очереди не превышала 2 заявок.

8. Требуется оценить работу АТС. Если все линий связи заняты, то абонент выбывает из системы. Звонки поступают с интенсивностью 2 вызов/мин.. Продолжительность разговоров распределена экспоненциально, и в среднем равна 1,5 мин. Определить предельные вероятности и показатели эффективности системы. Определить количество операторов, чтобы относительная пропускная способность АТС была не меньше 0,9.

9. В банке в зависимости от сложности запроса клиента кассиру требуется в среднем 10 минут. Клиенты подходят к нему в среднем через каждые 12 минут. Кассир зарабатывает 15000 ден.ед. за месяц. Очередь ограничена до 6 человек. Если в очереди больше 6 человек, клиент уходит, и потери за час составляют 200 ден.ед. Определить предельные вероятности состояний и характеристики обслуживания. Определить оптимальное количество кассиров.

10. В среднем на одну транзакцию у банкомата уходит 2 минуты. Клиенты подходят к нему в среднем через каждые 20 минут. Определить предельные вероятности состояний и характеристики обслуживания. Определить количество банкоматов, чтобы средняя длина очереди не превышала 2 человек.

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

12. В отделе заказов мебельной фабрики менеджеру по продажам в зависимости от заказа клиента требуется в среднем на оформление одного заказа 25 минут. Клиенты приходят в среднем через каждые 30 минут. Определить предельные вероятности состояний и характеристики обслуживания. Определить количество менеджеров, чтобы средняя длина очереди не превышала 3 человек.

Порядок выполнения работы

1.Рассчитайте в системе Excel показатели системы массового обслуживания по формулам, приведенным в методичке. Количество каналов обслуживания n=1, 2, 3...k перебирается для нахождения оптимального значения по варианту. Предполагается, что входные потоки и обслуживание соответствуют пуассоновскому распределению.

2.Проведите анализ полученных результатов.

3.Составьте отчет.

1) Цель работы;

2) постановка задачи;

3) результаты расчетов, проведенных в Excel;

4) выводы по выполнению работы.

Контрольные вопросы

1. Что включает в себя понятие система массового обслуживания?

2. Какие существуют виды систем массового обслуживания?

3. Что относится к основным характеристикам и показателям эффективности систем массового обслуживания?

4. Укажите основные свойства (характеристики) входящего потока требований?

5. Перечислите основные особенности и характеристики систем массового обслуживания с ожиданием?

6. Каковы основные характеристики СМО с отказами?

7. Приведите примеры различных видов СМО?

Библиографический список

1. Афанасьев М.Ю. Исследование операций в экономике: модели, задачи, решения. / М.Ю. Афанасьев, Б.П. Суворов.- М.:ИНФРА, 2003.-444с.

2. Вентцель Е.С. Исследование операций. Задачи, приниципы, методология./ Е.С. Вентцель.-М.: Высшая школа, 2001.-208с.

3. Зайченко Ю.П. Исследование операций./ Ю.П. Зайченко.- К.: Вища школа, 1975.-320с.

4. Конюховский П.В. Математические методы исследования операций. / П.В. Конюховский.- СПб.: Питер, 2001.-192с.

5. Кремер Н.Ш., Путко Б.А. Исследование операций в экономике./ Н.Ш. Кремер, Б.А. Бутко, И.М. Тришин.- М.:Банки и биржи, ЮНИТИ, 1997.-407с.

1. Кудрявцев Е.М. GPSS World.Основы имитационного моделирования различных систем.- М.: ДМК Пресс, 2004.- 320 с.

2. Советов В.Я., Яковлев С.А. Моделирование систем. - М.: Высшая школа, 1985

3. Советов В.Я., Яковлев С.А. Моделирование систем: курсовое проектирование. - М.: Высшая школа, 1989

Рассмотрим одноканальную систему массового обслуживания с ожиданием.

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

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

Заявка, поступившая в момент, когда канал занят, становится в очередь и ожидает обслуживания. Будем считать, что размер очереди ограничен и не может вместить более m заявок, т.е. заявка, заставшая в момент своего прихода в СМО m +1 заявок (m ожидающих в очереди и одну, находящуюся на обслуживании), покидает СМО.

Система уравнений, описывающих процесс в этой системе, имеет решение:

(0‑1)

Знаменатель первого выражения представляет собой геометрическую прогрессию с первым членом 1 и знаменателем ρ, откуда получаем

При ρ = 1 можно прибегнуть к прямому подсчету

(0‑8)

Среднее число находящихся в системе заявок.

Поскольку среднее число находящихся в системе заявок

(0‑9)

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

(0‑10)

и среднее число находящихся в системе заявок равно

(0‑11)

Среднее время ожидания заявки в очереди .

(0‑12)

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

Среднее время пребывания заявки в системе.

Время пребывания заявки в системе складывается из времени ожидания заявки в очереди и времени обслуживания. Если загрузка системы составляет 100%, то =1/μ, в противном случае = q / μ . Отсюда

(0‑13)

Содержание работы .

Подготовка инструментария эксперимента .

Выполняется аналогично в соответствии с общими правилами.

Расчет на аналитической модели .

1. В приложение Microsoft Excel подготовьте таблицу следующего вида.

2. В столбцах для параметров СМО таблицы запишите исходные данные, которые определяются по правилу:

m=1,2,3

(максимальная длина очереди).

Для каждого значения m необходимо найти теоретические и экспериментальные значения показателей СМО для таких пар значений:

= <порядковый номер в списке группы>

3. В столбцы с показателями аналитической модели впишите соответствующие формулы.

Эксперимент на имитационной модели .

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

2. Для каждой комбинации m , и осуществите запуск модели.

3. Результаты запусков внесите в таблицу.

4. Внесите в соответствующие столбцы таблицы формулы для расчета среднего значения показателя P отк , q и А.

Анализ результатов .

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

2. Для m=3 постройте на одной диаграмме графики зависимости P отк от на теоретически и экспериментально полученных данных.

Оптимизация параметров СМО .

Решите задачу оптимизации размера числа мест в очереди m для прибора со средним временем обслуживания = с точки зрения получения максимальной прибыли. В качестве условий задачи возьмите:

- доход от обслуживания одной заявки равным 80 у.е./час,

- стоимость содержания одного прибора равным 1у.е./час.

1. Для расчетов целесообразно создать таблицу:

Первый столбец заполняется значениями чисел натурального ряда (1,2,3…).

Все клетки второго и третьего столбцов заполняются значениями и.

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

В столбцы с исходными данными разделов Доход, Расход, Прибыль внесите значения (см. выше).

В столбцах с вычисляемыми значениями разделов Доход, Расход, Прибыль запишите расчетные формулы:

- число заявок в единицу времени

N r =A

- суммарный доход в единицу времени

I S = I r *N r

- суммарный расход в единицу времени

E S =E s + E q *(n-1)

- прибыль в единицу времени

P = I S - E S

где

I r - доход от одной заявки ,

E s - расход на эксплуатацию одного прибора ,

E q - расход на эксплуатацию одного места в очереди .

Графики для P отк ,

- таблицу с данными для нахождения наилучшего m и значение m опт,

- график зависимости прибыли в единицу времени от m .

Контрольные вопросы :

1) Дайте краткое описание одноканальной модели СМО с ограниченной очередью.

2) Какими показателями характеризуется функционирование одноканальной СМО с отказами?

3) Как рассчитывается вероятность p 0 ?

4) Как рассчитываются вероятности p i ?

5) Как найти вероятность отказа обслуживания заявки?

6) Как найти относительную пропускную способность?

7) Чему равна абсолютная пропускная способность?

8) Как подсчитывается среднее число заявок в системе?

9) Приведите примеры СМО с ограниченной очередью.

Задачи .

1) Порт имеет один грузовой причал для разгрузки судов. Интенсивность потока составляет 0,5 заходов в сутки. Среднее время разгрузки одного судна 2 суток. Если в очереди на разгрузку стоят 3 судна, то приходящее судно направляется для разгрузки на другой причал. Найти показатели эффективности работы причала.

2) В справочную железнодорожного вокзала поступают телефонные запросы с интенсивностью 80 заявок в час. Оператор справочной отвечает на поступивший звонок в среднем 0,7 мин. Если оператор занят, клиенту выдается сообщение "Ждите ответа", запрос становится в очередь, длина которой не превышает 4 запросов. Дайте оценку работы справочной и вариант ее реорганизации

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

Для СМО с отказами :

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

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

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

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

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

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

Рассмотрим примеры некоторых СМО.

2.5.1. Многоканальная СМО с отказами

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

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

Описание состояний:

Все инспекторы свободны;

Занят один инспектор;

Заняты два инспектора;

Заняты три инспектора.

Граф состояний системы приведен на рис. 2.11 .

Рис. 2.11.

На графе: - интенсивность потока грузовых автомобилей; - интенсивность проверок документов одним автоинспектором.

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

Решение

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

Пропускную способность этого поста автоинспекторов можно характеризовать относительной пропускной способностью :

Пример 2.6 . Для приема и обработки донесений от разведгруппы в разведотделе объединения назначена группа в составе трех офицеров. Ожидаемая интенсивность потока донесений - 15 донесений в час. Среднее время обработки одного донесения одним офицером - . Каждый офицер может принимать донесения от любой разведгруппы. Освободившийся офицер обрабатывает последнее из поступивших донесений. Поступающие донесения должны обрабатываться с вероятностью не менее 95 %.

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

Решение

Группа офицеров работает как СМО с отказами, состоящая из трех каналов.

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

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

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

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

В примере .

В рассматриваемой СМО отказ наступает при занятости всех трех каналов, то есть . Тогда:

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

Таким образом, только группа из шести офицеров сможет обрабатывать поступающие донесения с вероятностью 95 %.

2.5.2. Многоканальная СМО с ожиданием

Пример 2.7 . На участке форсирования реки имеются 15 однотипных переправочных средств. Поток поступления техники на переправу в среднем составляет 1 ед./мин, среднее время переправы одной единицы техники - 10 мин (с учетом возвращения назад переправочного средства).

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

Решение

Абсолютная пропускная способность , т. е. все, что подходит к переправе, тут же практически переправляется.

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

Коэффициенты использования и простоя переправы:

Для решения примера была также разработана программа. Интервалы времени поступления техники на переправу, время переправы приняты распределенными по экспоненциальному закону.

Коэффициенты использования переправы после 50 прогонов практически совпадают: .

Максимальная длина очереди 15 ед., среднее время пребывания в очереди около 10 мин.

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

Все каналы свободны, ;

Занят только один канал (любой), ;

- заняты только два канала (любых), ;
- заняты все каналов, .

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

Заняты все каналов и одна заявка стоит в очереди,

Заняты все каналов и две заявки стоят в очереди,

Заняты все каналов и все мест в очереди,

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

Запишем выражения для предельных вероятностей состояний:

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

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

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

Поэтому вероятность образования очереди равна: