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

УДК 004.422.635.5

ПРОГРАММНАЯ РЕАЛИЗАЦИЯ МАТЕМАТИЧЕСКИХ ОПЕРАЦИЙ МУЛЬТИПЛЕКСИРОВАНИЯ И ДЕМУЛЬТИПЛЕКСИРОВАНИЯ ПОТОКОВ ДЛЯ СЕТЕВЫХ МОДЕЛЕЙ

© 2009 Н. Ф. Бахарева Поволжский государственный университет телекоммуникаций и информатики, г. Самара

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

Характеристики распределения потоков - математические ожидания и дисперсии, мультиплексирование и демультиплексирование потоков, аппроксимация законов распределений и потоков, уравнения равновесия потоков.

Введение

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

Знание (прогнозирование) характеристик потоков важно также для оптимального или близкого к нему управления ими для ограничения загрузки буферов узлов коммутаций (УК), каналов связи и согласования скоростей передачи и приема информации между узлами источник - адресат и т.д.

В настоящее время не существует аналитических методов для точного определения характеристик распределений потоков в сетевых моделях, кроме их средних значений.

1. Постановка задачи и подход к ее решению

Пусть мы имеем открытую сетевую модель с матрицей вероятностей передач Р={р. }, (І І = 1,...,и), где р - вероятность

. ] У

того, что заявка, покидающая узел Б. , поступит в узел Б. . Для начала пусть узел представляет собой одноканальную систему 01/0/1 с бесконечной очередью, для которой определены числовые характеристики слу-

чайного времени обслуживания: т^ - среднее значение и П - дисперсия времени обслуживания. Для внешнего потока задана совокупность средних значений т0. и дисперсий По времени между соседними заявками

рекуррентного потока, входящего в узел Б.. В последующем узел может быть представлен как СМО с конечной очередью с потерями, а также с конечной очередью без потерь.

Для декомпозиции такой модели на отдельные узлы на уровне средних значений и дисперсий времен поступления и обслуживания заявок при произвольных законах их распределений не существует точных методов. Во многих случаях, например в [1, 2], пользуются только уравнениями равновесия

потоков на уровне их интенсивностей X .. Такой подход при произвольных потоках в сети МО означает описание случайного потока событий только его средним значением, т. е. математическим ожиданием без учета моментов высших порядков. Как известно, случайный поток событий на практике чаще всего определяется его такими характеристиками, как математическое ожидание, дисперсия и интервально-корреляционная функция. Поэтому учет дисперсий (вторых

моментов распределений) интервалов времен существенно может улучшить результаты расчетов. Описание потоков на уровне двух первых моментов распределений интервалов времен означает их аппроксимацию непрерывным гауссовским процессом с соответствующими характеристиками.

Если рассматривать структуру отдельного узла с номером i сетевой модели, то на входе /'-го узла сходятся (агрегируются) разреженные (демультиплексированные) потоки от других узлов с номерамиj=l,...,n сети.

Решением системы уравнений (1.1) равновесия потоков относительно интенсивностей 1. потоков на входе и выходе каждой СМО сети определяем средние значения интервалов времен между соседними заявками

т = 1- для каждого потока в сети:

4 = Лш + ЕМ-

j=1

где Л0/ - интенсивность потока извне в /-й узел.

Для вывода уравнений равновесия потоков относительно дисперсий времен между соседними заявками в потоках нужно знать дисперсии выходных потоков

Рвыхвсех узлов, а также формулы для определения дисперсий мультиплексированных и демультиплексированных потоков. Для этого предварительно докажем следующие утверждения.

Утверждение 1. Пусть мы имеем точку мультиплексирования потоков (точка А на рис. 1 или рис. 2), где сходятся два независимых потока заявок с параметрами: Т = 1/1 / (/ = 1, 2) - среднее время между соседними заявками в потоке /, РТ - дисперсия этого

же времени.

Тогда среднее значение и дисперсия времени в суммарном потоке между соседними заявками

*2 = Т^/

^ =

