Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Шнепс, М. А. Численные методы теории телетрафика

.pdf
Скачиваний:
8
Добавлен:
21.10.2023
Размер:
9.21 Mб
Скачать

1) взяты всевозможные наборы значений мощностей контакт­ ных множеств Ц, i = l , 2,..., v, при соблюдении условия i2 li = nd\

2) для каждого набора значений построена схема, удовлет­ воряющая принципам 2—4.

Таблица 5.3 иллюстрирует сложность выбора оптимальной схе­ мы в области небольших нагрузок. При Л = 0,0625 (и Л-Ю, как до­ казано ниже) оптимальной является схема ж, содержащая мак­ симальное число индивидуальных линий, однако при несколько большей нагрузке Л =0,25 не полностью равномерные схемы г и е являются лучше, чем равномерная схема а и ступенчатая схема ж.

3.Колебание нагрузки увеличивает область оптимальности равномерных схем

Подтверждение этой гипотезы покажем на примере схем, представленных на рис. 5.1. Предположим, что параметр К пред­ ставляет собой случайную величину, распределенную по нормаль­ ному закону с заданным средним значением и дисперсией, равной 0,01. Для вычисления вероятности потерь в этом случае берем дан­ ные табл. 5.1 и интегрируем их, умножая на весовую функцию плотности вероятностей нормального распределения. Получаем

табл. 5.4, 'из которой следует, что точ­ Т А Б Л И Ц А 5.4 ка пересечения для схем а и б 'сдви­ нулась 1в область 'меньших (потерь и кривые пересекаются приблизительно при Л=0,55. В целом ib области ма­ лых потерь при колебаниях нагрузки вероятность потерь растет, а в обла­ сти больших нагрузок падает (ср.

табл. 5.1 и 5.4).

 

 

 

4. О равномерных схемах при

20

0,7376

0,7317

повторных вызовах

 

12

0,6010

0,5878

Можно высказать гипотезу, что

8

0,4633

0,4424

соотношение

между равномерными и

4

0,2113

0,1897

ступенчатыми

схемами

сохраняется,

1

0,7288-10-2

0,6775-10- 2

если перейти

от

схем

с потерями к

