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

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

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

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

Добавлен: 17.10.2024

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

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

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

Зная вероятности Рь (t) (6 = 0, 1, .. ., п), можно проанализи­ ровать различные показатели эффективности функционирования системы массового обслуживания с отказами, часть из которых вве­ дена в § 12. В частности, характеристикой полного простоя си­ стемы в любой момент времени t является вероятность Po(t), а Рп(0 совпадает с вероятностью отказа в обслуживании. Матема­ тическое ожидание числа Y(t) приборов, занятых обслуживанием требований, в любой момент времени t находится по формуле

y { t ) = ' Z k P k{t). k=0

Коэффициент загрузки приборов

63аг (t) = 1 - у Ц)

( 1 4 . 8 )

( 1 4 . 9 )

Если время работы системы в два-три раза превосходит матема­ тическое ожидание времени обслуживания каждого требования

т

1

г^ ~ — , то, как показывают расчеты, можно считать, что переходIх

ной процесс закончился и происходит установившийся режим функ­ ционирования системы. Для вероятностей /?к (6 = 0, 1, .. ., п) на­ хождения системы в различных состояниях при этом справедливы

расчетные формулы (13.39) и (13.42)

при т — 0, а потому

P k = - j f P o ( 6 = 1 ,

2 , . . . , я ) ,

( 1 4 . 1 0 )

Ро

dd

а

 

( 1 4 . 1 1 )

 

V

6!

 

 

k=*

 

 

Выражения (14.10) и (14.11) называются формулами

Эрланга,

по имени впервые получившего их датского ученого А. К. Эрланга. При выводе этих формул предполагалось, что входной поток требо­ ваний простейший стационарный, а время обслуживания каждого требования распределено по показательному закону. Более деталь­ ные исследования, произведенные Б. А. Севастьяновым, показали, что эти формулы справедливы при любом законе распределения случайного времени обслуживания каждого требования с математи­

ческим ожиданием t[L, если принять р .= -= -. Формулами (14.10) и

(14.11) можно пользоваться при приближенных расчетах и в том случае, если входной поток требований не простейший. Замена дей­ ствительного стационарного потока требований простейшим с той же интенсивностью X обычно мало влияет на величины вероятно­ стей Рк(6 = 0, ! , . . . , « ) и на показатели эффективности, опреде-

101


ляемые через эти вероятности. Указанная замена возможна и при приближенных расчетах вероятностей (0 нахождения системы в различных состояниях для неустановившегося режима функцио­ нирования. Однако время обслуживания любого требования при этом должно быть распределено по показательному закону.

Вероятности (к = 0, 1, ... , п) могут быть рассчитаны по фор­ мулам (14.10) и (14.11) или найдены с помощью специальных таб­ лиц, составленных в зависимости от параметра а и числа п.

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

Л>бсл=1-/>п. (14.12)

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

У =

(14-13)

k=0

 

Так как согласно (12.17)

 

Р о бел = 4а- .

(14.14)

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

 

У = *(1 — Р п ) -

(14.15)

Коэффициент загрузки системы кзат, совпадающий с вероятностью того, что любой прибор обслуживания занят,

Аааг=“^------- ^ -(1— /»п) -

(14.16)

Обозначим через живания, т. е. Гпр — прибора до момента

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

a tnp = M (Tnр). Математическое ожидание времени обслуживания

каждого требования одним прибором равно tv. = — . Отношение Щ к

 

 

 

t1

 

£,1 + t

является вероятностью занятости прибора обслуживания

и потому совпадает с k3ZT, т.

е. справедливо равенство

 

 

 

__

У

(14.17)

 

tv-

^пр

п

 

 

102


Для математического ожидания времени простоя прибора обслужи­ вания при этом получается следующее выражение:

tnp= t j

~

- 1

- 1

(14.18)

 

 

а (1

- Р , i)

 

Пусть случайная

величина Тс — время

нахождения

требования

в системе массового обслуживания. Для определения функции рас­

пределения

Fc (t) = Р(Тс < t) этой случайной величины

восполь­

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

 

Р{ТС<Ь) =

Р {Н ,)Р {ТС< Ц Н , ) + Р { Н 2)Р

(14.19)

где гипотеза Н\ означает, что обслуживанием заняты все п прибо­ ров, а Н2— свободен хотя бы один из п приборов обслуживания. Вероятности этих гипотез следующие: Р(Н\) = рп\ Р(Н2) = \ — р^ Время обслуживания любого требования распределено по показа­ тельному закону с параметром р, поэтому Р(ТС<. t/H2) — 1 — е~ 11 при t > 0. Если все приборы обслуживания заняты, то очередное требование в систему не попадает, поэтому

Р(Тс < Ц Н г )^ в ^ )

0

при

£ < 0 ;

1

при

t > 0.

 

С учетом этих равенств соотношение (14.19) переписывается в’ виде

Fc V) = pnB{t) + (1 — />n)(l *-■“ )•

