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

книги из ГПНТБ / Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие

.pdf
Скачиваний:
27
Добавлен:
22.10.2023
Размер:
12.12 Mб
Скачать

Так как 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,

... ,

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

Соседние файлы в папке книги из ГПНТБ