■схемам с ожиданием (или (повторными

 

 

 

0,6

0,1469-10- 2

0,1417-10“ 2

вызовами.

Численный

пример

содер­

0,5

0,8376-10- 3

0,9007-10- 3

жит табл.

5.5. В

ней приведены ре­

зультаты (вычисления вероятностей по­

0,25 0,1141-10“ 3

0,1408-10—3

терь при наличии (повторных (вызовов

Ш

 

 

в схемах, изображенных на рис. 5-1.

 

 

 

Предполагается,

что

каждая

груп­

па абонентов имеет по одному месту ожидания, которое за­ нимает абонент, потерявший вызов. Ожидающий абонент произ­ водит повторные вызовы (с интенсивностью р,) с целью занятия

90

свободной линии. Подробнее о повторных вызовах рассказано в гл. 6. Таблица 5.5 показывает, что наличие повторных вызовов су­

щественно

не влияет на вероятность потерь и точка

пересечения

Т А Б Л И Ц А

5.5

 

 

 

 

 

и

X

 

 

а

 

6

 

 

Л,

яа

я,

|

ТСд

 

 

 

 

0,375

0,2699-10-1

0,3753-10- 2

0,2314-10—1

 

0,3133-10- 2

1

0,1875

0,2913-10- 2

0,1889-10—3

0,2864-10- 2

 

0,1846-Ю- 3

 

0,09375

0,2761-10- 3

0,8729-10- 5

0,3296-10- 3

 

0,1043-10—4

 

0,375

 

0,2769-10-1

0,3650-10“ 2

0,2355-10-1

 

0,3013-Ю- 2

2

0,1875

0,2956-10~2

0,1876-10—3

0,2892-10~2

 

0,1819-Ю- 3 '

 

0,09375

0,2780-10—3

0,8709-10- 5

0,3313-10—3

 

0,1037•10 4

кривых вероятности потерь я2 (отношение потерянных вызовов ко всем вызовам) практически остается на месте (ср. с табл. 5.1). То же можно сказать относительно щ — вероятности потери первич­ ного вызова.

5. Вероятность потерь при регулярном потоке

Замена пуассоновского потока регулярным потоком (длитель­ ность интервала фиксирована и равна единице) несколько умень­ шает область оптимальности равномерных схем. Таблица 5.6 пока-

Т А Б Л И Ц А

5.6

 

 

 

 

 

Л

а

б

Л

а

 

б

20

0,7334

0,7279

2

0 , 1246-Ю-1

0,1114-10-1

12

0,5875

0,5741

1

0,1419-10

3

0,1903-10—3

8

0,4354

0,4132

0,5

0,1315-10

6

0,2461-10- 6 '

4

0,1509

0,1303

 

 

 

 

зывает, что точка пересечения кривых вероятности потери для схем на рис. 5.1 переходит из точки Л =0,7, где они пересекались при пуассоновском потоке (см. табл. 5.1), в более высокую точку око­ ло Л = 2 при регулярном потоке. Заметим также, что при регуляр­

ном потоке, как и следовало ожидать, вероятность потерь

падает

особенно сильно при малых Л, например, в точке Л =0,5

она в

5000 раз меньше, чем при пуассоновском потоке.

 

91.

6. Схемы, для которых кривые вероятности потерь пересекаются дважды

Численные данные о схемах с небольшим числом линий (на БЭСМ-ЗМ можно вычислить вероятность потерь для схем с и ^ Ю ), а также приводимые ниже степенные разложения вероятности по­ терь развивают нашу интуицию в том направлении, что при упоря­ доченном искании свободной линии в области малых потерь опти­ мальными являются ступенчатые схемы, а в области средних и больших потерь — равномерные схемы. И хотелось бы думать, что кривые вероятности потерь любых двух схем пересекаются только

Т А Б Л И Ц А 5.7

О О

О—О

X

—О

О

о

о

 

У

2

0,3089

 

0,5

0,721 МО- 2

0,25

0,3770

-10~3

0,015625

0,2762

-10- 8

 

К

П

Ш

 

В)

0,3052

0,4178

0,7311 •10_ 2

0,2184-10-1

0,3951 -10—3

0,1468-10- 2

0,2750-10~8

0,8680-10~8

0—0 О—О

0,4151

0,2263-10-1

0,1677-10—2

0,7824-10- 8

однажды. Однако это не так, что сильно осложняет исследование неполнодоступных схем в области средних потерь. Результаты табл. 5.7 показывают, что для схем а и б, а также для схем б и г, соответственно, кривые вероятности потерь пересекаются дважды.

7.При случайном поиске оптимальными являются равномерные схемы

Данное утверждение хорошо иллюстрирует табл. 5.8, где приведены те же схемы, что в табл. 5.3. Там было показано, что при упорядоченном искании свободной линии преимущество рав­ номерных схем перед ступенчатыми по мере нарушения равномер­ ности схемы уменьшается и сохраняется только в области больших нагрузок, В табл. 5.8 можно четко проследить, что нарушение рав­ номерности в построении схемы равномерно увеличивает вероят­ ность потерь во всем диапазоне нагрузок, что равномерно лучшими являются полностью равномерные схемы (схема а в табл. 5.8) и самыми худшими являются ступенчатые (схема ж ). Преимущество схемы а перед схемой ж подтверждают и асимптотические разло­ жения 4) вероятности потерь по степеням X:

P a W =

^ 3+ . . . ;

/>*(*) =

+

I I

 

 

00

') Методика

получения асимптотических

разложений изложена в §§ 5.2

и 5.3.

92

откуда следует, что ори доста­ точно малых X схема а лучше схемы ж, а также асимптоти­ ческие 'разложения то степе-

О О О О

которые ‘подтверждают, что схема а лучше схемы ж и в области больших нагрузок. На­ ряду е численными данными о вероятности потерь inip-и конеч­ ных значениях X, ‘приведенны­ ми IB табл. 5.8, 'можно .вьюка- «ать нипотезу, что «три случай­ ном .иокании равномерные схе­ мы являются оптимальными во всем диапазоне нагрузок.

8.Изменение числа групп осложняет выбор оп­

тимальной схемы

 

В ьгше

расоматривалнсь

 

НС при

фиксированных

п, d,

 

v в 'зависимости от X. Однако

 

•для практики ‘большое

значе­

 

ние имеет также выбор оп­

 

тимального параметра п — чис­

 

ла групп абонентов, так как он

 

влияет на объем коммутацион­

 

ного

оборудования.

Молено

 

•высказать предположение, что

 

увеличение п приводит к умень­

 

шению потерь. В табл.

5.9

при­

со

ведены

результаты

расчетов

вероятностей

потерь

для

трех

<

НС с параметрами d = 2 ,

v =8,

Я

а п

меняется

от 2

до

4.

Так

К

как

Л

одно

in

тоже,

то

*=;

<

Х=А/п.

 

 

 

 

 

 

 

 

 

 

 

 

н

О—О 0—0 ^

О0—0 о

с•{ о—сГЪ

6

Q 6 р V

о—о о—с

со

о о о

—м

со т со со см

юСО 05

Tf Г- ■

О о о о

1

 

со

 

о

о о

со 1Л

00

to

СО

05

тр см

 

СО

О * о

о

о"

 

 

СО

о

о о

—1

■—1

00

т

о

СО

со

со

г— со

о

СО

 

СО

Оо" о о”

 

 

 

со

 

о

о

!

 

о

 

t-» о

—1

 

см

 

 

СО

оО

 

СО

см

 

 

см

о

о

о

о

 

 

 

со

 

о

о о

 

—■ — —*

 

СО

05

СО

 

со

со

 

СО

см

 

о

о

о

о

со

оо о

—1

CN СО

05 о с

00 о см

ю см

оо о о

сч со

о

о

о

—•

 

СО

о

см

СМ

о 00

а>

со

СО

 

00

*-*

о о о

о

ю

см ю со

Юсм о

см о о о

93

Таблица 5.9 показывает, что три малых нагрузках оптимальной является ступенчатая схема а, содержащая индивидуальные ли­ нии, а при больших нагрузках — равномерная схема б. Это соотно-

T А Б Л И Ц А 5.9

А

п

-

 

 

 

 

 

 

а)

 

У

 

40

0,9280

 

0,9278

0,9318

4

0,4911

 

0,4862

0,5329

1

0,1049

 

0,1064

0,1107

0,4

0,2065-1 0 - 1

0,2271-Ю -1

0,2234-10_|

0,04

0,2035-10- 3

0,2630-10_3

0,2481-10_3

0,004

0,2004-10- 5

0,2663-10- 5

0,2498-10-5

шение между схемами а и б сохраняется в области Я-»-0 и Л->-оо. Первые члены соответствующих асимптотических разложений мож­ но получить из аналитических выражений для вероятности потерь:

=

X2 (8Ь3 + 20Х2 +

 

16*,+ 3) .

 

Па

2(1 + Ь)3 (4Х2 +

4Ь + 3)

 

_

ЗА2 (1 +2Х )

 

 

.

 

 

 

Лб

6Л« + 9А.2 + 6А. + 2 ’

 

 

 

=

X2 (192Х3

+

256Х2

+

112?, +

15)

^ _

192Х,6 + 400Х4

 

+

256Х3

+

1Ш 2 +

49Х + 6

1

3

А,2-г... Подставляя

Например, при Л,-й) яа= — А2+ ...;

я б = —

К = А (2 в первом и Х=Л/3 во втором

случае,

имеем па — 8 Ла +

_3_ Л2+ ... .Отсюда следует, что схема а лучше схемы б три

18

малых значениях Л.

Однако высказанное в заглавии настоящего пункта предполо­ жение выполняется не всегда, что, ви­ димо, связано с дискретным характером НС: при увеличении числа групп п на 'единицу тип НС может измениться ка­ чественно. Такого рода противоречащий высказанному 'предположению пример содержит рис. 5.2, где приведены две

Рис. 5.2. Схемы, иллюстрирующие сложность вы­ бора оптимального числа групп п

94

схемы

с

параметрами d = 2, и = 5, но отличающихся по п:

а) п= 6;

б)

« = 4. Исходя из стремления заменить первую вертикаль

схемы на индивидуальные линии при малых значениях Л,, можно предположить, что схема б лучше схемы а. Однако, как показыва­ ют асимптотические разложения при зуалых А (к= А /п ), схема а лучше схемы б:

=-j- А2 + о (А2); л:7= -L- А2 + о(А2).

5.2.СТЕПЕННЫЕ РАЗЛОЖЕНИЯ ВЕРОЯТНОСТЕЙ СОСТОЯНИЙ НЕПОЛНОДОСТУПНОЙ СХЕМЫ

1.Уравнения для стационарных вероятностей состояний

При предположении, что на НС с параметрами п, d, v посту­ пают п взаимно независимых пуассоновских потоков (каждый с па­ раметром л), а время обслуживания подчинено экспоненциально­ му закону (с параметром 1), действие НС можно описать марков­ ским процессом с конечным множеством состояний 5 (число их равно 2и).

Как уже введено в § 1.3, в состоянии * 6 S через |х| обозначим число занятых линий. Обозначим через Ах множество тех состоя­ ний, из которых освобождением одной из |х| + 1 занятой линии можно перейти в х; Вх — множество тех состояний, которые полу­ чаются освобождением одной из |х| занятых линий в состоянии х. Пусть гху — число тех групп абонентов, от которых вызов, посту­ пивший в состоянии х, переводит схему в состояние у; s(x) — чис­ ло групп, которым доступна хотя бы одна свободная линия в со­

стоянии х, s(x ) = 2 гху.

иеАх

Стационарные вероятности рх являются решением системы алге­ браических уравнений, имеющей вид

[|*| + Я ,«(*)]/?,= £

Ру+Ъ Yi Ругух, * 6 5 .

(1)

УеАх

уев х

 

Система (1) является линейной однородной системой уравнений по отношению к {р^ х € 5 } , причем определитель системы равен ну­ лю (одно уравнение является следствием других), и, следовательно, система (1) имеет нетривиальное решение вида

Px = PxW,

(2)

где Я* (Я) — многочлены по X с целыми коэффициентами, и любое другое решение (1) пропорционально этому. При использовании нормирующего условия

в е р о я т н о с т и р х о п р е д е л я ю т с я о д н о з н а ч н о :

 

Рх(М

(4)

Рх

X Р*М

л-g S

т. е. рх получаются в .виде рациональных функций от А. Зная рх, нетрудно вычислить вероятность потерь

ле s I л |><г

или

X У(х)рх,

 

* =

6

 

X£ S

()

где у(х)

— условная вероятность потерь:

 

Hv) = i _ jlW .

(7)

2. Решение в виде разложений по степеням А,

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

Переходим к выводу разложений по степеням л. Ниже через гх обозначен элемент с индексами 0, х в |х|-й степени матрицы R, где R — матрица, составленная из элементов гху (остальные эле­ менты равны 0), т. е. первая строка в матрице Юх' . Легко видеть, что гх равно числу различных строго возрастающих путей (т. е. со­ стоящих из одних занятий), которые переводят схему из состояния Ов состояние х.

Будем искать решение системы (1)

в виде

 

 

рх -

(с„ (х) + сх (х) А +

с2 (х) А2 +

. . .),

х 6 S.

(8)

Подставляем выражение (8)

для рх в систему (1):

 

[ |х| +

A s(x)]

{?J ^(с„ (х) -|- Ci (х) A -J- Са (х) А2 +

• • •)} =

 

=

I

tix l + l {c0(y) + c1{ y ) X + c 2(y)W +

. . ,)}2+

 

 