(Т, °тх\

(т2,0Г1),/ 'А

Рис. 1. Мультиплексирование (агрегирование) потоков

Доказательство. Пусть М(/) означает число событий за время /. Тогда среднее значение потока N = ?/т , где г - среднее время между событиями в потоке М(/). Аппроксимируем дискретный поток М(/) непрерывным диффузионным процессом х(/). Процесс х(/) при больших / будет гауссовским процессом со средним значением / /т и дисперсией РТ  1/т [1]. Тогда сумма двух независимых потоков (/) = М1(/) + М2(/) так-

же будет гауссовским случайным процессом со средним / /т2 = / /т1 + / /Т2 и дисперсией

РтЕ " //Т2 = Р*1  / /*1 + Р*2  / / *2 . Отсюда можно записать следующие равенства: 1/тЕ = 1/Т1 + 1/ Т 2 - для среднего времени между соседними событиями в суммарном

потоке и От2 = (рп /т13 + Бт2 /т23)-Т2 - для дисперсии того же времени. Из последних равенств уже следует справедливость выражений (1.2) и (1.3). Таким образом, эти равенства получены на основе диффузионного приближения дискретных процессов на уровне двух первых моментов распределений интервалов времен. Утверждение 1 доказано.

Замечание 1. Допустимость такого приближения можно обосновать следующим образом. Теория сетей массового обслуживания основывается на двух основных допущениях: 1) статистическая независимость потоков в сети МО; 2) рекуррентность входных в СМО потоков. Заметим, что формула

(1.2) точна для любых независимых потоков, а формула (1.3) лишь приближенная вследствие допущения рекуррентности суммарно-

го потока [3, 4]. Это подтверждают и результаты имитационного моделирования.

Замечание 2. На основании полученных равенств легко доказывается справедливость утверждения о том, что сумма нескольких пуассоновских потоков даёт снова пуас-соновский поток.

Замечание 3. Для проверки справедливости выражений (1.2) и (1.3) применялось имитационное моделирование. Заявки генерировались по экспоненциальному, равномерному, нормальному и гиперэкспоненци-альному законам распределений времени между заявками. Результаты моделирования, приведенные в табл. 1, показывают, что формула (1.2) справедлива для любых законов распределений потоков, а формула (1.3) точна только для экспоненциального закона распределения. При этом, если коэффициенты вариаций потоков меньше 1, то формула (1.3) занижает дисперсию агрегированного потока. В случае, когда коэффициенты вариаций больше 1, формула (1.3) обеспечивает очень хорошие результаты. Это связано с тем, что моделирование потоков с таким коэффициентом вариации проводилось с использова-

нием гиперэкспоненциального распределения. В таблице 1 через т и Вт обозначены

—* о 2

теоретические моменты, а через т и дт -соответствующие статистические оценки по результатам имитации.

На рис. 2 приведены гистограммы распределений результирующего потока для некоторых законов распределений составляющих потоков, приведенных в табл. 1.

Результаты имитационного моделирования формул мультиплексирования потоков

(1.2) и (1.3), приведенные в табл. 1, показывают, что они точны для потоков, коэффициенты вариаций которых больше либо равны 1. Гистограммы распределений, приведенные выше на рис. 2-4, показывают, насколько сильно меняются законы распределений потоков при их мультиплексировании, когда они отличны от пуассоновских.

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

Таблица 1

* 1 t-i"1 l^-1 D1/ s’ — , —* Т2 / Т2 D 2/ S 2 — / —* ТХ 't'Z / Si

г 1аконы распределения - равномерный

0,5/0,499 0,083/0,083 0,5/0,50 0,083/0,083 0,25/0,250 0,021/0,037

Законы распределения - равномерный (0:1) и экспоненциальный с параметром X =2

0,5/0,499 0,083/0,084 0,5/0,497 0,25/0,247 0,25/0,248 0,042/0,046

Законы распределения - экспоненциальный с X =2 экспоненциальный с X =2

0,5/0,50 0,25/0,251 0,5/0,499 0,25/0,249 0,25/0,253 0,0625/0,0628

