
книги из ГПНТБ / Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие
.pdfТак как G(w; 0)— 1, то в результате интегрирования получаем
In О (и; |
i |
1 |
= a i 2 ^ s ( 7i — l)s — (е8М- I) |
||
|
S= |
1 |
Исключая из этого равенства параметр rj с помрщью (20.31), для искомой производящей функции получаем следующее выражение:
i
G(u\ t) = exp « i 2 ^ ( w - 1)5 (1— |
(20.4(5) |
S = 1 |
|
Зная производящую функцию, вероятность Pk(t) нахождения системы в момент t в состоянии Ск можно определить по формуле
|
|
|
|
|
1 |
|
dkG(u; |
t) |
|
(20.47) |
|
|
|
|
Р |
к |
k\ |
|
дик |
u=0 |
|
|
|
|
|
|
|
(k = |
0, 1, |
...). |
|
|
|
||
В частности, |
при k — О для вероятности |
P0(t) |
того, что в момент t |
||||||||
в системе нет ни одного требования, находим |
|
|
|||||||||
P0(t) = |
G (0; |
t) = |
exp |
|
|
|
( - D |
8 <7s(l - e ~ ^ ) |
. (20.48) |
||
|
|
|
|
|
|
|
S = 1 |
|
|
|
|
Дифференцируя выражение (20.46) по и, получаем |
|
||||||||||
dG{u.\ t)_ __ |
1 |
|
|
|
|
|
|
G(u; 0- |
(20.49) |
||
du |
L |
s=l |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
||||
Тогда |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
P i ( t ) = |
^ ^ G (u; |
*) |
|
^ |
|
( - |
i ) 5- 1?sd |
sM) |
PoV). |
||
|
|
|
u=0 |
|
S = 1 |
|
|
|
|
||
Согласно |
(20.41) и (20.44) |
|
|
|
|
|
|
|
|||
|
|
2 |
( - 1 ) 5- 1Чъ— Яо |
|
(0) ~ Qо = I) |
|
|||||
|
|
S = 1 |
|
|
|
|
|
|
|
|
|
поэтому окончательно |
|
|
|
|
|
|
|
|
|||
|
|
Л ( 0 |
= «1 |
1 + |
2 |
|
( - D |
'^ e - 4”* |
Po(t). |
(20.50) |
|
|
|
|
|
|
s~l |
|
|
|
|
|
Последнее выражение можно также получить с помощью первого уравнения системы (20.39), которое при известной функции Р0 (t) записывается в виде
(2 0 .5 1 )
Г1
150
Аналогично можно получить явные выражения для вероятно стей Pk(t) при любых значениях k. Для математического ожида
ния у (t) числа приборов, занятых обслуживанием требований, с по мощью (20.49) находим
(20.52)
где
4i = 2 /> j. |
(20.53) |
i-i |
|
Система дифференциальных уравнений (20.16) может быть ре шена в явном виде и в некоторых других частных случаях. Имея аналитические выражения для вероятностей Рк(£) (^ = 0, 1, ...
..., п -(- от), можно исследовать характер изменения различных показателей эффективности при неустановившемся режиме функ
ционирования |
системы |
массового обслуживания. Вероятности рk |
(k — 0, 1, .. ., |
п-\-т) |
состояний системы при установившемся ре |
жиме ее функционирования, если только такой режим существует, определяются как решение следующей системы алгебраических уравнений:
k
- («1 + k) рк- f {k + |
1 )рк+х + |
a, 2 rjPk-j = |
0 |
|||
|
|
|
|
|
j=l |
|
(k = 0, |
1,. .. , |
ti — 1); |
|
|||
(a! + Л + SpiJ/Vt-s + |
[n -f- (s -f- |
1) pj] /?n + s + l |
-f- |
|||
|
n-fs |
|
|
|
(20.54) |
|
+ |
a l 2 |
Гj/7 I1+s_ j = 0 |
|
|||
|
|
|||||
( s = |
0 , |
1 , . . . , |
о т — 1 ) ; |
|
где
|
V |
T — |
(20.55) |
Искомые вероятности связаны равенством
П+Ш
(20.56)
151
В частных |
случаях |
система (20.54) упрощается. |
Если, напри |
мер, т = 0, / > |
п и гг = |
1, то вместо (20.54) будет: |
|
|
(ai + |
k) рк — (k - f - 1) рк+1 |
(20.57) |
|
(k = |
0, 1, . . . , п — 1); |
=(1 — Р п ) .
Сучетом (20.56) из (20.57) находим:
Рп- |
|
; - Ро ■ |
п\ |
|
|
|
П ( * + |
ai) |
|
||
|
|
|
|
||
|
|
к—1 |
S-1 |
|
|
|
Роа1 |
a, |
п\ |
(20.58) |
|
|
|
||||
/V |
к\ |
П ( * + « ! ) = - I P |
ai) |
||
|
|
s = i |
k\ П |
(s + |
|
|
|
|
|
|
|
|
|
|
s=k |
|
|
|
(Л = 1, 2 , . . . , П— 1). |
|
|
Из (20.35) и (20.36) для этих вероятностей получаются другие выражения в виде
А — |
l)s С„ - к s _|_ £ _(_ а) |
(20.59) |
|
||
|
s=0 |
|
|
(й = 0, 1, . .. , п). |
|
Зная вероятности рк |
(к = 0, 1, ..., п -f- т), по аналогии с дру |
гими системами можно определить различные показатели эффек тивности установившегося режима функционирования рассматри ваемой системы массового обслуживания. Вероятность полной за грузки приборов обслуживания и вероятность наличия очереди оп ределяются формулами:
m n — 1
Л..э= 2 / » „ + . = ! - |
2 рк; |
(20.60) |
s-О |
к — 0 |
|
Рп+$ ^ ■ S a - |
( 20.61) |
|
|
к=0 |
|
Математическое ожидание числа приборов, занятых обслужива нием,
П
У = |
2 |
крк 4 - п Р оч. |
(20 .62) |
|
к=1 |
|
152
В частных случаях выражения для у могут быть получены с по мощью формул вида (20.38) и (20.52) при t оо. Коэффициент загрузки приборов, т. е. вероятность того, что прибор обслужива ния занят,
£3аг = “ - |
(20.63) |
Математическое ожидание числа требований, ожидающих начала
обслуживания,
m
z = |
2 spn+s. |
(20.64) |
|
|
S=1 |
|
|
Математическое ожидание числа требований в системе |
|
||
7 = |
г/+ 7 . |
(20.65) |
|
Обозначим через | число требований,поступающих |
в систему |
||
в единицу времени. Математическое |
ожиданиеэтойслучайной ве |
||
личины |
|
|
|
i |
|
|
|
И=7 2 |
/0 = |
4 - |
(20.66) |
j=i |
|
|
|
Произведение г/ц равно математическому ожиданию числа тре бований, обслуживаемых системой в единицу времени, поэтому от
ношение Уу к совпадает с вероятностью полного обслуживания любого требования, т. е.
|
|
^п.обсл = Y |
• |
(20.67) |
|
Вероятность того, что |
требование будет обслужено частично, |
||||
|
|
РЧа с т = Л . |
(20.68) |
||
Вероятность того, |
что |
требование покинет очередь, |
определяется |
||
с помощью равенства |
|
ZV |
|
|
|
|
|
|
|
|
|
|
|
Рп.о |
т |
|
(20.69) |
Вероятность отказа в обслуживании |
|
|
|||
Ротк'— 1 |
' (Рп.обсл -}" Рчаст) = 1 " |
(20.70) |
|||
|
|
|
|
* 1?1 |
|
Вероятность Рпеп |
непопадания |
требования в систему |
можно найти |
||
с помощью соотношения |
|
|
|
||
Рнеп |
1 |
(Рп.обсл |
| Рц |
ггЬ Рп.оч). |
(20.71) |
153
Для расчета этой вероятности можно также использовать равен ство
|
|
|
|
р |
— |
|
|
(20.72) |
|
|
|
|
А НР(1 |
|
|
||
|
|
|
|
|
6 |
|
|
|
где |
ц — математическое |
ожидание |
числа требований из группы, |
|||||
не попадающих в систему в единицу времени. |
|
|
||||||
|
Для определения ц воспользуемся формулой |
|
|
|||||
|
|
n +m |
|
|
|
|
||
|
fl — |
2 |
|
Рп+m -k^ (VCn+m-k) . |
|
(20.73) |
||
|
|
k=u |
|
|
|
|
|
|
где |
M(7i/Cn+m_ u) — математическое |
ожидание |
случайной |
вели |
||||
чины ц, вычисленное |
в |
|
предположении, что система находится |
|||||
в состоянии Сп+га_ к, т. |
е. |
|
что |
в системе имеется |
п-\-т — k |
требо |
ваний. При указанном условии в систему одновременно может по ступить не более k требований, а остальные требования из группы получают отказ. Следовательно,
М (7]/Cn+m_ k) = Ь |
2 |
(/ - к) г}. |
(20.74) |
|
Тогда |
|
j=i*+i |
|
|
|
l |
|
|
|
|
n-f m |
|
|
|
^ |
2 Рп+т-к |
2 |
а - к ) г }. |
(20.75) |
|
к=0 |
)=к+1 |
|
|
Пример 20.1. Предназначенная для обработки групповых сооб щений информационная логическая машина имеет пять одинако вых приборов обслуживания. Поток групп сообщений простейший с интенсивностью 2,5 группы в секунду. Каждая группа содержит три сообщения. Если все приборы заняты, то сообщение теряемся. Математическое ожидание времени обработки каждого сообщения равно 0,2 сек.
Определить показатели эффективности установившегося ре жима функционирования данной информационной логической ма
шины. |
|
|
По условию п,— 5; т = |
0; / = |
3; ri = |
r2 = |
0; |
r3 = 1; |
|||
Р е ш е н и е . |
|||||||||||
v = Y = |
0; |
Л = |
2,5 1/сек; |
p ,= -i |
= 5 |
1/сек. |
Тогда |
pi = |
p; |
щ = |
|
= а = |
— — 0,5; |
(?i = 0. Соотношения |
(20.54) |
для вероятностей р к |
|||||||
(& = 0, |
Iх |
... , |
5) |
различных состоянии записываются в виде: |
|||||||
1, |
|||||||||||
|
|
|
P i = 0,5p0; |
2р2 = |
l,5pi; |
Зр3=2,5р2; |
|
|
|
||
|
|
Арл — 3,5Рз — 0,5р0; |
5р5 = |
4,5р4 — 0,5pi; |
|
|
|
||||
|
|
|
|
5ps = |
0,5 (р2 + Рз -f- pi). |
|
|
|
|
154
Т о г д а |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Р\= 0,5 р 0; |
р2 = 0,375ро; |
Рз = |
0,3125ро; |
|
|
|
|||||||
|
|
|
|
Ра= ОД484р0; |
Рь = 0,0836р0. |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
Подставляя |
эти |
выражения в условие нормировки 2 |
Рк— 1> по- |
||||||||||||
лучаем 2,4195ро=1, |
а |
|
потому |
Ро = |
|
|
к = 0 |
|
0,2066; |
||||||
|
0,4133. Тогда |
pi = |
|||||||||||||
р2 = 0,1550; |
рз = |
0,1292; |
|
р4 = 0,0613; |
р5 = |
0,0346. |
|
|
|
||||||
|
Вероятность |
полной |
загрузки |
приборов |
/эп.з= ръ — 0,0346. По |
||||||||||
условию т = 0, |
а потому |
Роч= 0 . Математическое ожидание числа |
|||||||||||||
|
|
|
|
|
|
|
|
|
_ |
5 |
|
|
|
|
|
приборов, занятых обслуживанием, |
у = 2 |
кРь — 1,3224. |
Коэффи- |
||||||||||||
|
|
|
|
|
|
|
_ |
|
|
k-i |
|
|
|
|
|
циент загрузки приборов |
&заг = |
= |
0,2645. Так как z = |
0, |
то ма |
||||||||||
тематическое ожидание числа сообщений в системе х = р = |
1,3224. |
||||||||||||||
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
Имеем Ц\ = |
2 |
jf} = |
3, |
тогда |
= |
7,5. Вероятность |
полного об- |
|||||||
|
|
|
3=1 |
|
|
|
_ |
|
|
|
|
|
|
|
|
слуяшвания |
сообщения |
Ра.обсл = |
^=- |
=0,8816. |
При |
этом |
Рчаст = |
||||||||
= |
Рп.оч= 0, |
а потоку |
вероятность отказа в обслуживании |
Ротк = |
|||||||||||
= |
1 — Рп.обсл = |
0,1184. |
Вероятность |
непопадания |
сообщения |
в си |
|||||||||
стему |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рнеп — 1 |
|
Рп.обсл = Ротк== 0,1184. |
|
|
|
|
Математическое ожидание числа сообщений из группы, не попа дающих в систему в единицу времени,
г) = Я (Зр5 -f- 2Ра-f- Рз) = 0,8890.
§21. СМЕШАННЫЕ СИСТЕМЫ
СГРУППОВЫМ ОБСЛУЖИВАНИЕМ КАЖДОГО ТРЕБОВАНИЯ
При анализе функционирования различных систем массового обслуживания предполагалось, что обслуживание любого требова ния производится только одним прибором. Существуют также си стемы, в которых каждое принятое требование обслуживается сразу несколькими приборами. В зависимости от организации об служивание одного требования могут производить сразу все при боры пли часть из них. Когда число требований в системе меньше числа приборов, но все приборы заняты, на обслуживание вновь поступающего требования может переключаться часть приборов или оно становится в очередь на обслуживание. Как и в рассмот ренных выше сметанных системах, требование остается необслуженным, если в момент его поступления в систему очередь содер-
155
»
жит максимально возможное число т требовании. В одних систе мах момент окончания обслуживания любого требования совпадает с моментом завершения его обслуживания каким-либо одним при бором, а в других системах требование до конца обслуживает каж дый прибор независимо от других приборов. Освобождающиеся приборы могут участвовать в дообслуживании требований или ожидать поступления очередных требований.
Рассмотрим некоторые системы массового обслуживания указан ного типа, в которых при наличии хотя бы одного требования за няты все п приборов.
Система с полной информацией о результатах обслуживания
Данная система массового обслуживания имеет п не обяза тельно одинаковых приборов обслуживания и т мест ожидания. Время Гц. обслуживания любого требования /-м прибором имеет
показательное распределение с параметром pj ( / = 1 , 2,-..., п). Пусть работа системы организована так, что при поступлении в си стему одного требования все приборы приступают к его обслужи ванию независимо один от другого, причем обслуживание счи тается законченным при окончании обслуживания требования любым одним прибором (имеется полная информация о результа тах обслуживания каждым прибором). Если в момент поступления требования приборы обслуживания заняты, то требование стано вится в очередь. Время Т-, ожидания начала обслуживания слу чайное, имеющее показательное распределение с параметром v. Когда время ожидания больше 7\, требование покидает очередь н
потому остается необслуженным. Если в |
очередп уже имеется |
m требований, то очередное требование из |
потока получает отказ |
в обслуживании. |
|
Обозначим через Гр случайное время обслуживания любого тре бования. При указанной организации работы системы справедливо
равенство |
|
|
|
|
|
|
|
' |
P (Tv. > t ) = P { T ^ > t , |
Гр2> t, . . |
. , |
7Vn> |
t). |
(21.1) |
|
Так как случайные величины Гр^ |
( / = 1 , 2, |
... , |
п) |
взаимно неза |
|||
впсимы и имеют показательное распределение, то |
|
|
|||||
|
Р(Г„ > ^ ) = e- ^ |
+14+- +,ln,t == е~*\ |
|
(21.2) |
|||
где |
|
|
|
|
|
|
|
|
»* = |
£ |
h- |
|
1 |
|
(21-3) |
|
|
j=i |
|
|
|
|
|
Следовательно, время Гр обслуживания каждого требования яв ляется случайной величиной, имеющей показательное распределе ние с параметром р. Но тогда система массового обслуживания
156
с полной информацией о результатах обслуживания эквивалентна рассмотренной в § 18 смешанной системе с одним прибором обслу живания, производительность р. которого находится с помощью ра венства (21.3). Состояние Ск при этом означает, что в системе на ходится к требований (к— 0, 1, ... , т + 1), Предельные вероят ности рк (k~0, 1, ..., т -{- 1) состояний системы определяются формулами (18.1) — (18.3) при п = 1, т. е.
т + 1 |
а |
|
|
Ро— 1 + 2 |
|
||
|
|
||
s=i П О + / Р ) |
|
||
“1+Vo |
г=о |
(21.4) |
|
(s = 0, 1, ... , |
т), |
||
|
па + т
1 =0
где
V
|
|
|
|
(21.5) |
Вероятность |
полной |
загрузки приборов |
А | . з = 1 — Ро', |
вероят |
ность наличия |
очереди |
Роч = 1 — Ро — Рь |
Математическое |
ожида |
ние числа требований, обслуживаемых: в произвольный момент вре
мени при установившемся режиме, |
у = |
1 — Ро. Вероятность |
обслу |
|
живания любого требования |
|
|
|
|
Л ,бс«= -£- = - Ц |
к |
^ |
- 2 ft- |
(21.6) |
а |
|
j=1 |
|
Математическое ожидание числа требований, ожидающих начала обслуживания,
m |
|
z ~ 2 ^ 1+s = "h” (1 Pl+m) — "Ъ- • |
(21-7) |
Справедливы также формулы (18.14) — (18.19).
Пример 21.1. Для отражения нападения воздушных целей на корабле имеется пять зенитных комплексов двух типов. Математи ческое ожидание времени обстрела каждой цели любым из трех комплексов первого типа равно 0,3 мин, а каждым комплексом второго типа — 0,2 мин. Вероятность поражения цели за стрельбу любым комплексом равна 0,5. Управление огнем организовано так, что все зенитные комплексы начинают обстрел одной цели в мо мент ее входа в зону обстрела. Если в указанный момент хотя бы один комплекс занят обстрелом, то очередная воздушная цель остается необстрелянной.
157
Определить вероятность поражения каждой воздушной цели при
установившемся |
режиме, когда |
интенсивность |
поступления |
целей |
|||||||||||||
в зону обстрела в минуту равна 5. |
число |
приборов |
обслуживания |
||||||||||||||
Р е ше н и е . |
В |
данном |
случае |
||||||||||||||
п = |
5. |
Мест ожидания нет, |
т. е. т — 0. |
Известно, что |
t^ |
— 0,3 мин |
|||||||||||
( / = 1 , |
2, |
3); |
ti>.l = |
0,2 мин |
(/ = 4, |
5). |
Тогда |
= |
^ |
\/мин |
|||||||
(у = |
1, |
2, |
3); |
рг= 5 |
1 /мин |
|
|
|
|
5 |
|
|
|
|
Так |
||
(/ = |
4, |
5); |
р = |
2 Р] = 20 1 /мин. |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
j=i |
|
|
|
|
|
как |
А = 5 |
1/мин, то |
а = — = |
0,25. Согласно |
(21.4) |
имеем: |
Ро = |
||||||||||
— (1 + |
а) |
1= |
0,8; |
Pi =■ 0,2. |
При |
этом |
у — 1 — Ро = |
0,2; |
Робсл = |
||||||||
= —- = 0,8. |
Искомая |
вероятность |
поражения |
каждой |
|
цели за |
стрельбу
/>= Л,бслП-(1 — Р)5] = 0,8(1 - 0 ,5 5) = 0,8 -0,6875 = 0,55.
Система с отсутствием информации о результатах обслуживания
Рассмотрим |
смешанную |
систему |
массового обслуживания |
с . групповым |
обслуживанием |
каждого |
требования, в которой |
имеется п одинаковых приборов и т мест ожидания. Время 7\L обслуживания каждым прибором любого требования является слу чайной величиной, распределенной по показательному закону с па раметром р. При этом вероятность успешного обслуживания каж дого требования равна р. Информация о том, успешно или без успешно закончилось обслуживание того или иного требования, в систему не поступает, а потому дообслужпванпе требований не производится. При занятых k приборах (Р = 0, 1, ... , п — 1)
вновь поступающее в систему требование начинают обслуживать независимо один от другого сразу п — k свободных приборов. Когда все приборы заняты, вновь поступающее требование стано вится в очередь на обслуживание, где одновременно может нахо диться не более т требований. Время Т» ожидания начала обслу живания является случайной величиной, имеющей показательное распределение с параметром v.
Состояние Ck (k = 0, 1, ..., п) |
означает, |
что |
обслуживанием |
занято k приборов. Прп состоянии |
Cn+S (s = |
1, 2, |
... , m) обслу |
живанием заняты все п приборов, а в очереди на обслуживание находится s требований. Общее число требований в системе при
состоянии Ск (k = 1, 2, |
. .. , п) может быть любым от 1 |
до k, а при |
|||
состоянии Cn+S — любым от s + |
1 до s + |
п. |
в том |
слу |
|
Смена состояния Ск |
(k — 0, |
1, .. ., |
п) происходит |
||
чае, если в систему поступает требование или если один из k |
при |
158
боров |
(при £ > |
1) |
заканчивает обслуживание. |
Состояние Cn+S |
|
( 5 = |
1, 2, |
... , т) |
система изменяет в том случае, |
когда в систему |
|
поступает |
требование |
(при 5 < т — 1) или один |
из п приборов |
заканчивает обслуживание, или одно из s требований покидает очередь, не дождавшись начала обслуживания. Поэтому
Гк = |
Х + |
£р |
(£ = |
0 ,1 ........ я); |
|
|
|
Tn+s = |
X + |
Щ - f |
sv |
(s— 1, |
2 , . . . , |
т — 1); |
(21.8) |
Тп+ш = пу. + тч. . |
|
|
|
|
|||
Из состояния |
Ск |
( £ = |
1, 2, |
. .. , п) |
система переходит в состоя |
||
ние Ck_i, если освобождается |
один из k приборов. Следовательно, |
||||||
|
|
Tk,k-i = kY |
( k = i , 2, |
п). |
(21.9) |
Поступление в систему одного требования при исходном состоянии
Ск (к = 0, 1, ... , п — 1) |
приводит к переходу системы сразу в со |
||||
стояние Сп, поэтому |
|
|
|
|
|
Тк,п = х |
(А = |
0, 1, |
П - 1). |
(21.10) |
|
Справедливы также следующие равенства: |
|
|
|||
7n+s,n+s+i = ^ |
(s = 0, |
1 , . . . , |
т |
1); |
| |
fn + s .n + s - 1 =7=«p |
+ s v |
( s = 0 , |
1 , . . |
. , /га). |
I |
С учетом полученных выражений для отличных от нуля коэф фициентов Yk и Ykj система дифференциальных уравнений (13.1) для искомых вероятностей записывается в виде:
Рк (*) = |
- (к + ^ ) Рк (t) + (к + 1) аРк+1 (t) |
|||
|
(k = |
0, |
1......... га — 1); |
|
|
|
|
+ |
k=0 ^ ( о + |
|
+ |
(л**- + v) Pn+I (t) ; |
( 21. 12) |
|
-^П+S (^) — |
O' + nY+ |
sv) Pn+s (t) + |
^n+s-1 (^) + |
+ [n-Y- + (S+ 1) v] Pn+s+1 (*")
(s = 1, 2 ,. . . , /га — 1);
P n+m (0 = — («Р + |
Pn+m(t) + XPn+m_ 1 (t). |
|
|
Если в системе очереди быть не может, т. е. |
т = 0, то общее |
||
число дифференциальных уравнений равно га + 1. |
Первые га из них |
||
совпадают с уравнениями из (21.12), а |
последнее записывается |
||
в виде |
|
|
|
^ ( f ) = |
- « i ^ „ ( * ) + x 2 |
p,(t). |
(2i.i3) |
159