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

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

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

Для математического ожидания Wi (t) числа Wi (t) целей, обстре­ лянных за время t, согласно (14.37) имеем

 

« ■ = Т Т Т + ( Н ^ р - 11 - е"

=

 

 

 

 

 

 

= 0,25

+ 1 — е~‘).

 

 

 

 

Разность &wx (tk)=

Wi (tk) — ^i(^k-i)

равна математическому ожи­

данию числа целей, обстрелянных в промежутке времени от

i

до tk. Значения

функций Wi{t)

и

Ддо^)

при

tk— 4 - i = 1 мин

приведены в табл.

1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1

t

0

 

1

2

 

3

4

5

6

 

w ,(t)

0

 

0,408

0,716

0,988

1,245

1,498

1,749

 

Д®! (t )

 

0,408

0,308

0,272

0,257

0,253

0,251

 

Математическое ожидание числа обстрелянных воздушных це­

 

лей за 6 мин до (6)

1,75. При использовании приближенной фор­

 

мулы (14.40)

до =

1,5.

 

 

 

 

 

 

 

§

15. СИСТЕМА С НЕОГРАНИЧЕННЫМ ЧИСЛОМ

 

 

 

 

 

ПРИБОРОВ ОБСЛУЖИВАНИЯ

 

 

 

Когда число п приборов обслуживания достаточно велико, при­

 

ближенно можно считать, что в системе имеется неограниченное

 

число одинаковых приборов. В такой системе любое требование из

 

входного потока сразу начинает обслуживаться, причем вероятность

 

обслуживания Р 0бсл=1-

Число состояний системы п = со,

причем

 

состояние Ск означает, что обслуживанием требований занято к

 

приборов (к = 0, 1,

...).

Вероятности

(0

(к =

0, 1, ...)

нахож­

 

дения системы в указанных состояниях являются решением си­

 

стемы дифференциальных уравнений

(13.23)

при п = со, т. е.

 

Рk (t) — O' +

 

^Р)

(0 + ^ k —1 (t) + Р 1) Рк+1 (О

(1 5 .1 )

 

 

 

 

= 1, 2 , .

. . ).

 

 

 

 

по

Если в момент t — Q все приборы обслуживания свободны, то на­ чальные значения искомых функций следующие:

Л>(0) = 1; Рк(0) = 0 (Л= 1, 2, ...)•

(15-2)

Получим решение системы (15.1), состоящей из бесконечного числа связанных между собой дифференциальных уравнений, т. е. найдем явные выражения для вероятностей ЛДО (& = 0, 1, ...). Для этого воспользуемся возможностью преобразования системы (15.1) обыкновенных дифференциальных уравнений к одному диф­ ференциальному уравнению в частных производных относительно производящей функции, определяемой формулой

С ( » Л ) =

2 ^ к ( / ) .

(15.3)

 

к;—о

 

Умножив k-e уравнение из

(15.1) на ик (£ = 0, 1, ...)

и про­

суммировав результат по всем возможным значениям k, приходим к равенству

2

«к/>; (t)= х(« — 1) 2

и*рк (t) +

к =

0

к

= 0

 

 

+ ^

2(k + l)ukPk+1(t) -

2

ku*Pu(t)

 

к = 0

 

 

к = 1

 

Имеем

 

 

 

 

 

 

^ P k( t) = ^ 9 . S u’ Ъ ;

 

к - 0

&

 

 

2

(k + 1 ) икРк+1 ( t ) ~

2

ku'pk( 0 =

к =

0

 

к = 1

 

 

.= (1 - «

) 2

кик~1ръ ( о = о

-

« )

t] .

 

к = 0

 

 

 

 

Следовательно, система обыкновенных дифференциальных урав­ нений (15.1) эквивалентна следующему дифференциальному урав­

нению в частных производных:

 

 

+

 

(15.4)

С учетом (15.2) из (15.3) находим начальное

значение искомой

производящей функции G(«; t) при ^ =

0 в виде

 

G(u; 0 ) = 1.

 

(15.5 )

111

Метод решения уравнения (15.4) изложен в § 9. Из сравнения

(15.4) с (9.22) следует,

что координатами х, у, г для (15.4) яв­

ляются t, и, G соответственно, а 5 =

0. При этом х(|, У )~ G(u; 0) =

= 1, а потому а ( и )= 1 . Так как f(t, и, G) = ц — 1), g(t,

«, G) =

— Х(и— 1)G, то уравнения (9.24)

записываются в виде:

 

da

dG

X(u — \)G,

(15.6)

dt р («

— 1 ); dt

причем интегральная кривая должна проходить через точку (0, т], £),

где г} = ыД=о, а С =

ОД=о-

следует, что

 

Из первого уравнения (15.6)

 

 

 

 

da

 

,,

 

а потому

 

 

