книги / Прикладные задачи теории массового обслуживания
..pdfВероятность противоположного события T04^ t будет равна
Р ( Т 0ч < 0 = 1 - Р (Гоч > *). |
(5.1.27) |
Типичный график этой зависимости показан на рис. 5.1.6.
Функция распределения случайной величины Точ определяет
ся так:
ч V) = Р (Топ < ty= 1 - Р (Гоч > t ) - P (Гоч = *)• (5.1.28)
Вероятность P(T04 = t) имеет смысл только при ^= 0. График
функции распределения показан на рис. 5.1.7.
Мы видим, что случайная величина Точ является случайной
величиной смешанного типа, так как функция распределения имеет скачок и интервал непрерывного возрастания. Для вычис ления различных начальных моментов случайной величины Точ
можно воспользоваться следующей формулой:
«Л7’о . ,] = - ^ * -^ -P (T 04> t)dt. |
(5.1.29) |
О |
|
С помощью этой формулы при k= l можно получить выраже
ние для математического ожидания времени нахождения в оче реди:
^1^оч] = ^оч-
Для нахождения функции----—Р(Т0Ч> t) можно использо-
Ot
вать следующий прием. Время ожидания Ттданной заявки в оче
реди при условии, что перед ней стоит в очереди г заявок (см. рис. 5.1.5), равно сумме (г+1) времени обслуживания в предпо ложении, что работают все п каналов. Так как работают все ка
налы, то время, протекающее между вызовами двух очередных заявок из очереди, распределено по показательному закону с параметром n\i. Следовательно:
Г г = 2 Tt (г = 0, 1, . . . , т - 1 ) , |
(5.1.30) |
i= i |
|
181
где Ti(i=l, |
r + 1 ) — независимые случайные |
величины, рас |
|||||
|
|
|
пределенные одинаково по показательно |
||||
|
|
|
му закону с параметром лц. |
|
|
||
С подобной суммой мы сталкивались раньше в § |
1.5 при изу |
||||||
чении |
потока |
Эрланга. |
Там было показано, |
что |
сумма |
||
(5.1.30) |
подчинена |
закону |
Эрланга r-го порядка, т. |
е. [см. |
|||
(1.5.22)] |
плотность |
распределения величины Т,- |
будет иметь вид |
|
f r (t)= П|^ |
° Г е~^‘. |
(5.1.31) |
||
Теперь по формуле полной вероятности можно |
найти функ- |
||||
ци ю ---- ^ - Р ( Г оч> 0 : |
|
|
|
|
|
|
m — 1 |
|
m — 1 |
|
|
~ - ~ P ( T c, > t ) = X М 0 р» г= Х Л ^ И е- ^ ’р ,= |
|||||
|
г= 0 |
|
г = 0 |
|
|
= |
/7,)Л!хе-(л1-Х)</?(/л — 1, |
/i). |
(5.1.32) |
||
Заметим, что |
|
|
|
|
|
со |
|
ос т — 1 |
|
|
|
0 |
|
0 |
г—0 |
|
|
т — 1 |
|
оо |
|
т — 1 |
|
= Рп j |
(/г; Г |
1 \ tre~^dt= , j v//7„= |
|||
r = О |
’ |
0 |
|
r = О |
|
|
//2 — 1 |
|
I — У-т |
, , |
|
|
|
|
|
||
|
|
|
рп+г= рп 1 — - < \ . |
гS=0
Этот интеграл не равен единице, так как он не учитывает то го обстоятельства, что заявка вообще может не находиться в оче реди: получит немедленный отказ или немедленно будет приня та к обслуживанию.
Единице будет равно следующее выражение:
~ \ ~ t p {Т°" > tydt+ р {Т°ч= 0 ) = 1 ’
0
где
Р(Точ = 0) = \ - р п^ —
есть вероятность того, что время пребывания в очереди будет равно нулю.
182
Если ввести в рассмотрение дельта-функцию 6(х) (см. § 1.5), то можно записать выражение для плотности распределения вре мени ожидания в очереди /0ч(0- Напомним, что дельта-функция обладает следующими основными свойствами:
o(jt)=0 при хфО]
^<o(x)b(x)dx=<?(0),
если ф(л) непрерывна в точке дс=0; 8 (х) ф (л) = 0;
если ф (х) — нечетная функция.
Плотность распределения времени пребывания заявки в оче реди имеет вид
pnny.e~^-l)‘R (т — 1, >i) 4- о (г) X
/оч(*) = |
|
|
(5.1.33) |
|
|
|
|
||
P,;KR (т - 1 , |
х/) + |
8 (t) (1 - рпт) (•/.= 1), |
||
г д е |
Р (л, |
а) |
|
|
I |
1 — %т (у- т-1); |
|||
Р (л, |
|
|||
i R (п, I) + |
а) * —-------- |
|||
Рп=\ |
|
|
(5.1.34) |
|
|
Р(п, |
п) |
_____ |
|
R (я, |
п) + Р (л, |
л) т |
Если требуется определить лишь среднее время нахождения заявки в очереди ?0ч, то закон распределения времени 70ч отыс кивать не требуется. Покажем это (см. гл. 3). Допустим, что время ожидания данной заявки в очереди попало в элементарный ин тервал (t, t+dt). Вероятность этой гипотезы приближенно рав на f04{t)dt. За время пребывания заявки в очереди за этой за
явкой образуется очередь, в которой в среднем будет находить ся Kt заявок. Следовательно, полное математическое ожидание
числа заявок, находящихся в очереди, будет определяться из выражения
7 = ^ } J f04(t)dt—lJ04y |
(5.1.35) |
О |
|
откуда формула для определения среднего времени ожидания за явки в очереди примет вид
г
(5.1.36)
183
Среднее время нахождения заявки в системе
Т __ г Ч~ k |
__ |
I _ |
(5.1.37) |
||
~ |
л |
_ |
X ’ |
||
|
|||||
где 1 = г + к — среднее число |
заявок, находящихся в системе. |
||||
Таким образом, для отыскания |
среднего времени ожидания |
заявки в очереди требуется лишь знание среднего числа заявок, находящихся в очереди. Таковы основные характеристики време ни пребывания заявок в очереди.
Рассмотрим теперь случай, когда вероятность успешного обслуживания заявки, побывавшей на обслуживании, равна р< 1. Все характеристики системы, кроме вероятности Р0бС, оста
ются неизменными. Вероятность обслуживания заявки в этом случае определяется по формуле
Я о б с = р(1 — Р п + т )- |
(5.1.38) |
Абсолютная пропускная способность будет равна |
|
Х0 = }.Po6c= \р (1 — рп+ т). |
(5.1.39) |
Обратим внимание читателя на то, что при рассмотрении ста ционарного режима работы СМО с ожиданием мы ограничились лишь случаем ограниченного числа мест в очереди (/п<оо). При неограниченном числе мест в очереди (т = оо) стационарный ре
жим работы системы существует только при условии, что
-/. = — < 1 . |
(5.1.40) |
ГЦХ
Это условие указывает на то, что среднее число заявок, обслу живаемых в единицу времени всеми п каналами, должно быть
больше среднего числа заявок, поступающих в систему в единицу времени. Если это условие не выполняется и число мест в оче реди не ограничено (т = оо), то стационарного режима работы
наблюдаться не будет. Это можно объяснить следующим обра
зом. При |
х = — > |
1 система |
все время |
(неограниченно) будет |
|
|
П\х |
|
|
|
|
двигаться |
«вправо», |
т. е. перемещаться |
в сторону |
состояний с |
|
большим |
числом г, |
и очередь |
будет неограниченно |
возрастать. |
Однако стабилизации этого движения «вправо» не наступит, так как система не справляется с потоком заявок, что приводит к
постоянному |
увеличению |
очереди. Заметим, |
что если и ^ 1 , то |
для любого |
к о н е ч н о г о |
г |
|
|
Нш рп+г (/) = 0, |
(5.1.41) |
|
|
t-+oo |
|
так как система рано или поздно пройдет это состояние и практи чески в него уже не вернется.
184
Поэтому |
|
lim *P.”+'W = 0. |
(5.1.42) |
О д н а к о в э т о м с л у ч а е у р а в н е н и я |
( 5 . 1 . 2 ) и ( 5 . 1 . 3 ) у ж |
н е с п р а в е д л и в ы . |
|
Рассмотрим встречающийся на практике частный случай, когда имеется всего один обслуживающий канал и т мест в очереди. В этом случае величина к [см. (5.1.4) и (5.1.5)] будет равна величине а и вероятность того, что в системе будет k зая
вок (одна заявка обслуживается, а остальные ожидают в очере ди), будет равна
^ = - Г ^ а! |
+ |
(5.1.43) |
1 — а"* т 1 |
|
|
Э т о |
в ы р а |
ж е н и е |
п о л у ч а е т с я , |
е с л и в |
в |
ы (5р .а1ж.2)е нии (5я х.1.3) |
||
п о л о ж и т ь |
п = 1 |
и х = а . |
|
|
|
|
||
Вероятность обслуживания определяется из выражения |
||||||||
|
|
|
Л>бс= 1 —/ W |
i = |
|
|
(5.1.44) |
|
Среднее число занятых каналов (или вероятность того, что |
||||||||
единственный канал будет занят) |
равно |
|
|
|||||
|
|
k |
З.к |
обе |
|
l — gM-i-l |
|
(5.1.45) |
|
|
|
1 — <хт +2 |
|
||||
|
|
V- |
|
|
||||
|
|
|
|
|
|
|
||
Среднее число заявок, находящихся в очереди, будет |
||||||||
|
т |
|
т + 1 |
|
т -И |
|
rn + 1 |
|
|
|
|
*= 1 |
|
|
f - > |
^ |
Р к = |
|
|
|
|
*= 1 |
|
Л=1 |
||
= |
l - > |
~2 {[("* + 1)а” +а4-"-(-1~ ^ " ■->] - » ( ! - » « + « ) } (5.1.46) |
||||||
С р е д н е е |
в р е м я |
о ж и д а н и я |
в о ч е р е д и |
|
р а в н о |
|||
|
|
|
|
t04 |
Г |
|
|
(5.1.47) |
|
|
|
|
7 |
|
При неограниченном числе мест в очереди (ш = со) стацио нарный режим существует только при условии а< 1 . В этом случае получим следующее выражение для вероятности наличия k заявок в системе:
/7й = а * ( 1 — а ) |
( * = 0 , 1 , 2 , . . . ) . |
Такое распределение вероятностей называется распределени ем Паскаля. Оно уже встречалось у нас в гл. 1.
1 8 5
Вероятность обслуживания в этом случае будет равна едини це, вероятность того, что канал занят, будет [см. (5.1.45)]:
^зк = А = а. |
(5.1.48) |
Среднее число мест в очереди найдем по формуле
оо |
оо |
|
T = ^ r p 1+r = V ( k - l ) Pk= ^ — . |
(5.1.49) |
r=0 k=l
Среднее время ожидания в очереди равно
А)ч= - 7 - = - — т~ — |
(5-1.50) |
К |JL 1 — а |
|
Среднее время нахождения заявки в системенайдем по формуле
T r= L ± ± - = — L _ |
(5.1.51) |
ц— X
Закон распределения времени ожидания в очереди определим из выражения (5.1.27) при условии, что п=1; к = а и
HmA?(/tt — 1, Х/)=1; |
|
|
|
||
linrxm= |
Нт ат = 0 , |
|
|
|
|
TTL —*■оо |
Л 1-+ оо |
|
|
|
|
так как а < 1 . |
|
|
|
|
|
jjg—ц.(1—а)/ |
при |
/ > 0 ; |
|
||
Р(Т0ЧХ ) = |
1 |
при |
|
/ < 0 , |
(5.1.52) |
|
|
|
|||
|
откуда |
функция |
распределе |
||
|
ния F04(t) |
будет иметь вид |
|||
|
^оч(0 = |
1 - а е - '1<1- а)<( /> 0 ) . |
|||
|
|
|
|
|
(5.1.53) |
|
График |
функции |
распределе |
||
|
ния времени ожидания в очере |
||||
|
ди |
имеет |
вид, показанный на |
||
|
рис. 5.1.8. |
|
|
О с н о в н ы е р а с ч е т н ы е ф о р м у л ы
Вероятность того, что занято ровно k каналов, а очереди нет:
Р (k, а) |
■(k= 0, 1, 2........ п), |
Р ь = |
|
Р (л. а) + Р (л, а) V. |
1— х т |
1—х |
186
где
fi ** ли
Вероятность того, что все каналы заняты и в очереди имеется г заявок:
р „ + г = у ' р п ( r = 0 , 1, 2 , . . . , от).
При х = 1 (а = л) |
получим |
|
_________Я (£, л) |
(* = 0, 1 , . . га); |
|
Л" |
R(n, п) + тР (л, п) |
|
|
|
Рп+г=Р„ (Г—0, 1, 2,..., гаг).
Вероятность обслуживания равна
^обс = 1 —%тРп-
Среднее число занятых каналов (среднее число обслуживае мых заявок)
/г = а(1 — АтРп).
Среднее число заявок, находящихся в очереди:
Г — |
п |
(1 — ^.)2 |
4 Г П |
|
т (т+ |
1) |
|
|
|
|
, |
1 v |
||
|
Л.-----5----- |
(*= 1 ). |
||
Среднее время нахождения |
заявки в очереди |
Среднее время нахождения заявок в системе
t = ~ + k .
|
З а д а ч и и у п р а ж н е н и я |
|
|
||
5.1.1. |
Показать, |
что для любого |
т > 0 и любых |
параметров |
|
( п , Л , JLI) |
СМО с ожиданием имеет большую пропускную способ |
||||
ность, чем СМО с отказами и с теми же параметрами |
(я, |
Л, р). |
|||
Р е ш е н и е |
|
|
|
|
|
Вероятность обслуживания заявки для |
СМО с отказами |
равна |
|||
|
p i 1) _ R ( t l — 1, а) _ - __ Р ( п , а) |
|
|
||
|
°6С |
R(n , a) |
/?(л, 2 ) |
|
|
187
Вероятность обслуживания заявки для СМО с ожиданием равна
Яобс= =1 —*' |
Pin, а) |
1 — у."1 |
|
|
Я (П, я) + Р (л, я) х —-------- |
|
1 --X |
Рассмотрим первоначально случай, когда х = 1 (а =п), и по
кажем, что
P%L-P№c>0.
Разность вероятностей в этом случае равна
1 - |
р (/I, п)________ Г 1 |
Р (п, |
п) |
||
R (т, п) 4- тР (п, |
п) |
R (п, |
п) ] - |
||
|
|||||
= Р |
1 |
________1________ |
|||
R(n, п) |
R(n, п) + |
тР (л, п) |
|||
|
так как знаменатель первой дроби в квадратной скобке меньше знаменателя второй дроби.
Рассмотрим общий случай, когда к ф 1:
|
|
Я in, |
1 — |
У . п |
Г>(2) |
Pin, |
я) + Р (п, я) у. —--------— Я (л, я) у.т |
||
я) |
1 — |
У. |
||
‘ ож |
Я^бС- Я(п, |
я) |
|
1— у.т |
Я (л, я) + Р (Л, я) V.—--------
Для того чтобы эта разность была положительна, нужно, чтобы был положителен числитель последней формулы; покажем, что он положителен:
Я(п, |
|
1 _чт |
а)у.т — |
|
а)-\-Р{п, а)х------------ R(n, |
||||
1_ |
1— X |
|
||
[Я(/г, а)(1 —к) —чР(л, |
а)] = |
|||
= |
|
|||
= |
- Г^ |
Г [Я («, « ) - * / ? ( « - 1 , |
а)]. |
1_угп
Отношение —------ положительно при любом значении х. По-
кажем, что разность, стоящая в квадратных скобках, тоже поло жительна:
R ( n , a) — x R ( n — \ , a) = R ( n , |
а ) ------— R { n — 1, а) = |
|
|
|
П |
= Р ( 0. |
*) + J ]P (* . |
а) ( 1 — ~ ) > о, |
|
k=\ |
|
так как каждый член |
суммы .не. отрицателен |
188
Таким образом, мы показали, что при одинаковых парамет
р а |
(п, Ху р) |
система с ожиданием имеет большую пропускную |
||||
способность, |
чем |
система |
с отказами. |
Это достигается за счет |
||
увеличения времени нахождения заявки в системе, т. е. за счет |
||||||
того, что заявка будет ожидать в очереди. |
||||||
5.1.2. |
Рассматриваются опять две системы: 1) система с отка |
|||||
зами |
с параметрами п, X, pi и 2) система с ожиданием с пара |
|||||
метрами |
п, а, р2, тп. При этом p i> ji2, |
т. е. производительность |
||||
каждого |
канала |
системы |
с |
|
||
отказами |
больше, |
чем про y l t i y ( 2 i |
|
изводительность каждого ка |
л 0 ' А 0 |
\ { 1 ) |
|||
нала системы с ожиданием. |
7 1 } Х у |
Аа |
|||
Считая |
число |
каналов |
п |
|
|
и число мест в очереди т за |
|
|
|||
данными, |
определить, при |
/ х |
1 |
||
каких значениях |
параметра |
||||
X система с отказами будет |
г |
1 |
|||
1 |
|||||
иметь большую |
пропускную |
0 |
Л ’ |
||
способность, чем |
система |
с |
|
|
|
ожиданием, и наоборот. |
|
|
Рис. 5.1.2а |
Р е ш е н и е Рассмотрим характер изменения абсолютной пропускной спо
собности для СМО с отказами (Хс(1)) и для СМО с ожиданием (W2)) (рис. 5.1.2а) в зависимости от интенсивности потока зая вок X при условии, Ч Т О Р 1> Р 2 . При большой интенсивности по тока заявок X все каналы будут практически заняты и в этом
случае пропускная способность обеих систем будет определяться лишь общим числом каналов п и производительностью канала р. Так как по условию p i> p 2, а число каналов у них одинаковое, то при X— их> Я0(1)>^о(2). При X— *0 абсолютная пропускная способ
ность обеих систем будет тоже стремиться к нулю (заявок очень мало), а относительная пропускная способность для обеих си стем будет одинаковой:
lim Робе = |
lim Р^с = 1, |
Х-*0 |
Л -0 |
так как любая заявка будет принята к обслуживанию.
В примере 5.1.1 было показано, что в случае одинаковых па раметров (/г, X, х) пропускная способность СМО с ожиданием (при т > 0) всегда выше пропускной способности СМО с отка
зами. Следовательно, о г р а н и ч е н н о е (небольшое) увеличе ние производительности обслуживания канала СМО с отказами можно «компенсировать» увеличением числа мест т в очереди
для СМО с ожиданием, т. е. всегда можно подобрать такую про изводительность канала обслуживания p i> p 2 для системы с от казами, когда пропускная способность обеих систем будет оди наковой.
189
На рис. 5.1.26 показаны |
графики абсолютной |
пропускно1”1 |
|||||
способности СМО с отказами |
(Яо(|)) и СМО с ожиданием |
(>.о<20 |
|||||
при значении параметров п=5; т = 5, pi = 2, |
рг=1- |
Анализируя |
|||||
графики на рис. 5.1.26, приходим к выводу, |
что при Х<А .*«5 |
||||||
|
пропускная |
способ |
|||||
|
ность СМО с ожнданй* |
||||||
|
ем |
немного |
выше |
за |
|||
|
счет того, что у нее есть |
||||||
|
места для ожидания не |
||||||
|
смотря на то, |
что |
про |
||||
|
изводительность |
кана |
|||||
|
лов СМО с ожиданием |
||||||
|
в 2 раза ниже произво |
||||||
|
дительности |
каналов |
|||||
|
СМО |
с |
|
отказами: |
|||
|
PI/M2=2. Однако при |
||||||
|
л> л *. |
когда |
СМО |
с |
|||
|
ожиданием почти пол |
||||||
|
ностью |
|
загружается, |
начинает сказываться наличие большей производительности СМО с отказами. Так, например, при Л=8 вероятность того, что в системе с ожиданием будет очередь, равна [см. (5.1.17)]
/ > 0 " = ^ у- - Т |
Г Г = |
° ’ 8 6 5 , |
||
т. е. практически система бу |
||||
дет |
полностью |
загружена, |
||
достигнув своей |
максималь |
|||
ной производительности, рав |
||||
ной «(.12=5, в то |
время |
как |
||
СМО |
с |
ожиданием еще не |
||
достигла |
своей |
максималь |
||
ной |
|
производительности, |
||
равной n,ui = 10. |
|
|
||
5.1.3. |
Определить то число мест в очереди т*, при котором |
|||
абсолютная пропускная способность СМО с ожиданием Ko^2)(m*) |
||||
отличалась бы от предельной абсолютной пропускной способно |
||||
сти системы Хо(2) (при неограниченном увеличении числа мест в |
||||
очереди т— изо) |
на заданную величину А. |
Р е ш е н и е Абсолютная пропускная способность СМО с ожиданием
(А,0(2)(т )) есть неубывающая функция числа мест в очереди т
(рис. 5.1.3а), имеющая определенный предел присяг — ос:
1о2)— Пт Хо2) (ш).
/Я-*-«о
190