(14-20)

Производная от единичной функции e(t) равна дельта-функции б(^). Поэтому плотность распределения случайного времени Тс на­ хождения требования в системе

fc(t) = F'0( t ) = p ^ ( t ) + (\—pn)\ie-^ при t > 0. (14.21)

Математическое ожидание этой случайной величины

*с = | tfe( t)d t=

- ~ f - ,

(14.22)

о

 

 

Рассмотрим частный случай, когда в системе с отказами имеется только один прибор обслуживания, т. е. п = 1 . При этом уравнения (14.1) для вероятностей Po(t) и P\{t) записываются в виде

Я0(О =

- Х Р о (О +

рР1(^);

|

Р;(*) =

- рЛ(*) +

ЬРо(*).

I

103


Решение этой системы ищем в виде P0(t) = A0ert, Px(t) = А ^ .

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

в уравнения (14.23) и со­

кращения на еп приходим к равенствам:

 

- | - Х ) Л 0 — M i = 0 ;

(14.24)

 

 

 

 

— ХЛ0-f- (г + р) Л] = 0.

Решение этой системы

 

X

 

A l = l ± ± A

 

(14.25)

"

г +

Р

Р

05

т. е.

 

 

 

 

г +

Х

'

X

(14.26)

В, (г)

 

г

+

Р

 

Р

Характеристический определитель

 

 

 

r i~ X

— р

=

0 ,

 

 

X г +

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

 

 

 

г2-f- (X + р)г = 0.

 

 

 

Корни этого уравнения: г0 =

0; Г\ =

— (Я +

р ).

 

 

 

Согласно (14.6)

решение системы (14.23) записывается в виде:

 

P0(t) =

A00+ A 0le - P W ;

 

 

 

 

 

Pi ^ ) = - ^ - А ю - А пе - ^ .

 

 

(14.27)

 

 

 

 

 

 

 

Р

 

 

 

 

 

 

Если Исходным является состояние Со, то начальные условия

следующие: Р0( 0 )= 1 ;

Р](0) = 0. Тогда в соответствии с

(14.27)

 

^oo + ^oi — 1;

—~

-^оо

^oi — 0.

 

(14.28)

 

 

 

 

Р

 

 

 

 

 

Решение этих уравнений:

 

 

 

 

 

 

л ^

X

-

___L _ .

Л =

х

=

а

(14.29)

00

р

1 +

а ’

01

X —}—[i

 

1 — а

 

Следовательно, искомые вероятности:

p o ( ^ ) = - r b r l l + a e " (X+,l>tJ;

(14.30)

104


К такому же результату можно прийти другим способом, исклю­

чив, например, из первого уравнения системы (14.23)

функцию

P i { t ) = l — P0(t). Тогда для P0(t) получается уравнение

 

W =

04.31)

решение которого при начальном условии А>(0 ) = 1 определяется первой формулой из (14.30). Зная Ро{Ц, вероятность P\{t) нахо­ дим с помощью равенства Pi(t)— 1 — P&{t).

Если исходным является состояние Сi, т. е. момент ^ = 0 совпа­ дает с поступлением в систему требования, то вместо (14.28) будет

"4оо +

"4о1 — 0;

г ^оо

^oi ~ 1 *

 

(14.32)

При этом

 

 

 

 

 

 

Аоо —

1*

1

г 01

 

1

(14.33)

X -|- |1

1 + а ’

1 +

а ’

а искомые вероятности определяются формулами:

 

 

 

/>o(0 =--T ^ T

[ l - e - (X+rtt] ;

 

 

(14.34)

 

 

 

 

 

 

Функция Po(t) в любой момент времени t совпадает с вероят­ ностью простоя одноканальной системы массового обслуживания, а функция P\{t)— с вероятностью ее функционирования. Матема­ тическое ожидание y(t) числа Y(t) приборов, занятых обслужива нием в момент t, и коэффициент загрузки приборов &заг (0 в этом случае совпадают с P\{t). Если проинтегрировать функцию Р\(t) по t от 0 до Т, то получим математическое ожидание числа прибо­ ров, занятых обслуживанием в промежутке времени от 0 до Т. Ма­ тематическое ожидание времени обслуживания одного требования

равно

?!).= — . Поэтому

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

w(T) числа

W(Т)

I1

 

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

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

формулой

 

 

 

__

т

(14.35)

 

w ( T ) — р j P^{t)dt.

105