Научная статья на тему 'МЕТОДЫ ПРОГНОЗИРОВАНИЯ КАЧЕСТВА ОБСЛУЖИВАНИЯ САМОПОДОБНОГО ТРАФИКА В УСТРОЙСТВАХ КОММУТАЦИИ МУЛЬТИСЕРВИСНОЙ СЕТИ'

МЕТОДЫ ПРОГНОЗИРОВАНИЯ КАЧЕСТВА ОБСЛУЖИВАНИЯ САМОПОДОБНОГО ТРАФИКА В УСТРОЙСТВАХ КОММУТАЦИИ МУЛЬТИСЕРВИСНОЙ СЕТИ Текст научной статьи по специальности «Компьютерные и информационные науки»

CC BY
224
51
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МУЛЬТИСЕРВИСНАЯ СЕТЬ / САМОПОДОБНЫЙ ТРАФИК / СРЕДНЕЕ ВРЕМЯ ОЖИДАНИЯ

Аннотация научной статьи по компьютерным и информационным наукам, автор научной работы — Одоевский С. М., Хоборова В. П.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по компьютерным и информационным наукам , автор научной работы — Одоевский С. М., Хоборова В. П.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

METHODS OF PREDICTING SELF-SIMILAR TRAFFIC’S QUALITY OF SERVICE IN MULTISERVICE NETWORK’S SWITCHING DEVICES

Traffic in modern multiservice networks has features appropriate to self-similar random processes. Analytical and imitation methods, with help of which, data flows’ quality of service in multiservice network’s switching devices can be predicted by index of average waiting time for service’s beginning in different self-similar enter traffic conditions, are conducted in this article.

Текст научной работы на тему «МЕТОДЫ ПРОГНОЗИРОВАНИЯ КАЧЕСТВА ОБСЛУЖИВАНИЯ САМОПОДОБНОГО ТРАФИКА В УСТРОЙСТВАХ КОММУТАЦИИ МУЛЬТИСЕРВИСНОЙ СЕТИ»

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

С.М. Одоевский1*, В.П. Хоборова1

военная академия связи имени Маршала Советского Союза С.М. Буденного, Санкт-Петербург, 194064, Российская Федерация *Адрес для переписки: odse2017@mail.ru

Информация о статье

УДК 621.394.343 Язык статьи - русский

Ссылка для цитирования: Одоевский С.М., Хоборова В.П. Методы прогнозирования качества обслуживания самоподобного трафика в устройствах коммутации мультисервисной сети / / Труды учебных заведений связи. 2017. Т. 3. № 3. С. 86-92.

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

Ключевые слова: мультисервисная сеть, самоподобный трафик, среднее время ожидания.

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

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

В настоящее время качество телекоммуникационных услуг мультисервисных сетей нормируется в основном рекомендациями ITU-T, ETSI, 3GPP, IETF и отдельными отечественными руководящими документами. С учетом всех этапов (стадий) предоставления услуг (доступ к услуге, передача информации, завершение сеанса связи) и различ-

ных проявлений задержек, искажений и потерь информации в указанных документах приводится довольно большой перечень оцениваемых и нормируемых параметров (показателей) качества обслуживания. Например, в рекомендации П^-TG.1000 предлагается целая матрица параметров размером 11x7.

С точки зрения вклада в итоговое качество обслуживания со стороны отдельных устройств коммутации мультисервисной сети, реализующих, прежде всего, функции протоколов передачи и коммутации пакетов не выше третьего уровня ЭМВОС, можно ограничиться четырьмя показателями [1]: скорость передачи (V), задержка (т), джиттер (а) и вероятность потерь (р).

Учитывая случайный характер обслуживаемых потоков данных первые два (V и т) из перечисленных выше четырех показателей обычно оцениваются «в среднем» (V и т). Последние два показателя (а и р) сами по себе являются усредняющими. При этом в математических моделях все эти показатели представляют собой соответствующие моменты распределения вероятностей соответствующих случайных величин. Однако в нестационарных условиях с долговременными пульсациями интенсивности трафика все четыре показате-

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

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

Единственным адекватным показателем, не чувствительным (в явном виде) к абстрактным бесконечным всплескам случайной задержки, можно считать четвертый показатель - вероятность потери пакетов р. Именно такой показатель рассматривается в статье [3] для оценки качества обслуживания самоподобного потока данных в маршрутизаторе. Однако на практике для достоверной оценки (измерения) данного «слишком интегрального» и предельного (проявляющегося, когда уже начались потери) показателя требуется довольно много времени, что затрудняет его использование для динамического отслеживания и/или прогнозирования всплесков нагрузки и своевременного принятия упреждающих мер (в виде, например, перераспределения или ограничения входных потоков, выделения дополнительных ресурсов памяти устройств и т. д.).

