ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2013 Управление, вычислительная техника и информатика № 3(24)
УДК 681.324
В.В. Кокшенев, П.А. Михеев, С.П. Сущенко АНАЛИЗ СЕЛЕКТИВНОГО РЕЖИМА ОТКАЗА ТРАНСПОРТНОГО ПРОТОКОЛА В НАГРУЖЕННОМ ТРАКТЕ ПЕРЕДАЧЕ ДАННЫХ
Предложена модель асинхронной процедуры управления виртуальным соединением транспортного протокола в режиме селективного отказа в виде двумерной марковской цепи с дискретным временем, учитывающая влияние протокольных параметров размера окна и длительности тайм-аута ожидания сквозных квитанций, вероятности искажения пакетов в отдельных звеньях тракта передачи данных и длин очередей в транзитных узлах от «внешних» потоков на пропускную способность виртуального соединения. Проведен анализ зависимости пропускной способности управляющей процедуры от протокольных параметров, уровня ошибок в каналах связи, длины тракта передачи данных, распределения размеров очередей в транзитных узлах.
Ключевые слова: транспортный протокол, нагруженный тракт передачи данных, математическая модель, цепь Маркова, быстродействие виртуального соединения, размер окна, длительность сквозного тайм-аута.
В современных сетях, осуществляющих передачу компьютерных данных и мультимедийного абонентского трафика через единую сетевую инфраструктуру, существенно повышаются требования к наличию доступной полосы пропускания и, как следствие, к повышению эффективности ее использования. Важнейшей операционной характеристикой виртуального соединения, управляемого транспортным протоколом компьютерной сети, является его пропускная способность. Данный показатель определяется не только скоростью и достоверностью передачи данных в информационных каналах транспортного соединения, но и интенсивностью внешних по отношению к данному соединению потоков, имеющих с ним часть общего маршрута. Основным индикатором «внешней» нагрузки на тракт, в котором проложено исследуемое виртуальное соединение, являются размеры очередей перед единицами потока (сегментами, инкапсулируемыми на сетевом уровне в пакеты, которые на канальном уровне инкапсулируются в кадры) рассматриваемого соединения в транзитных узлах. Очевидно, что мониторинг такого индикатора позволяет оценить распределение длин очередей в транзитных узлах от внешних по отношению к анализируемому соединению сетевых потоков и использовать при расчете операционных характеристик соединения и выборе протокольных параметров на время сеанса связи между заданной парой абонентов. Известные модели асинхронных управляющих процедур отдельного звена передачи данных и транспортного протокола [1-7] не позволяют учитывать нагрузку на разделяемые сетевые ресурсы (пропускная способность отдельных межузловых каналов), обеспечиваемую соседством с другими виртуальными соединениями, агрегируемыми на различных участках пути в отдельных звеньях маршрута данного межабонентского соединения, и проявляющуюся в виде «внешних» очередей в транзитных узлах. В данной работе предложена математическая модель виртуального соединения, управляемого транспортным протоколом в режиме селективного отказа, учитывающая кроме фактора искажений в прямом и обратном
трактах передачи данных и механизмов повторных передач, обусловленных искажениями и истечением тайм-аута неприема ответа от получателя потока информации, еще и очереди от «внешних» межабонентских соединений.
1. Математическая модель транспортного протокола
Рассмотрим обмен данных между абонентами, соединенными многозвенным трактом передачи данных. Предположим, что выполняются следующие допущения. Узлы тракта соединены дуплексными каналами связи, имеющими одинаковые пропускные способности в обоих направлениях. Длина тракта передачи данных, выраженная в количестве участков переприема, равна Б. Обратный канал, по которому доставляются подтверждения отправителю о корректности приема последовательности сегментов данных, также имеет длину Б. Заданы вероятности искажения сегмента в канале связи для прямого Яп (ё), ё = 1, Б, и обратного -
Я0(ё), ё = 1,Б, направлений передачи каждого участка переприема. Тогда достоверности передачи сегментов данных вдоль тракта от источника до адресата и об-
Б Б
ратно составят ¥п = П( - К (ё)); =П( - К0(ё)). Время обработки сегментов
ё=1 ё =1
в узлах тракта одинаково. Взаимодействующие абоненты имеют неограниченный поток сегментов для передачи, а обмен выполняется сегментами одинаковой длины. Подтверждения получателя о корректности приема принимаемых данных переносятся в сегментах встречного потока. Полагаем, что повторная передача сегментов организована в соответствии с селективной процедурой отказа [1]. Считаем, кроме того, что потерь сегментов из-за отсутствия буферной памяти в узлах тракта не происходит. Задана функция вероятностей Ьп, п = 0, N, того, что каждый сегмент из потока анализируемого соединения в транзитном узле встретит очередь размером п < N, где N - максимальный размер очереди, определяемый емкостью буферных пулов транзитных узлов.
Тайм-аут длительности £ запускается перед началом передачи первого сегмента последовательности и фиксируется для всех сегментов в пределах ширины окна. Будем считать, что размер окна управляющего протокола определяется величиной Ш, а £>Ш - задает длительность тайм-аута ожидания подтверждения корректности доставки данных. После передачи очередного сегмента, протокол копирует его в очередь переданных, но не подтвержденных данных и запускает тайм-аут. Как только размер очереди становится равным ширине окна Ш, управляющий протокол приостанавливает передачу в ожидании получения квитанции или истечения тайм-аута £ ожидания подтверждения. При получении подтверждения из очереди удаляются сегменты, дошедшие до адресата без искажений. При истечении тайм-аута £ соответствующий сегмент передается повторно, и тайм-аут запускается вновь. Будем называть тактом время /, необходимое для вывода сегмента в линию. Такт определяется суммой времени вывода сегмента в линию, времени распространения сигнала в канале связи и времени обработки сегмента принимающим узлом. Тогда время получения отправителем сквозной квитанции распределено по геометрическому закону с параметром Е0 и длительностью такта дискретизации /.
Динамика очереди переданных, но не подтвержденных сегментов на узле-отправителе для различных режимов функционирования управляющего протоко-
ла в силу марковости дискретного процесса получения квитанций может быть описана двумерной цепью Маркова с дискретным временем и числом состояний по одному измерению равным длительности сквозного тайм-аута £, а по другому
- увеличенной на единицу максимальной длине очереди: N+1. Очевидно, что длительность тайм-аута должна быть достаточной для того, чтобы пакет с сегментом данных по прямому каналу достиг адресата и подтверждение получателя по обратному каналу было принято отправителем потока. Отсюда следует, что размер тайм-аута, выраженный в длительностях тактов (, должен быть не меньше суммы двойной длины пути и размера встреченной очереди в транзитном узле £ > 2Б+п. С учетом возможных повторных передач сегментов основного (прямого) потока и сегментов с подтверждениями во встречном потоке из-за искажений в отдельных звеньях тракта размер тайм-аута £ целесообразно выбирать с «запасом» на повторные передачи.
Квитанция на первый сегмент последовательности может поступить отправителю спустя время 5 > 2Б интервалов длительности t, необходимых для достижения первым сегментом адресата и возвращения отправителю подтверждения о корректности его приема. Если при передаче последовательности сегментов отправителем или переносе подтверждения пакет с сегментом в прямом тракте (или подтверждением в обратном) встретил очередь размера п, то время для получения квитанции (подтверждения) возрастает на размер встреченной очереди п и составит 5 > 2Б+п.
Процесс переноса информационного потока транспортным протоколом в однозвенном виртуальном канале моделируется марковской цепью [4]. Обобщение данной модели для пустого многозвенного тракта передачи данных выполнено в [5]. Функционирование виртуального соединения, управляемого транспортным протоколом, в нагруженном многозвенном тракте передачи данных с очередями сегментов перед отправляемыми данными или подтверждениями может быть описано марковизированным процессом, в котором размер очереди перед прямым или обратным потоком данных исследуемого соединения является дополнительной переменной марковского процесса. В состоянии цепи Маркова (1,п) источник отправил последовательность размера 1-п сегментов, которая в процессе переноса в одном из звеньев встретила очередь длиной п сегментов. Значениям координаты I = 0, Ж + п, п = 0, N, состояний цепи Маркова соответствует количество переданных, но не подтвержденных получателем сегментов и время от начала передачи последовательности, а значениям / = Ж + п +1, £ -1, п = 0, N, - время, в течение которого отправитель не активен и ожидает получение квитанции о корректности приема переданной последовательности из Ж сегментов. Обозначим через Р(/, п), I = 0, £ -1, п = 0, N, - вероятности состояний цепи Маркова. Тогда последовательность переданных, но не подтвержденных сегментов данных рассматриваемого виртуального соединения при очереди нулевой длины растет до состояния цепи Маркова с координатами (2Б -1,0) с вероятностью Ь0. Дальнейший
рост размера этой последовательности происходит с вероятностью Ь0 (1 - ^0). В состояниях (/, п), I = 2Б -1 + п, £ -1, п = 0, N, возможно получение отправителем квитанции, и в зависимости от результатов доставки отправитель передает новые сегменты (при положительной квитанции) либо повторно - искаженные. Поскольку отправленная последовательность сегментов исследуемого виртуального
соединения может встретить очередь ненулевой длины в любой момент процесса передачи (на пути последовательности до адресата или при переносе подтверждения отправителю информационного потока), то переход из состояния (/,0),
/ = 0, £ - 2, в состояние (/, п), / = 0, £ - 2, п = 1, N, происходит с вероятностью Ьп .
Обозначим через п-т переходные вероятности цепи Маркова, где (/,п) - координаты исходного, а (у,т) - измененного состояний цепи. Тогда динамику процесса передачи информационного потока в режиме селективного отказа можно задать следующими значениями переходных вероятностей:
(1)
Ь0, / = 0,2Б - 2, п = 0; у = / +1, т = 0;
Ь0 (1 - ^0), / = 2Б -1, £ - 2, п = 0; у = / +1, т = 0;
Ьт,/ = 0, £ - 2, п = 0; у = /,т = 1, N;
Ь0^0,/ = 2Б -1, Ж -1,п = 0;у = 2Б -1,т = 0;
Ь0 ^0, / = Ж ,Ж + 2Б - 2, п = 0; у = Ж + 2Б - 2 -/, т = 0;
Ь0^0,/ = Ж + 2Б -1,£ - 2,п = 0;у = 0,т = 0;
1, / = £ -1, п = 0, N; ] = 0, т = 0;
1, / = 0,2Б - 2 + п, п = 1, N; у = / +1, т = п;
1 - ^0, / = 2Б -1 + п, £ - 2, п = 1, N; у = / +1, т = п.
^0, / = 2Б -1 + п, Ж -1 + п, п = 1, N; у = 2Б -1, т = 0;
^0, / = Ж + п, Ж + п + 2Б - 2, п = 1, N; у = Ж + п + 2Б - 2 -/, т = 0;
^0, / = Ж + п + 2Б -1, £ - 2, п = ; у = 0, т = 0.
Пропускная способность виртуального соединения, управляемого транспортным протоколом, определяется как отношение среднего объема данных, передаваемых между двумя последовательными получениями квитанций, к среднему времени получения квитанции [4, 5]. Вклад в быстродействие виртуального соединения дают те состояния цепи Маркова, для которых возможно получение квитанции. Нормированная на единицу пропускная способность виртуального соединения в нагруженном тракте определяется отношением среднего количества сегментов данных, передаваемых отправителем между поступлениями двух последовательных квитанций, к среднему времени между поступлениями квитанций, выраженному в количестве интервалов длительности ^ 2 (Ж, £) = у/т . Поскольку квитанции переносятся в каждом сегменте независимо и поступают к отправителю каждый такт t при условии, что они не искажены на пути длины Б от получателя до отправителя информационного потока, то среднее время между приходами квитанций распределено по геометрическому закону с параметром Е0 и составит Т = 1/_Р0 . Средний объем передаваемых между поступлениями квитанций данных с учетом того, что каждый сегмент исследуемого соединения с вероятностью Ьп, п = 0, N, встречает очередь размера п и дает вклад в объем переданной информации обратно пропорциональный величине п +1, задается обобщением соотношения, приведенного в работе [4]:
1 ГЖ+2 Б-2+п_ £-1
£ 1р(1, п) + £ ЖР(1, п)
П +1
і=т-1+п
і =№+2Я-1+ п
Величины 1 и Ж определяются средним количеством сегментов, достигших адресата при селективной процедуре организации повторных передач искаженных сегментов: 1 = (1 - 2Б - п + 2)¥п , Ж = ЖЕп . Окончательно зависимость пропускной способности виртуального соединения от протокольных параметров ширины окна Ж и длительности сквозного тайм-аута £ примет вид
1
2 (Ж, £) = £ —1
п=0 п + 1
Ж+2 Б-2+п
£ (1 - 2Б + 2 - п)Р(1, п) + Ж £ Р(1, п)
1_ 1=2Б-1+ п 1=Ж+2Б-1+ п
. (2)
2. Анализ процесса передачи в однозвенном тракте
В стационарных условиях система уравнений равновесия, описывающая процесс переноса данных в виртуальном канале для тракта длины Б = 1, ширины окна Ж > 1 и длительности тайм-аута £ > Ж,£ > N + 2, согласно переходным вероятностям (1) имеет следующий вид:
Р(0,0) = *0
Ь0 £2 р(/,0) +£ £2 р(/,п)
_ /=Ж п=1 /=Ж+п
N
Р(1,0) = Ь0 Р(0,0) + *0
+ £ Р(£ -1, п);
п=0
N Ж-1+ п ~
Ь £ Р(/,0)+£ £ Р(/,п)
п=1 /=п+1
Р(/, 0) = Ь0(1 - *0)Р(/ -1,0), / = 2, £ -1; Р(0, п) = ЬпР(0,0), п = ;
Р( /, п) = ЬпР( /, 0) + Р( / -1, п), / = 1, п +1, п = 1, N;
(3)
(4)
(5)
(6)
(7)
(8)
Р(/, п) = ЬпР(/,0) + (1 - ^0)Р(/ -1, п), / = п + 2, £ - 2, п = 1, N;
Р(£ -1, п) = (1 - *0) Р(£ - 2, п), п = .
Найдем решение данной системы уравнений. Из уравнения (4) получаем р(/,0)=Р(1,0)[Ь0(1 - ^0)Г ', / =1, £ -1. Согласно (5), (6) и полученному выражению
1 -Ка - ]
для Р(/, 0), имеем Р(/, п) = Ьп
Р(0,0) + Р(1,0)-
/ = 0, п +1, п = 1, N.
1 - *0(1 - *0)
С учетом данного соотношения из (7) и (8) для произвольного п > 1 находим
Р(/, п) = Ьп (1 - *0)г-1 X 1 Ь0
Р(0,0) (1 - *0)п
Р(0,0) (1 - *0)п
+ Р(1,0)
*Ьп+1
г0°0
[1 - Ь0(1 - *0) ](1 - *0)п 1 - Ь0 (1 - Ь0)[1 - Ь0(1 - *0)]_
/ = п +1, £ - 2 :
+ Р(1,0)
Р(£ -1,п) = Ьп(1 - *0) 1
£-2 ,
, £-2
*0Ь0
п+1
[1 - Ь0(1 - *0) ](1 - *0)п 1 - Ь0 (1 - Ь0) [1 - Ь0(1 - *0)(
Подставляя найденные соотношения в (3), получаем выражение для Р(1,0):
Р(1,0) =
Р(0,0)Е
(1 - Е/
где Е = -
(1 - 60) (1 - *0(1 - Е>) ) - (1 - *0)(1 - -1 ]
N
(1 -Ь0)[1 -*0 + ^0ЬТ ] + *0^0 (1 -ЬТ-1)ьп (*0(1 -Е0))П
Из условия нормировки находим вероятность начального состояния (0,0):
ч5—п—1 V
Р(0,0) = (1- ^0)Т-1 М- ^0)Т-1
3 - 2(1+Ь - Е0) + *0(1 - Е0) N ґ
Е0
П=1
П —
+Е
Е0
2-(1+*0 -Е,)[1 -2*0(1-*0)]-*0(1-*0)[3 + *0(1 -*0)] (*0(1 --))5-2(1 -*0)
N
+£
П=1
Е0 (1-*0(1 - Е0))2 (1 - Е0)5-П-1 + (*0(1 - Е0) )
1-*0(1 - Е0)
*0П+1 (1
0 Vа 1 0 )
5-1
1 -
*0(1 - Е0) Е0 (1 - *0(1-Е0)) (1-*0(1 - Е0))2 (1-*0)(1-*0(1 - Е0))
Проанализируем полученное решение в ряде частных случаев. Пусть внешняя очередь в узлах отсутствует (Ь0 = 1). Тогда получаем известный результат [5]:
Р(0,0) =
Е0(1 - Е>)
Р(і,0) =
1 + (1 - Е0)Т-1 [Е0 - (1 - Е0)5-Т ] Е,(1 - Е,)і -1
1 + (1 - Е0)Т-1 [Е0 - (1 - Е0)5-Т ]
і = 1,5-1,
Р(і, п) = 0, і = 0,5 -1, п = 1, N .
Предположим, что поток исследуемого виртуального соединения всегда встречает очередь ненулевой длины (Ь0 = 0). При этом вероятность начального состояния и состояния (1,0) цепи Маркова принимает вид
Р(0,0) = -
Е)(1 - Е0)
1 + Е0 [1 + N ]+ (1 - Е0У
Е0-£ *П (1 - Е0)
П=1
5-Т-п
Р(1,0) =
Р(0,0) [1 - (1 - Е0)Т-1 ]
(1 - Е»)
(9)
где N = £ п*П . В случае очереди детерминированной длины (*П = 1, п > 1) для
П=1
начального состояния получаем
Р(0,0) =
1 + (п +1)Е0 + (1 - Е0)Т-1 [Е0 - (1 - Е0 )5-Т-П ] '
При неограниченной ширине окна (Ж =го), а следовательно, и длительности тайм-аута £ начальное состояние является невозвратным (Р(0,0) = 0), а вероятность состояния (1,0) принимает вид
Р(1,0) = F [1 -Ь0(1 -F0)]2
2-(1 + bo -Fo)[1 -2bo(1 -Fo)]-
N
(10)
-bo(1 - Fo) [3 + bo(1 - Fo) ] + NFo [1 - bo(1 - Fo) ] + Fo £ bn [bo(1 - Fo) ]n+
n=1
Если при этом соединение всегда соперничает за полосу пропускания с конкурентами (bo = o), то вероятность данного состояния преобразуется к
Fo
P(1,o) =-------(11)
v ' 1 + Fo (1 + N) v '
Для абсолютно надежного обратного канала связи (Fo = 1), длительности тайм-аута, превышающей максимальный размер очереди на удвоенную длину пути (S > N + 2), и старт-стопной управляющей процедуры (W = 1) множество вероятных состояний (i, n) определяется значениями индексов i = o, n +1, n = o, N и имеет вид
P(o,o) =----------г-1-------=, p(i, n) = bnP(o,o). (12)
3 - bo2 + (1 + bo) N
В тех же условиях при W > 2 значимыми будут только состояния
(i,n), i = 1,n +1, n = o, N :
Р(0,0) = 0, Р(1,0) =-=■, Р(/, п) = ЬпР(1,0), I = 1, п +1, п = 1, N. (13)
2 - Ь0 + N
3. Анализ процесса передачи в многозвенном тракте
В случае тракта произвольной длины аналитические зависимости для вероятностей состояний цепи Маркова при произвольной длительности тайм-аута £ > Ж получить не удается. Аналитическое решение удается найти в двух случаях: во-первых, при выполнении условия £ = Ж +1, Ж > 2Б + N и, во-вторых, если Ь0 = 0 . Система уравнений равновесия, описывающая процесс переноса данных в виртуальном канале для тракта произвольной длины Б для указанных в первом случае соотношений между протокольными параметрами £ и Ж , длиной тракта Б и максимальным размером очереди N, в соответствии с переходными вероятностями (1) принимает вид
N
Р(0,0) = £ Р(Ж, п);
P(i, o) = boP(i -1, o), i = 1,2D - 2;
W-1 N W-1
P(2D - 1,o) = bo P(2D - 2,o) + £ bo Fo P(i,o) + £ £ Fo P(i, n);
n=1 i=2D-1+n
n=o
Р(і,0) = Ь°(1 - -°)Р(і -1,0),і = 2Б,Ж; Р(0, п) = ЬпР(0,0), п = ЇЖ;
Р(і,п) = Р(і -1, п) + ЬпР(і, 0),і = 1,2Б -1 + п, п = 1, N;
Р(/, п) = (1 - )Р(/ -1,п) + ЬпР(/,0), / = 2Б + п,Ж -1,п = 1, N ;
Р(Ж, п) = (1 - )Р(Ж -1, п), п = Щ
Аналитическое решение данной системы уравнений равновесия определяется соотношениями
Р(і, 0) = Р(0,0)Ь°, і = 1,2Б - 2;
Р(і, 0) = Р(2Б -1,0) [Ь°(1 - -°)] 2, і = 2Б -1, Ж;
Р(і, п) = Р(0,0)Ьп
і+1
1 - Ь°
, і = 0,2 Б - 2, п = 1, N;
Р(і, п) = Р(0,0)Ь 1 Ь°
2Б-1 і г» (і г )1і-2Б+2
+ Р(2Б -1,0)Ьп 1 “[Ь0(1 - -°)]
1 - Ь0(1 - -°)
і = 2Б - 1,2Б -1 + п, п = 1, N;
Р(і, п) = Р(0,0)Ь,
1-Ь
1-Ь
Ь
’° (1 - —° )і-2 Б-п+' + Р(2 Б-1,0)Ь„(1-п
і-2 Б+1 ,
0
і-2 Б+2 0
Ьп+ -°° г°
_(1 - -°)п [1-Ь°(1--°)] 1-Ь° (1-Ь°)[1-Ь°(1--°)]_
і = 2Б+п,Ж-1, п = 1,=;
2 Б-1
Р(Ж, п) = Р(0,0)Ьп
1 - Ь0
- (1 - -°)Ж-2Б-п+1 + Р(2Б -1,0)Ьп (1 - - 0-2Б+1:
0
-2 Б+1
Ьп+1-у° г°
_(1 --0)п [1 -Ь0(1 --0)] 1 -Ь° (1 -Ь0)[1 -Ь0(1 --,)]_
Р(0,0)К
п = 1, N;
Р(2Б -1,0) =-
Ж-2Б+1
где К =
(1 - -0)
N Ь
[1 -Ь0С1 --°)]1 -Ь° -(1 -Ь°2Б-1)(1 --°)Ж-2Б+1
-1 (1 - -°)
N
ХЬ
п=1
1 Ь° + -Д^1 (1 --°)п ° °
Р(0,°) = (1- -°)
Ж-2 Б+1 ч
(1 - -°Г
1 - (1+-°)Ь°2Б-1 1 - Ь,
2Б+— ---------------+
-°
1-Ь
° п=1
-°
+K
2-(1+b0 -F0)[1 -2b0(1 -F0)]-b0(1 -F0)[3+b0(1-F0)] (^(1 --))W-2- -b0) +
F0 (1-b0(1-F0))2 1-b0(1-F0)
(1 - F0)W-2 D-n+2 + (b0(1 - F0)) b0”+1 (1 -F0f-2 D+2
N
+Xb
n=1
1-
b0(1 - F,) F0 (1 -b0(1-F,)) (1 - b0(1 - F0))2 (1-b0)(1-b0(1 - F0))
Рассмотрим найденное решение в ряде частных случаев. Нетрудно видеть, что данное решение при Б = 1 совпадает с решением, полученным при рассмотрении однозвенного тракта для протокольных параметров, связанных соотношением Ж = £ -1. При отсутствии внешней очереди в узлах (Ь0 = 1) получаем (известный результат [4]):
Р(0,0) =--------^о(1 ----------, Р(2Б -1,0) =----------------.
1 + (1 - F0 )№-2Б+1 [2БF0 -1] (1 - F0 )№-2Б+1
Если поток анализируемого виртуального соединения всегда встречает очередь ненулевой длины (Ь = 0), то вероятности состояний цепи Маркова определяются зависимостями
N
£ Ья (1 - Fо)f-2 Б-п+1 Р(0,0) = п=1 •
1 + N + -F- + ( 2D --F- ]£b„(1 - F0)
F0 V F0 ) n=1
P(0,0)
W-2 D-n+1
P(2D -1,0) =-
1 -X bn (1 - F0)
n =1
W-2 D-n+1
N
W - 2D-n+1
X bn (1 - F0)W
n=1
P(i, 0) = 0, i = 1,2D - 2,2D, W ; P(i, n) = bnP(0,0), i = 0,2D - 2, n = 1, N ;
P(i, n) = bn [P(0,0) + P(2D -1,0)], i = 2D -1,2D -1 + n, n = 1, N ;
P(i,n) = bn [P(0,0) + P(2D-1,0)](1 -F0)i-2D-n+l,i = 2D + n,W,n = 1N . (14)
При неограниченной ширине окна (Ж =го) состояния (/, п), / = 0,2Б - 2, п = 0, N, являются невозвратными (Р(/, п) = 0) и вероятность состояния (2Б -1,0) цепи Маркова принимает вид (10), что свидетельствует об инвариантности вероятности данного состояния к длине тракта передачи данных в указанных условиях. При этом для очереди детерминированной длины (Ьп = 1, п > 0 ) вероятности значимых состояний цепи Маркова, определяющих доступную межа-бонентскому соединению полосу пропускания тракта передачи данных, задаются
F F (1 - F )г'-2Б-п+1
зависимостями Р(2Б-1,0)=---------0----- и Р(г,п)=—-----------------0-, I>2Б-1+п .
1+F0(n + 1) 1+F0(n + 1)
Для детерминированного обратного тракта (F0 = 1) и ширины окна, превышающей максимальный размер очереди на удвоенную длину пути (Ж > N + 2Б), про-
странство значимых состояний (/, п) образует плоскость равнобедренного по ко-
ординатам I и п треугольника I = 2Б -1,2Б -1 + п, п = 0, N :
1
Р(2Б -1,0) =-----------=, Р(/,п) = ЬпР(2Б -1,0),/ = 2Б - 1,2Б -1 + п,п = 1,N . (15)
2 - Ь0 + N
Найдем решение для случая, когда поток анализируемого соединения всегда встречает очередь (Ь0 = 0 ),протокольные параметры связаны между собой, а также с максимальным размером очереди и длиной тракта передачи данных неравенствами £ > Ж + N и Ж > 2Б . Система уравнений равновесия при этом в соответствии с (1) имеет вид
Р(0,0) = £
п=1
Р(£ -1, п) + £ ^о Р(і, п)
і=Ж+2Б+ п-2
Р(і,0) = £ ^0Р(Ж + 2Б + п - 2 - і, п), і = 1,2Б - 2, п = 1, N;
п=1
N Ж-1
Р(2Б -1,0) = £ £ Е> Р(і, п);
п=1 і=2 Б-1+п
Р(0, п) = ЬпР(0,0), п = ;
Р(і, п) = Р(і -1, п) + ЬпР(і, 0), і = 1,2Б -1, п = 1, N;
Р(і, п) = Р(і -1, п), і = 2Б, 2Б -1 + п, п = 1, N;
Р(/, п) = (1 - F0)Р(1 -1, п),I = 2Б + п, £ -1, п = 1, N.
Для вероятностей состояний отсюда находим аналитическое решение:
Р(/,0) = Р(0,0)/7° ,/ = 1,2Б -2;
(1 - Fо)г
1 - (1 - F )Ж-2Б+1 Р(2Б -1,0) = Р(0,0)- 0^ •
(1 - Р0)Ш-1
Р(і, п) = Р(0,0)Ь» , і = 0,2Б - 2, п = 1, N; (1 - Р0)1
Р(0,0)Ь
Р(і,п) =—4 ’ 7 п ,і = 2Б -1,2Б + п -1,п = 1,N;
- — - Ж-1
(1 - Р0) Ъп (1 - ^,)і-2 Б-п+'
Р(і, п) = Р(0,0) ^^-, і = 2Б + п, £ -1, п = 1, N;
(1 - ^)Ж-1
Р(0,0) =-----------------------------------------------------N---------------------
1 + ^ ( N +1) + (1 - -2 Б+‘ [1 - (1 - ^ )2 Б-1 ] - £ Ь,, (1 - ^)£-2 Б-“+‘
п=1
При Б = 1 отсюда получаем (9), для детерминированного обратного тракта с учетом -0 = 0 - (15), а при Ж = го - приходим к
Р(1, п) = 0,I = 0,2Б - 2, п = 0, Ж,
р
Р(2Б -1,0) = 0
1 + Р0 (1 + Ж)
- р __________________________________ __________
Р(1,п) =------------------------" 0 _ч, I = 2Б -1,2Б + п -1, п = 1, Ж ,
1 + Р0 (1 + Ж) '
- р (1 - р )г-2Б-п+1 ______________________________
Р(1, п) = ------, I > 2Б + п, п = 1, Ж.
1 + Р0 (1 + Ж)
4. Анализ процесса передачи в многозвенном тракте с детерминированным обратным каналом связи
В предыдущих разделах найдены вероятностно-временные характеристики для виртуального соединения при Ж > 2Б . Остается открытым вопрос анализа быстродействия виртуального канала в случае Ж < 2Б , соответствующем «недогруженному» соединению. Уравнения равновесия, описывающие динамику очереди переданных, но не подтвержденных пакетов данных, при этом принимают вид
5-2 Ж 5-2 Ж ________
Р(0,0) = £ -0 р0 Р(1,0) + £ £ р0 Р(1,п) + £ Р(5 - 1,п),Ж = 1,2Б -1;
1=2 Б-2+Ж п=11=2 Б-2+п+Ж п=0
Ж
Р(1 , 0) = -0Р(1 -1,0) + -0р0Р(2Б - 2 + Ж -1 ,0) + £ р0Р(2Б - 2 + п + Ж -1, п),
п=1
| = 1,Ж -1, Ж = 2,2Б -1;
Р(1 ,0) = -0Р(| -1,0),I = Ж,2Б -1, Ж = 1,2Б -1;
Р(|, 0) = -0(1 - р0) Р(1 -1,0), I = 2Б, 5 -1;
Р(0, п) = -пР(0,0), п = 1Ж;
Р(| ,п) = Р(1 -1,п) + -пР(1,0),I = 1,2Б -1 + п,п = 1Ж;
Р(|, п) = (1 - р0) Р(1 -1, п) + -пР(1,0), I = 2 Б + п, 5 - 2, п = 1Ж;
Р(5 -1, п) = (1 - р0) Р(5 - 2, п), п = 1Ж.
Найдем решение данной системы уравнений для абсолютно надежного обратного тракта (р0 = 1), определяющего потенциально достижимый виртуальным каналом уровень быстродействия при недогрузке виртуального соединения, обусловленной малым размером окна (Ж < 2Б). Из приведенных уравнений в данных условиях нетрудно видеть, что состояния (I, п), I = 0, Ж - 2, п = 0, Ж являются невозвратными (Р(!, п) = 0). Для значимых состояний из системы уравнений равновесия получаем
Р(і, 0) = Р(Ш -1,0)ЬІ~Ш+1, і = Ш -1,2Б -1, Г = 1,2В -1;
і-Ш+1
Р(і, п) = Р(Ш - 1,0)й„ £ Ц = Р(Ш - 1,0)йп —-і=о 1
і = Ш -1,2В -1,п = ЇЖ,Ш = 1,2В -1;
- Ь0
(16)
Р(і, п) = Р(Ш - 1,0)Ьп
1 - ь02 В-Ш+1
1 - Ь0
1 = 2В - 1,2В -1 + п,п = 1,Ж ,Ш = 1,2В -1. Из условия нормировки находим
1
Р(Ш -1,0) = -
2 В-Ш+1
'0
(17)
1 - Ь0
N
Нетрудно убедиться, что при Ш = 1, В = 1 отсюда получаем найденное ранее решение (12).
5. Исследование доступной полосы пропускания
Анализ процесса передачи информационного потока в виртуальном канале, управляемого транспортным протоколом, показывает, что индекс быстродействия виртуального соединения при абсолютно надежных каналах связи в отдельных звеньях обратного тракта (= 1) и достаточной длительности тайм-аута согласно (2) и найденных вероятностей состояний цепи Маркова (12) для однозвенного тракта (В = 1) при Ш = 1 определится соотношением
~ Ы Ь "
Ь + (1 + Ьо) X -+7 _ П=1п+1_
2 (1,5) - 2
3 - Ьо2 + (1 + Ьо)N
а при Ш > 2 в соответствии с (13) составит
2 - Ь0 + N
" Ь 1 + £ ——
. п=1 п + 1,
(18)
Для случая, когда поток исследуемого однозвенного виртуального соединения всегда встречает очередь ненулевой длины (Ь0 = 0), его быстродействие согласно (2) и (9) определится выражением
2 (Ш, £) = ¥п
N и
Е2 [1 -(1 -е)Ш-1 ]+£[1 -(1 -^)Ш -ШЕ(1 -е)£1 ]
_ N '
1 + Е>(1 + N) + (1 - Е0)Ш-1 - (1 - Е0)Ш-£ Ьп (1 - Е0) ^
п=1
В случае многозвенного тракта передачи данных (В > 1) при Ь0 = 0, 5 = Ш +1 и Ш > 2В + N в соответствии с найденными вероятностями состояний цепи Маркова (14) из (2) получаем
Р2 г0
= Р„
г (V V+1)=
N
V-2 Б-п+1
1-£Ь (1-Ро)
П=1
N ь
+ £-+7{1 - (1 - Ро)№-2°-”+2 [1 - Ро( V - 2Б - п + 2)]}
п=1П +1
_ N
1+Ро(1 + N) + (2БРо -1)£Ь,(1 -Ро)И п=1
Если протокольные параметры удовлетворяют неравенствам £ > V + N , V > 2Б , то доступная полоса пропускания составит
N Ь
Ро2 [1 - (1 -РоТ-2Б+1 ] + £[1 - (1 -Ро)Ш -№Ро(1 -Ро)£-2Б-П+1 ]
г (V, £) = рп------------------------------------------2="+-N------------------.
1 + Ро(1 + N) + (1 - Ро)Ж-2Б+1 - (1 - Ро)Ж - £ Ь (1 - Ро)£-2Б-п+1
п=1
Если V > 2Б и Ро = 1, то вероятности состояний марковской цепи задаются выражением (15) и доступная виртуальному соединению доля полосы пропускания оказывается инвариантной к длине тракта и совпадает с зависимостью (18). Если же при этом еще и очередь имеет детерминированную длину (Ьп = 1, п > о), то доля пропускной способности многозвенного тракта (2), доступная виртуаль-
Рп
ному соединению составит 2 (V, £) =—— . Согласно найденному решению (16),
п +1
(17) при Ро = 1 и V < 2Б быстродействие виртуального канала (2) определяется соотношением
1 Ь2Б-V+1 N ь
Ь2Б-Ш + 1_^---------£_ Ьп
о
1 - Ьо п= п +1
2^,£) = Рп------------о . ^Б-V+1_. (19)
2Б - V + 2 - Ьо2Б-ж+1 +---о-------N
о 1 - Ьо
Для постоянной очереди ненулевой длины (Ьп = 1,п > о) получаем 2(V,£) =
Рп
=------------п----------, а при V = 1, Б = 1, Ьо = 1 приходим к известному ре-
(п +1)[2Б - V + 2 + п] о
Рп
зультату [3]: 2(1, £) = . Таким образом, потенциально достижимая скорость
передачи в виртуальном канале, проложенном в тракте передачи данных, за ресурсы которого соперничают множество сетевых абонентов, определяется соотношениями (18) и (19) при V > 2Б и V < 2Б соответственно. Численные исследования показывают, что доступная виртуальному соединению полоса пропускания для V > 2Б и прочих равных условиях практически инвариантна к длине тракта передачи данных (см. рис 1), ощутимо снижаясь от области насыщения лишь при V = 2Б и Ро < 1. В случае же V < 2Б доступная полоса пропускания недогружена, и эффективная скорость передачи данных значительно снижается (см. рис. 2 и 3). С ростом конкуренции между абонентами за полосу пропускания тракта передачи данных размер очереди увеличивается и скорость информационного переноса быстро падает (см. рис. 4). Зависимость быстродействия от достоверности передачи подтверждений в обратном канале связи приведена на рис. 5.
От достоверности доставки данных в прямом канале связи индекс быстродействия имеет линейную зависимость.
Рис. 1. Зависимость доступной полосы пропускания от длины тракта передачи данных при £=35, №=16, Ь1 = 1, Рп = 1
Рис. 2. Зависимость быстродействия виртуального соединения от ширины окна при £=3о, Б=1, Ь5 = 1, Рп = 1
Рис. 3. Зависимость быстродействия виртуального соединения от ширины окна при £=3о, Б=2, Ь5 = 1, Рп = 1
Рис. 4. Зависимость быстродействия виртуального соединения от длины очереди при £=3о, W=1, Б=4, ЬN = 1, Рп = 1
о о,2 о,4 о,6 о,8 1 Р0
Рис. 5. Зависимость доступной полосы пропускания от достоверности передачи данных в обратном канале при £=3о, Б=4, Ь1 = 1, Рп = 1
При неограниченном росте размера окна, а следовательно, и тайм-аута доступная абонентскому соединению доля пропускной способности тракта (2) имеет вид
ч Рп [1+ (п +1)Ро2]
2 (да, ^) =---- -------------— .
(п +1)[1 + (п + 1) Ро ]
Отсюда для Ро = о и Ро = 1 получаем одинаковые значения доступной полосы
2 ( ) Рп
2 (да, да) = -
п+1
г л/п+2-1 „
а в точке Ро =---------- показатель доступной полосы пропускания достигает
п +1
минимума:
Z (да, да) = Fn
2V n + 2 - 2 (n +1)2
Заключение
Проведенный анализ направлен на исследование доступной полосы пропускания виртуального соединения в тракте передачи данных, разделяемом конкурирующими абонентами. Предложена математическая модель нагруженного тракта при заданном распределении длин очередей протокольных блоков данных або-нентов-соперников в виде цепи Маркова, описывающей динамику объема переданных, но не подтвержденных данных исследуемого виртуального соединения. Найдены распределения состояний марковской цепи при различных соотношениях между протокольными параметрами и характеристиками тракта. Получены аналитические зависимости быстродействия виртуального соединения для различных условий его функционирования. Численные исследования быстродействия виртуального канала в селективном режиме повторной передачи показали, что скорость передачи в виртуальном соединении в значительной мере определяется интенсивностью потерь данных, распределением размера очередей в транзитных узлах, а также соотношением между длиной тракта и размером окна. Очевидно при этом, что длительность сквозного тайм-аута должна превосходить сумму удвоенной длины тракта и совокупной длины очередей перед информационным потоком взаимодействующих абонентов виртуального соединения. Полученные результаты позволяют утверждать, что при заданном размере окна показатель пропускной способности виртуального соединения возрастает с увеличением длительности тайм-аута ожидания сквозного подтверждения и практически достигает теоретического предела при насыщении по протокольному параметру W для значений S, превосходящих ширину окна на 3-5 тактов длительности t.
ЛИТЕРАТУРА
1. Богуславский Л.Б. Управление потоками данных в сетях ЭВМ. М.: Энергоатомиздат, 1984. 168 с.
2. Gelenbe E., Labetoulle J., Pugolle G. Performance evaluation of the HDLC protocol // Comput. Networks. 1978. V. 2. P. 409-415.
3. Боровихин Е.А., Коротаев И.А. Анализ функционирования и оптимизация протокола HDLC // Автоматика и вычисл. техника. 1993. № 2. C. 47-51.
4. Сущенко С.П. Аналитические модели асинхронных процедур управления звеном передачи данных // Автоматика и вычисл. техника. 1988. № 2. C. 32-40.
5. Кокшенев В.В., Сущенко С.П. Анализ быстродействия асинхронной процедуры управления звеном передачи данных // Вычислительные технологии. 2008. Т. 13. № 5. С. 61-65.
6. Nimbe L. Ewald, Andrew H. Kemp. Analytical model of TCP new reno through CTMC // Computer Performance Engineering, Proc. 6th European Performance Engineering Workshop, EPEW 2009; London, UK, July 9-10, 2009. P. 183-196.
7. Padhye J., Firoiu V., Towsley D.F., et al. Modeling TCP Reno Performance: A Simple Model and Its Empirical Validation // IEEE/ACM Transactions on Networking, April 2000. V. 8. №. 2. P. 133-145.
Михеев Павел Андреевич Кокшенев Владимир Владимирович Сущенко Сергей Петрович Томский государственный университет E-mail: [email protected]; [email protected]; [email protected]
Поступила в редакцию 12 февраля 2013 г.
Mikheev Pavel A., Kokshenev Vladimir V., Sushchenko Sergey P. (Tomsk State University). Transport protocol selective acknowledgements analysis in loaded transmission data path.
Keywords: transport protocol, transmission data path with queues, mathematical model, Markov chain, transport connection throughput, window size, through time out, packet loss.
A mathematical model of the transport protocol with selective acknowledgements for multihop transmission data path with queues is proposed. The protocol behavior is modeled through the two-dimensional Discrete-Time Markov chain which is used for estimation of transport connection throughput. The model takes into account the influence of window size, retransmission time out, probability of packet loss and queue size on transport connection throughput. Probabilities of Markov chain states are found. Analytical expressions of bandwidth available for transport connection in different load conditions were obtained. Analysis of transport protocol throughput for different protocol parameters and transmission data path conditions has been conducted.