УДК 004.057.4
С. Ю. Шаров
Московский физико-технический институт (государственный университет)
Аналитическая модель гибридного канала связи с холодным резервированием
Рассматривается гибридный канал передачи данных, построенный на базе лазерного атмосферного канала, который резервируется радиоканалом. Строится математическая модель, позволяющая описать предлагаемое устройство. Кроме того, в работе предлагается комплексный анализ построенной модели, позволяющий получить основные характеристики гибридного канала с холодным резервированием. Находятся условия существования стационарного режима.
Ключевые слова: гибридный канал, лазерный атмосферный канал, радиоканал миллиметрового диапазона радиоволн, математическое моделирование.
1. Введение
Лазерный атмосферный канал (англ. FSO - Free Space Optics) все активнее завоевывает рынок беспроводных устройств связи, обеспечивающих высокую пропускную способность. На больших расстояниях оптоволоконные кабели по-прежнему предпочтительнее для обеспечения высокой скорости передачи данных. Однако благодаря тому, что системы FSO значительно дешевле, а время развертывания системы намного меньше, применение этой технологии вместо волоконной оптики считается эффективным на расстояниях до 5 км.
Благодаря высокой несущей частоте, находящейся в диапазоне 300 ТГц, технология лазерных атмосферных каналов позволяет осуществлять высокоскоростную передачу данных, необходимую для таких приложений, как потоковое видео и звук, видео по запросу, конференцсвязь с использованием телевизионных каналов и т.д. Лазерный атмосферный канал может быть также широко использован для организации корпоративных сетей. Однако подобные оптические беспроводные линии подвержены сильному влиянию со стороны среды, по которой распространяется сигнал, которая в абсолютном большинстве случаев является земной атмосферой. Среди всех факторов, ослабляющих сигнал в таком канале, туман вносит наибольший вклад, т.к. размер капелек тумана сравним с длиной волны в используемом оптическом диапазоне; в то же время затухание из-за дождя относительно менее значимо для таких беспроводных оптических линий. Сравнимую высокую скорость передачи данных может обеспечить и радиоканал в миллиметровом диапазоне, работающий в диапазоне частот от 30 до 300 ГГц. Для таких радиоканалов наибольшую проблему представляют дождь, град и мокрый снег, значительно уменьшая мощность сигнала в радиоканале. Особенно важно отметить, что туман не оказывает значительного влияния на функционирование миллиметрового радиоканала, т.к. повышенная влажность вызывает затухание менее чем 5 дБ/км.
Гибридные системы связи находят широкое применение для решения проблемы последней мили и в ad-hoc сетях. Применение различных комбинаций лазерного атмосферного канала и радиоканала предлагалось во множестве областей, например, было предложено использовать гибридный канал в мобильных ad-hoc сетях [1], совместно с другими наземными и спутниковыми линиями связи, например, [2, 3]. Проведён ряд экспериментов по сбору статистики работы гибридных систем. Так, в работе [4] в течение 14 месяцев собиралась статистика для двух дублирующих друг друга каналов: FSO и миллиметрового радиоканала. «Дублирующие» в данном случае означает, что одни и те же данные одновременно передавались по двум параллельным каналам и измерялся коэффициент доступности как
для каждой линии в отдельности, так и для всего гибридного канала в целом. Согласно собранной статистике время доступности гибридного канала составило 99.93%. Это значение в большинстве случаев приемлемо для построения корпоративных сетей.
В целом рассматриваются два основных подхода для построения гибридных систем. Первый заключается в том, что все время параллельно работают обе линии, тем самым приводя к потере 50% доступной пропускной способности, вызванной дублирующей передачей данных. Во втором подходе используется механизм переключения таким образом, что большую часть времени функционирует только одна линия, тем самым позволяя уменьшить объем излишне передаваемых данных и использовать резервный канал только по мере необходимости, в случае, когда в основном канале возникают ошибки. Механизм переключения в гибридных каналах, используемый во втором подходе, обсуждался, например, в [5]. Однако во всех этих работах предложенный механизм переключения не использует пропускную способность резервного канала, если доступен основной.
В последнее время стали рассматривать и третий вариант. Предлагается механизм переключения, позволяющий увеличить полезное использование пропускной способности каналов, применяя механизм распределения нагрузки, когда оба канала доступны. Кроме того, если линия недоступна, то механизм продолжает следить за каналом до тех пор, пока связь не восстановится. В большинстве работ механизм переключения основывается на сравнении уровня полученного сигнала с пороговым значением. Однако такой механизм переключения обладает большим количеством недостатков, поэтому в последнее время появились работы, предлагающие альтернативные способы переключения.
В данной работе рассматривается гибридная система, в которой лазерный атмосферный канал (ЛАК) резервируется с помощью радиоканала. Это позволяет приблизить систему к классу операторской надежности. В работе представлена аналитическая модель подобного гибридного канала. Разработанная модель позволяет проанализировать большинство интересующих характеристик и получить комплексную картину работы системы в целом, при этом нет необходимости иметь в наличии реальное устройство.
Остальная часть работы организована следующим образом. Следующий раздел посвящен описанию предлагаемой математической модели. Далее проводится ее комплексный анализ, рассчитываются основные характеристики системы. В конце работы приводятся полученные численные результаты и дается их краткий анализ.
2. Математическая модель гибридного беспроводного канала связи с холодным резервом
В данном разделе построена и исследована стохастическая модель, описывающая функционирование гибридного канала связи. Модель представляет собой систему массового обслуживания типа М/М/1 с двумя возможными скоростями обслуживания.
Входящий поток заявок в систему - стационарный пуассоновский с параметром Л. Число мест для ожидания не ограничено.
При использовании к-т скорости (к-то режима) время обслуживания заявки распределено экспоненциально с параметром к = 1, 2. При этом полагаем, что > ц,2- Время использования к-и скорости ограничено и характеризуется абсолютно непрерывной случайной величиной имеющей плотность гиперэкспоненциального распределения второго порядка:
^ (*) = рк е-^4 + (1 - Рк )7*
Если время первого режима завершается (ЛАК становится недоступным), текущее обслуживание заявки (передача пакета) прерывается и система переходит на второй режим работы (использование радиоканала). Сразу после перехода на второй режим обслуживание пакетов не производится. Система должна выждать время Q2, и если за это время второй режим работы не завершился (£2 > $2) то по истечении времени начинается обслуживание заявки (передача пакета) на второй скорости (по радиоканалу). При этом
предполагаем, что заявка, обслуживание которой по оптическому каналу было прервано, обслуживается заново.
По истечении времени работы второго режима (когда ЛАК-канал вновь становится доступным) система, продолжая обслуживать заявки во втором режиме, наблюдает доступность оптического канала по времени. И если ЛАК доступен в течение определенного времени, то по истечении этого времени система переходит на первый режим работы (оптический канал). При этом заявка, в ходе обслуживания которой произошла смена режима, обслуживается заново на новой скорости. Время, которое должен превысить период доступности ЛАК, определяется всякий раз, когда ЛАК становится доступным, из экспоненциального распределения с параметром Ql.
3. Стационарное распределение вероятностей состояний системы
Предположим, что система функционирует в стационарном режиме. Условия существования стационарного режима будут получены позже.
Предполагаем, что система может находиться в следующих четырех состояниях:
1 - работа в первом режиме (использование ЛАК),
1' - переключение с первого режима на второй (с оптического на радиоканал),
2 - работа во втором режиме (использование радиоканала),
2' - переключение со второго на первый режим (с радиоканала на оптический). Будем наблюдать поведение системы в моменты завершения обслуживания и моменты окончания и смены режимов работы. Обозначим через Ьп п-й такой момент, п ^ 1. Под состоянием системы в момент Ьп будем понимать вектор
Хп = {гп,тп,т'п) ,
где гп — число заявок в системе в момент времени Ьп + 0 тп — состояние обслуживающего прибора в момент Ьп — 0 т'п — состояние обслуживающего прибора в момент Ьп + 0. Здесь мы полагаем, что состоянию (г, т, т) соответствует завершение обслуживания заявки, т € {1, 2, 2'}.
Введем стационарные вероятности состояний системы:
п(г, т, т') = Нш Р {гп = г, тп = т, т'п = т'\ , г ^ 0, т,т' € {1,1', 2,2'} .
п^те
Пусть также
те
= / ^е-^е-^ (Рке-^ + (1 — Рк)е-^) цкМ о
- вероятность того, что за время обслуживания заявки в режиме к в систему поступит г заявок,
те
^ = I ^е-^* (рк^е-^ + (1 — РкЬГе-^) <И
о
- вероятность того, что за время обслуживания заявки в режиме к время использования режима истекает (оптический канал становится недоступным, если использовался режим 1 (к = 1), или доступным, если использовался режим 2 (к = 2) и за неполное время обслуживания в систему поступит г заявок,
Я*
^ = / ^е-* (р^е-^ + (1 — Р2)722)е-^) <И о
- вероятность того, что время переключения с первого режима на второй меньше, чем ^2, и за время переключения в систему поступит г заявок,
2) = Й^!(>-*Я* ^-^Я* + (1 — р^е-^Я*)
- вероятность того, что время переключения с первого режима на второй больше, чем ^2, и за время переключения в систему поступит г заявок,
те
п = 1 М. е-^е-^^е-^
о
- вероятность того, что за время обслуживания заявки при переключении со второго режима на первый в систему поступит г заявок,
те
№ = ■ (А()'--
- вероятность того, что за время обслуживания заявки при переключении со второго режима на первый это переключение завершается и за неполное время обслуживания в систему поступит г заявок.
Заметим, что за время переключения со второго режима на первый заявки как поступают в систему, так и уходят из нее. Поэтому для удобства будем полагать, что в каждый момент завершения такого переключения с вероятностью
(1) (2) л л а = -+ (1 - Р^^Г-
7( + $1 71 + $1
система начинает использовать первый режим работы, а с дополнительной вероятностью 1 — а остается работать во втором режиме. Уточним, что величина а есть вероятность того, что начавшееся время первого режима больше, чем время переключения, и, таким образом, с этой вероятностью система может переключиться на первый режим работы.
Далее, используя введенные выше вспомогательные вероятности, составим уравнения равновесия для стационарных вероятностей п(г, т, т') состояний системы для г ^ 0:
г+1
-к(г, 1,1) = ^(п(к, 1,1) + ф, 2', 1) + ф, 1', 1)^+1+ к=1
+(эт(0,1,1) + я-(0, 2', 1) + эт(0,1', 1))51/^(1), ж(г, 1,1') = ^(ж(к, 1,1) + ж(к, 2', 1) + ж(к, 1', 1))^-+
(1)
(1) + г-к
к=1
+(эт(0,1,1) + Я"(0, 2', 1) + ^(0,1', 1)Я1),
(2)
(г, 1,г) = £ п(к, 1, М-к, г = 1,2, (3)
■К(г, 1,Г) = П(К, 1, 1 г
к=0
+1
п(г, 2, 2) = ^2(Ф, 2, 2) + п(к, 2', 2) + п(к, 1', 2))^-^+ к=1
+(эт(0, 2, 2) + я-(0, 2', 2) + эт(0,1', 2))82^2), ж (г, 2, 2') = 2, 2) + ж(к, 2', 2) + ж(к, 1', 2))д^к+
(2) + г-к
к=1
+(эт(0, 2, 2) + я-(0, 2', 2) + ^(0,1', 2))у\2),
(4)
(5)
+1
ж(1, 2', 2') = ^(ъ(к, 2', 2') + ж(к, 2, 2'))гг-Л+1+
к=1 (6)
+(*(0, 2', 2') + *(0, 2, 2'))
л +
о
ж (г, 2', 1) = ^(ж(к, 2', 2') + ж(к, 2, 2'))mi-ka+
k=i
+(ж(0, 2', 2') +ж(0, 2, 2')) ^
Xa
X + Qi
mi-1 +
Qia X + Qi
1 {i=0}
(7)
ж(г, 2', 2) = ^(ж(к, 2', 2') + ж(к, 2, 2'))mi-k(1 - a) +
k=1
+(ж(0, 2', 2') + ж(0, 2, 2')) ( Л,(1, ^ Шг-х +
V ^ + 41
где 1{А] — функция-индикатор события А,
Qi(1 -а) X + Qi
{ =0}
(8)
S m — pm
x+
+ (1 - Pm)
x+7m
, m — 1, 2,
(m) (m) т ,л \T • ^ n
Vi — Sm5i-1 7{i>0} + (1 - «mjI{i=0}, ^ ^ 0.
Заметим, что из системы (1) - (8) можно непосредственно исключить равенства (3). Для вычисления стационарных вероятностей ж (г ,т,т'), г ^ 0, т,т' € {1,1', 2, 2'} (за исключением ж (г, 1', г), которые могут быть вычислены по формулам (3)) применим матрично-аналитический подход, который основан на идее сенсорных цепей Маркова (см., например, [6]). Близкие алгоритмы изложены в работе [7].
Введем в рассмотрение векторы стационарных вероятностей
щ = (ж (г, 1,1), ж(г, 1,1'),ж(г, 2,2),ж(г, 2, 2'),ж(г, 2', 2'),ж(г, 2', 1),ж(г, 2', 2)). Тогда систему уравнений (1) - (8) можно переписать в векторно-матричной форме
i+1
TTi
J2*kPk,i, г^ о, k=0
где матрицы Рц, 0 l ^ min {г - 1, 0}, имеют вид
P0
0, —
' Л
^2x4
C
Оэхэ' Вг
02хз
Pi,i — i+i для l ^ i - 1, где
Ъ —
(i) J1 У i-i 0 0 0 0 0
f(i) J1 M) Уг-i f(2) J1 (2) 0 0 0
0 0 f( 2) J1 (-2) 0 0 0
0 0 0 0 n ami-i (1 - a)mi-i
0 0 0 0 n ami-i (1 - a)mi-i
f(i) J1 (1) У i-i 0 0 0 0 0
0 0 f(2) J1 (2) 0 0 0
l > 0,
(9)
X
X
Л =
л
«1Ъ
(1)
,(1)
'(1) 9((1) + ^ - 9((1))
0
,(2)
8 2 /(
(2)
,((2)+^((2) - ,;2))
(2)
В( _
л + Я1 л
(
Ла
Л + Я1 л а
т(-1 +
(
Л + Я1 Л + Я1
С( _
т~ 1 +
Я1а
л + Я1
Я 1а Л + Я1
{(=0}
Л(1 - а)
т(-1 +
Я1 (1 -а)
{(=0}
{(=0}
Л + Я1 ( Л + Я1
Л(1 -а) Я1 (1 - а)
-Ш(-1 + х , ^-1{(=0}
Л + Я1
Л + Я1
1 ( 0
(1) (1)
0
0 821]
0
(2) (2)
I > 0,
У
(т) _ (е Ат) — ( о 1
'т J г
е (т)\и(т) <1 г+1 -
- Г(т>0т) + п
(т) <1 г+1 ,
т _ 1,2, г^ 0.
о (<л\ (1) (2) -(1) -(2) п
1 _ 1 _ 1 _ 1 _ т-1 _ 0
Отметим, что известные алгоритмы, основанные на идее сенсорных цепей Маркова, предполагают, что счетная компонента цепи Маркова, описывающей поведение системы (в нашем случае компонента гп цепи Хп, п ^ 1), совершает переход не более чем на единицу влево за один шаг.
Сенсорная цепь Маркова хп> п ^ 1 для цепи Хп, п ^ 1 строится следующим обра-
п
N * ^ 0. Состояние сенсорной цепи совпадает с состоянием цепи Хп, п ^ 1, если гп _ 0, N *. Когда состояние компоненты г п цеп и Хп, п ^ 1 становится больше, чем N *, компонента гп сенсорной цепи принимает значение N*, а конечные компоненты этой цепи изменяют состояние так же, как и соответствующие компоненты цепи Хп, п ^ 1.
Зафиксируем целое число N * ^ 0 и для цепи Маркова Хп,п ^ 1 построим сенсорную цепь хп-1 п ^ 1 с пространством состояний {0,1,..., N*} счетной компоненты гп.
Обозначим через Р(М \ 1,1 — 0, N * матрицы одношаговых вероятностей переходов сенсорной цепи хп1 п ^ 1. Нетрудно видеть, что
Р
(м *)
Р,
1 — 0^ *, I _ 0, N * - 1.
(м *)
Определим матрицы Р(1 ) при I _ N*. Для этого необходимо проследить переходы
исходной цепи Хп, п ^ 1 в области гп > N * от момента, когда значение компоненты
гп стало больше, чем N*, и до момента, когда гп вновь вернется в состояние N*. Решив
,(м *)
систему матричных уравнений для Р( п(г, т, т).
можно вычислить стационарные вероятности
4. Стационарное распределение вероятностей состояний системы в произвольный момент времени
Под состоянием системы в произвольный момент времени £ ^ 0 будем понимать состояние процесса Хг _ {ц, тг}, £ ^ 0 гДе Ч — чисто заявок в системе в момент времени тг — состояние обслуживающего прибора, тг € {1,1', 2, 2'}.
Стационарное распределение вероятностей состояний системы в произвольный момент времени определяется как предельное распределение процесса Хг, £ ^ 0. Процесс Хг, £ ^ 0, является полурегенерирующим процессом с вложенным процессом марковского восстановления {гп, Ьп}, п ^ 1 (см. [7]). На основании результатов работы [7] достаточным условием существования стационарных вероятностей
р(г, т) _ Нш Р {и _ г,ть _ т} , г^ 0,
т
€ {1,1', 2, 2'}
Хп
п ^ 1, а стационарные вероятности р(г,т), г ^ 0, состояний системы в произвольный момент времени могут быть выражены через п(г,т,т').
0
0
0
5. Характеристики производительности системы
Вычислив стационарное распределение вероятностей состояний системы в произвольный момент времени, находим характеристики производительности:
те
дг = ^р(г, г), г€ |1,1', 2, 2'|,
г=0 те
и = ^ гр(г, г), Г€ |1,1', 2,2'|,
г=0
• среднее число заявок в системе в произвольный момент времени
Ь = Ьхдх + Ьудх> + + Ь^д2>,
• среднее время пребывания заявки в системе (вычисляется по формуле Литтла)
Т = Ь/Л.
6. Численные результаты математического моделирования
В данном разделе приведены численные результаты аналитического моделирования, проведенного с помощью программного комплекса, разработанного автором, а также представлен анализ полученных результатов.
Моделирование гибридного канала, работающего в режиме холодного резерва, проводилось для следующих входных параметров:
Таблица1
Значения параметров математической модели, использовавшихся для получения численных результатов
Параметр Значение
Пропускная способность оптического канала 1024 МЬрз
Пропускная способность радиоканала 512 МЬрз
Время переключения 10 с
Время переключения 1 с
Средняя длина пакета 1024 байта
Средний интервал между пакетами 0.5-2 мс
На рис. 1 представлен график зависимости среднего времени ожидания пакета в очере-Л
жается к средней пропускной способности системы, то задержка в передаче пакетов резко возрастает. В то же время если интенсивность входного потока относительно невелика, то
Л
7. Заключение
В работе проведено комплексное математическое исследование беспроводной гибридной системы связи, включающей лазерный атмосферный канал и резервный радиоканал миллиметрового диапазона радиоволн. Приведены численные результаты, полученные с помощью разработанной модели, а также проведен их анализ. Показано, что адекватной моделью гибридного канала является система массового обслуживания с двумя обслуживающими приборами, обладающими разными скоростями обслуживания. Исследован стохастический процесс, описывающий состояния системы и условия его эргодичности. Разработан
Рис. 1. Зависимость среднего времени ожидания пакета в очереди от Л
а.;п'оритм вычисления стационарных вероятностей состоянии марковскхнх) процесса и характеристик производительности гибридной системы.
Литература
1. Dc.rc.nick 3., Thome С., Spletzer J. On the deployment of a hybrid free-space optic/radio frequency (FSO/RF) mobile ad-hoe network /7 IEEE/RSJ International Conference on Intelligent Robots and Systems. 2005. P. 3990 3996.
2. Вишневский B.M., Кузнецов С., Лаконцев Д., Пом,яков С. Гибридное оборудование на базе радио- и лазерной технологии. Первая миля. http://\v\v\v.moctkom.ru/articlcs/articlc7/articlc7.htm
3. Bloom S., Hartley W.S. The Last-Mile Solution: Hybrid FSO Radio. http:/\/4vww.systemsupportsohitions.(X)m/'whitepapers.html
4. Kim I., Kovervaar E. Availability of Free Space Optics (FSO) and hybrid FSO/RF systems. /7 Optical wireless communication. 2005. V. IV.
5. N adeem F., Geiger В., Le.itge.b E., Awan M.S., Kandus G. Evaluation of Switch over Algorithms for Hybrid FSO-WLAN systems /7 Wireless VITAE. 2009. P. 565 570.
6. Grassman W., Heyman D. Equlibrium distribution of block-structured Markov chains with repeated rows. /7 .Journal of applied probability. 1990. V. 27. P. 57 576.
7. Семенова О. Устойчивый а.;п'оритм расчета етационаржнх) распределения системы обслуживания BMAP/SM/1 с марковским потоком сбоев и двумя режимами работы /7 Автоматика и вычислительная техника. 2004. № 1. С. 75 84.
Поступила в редакцию 05.03.2014-