Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Прикладные задачи теории массового обслуживания

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
12.79 Mб
Скачать

Вероятность противоположного события 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

Соседние файлы в папке книги