Самоподобный (self-similar) входной поток в системах M/M/1 и G/M/1
Ильницкий С.В.
Ссылка на оригинал.Последние исследования различных типов сетевого трафика убедительно доказывают, что сетевой трафик является самоподобным (self-similar) или фрактальным (fractal) по своей природе, т.е. в нем присутствуют так называемые вспышки или пачки (burst) пакетов, наблюдаемые в различных временных интервалах (от милисекунд до минут или даже часов), а также корреляция между пакетами. Из этого следует, что широко используемые в настоящее время методы моделирования и расчета сетевых систем, основанные на традиционных предположениях, не дают адекватной картины происходящего в сети. Было бы логично предполагать, что рабочие параметры сетевых систем с самоподобным трафиком отличаются по своим значениям от систем с простейшим (Пуассоновским) потоком на входе.
В данной публикации проведено cравнение рабочих параметров сетевых систем, представленных моделью M/M/1 (с Пуассоновским входным потоком (Poison input flow)) и моделью G/M/1 с самоподобным входным потоком, а также анализ средней длины очереди, среднего времени ожидания в очереди и определение максимальной длины очереди в случае различных входных потоков.
КОМПЬЮТЕРНЫЕ СЕТИ, ТРАФИК, МОДЕЛИРОВАНИЕ, МЕТОДЫ РАСЧЕТА, САМОПОДОБНОСТЬ.
Keywords: network traffic, self-similar traffic, traffic analysis.
Последние исследования различных типов сетевого трафика убедительно доказывают, что сетевой трафик является самоподобным (self-similar) или фрактальным (fractal) по своей природе, т.е. в нем присутствуют так называемые вспышки или пачки (burst) пакетов, наблюдаемые в различных временных интервалах (от милисекунд до минут или даже часов). Из этого следует, что широко используемые в настоящее время методы моделирования и расчета сетевых систем, основанные на традиционных предположениях, не дают адекватной картины происходящего в сети. Было бы логично предполагать, что рабочие параметры систем с самоподобным трафиком отличаются по своим значениям от традиционных, рассчитанных для систем с простейшим (Пуассоновским) потоком на входе. Сравнение удобно провести на моделях систем массового обслуживания. Для этого проанализируем поведение моделей M/M/1 (с Пуассоновским входным потоком (Poison input flow)) и G/M/1 с самоподобным входным потоком, а также сделаем анализ средней длины очереди, среднего времени ожидания в очереди и определение максимальной длины очереди в случае различных входных потоков.
2. Сравнение рабочих параметров для систем М/М/1 и С/М/1 с самоподобным входным потоком
Проанализируем основные характеристики (среднее время ожидания, среднюю длину очереди и максимальную длину очереди) для системы G/M/1. В качестве входного потока будем использовать поток, сгенерированный с помощью модели, описанной в работе [7]. Время ожидания измеряется встроенной подсистемой GPSS в зависимости от загрузки системы:
где
X - средняя интенсивность входного потока, р - средняя интенсивность обслуживания.
Время обслуживания транзакции имеет экспоненциальный закон распределения.
В начале выполнения программы происходит инициализация параметров. Количество ON-OFF источников выбрано согласно топологии исследуемой сети [7] и подразумевает следующее: каждая рабочая станция в сети является источником запросов к серверу. Общее количество рабочих станций в сети равно 350. В 1-м блоке программы происходит прозрачное включение 350 ON-OFF источников.
Во 2-м блоке, в зависимости от количества включенных источников, генерируется поток транзакций, который затем будет передан в 3-й блок программы.
В 3-м блоке транзакции поступают на обслуживающий прибор GPSS, названный PRIB, и задерживаются в нем на время, которое распределено по экспоненциальному закону и имеющим среднее значение TM2. Когда устройство обслуживания занято, транзакции поступают в очередь и ожидают там до тех пор, пока обслуживающее устройство не будет свободным. Длина очереди в нашем случае бесконечна, а приоритеты не используются.
В 4-м блоке программы определяется время моделирования.
Для получения точных результатов необходимо повторять эксперименты много раз, изменяя только случайные последовательности, инициализируя генератор и вычисляя средние значения результатов. Исходные коды программы находятся у автора данной статьи.
Для сравнения результатов с классической системой с экспоненциальным временем распределения необходимо провести дополнительные расчеты, сохраняя без изменения время моделирования и другие параметры. Входной поток имеет экспоненциальное время распределения со средним значением а процесс обслуживания происходит в том же самом блоке программы, что и в случае самоподобного входного трафика в предыдущем эксперименте.
3. Оценка результатов экспериментов
Результаты моделирования систем G/M/1 и M/M/1показаны в
Таблице 1 и 2.
р |
Средняя длина очереди |
Максимальная длина очереди |
Ср. время ожидания в очереди |
0.179579 |
1.574 |
42 |
12.87 |
0.236548 |
2.587 |
39 |
14.275 |
0.279621 |
3.165 |
41 |
15.904 |
0.348356 |
4.492 |
42 |
19.443 |
0.410023 |
4.949 |
44 |
21.747 |
0.447126 |
5.094 |
42 |
25.465 |
0.490738 |
5.792 |
42 |
26.549 |
0.53174 |
6.184 |
47 |
27.107 |
0.581762 |
6.845 |
48 |
30.314 |
0.609915 |
7.112 |
50 |
31.599 |
0.621431 |
7.078 |
54 |
32.919 |
0.649165 |
7.534 |
59 |
34.829 |
0.681342 |
7.993 |
56 |
37.032 |
0.72452 |
8.892 |
76 |
40.846 |
0.749358 |
9.438 |
93 |
43.309 |
0.781564 |
10.517 |
102 |
49.282 |
0.79131 |
10.752 |
107 |
52.456 |
0.79782 |
11.849 |
115 |
55.693 |
0.816593 |
12.994 |
114 |
62.856 |
0.84527 |
14.528 |
145 |
71.672 |
0.865794 |
18.319 |
171 |
89.568 |
0.881237 |
24.966 |
284 |
121.432 |
0.897672 |
41.476 |
422 |
212.654 |
0.931744 |
56.362 |
678 |
345.615 |
0.961234 |
212.26 |
1054 |
1023.312 |
р |
Средняя длина очереди |
Максимальная длина очереди |
Ср. время ожидания в очереди |
0.975936 |
524.23 |
2952 |
2563.912 |
0.983675 |
914.96 |
3447 |
5486.7817 |
0.998179 |
2107.48 |
4589 |
9901.5832 |
Таблица 1. Основные параметры для системы G/M/1.
Р |
Средняя длина очереди |
Максимальная длина очереди |
Ср. время ожидания в очереди |
0.1682 |
0. 034 |
3 |
0. 154 |
0.227093 |
0. 063 |
5 |
0.263 |
0.267587 |
0. 953 |
7 |
0.424 |
0.323471 |
0. 142 |
8 |
0.712 |
0.38683 |
0.241 |
8 |
1.203 |
0.415612 |
0.301 |
11 |
1.487 |
0.455437 |
0.324 |
12 |
1. 816 |
0.50124 |
0.398 |
13 |
2 . 356 |
0.525853 |
0.586 |
14 |
3. 132 |
0.583696 |
0. 614 |
18 |
3. 812 |
0.60468 |
0.763 |
18 |
4.221 |
0.626788 |
1. 002 |
19 |
5.487 |
0.646423 |
1. 119 |
20 |
6. 021 |
0.663456 |
1.382 |
23 |
6. 925 |
0.701446 |
1.409 |
22 |
7 . 826 |
0.704313 |
1. 634 |
26 |
8 . 993 |
0.724567 |
1. 893 |
29 |
9. 978 |
0.746378 |
1. 065 |
31 |
10.52 |
0.774567 |
1.255 |
35 |
12.260 |
0.783567 |
2 . 043 |
39 |
14.016 |
0.816843 |
2 . 472 |
45 |
15.402 |
0.840467 |
2 . 989 |
50 |
19.984 |
0.847783 |
3. 172 |
52 |
24.187 |
0.868433 |
5. 154 |
60 |
31.116 |
0.89286 |
6. 135 |
71 |
35.136 |
0.922677 |
7 . 747 |
72 |
42.858 |
0.955677 |
16.634 |
104 |
97.865 |
0.983567 |
29. 984 |
247 |
202.537 |
0.99853 |
125. 183 |
496 |
595.173 |
Таблица 2. Основные параметры для системы M/M/1.
Покажем отличия в значениях рабочих параметров для случая классического входного потока (система M/M/1) и самоподобного входного потока (система G/M/1). Для этого построим графики средних значений длины очереди, максимальной длины очереди и среднего времени ожидания в очереди в зависимости от загруженности системы.
Используя данные из Таблицы 1 и 2, строим графики средней длины очереди в зависимости от загруженности системы.
Рисунок 1. Средняя длина очереди для систем G/M/1 и M/M/1.
Так как результаты для различных р отличаются в тысячи раз, разделим график на две части: диапазоны значений 0.2< р<0.8 и 0.8 < р< 1.0, затем проанализируем их независимо.
На Рисунке 2 изображен график для диапазона 0.2< р<0.8. Как видно из графика, в случае самоподобного входного трафика средняя длина очереди намного больше по сравнению с простешим (Пуассоновским) входным трафиком. Разница составляет 2-6 раза даже в случае, когда система не была нагружена. Очередь растет даже тогда, когда интенсивность на входе намного меньше интенсивности обслуживания. Это может быть объяснено особой структурой пакетов самоподобного трафика (в нем присутствуют так называемые вспышки или пачки (burst) пакетов).
Рисунок 2. Средняя длина очереди для О/М/1 и М/М/1 систем при 0.2< р<0.8.
На Рисунке 3 показан график для случая 0.8< р<1.0. нагрузка на систему возросла и разница между значениями средней длины очереди для различных законов распределения тоже возросла. В начале диапазона (р=0.8) разница между двумя графиками составляет 6 раз. Когда р=0.95среднее значение длины очереди для систем О/М/1 в 17 раз больше, чем для систем М/М/1.
Рисунок 3. Средняя длина очереди для С/М/1 и М/М/1 систем при 0.8 < р<1.0
С помощью полученных в ходе экспериментов данных [7] (Таблица 1 и 2) строим график максимальной длины очереди в зависимости от нагрузки на систему (Рисунок 4).
Рисунок 4. Максимальная длина очереди для систем О/М/1 и М/М/1
Так как результаты для различных р отличаются в тысячи раз, разделим график на две части: диапазоны значений 0.2< р<0.8 и 0.8 < р< 1.0, затем проанализируем их независимо.
На Рисунке 5 изображен график для диапазона 0.2< р<0.8. Как видно из графика, в случае самоподобного входного трафика максимальная длина очереди больше по сравнению со случаем простейшего входного потока. В начале диапазона (р=0.2) разница между двумя графиками составляет 5 раз. Когда р возрастает до значения 0.5, разница сокращается до 2-х раз, но в дальнейшем, когда р продолжает расти, разница между значениями для О/М/1 и М/М/1 системами вновь возрастает и в момент, когда р=0.8 отличие составляет уже 3 раза.
Данные различия на графике можно объяснить характерной для самоподобного трафика структурой пакетов. (Структура пакетов в этом случае характеризуется наличием множественных коллизий даже если средняя интенсивность входного трафика меньше интенсивности обслуживания.) В результате получаем, что для малых значений р разница максимальной длины очереди для систем О/М/1 и М/М/1 выражена больше, чем для средних значений р.
Рисунок 5. Максимальная длина очереди для О/М/1 и М/М/1 систем, 0.2< р<0.8
С помощью полученных в ходе экспериментов данных [7] (Таблица 1 и 2) строим график максимальной длины очереди в зависимости от нагрузки на систему (Рисунок 6).
По аналогии с предыдущими разделами, так как результаты для различных р отличаются в тысячи раз, разделим график на две части: диапазоны значений 0.2< р<0.8 и 0.8 < р< 1.0, затем проанализируем их независимо.
Рисунок 6. Среднее время ожидания для систем О/М/1 и М/М/1.
На Рисунке 7 показан график для среднего времени ожидания в случае 0.2<р<0.8. Как можно видеть на графике, в случае самоподобного трафика среднее время ожидания в очереди больше по сравнению со случаем простейшего входного потока. В начале диапазона (р=0.2) разница между графиками составляет 15 раз. Когда р возрастает до значения 0.5, разница снижается до 10 раз. В дальнейшем, когда р продолжает расти, разница между системами О/М/1 и М/М/1 уменьшается и в момент, когда р=0.8 составляет всего 4.5 раза.
Объяснение этому факту следует искать в особенностях структуры пакетов самоподобного трафика. Даже в случае, когда интенсивность входного потока меньше интенсивности обслуживания, имеются различия между системами с самоподобным и простейшим входными потоками.
Рисунок 7. Среднее время ожидания для систем О/М/1 и М/М/1, 0.2< р<0.8
После анализа рабочих параметров систем О/М/1 и М/М/1 можем сделать следующие выводы:
• Все основные параметры (средняя длина очереди, максимальная длина очереди а также среднее время ожидания в очереди) для системы О/М/1 имеют большие значения по сравнению с системой М/М/1 с аналогичными параметрами (за исключением типа входного потока)
• Когда система не нагружена (р=0...3), разница составляет 2 раза для средней длины очереди, 5 раз для максимальной длины очереди и 15 раз для среднего времени ожидания в очереди. На это сильно влияет структура пакетов самоподобного трафика, имеющего так называемые вспышки или пачки (burst) пакетов.
• Когда система работает в нормальном режиме (р=0.4...0.7), разница между значениями рабочих параметров для систем G/M/1 и М/М/1 не столь велика, как для небольших значений р, но с повышением нагрузки на систему тоже постепенно возрастает: 6 раз для средней длины очереди, 2-3 раза для максимальной длины очереди и 3 раза для среднего времени ожидания.
• Когда система работает в нагруженном режиме (р=0.75.. .0.99), разница в значениях рабочих параметров между двумя системами продолжает возрастать и когда р превышает значение 0.86, это возрастание становится очень стремительным: 17 раз для средней длины очереди, 10 раз для максимальной длины очереди. (Стремительное возрастание разницы значений параметров для максимальной длины очереди начинается даже раньше, когда р=0.75.)
• Длина очереди в нашем эксперименте была бесконечной и не были использованы приоритеты, но в реальной жизни длина очереди имеет определенное значение и поэтому влияние самоподобного входного потока на производительность системы еще более заметно и иногда приводит к нестабильности в ее работе или большим задержкам.
• Влияние эффекта самопободности необходимо учитывать при проектировании и анализе систем.
[1] Leland W.E., Taqqu M.S., Willinger W. and Wilson D.V. On The Self-Similar Nature Of Ethernet Traffic.- Proc. АСМ SIGCOMM'93, San-Fransisco, 1993, p 183-193.
[2] Walter Willinger, Murad S. Taqqu, Robert Sherman and Daniel V. Wilson, "Self-Similarity Through High-Variability: Statistical Analysis of Ethernet LAN Traffic at the Source Level". IEEE/ACM Transactions on Networking, Vol. 5, No. 1, 1997
[3] Adrian Popescu, Traffic Self-Similarity, Department of Telecommunications and Signal Processing, University of Karlskrona/Ronneby
[4] Paxson V., Fast Approximation of Self-Similar Network Traffic, preprint Lawrence Berkely Laboratory and EECS Disvision, University of California, 1996
[5] P.Ulanovs, E.Petersons, Modeling of Self-similar traffic in high-performance Computer Networks, RTU, report collection for 42 International Conference, Riga, 2001
[6] Цыбаков Б.С. Модель телетрафика на основе самоподобного случайного процесса. "Радиотехника", 1999 г., № 5.
[7] S.Ilnickis, Research of the network server in self-similar traffic environment, RTU, Riga, 2003.
Об авторе:
Ильницкий Сергей Васильевич, докторант Рижского Технического Университета, факультет „Транспортная Телематика".
Служ. адрес: Латвия, Рига, Ломоносова l, Рижский Технический Университет. Тел: (+371) 7089984,
Дом. адрес: Латвия, Олайне, LV-2114, ул. Земгалес, 34-48. Моб. тел: (+371) 9441875 E-mail: Sergejs.Ilnickis@dati.lv