УеАх

 

 

 

 

 

+

X гух^ Х] (со (У) + ci (У) АД- • • •)•

 

 

 

уевх

 

 

 

 

 

96

Приравняем нулю коэффициенты при одинаковых степенях X. При А,1-'1 имеем

|а'|с0( * )=

Y. гухс0(у),

 

 

 

 

(9)

приА,Л'|+ 1

i/e в х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\х\С1 (х) +

.?(х) с0(х) = 2

Ф ) +

2

V

i (У)

 

и т. д. В общем случае

 

уе лх

уе вх

 

 

 

 

 

 

 

 

 

|Л’ I ст (х) + 9 (х) сш_ 1 (х)

=

V

сш_,(г/)+

2

rvxcm{y),

т = 1,2, . . .

 

 

 

УеАх

 

 

Уевх

 

(19)

Так как из системы (1)

рх определяются с точностью до нормирую­

щего множителя, то положим:

 

 

 

 

 

со (0) = 1,

ст(0) == 0,

щ ф 0,

 

 

 

(11)

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

рекуррентной

системы (L0).

 

 

 

 

 

 

 

 

Выводим второе граничное условие. Для хб L r (т. е. |х] = 1) из

(9) получаем

 

 

 

 

 

 

 

 

=

гУхс0 (У) =

