Модель управления очередями на маршрутизаторах

УДК 004.021, 519.2

Модель управления очередями на маршрутизаторах

Т. Р. Велиева, А. В. Королькова, Д. С. Кулябов, Б. А. Сантуш

Кафедра систем телекоммуникаций Российский университет дружбы народов ул. Миклухо-Маклая, д.6, Москва, 117198, Россия

Проблемы моделирования активного управления очередью (Active Queue Management, AQM) давно находились в сфере интересов авторов. Одно из направлений работ было связано с динамической моделью управляющего модуля типа Random Early Detection (RED) на основе стохастических дифференциальных уравнений с пуассоновским процессом. Данные уравнения применяются в теории очередей достаточно недавно и не очень хорошо изучены. В качестве недостатков изученного ранее подхода авторы выделяли его частный характер. Было описано взаимодействие модуля RED и протокола TCP Reno, но его расширение на другие варианты протокола TCP и управляющего модуля не представлялось возможным.

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

Ключевые слова: стохастические дифференциальные уравнения, основное кинетическое уравнения, уравнение Фоккера—Планка, активное управление очередями, алгоритм RED.

1. Введение

В работах [1-5] рассмотрена модель управляющего модуля типа Random Early Detection (RED). Представленная модель базируется на теории автоматического управления и является неформализованной. Основная цель данной статьи — построение более реалистичной модели управляющего модуля маршрутизатора. Для этого используется теория стохастических дифференциальных уравнений с винеровским процессом, методы стохастических дифференциальных уравнений с пуассоновским процессом, а также применяется разработанная в [6-8] методика построения одношаговых моделей на базе основного кинетического уравнения.

Приведём основные понятия и принципы работы протокола TCP, а точнее его реализации — TCP Reno, а также рассмотрим алгоритм RED [9], который лежит в основе ряда механизмов предотвращения и контроля перегрузок в очередях маршрутизаторов. Основной акцент сделаем на методике построения одношаго-вых процессов на базе основного кинетического уравнения. Построим стохастическую модель управляющего модуля RED маршрутизатора на основе разработанной в [6-8] методики одношаговых процессов. За основу возьмём ранее рассмотренную детерминистическую модель [1-3]. Её же используем для верификации построенной стохастической модели [10].

2. Адаптивное управление перегрузками 2.1. Механизм управления перегрузками в TCP

В протоколе TCP используется механизм скользящего окна для борьбы с перегрузками. Реализация данного механизма зависит от конкретного стандарта

Статья поступила в редакцию 30 декабря 2013 г.

протокола TCP. Вслед за оригинальной статьёй [1,2] будем рассматривать протокол TCP Reno.

В TCP Reno механизм управления перегрузками состоит из следующих фаз: медленный старт, предотвращение перегрузок, быстрая передача и быстрое восстановление. Динамика изменения размера окна перегрузки (Congestion Window, CWND) зависит от конкретной фазы (см. рис. 1).

Рис. 1. Фазы функционирования TCP

В фазе медленного старта происходит увеличение окна перегрузки каждый раз при получении источником сообщения о доставке пакета (Acknowledge, ACK), т.е. источник увеличивает размер окна перегрузки в зависимости от количества подтверждённых сегментов (Segment Size, SS): cwnd = cwnd +1 для каждого переданного ACK. Первоначальный размер окна перегрузки (Minimum Segment Size, MSS) может принимать значение 1, 2 или 10 сегментов. Приёмник посылает ACK для каждого пакета, однако в реальности можно считать, что подтверждения приходят все вместе по истечении времени двойного оборота (Round-Trip Time, RTT). Таким образом окно перегрузки удваивается по истечении времени двойного оборота.

При достижении определённого размера окна TCP Reno входит в фазу предотвращения перегрузки. В этой фазе окно перегрузки увеличивается на величину 1 /cwnd для каждого подтверждения ACK, что эквивалентно увеличению окна на один пакет за время двойного оборота.

Протокол TCP Reno отслеживает два варианта потери пакетов:

- Тройное дублирование подтверждения (Triple Duplicate ACK, TD). Пусть n-й пакет не доставлен, а последующие пакеты (п + 1, п + 2 и т.д.) доставлены. Для каждого пакета, доставленного в нарушении очерёдности (для п+1, п+2 и т.д.) получатель отсылает сообщение ACK для последнего недоставленного (n-го) пакета. При получении трёх таких пакетов источник перепосылает п-й пакет. Кроме того, размер окна уменьшается в 2 раза cwnd ^ cwnd/2.