Ч Г = Т =

 

 

 

 

 

 

 

 

 

 

In — 1 ) = pi + С ,

 

причем С =

1п (т) — 1). Тогда

 

 

 

т. е.

 

 

In (и - 1 ) = {Л + 1п(т]— 1),

 

 

 

м = (г| — 1 )

+ 1 .

(15.7)

 

 

 

Согласно

(15.6)

имеем

 

 

 

 

 

 

dG

XG

 

 

поэтому

 

 

da

 

 

 

 

 

 

 

 

 

 

 

 

G — Схе 11

,

 

причем 6Д — Се

11

.Тогда

 

 

 

 

 

 

 

■ (и— Г|)

 

(15.8)

 

 

 

G = Се ^

 

Согласно (9,24)

в (15.8) следует заменить £ на <й(т) ) = 1.

Следова­

тельно,

 

 

 

 

 

 

 

 

 

G =

( U- Ti)

(15.9)

 

 

 

e *

.

Исключив из (15.7) и (15.9) параметр ц, получаем следующее выражение для искомой производящей функции:

G(u; t) — exp — (и - 1)(1 е-*)

(15.10)

Известно, что производящая функция случайной величины X, рас­ пределенной по закону Пуассона, имеет вид

0 ( и ) = ехр [(и — 1 )л ]

( 1 5 .1 1 )

112

Из сравнения (15.10) с (15.11) следует, что случайное число X(t) приборов, занятых обслуживанием требований в момент вре­ мени t, имеет распределение Пуассона. Математическое ожидание этой случайной функции

 

x{t) = — ( Г - е - И ) .

 

(15.12)

 

 

t1

 

 

 

Вероятность Рк (t) того,

что в момент t обслуживанием требований

занято k приборов, определяется формулой

 

 

Pk(f) = M

2

] l e-7(t)

(6 = 0 , 1 , . . . ) .

(15.13)

Таким образом, получено решение (15.13)

системы

дифферен­

циальных уравнений (15.1)

при начальных условиях (15.2).

Предельные вероятности

рк (6 =

0, 1, ...)

состояний системы

массового обслуживания с неограниченным числом приборов можно

найти с помощью

(15.13), полагая

t = c o . Для математического

ожидания у ~ х ( оо)

числа приборов,

занятых обслуживанием при

установившемся режиме функционирования системы, из (15.12) на­ ходим

у = — = о.

(15.14)

I1

 

Вероятность рк того, что обслуживанием занято 6 приборов, при этом определяется формулой

(* = 0 ,1 ,...) . (15.15)

Пример 15.1. При работе электронной вычислительной машины

ее отдельные блоки выходят из строя с интенсивностью два блока

всутки. Эти блоки сразу поступают в ремонт, а на их место ста­ вятся резервные. Математическое ожидание времени ремонта каж­

дого блока равно двум суткам.

 

_

 

 

 

 

Найти

математическое

ожидание

x(t)

числа

блоков,

находя­

щихся в ремонте в момент t, вероятности

Ps{t)

(s =

0, 1 ,

2 ) того,

что в момент t имеется s неисправных блоков,

и вероятности рк

(k = 0, 1 ,

...)

того, что при установившемся режиме работы вычис­

лительной

машины

в ремонте

будет находиться

одновременно

6 блоков.

 

 

 

 

 

 

 

 

 

 

Р е ш е н и е .

Известна

интенсивность поступления

неисправных

блоков: 6 = 2

l/сут.

Математическое

ожидание

времени

ремонта

каждого блока 1^ = 2 суг,

поэтому

р =

4 г - = 0,5 l/сут, а =

-----= 4.

 

 

 

 

 

 

tf.

 

 

 

1*

8

113

Согласно (15.12) математическое ожидание числа блоков, находя­ щихся в ремонте в момент t (t в сутках),

