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

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

П.А. Трещановский

ЗАО “Ангстрем-Телеком”

Московский институт электронной техники

Передача мультимедийного трафика реального времени с соблюдением требований по качеству обслуживания (QoS) является одной из наиболее актуальных задач в современных телекоммуникациях. Ограниченная пропускная способность сетей, большое число абонентов и жесткие требований к параметрам QoS могут становиться причинами падения качества предоставляемой мультимедийной услуги в многоканальных абонентских сетях. Для решения данной проблемы необходимо динамическое управление ресурсами сети с резервирующее необходимой пропускной способности для каждого мультимедийного потока. Широко известными являются применяемые в Интернете методы интегрированного обслуживания (IntServ) [1] и дифференцированного обслуживания (DiffServ) [2]. Эффективность управления сетевыми ресурсами, проявляющаяся в максимально возможном количестве обслуживаемых потоков, во многом зависит от математической модели трафика, на основе которой производится резервирование. В данной работе предложена новая параметризованная модель, позволяющая за счет более полного учета вероятностных свойств трафика добиться повышения коэффициента использования сети. Получение оптимальных параметров модели возможно с помощью предложенного метода расчета на основе стохастической глобальной оптимизации. Одной из наиболее распространенных моделей трафика для управления сетевыми ресурсами и контроля доступа является кривая поступления [3]. Эта модель представляет собой функцию времени а(т), определяющую максимальный объем данных, передаваемый источником трафика за время t. Данная модель позволяет легко характеризовать агрегированные мультимедийные потоки с помощью суммы функций отдельных кривых , а также осуществлять профилирование потоков на входе в сеть, например, с помощью алгоритма “дырявого ведра” [ 4]. Для определения параметров QoS используется дополнительная модель, характеризующая узел коммутации (УК) сети и называемая кривой обслуживания р(т). Кумулятивный выходной поток R(t) такого узла связан с входным потоком A(t) неравенством R(t) > A(t)  р(т), где  обозначает операцию свертки в мин-плюс алгебре. Кривая обслуживания обладает рядом полезных свойств. Во-первых, важнейшие параметры QoS могут быть вычислены на основе кривых поступления и обслуживания с помощью теории сетевого исчисления. В частности, максимальная задержка передачи определяется равенством

dmax = Ка(т),Р(т)Х

где h(-) обозначает максимальное горизонтальное расстояние. Во-вторых, данному определению кривой обслуживания удовлетворяют многие популярные алгоритмы выборки пакетов, в числе которых приоритетная выборка и взвешенная “справедливая” выборка (Weighted Fair Queueing) [5], что позволяет характеризовать реальные узлы коммутации абонентских сетей. Наконец, возможность конкатенации кривых обслуживания позволяет легко переходить от исследования отдельных узлов коммутации к расчету цепочек УК. Для этого совокупная кривая обслуживания вычисляется с помощью мин-плюс свертки отдельных кривых: р(т) = р0 (т) р! (т) ... рд, (т). Перечисленные свойства позволяют использовать кривую поступления при реализации интегрированного или дифференцированного подхода к управлению трафиком. В частности, в первом случае используется кусочно-линейная кривая поступления, определенная на двух интервалах времени и носящая название TSPEC. Несмотря на вышесказанное, необходимо

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

Исправить недостатки детерминированной кривой поступления призвано семейство моделей, называемых стохастическими кривыми поступления. В дополнение к функции а(т), характеризующей поток данных на произвольном интервале т, каждая СКП включает в себя функцию f (с), ограничивающую вероятность превышения некоторого порога, заданного функцией а(т). Одной из наиболее общих СКП является модель с экспоненциальной ограничивающей функцией gSBB [6]. В данной модели кривая поступления а(т)= гт интерпретируется как кривая обслуживания виртуального УК, передающего исследуемый поток. В моменты превышения потоком пропускной способности сервера, в виртуальном УК образуется очередь пакетов. Закон распределения длины этой очереди F(c) отражает вероятностные свойства передаваемого потока, и ограничивающая функция f (с) выбирается таким образом, чтобы выполнялось неравенство f (с) > F(с). Теория стохастического сетевого исчисления позволяет вычислять вероятностные параметры QoS для потока с СКП( f, а), передаваемого через УК с кривой обслуживания р(т). Как правило, практическое значение имеет максимально допустимое значение некоторого параметра QoS (например, задержки), а также вероятность превышения этого значения. Для максимальной задержки справедливо следующее равенство:

d 1 = h(ri + cv,oi, Р(т))>

где с^ы = inf {с: f (с) < pviol}, а pviol - вероятность превышения максимальной задержки. Преимущества модели gSBB проявляются при агрегированном обслуживании группы потоков. Пусть каждый i-й поток характеризуется СКП f, гг). Тогда совокупный поток также

характеризуется некоторой СКП (fagg, rag, и выполняются следующие равенства:

ragg = Г0 + - + rN,

fagg = 1 - (1 - fo) - (1 - fN ).

В случае независимости составляющих потоков мин-плюс свертку во втором равенстве можно заменить на свертку Стилтьеса:

fagg = 1 - (1 - f0) *•••* (1 - fN ).

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

одинаковых микропотоков)

P{A PLO (t) - R(t) > с}< rLO (с),

где pLO (т) = р(т) - (N -1) • гт и TLO (с) = 1 - (1 - f)(N-1)* (с), а f”* (x) обозначает n-кратную свертку Стилтьеса. Тогда максимальная задержка определяется равенством

dГГ = h(rT+ сvol, PLO COX

где с^ы = inf {с: f TLO (с) < pviol}. Аналогично могут быть рассчитаны и другие параметры QoS, например, максимальная длина очереди в УК. Приведенные результаты могут быть использованы для создания системы управления сетевыми ресурсами на основе стохастической модели gSBB. Однако решающим критерием оценки является коэффициент использования, определяемый адекватнгостью модели реальному трафику. Представляется, что использование линейной функции в качестве кривой поступления может существенно ограничивать точность модели, если свойства исследуемого трафика меняются в зависимости от выбранного масштаба времени. Так, например, видео-трафик может демонстрировать различные вероятностные свойства на уровне регулярной группы видео-кадров (group of pictures) и на уровне законченной сцены. Для устранения данного недостатка предлагается ввести дополнительные параметры модели и перейти к кусочно-линейной кривой поступления.

Кусочно-линейная СКП обобщает модель gSBB за счет использования кусочно-линейной кривой поступления а(т) для характеризации трафика:

а1т

(т) =

Ь0 + гт , если 0 < т < т0;

А-1 + гк-\т , если тк-2 <т.

Будем считать, что поток данных A(t) ограничен кусочно-линейной СКП (f, а}, если Vt выполняется следующее неравенство:

P{Q( A, а, t) >с}< e“to,

где Q(A, а, t) - размер очереди виртуального УК с кривой обслуживания а(т) в момент t, а в качестве ограничивающей функции f (с) выбрана экспоненциальная функция. По аналогии с моделью gSBB можно показать, что максимальная задержка для кусочно-линейной СКП определяется равенством

d!r = ^а

при некоторой вероятности превышения pviol . Из определения предлагаемой модели видно, что для любого произвольного потока данных существует бесконечное множество комбинаций параметров Ь0,...,bk-1,r0,...,гк-1, лишь часть из которых дает задержку, удовлетворяющую требованиям к QoS. Из числа подходящих комбинаций параметров наиболее предпочтительными являются те, которые обеспечивают максимальное количество

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

rk-1

Таким образом, искомой комбинацией параметров кусочно-линейной СКП является решение следующей оптимизационной задачи.

Пусть поток данных данных задан ступенчатой функцией Atrace(t), где высота каждой “ступеньки” равна размеру переданного кадра, а расстояние между любыми соседними “ступеньками” одинаково и равно T . Допустим, что N таких потоков передаются по маршруту с совокупной кривой обслуживания р(т) = min(0, R • (т - L)). Для каждого потока задана максимальная допустимая задержка передачи D и вероятность превышения этой задержки риЫ. С учетом этого необходимо найти набор действительных параметров (b0,...,bk-1,r0,...,rk-1,N),

R

который максимизирует целевую функцию U(Ь0,..., , г, •, r_i, N) =----- с соблюдением

Nrk-1

следующих условий.

1. 1 - Fq (с) < e , где Fq (с) - эмпирическая функция распределения длины очереди в

виртуальном УК с кривой обслуживания а(т).

2. dmax = ^а(т) + сую1 , pLO (т)) < D , где