Законы распределения - гиперэксп-ый и гиперэкспон-ый с функцией плотности /(0 = (2р2 /т )ехр{-2pt/т} + [2(1 - р)2 /т]ехр{-2(1 - рУ/т}, гдер=0,8873, т =1,0.

1,0/1,023 4,0/4,026 1,0/1,011 4,0/4,012 0,5/0,502 1,0/1,013

Законы распределения - гиперэксп-ый и гиперэкспон-ый с функцией плотности f (t) = (2p2 /t )exp{-2pt/t} + [2(1 -p)2 /t]exp{-2(1 -p)t/t}, гдеp=0,9472, t =1,0.

1,0/0,994 9,0/8,918 1,0/1,014 9,0/9,442 0,5/0,496 2,25/2,351

Рис. 2. Гистограмма распределения интервалов времени в суммарном потоке

Рис. 3. Гистограмма распределения интервалов времени в суммарном потоке

2,872 I Т7Г[

1 Я?1 I

О 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

Рис. 4. Гистограмма распределения интервалов времени в суммарном поток.

Рис. 5. Демультиплексирование потока

и дисперсия времени между соседними событиями в просеянном потоке

т р =т / р

°тр = Вг/р + 7-2(1 - р)/р2 •

г-1

Определим математическое ожидание величины тр с учетом того факта, что поток рекуррентный, т.е.

М(тг) = М(т), М(тг2) = М(т2):

^ 1

М(тр) = М(т) • рЕг • ц1 = рМ(т) х

1=1

х (ц+ч1+^ + ...)Ц = рм (т)(т^)Ц =

1 - ц

= рМ

(1 - ч)2

= М (т)/ р.

Доказательство. Представим поток событий с параметрами ( т, ), где т, -

среднее значение и дисперсия времени между соседними событиями как последовательность случайных точек {.} на оси времени 0^. Случайные интервалы времени между ними обозначим соответственно через т1,т2,т3,.... При разрежении потока 1 случайный интервал времени в потоке 2 - т р равен: т - с вероятностью р; (т +12) - с вероятностью рц; (т +Г2 +тз)- с вероятностью рц2 и т.д. Запишем выражение для величины

тр следующим образом:

т с вероятностью р,

(т1 + т2) с вероятностью рц,

(т1 +т2 +т3) с вероятностью рц

(т1 + т2 +... + тг) с вероятностью рц

Для определения дисперсии вычис-

2

лим второй начальный момент М( р):

М(тр) = рЕ км(т2) + (к2 -к)(М(т))2

к=1

X цк-1 = р Е \кВт + к 2(М (т))2]цк-1 =

к=1

= рВт Е кцк-1 + р(М (т))2 ЕЕ к У-1 =

к=1

к=1

Здесь вероятность ц=1-р. Заметим, что Р - преобразование потока является случайным его разрежением в отличие от детерминированного просеивания при получении потока Эрланга из простейшего потока.

= От /р + (1 + ц)(М(т))2/р2.

Из полученных равенств для М (т2) и

М (т р ) следует справедливость выражений

(1.4) и (1.5). Утверждение 2 доказано. Следствие. Полученное из формулы

(1.5) выражение для квадрата коэффициента

2 2

вариации просеянного потока ср = р • с + ц

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

Замечание. Для проверки справедливости выражений (1.4) и (1.5) также использовалось имитационное моделирование. Заявки генерировались по экспоненциальному, равномерному и нормальному законам распределений времени между заявками. Для каждого закона распределения генерировалось по 10000 заявок. Результаты имитаци-

1

х

онного моделирования, приведенные в табл. 2, полностью подтверждают справедливость утверждения 2 для любого закона распределения.

На рис. 6-8 приведены гистограммы распределений разреженных потоков для некоторых законов распределений, приведенных в табл. 2.

Результаты имитационного моделирования формул демультиплексирования потоков (1.4) и (1.5), приведенные в табл. 2, показывают, что они являются точными для любых законов распределений. Гистограммы распределений разреженных потоков, приведенные на рис. 6-8, показывают, что демультиплексирование потока также существенно влияет на характер распределения исходного потока, кроме пуассоновского.

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

2. Повышение точности математической модели мультиплексирования потоков

В работе [5] приведено выражение для функции распределения интервала времени

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

^1 и 2 •

Таблица 2

Закон распределения - экспоненциальный с параметром 2 =2; р=0,2; ^=0,8

т/ т* О / — / —* тр/т р Ор / 52р — /—* ^ / tq Ъщ / 5 2

0,5/0,498 0,25/0,249 2,5/2,491 6,25/6,121 0,625/0,622 0,39/0,389

Закон распределения - экспоненциальный с параметром 2 =0,5; р=0,2; ^=0,8

2,0/2,01 4,0/4,01 10,0/9,96 100,0/99,33 2,5/2,51 6,25/6,308

Закон распределения - равномерный с параметрами о=1, Ъ=3; р=0,2; д=0,8

2,0/2,001 0,333/0,333 10,0/9,99 81,666/84,23 2,5/2,51 1,666/1,703

Закон распределения - нормальный с параметрами тт =2, ат =0,5; р=0,2; д=0,8

2,0/1,993 0,25/0,249 10,0/9,98 80,25/80,59 2,5/2,497 38,75/38,26

¥

т

+В2 (г) | В1 (и)ёи}

где (г) = 1 - (г), (=1, 2), а В(г) - функ-

ция распределения интервалов времени между событиями в потоке.. Используя выражение (2.1), определим функцию плотности для результирующего потока для последующего вычисления первых двух моментов его распределения. Для этого введем обозначения: ¥¥

81 (г) = I [1 - В (и)]<^и, 8 2(г) = Ш - В2 (и)] <^и .

