Файл: Шнепс, М. А. Численные методы теории телетрафика.pdf

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

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

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

Добавлен: 19.10.2024

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

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

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
Диаграмма переходов этого процесса изображена на рис. 1.8. Для процесса размножения и гибели каждое подпространство L>, содержит одно единственное со­
стояние.
При помощи чисел аху для х ф у задаются вероятности пере­ хода; вероятность перехода из х в у за промежуток времени At
равна axyAt+o^At) ери At-^О. На рис. 1.12 указаны интенсивности переходов аху, х ф у , для НС, (при­ веденной на рис. 1.2 и 1.9.
Для дальнейшего полезно иметь другую интерпретацию ин­ тенсивностей перехода, раскры­ вающую простой вероятностный смысл траектории процесса x(t). Как известно из теории процес­ сов Маркова, почти все траекто­ рии однородного марковского
Рис. 1.12. Диаграмма Хассе трехли­ процесса с непрерывным време­
нейной НС, дополненная интенсивно­ нем и дискретным множеством стями переходов
состояний имеют следующие свой-
22

3. Матрица интенсивностей перехода марковского процесса

Если для НС даны диаграмма Хассе и матрица занятия Я, а

также

дано, что

потоки вызовов

пуассоновские

(с интенсивно­

стью

X) и время

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

экспоненциально

распределенное

(с параметром, равным 1), то можем построить

марковский про­

цесс {x(t), t^O }

с конечным множеством состояний 5 и непрерыв­

ным временем t, описывающий деятельность данной коммутацион­ ной системы.

Матрица интенсивностей перехода А = {a,j,} марковского про­

цесса x(t) имеет элементы:

 

1.

у е в ,;

 

ьгху>

уе 4 ;

(37)

— I X ] — X S (х), у = х\

 

0для остальных случаев.

Вслучае процесса размножения и гибели, заданного на множе­ стве чисел 0, 1, 2, 3... (конечном или счетном), матрица А имеет элементы:

г-Н

ai, г-i

Рч) ^Zq,—1 0

 

(38)