лг(^) = 4(1 — e_0'5t).

Искомые вероятности Ps(0 согласно (15.13) рассчитываются с по-‘ мощью равенств:

 

Р0(*)=

 

P l(t) = x ( t ) e - ^ m = x {t)P 0(t);

 

 

 

р 2 (t)= Y

Й012^-7(t) = о,5*(t) Pi (t ).

 

 

Значения функций x(t)

и

Ps{t)

(s = 0,

1 , 2)

при различных зна­

чениях t приведены в табл. 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2

t, сут

0

1

2

3

4

5

6

7

8

9

10

x(t)

0

1,57

2,53

3,11

3,46

3,67

3,80

3,88

3,93

3,96

3,97

PoV)

1

0,208

0,080

0,045

0,031

0,026

0,022

0,021

0,020

0,019

0,019

Л ( / )

0

0,326 0,202

0,140 0,107

0,095

0,084

0,081

0,079

0,075

0,075

А (0

0

0,256

0,256

0,218

0,186

0,175

0,159

0,158

0,154

0,149

0,150

Вероятность рк того, что при установившемся режиме работы вычислительной машины в ремонте будет находиться одновременно k блоков, определяется с помощью равенства

 

 

 

 

 

 

 

(£ = 0, 1 , . . . ) .

 

 

 

 

Значения этих вероятностей приведены в табл. 3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3

k

0

1

2

3

4

5

6

7

8

9

10

11

12

Ръ

0,018 0,073

0 ,146

0,195

0 ,195

0,156

0,104

0,060

0,030

0 ,013 jo,005

0,002

0,001

114

§ 16. СИСТЕМА С ОЖИДАНИЕМ

Система с ожиданием имеет конечное число п одинаковых при­ боров обслуживания и неограниченное число место ожидания, т. е. т = со. Если требование поступает в систему, когда свободен хотя бы один из приборов обслуживания, то оно сразу начинает обслу­ живаться. Когда все п приборов обслуживания заняты, поступаю­ щие требования становятся в очередь и обслуживаются в порядке их поступления по мере освобождения обслуживающих приборов. Требования не покидают очереди, и потому все они обслуживаются.

Число состояний системы с ожиданием бесконечное. Состояние

Ck

(k — 0, 1, .. .,

п) означает, что обслуживанием занято k прибо­

ров,

а состояние

Cn+S (s = l, 2, ...) означает, что обслуживанием

занято п приборов и ожидает начала обслуживания s требований. Вероятности Рк (t) (k — 0, 1, ...) нахождения системы в различных состояниях являются решением системы обыкновенных дифферен­

циальных уравнений, которая получается из

(13.23) при т = со и

v = 0 (допустимое время ожидания начала

обслуживания равно

бесконечности). Эта система, состоящая из бесконечного числа уравнений, записывается в виде:

P'0(t) =

- \ P 0(t) + v.P1(t);

 

 

 

Рк (t) =

-

(X +

kv) Рк (t) +

1Рк- , (t) +

v. (k + 1) Рк+1 (t )

 

 

 

{k = 1, 2,...,

n— 1);

(16.1)

^ n + s 0 ) —

0 - +

raH-) P n+s ( 0 +

^

n + s — 1 0 )

^ ^ ^ n + s + i 0 )

(s — 0, l , . . . ) .

Если при t = 0 все приборы обслуживания свободны, то начальные значения искомых вероятностей Рк (/) следующие:

p о (0) = 1 ; pk (0) =

0 (£ = 1 , 2 , . . . ) .

(16.2)

Предельные вероятности рк =

Нш Рк (t) (k = 0, 1 , ...),

если они

 

t-*- оо

 

существуют, определяются с помощью равенств (13.39), (13.40) п

(13.42) при т — со и

0. Следовательно,

 

 

Cl

(16.3)

Pk= ffiPo ( * = 1 , 2, . . . , Л);

(16.5)

115

Имеем

 

 

 

 

а

-----------

при

а <

я;

1 - -п

 

 

 

2 п

 

 

(16.6)

s=0

с о

при

а >

п.

 

Поэтому при а < п расчетная формула для вероятности Ро записы­ вается в виде

Ро =

S

(16.7)

 

+ (я —<»)(л — 1)!

Если а > п, то ро =

0, а потому рк= 0 (k = 0 , 1, ...), т. е. в этом

случае не существует стационарный режим функционирования си­ стемы массового обслуживания с ожиданием. Условие а > п экви­ валентно условию к пц, согласно которому математическое ожи­ дание числа требований, поступающих в систему в единицу вре­ мени, не меньше числа требований, обслуживаемых всеми п приборами в единицу времени. При этом система массового обслу­ живания не имеет возможности обслужить поступающий поток тре­ бований, вследствие чего очередь на обслуживание непрерывно растет и при t-* со становится неограниченной.

При установившемся стационарном режиме функционирования системы, который существует при а < я, любое поступающее в си­ стему требование обязательно обслуживается, а потому вероятность обслуживания Р0бсл— 1- Для математического ожидания числа Y

приборов, занятых обслуживанием, согласно (12 .2 1 )

получаем

У — а.

(16.8)

Коэффициент загрузки приборов (вероятность того,

что прибор об­

служивания занят)

 

£ з а г = ^ .

06.9)

Пусть Рп.з — вероятность полной загрузки приборов обслужива­ ния при установившемся режиме функционирования системы. Дан­ ная вероятность определяется с помощью равенства

ею

(16.10)

8=0

Учитывая (16.6) при ос < я, получаем

р = ДРп

(1 6 .1 1 )

116

где

Рп = ^ Р о -

(16.12)

Обозначим через Тн.3 случайное время неполной загрузки прибо­ ров обслуживания, т. е. Тн.3— длительность промежутка времени от момента перехода системы из состояния Сп в состояние Cn_j до момента возвращения системы в состояние Сп. Числовые характе­ ристики этой случайной величины, в частности ее математическое