С другой стороны, поскольку абстрактная бесконечность в известных моделях самоподобного трафика «грозит» в основном только дисперсии (третьему показателю а), то в качестве основного оцениваемого ПКО самоподобного трафика в устройствах коммутации мультисервисной сети целесообразно использовать второй показатель -среднее время задержки т. При этом следует отметить, что первый показатель - скорость передачи пакетов V на выходе устройства коммутации фактически является обратной величиной времени передачи (обслуживания) тоб ~ 1/v, являющейся составной частью времени задержки т = тож + тоб, дополнительно включающей еще время ожидания тож начала передачи (обслуживания) в очереди на входе устройства коммутации.

В общем случае время передачи (обслуживания) Тоб является случайной величиной, которая теоретически может иметь произвольное распределение, в том числе, «с длинным хвостом». Однако на практике в связи с подведомственной ролью данного показателя в конкретной мультисервис-ной сети, которая может им управлять или, по крайней мере, достаточно точно его прогнозировать, в большинстве математических моделей даже с самоподобным входным трафиком допускается время обслуживания считать детерминированным (D - Deterministic) или распределенным по экспоненциальному (характерному для марковских случайных процессов) закону (M - Markovian).

Другая составная часть времени задержки -время ожидания начала обслуживания Тож, с одной стороны, зависит от указанного выше времени обслуживания Тоб (непосредственно связанного с первым показателем - скоростью передачи v), а с другой стороны, непосредственно связана с характеристиками входного потока пакетов, а, следовательно, в первую очередь «ощущает на себе» влияние нестационарности трафика и всех проявлений его свойства самоподобия. Таким образом, наиболее подходящим как с практической (измерительной), так и с теоретической (вычислительной) точек зрения обобщенным показателем качества обслуживания самоподобного трафика в устройствах коммутации мультисервисной сети является среднее время ожидания Тож начала обслуживания.

Для упрощения и большей общности вместо среднего время ожидания Тож удобнее использовать относительную величину этого времени q, нормированную по отношению к среднему времени обслуживания Тоб = 1/^, где ц - интенсивность обслуживания, т. е.:

(1)

_ 1ож

Тоб

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

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

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

_ _ тоб

= я • тож = Я • тоб • q = — • q = р • q. (2)

^вх

При заданной максимальной допустимой длине очереди nq.max, определяемой объемом памяти, выделенном для этой цели в устройстве коммутации, и зная величину щ = р . q (2), можно оценить четвертый ПКО - вероятность потери пакетов p из-за переполнения памяти (даже при неизвестном законе распределения, например, с использованием неравенства Чебышева). Это является дополнительным подтверждением целесообразности выбора в качестве обобщенного ПКО относительного среднего времени ожидания q.

Очевидно, что в «нормальных» условиях, когда очередь на входе устройства коммутации «в среднем» не переполняется, т. е. когда щ < nq.max, должно выполняться условие р < 1. Следует отметить, что именно контроль за возрастанием величины nq (средней длины очереди) позволяет при колебаниях входной нагрузки прогнозировать близость перегрузок. Для ее количественной оценки (с целью упреждающего увеличения nq.max или снижения р) необходимо знать функциональную связь величины q с величиной р, т. е. нужна математическая модель функциональной зависимости q(p), в том числе, с учетом ожидаемых проявлений свойств самоподобия входного трафика.

Таким образом, под методом прогнозирования качества обслуживания самоподобного трафика в устройствах коммутации мультисервисной сети в данной статье понимается метод выбора и использования функциональной зависимости q(p) прогнозируемого относительного среднего времени ожидания q (однозначно связанного со средней длиной очереди nq = р . q) от контролируемой (ожидаемой) интенсивности входной нагрузки р с учетом дополнительной (доступной или предполагаемой) информации об особенностях входного, в частности, самоподобного трафика.

Традиционным дополнительным параметром самоподобного трафика является коэффициент (показатель) Хёрста H, который может принимать значения в диапазоне 0,5 < H < 1. При этом значению H = 0,5 соответствует простейший (несамопо-добный) поток с экспоненциальным распределением интервалов времени Твх между моментами поступления очередных пакетов.

