- •Набережные Челны 1999
- •2. Теория массового обслуживания
- •2.1 Задачи теории массового обслуживания
- •2.2. Классификация систем массового обслуживания и их основные характеристики
- •2.4. Многоканальная смо с отказами
- •2.5. Одноканальная смо с ожиданием
- •2.8. Многоканальная смо с ожиданием
- •2.9. Смо с ограниченным временем ожидания
- •2.10. Замкнутые системы массового обслуживания
- •2.11. Системы массового обслуживния с не пуассоновскими потоками событий
- •3.3. Оценка количества испытаний,
- •3.4. Генерирование, равномерно
- •3.6. Имитация случайных величин
- •3.7. Имитация случайных потоков
- •3.8. Имитация процессов функционирования
2.5. Одноканальная смо с ожиданием
(число мест в очереди ограничено)
Рассмотрим сначала простейшую из всех возможных СМО с ожиданием - одноканалъную систему (п = 1) , на которую поступает поток заявок с интенсивностью λ ; интенсивность обслуживания µ (т.е. в среднем непрерывно занятый канал будет выдавать µ обслуженных заявок в единицу времени). Заявка, поступившая в момент, когда канал занят, становится в очередь и ожидает обслуживания.
Предположим сначала, что количество мест в очереди ограничено числом т, т.е. если заявка пришла в момент, когда в очереди уже стоят т заявок, она покидает систему необслуженной. В дальнейшем, устремив т к бесконечности, получим характеристики одноканальной СМО без ограничений по длине очереди.
Будем нумеровать состояния СМО по числу заявок, находящихся в системе ( как обслуживаемых, так и ожидающих обслуживания):
S0 - канал свободен,
S1 - канал занят, очереди нет,
S2 - канал занят, одна заявка стоит в очереди,
............................................................................
Sm+1 - канал занят, т заявок стоят в очереди.
Граф состояний СМО показан на рис.23. Интенсивности потоков событий, переводящих систему по стрелкам слева направо, все равны λ, а справа налево - µ . Действительно, по стрелкам слева направо систему переводит поток заявок (как только придет заявка, система переходит в следующее состояние), справа же налево - поток " освобождений " занятого канала, имеющий интенсивность µ (как только будет обслужена очередная заявка, канал либо освободится, либо уменьшится число заявок в очереди).
Изображенная на рис.23 схема представляет собой схему гибели и размножения. Пользуясь общим решением, данным для схемы гибели и размножения, напишем выражения предельных вероятностей состояний:
(2.19)
Вводя обозначение перепишем формулы (2.19) в виде
(2.20)
Заметим, что в знаменателе последней формулы (2.20) стоит геометрическая прогрессия с первым членом 1 и знаменателем р ; суммируя эту прогрессию, находим:
(2.21)
Обратим внимание на то, что
формула (2.21) справедливо только при
(при
р=1 она
дает неопределенность вида 0/0). Но сумму
геометрической прогрессии со знаменателем
р = 1
найти еще проще, чем по формуле (2.21): она
равна m+2,
и в этом случае P0
= 1 /(m
+ 2). Заметим, что тот же результат
мы могли бы получить более сложным
способом, раскрывая
неопределенность (2.21) по правилу Лопиталя.
Определим характеристики
СМО: вероятность отказа Ротк,
относительную пропускную способность
q,
абсолютную пропускную способность А
, среднюю длину очереди
,
среднее число заявок, связанных с
системой
.
О
чевидно,
заявка получает отказ только в случае,
когда канал занят
и все т мест
в очереди - тоже:
(2.22)
Находим относительную пропускную способность:
(2.23)
Абсолютная пропускная способность:
A
= λq.
Найдем среднее число заявок, находящихся в очереди. Определим эту величину как математическое ожидание дискретной случайной величины R - числа заявок, находящихся в очереди:
С вероятностью Р2 в очереди стоит одна заявка, с вероятностью Р3 две заявки, вообще с вероятностью Рk в очереди стоят к — 1 заявок, наконец, с вероятностью Рт+1 в очереди стоят т заявок. Среднее число заявок в очереди получим, умножая число заявок в очереди на соответствующую вероятность и складывая результаты:
(2.24)
Выведем формулу для суммы, стоящей в скобках. Очевидно, рассматриваемая сумма представляет собой не что иное, как производную пo p суммы
а для этого выражения можем воспользоваться формулой суммы геометрической прогрессии:
(2.25)
Продифференцируем (2.25) по р и в результате получим
(2.26)
Подставляя его в (2.24), получим:
Учитывая выражение (2.21) для Р0 , имеем
и
ли,
окончательно
(2.27)
Таким образом, получено выражение для среднего числа заявок, ожидающих обслуживания в очереди. Выведем теперь формулу для среднего числа заявок, связанных с системой (как стоящих в очереди, так и находящихся под обслуживанием). Будем решать эту задачу следующим образом: рассмотрим общее число заявок К, связанных с системой, как сумму двух случайных величин: числа заявок, стоящих в очереди, и числа заявок, находящихся под обслуживанием:
K=R+Ω
где R - число заявок в очереди, Ω - число заявок под обслуживанием. Имеем
где
-
среднее число заявок
в очереди,
-
среднее число заявок
под обслуживанием..
Величина найдена, найдем величину . Так как канал у нас один, то случайная величина Ω может принимать только два значения: 0 или 1. Значение 0 она принимает, если канал свободен; вероятность этого равна
Значение 1 она принимает, если канал занят; вероятность этого равна
Отсюда находим математическое ожидание числа заявок, находящихся под обслуживанием:
Таким образом, среднее число заявок, связанных с СМО, будет
(2.28)
где величина r определяется по формуле (2.27).
Выведем
выражение еще для одной существенной
характеристики СМО с ожиданием;
среднего времени ожидания заявки
в очереди. Обозначим его
.
Пусть заявка приходит в систему
в какой-то момент времени. С вероятностью
P0
канал обслуживания
не будет занят и заявке не придется
стоять в очереди
(время ожидания равно нулю). С вероятностью
Р1
она
придет
в систему во время обслуживания какой-то
заявки, но перед
ней не будет очереди, а заявка будет
ждать начала своего обслуживания
в течение времени
(среднее время
обслуживания
одной заявки). С вероятностью Р2
в
очереди перед рассматриваемой
заявкой будет стоять еще одна, и время
ожидания в среднем
будет равно
,
и т.д.. Вообще, с
вероятностью Pk
пришедшая заявка
застанет в системе k
заявок и будет
ждать в среднем
единиц времени; здесь к
может
быть любым целым числом до
m
. Что же касается к = т
+ 1, т.е.
случая, когда вновь приходящая заявка
застает канал обслуживания занятым
и еще т
заявок
в очереди (вероятность этого
Рт+1),
то время ожидания в этом случае также
равно нулю,
потому что заявка не становится в очередь
(и не обслуживается).
Поэтому среднее время ожидания будет
П
одставляя
сюда выражения для Р1,...,
Рт,
получаем
Преобразуем сумму в скобках, пользуясь формулой (2.26):
или выражая Р0 через р :
(2.29)
Сравнивая это выражение с формулой (2.27), замечаем, что
(2.30)
т.е. среднее время ожидания равно среднему числу заявок в очереди, деленному на интенсивность потока заявок.
Выведем еще формулу для среднего времени пребывания заявки в системе. Обозначим Тсист случайную величину - время пребывания заявки в СМО. Эта случайная величина складывается из двух слагаемых (тоже случайных):
Tсист=Tож+θ
где Тож - время ожидания заявки в очереди, θ – случайная величина, равная времени обслуживания Тоб, если заявка обслуживается, и нулю, если она не обслуживается (получает отказ).
Имеем:
но в наших обозначениях,
a
О
тсюда
находим:
или, с учетом формулы (2.30), (2.31)
Пример. Автозаправочная станция (АЗС) представляет собой СМО с одним каналом обслуживания (одной колонкой). Площадка при станции допускает пребывание в очереди на заправку не более трех машин одновременно (m =3). Если в очереди уже находятся три машины, очередная машина, прибывшая к станции, в очередь не становится, а проезжает мимо. Поток машин, прибывающих для заправки, имеет интенсивность λ = 1 (машина в минуту). Процесс заправки продолжается в среднем 1,25 мин. Определить:
вероятность отказа;
относительную и абсолютную пропускные способности СМО;
среднее число машин, ожидающих заправки;
среднее число машин, находящихся на АЗС (включая и обслуживаемую);
среднее время ожидания машины в очереди;
среднее время пребывания машины на АЗС (включая обслуживание).
Решение. Находим приведенную интенсивность потока заявок:
По формулам (2.21), (2.20) вычисляем предельные вероятности состояний:
Относительная пропускная способность СМО
Абсолютная пропускная способность СМО
A = λq =0,703.
Среднее число машин в очереди находим по формуле (2.27):
т.е. среднее число машин, ожидающих в очереди на заправку, равно 1,56.
П
рибавляя
к этой величине среднее число машин,
находящихся под
обслуживанием
получим среднее число машин, связанных с АЗС:
С
реднее
время ожидания машины в очереди, по
формуле (2.30),
равно:
Прибавляя к этой величине
получим среднее время, которое машина проводит на АЗС:
2.6. ОДНОКАНАЛЬНАЯ СМО С ОЖИДАНИЕМ (число мест в очереди не ограничено)
В предыдущем 2.5 рассматривали работу одноканальной СМО с ожиданием при ограниченном числе т мест в очереди.
Теперь снимем это ограничение, т.е. устремим т к бесконечности. При этом число возможных состояний системы станет бесконечным, и граф состояний примет вид, показанный на рис.24.
λ λ
λ λ λ
µ
µ µ µ
µ
Попытаемся
получить вероятности состояний СМО с
неограниченной очередью путем предельного
перехода (при т
)
из формул (2.20), (2.21).
Заметим, что при этом знаменатель в последней формуле (2.20) представляет собой сумму бесконечного числа членов геометрической прогрессии. Эта сумма сходится только, когда прогрессия бесконечно убывающая, т.е. р < 1. Можно совершенно строго доказать, что р < 1 есть условие, при котором в СМО с ожиданием существует предельный установившийся режим; при р>1 такого режима не существует, и очередь при t растет до бесконечности. Предположим, что
т.е. предельный режим
существует. Устремим в формулах (2.21),
(2.20) т к
и
выведем формулы для предельных
вероятностей
состояний в СМО без ограничений по длине
очереди.
Получим:
(2.32)
При отсутствии ограничений по длине очереди каждая заявка, пришедшая в систему, будет обслужена, поэтому
q = 1, А = λq = λ.
Среднее число заявок в очереди получим из (2.27) при т :
(2.33)
Среднее число заявок в системе по формуле (2.28) при т будет равно:
(2.34)
Среднее время ожидания
также получим из
формулы (2.30) при т
(2.35)
Среднее время пребывания заявки в СМО равно среднему времени ожидания плюс среднее время обслуживания
(2.36)
Пример. На железнодорожную сортировочную горку прибывают составы с интенсивностью λ = 2 (состава в час). Среднее время, в течение которого горка обрабатывает состав, равно 0,4 часа, причем одновременно обрабатывается лишь один состав. Составы, прибывшие в момент, когда горка занята, становятся в очередь и ожидают в парке прибытия, где имеются три запасных пути, на каждом из которых может ожидать один состав. Состав, прибывший в момент, когда все три запасных пути в парке прибытия заняты, становится в очередь на внешний путь, где число мест для ожидания не ограничено. Все потоки событий - простейшие. Найти:
среднее число составов, ожидающих в очереди (как в парке прибытия, так и вне его);
среднее число составов, ожидающих на внешних путях;
среднее время ожидания состава в парке прибытия и на внешних путях;
среднее время нахождения состава на сортировочной станции (включая ожидание и обслуживание);
вероятность того, что прибывший состав займет место на внешних путях.
Решение.
В нашем случае λ
= 2, µ
=
=
2,5
,
и СМО в среднем "справляется"
с поступившим на нее потоком
заявок.
Среднее число составов, ожидающих очереди (как в парке прибытия, так и вне его), найдем по формуле (2.33):
Среднее число составов, ожидающих очереди на внешних путях, подсчитываем так: с вероятностью Р5 вне парка прибытия будет ожидать один состав, с вероятностью
Р6 - два состава, и т.д. с вероятностью Рk {к > 5) – (к — 4) состава.
С
реднее
число составов, ожидающих вне парка,
будет:
Формулу для бесконечной суммы в скобках получаем
предельным переходом.(при т ) из формулы (2.26):
Подставляя сюда р
= 0,8 , получим:
=
1,64.
Вероятность того, что прибывший состав займет место на внешних путях, определяется еще проще: она равна вероятности того, что длина очереди будет не меньше трех, т.е.
Среднее время ожидания в парке прибытия определяем, рассматривая различные гипотезы о числе составов, находящихся в системе:
Для ρ= 0,8 , µ = 2,5 , получаем, что среднее время ожидания в парке прибытия равно
Что же касается среднего времени ожидания на внешних путях, то оно равно:
т.е. для наших численных данных
Среднее время пребывания состава на сортировочной станции (считая ожидание и обслуживание) будет равно:
2.7. МНОГОКАНАЛЬНАЯ СМО С ОЖИДАНИЕМ (число мест в очереди ограничено)
Рассмотрим п - канальную СМО с ожиданием, на которую поступает поток заявок с интенсивностью λ ; интенсивность обслуживания (для одного канала) µ число мест очереди т .
Состояние системы будем нумеровать по числу заявок, связанных с системой:
S0 - все каналы свободны,
S1 - занят один канал, остальные свободны,
.............................................................................
Sk -заняты к каналов, остальные свободны,
..........................................................................
Sn - заняты все п каналов, очереди нет,
Sn+1 - заняты все п каналов, одна заявка стоит в очереди,
.................................................................................
Sn+r - заняты все п каналов, r заявок стоит в очереди,
..........................................................................................
Sn+m - заняты все п каналов, т заявок стоит в очереди.
Граф состояний приведен на рис.25. У каждой стрелки поставлены соответствующие интенсивности потоков событий. Действительно, по стрелкам слева направо систему переводит всегда один и тот же поток заявок с интенсивностью λ; по стрелкам справа налево систему переводит поток обслуживании, интенсивность которого равна µ умноженному на число занятых каналов.
Граф на рис.25 представляет собой схему гибели и размножения, для которой решение в общем виде уже получено. Напишем выражения для предельных вероятностей состояний, а сразу же обозначая
(2.37)
или,
суммируя геометрическую прогрессию с
знаменателем p/n
и первым членом
в
последней сумме
(2.38)
Таким образом, все вероятности состояний найдены.
Найдем некоторые характеристики эффективности обслуживания. Поступившая заявка получает отказ, если заняты все п каналов и все т мест в очереди:
(2.40)
Относительная пропускная способность, как всегда, дополняет вероятность отказа до единицы:
Абсолютная пропускная способность СMO будет равна
(2.40)
Найдем
среднее
число занятых каналов,
которое
обозначим
.
Каждый занятый канал
обслуживает в среднем µ
заявок
в единицу времени; вся же СМО обслуживает
в среднем А
заявок в единицу
времени. Деля одно на другое, получим:
(2.41)
Среднее число заявок в очереди можно вычислить непосредственно, как математическое ожидание дискретной случайной величины, умножая любое возможное число заявок на вероятность того, что именно это число заявок будет в очереди, и складывая результаты:
(2.42)
Введем обозначение р/п =χ и перепишем (2.42) в виде :
(2.43)
Заметим, что выражение в скобках есть не что иное, как уже вычисленная сумма (2.26 ), где вместо р поставлено χ .Пользуясь этой формулой и подставляя результат в (2.43), получим
(2.44)
Складывая среднее число заявок в очереди r и среднее число занятых каналов z , получим среднее число заявок, связанных с системой к:
к =r+z (2.45)
Найдем среднее время ожидания заявки в очереди .
Сделаем ряд гипотез о том, в каком состоянии застанет систему вновь пришедшая заявка и сколько времени ей придется ждать обслуживания.
Если заявка застанет не все каналы занятыми, ей вообще не придется ждать (соответствующие члены в математическом ожидании отбросим, как равные нулю). Если заявка придет в момент, когда заняты все п каналов, а очереди нет, ей придется ждать в среднем время, равное 1/nµ (потому что поток освобождении п каналов имеет интенсивность nµ ). Если заявка застанет все каналы занятыми и одну заявку перед собой в очереди, ей придется в среднем ждать время 2/nµ ( по 1/nµ на каждую впереди стоящую заявку) и т.д. Если заявка застанет в очереди r заявок, ей придется ждать в среднем время r/пµ. Если вновь пришедшая заявка застанет в очереди уже т заявок, но она вообще не будет ждать (но и не будет обслуживаться). Среднее время ожидания найдем, умножая каждое из этих значений на соответствующую вероятность:
Так же, как и в случае одноканальной СМО с ожиданием, замечаем, что это выражение отличается от выражения для средней длины очереди (2.42) только множителем 1/pµ=1/λ, т.е.
(2.46)
Среднее время пребывания заявки в системе, так же, как и для одноканальной СМО, отличается от среднего времени ожидания на среднее время обслуживания, умноженное на относительную пропускную способность:
Пример. Автозаправочная станция (АЗС) с двумя колонками (п = 2) предназначена для обслуживания машин. Поток автомашин, прибывающих на АЗС, имеет интенсивность λ = 2 (машин в минуту); среднее время обслуживания одной машины
Площадка у АЗС может вместить очередь не более m=3 (машин). Машина, прибывшая в момент, когда все три места в очереди заняты, покидает АЗС (получает отказ). Найти характеристики СМО:
вероятность отказа,
относительную и абсолютную пропускную способность,
среднее число занятых колонок,
среднее число машин в очереди,
среднее время ожидания и пребывания машины на АЗС,