ожидание tH.3, не зависят от наличия очереди'в системе массового обслуживания с ожиданием и потому могут быть определены, как для эквивалентной системы с отказами с п такими же приборами

обслуживания. Чтобы найти tH.3, введем случайное время Гй.з пол­ ной загрузки приборов обслуживания эквивалентной системы с от­ казами. Случайная величина Тп.з является длительностью проме­ жутка времени от момента занятия всех п приборов системы с от­ казами до момента освобождения хотя бы одного прибора. Эта случайная величина будет больше фиксированного значения t, если оставшееся время обслуживания требования любым из п приборов больше t. Так как приборы обслуживания функционируют незави­ симо один от другого, время Тц обслуживания любого требования имеет показательное распределение с параметром р, а время дообслуживания имеет такое же распределение, как и Тti. , то

П

 

P ( T L > t ) = П Я ( Г , > 0 = г ”^ .

(16.13)

j=i

 

Из этого выражения следует, что случайная величина Т„.3 имеет показательное распределение с параметром ftp. Следовательно, ма­ тематическое ожидание времени полной загрузки эквивалентной си­ стемы с отказами

 

tП*

___ 1_

(16.14)

 

 

ftp

 

Отношение tn.з к сумме

t$.3 +

tn.3

равно вероятности рп полной

загрузки приборов обслуживания в системе с отказами, т. е.

 

 

/*

 

Рп,

(16.15)

 

П.З = ~ =

^п.з

4 " ^н.:

 

 

где

 

 

 

 

Рп

 

 

( 1 6 . 1 6 )

k=0 R-

117

Разрешая (16.15) относительно ^н.з, с учетом (16.14) получаем

tн.з — - Ч 4 - 1

(16.17)

ППРп

 

Для системы с ожиданием случайное время Т„.3 полной за­ грузки приборов обслуживания не совпадает с Гп.з, так как при на­ личии требований в очереди окончание обслуживания любого тре­ бования не означает освобождения прибора обслуживания. Отно­

шение математического ожидания ^п.з случайной величины Т„.3

к сумме ^п.з+^н.з равно вероятности Рп.з полной загрузки прибо­ ров обслуживания системы с ожиданием, т. е. справедливо равенство

 

tn.3

Г>

 

(16Л8)

 

 

* П.З•

1

in.3 Ц- Гн.З

 

 

 

Тогда

 

 

 

 

 

 

 

Т

I

Рп'3

 

(16.19)

 

^П.З Ь'Н.З*1

р

 

-- * П.З

 

Математическое ожидание числа Z требований, ожидающих на­ чала обслуживания, находится с помощью равенства

 

Z

SPn+s —Pn

s=l

(16.20)

 

 

s=1

 

 

 

При |л: |<

1 имеем

 

 

 

 

 

2

- - = £

2

x

d

1

1

dx 1 —

( I - - * ) 2

S=1

dx s=0

Следовательно, при a <

n

 

 

 

 

 

 

 

anpn

_ аРп.з

(16.21)

 

Z

(n — a)2

tl — a

 

 

Отношение z к математическому ожиданию числа требований, поступающих в систему в единицу времени, т. е. к к, равно мате­

матическому ожиданию t0)K времени Тожожидания начала обслужи­ вания (времени пребывания требования в очереди), ’г. е.

Т

— z —

ПРп

Рп-3

(16-22)

ож

X

р (га — а)2

~ р (га — a)

 

Для, длительности ожидания Гож, являющейся основной харак­ теристикой функционирования системы с ожиданием, определим функцию распределения

F{t) = P{T0&<t).

(16.23)

118

Обозначим через Ps( Гож > t) условную вероятность того, что время ожидания начала обслуживания будет больше t, если в моч мент поступления требования в систему там была очередь из s тре­

бований (s = 0,

1, ...)• Тогда в соответствии с формулой полной ве­

роятности

 

 

 

 

Р(Тож > t ) =

Hp„+sPAPoM>t).

(16-24)

 

 

s=0

 

Пусть Qi(t)

— вероятность

того, что за время

t после поступ­

ления в систему рассматриваемого требования обслужено I требо­ ваний. По условию требования обслуживаются в порядке их по­ ступления в систему. При наличии очереди из s требований время ожидания начала обслуживания вновь поступающего требования будет больше t, если за это время обслужено нс более s требова­ ний. Следовательно,.

 

S

(16.25)

P,(T0A> t )

= '2Ql(t)

 

1=0

 

(s==0,

1, ...).

 

При наличии очереди поток обслуженных требований

простейший

с интенсивностью п\х, поэтому

 

 

 

 

(16.26)

Подставляя (16.26) в (16.25), находим

(16.27)

Тогда

P(T0* > t )

Меняя порядок, суммирования, получаем

п

119

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