Одной из самых популярных моделей самоподобного трафика является модель типа fbm (от англ. fractal brownian motion - фрактальное броуновское движение) [4]. Для данной модели известны аналитические функциональные зависи-

мости q(p, Н) при детерминированном (О) и экспоненциальном (М) распределении времени обслуживания в одноканальных системах массового обслуживания (СМО) вида /Ьт/О/1 и /Ьш/М/1 [1, 2]. Данные функциональные зависимости приведены в таблице 1.

ТАБЛИЦА 1. Аналитические зависимости q(p,H) и q(p)

fbm/D/1 fbm/M/1 M/D/1 M/M/1

р1/[2(1-Н)] р1/[2(1-Н)] р р

Ч_ 2(1-р)«/(1-н) q = (1 - р)Н/(1-Н) * 2(1-р) * (1-Р)

В той же таблице 1 приведены также аналогичные известные функциональные зависимости q(p) для СМО вида М/О/1 и М/М/1. Нетрудно видеть, что при Н = 0,5 зависимости q(p, Н) для СМО /Ьт/О/1 и /Ьш/М/1 полностью совпадают с соответствующими зависимостями q(p) для СМО М/О/1 и м/м/1.

На рисунке 1 приведены графики зависимостей q(p, Н), рассчитанные по формулам из таблицы 1 при трех значениях Н = 0,5, 0,65 и 0,8.

Я=0.8 Я= 0.65 #=0.5

Рис. 1. Графики зависимости относительного среднего времени ожидания q от нагрузки р при различных значениях коэффициента Хёрста Н, полученные с помощью аналитических моделей

Как видно из рисунка 1, при Н > 0,5 самоподобный трафик в области большой нагрузки, превышающей некоторую величину р^, приводит к значительному увеличению относительного среднего времени ожидания q (в пределе q ^ го при Н ^ 1), а в области малой нагрузки - к уменьшению (в пределе q ^ 0 при Н ^ 1).

Первое явление соответствует известному свойству коррелированного пульсирующего трафика «с длинными хвостами», из-за которого приходится предусматривать в устройствах коммутации накопители гораздо большей емкости qmax, чем это требуется на основании расчетов для традиционных моделей М/О/1 и М/М/1.

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

В [8] получена формула для расчета значения пороговой нагрузки рл на основании совпадения выражений для д(р, Н) и д(р) из таблицы 1 для СМО /Ьт/М/1 и м/м/1 при р = рл :

= (3 -Т5)/2 ~ 0.382. (3)

Выражение (3) соответствует известной «Золотой пропорции», характеризуемой числом Фибия: Ф = (1 + )/2 « 1,618, при этом рь = 1/(1 + Ф).

Приведенные в таблице 1 аналитические выражения для расчета зависимостей д(р, Н) можно считать эталонными для СМО с самоподобным входным трафиком при известном коэффициенте Хёрста Н. Однако идеальный самоподобный трафик типа/Ьт на основе модели фрактального броуновского движения для мультисервисных сетей все же следует считать некоторой математической абстракцией [1]. Да и коэффициент Хёрста для реальных потоков на конечных отрезках времени оценить весьма проблематично [6, 7]. Реальные нестационарные потоки влияют на ПКО более сложным и менее предсказуемым образом.

Как показано в [9], при любой заданной (известной, измеренной или контролируемой) интенсивности входной нагрузки р (0 < р < 1) можно подобрать такое распределение интервалов времени между очередными пакетами, при котором относительное среднее время ожидания д(р) будет иметь любое наперед заданное значение р[0, Отчасти это подтверждает отмеченное выше возможное отсутствие конечного значения математического ожидания у распределения «с длинным хвостом» [7].

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

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

входного потока для корректной оценки конечной величины ПКО в виде относительного среднего времени ожидания q знания только интенсивности входной нагрузки р явно недостаточно. Дополнительным параметром входного потока, позволяющим получить более определенные оценки ПКО, может служить коэффициент Хёрста H. Однако и его знание не позволяет однозначно оценить величину q. Точный вид аналитической зависимости q(p, H) в таблице 1 справедлив только для входного потока типа fbm и СМО fbm/D/1 и fbm/M/1. Но входной поток, как отмечено выше, может не соответствовать модели fbm, хотя может иметь те же значения параметров (р, H), зависимость от которых ПКО q(p, H) будет не такой, как в таблице 1 [5].

Известно, что для моделирования самоподобного трафика кроме модели fbm могут использоваться многие другие распределения с «длинными хвостами» [1-3, 5, 6]. Наиболее популярными являются распределения Парето (P - Pareto) и Вей-булла (W - Weibull) [3, 5]. В таблице 2 представлены основные математические выражения, характеризующие данные распределения, где: F(x) -функция распределения, P(x) - его плотность, Mo -математическое ожидание, Do - дисперсия. В этой же таблице (для сравнения) приведены аналогичные математические выражения для экспоненциального распределения (E - exponential).

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

В частности, представление всех трех функций распределения в виде функций относительно нормированной величину x/k позволяет более явно наблюдать влияние параметра а на величину «хвостов» распределений в области x/k >> 1.

ТАБЛИЦА 2. Характеристики функций распределения интервалов времени между моментами поступления пакетов входного потока

Распределение

E P W

F(x) 1-ехр(-£) fx\-a 1-U 1-exp(-©")

PX) a/Xx-O+i) k\k) КГехрНкЛ

