ОБ ОДНОЙ МОДЕЛИ ОЦЕНКИ ПРОИЗВОДИТЕЛЬНОСТИ ШИРОКОПОЛОСНОГО ГИБРИДНОГО КАНАЛА СВЯЗИ НА ОСНОВЕ ЛАЗЕРНОЙ И РАДИОТЕХНОЛОГИЙ
В. М. Вишневский, О. В. Семенова
Закрытое акционерное общество «Научно-производственная фирма "Информационные и сетевые технологии"», 129626, Москва, Россия
УДК 004.057.4, 519.872
Рассмотрена система массового обслуживания для моделирования широкополосного гибридного канала связи (атмосферного оптического и радиоканала). Радиоканал является резервным и подключается в тех случаях, когда оптический канал недоступен. Время доступности (недоступности) атмосферного оптического канала аппроксимируется гиперэкспоненциальным распределением второго порядка на основе статистического анализа метеоданных. Для рассмотренной модели получено стационарное распределение состояний системы, а также основные характеристики производительности, в том числе среднее время ожидания и вероятность доступности гибридного канала. Приведены численные примеры.
Ключевые слова: широкополосные беспроводные сети передачи данных, гибридный канал, система массового обслуживания, матрично-аналитический подход.
In the paper, we present a queueing system to model a free space optics-radio hybrid channel. The radio-frequency channel acts as a backup when the free-space optical link is down. Time when optical link is/is not used mainly depends on the weather conditions and is approximated by the hyperexponential distributions of order two. For the model proposed, the stationary state probabilities are obtained and performance characteristics are derived including the mean waiting time and the channel accessibility. The results are illustrated by numerical examples.
Key words: broadband wireless networks, hybrid channel, queueing system, matrix-analytical approach.
Введение. Быстрое и непрерывное увеличение количества пользователей Интернетом и соответствующее увеличение объемов передаваемой мультимедийной информации выдвигает в ряд первоочередных задач повышение производительности и надежности каналов связи для решения проблемы "последней мили". Работы по увеличению быстродействия и качества широкополосных беспроводных каналов связи ведутся во многих исследовательских центрах и стандартизирующих организациях всего мира [1]. В конце 2009 г. Институт инженеров по электротехнике и электронике (Institute of electrical and electronics engineers (IEEE)) опубликовал высокоскоростной стандарт IEEE 802.11n [2]. На базе этого стандарта разработана аппаратура, обеспечивающая номинальную скорость передачи информации в беспроводном канале связи до 300 Мбит/с, что значительно превышает быстродействие существующих каналов WiFi и WiMAX (номинальная скорость существующих беспроводных каналов составляет 54 Мбит/с). Разрабатываются также следующие новые перспективные стандарты беспроводной связи: IEEE 802.11n VHT (very high throughput), на базе
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (код проекта № 10-07-90006-Бел).
которого будут реализовываться беспроводные сети со скоростью передачи информации до 500 Мбит/с, завершение планируется в 2012 г.; LTE Advanced (релиз 10) - стандарт сотовой связи 4G со скоростью передачи мультимедийной информации до 1000 Мбит/с в нисходящем канале и 500 Мбит/с в восходящем канале, завершение планируется в 2012 г.; IEEE 802.11m (релиз 2.0) - мобильный WiMAX со скоростными характеристиками, близкими к характеристикам LTE, завершение планируется в 2011 г.
Одним из перспективных направлений повышения производительности и надежности беспроводного канала связи является создание гибридных каналов на базе лазерной и радиотехнологий. Оборудование, реализующее гибридный канал, объединяет преимущества лазерных атмосферных каналов связи (высокая скорость передачи информации в дуплексном режиме) и широкополосных радиосредств (возможность эффективной работы в неблагоприятных погодных условиях: туман, снег и т. д.). Основным режимом работы комбинированных приемопередатчиков является лазерный режим. При ухудшении соотношения сигнал - шум осуществляется переход на широкополосный радиорежим, а после восстановления качества атмосферного оптического канала происходит обратный переход. При этом радиоканал может либо прекращать работу при использовании лазерного канала, либо параллельно использоваться для передачи данных. Большинство работ, в которых исследуются модели гибридного канала, связано главным образом с имитационным моделированием [3-6].
Среди работ, посвященных построению математических моделей и их анализу, следует отметить [7], однако и в этой работе авторы ограничиваются лишь нахождением вероятности нарушения связи в оптическом канале. Заметим также, что в отличие от моделей, рассматриваемых в [3-6], в модели [7] предполагается, что оба канала (оптический и радиоканал) работают параллельно.
В данной работе рассмотрен случай, когда основным каналом связи является оптический канал. Когда его работа нарушается и не возобновляется дольше определенного периода времени, передача данных осуществляется по радиоканалу, до тех пор пока работа оптического канала не возобновится. Смена каналов передачи данных описывается с помощью модели массового обслуживания с чередующимися скоростями обслуживания (для оптического и радиоканала). Использование одной скорости обслуживания (время работы оптического или радиоканала) ограничено случайной величиной, распределение которой может быть определено путем статистической обработки метеорологических данных той местности, в которой предполагается функционирование гибридного канала. Рассматривается пример, когда длительности использования каждого из каналов могут быть аппроксимированы гиперэкспоненциальным распределением второго порядка.
1. Статистическая обработка метеоданных и определение функции распределения периодов доступности и недоступности оптического канала. В опубликованных к настоящему времени работах, посвященных исследованию качества открытого оптического атмосферного канала (ОАК) и характеристик погодных явлений, влияющих на доступность оптического канала, оценивается лишь вероятность того, что метеорологическая дальность видимости (МДВ) меньше некоторого значения. Однако в действительности важно не только то, в течение какого промежутка времени МДВ ниже некоторого показателя, но и то, насколько часто происходит смена погодных условий, приводящая к переходу между оптическим и радиорежимами работы гибридного канала связи. Известно также, что на доступность оптического канала оказывает влияние большое количество других, не связанных с МДВ, факторов, таких как сцинтилляция, юстировка, колебания опоры, влия-
гиперэкспонента (р = 0,8, ^\ = 0,223, = 0,01): ' 1=23,37,т2=3949,т З=1163216
экспериментальные данные:
т= 23,37, т = 3949, тъ= 987978
^ экспонента: т] = 23,37, те2 =1092, т3 = 76588
10 20 30 40 50 60 70 Длительность работы резервного канала без учета переключений, мин
80
а
гиперэкспонента (р= 0,65, ^ ° = 0,345, "Ц ° = 0,014): т =26,08, т2= 3355, т3= 693351
экспериментальные данные'
т= 26.08 , т2 = 3355 , т 3= 592961
«»л у' экспонента: т= 26,08, т= 1360, тъ= 106377
Г _
10 20 30 40 50 60 70 80 Длительность работы ОАК без учета переключений, ч
Рис. 1. Распределения периодов недоступности (а) и доступности (б) ОАК
ние солнца и др. Помимо перечисленных факторов, связанных с атмосферой и погодными явлениями, большое значение имеет непосредственно конструкция приемопередатчика: длина оптической волны, количество лучей, расходимость луча, чувствительность приемника и мощность передатчика, наличие системы автоматического наведения и др. Таким образом, очевидно, что невозможно создать модель, достаточно точно учитывающую все перечисленные характеристики гибридного канала связи, поэтому для оценки характеристик разрабатываемого гибридного канала был проведен продолжительный натурный эксперимент, позволивший получить данные для оценки времен доступности и недоступности атмосферного оптического канала связи для конкретного оборудования при протяженности канала связи, расположенного в Рязанской области, 4,5 км.
На рис. 1 приведены полученные в натурном эксперименте длительности работы резервного и оптического атмосферного каналов, а также две функции, аппроксимирующие данное распределение. Длина гибридного канала 4,5 км, натурный эксперимент проводился в Рязанской области. Из рис. 1 следует, что экспоненциальное распределение, построенное по среднему значению функции, является недостаточной аппроксимацией времен доступности канала, в то время как гиперэкспоненциальное распределение, построенное по первому и второму начальным моментам, хорошо согласуется с экспериментальными данными.
Для построения гиперэкспоненциального распределения второго порядка для выборки данных {Х1, • ••,Х*} (например, периодов доступности) используем следующий подход [8].
Вычисляются первые три начальных момента ш\, т2 и т3 по формулам
1 *
т* = , * = 1, 2 3
к=1
Определяются величины
2
т3 — т1т2 т1т3 — т2 /1 = --17Х", /2 = —-"Х" •
т2 — т1 т2 — т1
Если выполнены неравенства
т2 > т\, т3 > т1т2, т1т3 > т^, /1 > 4/2,
то параметры гиперэкспоненциального распределения с плотностью вида
F(t) = pYie-Ylt + (1 - p)72e-724 определяются следующим образом:
=_2__=_2_ = Yi(l2mi - 1)
71 = fi + VW-4T2' 72 = fi — f—f' P = Y2 - 7i "
В противном случае (если хотя бы одно из неравенств не выполнено) параметры Y1, y2 и p находятся из условий p = m\/m2, 0 < p < 1 и m1Y2 — 1+p > 0, где j1 = pj2/(m1 y2 — 1+ p); Y2 = {ш1(1 — p) — \J(1 — p)pi(l2 — m<2)^ /(m2 — m2p); величина p выбирается таким образом,
чтобы минимизировать разность \p/m1 + (1 — p)/m\ — m3\ .
Таким образом, считаем, что для выбранного устройства и канала протяженностью 4,5 км, расположенного в Рязанской области, времена доступности и недоступности ОАК могут быть приближены гиперэкспоненциальным распределением с плотностью
F (t) = p^1e-Ylt + (1 — p)Y2e-Y2 4
(параметры Y1, y2 и p для интервалов работы оптического и резервного каналов приведены на рис. 1).
При построении функций распределения для канала произвольной длины можно использовать формулу W(L) = aiLbi, где W(L) - вероятность наступления погодных условий, при которых МДВ меньше расстояния L (в километрах); ai,bi - константы для конкретной географической точки, предложенные в работе [9]. Исходя из пропорциональности длины канала и МДВ, при которой оптический канал становится недоступным, можно вывести формулу для параметров распределения длительности работы резервного канала произвольной длины:
YP = Y[(/o)(!f) Ь .
Здесь i =1, 2; b =1,46 для Рязанской области; !0 =4,5 км; ! - протяженность канала, км; Yp -параметры гиперэкспоненциального распределения длительности работы резервного канала (радиоканала) протяженностью 4,5 км, которые приведены на рис. 1,а в минутах.
Поскольку при изменении длительности работы резервного канала изменяется длительность работы оптического канала, будем считать, что параметры распределения длительности работы ОАК определяются следующим образом:
* = *<«( йЬ-
Здесь yO(!0) - параметры гиперэкспоненциального распределения длительности работы ОАК протяженностью 4,5 км, которые приведены на рис. 1,б в часах.
Предложенная методика позволяет определить искомые характеристики канала в произвольном регионе на основе численных данных, полученных в натурном эксперименте для канала фиксированной длины.
2. Математическая модель гибридного беспроводного канала связи. Модель гибридного канала представляет собой однолинейную систему массового обслуживания с двумя возможными скоростями обслуживания и ограниченным временем их использования.
Входящий поток заявок в систему - стационарный пуассоновский с параметром Л. Число мест для ожидания не ограничено.
При использовании k-й скорости (k-го режима) время обслуживания заявки распределено экспоненциально с параметром ¡k, k = 1, 2. При этом полагаем, что ¡1 > ¡2. Время использования k-й скорости ограничено и характеризуется абсолютно непрерывной случайной величиной ^k, имеющей плотность гиперэкспоненциального распределения второго порядка:
Fk (t) = pk Y^e-^4 + (1 — pk bfe-^,
параметры которого определены с помощью подхода, описанного в п. 1.
Если время первого режима завершается (ОАК становится недоступным), то текущее обслуживание заявки (передача пакета) прерывается и система переходит на второй режим работы (использование радиоканала). Сразу после перехода на второй режим обслуживания пакетов не производится. Система должна выждать время Q2. Если за это время второй режим работы не завершился (£2 > Q2), то по истечении времени Q2 начинается обслуживание заявки (передача пакета) на второй скорости (по радиоканалу). При этом предполагается, что заявка, обслуживание которой по оптическому каналу было прервано, обслуживается заново.
По истечении времени работы второго режима (когда ОАК вновь становится доступным) система, продолжая обслуживать заявки во втором режиме, наблюдает за доступностью оптического канала по времени. Если в течение определенного времени ОАК доступен, то по истечении этого времени система переходит на первый режим работы (оптический канал). При этом заявка, в ходе обслуживания которой произошла смена режима, обслуживается заново на новой скорости. Время, которое должен превысить период доступности ОАК, определяется каждый раз, когда ОАК становится доступным, из экспоненциального распределения с параметром Q1.
В п. 3 данная модель исследуется методом вложенных цепей Маркова с целью нахождения стационарных вероятностей состояний системы, которые вычисляются в п. 4 с помощью матрично-аналитического подхода.
3. Вложенная цепь Маркова. Предположим, что система функционирует в стационарном режиме. Условия существования стационарного режима будут получены ниже.
Предполагаем, что система может находиться в следующих четырех состояниях: 1 - работа в первом режиме (использование ОАК); 1' - переключение с первого режима на второй (с оптического на радиоканал); 2 - работа во втором режиме (использование радиоканала); 2' - переключение со второго режима на первый (с радиоканала на оптический).
Отметим на временной оси моменты завершения обслуживания заявок, моменты начала и моменты окончания использования каждого из режимов работы. Обозначим через tn, n > 1 последовательность таких моментов. Под состоянием системы в момент tn будем понимать вектор
Xn = (in,mn,m'n),
где in - число заявок в системе в момент времени tn + 0; mn - состояние обслуживающего прибора в момент tn — 0; m'n - состояние обслуживающего прибора в момент tn + 0. В данном случае будем полагать, что состоянию (i, m, m) соответствует завершение обслуживания заявки m е{1, 2, 2'}.
Введем стационарные вероятности состояний системы:
п(г,т, т') = Иш Р{гп = г,тп = т, т'п = т'}, г > 0, т,т' €{1,1', 2, 2'}.
п—
Пусть также
о
/(к) = 1 ^е-Л'е-^' (Рке-^ + (1 — рк)е-^') Цк¿1 — о
вероятность того, что за время обслуживания заявки в режиме к в систему поступят г заявок,
оо
д(к) = 1 ^е-Л'е--' (рк^е-^ + (1 — рк¿1 — о
вероятность того, что за время обслуживания заявки в режиме к время использования режима истекает (оптический канал становится недоступным, если использовался режим 1 (к = 1), или доступным, если использовался режим 2 (к = 2)) и за неполное время обслуживания в систему поступят г заявок,
<2
= / ^Ге'" + (1 — Р2ИГе-221') А —
о
вероятность того, что время переключения с первого режима на второй меньше, чем Q2, и за время переключения в систему поступят г заявок,
(1)< — ^ -42)52] _
^(2) = е'Л< ре'^5 + (1 — р2)е'Л)
вероятность того, что время переключения с первого режима на второй больше, чем Q2, и за время переключения в систему поступят г заявок,
^ е-Л'е-51^2е-№ '¿1 —
вероятность того, что за время обслуживания заявки при переключении со второго режима на первый в систему поступят г заявок,
^ е-Л'е-<1' Qle-»2 '¿1 —
т* = ~Ц~е
о
вероятность того, что за время обслуживания заявки при переключении со второго режима на первый это переключение завершается и за неполное время обслуживания в систему поступят г заявок.
Данные вероятности (кроме которые заданы явно) могут быть вычислены следующим образом:
Ак) Цк А ,л ч Цк А , 0
) = рк ---тг+т + (1 — рк) —-п+т, к = 1,2,
(7(1) + Цк + а) (т(2) + Ц к + х)
п,(1) \г ~(2) хг
д(к) = рк-ЖЛ-+ (1 — рк)-ЖЛ-, к =1, 2,
»1 , (1) чт 1 ^ ^к> , (2) \г+1' ' '
(7(1) + Цк + а) (т(2) + Ц к + х)
г
^ = ^ ( *--^^ -- ) +
4 + М (т2Ч + А)4 ^=0 т!^ + А)-
т22)Аг(1 - р^) ( 1_ „-(7(2>+аЮ2 ^ ят
-А)4
+ 72 А - Р2) 1 _ е-(7(2)+А)д2 у^ Я2
т22) + А Ы2) + А)4 т=0 т!(722) + А)4
Г 4 =
(Я + ^2 + А)4+1 (Я + ^ + А)4
т = ——:-"ттт+г, г ^
Заметим, что за время переключения со второго режима на первый заявки не только поступают в систему, но и уходят из нее. Поэтому для удобства будем полагать, что в каждый момент завершения такого переключения с вероятностью
(1) (2) 71 т1
а = Р1~ш—;т +(1 - —;т
71 + Я1 71 + Я1
система начинает использовать первый режим работы, а с дополнительной вероятностью 1 — а остается работать во втором режиме. Величина а есть вероятность того, что начавшееся время первого режима больше времени переключения, и, таким образом, с этой вероятностью система может переключиться на первый режим работы.
Используя введенные выше вспомогательные вероятности, составим уравнения равновесия для стационарных вероятностей п(г,т,т') состояний системы
4+1
(1)
п(г, 1,1) = £ (п(к, 1,1) + п(к, 2', 1) + п(к, 1', 1)) /—+1+ к=1
+ (п(0,1,1) + п(0, 2', 1) + п(0,1', 1)) /; (1)
г
п(г, 1,1') = £ (п(к, 1,1) + п(к, 2', 1) + п(к, 1', 1)) д-к+ к=1
+ (п(0,1,1) + п(0, 2', 1) + п(0,1', 1)) ^(1); (2)
г
(г, 1',г) = £ п(к, 1,1>-к, г = 1, 2; (3)
п(г, 1 ,г) = п(
к=0
4+1
п(г, 2, 2) = £ (п(к, 2, 2) + п(к, 2', 2) + п(к, 1', 2)) /Д+1+ к=1
+ (п(0, 2, 2) + п(0, 2', 2) + п(0,1', 2)) /(1); (4)
г
п(г, 2, 2') (п(к, 2, 2) + п(к, 2', 2) + п(к, 1', 2)) д-к+
к=1
+ (п(0, 2, 2) + п(0, 2', 2) + п(0,1', 2)) ^(2); (5)
4+1
п(г, 2', 2') =5>(к, 2', 2') + п(к, 2, 2'))г-к+1 +
+ (п(0, 2', 2')+ п(0, 2, 2')) А+А^тг; (6)
г
п(г, 2', 1) = ^(п(к, 2', 2') + п(к, 2, 2'))т-к а+ к=1
+ (п(0, 2', 2') + п(0, 2, 20) (^ХО-1тг-1 + ; (7)
г
п(г, 2', 2) = ^(п(к, 2', 2') + п(к, 2, 2'))т-к(1 — а) + к=1
+ (п(0, 2', 2') + п(0, 2, 2')) (тг-1 + QA+:Qf)!{г=о}) , (8)
где 1{А} - функция-индикатор события А,
рт ^ (1) + (1 рт) (2) , т 1, 2,
т л Т т
((т) = Зтдг—г 1{г>о} + (1 — 8т)!{г=0}, % > 0.
Vг $тдг-
Заметим, что из системы (1)-(7) можно непосредственно исключить равенства (3). В результате получим
г+1 г+1
К1) , 1 1Л /(1)
п(г, 1,1) = £ (п(к, 1,1) + п(к, 2', 1)) /-+1 + ^п(к, 1,1'+ к=1 к=1
+ (п(0,1,1) + п(0, 2', 1)) / + п(0,1,1') [(81/(1) — /г+1)^ + ; (9)
г г
п(г, 1,1') = £ (п(к, 1,1) + п(к, 2', 1)) д(--к + ^ п(к, 1,1')/- + к=1 к=1
(1)
+ (п(0,1,1) + п(0,2', 1)) ^ + п(0,1,1') (д(1) — к01)д(1) + ; (10)
г+1 г+1
п(г, 2, 2) = £ (п(к, 2, 2) + п(к, 2', 2)) к+1 + ^ п(к, 1,1'+ к=1 к=1
+ (п(0, 2, 2) + п(0, 2', 2)) / + п(0,1,1') ((/ — /1+)1)^02) + /й) ; (11)
п(г, 2,1') = £ (п(к, 2, 2) + п(к, 2', 2)) д-к + ^ п(к, 1,1')д— +
& + £ п(к, 1,1')д(-к к=1 к=1
+ (п(0, 2, 2) + п(0, 2', 2)) *<2) + п(0,1,1') (д(2) — к^д(2) + ^(2)к02)) ; (12)
г+1 А
п(г, 2', 2') = У>(к, 2', 2') + п(к, 2, 2'))гг-к+1 + (п(0, 2', 2') + п(0, 2, 2'))-—-тг; (13) к=1 А + Ql
п(г, 2', 1) = ^(п(к, 2', 2') + п(к, 2, 2'))т-ка+
+ (п(0, 2', 2') + п(0, 2, 2'))
Аа
Я\<а
А + Q1 т-1 + А + Q1
{4=0}
п(г, 2', 1) = ^(п(к, 2', 2') + п(к, 2, 2'))ш-к(1 - а) +
к=1
+ (п(0,2,2) + п(0,2,2)) ( Аг1-0)ш- + ^ -а) 1{г=о}
А + Ql
А + Ql
где
(14)
(15)
/
(т)
£4/1, !й = £ ^зй, г > 0, ш =1, 2.
к=0
к=0
4. Матрично-аналитический метод вычисления стационарных вероятностей.
Для вычисления стационарных вероятностей п(г, ш, ш'), г > 0, ш,ш' Е{1,1', 2, 2'} (за исключением п(г, 1',г), которые могут быть вычислены по формулам (3)) используем матрично-аналитический подход, основанный на идее сенсорных цепей Маркова (см., например, [1011]). Близкие алгоритмы изложены в работах [12-15].
Введем в рассмотрение векторы стационарных вероятностей
П = (п(г, 1,1),п(г, 1,1'),п(г, 2, 2),п(г, 2, 2'),п(г, 2', 2'),п(г, 2', 1),п(г, 2', 2)).
Тогда систему уравнений (9)-(15) можно записать в векторно-матричной форме
4+1
П = ^ ПкРк,4, г > 0, к=0
где матрицы Р4г (г > 0, I > шт{г — 1, 0}) имеют вид
Р
0,1 =
Аг Оэхэ"
О 2x4 Вг
Сг О2х3
Р4,1 = Уг-г+1 для I > г — 1,
/(1) /г зй 0 0 0 0 0
/,(1) зЙ /!(2) (2) 3-1 0 0 0
0 0 /(2) (2) 3((_1 0 0 0
Уг = 0 0 0 0 гг ашг-1 (1 — а)шг-1
0 0 0 0 гг ашг-1 (1 — а)шг-1
/(1) г! 0 0 0 0 0
0 0 /!(2) (2) 3г-1 0 0 0
I > 0,
(16)
А,
г =
(1)
Уг
(1)
(1)
+ ^
г(1)
(^—3((1 0
0
,(2)
0
) У((2) !((2) + ^0>((2) — зГ)
(2)
,(2)
B,
i —
A r \a m I Qia т X(l-a) m + Ql(l-a) г '
X+Qi ri XTQTmi-1 + A+QT T{i=0} "A+Q7"mi-1 + A+Qi T{i=0}
A
Aa I QTa t A(1-a) . Qi(1-a) T
A+Qi ri a+QTmi-i + A+QTT{i=0} ЛОГmi-i + QA+Q1-2T(i=0}
C, —
f V((1) 0 0
(2) „,(2)
0 sif
v;
/ > 0,
yim) = (/ - /¿i )4m) + f№, m = 1, 2, i > 0.
В (16) g-\ = g— = g-1) = g-l = m_i = 0.
В известных алгоритмах, основанных на идее сенсорных цепей Маркова, предполагается, что счетная компонента цепи Маркова, описывающей поведение системы (в данном случае компонента in цепи Xn, n > 1), за один шаг совершает переход не более чем на единицу влево.
Пусть G(k) - матрица вероятностей переходов конечных компонент цепи Маркова Xk,k > 1 за время, в течение которого значение счетной компоненты уменьшится с k + 1 до k (k > 0). Матрицы G(k) (k > 0) удовлетворяют соотношениям
G(k) — Pk+i,k +Y, Pk+i,гО(г-1)С(г-2) • • • G(k), k > 0. i=k+1
(17)
В силу того что матрицы Рг, I зависят лишь от разности I — г, матрицы С(к) в (17) не зависят от к (будем полагать, что все С(к) = С) и определяются как решение матричного уравнения
g — PiG
i-i+1
(18)
i=i- i
Уравнение (18) решается итерационно, например с помощью итерационной схемы
-1
Go — I,
Gn — II Pi i iG-W Рг i i-1, n > 1,
(19)
а матрица G определяется как G = lim Gn.
Согласно рассуждениям, приведенным в [10, 11], векторы стационарных вероятностей п (i > 0) вычисляются следующим образом:
Пк — ПоФк, к > 1.
(20)
Здесь Ф0 — I,
к-1
Фк — (I - рк,к)-1, к > 1
i=0
матрицы Piкк (к > 0, i — 0,к) определяются равенствами
Р,к = Pi,iGl-k, г = 0,k, k > 0,
-ii-k
Pi, l G l=k
вектор по есть решение системы уравнений
По(1 - Ро ,о) = 0; (21)
те
По Фг12 = 1, (22)
г=0
12 = (1, 2,1,1,1,1,1)4.
Заметим, что уравнение (22) получено из условия нормировки
те 2
п(г, 1,1') + п(г, 2, 2') + п(г, 2', 2') + (п(г, r, r) + п(г, 2',r) + п(г, 1',r))
г=о
r=1
с учетом равенств (3), (20).
5. Стационарное распределение вероятностей состояний системы в произвольный момент времени. Под состоянием системы в произвольный момент времени t > 0 будем понимать состояние процесса Xt = {it,mt},t > 0 (it - число заявок в системе в момент времени t; mt - состояние обслуживающего прибора, m Е {1,1', 2, 2'}).
Стационарное распределение вероятностей состояний системы в произвольный момент времени определяется как предельное распределение процесса Xt (t > 0). Процесс Xt (t > 0) является полурегенерирующим процессом с вложенным процессом марковского восстановления {in,tn},n > 1 (см. [16]). Из анализа результатов работы [16] следует, что достаточным условием существования стационарных вероятностей
p(i, m) = lim Р{it = i,mt = m}, г > 0, m g{1, 1', 2, 2'} t^tt
является существование стационарного распределения вложенной цепи Маркова Xn (n > 1), а стационарные вероятности p(i,m) (г > 0) состояний системы в произвольный момент времени определяются равенствами
р(г, r) = Ä(n(0,r, r) + п(0,1 ,r) + п(0, 2,r)) х ( вгh-i + + (1 - Pr^ ) +
Л + 7r Ä + 7r
+äJ](n(k,r,r)+ п(к, 1',r)+ п(к, 2',r)) fi-k, г > 0, r = 1, 2, k=1
i
р(г, 1') = Ä^>(k, 1,1')hi-k, г > 0,
к=о
р(г, 2') = Л(п(0,2', 2') + п(0,2,2' S-i + Ä+Ql) +
i
+Л (п(к, 2', 2') + п(к, 2, 2'))дг-к, г > 0,
где
Г = / Ме'Л*е'^г' ие'7(1)' + (1 — р„)е-^') Л =-Р^-^ + (1 — РгА ,
] %! ^ ) (А + ЦГ + 7Г1))г+1 (А + ЦГ + 7Г2))г+1
к = 7 (М-е'Л' (р2е'^(1)' + (1 — р2)е'7(2)^ ¿1
о
р2Аг р2АгТ_^_е-л<2 е'^(1)<2 +
Р2А , Шч г—т+1 е е +
9г
(А + 721))г+1 2 т=0 т!(А + 721))г-т+1
+ (1 — р2)Аг _ (1 _ р2)хг V_^_е-л<2е-^(2)<2
(А + 722))г+1 т=о т!(А + 7(2))г-т+1
о
= [ ^е-Л'е-^е-51'М = --А* . + , г > 0.
У %! (А + Ц2 + Ql )г+1
о
6. Условия существования стационарного режима. Для определения условий существования стационарного режима рассматриваемую систему представим как систему с отдыхами обслуживающего прибора (под отдыхом прибора понимается время между моментом, когда завершился первый режим работы, и моментом, когда система переключилась на первый режим). Заметим, что дисциплина обслуживания заявок в системе является ограниченной (время обслуживания очереди ограничено значением случайной величины ).
Для систем с отдыхами обслуживающего прибора необходимое и достаточное условие существования стационарного режима получено в работе [17] и имеет вид
р1 + АУ/Ы < 1, (23)
где р1 = А/ц1; У - средняя длительность отдыха; Ы - среднее максимальное число заявок, которые могут быть обслужены в очереди между соседними отдыхами. Заметим, что, поскольку для системы в [17] обслуживание заявок во время отдыха прибора не проводится (а в данном случае - наоборот), условие (23) является достаточным для существования стационарного режима рассматриваемой системы. Величина Ы определяется по формуле
лж ( Р1 I 1 — Р1 Ы = Ц1 — +
\Т1 ) 71 ) Найдем среднюю длительность отдыха У:
У = У1>, 1 + (1 — щ^У, 2 + У2 + У2' )•
Здесь
Щ1 = 1 — Р2е-7(1)<2 — (1 — р2)е-7(2)<2 —
вероятность того, что переключение с первого режима на второй завершилось переходом на второй режим;
V = Р2 | 1 ~ Р2
'1 = "НУ + ~12Т
средняя длительность переключения с первого режима на второй, которое завершилось переходом на первый режим; Уу,2 = - средняя длительность переключения с первого режима на второй, которое завершилось переходом на второй режим;
у = Р2 + 1 - Р2
ъ' 7(
средняя длительность второго режима; У2/ - средняя длительность переключения со второго режима на первый (с учетом возможных "возвратов" на второй режим), которое завершилось
переходом на первый режим.
Найдем величину У2/. С вероятностью а имеем У2/ = 1/^1 (а - вероятность того, что
переключение со второго режима на первый завершается переходом на первый режим).
С вероятностью (1 — а)т_1 а за время с момента завершения второго режима до момента
перехода на первый последует т — 1 "возвратов" на второй режим, и на т-м переключении
произойдет переход на первый режим (т > 2). Средняя длительность такого периода равна
¿(т + (т — 1)У^ (1 — а)т_1а,
и, таким образом, 1 , 1 ч
У2/ = -(тГ + У2 (1 — а)) . а )
Средняя длительность отдыха У определяется равенством
У = Ч % + + <1 — -К * + ак + У2
Таким образом, получаем следующий результат.
Утверждение. Для существования стационарного режима в рассматриваемой системе достаточно, чтобы выполнялось неравенство
, лА Р1 , 1 — рЛ ^ 1
Р1 + Р1М — + <1
7. Характеристики производительности системы. Вычислив стационарное распределение вероятностей состояний системы в произвольный момент времени, находим характеристики производительности:
1) время, в течение которого система работает в г-м режиме:
те
Чт = ^p(i,г), г = 12;
г=0
2) среднее число заявок в системе при ее работе в г-м режиме:
те
^ — / гр(г,г)
г=0
3) среднее число заявок в системе в произвольный момент времени:
Ь = Ь
1 + Ь1/ + Ь2 + Ь2/;
Ьт = ^ гР(г,г), г е{1,1', 2, 2'};
Рис. 2. Зависимость недоступности гибридного канала от времени переключения из оптического в радиорежим: а - при мгновенном обратном переключении; б - при пятиминутном обратном переключении
4) среднее время пребывания заявки в системе (вычисляется по формуле Литтла)
T = L/А.
8. Численные примеры. Рассмотрим систему А = 2000, ц1 = 15 258, ц2 = 2441. Параметры распределений, аппроксимирующих периоды доступности (недоступности) ОАК, получены в п. 1 и имеют следующие значения: y(1) = 9,57 • 10-5, y(2) = 4,019 • 10-6, y(1) = 3,72 • 10-3, y22) = 1,684 • 10-4 c-1, pi = 0,65, p2 = 0,8.
Основной характеристикой производительности гибридного канала является его доступность (среднее время, в течение которого канал доступен для передачи пакетов). Доступность гибридного канала может быть вычислена по формуле
те
R =1 — ^ P(i, 1').
г=0
На рис. 2 показана зависимость недоступности канала связи (величины 1 — R) от времени Q2 переключения на резервный радиоканал. Из рис. 2 следует, что доступность канала 99,99 % обеспечивается при значении Q2, не превышающем 9,5 с, а при времени переключения до 1 с может быть обеспечена доступность 99,999 %. Заметим, что среднее время переключения с радиоканала на ОАК практически не оказывает влияния на недоступность канала, поэтому может принимать различные значения в широком диапазоне.
На рис. 3 представлены результаты для среднего числа переключений с ОАК на радиоканал в единицу времени при различных значениях Q2. Среднее число переключений определяется формулой
- те те -|
£n(i, 1,1') ^е"^2 + (1 — p2)e~Y22)Q2) + J] n(i, 2', 1) ,
-г=0 г=0
Nor = (а +-2-)
V mo + mR>
2 4 6 8 10
Время переключения из оптического в радиорежим, с
2 4 6 8 10
Время переключения из оптического в радиорежим, с
Рис. 3. Зависимость числа переключений между режимами от времени переключения из оптического в радиорежим: при мгновенном обратном переключении; б - при пятиминутном обратном переключении
a
где mo = pi/7(1) + (1 - Pi)/Yi2),
mR = P2/Y21) + (1 - P2)/j¥) - средние длительности доступности (недоступности) ОАК.
На рис. 4 приведена зависимость среднего числа заявок в системе от интенсивности входного потока.
„(2)
3,0
12'5 &
£
2,0
! 1.:
U
1,0
0,5
0,5 1,0 1,5 2,0 2,5 3,0 3,5 4,0
(1/Х)-10 3, с
Рис. 4. Зависимость среднего числа заявок в системе от интенсивности входного потока
Список литературы
1. Вишневский В., Портной С., Шлхнович И. Энциклопедия WiMAX. Путь к 4G. М.: Техносфера, 2009. 470 с.
2. Part 11: Wireless LAN Medium Access Control (MAC) and Physical Layer (PHY) Specifications Amendment: Enhancements for Higher Throughput // IEEE 802.11n-2009. 2009. Oct. 502 p.
3. Akbulut A., Gokhan H., Ari F. Design, availability and reliability analysis on an experimental outdoor FSO/RF communication system // Proc. of the 7th Intern. conf. on transparent optical networks. S. l., 2005. P. 403-406.
4. Derenick J., Thorne C., Spletzer J. On the deployment of a hybrid free-space optic/radio frequency (FSO/RF) mobile ad-hoc networks // IEEE/RSJ Intern. conf. on intelligent robots and syst., comp. sci. and engng. Bethlehem (USA): Lehigh Univ., 2005. P. 3990-3996.
5. Wang D., Abouzeid A. A. Throughput of hybrid radio-frequency and free-spaceoptical (RF/FSO) multi-hop networks // Inform. theory and applications workshop. USA, 2007. P. 1-8.
6. Nadeem F., Leitgeb E., Kvicera V., et AL. Simulation and analysis of FSO/RF switch over for different atmospheric effects // Proc. of the Intern. conf. on telecomm. (ConTEL 2009). Graz (Austria): Inst. of Broadband comm., 2009. P. 39-43.
7. Letzepis N., Nguyen K. D., Guillen i Fabregas A., Cowley W. G. Outage analysis of the hybrid free-space optical and radio-frequency channel // IEEE J. Select. Areas Comm. 2009. V. 27. P. 1709-1719.
8. Клзимирский А. В. Методика аппроксимации MAP-потоков при помощи MAP-потоков второго порядка // Информационные системы и технологии (IST'2002): Материалы 1-й Междунар. конф. Дополнение. Минск: БГУ, 2002. С. 27-31.
9. Зеленюк Ю. И., Огнев И. В., Поляков С. Ю., Широбакин С. Е. Влияние погодных условий на надежность атмосферной оптической связи // Вестн. связи. 2002. № 4. С. 136-138.
10. Grassmann W., Heyman D. Equlibrium distribution of block-structured Markov chains with repeated rows //J. Appl. Probability. 1990. V. 27. P. 557-576.
11. Kemeni J., Snell J., Knapp A. Denumerable Markov chains. N. Y.: Van Nostrand, 1966. 348 p.
12. Breuer L., Dudin A. N., Klimenok V. I. A retrial BMAP/PH/1 system // Queueing Systems. 2002. V. 40. P. 433-457.
13. Ramaswami V. A. A stable recursion for the steady state vector in Markov chains of M/G/1 type // Comm. Statist.-Stochastic Models. 1988. V. 4. P. 183-188.
14. Dudin A., Semenova O. Stable algorithm for stationary distribution calculation for a BMAP/SM/ 1 queueing system with Markovian arrival input of disasters //J. Appl. Probability. 2004. V. 42, N 2. P. 547-556.
15. Семенова О. В. Устойчивый алгоритм расчета стационарного распределения системы обслуживания BMAP|SM|1 с марковским потоком сбоев и двумя режимами работы // Автоматика и вычисл. техника. 2004. № 1. С. 75-84.
16. Cinlar E. Introduction to stochastic processes. New Jersey: Prentice-Hall, 1975. 400 p.
17. Fricker C., Jaibi R. Monotonicity and stability of periodic polling models // Queueing Systems. 1994. V. 15. P. 211-238.
Вишневский Владимир Миронович - д-р техн. наук, проф., ген. директор ЗАО НПФ "Информационные и сетевые технологии";
тел. (495) 643-57-72, e-mail: vishn@inbox.ru Семенова Ольга Валерьевна - канд. физ.-мат. наук, ст. науч. сотр. ЗАО НПФ "Информационные и сетевые технологии";
тел. (499) 789-96-91, e-mail: olgasmnv@gmail.com Дата поступления - 20.03.10