Файл: Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 17.10.2024

Просмотров: 200

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

требования является случайной величиной, имеющей показательное распределение с параметром ц. При занятых приборах обслужива­ ния вновь поступающее в систему требование становится в очередь на обслуживание. Однако очередь не может содержать более т тре­ бований. Когда в системе имеется п -f- т требований, все очередные требования из входного потока получают отказ и не обслуживаются. Пусть, кроме того, любое требование, ожидающее начала обслужи­ вания, покидает систему и потому остается необслуженным, если время его пребывания в очереди превосходит допустимое значение Г». Время Тч также будем считать случайной величиной, имеющей

1

т

показательное распределение с параметром v — -=r где

t v =

= М{ Тч) — математическое ожидание времени пребывания требо­ вания в очереди.

Частными случаями введенной системы являются следующие си­ стемы массового обслуживания:

система с отказами;

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

система с ожиданием;

система с ограниченным числом мест ожидания;

— система с ограниченным временем пребывания требования

вочереди.

Всистеме с отказами (потерями) ограниченное число п прибо­ ров обслуживания и нет мест ожидания, т. е. пг — 0. В системе вто­ рого типа число п приборов обслуживания настолько велико, что

можно считать п = со. В системе с ожиданием ограниченное число п приборов обслуживания и неограниченное число мест ожидания, т. е. tn = оо, а допустимое время ожидания начала' обслуживания

каждым требованием из очереди бесконечное,

т. е. v =

0.

Система

с

ограниченным числом мест

ожидания отличается

от

системы

с

ожиданием только тем, что

в этом случае

m — конечное число.

В системе с ограниченным временем пребывания требования в оче­

реди может не быть

ограничения

по числу мест ожидания,

т. е.

в этом случае m —- любое целое число или пг = со.

 

Дифференциальные

уравнения

