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

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

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

4. Локальное укрупнение состояний

Иногда описание коммутационной системы можно упростить укрупнением состоянии системы. Изложим алгоритм укрупнения, основанный па применении метода исключения Гаусса.

Систему уравнений для расчета стационарных вероятностей со­ стояний

£ a }ipj = 0, / =■ 1. . . ., п,

(38)

i=i

 

с условием нормировки

 

V ру = 1

(39)

;'=1

можно преобразовать следующим образом. Из первого ур-ния (38) выражаем рх как функцию от остальных неизвестных

 

(40)

и подставляем это значение в (38), что дает

 

= 1.

 

/>2

(39), а именно:

Это соотношение получает вид, аналогичный

2 < i, = i.

н о

/>2

 

при введении новых неизвестных

 

 

 

 

; >

2 -

 

т

Неиспользованная часть ур-ний (38) после подстановок (40) и

(42) принимает вид

 

 

 

Yi

bjfli = °.

i = 2. • • -Л,

 

(43)

i>з

 

 

 

 

 

где

 

 

 

 

 

b

= а и а 1 + а и а и

^

.2>

(44)

 

а-х +

a i i

 

 

 

 

 

 

Очевидно, если а,ц 0, то

4i = Pj, bjc = ац. (45)

Полученная новая система ур-ний (43) и (41) имеет ту же струк­ туру, что и исходная система (38) и (39), только на одну неизвест­ ную меньше. При этом изменяются значения лишь тех неизвестных и могут изменяться значения лишь тех коэффициентов, которым в

80

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

Приведем соотношения для одновременного исключения состоя­ ний с номерами 1, 2,..., k, /г<п. Получаем систему вида (43) и (41), где j ^ k + 1, а формулы преобразования имеют вид:

Я 1 = -^ Р Г ,

(46)

Ьц = ^ ~ ,

(47)

ci

 

--- ^21 • .

— 0*1

д =

 

---a l k ---- °2* ‘

• •а к

a Cj н dji — соответствующим образом окаймленное Д:

1

1 .

. 1

 

an

alt .

■ -a ki

 

Cj = —Яд

аг . .

Oftl

, j

k , dji ~ а ц

<h ■

■—akl

 

ftjk <hk • • •ak

 

—ajk

aik

ak

 

При этом расчеты следует проводить только для тех j > k ,

для

которых имеется

хотя

бы

одно а ^ > 0 о i ^ /г. Для всех других /

имеет место

(45).

 

 

 

 

 

 

В случае больших k для расчета значений определителей можно

использовать их приведение к наддиагональному

виду. При этом

первые строки всех ёц

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

как

и первая строка Cj, и значения Ьц равны отношениям конечных пер­ вых элементов этих строк. Легко проверить, что коэффициенты Ьц имеют те же свойства, что и величины ад, только с учетом измене­

ния множества значений

индексов. Очевидно, с: > 0 и

0 при

в силу положительности и неотрицательности значений эле­

ментов первой строки соответственно.

 

Соотношению

 

 

— 6, / — 1, .

. ., П,

(48)

L

 

 

соответствует

 

 

2 &д = 0* i = k + \, . . ., п,

(49)

i>k

 

 

81

а последнее будет иметь место одновременно с таким же соотно­ шением для dji. Сумма же последних равна

S a /(

V flli .

*

V aki

i>k

i>k

i>k

а.ц

аг ■

 

-Oki

&jk

 

 

ak

В этом определителе в силу (48) элементы первой строки рав­ ны сумме всех других элементов соответствующего столбца. Следо­

вательно, имеет место (49) и Ьц = —БЬд. i'ti

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

5. Укрупнение посредством производящих функций

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

мерованными от 1 до п, причем при малых

значениях

номера t

интенсивности Аа = а,-,г--н1 и ц,г =

i_t имеют произвольные положи­

тельные значения и, понятно,

p,i= 0, а при i ^ k

p.j='[x (но

Х„ = 0). Такой процесс описывает, например,

^-линейную систему с

ожиданием.

При расчете стационарных вероятностей р, вводим производя­

щие функции

 

 

 

Р (*) =

pk +

pk+xx + . .

. + pnxn~k t

(50)

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

 

 

Х Pk+j

