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

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

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

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

Добавлен: 17.10.2024

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

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

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

боров

(при £ >

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



При этом

система

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

с (20.21), когда fxi =

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

искомые вероятности опре-,

деляются формулами (20.35) и (20.36), т. е.

 

P0{t) =

+

( -

1)5Сп

[1 _ e-(x+^)t];

(21.14)

 

 

 

s=0

 

 

 

рк (t) =

Сп " | V

I)3CiU

 

[ l - ^ + » + « ) M t ]

(21.15)

 

 

 

( * =

1, 2, ....

л),

 

где а == — . • (j.

Предельные вероятности определяются формулами (20 58) или

(20.59), т. е.

п\

Ро-

f l( s

+ «)

/V

Л + а ’

 

 

 

 

 

S— 1

 

 

 

 

 

k—1

 

 

are!

(21.16)

Рк =

П

(s +

°0:

 

 

 

s=0

 

А |П (* +

«)

 

 

 

 

 

 

 

s=k

 

или

( f t = l , 2 ,. .. , л - 1 )

 

 

 

 

 

 

 

 

 

 

а

(21.17)

 

 

 

 

 

 

(k =

0, 1,

...,

re).

 

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

k) Рк— (k + 1) Рк+i

— 0, 1, . . . , ге — 1 );

 

П-1

 

 

(я + л) ра — а к-0 Рк +

(л + Р) Pn+i ;

(а -j- re -f- sp) jPn+s — ^Pn+s—i

[л ~Ь (s ~Ь 1 ) Р] ft+s+i

 

(s = 1, 2, . . •, /ге — 1);

 

(л -)■' Л/Р)/ ’n+m — аРп+т—11

где Р =

— .

 

 

V-

 

160


Из первых п уравнений следует, что

 

Ро

к — 1

 

 

 

 

 

 

 

П ( И - а )

(к— \, 2 ..........п).

(21.19)

Положим

k\

1=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я-S

aPn+s—1

 

ip +

5р) Pn+s

(21.20)

 

 

 

( s = 1 , 2,

т).

 

 

Тогда последние т уравнений из

(21.18)

принимают вид

as+1 =

as

( s = 1,

2 , . . . ,

т -

1); ат= 0.

 

Но тогда as =

О (s =

1 , 2, ..., т) , а потому

 

 

 

 

я

__aA i+ s —1

 

 

т. е.

 

 

Pn+S

n +

s$ ’

 

 

 

 

 

 

 

 

 

 

Р n+s

 

«Л .

(S = l, 2,

т).

( 21. 21)

 

 

 

 

П ( « + /Р)

 

 

 

 

 

 

j-i

 

 

 

 

 

 

С учетом

равенства £0 = 0

из

(и +

1)-го уравнения

(21.18) на­

ходим

 

 

П—1

 

 

 

 

 

 

 

 

пра

 

 

 

 

 

 

_

 

 

( 21.22)

 

 

 

^ Р к :

 

а

 

 

 

 

 

к=0

 

 

 

 

 

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

и (21.21) в соотношение

 

 

 

 

n -fm

 

 

 

 

 

 

 

 

2

а

= 1,

 

 

(21.23)

 

 

 

к—О

 

 

 

 

приходим к равенству

+т+2-

Рп= 1.

8=1 П (« +

/Р)

j-1

 

Тогда

 

Рп = i + i + 2

 

j=i

 

п\рп

 

Ро ~ 1Г=Т------

 

П (1+ а)

 

/=о

 

( 2 1 . 2 4 )

(2 1 . 2 5 )

11

161


Зная предельные вероятности рк (k = 0, 1, ..., п + т), можно рассчитать различные показатели эффективности рассматриваемой системы массового обслуживания. В частности, математическое ожидание числа приборов, занятых обслуживанием,

п

m

п

+

Рп

(21.26)

у = 2

л 2

Pn+s = 2 kPk + п 1 - 1

к-1

s=l

к=1

 

 

 

Математическое ожидание числа требований в

очереди

 

 

 

ГП

 

 

 

 

 

* = 2 « / W

 

 

(21.27)

 

 

S=1

 

 

 

Вероятность непопадания требования в систему совпадает с ра+т .

Вероятность ухода требования из очереди равна - у Суммируя эти

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

Л™ = />п+т + Т .

'■

(21.28)

Чтобы найти вероятность Рп.обсл полного обслуживания любого

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

 

Рп.обсл = 75 ( Л ) = 2 Р(Нк)Р(А/Нк).

 

(21.29)

к=0

 

 

Случайное событие А означает полное обслуживание требования. Гипотеза Нк (k = 0, 1, . .. , п — 1) означает, что в момент поступ­ ления требования в систему было занято k приборов, а потому это требование обслуживалось п k приборами. Гипотеза Нп— что при поступлении требования в систему все приборы были заняты,

но это

требование

обслуживалось

(одним прибором), а гипотеза

PIп+1 — что требование получило

отказ в обслуживании. Вероятно­

сти этих гипотез:

 

 

 

 

 

 

 

 

 

/3 (^ k )= P k

(* =

0,

1 ,

..., л - 1 ) ;

 

 

(21.30)

 

 

n+m

 

 

 

 

n—1

 

 

 

Р (Яп) =

s2= 0 Pn+s

Pqik1

 

k=02P\/L ^отк!

P (#„+,) = P0TK

Кроме

того,

P(A/Hn+l)= 0,

Р (Л ///п) = р, a

P(AfНк) =

1 - q"~k

(* = 0, 1 , . .. ,

n — 1 ), где <7=

1 — p.

 

 

 

Искомая вероятность полного обслуживания требования

 

 

 

п—1

 

 

 

 

/

п—1

\

 

 

р„.обсЛ= 2

p k ( i -< ?"-“) +

1 - Л ™ -

2 р к )р-

(21-31)

 

 

к = 0

 

 

 

 

\

к = 0

/

 

162