г г

Заметим, что %1(0) = т1з %2(0) = т2, т.е. эти функции в точке 0 равны соответствующим средним значениям интервалов времен в потоках. Несложно показать, что тогда функция плотности

2 2

/(г) = *1 (г) = —[81 (г)• %2{г)]'.

2*1 + ^2

Математическое ожидание, т. е. среднее значение интервала между событиями в результирующем потоке

¥ 2 2 ¥

Ч = I г/1(г¥г = - 1 2- Iг-[81(гV 82 (г)]"^г =

0 \ + ^2 0

= 1112 [г ( (г) (г))' ¥ ¥

[г •(81(г) • 82(г))' ¥ -1 [81(г) • 82(гШ.

Рис. 6. Гистограмма распределения интервалов времени в разреженном потоке

Рис. 7. Гистограмма распределения интервалов времени в разреженном потоке

Рис. 8. Гистограмма распределения интервалов времени в разреженном потоке

Здесь первое слагаемое в правой части равно нулю при определенных ограничениях на функции распределения В.(г), а именно, для функции из семейства экспоненциальных распределений. Тогда

Ту = -

X1X2 Х1 + Х2

[ gl(t )• g 2 (t )]|¥ =

= о+

1

1

Х1 + Х2

Т1Т2 =

Х1 + Х2

что полностью подтверждает справедливость формулы (2.1). Результаты, полученные по формуле (2.2), совпадают с результатами аналогичной формулы (1.2).

Определим теперь второй начальный момент распределения интервала т 2 для вычисления дисперсии этой случайной величины:

¥ 2 2 ¥

М(Т!) = 1 г2/! (г)ёг=-12-1 г ^(гЫг)]'*=

0 2 + 2 0

2 2 ¥

= т-12- {г 2[81(г)82(г)]' |¥-21 г[ а(г) %2(г)]'^г}.

21 1 22 0

Первое слагаемое в правой части последнего выражения при тех же ограничениях на функции В(г) опять будет равно нулю. Тогда, вычисляя второй интеграл по частям, получим

М(т2) =-2[г%1(г)%2(г) |¥ -¥%1(г)%2(г)^г] .

21 1 22 0

Здесь первое слагаемое снова равно нулю. Тогда окончательно имеем:

22 ¥

2) = 2 ^ I g1(t) g 2 (t) dt,

M (T s2) = 2

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

разом, в общем случае дисперсию величины ts -интервала времени между событиями результирующего потока нельзя будет выразить в виде элементарной функции от дисперсий и математических ожиданий составляющих. Этот интеграл можно вычислить только при конкретных выражениях функций распределений F(t). Учитывая тот факт, что в сетевых моделях и в самих реальных сетях мы не знаем точных законов распределений потоков, остается единственно возможный путь для вычисления этого интеграла через элементарные функции - это аппроксимация функций g(t) (/=1,2) на уровне двух первых моментов распределений интервалов времени. Эти моменты можно на практике определить путем съёма трафика с помощью программно-аппаратных средств измерения трафика в узлах сети.

Рассмотрим конкретные примеры распределений. В качестве первого примера возьмем два экспоненциально распределенных потока с параметрами:

F(t) = 1 - е-2 ',F2 (t) = 1 - е

Тогда по формуле (2.3) получим

М(Т 2) = 2 /(2 + 2 А

следовательно, дисперсия величины т2 -

О = 1 /(\+ 22 )2, что полностью совпадает с выражением (1.3), если в него вместо дисперсий О подставить их значения 1 / 22.

Это означает, что при мультиплексировании потоков, распределенных по экспоненциальному закону, снова получается пуассоновс-кий поток.

В качестве следующего примера рассмотрим два независимых потока событий, распределенных по равномерному закону на интервале (0;1). Тогда дисперсия величины

т2 по формуле

2 1

D(тs) = M(ts) - —.

(2.4) D =

11 1 (1Y 3

21 [ |( 1 - u)du |( 1 - u)du]dt -|^ 4 J = 80.

Под интегралом в выражении (2.3) сто-

Формула же (1.3) в этом случае дает

ит произведение двух функций, и, таким об- результат 1/48, следовательно, она занижает

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

С2 . <1 и С2 . >1.

3. Определение параметров аппроксимирующих функций распределений

Как известно из теории массового обслуживания, такими функциями распределения в первом случае является гипоэкспонен-циальное распределение, а во втором - ги-перэкспоненциальное [1, 5].

Функция распределения в первом случае

Если же при этом один поток будет иметь коэффициент вариации меньше 1, а другой - больше 1, то в таком случае функции 8.(г), очевидно, будут скомбинированы из выражений (3.3) и (3.4). Параметры иско -мых аппроксимирующих функций распределений (3.1) и (3.2) подберем, используя метод моментов, приравняв первые два момента данных распределений. Математическое ожидание случайной величины, распределенной по закону (3.1), равно

¥

Т* = I гехр{-(г-т 1 )/Т}2}/т}2Л =

т11

= - г ехр{-(г - Тп )/ Т12 } |¥ +

¥

+ I еХР{-(г - Тп )/ Т. 2 }Лг = Т. + Т. 2 .

т11

Проинтегрировав дважды по частям, найдем дисперсию распределения (3.1):

* I 0, і < т;1

Р} [1-ехр{-

Б* = 7 (7 = 1,2).

Используя метод моментов, запишем:

а во втором

рі(і) =1 - Р,ехР(~ 2РІ /Ті) -

-

Возникает задача выбора параметров распределений (3.1) и (3.2). Для этого опре-

¥

делим функции Еі(і) = |[ 1 -¥і(і)]йі, под-

і

ставив в их выражения (3.1) в случае Сл< 1, 7=1, 2:

81(і) = -

Т +Т -і і < Т

1]1^ 2 1 - 1]1

71-7 2.

-71-

Т71 + Т7 2 = Т7

.Т? 2 = А/

Отсюда параметры функции распределения (3.1) равны:

Т7 2 =4Вт,

Аналогично проделаем те же операции с функцией распределения (3.2). В этом случае функция плотности

2 Р

2 РІ

/(і) = ехР{ --7} +

Т7 Т 7

В случае Сл 7 >1 подставим (3.2) в выражение для функции g7(і):

%і(і) = Т7 {еХР[ -2 Р/ /т7] + еХР[ -2( 1 - Р7 )г / Т7 ]} / 2.

2( 1 - Р,)2

єхр{ -

2( 1 - Р1)і

}.

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

равно т* = тг Дисперсию этой величины найдем, дважды проинтегрировав по частям

М(т I) = 2-р+у- Г gl(t) • (t)dt-к, + ^0

интеграл

си

/1 ги 0)Л :