* pLO

* с^1 = inf : (е_Я° )@ rLO (с) < Pviol }

* Tlo (с) =1 - (1 - e_to }N 1)1 (с).

Для упрощения задачи второе условие может быть отброшено, при этом выражение min(0, dmx - D) вычитается из целевой функции в качестве штрафа за превышение допустимой задержки. Благодаря выбору f (с) ограничивающая функция агрегированного потока rLO (с) представляет собой функцию распределения Эрланга и рассчитывается в любой точке стандартными методами. Нахождение вероятностной составляющей сviol для расчета dmax подразумевает решение двух подзадач. Во-первых, вычисление свертки в мин-плюс алгебре требует нахождения минимума суммы в каждой точке: e-to rLo(с) = inf [e-to +rLo(с-р)]

0<р<с

Поскольку сумма двух функций в общем случае является многоэкстремальной, для нахождения минимума используется алгоритм глобальной детерминированной оптимизации на основе липшецовской модели целевой функции [7]. Вторая подзадача заключается в нахождении наименьшего решения неравенства e~yss rLO (с) < риоv Так как левая часть является убывающей функцией, данная задача эквивалентна поиску корня уравнения e~to rLO (с) = рио/ и может быть решена с помощью алгоритма Брента. При известном значении с^о1 целевая функция U может быть получена с помощью простых алгебраических операций.

Так как функция ^ace(t), описывающая исследуемый поток данных, является неизвестной и в значительной степени произвольной, зависящая от нее целевая функция U является многоэкстремальной, а ее производная не может быть определена аналитически. При данных условиях для поиска решения задачи может быть использован алгоритм глобальной стохастической оптимизации. Поскольку вычисление целевой функции требует значительных вычислительных ресурсов, алгоритм поиска должен включать элемент адаптации. Данному требованию удовлетворяет семейство марковских алгоритмов оптимизации под названием имитация затвердевания [8]. В данной работе для выбора новой точки на каждой итерации используется правило Метрополиса-Гастингса, и условием принятия этой точки является выполнение неравенства AU'

R > exp -

\ T

где R - это выборка равномерного распределения на [0,1], а T интерпретируется как “температура”. Полученное таким образом множество точек подчиняется распределению Больцмана. С уменьшением температуры T это распределение концентрируется вокруг решения оптимизационной задачи, что гарантирует сходимость алгоритма при достаточно большом числе итераций. Таким образом, рассмотренный алгоритм позволяет получать комбинацию параметров кусочно-линейной СКП, обеспечивающую соблюдение требований к качеству обслуживания и возможность передачи такого числа потоков, которое является сколь угодно близким к максимальному для данной сети. Заметим, что в частном случае алгоритм может находить параметры модели gSBB и кусочно-линейной детерминированной кривой поступления.