для вероятностей Рк (0

{& =

— 0, 1 , .. ., п + т) применительно к частным видам систем массо­ вого обслуживания получаются при указанных значениях чисел я, т и коэффициента v из дифференциальных уравнений для введен­ ной системы. Получим дифференциальные уравнения, описываю­ щие функционирование этой общей системы массового обслужи­

вания.

Смена состояний рассматриваемой системы в любой момент мо­ жет произойти из-за поступления очередного требования, при окон­ чании обслуживания любого требования и (при наличии очереди) вследствие ухода требования из очереди. Так как номер k состоя­

ния Ск

совпадает с числом требований в системе, то переход из со­

стояния

Ск в C k + j при j

1 за малое время At возможен только

92


при поступлении в систему не менее j требований. Вероятность по­ ступления ровно / требовании при стационарном потоке опреде­ ляется формулой

Pj {t, t + At) = - Z f - е~ш .

(13.4)

гГак как

 

е ■*«=-. 1

(13.5)

то вероятность (13.4) пропорциональна /-й степени малого проме­

жутка Д(. Переход из состояния Ск в Ck+J-

за время At возможен

и другим способом, когда,

например,

за

указанный

промежуток

времени закончат обслуживание I пз

к (k

п) приборов,

а

в си­

стему поступит не /, а / + j требований. Вероятность

такого

псхо

да пропорциональна At в степени

 

В общем случае вероят­

ность Pk.k+j (t, t -f- At)

перехода системы пз состояния

Ск

в

Qi+j

за малый промежуток времени -At можно записать в виде

 

 

p k.u+j(P

* +

д о = -

^

Ч

1 + о (ло ] .

 

 

(13.6)

С помощью этого выражения и (13.2) получаем, что временная

плотность вероятности 7;k

равна пулю при £ > / - ) -

1 ,

т. е.

 

 

 

 

' Т/k =

0

 

 

 

 

(13.7)

(k = I —j—2, l —f-3, ... ,

и tti]

l

0, 1 ,

..., n

in — 2 ).

 

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

Чк.к+i находим

 

 

 

 

 

"(к,к+1 =

А

 

 

 

 

(13.8)

(к — 0, 1 , .. ., я - f т — 1 ).

 

 

 

 

Коэффициенты 7к,к-н

при

1,

характеризующие переходы

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

Ск в Cj.+j, зависят только от интенсивности

потока А, а от параметров ц и v не зависят. Аналогично доказы­ вается, что временные плотности вероятности перехода системы из состояния Ск в состояния с меньшим номером не зависят от пара­ метра А, а потому при получении соответствующих выражений для Тk,k—j при 1 можно просто не учитывать входной поток требо­ ваний.

Пусть номер k исходного состояния Ск меньше п, т. е. в системе обслуживанием занято k приборов, и очереди на обслуживание нет. С учетом сделанного выше замечания можно считать, что переход системы из состояния С к в C k _ j при 1 k за малое время возможен только в том случае, если за это время / приборов из к закончат обслуживание. Время обслуживания имеет показательное распределение с параметром д, а закон распределения оставшегося времени обслуживания при этом не зависит от того, сколько вре­

93


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

Р ( П < At) = F^(At)=

(13.9)

Так как требования обслуживаются приборами независимо одно от другого, то вероятность перехода системы из состояния Ск в C k_ j за малое время At может быть записана в виде

 

PKk-i(t,

=

 

 

 

 

+

0

( 1 3 . 1 0 )

и л и ,

ч т о т о ж е с а м о е ,

 

 

 

 

 

 

 

 

 

Pk,k-j(t, t -f- Л() = Cl (рЛ()] -f- 0[(A^jJ].

(13.11)

Подставляя (13.11)

в (13.2), находим

 

 

 

 

 

 

 

 

T/k =

0

 

 

 

 

(13.12)

 

(* =

0,

1, ...,

I - 2;

/ — 2,

3,

... ,

n);

 

 

 

 

 

Tk.k-l —

 

 

 

 

(13.13)

 

 

 

(6 =

1 , 2,

...,

n).

 

 

 

 

П у с т ь и с х о д н ы м

я в л я е т с я

с о с т о я н и е

C n+S

( s =

1, 2 , . . . ,

m ) . П р и

э т о м з а м а л о е в р е м я At п р а к т и ч е с к и н е в о з м о ж е н п е р е х о д в с о с т о я ­

н и е

C n+S _ j п р и / >

1,

а п о т о м у

 

 

 

 

 

 

 

 

 

Tn+s,n+s-i == о

 

 

 

(13 -14)

 

(/ =

2, 3, .. ..

л + s;

s = l ,

2 ,

... ,

т).

 

Переход из состояния Cn+S в Сп+8_! произойдет в том случае, если один из п приборов закончит обслуживание или одно из s требова­ ний, не дождавшись начала обслуживания, уйдет из очереди. По аналогии с (13.13) получаем

Tn+s,n+s-i —

Щ + sv

(1 3 .15 )

(s = l,

2 , ... , т).

 

Вероятность Ркк (t, t-\-At)

того,

что за малое время At

система

не изменит исходного состояния Ск при k О я, с точностью до сла­ гаемых, пропорциональных (At)2, равна произведению вероятностей

e-XAt и (g-iJ.At)k _

Первая из них совпадает с вероятностью того, что

за время At в

систему не поступит ни одного требования, а вто­

рая — вероятность того, что за время At ни один из 'k приборов не закончит обслуживания. Следовательно,

Ркк (t, t + At) = о (Д*)= 1 — (X - f fcji) At + 0 (At). (13.16)

94


Подставляя эту вероятность в (13.3), находим

 

 

 

 

7к =

^ +

 

 

(13.17)

 

 

(* = 0, 1,

П ) .

 

 

При исходном состоянии Cn+S

аналогично

получаем

P n +S,n + S (t,

t +

At) = e - W te -n |iA tg -* » 4 t

О

s

=

1

(A. -f- wp - f -

sy) At + О (Д^),

( 1 3 . 1 8 )

поэтому

 

7n+s — Я -b np, -|-sv

 

(13.19)

 

 

 

 

( s = l , 2,

.... m — 1 ).

 

 

Если в системе находится п + т требований, то очередное требова­ ние получает отказ в обслуживании и потому номер состояния Сп+т не изменяется. Следовательно,

Рп+ra,n+ m (t, t + At) =

e-n^tg-mvM + Q

=

= 1 — (Л|1 + л м )Д * + 0 ( Д 0 ;

(1 3 .2 0 )

Tn+ш

Пр -)- lfl\.

(13.21)

Выражения (13.17), (13.19) и (13.21) можно было также полу­ чить с помощью соотношений

П+

ГП

 

 

7к = 2

Tkj (k = 0, I,

п + т)

(13.22)

j=0

 

 

\

по найденным ранее коэффициентам fkj-

Подставляя коэффициенты 7к и 7fk в (13.1), получаем диффе­ ренциальные уравнения для вероятностей состояний системы мас­ сового обслуживания в виде:

Р '( 0 = - м > о ( * ) + !

Рк (t) = -

(X +

Рк (t) +

(0 + v- {k 4- 1) Рк+1 (0

 

 

 

( k = l , 2 , . . . , л - 1);

 

Рп' +s (t) =

- (X +

Лр +

ЗУ) Pn+S (*) + XPn+s_ 1 (*) +

(13-23)

 

4~ \пР +

(s +

1) v] Лы-S+I (t)

 

 

(s = 0, 1,. . . , от — 1 );

 

■^П+Ш(^) ~

+ mv) ^ П+Ш(0 + ^n+m -l (t).

 

95


Начальными условиями для этих уравнений являются значения искомых функций в начальный момент времени t = 0. Если, напри­ мер, при t = 0 в системе не было требований, то исходное состоя­ ние СоТогда начальные условия для системы дифференциальных уравнений (13.23) следующие:

Р0(0) = 1 ;

Pk(0) = 0 (k = \ , 2 , ... , л +

m).

(13.24)

Нетрудно видеть,

что система уравнений (13.23)

совпадает с си

етемой дифференциальных уравнений (9.8) для процесса размно­

жения н гибели, когда N =

п

m:

 

 

Xk =

/,

(k — 0,

1, • . . ,

п + т — 1)

;

Рк —

^

 

 

2............ и );

(13.25)

P„+s =

« р - f

sv

(s 1,

2, •. . , т ) .

 

В результате решения системы дифференциальных уравнений (13.23) находятся вероятности Рк (t) 1 , ..., я -f-m), пер­ вые п + 1 из которых характеризуют загрузку системы массового обслуживания, а т остальных — очередь на обслуживание в любой момент времени t. Данные функции при любом / связаны условием нормировки

П+Ш

 

2 ^ ( 9 = 1 .

(13.26)

к=0

 

Если п-\-т= со, то при любых параметрах К, р и v равенство (13.26) может не выполняться. Достаточное условие существования единственного решения системы (13.23), которое удовлетворяет со­ отношению (13.26), записывается в виде (9.58), т. е.

1 Ё т ) = “ -

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

(яр + sv)

(13.28)

Когда бесконечным является число приборов обслуживания! т. е.

п — оо, условие (13.27) выполняется, так как 2

^' ( Т” ) =С0, ^сли

к=1

\ ' '

число т мест ожидания не ограничено, то при v Ф 0 условие (13.27) также выполняется, так как

 

+

= со .

(13.29)

i = i

s = i

i = i V 1

 

96