= 2-

к1к2 к, + к

11

{ | (Т11 + Т12 - t)(T2l +'

б: = 7 1 / 2 Р1 + 1/2(1-р1)]-т

Параметр распределения

Ь21

+ | Ь(т2! + т22 - t)exp[-(t - тц)/т12^ + т„

¥

+ Г т т exp[-(t-т )/т -(t-т )/т ]dt}.

Л 12 22 11 12 2Г 22-' '

21

Из (3.6) определим значение вероятности pJ:

В случае т > т

21

1

(3 8) М(т2) = 2кк-{ Г (тп + т12 - 0^ + т22 - (№ +

4 ’ к1 + к2 0

тп

+ 1 т22 (тц + ти-1)ехр[-(г - т21)/т 22]dt +

т 21 ¥

+ Г т т ехр[-(?-т )/т -(I-т )/т ]dt}.

3 12 22 11 12 21 22 >

4 2

Таким образом, параметры функций распределений ¥* ^), аппроксимирующих

законы распределений Р(р) составляющих результирующего потока, полностью определены для всех случаев С к J £ 1 и Ск J >1. Тогда, подставив функции gj(t) 7=1,2 с одно-

11

т < т :

11 21 •

В случае же равенства значений

т =т второй интеграл в выражениях (4.1) значно определенными их параметрами в вы- 11 21

ражения (2.3), (2.4), после вычисления всех и (4.2) будет равен 0. Обозначив интегралы в

интегралов можем определить дисперсию правых частях выражений (4.1) и (4.2) через

интервала времени мультиплексированного I I I запишем их значения в случае

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

4. Определение характеристик распределения мультиплексированного потока

Как было показано выше, среднее значение интервала времени в результирующем

потоке легко определяется по формуле (2.2) 12 = т^(т + т - т - т ) +

независимо от значений коэффициентов ва-

риаций распределений составляющих пото /_ ” ' "г/_ ~ '/_ 1

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

чае Ск , £ 1. Для этого в формулы (2.3) и (2.4) I = т2 т2 ехр[(т -т )/т ]/(т +т ).

'> 3 1Л00-*|-'77 07'/70-/'70

подставим выражение (3.3) с найденными ранее параметрами распределения (3.5). Тогда

I = (т / 2)(т + т -т /3-т ) + 1 '11 21 22 11 12

+т т (т + т ),

11 1Г 21 227

+т2 (т -т )ехр[-(т -т )/т ],

1Г 12 22 ' 21 1Г 12

12 22

11 2Г 1Г ' 12 22'

В случае 11 > 21 эти интегралы рав-

в случае <

11 21

Г

2

21

г

її = (т221/2)(т„ + т12 - т21 /3 -+

+Т21 Т22(Т11 + Т12 )’

ї2 = т2 (т +т -т -т ) +

2 22 11 12 21 22

+т2 (т -т )ехр[-(т -т )/т ],

22 22 1Г ^ ' 11 2Ґ 22"'

ї3 = т2 т2 ехр[(т -т )/т ]/(т + т ) .

3 12 22 21 1Ґ 22 ' 12 22у

Здесь параметры распределений т]к

(/=к=1,2) определяются выражениями (3.5). Подставив полученные значения интегралов (4.3) в (4.1), а - (4.4) в (4.2) для второго начального момента распределения интервала времени между событиями результирующего потока, определяем дисперсию распределения этого интервала:

Б(т, ) = М(ті )-

1

+1 ехр[-2(— + 2(1 р 2) )і ]йі +

о

¥

Ті То

+ І ехр[-2( 2(1 Р|) + Щі ]йі +

о

¥

Т1 Т2

т2т2

I '= т1 т2______________ I '=_

1 1 * 1 ->

т1 т2

1 8(р1т2 + р2т1) ’ 2 8[р1т2 + (1 - р2)т!] ;

13 =

т1 т2

8[—2т1 + (1 - —1)т2 ]

-2 2

14 =

8[(1-—1 )т2 + (1 - —2)т1 ] '

Тогда дисперсия величины т2 при ги-

перэкспоненциальном распределении составляющих результирующего потока будет равна

і і 1

і = 2-^Мї/ + г2 +13 +14) — 1

ті

= 2^^

V Vі

Теперь те же операции выполним для случая гиперэкспоненциального распределения составляющих. Для этого функции

gJ(t), определяемые выражением (3.4) с параметрами распределения (3.7) и (3.8), подставим в (2.3).

к к ¥

М (т2)=2.12 I ^(t) • g2(t) ^ =

к + 0

= 2_к^тт^{¥ ехр[-2(Р+-Р2)t ]dt +

к1 + к2 4 о :1 :2

+ | ехр[-2( 2(1 —1) + 2(1 —2) )і ]^і}.

Обозначим через I/ , 1'2, !3' и Г4 интегралы, стоящие в фигурных скобках, и определим их значения:

(к + к2)2'

Далее рассмотрим последний случай, когда одна составляющая результирующего

потока имеет коэффициент вариации СкJ. < 1, а вторая - Ск>1. Не умаляя общности, в качестве функции g1(t) возьмем функцию

(3.1), а в качестве функции g2(t) - функцию

(3.2) с известными уже параметрами. Тогда, подставив эти функции в (2.3), получим

л л ¥

М(т2) = 2-кк. Гgl(t)• g2(^Л =

к1 + 2 О

т11

= | (т11 + т12 - t)[eXР(-2Р2t / т2 ) +

0

+ ехр(-2( 1 -р2^/т2)](т2 /2^ +

¥

+ | т-12 ехр[-(t-т„)/т-12]х

ти

х[ехр(-2р^ / т2 ) + ехр(-2( 1 - р2 ^/ т2)] х х(:2 / 2 ^.

Обозначим первый интеграл через I'',

второй - 12'.

Выражение I' ' найдем после интегрирования по частям:

2

К - (t11 + t12 >22 /(4P2 ) - f Т12Т22 /(4P2 ) --t23 /(8P2 )] eXP(-2P2t\ 1 / 12 ) - 13 /(ВP22 J + +(t11 + Ті2 ^ /f 4( 1 - P2 )] - Т12Т22 /f 4( 1 - P2 )] --t23 /f В( 1 - P2 )2 ] eXPf -2( 1 - P2 )Тіі / 12 ] --т23 /f В( 1 - p2 )2 ],

Т\2Т2 reXp

Т2 + 2 P 2і 12

].

1

+ еХР[—2( 1 ~ Р2 )Т11 7 Т2 ]

Т 2 + 2(1 — Р2 )Т12 Тогда

д Е= 2-^^ (/;'+ г2' ) — 2

тЕ V 1 2 (\+ х2)2 •

Схема алгоритма операции мультиплексирования двух потоков, реализованная по вышеполученным формулам, представлена на рис. 9.

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

чений (2.1) можем записать уравнения равновесия относительно их дисперсий:

Д « = Д

Здесь Д - операция вычисления дисперсии, а (П — и * П}1) означает попарное мультиплексирование выходных потоков от (/-1)-го ( П ^-и ) и у-го узлов ( П}1), поступающих на вход /-го узла после их р преобразования. Тогда дисперсии этих потоков равны:

1

1 - Pj

) j=1,2,---,n);

