Материал: Прикладные задачи теории массового обслуживания.

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Соответствующая этому графу система дифференциальных уравнений имеет вид

^ % L = _X p0(/)

dp]}P~— ~ ~Ь нОР\ (О4- ^Ро(О-

at

Преобразование Лапласа для этих уравнений при отыскании за­

кона распределения величины Г0(ро(0) = 1; P i(0)= 0) запишет­ ся так:

($

Ро(5)~~ \**Pi (5) —55

(5~Ь ^ ~Ь Iх) Р\ (5) 1Ро(5)~

Решая эту систему относительно изображения p\(s), получим

Рi(s )= S2 +

\s

\s

5 (2Х + (1) +

Х2 (S + a) (s + Ь)

откуда

e - b t_ e-at

 

px[t)=\

а — b

 

 

где

 

 

а--

21 + ц — у

4Х<л +

(2.7.4)

^__2Х + [I -Ь \' 4X;J. 4- JJ.2

— 2

корни уравнения х2+ (2Л+ц)х + л2 = 0, взятые с обратным знаком.

Тогда искомая плотность распределения случайной величины Т0 будет

f о(0= ~ ~ ~ ~ —^Pi

e-bt —e~at

(*> П).

а — b

dt

 

Нетрудно убедиться в том, что условие

j u t ) d t = \

о

выполняется.

Действительно:

81

так как произведение корней равно свободному члену. Найдем математическое ожидание Т0:

-

>-2

^ te~b,dt — ^ te~atdt

M[T0] = t{

а — Л

 

 

 

 

м

 

 

к-

r j ____i_l=>2 JL+jL— J

I

Л

а — Ь

[ *2

a1- j ' a-b-

>. '

>.2

При отыскании закона распределения случайной величины Т\

начальными условиями будут ро(0)=0; pi(0) = l. В этом случае преобразование Лапласа для дифференциальных уравнений при­ мет вид

(s-г'-) AjOO-tJ-Pi (s) = 0;

(s -j- >'+!*) PI (s)— Xp0 (s) = s.

Решая эту систему относительно Pi{s), получим

s(s + l)______

Р\ (s) = - s2 + S (2), + jx) + X2

откуда

a— b

где а и b находятся из формул (2.7.4).

Следовательно:

/ i( 0 = x ? ,

l a — b

 

+ ± z L e-bi]

 

 

 

b — a

J

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

Tx будет равно

М [7^]—tx— ^ t f x (t) dt = X

 

 

oo

 

a— b

\ te~atdt 4-

 

 

о

J

 

 

 

 

 

о

 

+ - — -

^ te~bidt

[ a- ь op-

 

b — l

J _ l_ J_

b — a

.)

 

b — a b2 J

X 1 12 *

 

о

 

 

 

 

 

Сравнивая это выражение с выражением для математического ожидания случайной величины Г0, убеждаемся в том, что для любых К>0 и р > 0 значение ?0>Г Ь

82

Это объясняется тем, что система будет блуждать по состоя­ ниям х0 и Х\ в среднем больше, если она в начальный момент

находится

в

состоянии

х0. На

рис. 2.7.4в показаны

графики

плотностей распределения f0(t)

и f\(t) при Я = р=1.

 

 

2.7.5.

 

Через систему ПВО прорывается группа самолетов,

состоящая из т постановщиков помехи одного бомбардировщика.

Огонь сначала ведется по постановщикам помех, а затем по бом­

бардировщику.

 

Поток

пусков

 

 

 

 

ракет

можно

считать

пуассо­

 

 

 

 

новским

с

интенсивностью Я;

 

 

 

 

вероятность

 

поражения

одной

 

 

 

 

ракетой одного

 

самолета (лю­

 

 

 

 

бого)

равна

 

р.

 

После пораже­

 

 

 

 

ния самолета огонь

мгновенно

 

 

 

 

переносится на следующий са­

 

 

 

 

молет. Время

нахождения са­

 

 

 

 

молетов в зоне обстрела ПВО

 

 

 

 

равно

т. Определить

вероят­

 

 

 

 

ность того, что за время обстре­

 

 

 

 

ла т бомбардировщик будет по­

 

 

 

 

ражен.

 

 

 

 

 

 

 

 

 

 

 

Р е ш е н и е .

 

 

 

 

 

 

 

 

Группа

налетающих само­

 

 

 

 

летов может быть в следующих

 

 

 

 

состояниях: Xi — сбито

i само­

 

 

 

 

летов

(/ = 0,

 

1,

2,..., т + 1 ) .

 

 

 

 

Граф состояний системы по­

 

 

 

 

казан

на

рис. 2.7.5а.

 

 

 

 

 

Величина

/ч

определяется

 

 

 

 

следующим

образом: К\ = ЯР.

 

 

 

 

Этот граф с точностью до обозначений совпадает с графом,

изображенным на рис. 2.7.36*

 

 

 

 

 

 

 

 

^ 7

I

| ^ 7

\

1

Следовательно:

 

 

 

 

 

'V

dpm-t-1 (') h (V)"1

 

 

 

 

 

 

 

 

 

**

 

.-М

 

 

 

 

 

 

 

 

 

 

dt

т\

 

 

 

 

 

 

Рис. 2.7.5а

 

 

 

 

 

 

 

 

 

 

 

 

где pm+\(t) — есть вероятность того, что к моменту t бомбарди­

ровщик будет поражен.

Тогда искомую вероятность найдем по формуле

ут-ы

dPт+1(О

d t = \ —R(jn, >ux),

dt

 

где R(m, а) — табулированная функция (см. приложение).

83

Лг.пример, при Х=1 ------;

т=16 мин; /7 = 0,5; т = 10 получим

 

 

мин

 

 

 

(х) =

1—/?( 10,8)==1 —(1 —0, 184)=

0Э184.

2.7.6.

Граф

состояний

системы приведен

на рис. 2.7.6а. Из­

вестна интенсивность пуассоновского потока событий 'k(t), пере­

водящего систему из состояния

х$. Известно

также,

что при

 

выходе

системы

из со­

 

стояния х0 она

с

веро­

 

ятностью

 

оказывает­

 

ся

в С О С Т О Я Н И И

X i

(/= 1,

 

2,

...,

п).

 

Определить

 

а)

интенсивность

пуас­

 

соновского

потока

со­

Рис. 2.7.6а

бытий,

 

переводящего

систему

из

состояния

 

х0

в состояние

х\

и

б) вероятность pi(t) того,

что к моменту

t

система

будет в со­

стоянии x it если в момент t = 0 система была в состоянии х0.

О т в е т

 

 

 

 

 

t

-

| X (х) dx

 

а) >-|(0 = Ч 0 Йо

б)/7,(/)= ( >ч(~)е

0

d-.

 

О

 

 

 

При л(/) = \ = const получим:

 

a)

\ l= \ql) б)

Pl(i) = q,[ 1— е~и ].

2.7.7.

Граф

состояний

системы имеет вид, показанный на

рис. 2.7.7а. Требуется найти закон распределения времени Г0,1 нахождения системы в состояниях х0, хх в стационарном режиме,

если интенсивности всех потоков равны 1----

сек

Ре ш е н и е

Всоответствии с графом состояний система из группы со­ стояний х0ч Х\ может попасть только в состояние х2. Попав в со­

стояние х2, система с вероятностью 7г перейдет в состояние х0,

либо с вероятностью 72 — в состояние х3 (а оттуда в состояние *i). Следовательно, блуждание в группе состояний х0, Х\ с равной

вероятностью 72 может начаться из состояния х0 или из состоя­

ния Х\ . Граф состояний подсистемы X, с помощью которого мож­ но найти искомый закон распределения времени Г0,1 блуждания по состояниям х0, хи представлен на рис. 2.7.76.

84

Этот граф с точностью до обозначений интенсивностей пото­ ков событий совпадает с графом, изображенным на рис. 2.7.46. Следовательно, и решение с точностью до обозначений будет таким же, как и в примере 2.7.4.

Таким образом, если блуждание начинается из состояния л*0,. то плотность распределения времени Тол будет равна

_ з - УГ t

_ з + /5~

t

m = e- — -—

— -—

( < > о ) .

V 5

 

Если блуждание начинается из состояния х\, то плотность рас­ пределения времени Тол будет равна

3 + V 5

3 - / 5

К — Ч - 3- # )

т = -

{t > 0).

Уъ

Сучетом того, что блуждание в состояниях х0\ х{ может с веро­ ятностью 72 начаться из состояния х0 и с вероятностью 7г — из состояния Х\, получим следующее выражение для плотности, распределения времени блуждания в состояниях л:0; Х \ \

 

/о,1 (0 = 4" 1/о (*) +

f i

=

 

 

_

_ 3 +

t

_

 

(t> 0).

= —

( / 5 - 1 ) 6

2

- f ( / 5

+

l)e'

4 У

5 L

 

 

 

 

 

8S

График этой функции представлен на рис. 2.7.7в. Математическое ожидание времени пребывания в состояниях

-т0, Xj будет равно

/о,1= \ i/o.i{t)dt=-Y сек.

О

2.7.8. Автомашина при ее эксплуатации может находиться в следующих состояниях:

 

— исправна;

 

A'I — неисправна, проходит осмотр, ко­

 

торый проводится с целью опре­

 

деления вида ремонта;

 

а'2 — неисправна, проходит капиталь­

 

ный ремонт;

 

а'з — неисправна, проходит средний ре­

Рис. 2.7.8а

монт;

лг4 — неисправна, проходит текущий ре­

 

монт.

Граф состояний автомашины показан на рис. 2.7.8а. Среднее время межремонтного пробега равно 10. Среднее время осмотра машины равно t\. После осмотра автомашина подвергается капи­ тальному ремонту с вероятностью q2i среднему ремонту с веро­ ятностью <7з и текущему ремонту с вероятностью (<72+<7з4-<74= = 1). Среднее время проведения капитального ремонта равно t2,

среднее время проведения среднего ремонта равно /з, среднее

время проведения текущего ремонта равно /4. Определить веро­ ятность того, что машина будет исправна (для стационарного режима), вероятность того, что время простоя машины будет не более t и среднее время простоя машины.

Р е ш е н и е

Найдем интенсивности потоков, переводящих автомашину из состояния в состояние:

/ OA =

~ Zr ;

 

<7? .

__

Яз .

\

Я 4 .

м , 2 =

__

У

/•1,3 =

------- I

Л 1 ,4 =

— ,

 

*0

 

и

 

 

/ 1

 

h

л2,0 =

1

 

1

 

 

 

 

 

— ;

^‘З . о =

__

У

' • 4 , 0 =

—

 

 

 

h

 

h

 

 

и

 

 

Считая, что все потоки простейшие с выше указанными интенсивностями, найдем вероятность р0 того, что машина будет

исправна для стационарного режима. Решим соответствующую систему однородных алгебраических уравнений:

о = — ^0,1Ро+ *'2,0Р2“Ь ^3.0Pz ^4,0^1>

86

Источник: https://tut-files.ru/previewfile/83058