Mo к ak a — 1 kr(i+a)

Do k2 ak2 (a — 1)2(a — 2) kfHHK))

Примечание. Для распределения Парето формулы расчета F(x) и P(x) справедливы только для x > к. Иначе (при x < к) F(x) = 0 и P(x) = 0.

Очевидно, чем больше величина параметра а (в области значений а > 1), тем быстрее убывают «хвосты» распределений Парето и Вейбулла. Причем у распределения Вейбулла а > 1 «хвост» убывает даже быстрее, чем у экспоненциального распределения, которое, как видно из формул в таблице 2, полностью совпадает с распределением Вейбулла при а = 1. Соответственно, при этом совпадают и выражения для расчета математического ожидания Mo = k и дисперсии Do = к2.

При уменьшении параметра а < 1, «хвост» распределения Вейбулла начинает расти. При этом начинают возрастать значения математического ожидания Mo и дисперсии Do, которые в пределе при а^ 0 стремятся к бесконечности, что асимптотически характерно для самоподобных случайных процессов. Для Mo это видно из монотонно возрастающей зависимости гамма функции Г(1+1/а) при увеличении аргумента (1+1/а) ^ го, из-за а^ 0, а для Do это следует из более сильного возрастания первого слагаемого Г(1+2/а) ~ [2/а]! по сравнению с возрастанием второго Г2(1+1/а) ~ ([1/а]!)2 при а^ 0. Но математическое ожидание Mo распределения Вейбулла можно «удержать» от возрастания «на месте» (например, приравняв его математическому ожиданию экспоненциального распределения) при а ^ 0 за счет компенсирующего уменьшения параметра к. Но дисперсия Do при а ^ 0 все равно будет возрастать и стремиться к бесконечности, что асимптотически присуще самоподобным процессам.

Распределение Парето, как видно из таблицы 2, по сравнению с распределением Вейбулла при любых значениях параметра а имеет более длинный «хвост», так как характеризуется обратной гиперболической зависимостью от (х/к)а, а не экспоненциальной. При этом, как видно из таблицы 2, дисперсия распределения Парето стремится к бесконечности при уменьшении значения а > 2 и приближении а ^ 2. А при а < 2 согласно приведенной формуле дисперсия распределения Парето становится отрицательной, т.е. практически неопределенной. В то же время в области значений а> 1 математическое ожидание распределения Парето имеет конечную величину и стремится к бесконечности только при а ^ 1. При а < 1 и дисперсия, и математическое ожидание распределения Парето становятся отрицательными (практически неопределенными).

Из вышесказанного следует, что распределение Парето является классическим представителем семейства распределений, которые могут иметь бесконечные или неопределенные значения дисперсии и математического ожидания [7] и в наибольшей степени подходят для имитации самоподобных процессов. В то же время очевидно, что для моделирования потоков данных с конечной определенной средней интенсивностью Я, обрат-

ной величине среднего интервала времени Твх между моментами поступления очередных пакетов допустимо использование распределения Па-рето только с а > 1. Причем наиболее явное проявления свойства самоподобия следует ожидать в области значений 1 < а < 2, при которых дисперсия имеет неопределенное (отрицательное) значение.

Отмеченные выше диапазоны значений параметра а, при которых у распределений Вейбулла (а е (0, 1]) и Парето (а е (1, 2]) наиболее сильно изменяется «длина хвостов», согласно многим публикациям (в частности, [3, 5, 6]) непосредственно связаны (но инверсно) с диапазоном изменения коэффициента Хёрста Н, отражающего степень самоподобия трафика H е (0,5, 1], описываемого данными распределениями:

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