DП ji (Di вых ' л 2

' Pi ' Pi •

П0І - дисперсия потока П0І от внешнего источника, поступающего на вход /-го узла. Решение уравнений (1.1) и (4.10) позволяет декомпозировать сетевую модель на отдельные узлы на уровне двух первых моментов распределений потоков для последующего расчета их характеристик.

Таблица 3

_ * Т/Ті D і/ Si! — / —* 12 / X 2 Dx 2/ S22 — / —* Dffi /

г іаконьї распределения - равномерный

0,5/0,499 0,083/0,083 0,5/0,50 0,083/0,083 0,25/0,250 0,029/0,037

Законы распределения - равномерный (0:1) и экспоненциальный с параметром X =2

0,5/0,499 0,083/0,084 0,5/0,497 0,25/0,247 0,25/0,248 0,041/0,046

Законы распределения - экспоненциальный с X =2 экспоненциальный с X =2

0,5/0,50 0,25/0,251 0,5/0,499 0,25/0,249 0,25/0,253 0,0625/0,0628

Законы распределения - гиперэксп-ый и гиперэкспон-ый с функцией плотности f (t) = (2p2 / t )exp{-2pt /t} + [2(1 - p)2 /r]exp{-2(1 - p)t /t}, гдеp=0,8873, t =1,0.

1,0/1,023 4,0/4,026 1,0/1,011 4,0/4,012 0,5/0,502 0,625/0,999

Законы распределения - гиперэксп-ый и гиперэкспон-ый с функцией плотности /(0 = (2р2 /т )ехр{-2рґ/т} + [2(1 - р)2 /т]ехр{-2(1 - р)ґ/т}, гдер=0,9472, т =1,0.

1,0/0,994 9,0/8,918 1,0/1,014 9,0/9,442 0,5/0,496 1,25/2,351

Законы распределения - гиперэксп-ый с функцией плотности /

1,0/1,001 4,0/3,996 1,0/0,996 0,333/0,333 0,5/0,500 0,204/0,246

т12=е!р/е!з;

ег=0;

s1=sqrt(Dv1); М1=1/є!1-з1; t12=s1;

s2=sqrt(Dv2); t21 =1/e!2-s2; t22=s2;

Э11={т225*

^21+122-411/3-

П2Ж11Ч12*д2^22);

Si2=t122*(t21+t22-t11-

t12)+t122*(t12-t22)*

exp(-(t21-t11)/t12);

Si3=t122*t222*exp((t11-

t21)/t12)/(t12+t22);

Si1=(t112/2)*(t21+t22- П1/3-П2) +П1Ч12*^2^22);

Si2 =0;

Si3=t122*t222/

(t12+t22);

Si1=(t212/2)*(t11 + ^2^21/3^22) +t21*t22*(t11+t12);

Si2=t222*(t11+t12-t21-t22)+t222*(t22-t12)* exp(-(t11-t21 )/t22);

Si3=t122*t222*exp((t21 -t11 )/t22)/(t12+t22);

Ss=si1+si2+si3;

s1=sqrt(Dv12-1/e!14) s2=sqrt(Dv22-1/e!24) s3=Dv1+1/el12 s4=Dv2+1/e!2•

t110=e!1*(s3-s1)/2;

t210=e!2*(s4-s2)/2;

t120=e!1*(s3+s1)/2;

^220=e!2*(s4+s2)/2;

p11=1/(2*e!1*t110); p12=1/(2*e!1 *t120);

p21=1/(2*el2*t210);

p22=1/(2*e!2*t220);

Si10=t1102*t2102*p11

*p21/(t110+t210);

Si20=t1102*t2202*p11

*p22/(t110+t220);

Si30=t1202*t2102*p12

*p21/(t120+t210);

Si40=t1202*t2202*p12

*p22/(t120+t220);

2

s4=Dv2+1/el2

t11=1/e!1-sqrt(Dv1);

t12=sqrt(Dv1);

t210=e!2*(s4-s2)/2;

t220=e!2*(s4+s2)/2;

p21=1/(2*e!2*t210);

p22=1/(2*e!2*t220);

Si11=(t11+t12)(p21* t2102+p22*t22CP)-p21* t2103-p22*t2203+p21 *t21C2*(t210-t12)* exp(-t11/t210) +p22*t2202*(t220-t12) *exp(-t11/t220);

Si21=p21*t122*t2102/ (t210+t12)*exp(-t11/ t210)+p22*t122*t2202/ (t220+t12)*exp(-t11/ t220);

Ss=si10+si2C+si30

+si40;

Ss=si11+si21;

ep=e!1; є!1=є!2; e!2=ep;

ep=Dv1; Dv1=Dv2; Dv2=ep;

+

Возврат

Рис. 9. Схема алгоритма процедуры мультиплексирования потоков

183

Заключение

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

Полученные математические результаты по мультиплексированию двух различных потоков реализованы в виде процедуры МиШрІ с соответствующими параметрами. Схема алгоритма этой процедуры приведена на рис. 9.

В табл. 3 приведены результаты имитационного моделирования по мультиплексированию потоков в сравнении с результатами расчетов по полученным математическим моделям во второй части статьи. Анализ результатов таблиц 1 и 3 показывает, что формула мультиплексирования потоков (1.3) точнее для потоков, коэффициенты вариаций которых больше или равны 1, чем формула

(4.7). И наоборот, формула для вычисления дисперсии результирующего потока в случае, когда коэффициенты вариаций составляющих потоков меньше 1 - (4.5), а также в смешанном случае (4.9), точнее формулы (1.3).

Библиографический список

1.Клейнрок Л. Вычислительные системы с очередями: Пер. с англ./ Под ред. Б.С. Цыбакова. - М.: Мир, 1979. - 597 с.

2. Вишневский В. М. Теоретические основы проектирования компьютерных сетей. - М.: Техносфера, 2003. - 512 с.

3.Тарасов В. Н. Вероятностное компьютерное моделирование сложных систем. -Самара: Самарский научный центр РАН,

2002. - 194 с.

4.Тарасов В. Н., Бахарева Н. Ф. Организация интерактивной системы вероятностного моделирования стохастических систем // Известия Самарского научного центра РАН,

2003. - №1 - С. 119 - 126.

5. Шнепс М. А. Системы распределения информации. Методы расчета. Справочное пособие. - М.:Связь, 1979. - 342 с.

6. Овчаров Л. А. Прикладные задачи теории массового обслуживания. - М.: Машиностроение, 1969. - 324 с.

References

1. Kleinrock L. Computing systems with queues: the Lane with English / Under the editorship of B.S.Tsybakova. - М: the World, 1979.

2. Vishnevsky V.M. The theoretical of a basis of designing of computer networks. - М: the Technosphere, 2003.

3. Tarasov V. N. Probability computer modelling of difficult systems. - Samara: the Samara centre of science of the Russian Academy of Sciences, 2002.

4. Tarasov V. N., Bahareva N.F.Organisation of interactive system of stochastic modelling of stochastic systems//News of the Samara centre of science of the Russian Academy of Sciences, 2003, №1, With. 119 -126.

5. Shneps M. A. Systems of allocation of the information. Calculation methods. The handbook. - M.:Link, 1979.

PROGRAMME REALIZATION OF MATHEMATICAL OPERATIONS OF FLOW MULTIPLEXING AND DEMULTIPLEXING FOR NETWORK MODELS

© 2009 N. F. Bakhareva Povolzhsky State University of Telecommunications and Information Science

The paper deals with the problem of constructing mathematical models of event flow multiplexing and demultiplexing operations. The flows are described at the level of the first two moments of time distribution between events. The results obtained make it possible to write the equilibrium equations for the flows circulating in network models and to decompose the latter into individual nodes to calculate their characteristics.

Flow distribution characteristics - mathematical expectations and dispersions, flow multiplexing and demultiplexing, approximation of distribution laws and flows, flow equilibrium equations.

Информация об авторе Бахарева Надежда Федоровна, кандидат технических наук, доцент кафедры «Программное обеспечение и управление в технических системах», Государственное образовательное учреждение высшего профессионального образования «Поволжский государственный университет телекоммуникаций и информатики». Область научных интересов: программирование, моделирование компьютерных систем. E-mail: vt@ist.psati.ru.

Bakhareva, Nadezhda Fyodorovna, candidate of technical science, associate professor of the department “Software and control in engineering systems”, Povolzhsky State University of Telecommunications and Information Science. Area of research: programming, modeling of computer systems. E-mail: vt@ist.psati.ru.