(Ц +

Ц Pk+j+l +

Г Pk+ л -2 ~ 0

 

соответственно на х*+2, уравнение

Рп-1 + И- Рп = 0

на xn~k+i и суммируем. В результате получаем

(X1)(A.JC — |Х) Р (*) + (* + 1 ) [|ЛPk — xn~“+l Х Рп\+Х\х Рк—т + 1 1= °-

82

Подставляя х = \ и х —ц/Х, а также полагая а=Я/ц, находим

Pk+l

~

а Pk

|

(51)

 

n—k

Г

Рп — а

Pk

I

 

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

 

 

Р{х)

=

1 — ах

 

 

 

 

 

откуда

 

 

 

п—*+1

Р(1)

 

: Pk +

• • •+

 

Рп== Pk - — ----

 

 

 

 

1 — а

Для определения ри Рь—, Pk имеем систему, состоящую из неис­ пользованных исходных уравнений, и условия нормировки

Pk

1 _

t/'-ft+i

 

(52)

 

1— о

Эту систему можем решить рассмотренным выше способом. Если

же нас интересуют только pi с i^ k ,

то требуется определить толь­

ко Ph, а дальнейшие

находим из

(51); отсюда же можем полу­

чить среднее число ожидающих в очереди и т. д.

Аналогичным образом, если требуется определить среднее вре­ мя достижения 7\- некоторого состояния Xj, j< .k , полагаем

Т (*) = Tk +

хТk+\ + . . .+ * " k Tk.

Действуя, как выше с уравнениями, где в коэффициентах фигу­

рируют только Яиц,

получаем

 

 

(х — 1) (цх — Х)Т{х) — \ixTk-i

+ h T k +

([Iх— Я) х

"1" Тп — х

х? -)- .

. . -f- хп.

Здесь подстановки х = а

и х = 1

дают соответственно:

Tk—Tk-1 +

1

1 _

g» -fe+i

 

 

 

1— с

 

(53)

 

 

 

 

 

1

 

 

Т „ = 7 * - , +

 

 

 

Р ( l — o n ~ k )

ц'Д — о) п k + 1

 

1— а

Для определения Ти, Th-i и Д, i</e— 1, ‘имеем 'неиспользованные исходные уравнения и (53). Затем получаем Т(х) и т. д.

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

а) использование части исходных уравнений при выводе соот­ ношений для производящей функции;

83

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

случае соотношения для величин />,•и Тп);

в)

решение систем уравнений типа (53), содержащих соотноше­

ния,

полученные в б), п исходные

уравнения, не использован­

ные в а );

 

г)

в случае надобности — расчет

величин, определяемых про­

изводя щимн фуикцнями.

Соотношения для производящих функций, естественно, позво­

ляют определить поведение процесса при п-*-оо. В

рассмотренном

примере как для того, чтобы р„->-0, так и для того,

чтобы Г/, оста­

валось

конечным, необходимым и достаточным условием являет­

ся а <

1.

 

Использование частичных производящих функций может быть полезным и для более сложных систем обслуживания. Если обра­ ботка уравнений для производящей функции связана с громоздки­ ми математическими выкладками, то последнее можно поручить ЭВМ при использовании соответствующих программ.

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

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

ссылки

В

настоящей

главе использованы статьи Гринберга [42,

43],

Гринберга

и Дам'бита

[44] и другие 1мате-риалы, любез­

но предоставленные Гринбергом. При описании марковского про­ цесса в виде потоков в графе автор следует монографии Бержа [25]. Комбинаторный метод классификации состояний, рассматриваемый в § 4.1.2, более простой, чем алгоритмы, ко­ торые обычно рекомендуются для анализа марковских процессов, например, действия со степенями матрицы А [137, 145] или других матриц того же порядка либо с собственными значениями и векто­ рами таких матриц [31]. Оба эти способа весьма трудоемки, если число состояний достигает нескольких сотен или больше. Для клас­ сификации состояний марковских цепей Сарымсаков [125] предла­ гает еще более трудоемкий прием, требующий как возведения матрицы типа А в степени, так и исследования всех траекторий, ис­ ходящих из отдельных состояний.

Теорема 4.1 обобщает правила, данные в свое время Кирхго­ фом [239] для расчета электрических цепей и широко реклами­ руемые в последнее время под различными другими названиями в качестве метода расчета электрических цепей [109]. Такой способ расчета предлагался также и для марковских систем [26].

84

Изложение итерационного алгоритма в § 4.3.2 следует статье Седола [129]. Этот алгоритм применяется в теории телетрафика давно (например, Карлсон и Эллдин [176]).

Среди алгоритмов решения системы (38) необходимо упомянуть решение в виде степенных рядов, т. е. поиск решения в виде беско­ нечных разложений по степеням А или А,-1:

Pi =

А * (с/о +

Сд А + С/гА2 - ) - , . . )

(54)

или

 

 

 

Pi =

~Т~ (d

+ йц-^- + dj2-^ -+ . . .j .

(55).

 

А.

 

 

Коэффициенты таких разложений вычисляются по

рекуррентным

соотношениям, определяемым системой (38) при .подстановке туда (54) или (55). Таким образом, можно получить несколько первых членов разложения, которые характеризуют систему при малых или больших нагрузках. В следующей главе на примере изучения прин­ ципов построения оптимальных неполнодоступных схем будет по­ казано, как работает этот метод.

Чтобы сравнить возможности различных алгоритмов, приведем некоторые данные о порядке систем линейных алгебраических уравнений, решаемых на ЭВМ, имеющей быстродействующую фер­ ритовую память на 64 00,0 слов. Методом Гаусса можно решать си­ стемы 200-го порядка, а при использовании памяти на магнитных дисках до 2500-го порядка. Если учесть, что матрица А является квазиякобиевой (ненулевые элементы находятся в полосе вокруг диагонали), то порядок решаемой системы можно увеличить. На­ пример, если полоса имеет ширину 600 членов (по 300 членов с каждой стороны диагонали), то можно решать системы 6000-го по­ рядка. Итерационным методом удается решить системы с 40 000 уравнениями. Заметим, что это не так уж много, если учесть, что неполнодоступная система с 15 линиями описывается системой

215 = 32 768-го порядка.

Г л а в а

5

Исследование неполнодоступных схем

Внастоящей главе подробно рассмотрена од­ на из классических проблем теории теле­ трафика — проблема выбора оптимальных

неполнодоступных схем (НС). В § 5.1 приводятся численные при­ меры, решенные на ЭВМ и иллюстрирующие различные принципы построения НС при конечных значениях нагрузки. Названия пунк­ тов в § 5.1 представляют собой перечисление этих принципов. Ос­ новной практический вывод — это целесообразность широкого при­ менения равномерных схем. Параграф 5.2 содержит математичес­ кий аппарат степенных разложений для оценки вероятности потерь при предельных больших и малых нагрузках. На основе этого ма­ тематического аппарата в § 5.3 доказаны важные предельные тео­ ремы: 1) об оптимальности равномерных схем при больших на­ грузках и 2) об оптимальности схем, содержащих максимальное число индивидуальных линий, при малых нагрузках. В § 5.4 изла­ гаются приемы построения равномерных НС.

5.1.ЧИСЛЕННЫЕ ПРИМЕРЫ АНАЛИЗА НЕПОЛНОДОСТУПНЫХ СХЕМ1

1.Существуют ли неполнодоступные схемы, равномерно лучшие по %

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

На рис. 5.1 приведены две НС, имеющие одни и те же парамет­ ры: число входных потоков п = 4, доступность d = 3, число линий и = 6. Приведенные НС отличаются только видом включения, т. е.

86

видом распределения контактов коммутационной схемы

(всего их

n d= 12) и по и = 6 линиям.

одной из-

Схема, представленная на рис. 5.1а, — это пример

ступенчатых схем, которые еще в 60-х годах рекомендовались для.

а)

Ю

1 3

5

 

X

X

 

 

 

X

X

 

 

Рис. 5.1. Две простые

X

Л

И

 

неполнодоступные

 

схемы:

 

Л

 

а) ступенчатая; б)

Л

сГ

сГ

равномерная

практического применения1). В табл. 5.1 приведены результаты вы­ числения вероятности потерь рассматриваемых схем, которые ста­ вят под сомнение целесообразность такой практической рекоменда-

Т А Б Л И Ц А

5.1

 

 

 

 

л

(Л) для схемы на

 

 

рис.

5.1 о

рис. 5.16

 

28

0,80657

0,80132

4

24

0,77800

0,77129

2

20

0,73956

0,73074

1

16

0,68517

0,67315

0 ,5

12

0,60264

0,58568

0,25

10

0,54335

0,52310

 

8

0,46459

0,44079

0,125

6

0,35720

0,33097

0,0625

я (Л)

для схемы н а

 

рис. 5.Ю

 

 

рис.

5.16

 

0,21187

 

0,18903

 

0 ,5 1 8 8 - 10~ 1

0,4489

-10

11

0 ,7 0 8 8 - 10~2

0,6614 -10 —г