а = 3 - 2H, (4)

для распределения Вейбулла

а = 2 - 2H. (5)

С учетом выражений (4), (5) и формул в таблице 2 (для расчета математического ожидания Mo = Твх = 1/Я) можно получить аналитические зависимости параметров а и k для распределений Парето и Вейбулла от параметров моделируемого самоподобного потока данных Я и H (таблица 3).

ТАБЛИЦА 3. Взаимосвязь параметров потоков данных и параметров распределений, используемых для моделирования этих потоков

Распределение

Е Р IV

Функция распределения 1-ехр(-£> /Х\-и 1 — (х) 1—ехфК1Л

Параметры потока и распределения к II а = 3 — 2Я; а — 1 аЯ а = 2 — 2Я; 1 /• —

Обратная функция х = Р-1(у) х = — к • /п(у) к х = —;= ТУ х = кТ- 1п(у)

у - случайная величина с равномерным распределением Р(у) = у, уе(0,1]

К сожалению, точные аналитические выражения для расчета прогнозируемого ПКО д(р, И), р = Я / ц, для СМО типа Р/М/1 или P/D/1, а также Ш/М/1 или W/D/1, т. е. при поступлении на вход устройства коммутации потоков данных, описываемых распределениями Парето (Р) и Вейбулла (W), не известны. Поэтому в настоящей работе для оценки указанной зависимости по аналогии с [3, 5] предлагается использовать имитационное моделирование.

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

Пример результатов расчета зависимостей д(р, Н) с использованием имитационных моделей СМО Р/М/1, Р/й/1, Ш/М/1, Ш/й/1, М/М/1,М/й/1 (с использованием формул из таблицы 3) приведены на рисунке 2 совместно с результатами расчета с использованием аналитических моделей СМО /Ьт/М/1,/Ьт/й/1, М/М/1, М/й/1 (с использованием формул из таблицы 2).

б)

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

Корректность разработанных моделей и их программной реализации подтверждается ожидаемым совпадением результатов расчета с использованием аналитических моделей для СМО М/М/1, М/й/1 и результатов расчета с использованием имитационных моделей для аналогичных СМО М/М/1, М/й/1, а также для СМО Ш/М/1, Ш/й/1 при Н=0.5 (см. рисунок 2б) во всем диапазоне (выбранном для вычислительного эксперимента) изменения нагрузки р е [0,1; 0,9] с шагом Др = 0,02.

Распределение Парето подобной сходимостью к экспоненциальному распределению при Н = 0,5 не обладает, что подтверждается несовпадением аналогичных графиков для М/М/1 и Р/М/1, а также М/й/1 и Р/й/1 на рисунке 2а. В то же время при

большом значении коэффициента Хёрста Н = 0,8 распределение Парето приводит к зависимости времени ожидания от нагрузки д(р) более близкой к аналогичной зависимости для потока типа /Ьт (эталона самоподобного процесса), чем распределение Вейбулла.

Следует отметить, что при имитационном моделировании распределений с теоретически бесконечной или неопределенной дисперсией выборочная дисперсия, естественно, все равно оказывается конечной и определенной. Потенциальная бесконечность отдельных реализаций моделируемых случайных интервалов времени х согласно формулам обратного преобразования х = ^ _1(у) в таблице 3 кроется в предельном бесконечном значении функций ^ _1(у) ^ ю при у ^ 0. Однако конечность численного представления значений данных функций не позволяет генерировать числа больше некоторой величины хтах (максимальное число Хтах, используемое в MathCad в качестве условной бесконечности, равно 10307). Это ограничение не позволяет в полной мере воспроизвести самоподобные (фрактальные) процессы с существенным влиянием дальней части «длинных хвостов», уходящих в бесконечность при любом реально возможном объеме выборки N.

При отмеченном конечном диапазоне изменения случайных чисел при имитационном моделировании с ростом объема выборки N наблюдается неустойчивое значение выборочной дисперсии, что характерно для самоподобных процессов и теоретически, и практически. А устойчивое значение выборочного математического ожидания оцениваемой величины ПКО <|(р, Н) удается получить только при очень большом значении объема выборки (особенно при больших значениях коэффициента самоподобия Н). Приведенные на рисунке 2 зависимости ^(р), для СМО Р/М/1, Р/й/1, Ш/М/1, Ш/й/1 при Н = 0,8 удалось «сгладить» (и то, как видно из графиков, не полностью) лишь при N = 2105.