— Тайм-аут (Timeout, TO). При отправке пакета запускается таймер тайм-аута. При получении каждого подтверждения таймер перезапускается. Окно при этом устанавливается в начальное значение окна перегрузки. Первый потерянный пакет перепосылается. Протокол переходит в фазу медленного старта.

Общий алгоритм управления перегрузкой относится к типу AIMD (Additive Increase, Multiplicative Decrease) — аддитивное увеличение размера окна и мультипликативное его уменьшение.

2.2. Механизм адаптивного управления перегрузками RED

Для улучшения рабочих характеристик канала необходимо оптимизировать управление очередями на маршрутизаторах. Одним из возможных подходов является применение алгоритма случайного раннего обнаружения (Random Early Detection, RED) [9,11].

Алгоритм RED использует взвешенное значение длины очереди в качестве фактора, определяющего вероятность отбрасывания пакета. По мере роста средней длины очереди растёт вероятность отбрасывания пакетов (см. (1)). Для управления функцией сброса используются два пороговых значения средневзвешенной длины очереди, влияющих на функционирование алгоритма (рис. 2):

( 0,

0 < Q <Qn

Р(Q) Н У Qmin Pmax, Qmin < Q < Qm

Qmax Qm 1,

Q > Qn

Здесь p((Q)) — очереди, Qmin очереди, max

- функция сброса пакета, ( — средневзвешенное значение длины и (шах — пороговые значения средневзвешенного значения длины

— максимальный уровень сброса пакетов.

Ртах

Рис. 2. Функция сброса RED

3. Стохастические дифференциальные уравнения

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

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

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

3.1. Обозначения и соглашения

1. Будем использовать нотацию абстрактных индексов. В данной нотации тензор как целостный объект обозначается просто индексом (например, жг), компоненты обозначаются подчёркнутым индексом (например, ж-).

2. Запятой в индексе обозначается частная производная по соответствующей координате.

3.2. СДУ с винеровским процессом

Обычно под стохастическим дифференциальным уравнением понимают уравнение в форме Ланжевена, имеющего вид:

dx = adt + bdW,

где х := x(t) £ R — некоторый случайный процесс, а := а(х, t) £ R, b := b(x, t) £ R, а W := W(t) £ R — винеровских процесс.

В многомерном случае уравнение (2) принимает вид:

dx1 = aidt + b^dWа,

где хг := xl(t) £ R" — n-мерный случайный процесс, аг := аг(хк,t) £ R" ^ R™, К := Ьга(хк,t) £ R" х Rs, а Wа := Wa(t) £ Rs — s-мерный винеровский процесс.

При интегрировании уравнения Ланжевена с непостоянными коэффициентами (Ьга = const) у винеровского процесса возникает неопределённость учёта скачка, задаваемого винеровским процессом, при интегрировании уравнения. Данная проблема снимается с введением определённой интерпретации стохастического уравнения (чаще всего используют интерпретации Ито или Стратановича). Мы будем использовать интерпретацию Ито.

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

Пусть f := f (xk,t) — функция от случайного процесса xk(t), f £ C2. Тогда формула дифференциала будет выглядеть следующим образом [12]:

df =

dtf + a'fr + - bibiaftij

dt + blf,i dWa,

где f := f (xk,t), ai := ai(xk,t), Ьa := bi(xk,t) и d^a := d^a(t).

Стохастическому дифференциальному уравнению (3) соответствует уравнение Фоккера-Планка. Оно представляет собой уравнение в частных производных и описывает изменение плотности вероятности во времени:

dtp = -дг[Агр]+ didj [Вгэ р],

где р := р(хк,t), Аг := Al(xk,t), ВгJ := ВгJ(xk,t), р имеет смысл плотности распределения случайной величины хк, Аг — вектор сноса, ВгJ — вектор диффузии.

3.3. СДУ с пуассоновским процессом

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

dx = adt + bdN,

где х := х(1) £ R — некоторый случайный процесс, а := а(х, £) £ R, Ь := Ъ(х, £) £ R, а N := N(£) £ R — пуассоновский процесс.

Для многомерного случая уравнение (6) принимает вид:

dxi = а^ + Ь^Ма,

где хг := хг(1) £ Rn — п-мерный случайный процесс, аг := аг(хк,1) £ Rn, Ьга := Ъга(хк,1) £ Rn х Rs, а Ма := Nа(1) £ Rs — 8-мерный пуассоновский процесс.

Также запишем правило Ито для пуассоновского процесса:

d f = [dtf + aif,i]dt + [f

4. Построение модели модуля RED

Рассмотрим модель, содержащую два основных элемента: источник и получатель. В качестве получателя представлена очередь. Источник отправляет пакеты, получатель их обрабатывает и посылает подтверждение о принятии пакета.

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

4.1. Модельные допущения

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

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

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

4.2. Переход к непрерывной модели

Для применения метода стохастизации одношаговых процессов (см. [6-8]) мы должны перейти к модели с непрерывным временем. Вектор состояния получившейся системы также должен быть непрерывным. В качестве вектора состояний будем использовать ( Ш, ()т, где Ш := Ш(£) — окно, ( := ((£) — размер очереди, ( :=((£) — экспоненциально взвешенное скользящее среднее длины очереди. Строчные буквы мы зарезервируем для плотностей распределения соответствующих величин.

4.2.1. Уравнение для окна перегрузки

Используя результаты раздела 2.1, формализуем поведение нашей модели. Изменение размера TCP-окна связано с элементарным событием, которому соответствует либо приход одного подтверждения, либо приход всех подтверждений. Примем за элементарное событие приход всех подтверждений, который происходит за время двойного оборота (RTT).

В фазе медленного старта (Slow Start) размер окна увеличивается при каждом приходе подтверждения ACK:

W (tfJK + A tACK) = W (tiCK) + 1.

Перепишем теперь (9) относительно времени двойного оборота Т (за время двойного оборота приходят все подтверждения на отправленное окно):

—

Аг аь '

Считая, что А! = Т, получим

ё — —

— = —,

ё1п — =М, 1пШ = — = ехр Т.

Таким образом, окно растёт экспоненциально, как и должно быть в стадии медленного старта в соответствии с описанием ТСР.

Аналогично рассмотрим фазу избежания перегрузок. При каждом приходе подтверждения АСК размер окна увеличивается:

—

Перепишем для времени двойного оборота (т.е. при получении всех подтверждений):

—(1п + А€) = — (1п) + ^п);

— п)

—(1п + А€) - — (1п) 1

Аг Аг'

Считая, что А! = Т, имеем

— = -

ё = г

ё — = — = 1. (16) В результате получим линейный рост окна, как и описано в спецификации ТСР.

4.2.2. Уравнение для средней длины очереди

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

Пусть — параметр скользящего среднего [3,9]. Исходя из формулы для экспоненциально взвешенного скользящего среднего, запишем:

(Ш п + А = % т п) -да.)).

Положим ё = А£. Придадим ё смысл времени нахождения одного пакета в очереди. Считая, что интенсивность обслуживания есть С, можно записать ё = -—. Тогда уравнение для экспоненциально взвешенной скользящей средней длины

очереди принимает вид:

= ^

4.3. Кинетические уравнения

4.3.1. Уравнение для окна перегрузки

Рассмотрим поведение пакетов в системе. Количество пакетов будет задаваться окном TCP. Возможны два процесса: с интенсивностью k} пакеты из бесконечного резервуара поступают в систему, а с интенсивностью k} пакеты покидают систему. Запишем это в виде кинетических уравнений:

к1

0 %W,

ki , (19) W % 0.

Первое соотношение описывает появление пакетов, второе — вывод пакетов из системы.

Распишем получение стохастических дифференциальных уравнений. Из кинетических уравнений (19) получаем операторы состояний. Число пакетов до взаимодействия:

N? =

число пакетов после взаимодействия:

М* =

Оператор перехода г* = М* — М* имеет вид:

г} =

Общая форма для удельной вероятности перехода в прямом направлении задаётся формулой:

% !

= к*П

Учитывая (23) и (20), получим:

= к\, = к12~Щ.

Используя (24) и (22), получим коэффициенты уравнения Фоккера-Планка:

А{ = Б+г * = к! —

* 1 2 ,

Вij = S+r airaj = k} + k}W.

Придадим коэффициентам конкретные значения. Из (15) получим, что k

W

„™nw„OD„a ™ 1,1 -__

2 di

щ. Поскольку мы учитываем потерю пакетов только по тройному дублированию подтверждения, то положим к^ = ^^Т, где — винеровский процесс, аналогичный введённому в [2].

Запишем уравнение Фоккера-Планка:

дю

д_ д—

1 — ёЛТ\

1 д2 + х-

2дШ2

1 — ёЛТ\

где ю := ю(1) — плотность распределения случайного процесса — (¿). Соответствующее уравнение Ланжевена имеет вид:

1 , — „г 1 — ёМ .

где ё У1 — винеровский процесс, соответствующий случайному процессу — (Ь).

4.3.2. Уравнение для мгновенной длины очереди

Аналогично описывается поведение очереди. Возможны два процесса: с интенсивностью А;2 пакеты из бесконечного резервуара поступают в очередь, а с интенсивностью - пакеты покидают систему (мы рассматриваем это как поступление в систему отрицательного количества пакетов). Запишем это в виде кинетических уравнений:

к2

0 ^ (28) (.

0 к22

Из (28) получим операторы состояний:

МГ =

(0) • м

Оператор перехода г" = МГ - МГ имеет вид:

Г1 =(1), Г2 =(1).

Удельная вероятность перехода в прямом направлении:

= , = А2 .

Используя (31) и (30), получим коэффициенты уравнения Фоккера-Планка:

А* = Б+гг = к\ + к%,

В* = Б+г ггг = к2 + к2.

Придадим коэффициентам конкретные значения. Интенсивность входящего потока есть количество пакетов за единицу времени: к2 = . Интенсивность обслуживания С, соответственно, имеет вид к| = -С.

Запишем уравнение Фоккера-Планка:

д д

д_ д(

(Т-с >

1 д2

+

2д(2

(— -с).

где д := — плотность распределения случайного процесса <3- Тогда уравнение Ланжевена имеет вид:

<1<э =

где (IV2 — винеровский процесс, соответствующий случайному процессу <5-

5. Результаты

Учитывая полученные уравнения (27), (34) и (18), запишем результирующую систему:

< dQ=

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

В этой работе мы не проводим тщательный анализ влияния винеровского процесса на полученную модель, однако для иллюстрации такого влияния провели численный эксперимент. На приведённых графиках демонстрируется поведение величин, описывающих систему, а именно изменение окна перегрузки без учёта (рис. 3) и с учётом (рис. 4) винеровского процесса. Аналогично приводятся графики для изменения мгновенной (рис. 5 и 6) и средней (рис. 7 и 8) длины очереди. В конце приводится фазовый портрет (рис. 9 и 10) без учёта и с учётом винеровского процесса соответственно.

Рис. 3. Окно перегрузки без учёта винеровского процесса

Рис. 4. Окно перегрузки с учётом винеровского процесса

Рис. 5. Мгновенная длина очереди Рис. 6. Мгновенная длина очереди с без учёта винеровского процесса учётом винеровского процесса

Детерминированная модвгъ НЕР

ЯП)

20 40 60 80

4 — время

1 ао

5

т 60

&

1 40

0 1

СУ 20

Стохастическая модель ЯЕР

10 20

40 50

Рис. 7. Средняя длина очереди без Рис. 8. Средняя длина очереди с учёта винеровского процесса учётом винеровского процесса

Рис. 9. Фазовый портрет без учёта Рис. 10. Фазовый портрет с учётом винеровского процесса винеровского процесса

6. Заключение

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

Литература

1. Misra V., Gong W.-B., Towsley D. Stochastic Differential Equation Modeling and Analysis of TCP-windowsize Behavior // Proceedings of IFIP WG 7.3 Performance. — 1999. — Vol. 99. — http://dna-pubs.cs.columbia.edu/citation/ paperfile/24/Misra99-TCP-Stochastic.pdf.

2. Misra V., Gong W.-B., Towsley D. Fluid-Based Analysis of a Network of AQM Routers Supporting TCP Flows with an Application to RED // ACM SIGCOMM Computer Communication Review. — 2000. — Vol. 30, No 4. — Pp. 151-160.

3. Королькова А. В., Кулябов Д. С. Математическая модель динамики поведения параметров систем типа RED // Вестник РУДН. Серия «Математика. Информатика. Физика». — 2010. — № 1. — С. 68-76. [Korolkova A. V., Kulyabov D.S. Mathematical Model of the Dynamic Behavior of RED-Like System Parameters // Bulletin of Peoples' Friendship University of Russia. Series "Mathematics. Information Sciences. Physics". — 2010. — No. 1. — Pp. 68-76. — (in russian). ]

4. Королькова А, Черноиванов А. Модификация модели процесса передачи с регулированием алгоритмом типа RED интенсивности потока для случая TCP-NewReno трафика // International Workshop "Distributed Computer and Communication Networks (DCCN-2010)". — М.: R&D Company "Information and Networking Technologies", 2010. — С. 262-267. [Korolkova A., Tchernoivanov A. Modification of Transmission Process with Intensity Regulated by RED for the TCP-NewReno Traffic // International Workshop "Distributed Computer and Communication Networks (DCCN-2010)". — Moscow: R&D Company "Information and Networking Technologies", 2010. — P. 262-267. — (in russian). ]

5. Королькова А. В, Черноиванов А. И. Использование СДУ для моделирования поведения TCP-трафика при взаимодействии с узлом, работающим по алгоритму RED. Определение области возникновения автоколебаний на примере алгоритмов RED, ARED, RARED, POWARED // «Математика. Компьютер. Образование». C6. научных трудов. / под ред. Г. Ю. Ризниченко. — М.-Ижевск: НИЦ «Регулярная и хаотическая динамика», 2010. — Т. 1. — С. 270-282. — http://www.mce.su/rus/books/book72814/. [Korolkova A.V., Tchernoivanov A. I. The TCP-Traffic and RED-Like Node Interaction Model Based on Sthchastic Differential Equations. Region of Self-Oscillations for RED, ARED, RARED, POWARED Algorithms // Proceedings of "Mathematics. Computer. Education" / ed. by G. Yu. Riznichenko. — Moscow-Izhevsk: Research Centre "Regular and Chaotic Dynamics", 2010. — Vol. 1. — Pp. 270-282. — (in russian). ]

6. Демидова А. В. Метод стохастизации математических моделей на примере системы «хищник-жертва» // Научная сессия НИЯУ МИФИ-2013. — 2013. — С. 127. [Demidova A. V. The Method of Stochastization of Mathematical Models for the Example of the "Predator-Prey" // A Scientific Session NRNU MEPHI-2013. — 2013. — P. 127. — (in russian). ]

7. The Method of Stochastization of One-Step Processes / A. V. Demidova, A. V. Korolkova, D. S. Kulyabov, L. A. Sevastianov // Mathematical Modeling and Computational Physics. — Dubna: JINR, 2013. — P. 67. — http://mmcp2013.jinr.

ru/prog.php.

8. Влияние стохастизации на одношаговые модели / А. В. Демидова, М. Н. Геворкян, А. Д. Егоров и др. // Вестник РУДН. Серия «Математика. Информатика. Физика». — 2014. — № 1. — С. 71-85. [Influence of Stochastization on One-Step Models / A. V. Demidova, M.N. Gevorkyan, A.D. Egorov et al. // Bulletin of Peoples' Friendship University of Russia. Series "Mathematics. Information Sciences. Physics". — 2014. — No. 1. — Pp. 71-85. — (in russian). ]

9. Floyd S., Jacobson V. Random Early Detection Gateways for Congestion Avoidance // IEEE/ACM Transactions on Networking. — 1993. — Vol. 1, No 4. — Pp. 397-413.

10. Кулябов Д. С., Королькова А. В., Зарядов И. С. Обзор подходов к моделированию модуля управления трафиком // T-Comm — Телекоммуникации и транспорт. — 2012. — № 7. — С. 122-125. [Kulyabov D. S., Korolkova A. V., Zaryadov I. S. Traffic Management Module Modeling Approaches Review // T-Comm. — 2012. — No 7. — Pp. 122-125. — (in russian). ]

11. Королькова А. В., Кулябов Д. С., Черноиванов А. И. К вопросу о классификации алгоритмов RED // Вестник РУДН. Серия «Математика. Информатика. Физика». — 2009. — № 3. — С. 34-46. [Korolkova A.V., Kulyabov D.S., Tchernoivanov A. I. On the Classification of RED Algorithms // Bulletin of Peoples' Friendship University of Russia. Series "Mathematics. Information Sciences. Physics". — 2009. — No. 3. — Pp. 34-46. — (in russian). ]

12. 0ksendal B. K. Stochastic Differential Equations: An Introduction with Applications. — Berlin: Springer, 2003.

UDC 004.021, 519.2

Model Queue Management on Routers

T.R. Velieva, A. V. Korolkova, D. S. Kulyabov, B.A. Dos Santos

Telecommunication Systems Department Peoples' Friendship University of Russia 6, Miklukho-Maklaya str., Moscow, Russia, 117198

Problems of modeling of active queue management (AQM) systems have been for a long time in the sphere of interests of authors. One of the areas of work was associated with a dynamic model of the control module Random Early Detection (RED) based on Poisson process driven stochastic differential equations. These equations are used in queuing theory quite recently, and not very well understood. As disadvantages of this approach the authors underlined its non-generic character. We describe the interaction between module RED and protocol TCP Reno. But its extension to other variants of the TCP protocol and the control module is not possible.

Our group studied common approaches to modeling of such phenomena. As a result a method for randomization of one-step processes, allowing to obtain new models in a universal manner was developed. In this paper, the authors use this technique to model previously investigated RED module and protocol TCP Reno to demonstrate its applicability to this kind of problems. As a result an extended model of control module type RED for traffic type TCP Reno, was created. It contains previously studied model as a special case.

Key words and phrases: stochastic differential equations, master equation, Fokker-Planck equation, AQM, RED.