Одним из возможных применений предложенной стохастической модели трафика является система динамического контроля доступа в многоканальных мультисервисных сетях доступа. Рассмотрим более подробно этот вариант применения. Типичная абонентская сеть состоит из узлов коммутации, соединенных в древовидную топологию. К листьям дерева подключены абоненты сети, к корневому узлу - оператор. Каждый УК обслуживает несколько классов трафика с соответствующими очередями, и планировщик пакетов осуществляет строгую или “справедливую” приоритетную выборку из этих очередей. По запросам от оператора УК производит резервирование и освобождение ресурсов, в зависимости от чего меняется кривая обслуживания, предоставляемая соответствующему классу трафика. Для управления резервированием может быть использован традиционный протокол, например, RSVP [ 9]. Сеть доступа осуществляет доставку динамически меняющегося набора мультимедийных потоков (например, видео по запросу). Перед установкой нового соединения поток характеризуется стохастической кривой поступления, параметры которой определяются с помощью глобальной стохастической оптимизации и с учетом пропускной способности сети. На основе полученных параметров формируется дескриптор потока, передаваемый в запросе на резервирование. Для этого можно воспользоваться имеющимся (при использовании RSVP) дескриптором трафика T -SPEC, расширив его полями для доставки параметров стохастической кривой поступления. Получившийся дескриптор предлагается называть ST-SPEC (Stochastic Traffic Specification). Каждый УК на пути передачи потока определяет объем резервируемых ресурсов для полученного ST-SPEC с помощью формулы (1). Новое соединение устанавливается в случае наличия необходимой пропускной способности при условии соблюдения вероятностных гарантий качества обслуживания для всех потоков. На рис. 1 показаны графики зависимостей максимального количество типичных видео-потоков от допустимой задержки для сети с предложенной системой контроля доступа по сравнению с традиционными сетями DiffServ или IntServ. Аналогичная зависимость при меняющейся пропускной способности показана на рис. 2. Из графиков видно, что предложенный метод обеспечивает более эффективное использование сетевых ресурсов по сравнению с традиционными методами, причем наибольший выигрыш за счет статистического мультиплексирования достигается при передаче большого количества потоков.

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

скорость канала передачи, Мбит/с

Рис. 1. Зависимость количества потоков от максимально допустимой задержки для видеофрагментов “Terminator I” и “Silence of the Lambs”

Рис. 2. Зависимость количества потоков от пропускной способности для видео-фрагментов “Terminator I” и “Silence of the Lambs”

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

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

Литература

1. Braden, R. Integrated services in the Internet architecture: an overview / R. Braden, D. Clark, S. Shenker // RFC 1633. - 1994.

2. An architecture for differentiated services / S. Blake, D. Black, M. Carlson et al. // RFC 2475. - 1998.

3. Le Boudec, J.-Y. Network Calculus: A Theory of Deterministic Queuing Systems for the Internet / J.-Y. Le Boudec, P. Thiran. Lecture Notes in Computer Science. - Springer, 2004.

4. Konstantopoulos, T. Optimal flow control schemes that regulate the burstiness of traffic / T. Konstantopoulos, V. Anatharam // IEEE/ACM Transactions on Networking. - 1995. - Vol. 3, no. 4. -Pp. 423-432.

5. Demers, A. Analysis and simulation of a fair queueing algorithm / A. Demers, S. Keshav, S. Shenker // Journal of Internetworking Research and Experience. - 1990.- Pp. 3-26.

6. Liu, Y. A calculus for stochastic QoS analysis / Y. Liu, C.-K. Tham, Y. Jiang // Perform. Eval. -2007. - Vol. 64, no. 6. - Pp. 547-572.

7. Jones, D. Lipschitzian optimization without the Lipschitz constatnt / D. Jones, C. Pertitunen, B. Stuckman // Journal of Optimization Theory and Applications. - 1993. - Vol. 79, no. 1. - Pp. 157_181.

8. Kirkpatrick, S. Optimization by simulated annealing / S. Kirkpatrick, C. D. Gelatt, M. P. Vecchi // Science. - 1983. - Vol. 220, no. 4598. - Pp. 671-680.

9. Resource reservation protocol / R. Braden, L. Zhang, S. Berson, S. Herzog // RFC 2205. - 1997.