Достаточно гладкий вид полученных путем имитационного моделирования графических зависимостей д(р, Н) для СМО Р/М/1, Р/й/1, Ш/М/1, Ш/й/1 позволяет рассчитывать на возможность их аналитического представления теми или иными эмпирическими формулами. В частности, в [1] для аналитического представления зависимости средней длины очереди от нагрузки при описании входного самоподобного потока данных распределением Парето предлагается использовать формулу Поллачека-Хинчина с учетом предварительно вычисляемой конечной дисперсии, возникающей при ограничении максимальных выбросов случайных интервалов времени некоторой величиной Хтах, отражающей конечность реальных всплесков трафика в мультисервисных сетях.

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

него времени ожидания начала обслуживания (однозначно определяющего среднюю длину очереди) при небольшом значении коэффициента Хёрста (0,5 < Н < 0,75) использовать модели СМО с входным потоком, описываемым распределением Вейбулла, а при большом (0,75 < Н < 1) - распределением Парето.

Список используемых источников

1. Назаров А.Н., Сычев К.И. Модели и методы расчета показателей качества функционирования узлового оборудования и структурно-сетевых параметров сетей связи следующего поколения. Красноярск: Изд-во ООО «Поликом», 2010. 389 с.

2. Крылов В.В., Самохвалова С.С. Теория телетрафика и ее приложения. СПб.: БХВ-Петербург. 2005. 288 с.

3. Лосев Ю.И., Руккас К.М. Анализ моделей вероятности потери пакетов в буфере маршрутизатора с учетом фрак-тальности трафика // Вестник ХНУРЭ. Серия «Математическое моделирование. Информационные технологии. Автоматизированные системы управления. 2008. № 833. С. 163-169.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4. Norros I. On the Use of Fractional Brownian Motion in the Theory of Connectionless Networks //IEEE Journal on Selected Areas in Communications. 1995. Vol. 13. PP. 953-962.

5. Ложковский А.Г., Вербанов О.В. Моделирование трафика мультисервисных пакетных сетей с оценкой его коэффициента самоподобности // Научные труды ОНАС им. А.С. Попова. 2008. № 1. С. 57-62.

6. Агеев Д.В., Игнатенко А.А., Копылев А.Н. Методика определения параметров потоков на разных участках муль-тисервисной телекоммуникационной сети с учетом эффекта самоподобия // Электронное научное специализированное издание - журнал «Проблемы телекоммуникаций». 2011. № 3 (5). URL: http://pt.journal.kh.ua.

7. Колядин В.Л. Распределения с бесконечной дисперсией и ограниченность классической статистики // Радиотехника. 2002. № 2. С. 4-11.

8. Одоевский С.М., Сорокина Е.А., Хоборова В.П. Варианты расчета показателей качества обслуживания информационных потоков в узлах коммутации мультисервисной сети военного назначения // Труды учебных заведений связи. 2016. Т. 2. № 4. С. 92-99.

9. Одоевский С.М., Хоборова В.П. Моделирование пульсирующего трафика на входе инфокоммуникационной системы с заданными выходными характеристиками процесса обслуживания // Современные проблемы создания и эксплуатации вооружения, военной и специальной техники. III Всероссийская НПК. СПб.: ВКА. 2016. С. 293-295.

11. Одоевский С.М., Калюка В.И. Адаптивно-игровое моделирование военных сетей беспроводного абонентского доступа. СПб.: Изд-во Политехн. ун-та, 2017. 342 с.

* * *

METHODS OF PREDICTING SELF-SIMILAR TRAFFIC'S QUALITY OF SERVICE IN MULTISERVICE NETWORK'S

SWITCHING DEVICES

S. Odoevskiy1, V. Khoborova1

telecommunications Military Academy, St. Petersburg, 194064, Russian Federation

Article info:

Article in Russian

For citation: Odoevskiy S., Khoborova V. Methods of Predicting Self-Similar Traffic's Quality of Service in Multiservice Network's Switching Devices // Proceedings of Telecommunication Univercities. 2017. Vol. 3. Iss. 3. PP. 86-92.

Annotation: Traffic in modern multiservice networks has features appropriate to self-similar random processes. Analytical and imitation methods, with help of which, data flows' quality of service in multiservice network's switching devices can be predicted by index of average waiting time for service's beginning in different self-similar enter traffic conditions, are conducted in this article.

Keywords: multiservice network, self-similar traffic, average waiting time.

i Надоели баннеры? Вы всегда можете отключить рекламу.