г0хс0(0) -

г0х\

 

(12)

 

У & вх с . ц

 

 

 

 

 

 

 

длкх(ЕТ2 из (9)

и (11)

следует

 

 

 

 

 

со (*) = -£-

2 ] V »

=

2 ]

Г°'/ух =

2\Гх

 

 

yeBxcLt

 

 

.ves^ct,

 

 

 

и т. д. В общем случае получаем второе граничное условие для ре­ куррентной системы (10):

СоМ = г £г;> x€ s -

(13)

I* И

 

Таким образом, получаем, что на основе рекуррентного соотно­

шения (10), переписанного в виде

 

CraW =

W cm-1W +Cm2]- 1 (У ) + 2j Ф

т Ц ’

:(14)

 

yeAx

у&вх

 

 

совместно с граничными условиями (11)

и (13) можем найти раз­

ложения вероятностей состояний по степеням X в виде

(8).

 

Граничные условия

(11) можно видоизменить так,

чтобы найти

степенные разложения рх, удовлетворяющие нормирующему усло­ вию, а именно, следует брать

с0 (0) =

1;

 

Ci (0) +

2

Со (х) — 0;

 

-veL,

(15)

са (0) +

2

Ci (х) + 2 с0 W = 0.

 

хев,

xql,

А— 264