0 ,7205 -10- 3

0 ,7 9 5 8 - 10-3,

0 ,6 9 1 2 - 10- 4

0,9253

- 10~4

0,6973 -10

5

0,1094•10

*

0,7572 -10—6

О

со to

о1

ЧУ

ции. А именно, оказывается, что ступенчатая схема на рис. 5.1а является оптимальной только в области малых потерь — приблизи­ тельно до я< 0,001 или, другими словами, до A fa 0,.7 (при этом А = 4Я). Для всех остальных значений параметров к более выгод­ но применять равномерную схему, приведенную на рис. 5.16. Тем самым на вопрос, поставленный в заголовке настоящего раздела, следует ответить отрицательно, так как при данном наборе пара­ метров п, d, v не существует схемы, равномерно лучшей при всех значениях интенсивности потока к (это доказано в § 5.3).

Заметим, что для вычисления вероятности потерь любой из двух, схем, изображенных на рис. 5.1, необходимо решить систему линей­ ных алгебраических уравнений порядка 26 = 64. Учитывая симмет­

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

87

рию схем и потоков, состояния удалось укрупнить и тем самым по­ низить порядок системы. Минимальное марковское разбиение ис­ ходного пространства состояний для схемы на рис. 5.1а имеет ^ с о ­ стоянии, для схемы на рис. 5.16 — 28. Методика составления систе­ мы уравнений подробно изложена в § 1.3.

2. Принципы построения равномерных схем и их численное подтверждение

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

1. Контактное поле с параметрами п, d, v разбивается на кон­ тактные множества, содержащие по г или по г+1 контактов, где

г — д г (здесь квадратная скобка — знак целой части). 1

2.Контактные множества образуются так, что лднии, доступ­ ные г группам, размещены левее линий, доступных (г+1)-й группе (направление искания слева направо).

3.Каждая пара линий имеет число общих групп (по возмож­ ности), одинаковое с числом -общих групп любой другой пары линии.

4.Каждое отдельное контактное множество размещено (по воз­ можности) в минимальном количестве вертикалей контактного поля.

Приводимые ниже примеры показывают, что принципа 1 целесо­ образно придерживаться при средних и больших нагрузках, в то время как нарушение любого из принципов 2—4 приводит к увели­ чению потерь при всех нагрузках.

В качестве численного подтверждения применимости равномер­ ных схем рассмотрим два набора схем с одними и теми же пара­ метрами /2= 4, d = 3, а = 6. Каждый набор включает две схемы, рас­ смотренные выше. Данные по первому набору сведены в табл. 5.2.

Т А Б Л И ЦА 5.2

1

“ П

Xd

in

о 0 о о

О О О

0 0 °

QО. О

ОО О

8)

0

? 0

0—

0—

0—0

<* О О—О

4

0,6731

0,6852

0,6891

0,6738

0.6732

1

0,1890

0,2119

0,2500

0,1930

0,1928

0,25

0,6614-10-2

0,7088-10-2

0,2041 -10_|

0,7138-10“ 2

0,9147-10“ 2

0,03125

0,1094-10~1

0,6976-10-5

0,7513-10-4

0,1128-10—4

0,2442-10-4

 

88

Принципам 1—4 удовлетворяет только схема а, для схем б, в, г, (3 нарушены принципы 1,2,3, 4 'соответственно с соблюдени­ ем остальных принципов. Ис­ ключение составляет схема в, для которой, кроме принци­ па 2, 'нарушен также принцип 1, 'Что приводит, как видно из численных данных, к дополни­ тельным потерям по сравне­ нию со схемой б во всем диа­ пазоне X. Схема г хуже схе­ мы а из-за нарушения прин­ ципа 3: труппы имеют неоди­ наковое число общих линий, .в отличие от схемы а. Схема д хуже схемы а из-за наруше­ ния принципа 4: в ней линии размещены не в одной верти­ кали, как в схеме а, а со сдви­ гом, что порождает добавоч­ ные потери.

В качестве вывода по это­ му численному примеру мож­ но высказать гипотезу, что равномерные схемы, т. е. пос­ троенные согласно принципам 4—4, оптимальны в диапазоне больших и средних значений параметра X, а в диапазоне ма­ лых значений X может быть выгодным иарушение принци­ па 1. Что же касается прин­ ципов 2—4, то их нарушение при любых значениях X увели­ чивает вероятность потерь.

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

Схемы, приведенные в табл. 5.3, построены следующим об­ разом:

89)

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