a it[ — — (^г + Рг)

ai,i = 0

i — i I > 1


ства. Если в некоторый момент t процесс x(t) находится в состоя­ нии х, то время пребывания tx в этом состоянии является случай­ ной величиной, распределенной по экспоненциальному закону с параметром —ахх, т. е.

Р { ^ > 0 = е а* Л

(39)

н в момент t + tx процесс переходит в состояние у с вероятностью

х ф у .

(40)

--- вхх

Так как

&ХХ 2.J ахУ<

U&S, хфу

то эти соотношения определяют цепь Маркова, вложенную в про­ цесс Маркова над множеством моментов переходов.

4. Уравнения равновесия и стационарные вероятности

При помощи матрицы интенсивностей перехода А задается стационарный марковский процесс x(t), 0, принимающий значе­ ния из множества состояний S. Матрица переходных вероятностей

P(t) = {Pxv(t)}, где pxy(< t)= P {x(t)= y/x(0 )= х }, процесса x(t) удов­ летворяет уравнениям Колмогорова—Чепмена:

4at р (*)= ATP(t) = P(t)A

с начальным условием Р(0) =/ (I — единичная матрица). Решени­ ем этой системы является P ( t ) = e At.

Если множество 5 образует один существенный класс сбстояний без циклических подклассов, что имеет место в случае НС и ком­ мутационных систем вообще, то существует распределение стацио­ нарных вероятностей р = { р х}- Вектор финальных вероятностей яв­ ляется решением однородной системы

Атр = 0,

(41)

где Т обозначает транспонирование матрицы.

Для окончательного определения стационарного распределения требуется добавить лишь нормирующее условие 2 p i = l . С учетом

конкретного вида коэффициентов матрицы

I

(37) имеем

согласно

систему

 

 

 

 

[|х| + М *)]р х = £ Р у +

£

Pybryx,

x £ S .

(42)

у s ах

у е

в х

 

 

Это так называемые уравнения равновесия, имеющие простой физический смысл. Левая часть уравнения выражает среднюю ско­

23


рость выхода из состояния х, а правая часть — среднюю скорость входа в состояние х при условии, что коммутационная система на­ ходится в равновесии (в стационарном режиме).

5. Численный пример

Составим систему уравнений равновесия (42) для схемы, представленной на рис. 1.2, пользуясь диавраадмой переходов, изо­ браженной на рис. 1.12. Заметим, что из-за симметрии схемы мож­ но попарно объединить состояния, в которых занята одна из инди­ видуальных линий. При этом получим систему из 6 состояний вме­ сто 8 (23), что изображено на рис. 1.12. Введем обозначения со­ стояний (в фигурных скобках указаны номера занятых линий):

*1 =

{ }.

*2 = {1 или 2},

х3 = {3},

х4 = { 1, 2},

а-5 =

{(1,

3) или (2, 3)},

х6 = { 1, 2,

3}.

Тогда система (42) принимает вид:

2Хрх — р2+ р3

(2Х + \) ръ—2'крх -f- 2р4+ ръ

(2Х + 13 = р-а

(2Я, -{- 2) р\ = Xр2 Ре

(Я,-)- 2)ръ—Хр2+ р3 + 2р6

Зрв = 2X рх X р5

Одно из ее решений следующее (без учета нормирующего усло­ вия) :

Рх =

3 -(- 5Я,2;

р2 =

6 Х +

11Я,2 + 6Я.2;

р3 =

ЗЯ,2 +

4Я,3

р4 = ЗХ2 + 4Х3+ 2Х“;

р & = ЗЯ,2 + ЮЯ3 + 8Я4;*

р6 = ЗЯ,3 + 6Я4 + 4ЯЛ

Отсюда получаем вероятность потерь с учетом того, что по (35) У5= 1/2, у6= 1 и остальные уг= 0:

4 Р о + Ре

4 Я2+ 8Я3+ 10Я4+ 4Х6

 

л = ---------------

=

------------------------------------------------ .

(44)

Pi - r - ..

+Ре

3 + 13Я + 25Х,2 + 27Я3 + 16А.1 + 4ЯБ

 

24


Замечания и литературные ссылки

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

Описанию произвольных систем коммутации большое внимание уделяется в работах Бенеша [24], Башарина [12, 18]. Особенностью работ Беиеша является предположение о том, что матрица интен­ сивностей перехода марковского процесса, описывающего коммута­ ционную систему, является симметризуемой (об этом подробнее в § 11.3). Это предположение означает, что занятие любого свобод­ ного пути (пары контактов «вход—выход» и соединительного пути между ними) является равновероятным, что, конечно, на практике далеко не так. Башарин [11, 18] провел более глубокое теоретико­ множественное изучение свойств коммутационных систем, чем это делает Бенеш. Для построения уравнений равновесия Бенеш ис­ пользует две системы множеств (для каждого х 6 S строит множе­ ство Ах и множество В х), как это сделано в настоящей главе. Ба­ шарин для каждого состояния х строит четыре множества состоя­ ний: Г х, Г ~ х, Нх, Н ~1, которые дают более детальный анализ

структуры коммутационной системы (включая ее алгоритм управ­ ления). Заметим, что Н~1 совпадает с Ах и T j 1 совпадает с Вх.

Этих двух множеств достаточно для составления уравнений и вы­ числения вероятностей состояний системы, как это сделано в § 1.3.

Подробнее об условиях существования стационарного распреде­ ления процессов размножения и гибели (33) можно просесть в различных источниках, например [36].

Терминология теории графов (диаграмма Хассе и т. п.,) исполь­ зуемая при анализе состояний коммутационной системы, подробно изложена в [25, 114].

Г л а в а

2

Основные формулы теории телетрафика

Внастоящей главе на основе теории процес­сов размножения м гибели получены про­

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

2.1.ПРОСТЕЙШИЕ ФОРМУЛЫ

1.Формула Эрланга для полнодоступного пучка с потерями

Пусть дан о-линейный полнодоступный пучок с потерями при пуассоновском потоке вызовов интенсивности К и экспоненциально распределенной длительности разговора со средним значением, равным единице. Тогда вместо (1.38) имеем:

а,-,ж

=

^,

0

 

ас<

=

i,

0 ^

v,

и а,-j равны нулю при остальных значениях индексов, г ф j. Формулой Эрланга E V(X) называют вероятность занятия всех

линий, которая согласно (1.34) имеет вид

 

Е0(Х) =

Ы

( 1)

 

Это так называемая первая формула Эрланга. Введем еще распределение Эрланга

 

А/

 

 

E i,v(X)

П

 

* = 0, 1,

.v.

 

А*

 

1

 

(2)

 

k\

 

 

 

k=0

 

 

26


2.Формула Энгсета

Вотличие от предыдущего случая — вывода формулы Эрлан­ га, изменим только одно условие: вместо пуассоновского потока, что в математическом смысле возможно при бесконечном числе абонентов, 1иред(положи1М, что ноток создается конечным числом абонентов N. В состоянии с i занятыми линиями интенсивность

перехода вверх (в состояние £+1)

a iti+i = (Ni)%, £ = 0,

1,..., у . Сле­

довательно, имеем

процесс размножения

и гибели с

интенсив­

ностями:

0 < £ < у ;

 

 

 

Х[ = {N — 1)1,

 

 

 

р,- = £,

0 < £ < т ,

 

 

(3)

и равными нулю при остальных

значениях £. Подставляя (3) в

(1.34), имеем

 

 

 

 

Л =

+

V Ро =

Г ) Vpo,

(4>

где = С1'N — обозначение биномиального коэффициента.

С учетом нормирующего условия из (4) получаем формулу Энгсета — вероятность занятия у линий:

1 = 0

Введем два определения нагрузок, связанные с ф-лами (1) и (5): 1) пуассоновская нагрузка первого рода (та, при которой вы­

ведена формула Эрланга (1)); ‘ 2) пуассоновская нагрузка второго рода (нагрузка, создавае­

мая конечным числом источников нагрузки, как в случае вывода формулы Энгсета (5)).

Если при пуассоновской нагрузке первого или второго рода ин­ тенсивность обслуживания равна р, то в ф-лах (1) и (5) в качест­ ве X следует подставлять Х/р.

3.Формула Эрланга для полнодоступного пучка с ожиданием

Пусть поступает пуассоновская нагрузка первого рода на у-линейный полнодоступный пучок. Вызовы, поступившие в состоя­ нии с у занятыми линиями, не теряются, а ожидают. Соответствую­ щий процесс размножения и гибели имеет интенсивности переходов:

%i = К

 

i = 0, 1, 2,

 

P i =

i,

0 <

£ < у;

 

v,

£ >

у .

( 6 >

 

 

 

27