9 7

Соотношения (13)— (15) также определяют решение системы (1).

П р и м ер . Рассмотрим трехлинейную НС, изображенную на рис. 5.3 и подробно исследованную в § 1.3.5. Вычислим для нее пер­

/

 

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

разложения mo X,

удов-

Д

летворяющих

условию

 

нормированности

(15).

О

 

у

Для анализа схемы .введем обозначения состоя -

2

I

ний x = ( i j k ) ,

где i, /,

k

относятся

к 1, 2

и 3-й

О

Л.ШНЧ1Я1М соответственто

и

равны 0,

если

линия

О

АВ .свободна, и 1, если линия занята. Матрица и«-

Р и с .

5 .3 .

Трех-

течтсивностей

терехода,

как

ои а постр

$ 1.3.5, 'имеет .вид

 

 

 

линейная НС

 

 

 

 

— 2 Х

X

X

0

0

0

0

0

/ 1

— 1— 2 Х

0

0

X

X

0

0

 

1

0

— 1 — 2Ь

0

X

0

X

0

А ==

1

0

0

-1 — 2 Х 0

X

X

0

0

1

1

0

— 2— 2 Х

0

0

2 Х

 

 

0

1

0

1

0

- 2 — X

0

X

 

0

0

1

1

0

0

— 2 - -X

X

'

0

0

0

0

1

1

1

— 3

Перед вычислениями отметим,

что

 

 

 

0

1

1

0

0

0

0

0

0

0

0

1

1

0

0

0

0

0

1

0

1

0

0

0

0

0

1

1

о

о

0

0

0

0

0

0

0

0

0

2

 

0

0

0

0

0

0

0

1

0

0

0

0

0

0

0

I

1

Первая строка матрицы Д2= (0 0 0 0 2 1 1 0), первая строка матрицы R3= (0 0 0 0 0 0 0 6).

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

ной записи матрицы А:

1) составляем таблицу коэффициентов, взяв матрицу А и опус­ кая в ней К;

2)справа от нее пишем компоненты векторов-столбцов Со, сь сг,

3)вычисление коэффициентов начинаем с того, что пишем

с0(000) = 1

согласно (15);

условия

в (13)

имеем Со(100) = 1,

4) из второго граничного

fn(€10) =1,

Со(001) =0, с0(110) =

I

=1

(г1ю=2,

потому

|*|!|*=(1Ю )

что имеем две возможности непосредственного перехода из (000) в (110)), с0(101)=Со(011) =il/2, со(111) =6/31 = 1;

эк

5) из BTioipoiro условия iB (15) ct (000) = —[со( 100) + cO'(010) +

+Co(001)] = —2;

6)из (14) находим Ci(100) = —2-1 +1 + 1/2+1 •(—2) = —5/2итак

весь столбец сп

7)

из третьего уравнения в (15)

 

с%(000) = — 5/2 — 5/2 + 1 + 1 + 1/2 + 1/2 = — 2;

8)

из (14)

 

св(100) = — 2(— 5/2) + (— 3 — 1/2)+ 1-2 = 7/2

и так весь столбец с% Результаты вычислений запишем в виде таблицы:

 

Со

Cl

Со

(000)

1

—2

2

(100)

1

—5/2

7/2

(010)

1

—5/2

7/2

(001)

0

1

—3

(НО)

1

—3

16/3

(101)

1/2

-1/2

—2/3

(011)

1/2

-1/2

—-2/3

(111)

1

—7/3

28/9

Подставляя результаты вычислений в (8), получаем

р(000) =

1 — 2Я + 2Я2 +

* »

 

 

р(Ю0) =

Я — 5/2Я2 +

7/2Я3 +

• • •;

 

Р(010)=ь

Я — 5/2Я2 +

7/2Я3 +

• - •;

 

Р(001) =

Я2 — ЗЯ3 +

• - •;

 

Р (1Ю) =

Я2 — ЗЯ3 +16/ЗЯ4+

• ■ •;

р(101)

=

1 /2Я2 — 1 /2Я3

— 2/ЗЯ4 +

• ■

Р(011)

=

to 1

■1/2Я3

— 2/ЗЯ4+

• ■

То >>

р(111) =

 

Я3

— 7/ЗЯ4 +

28/9Я5 +

Так как вероятность потерь рассматриваемой схемы

* = ^-[р(Ю1)+р(011)] + р(111),

то получаем первые члены разложения л по степеням Я:

л = 1/2Я2 + 1/2Я3 — ЗЯ4 + • • ;

3.Решение в виде разложений по степеням Я-3

Выше было рассмотрено, решение системы (1) в виде разло^ жений по степеням' Я.' Такое решение дает представление о поведе­ нии НС при малых Я, т. е. при малых нагрузках. Для изучения НС

Соседние файлы в папке книги из ГПНТБ