Имитационное моделирование межсотовых перемещений И отслеживание локальных перегрузок в системах связи с подвижными объектами
Разработана имитационная модель межсотовых перемещений пользователей сети подвижной радиотелефонной связи, позволяющая моделировать поведение мобильного потока абонентов во времени. При этом учитываются наличие двух типов абонентов: находящихся в режиме эстафетной передачи (HO) и абоненты входящие в соту в пассивном состоянии (OR). В качестве среды моделирования используется среда GPSS (General Puipose Simulation System — общецелевая система моделирования) — язык программирования, используемый для имитационного моделирования различных систем, в основном систем массового обслуживания. Одновременной с моделированием и сбором информации о межсо-товых перемещениях абонентов, в данной модели реализованы алгоритмы обработки полученной выборки, с целью предсказания локальных перегрузок в выбранной соте, в режиме реального времени. В качестве примера работы алгоритма приведены результаты моделирования для макросоты Синтезированную программную модель можно использовать на сетях подвижной радиосвязи для предсказания перегрузок и прогнозирования параметров оптимизации. Вычислительная часть полученного алгоритма может быть экспортирована в практически любую программную среду и использована в качестве Ключевые слова: моделирование, GPSS, программного продукта на системе OSS оператора. Интеллектуальные системы данного класса
SON, предсказания, оптимизация. являются частью концепции SON (Self Optimizing Network) разрабатываемой консорциумом 3GPP.
Щучкин В.М.,
аспирант МТУСИ, [email protected]
Введение
Существующие сети подвижной радиосвязи (СПРС) испытывают частые локальные перегрузки вследствие резкого увеличения потребления услуг со стороны абонентов, спрос на которые растет с каждым днем. Величину потребления услуг обычно принято выражать в максимально обслуживаемой нагрузки, при условии обеспечения качества связи (QoS) на заданном уровне. В общем случае, нагрузку можно определить, воспользовавшись соотношением:
Pz=NPy4 (1)
где Ру - суммарная нагрузка [Эрл], N - число абонентов на заданной территории (соте), р^ — удельная нагрузка на абонента [Эрл]. В то же время, одним из основных критериев качества работы СПРС является доля блокировок абонентов, обычно данный показатель не превышает 2-5% от общего числа попыток установления соединения. С другой стороны, в простевшем случае, вероятность блокировки можно определить воспользовавшись формулой Эрганга В, для этого необходимо задаться числом каналов п и нагрузкой л.. Очевидно, что вероятность блокировки зависит только от суммарной нагрузке в соте, в случае лимитированного числа трафиковых каналов РгЬ(Рг) •
Отсюда следует, что увеличение вероятности блокировки может произойти из-за увеличения удельной нагрузки или
вследствие роста числа абонентов на заданной территории (соте). Удельная нагрузка зависит от потребления абонентом услуг связи, например средняя длительность разговора при голосовом вызове. Колебания среднего значения данного параметра возможны, однако их легко предугадать, так как обычно они зависят от активности абонентов в определенные часы суток. В случае изменения удельной нагрузки у небольшого количества абонентов, отклонение суммарной нагрузки от своего среднего значения будут незначительны.
Число абонентов напротив, может претерпевать значительные изменения в короткие интервалы времени, например образование затора на автостраде, проведение массовых мероприятий и т.д., что говорит об актуальности исследования процессов межсотовых перемещений абонентов.
Регулирование нагрузки в соте
Из вышесказанного следует, что поддержание QoS в заданной соте можно осуществлять путем регулирования числа абонентов на обслуживаемой территории. Наиболее привлекательным, с точки зрения простоты реализации и эффективности являются методы регулировки площади при использовании интеллектуальных антенн [1], [5].
Интеллектуальные антенные системы (ИАС) позволяют управлять площадью обслуживания благодаря изменению параметров диаграммы направленности, что является частью концепции SON (Self Optimizing Network) разрабатываемой консорциумом 3GPP.
Наиболее распространенными методами регулировки диаграммы направленности (ДН) являются:
• изменение угла наклона максимума ДН, при этом выделяются случаи электрического и механического регулирования;
• изменение ширины диаграммы направленности, при использовании алгоритмов адаптивного формирования луча.
В работе [ 1 ], был проведен анализ ИАС из которого следует, что использование антенных систем с возможностью удаленного управление углом наклона ДН наиболее эффективно с точки зрения максимально обслуживаемо трафика и простоты реализации. Антенные системы подобного класса получили название RET (Remount Electrical Tilt) в англоязычной литературе [6-8]. Современные RET системы позволяют изменять угол наклона ДН со скоростью до 1 градус на 5-10 секунд.
В [1] была получена аналитическая зависимость, позволяющая определить площадь обслуживаемой соты при различных значениях угла наклона ДН:
5(a), °^Т «>
где М — число шагов отсчета в горизонтальной плоскости, S = 360/M, r(i,а) - функция определяющая дальность связи в направлении / и угле наклона ДН а . Дальность связи это максимально расстояние между РЭС, на котором еще возможно установить соединение.
Если положить, что абоненты распределены по обслуживаемой территории равномерно, то число абонентов на территории соты будет определяться, как
1(a) = S(cc)cj (3)
где а плотность абонентов аб I м2.
ших на интервале времени [О,/„J, К — число запросов на выход абонента из соты, поступивших на интервале времени [0,rj, =;(/„), Е'(/0), Е}(/0), г'р(/0), £>;(О-усредненные моменты входящего и выходящего потоков абонентов [2].
Моделирование межсетевых перемещений абонента
Используя алгоритмы управления площадью обслуживания соты совместно с методами предсказания перегрузок возможно управлять оборудованием СПРС в масштабе реального времени, с целью получения максимального QoS. Для этого необходима разработка имитационной модели, на основе машинного кода, позволяющей в масштабе реального времени оценить параметры адаптации системы.
В данной работе в качестве программной среды моделирования используется GPSS (General Purpose Simulation System — общецелевая система моделирования) — язык программирования, используемый для имитационного моделирования различных систем, в основном систем массового обслуживания.
Блок схема описывающая процесс обслуживания абонента (транзакта) на территории одной соты СПРС (в многоканальном устройстве - МКУ), приведена на рис. 1 и 2. Абоненты условно разделены на абонентов находящихся в режиме эстафетной передачи (НО) и абоненты входящие в соту в пассивном состоянии (OR).
Предсказание локальных перегрузок
Наиболее эффективными, с точки зрения поддержания качества обслуживания на сети, являются методы, позволяющие предотвратить перегрузку в соте, так как это позволит уменьшить коэффициент блокированных пользователей и увеличить обслуживаемых трафик. В работе [1-3], рассмотрены алгоритмы предсказания локальных перегрузок, основанные на контроле межсотовых перемещений абонентов. В данной работе ограничимся рассмотрением алгоритмов обнаружения скачка интенсивности входящего и исходящего потока абонентов, т.к. это наиболее распространенная причина возникновения локальных перегрузок.
Алгоритм обнаружения скачка интенсивности входящего на территорию соты потока абонентов имеет следующий вид:
(Я
",
>
<
П{Р»)
(4)
где I — параметр предполагаемого момента возникновения скачка, и. (О — число запросов, поступивших на интервале времени [о,/], N. - число запросов, поступивших на интервале времени [о, 7"]. Для обеспечения заданной вероятности ложных тревог П выбирается согласно требуемому значению Ра — вероятности ложного обнаружения [2]. В данном
случае в качестве критерия оптимальности используется критерий Неймана-Пирсона.
Для детектирования скачкообразного изменения интенсивности исходящего потока, так же можно воспользоваться соотношением полученным в [2]:
[=#о)-2^гН!(<о)] ".
>ПМ (5)
T-L
'(т-1,Г
где 10 — параметр предполагаемого момента возникновения скачка, К(/0) — число запросов на выход из соты, поступив-
На рис. 1 и 2 приняты следующие обозначения: Р1- Количество сек. до совершения первой попытки звонка, Р2-Длительность разговора, РЗ- Длительность пребывания абонента в соте, Р4- Длительность ожидания канала в случае блокировки, М1- Длительность пребывания абонента в системе. Параметры Р1-Р4 определяются моделью автоматически для каждого транзакта отдельно, на основе псевдослучайных счетчиков.
Данная блок-схема несколько отличается от схемы, отражающей процессы обработки транзакта в среде ОР85. Данные упрощения приведены улучшения восприятия материала и обеспечения прозрачности работы модели. Модель хорошо согласуется с теоретическими формулами, которые описывают соту СПРС как СМО, например Эрланга В.
/"’ТйжбывТтье N сот» Р1_________У
н*т
/^ЗЛНЯТЪ АЛИЧ~""'\
У^. щрг ^
Рис. 2. Модель обслуживания абонента (Ж в соте
Для отслеживания локальных перегрузок в синтезированную модель соты введены два дополнительных блока:
• блок сбора и накопления статистики межсотовых перемещений абонента (позволяет накапливать данные о времени входа и выхода абонента из соты);
• блок обработки статистики (математический аппарат, который используя накопленную статистику, определяет возможность наступления перегрузки и оценивает среднее значение параметров мобильности абонента).
Блок сбора и накопления статистики реализуется при помощи стандартных средств ОРРБ. В данном блоке фиксируются и хранятся основные параметры мобильности абонента: время входа абонента на территорию соты, время выхода абонента с территории соты.
Блок обработки статистики извлекает необходимые данные о межсотовых перемещениях абонентов и прогнозирует наступление перегрузки, а так же оценку параметров мобильного потока. Основные параметры для оценки межсотовых перемещений абонентов являются, для входящего потока: 1а — оценка момента возникновения скачка. А* и —
оценка интенсивности входящего потока до и после момента наступления скачка мобильности абонентов. Для исходящего потока: 10 - оценка момента возникновения скачка, и
“ — оценка интенсивности исходящего потока до и после момента наступления скачка мобильности абонентов. Ниже приведены результаты имитационного моделирования детектирования скачка входящего/исходящего потоков, а так же произведем оценку параметров потоков.
В качестве исходных данных для создания модели соты и оценки мобильности абонентов воспользуемся данными, приведенными в табл. 1.
Данные значения использовались для определения параметров транзакта в системе вР88, определяя, таким образом, средний интервал времени между двумя последовательными входами и выходами абонентов на/с территорию/и соты. Реальные значения параметров, определялись при помощи генераторов псевдослучайных чисел встроенных в ОРБ8 по экспоненциальному закону. В качестве параметра, позво-
ляющего менять интенсивность поступления абонентов в соту (транзактов в многоканальное устройство - в терминологии ОРБЭ), использовался параметр а — плотность распределения абонентов.
Таблица 1
Исходные данные
Параметр Значе- ние Размер- ность
Процент абонентов совершающих handover от общего числа абонентов 1 %
Интенсивность входа в соту абонентов, X* 0,06724 с
Интенсивность выхода абонетх» из сопы, ц* 121 с
Удельная нагрузка на абонента 0.015 Эрл
* - параметры рассчитаны из условия движения абонентов 50 к.и/ч. площадь соты 3600000 м:, плотность абонентов вне соты 0.0005 аб/ м1. Указания по расчету параметров приведены в [2].
Условия проведения эксперимента для входящего потока абонентов:
• Эксперимент проводится на Зх различных интервалах анализа То: 100с/300с/500с;
• Детерминированным временем появления скачка интенсивности на выбранных интервалах будет временная отметка 70с/200с/400с (детерминированное время скачка необходимо для оценки параметров алгоритма обнаружения);
• Интенсивность входящего потока меняется путем увеличения/уменьшения относительной плотности абонентов оскач/о. Рассматриваемые значения относительной плотности: 0,5; 0,6; 0,7; 0,8; 0,9; 1; 1,1; 1.2; 1,3; 1.4; 1,5.
Результаты моделирования скачка интенсивности входящего потока абонентов и оценка его параметров приведены на рис. 3-5.
'4нисм1с порогаобкфудсмшгАкнснп интенсивное т входящегогн>ток*
ОД 0л о fi О,в ! |Д 14
OIHOC1ПСЛЫККI|>СДЖС «наЧОЯ 1C плотное* I .Гюнсню» ОСК1Ч 3
■ Тл-irwv То-300< -•-То-50<к
Рис. 3. Зависимость порога обнаружения изменения интенсивности входящего потока при различных значениях относительной плотности абонентов
Отн(ХИ1сльнэсотклонение оценки момем'д сьачка
£ Ш
Н I»)
£ ол» с
•• V’C'4>’*-
* ■* ■ ■
—Тс» 100с —Тс-300с “Тс—500с
отмоапсльнос е|:еднсс яичлк плотное иI кнчктов оекач о
Рис. 4. Относительное отклонение оценки момента скачка входящего потока при различных значениях относительной плотности абонентов
1 ^тосшгаьнае отклонение <>цгни1 интгмавмогт входящггопггок»
-ІО-ЮОС
“То^ЗООС
-Т«-ЗООе
отнштслыкке рсдисс шопе гмотжчти лСн.налов лклч о
Рис. 5. Относительное отклонение оценки интенсивности скачка входящего потока при различных значениях относительной плотности абонентов
Из полученных результатов можно сделать вывод, что максимальное отклонение от действительного значения параметра скачка интенсивности наблюдается в моменты отсутствия самого скачка. Данный эффект не влияет на работоспособность алгоритма при правильном выборе порога обнаружения изменения интенсивности, что подтверждается результатами, см. рис. 3. При отсутствии изменения интенсив-
У
ности входящего потока оценка параметров не производится, при этом малые колебания числа абонентов в соте, которые не будут отслеживаться, незначительно скажутся на общее качество обслуживания.
Па рис. 6-8, для полной оценки работы алгоритма, показан результат ряда дополнительных экспериментов. Для этого зафиксировали величину скачка интенсивности (1,3) и отследим ошибки оценки параметров при различных моментах времени возникновения изменения интенсивности. Результаты моделирования представлены на рис. 6-7.
Отоягтслыюсспоюнсмк оиснкнмлмош скачка при рагшгтых моментах I именем !Я П.1| лмстр о В М > *>1 П ЬНОС ПI абонентов
“|Т<к- Тск |Тск
Рис. 6. Относительное отклонение оценки момента скачка, при различных моментах изменения параметров мобильности абонентов
Отнои пельносотклсжн! к оценки ютсмснвносш скачка. р.тмчмычмоментах
I именем н параметров мобтьносиI абонентов
* О 1
-ц.<х>кх
Тс
Рис. 7. Относительное отклонение оценки интенсивности скачка, при различных моментах изменения параметров мобильности абонентов
Дополнительно приведем полученную в результате эксперимента зависимость значение порога обнаружения изменения интенсивности исходящего потока, при различных моментах изменения параметров мобильности абонентов
Чюпение П&рОП».• бшружвпніпменеіпи твгтенпшюсп! ікхлдщеглПОТОКА гри
І».тліпнкг:моментахі пмснснпя параметров моЛіиьносшябснентов
Рис. 8. Значение порога обнаружения изменения интенсивности исходящего потока, при различных моментах изменения параметров мобильности абонентов
На рис. 6-8 можно наблюдать наличие ошибки в случае возникновения изменения интенсивности на границе интервала наблюдения. Для устранения данного эффекта, проявляющегося на границах интервала наблюдения необходимо использовать метод параллельных вычислений с перекрывающимися интервалами наблюдений, технически пригодный вариант которых приведен на рис. 9.
Рабочая область интервала для формирования опенки параметра момента /0
Рис. 9. Организация разделения времени на интервалы анализа
Для детектирования скачка исходящего потока абонентов определим следующие начальные условия:
• Эксперимент проводится на трех различных интервалах анализа То: 100с/300с/500 с;
• Детерминированным временем появления скачка интенсивности на выбранных интервалах будет временная отметка 30с/150с/250с (детерминированное время скачка необходимо для оценки параметров алгоритма обнаружения);
• Интенсивность входящего потока меняется путем увеличения/уменьшения относительного пребывания абонента в соте Ат. Рассматриваемые значения: 0,5; 0,6; 0,7; 0,8; 0,9; 1; 1,1; 1,2; 1,3; 1,4; 1,5;
• Интенсивность входящего потока на интервале постоянна.
Результаты моделирования скачка интенсивности исходящего потока абонентов и оценка его параметров приведены на рис. 10-15.
*ОЧ€ИК порог* ОСНП|>уЛ.(Н1Я!ГМеНеН11Я ЮТСНСИВНОСЛ! исходящего потока
-*-То=100< -*-То=300с -А-То=500с
О 0,2 04 06 04 1 и 1.4 1,6
отоагтеяыюе еднеетндчемя дигтетьносшпреоыв.тння ложенгз б соте Л;
Рис. 10. Зависимость значение порога обнаружения изменения интенсивности исходящего потока от Ат для различных интервалов наблюдения
Стхігтмьное .іікілкнк оценки помет мчи
їло К шо ?омо Ьомо У оюо £0.200 “о 000
02 0.4 )Л 0 Я
и
Тс—100с Тс-300с
Те=500с
отхкгеинм среднее значеим дтпельнесипфебимия абонента всоте Д1
Рис. 11. Зависимость относительного отклонения оценки момента скачка от Ат для различных интервалов наблюдения
Относ] пгяшое отклоняя |р одоыияпснснвностн исходящего потока до момента шменени параметров ‘цэержкн абонентов в сок
10.*»
Г ММ
2 6 <к>) Г 4.00)
2 2.00)
— ООО)
г
А й А— .
0.2 04 Об (.8 1 и 14 16
оі№чіпсіьиосі|і(цмсе тінта д иисіин* иіі4>счік.ми .пЧчммііі щоїс Лі
-Т(-100с
-Тс“}00с
-Тс=500с
Рис. 12. Зависимость относительного отклонения оценки интенсивности исходящего потока, до момента изменения параметров задержки абонентов в соте от Дт для различных интервалов наблюдения
Полученный результат аналогичен предыдущему случаю. Однако необходимо отметить еще одно важное свойство: так как точность результатов зависит от величины выборки [2], следовательно она зависит и от интервала наблюдения То, что подтверждается результатами эксперимента. Например, в случае интервала наблюдения 500 с, ошибка детектирования скачка и оценки параметров межсотовых перемещений абонентов меньше, чем при То = 100, 300 с, см. рис. 4-7, 11-12, при идентичных начальных условиях.
Выводы
1. Полученные результаты подтверждают целесообразность использования синтезированной модели при оптимизации сети СПРС. При этом интервал наблюдения, время накопления статистики, необходимо выбирать в зависимости от необходимой точности. В рассматриваемом случае площадь соты 3,6км2, средняя скорость движения абонентов V = 50лсу// •/, плотность абонентов на обслуживаемой терри-
тории а = 0,0005аб/лГ интервал наблюдения целесообразно выбирать не менее чем 500 с., при этом интервал анализа, на котором происходит формирование оценки, целесообразно выбирать не менее 100 с., что позволит достаточно точно детектировать наличие скачка и иметь запас времени на применение алгоритмов оптимизации.
2. Определена возможность использования алгоритма оптимизации покрытия соты при помощи ИЕТ систем, при условии правильного детектирования перегрузки и корректной оценки параметров. Например, если вследствие оценки параметров определена возможность перегрузки и оценка интенсивности входящего потока = кА*, то учитывая, что Я“ = 2/? • а ■ V или при пересчете /? в эквивалентную площадь Xм = 2^1 Б / п ■ а ■ V (см. [2]). Очевидно, что площадь соты нужно уменьшить в к ~ раз. При этом угол наклона выбирается в соответствии с диаграммой направленности антенны БС. Аналогично для случая со скачком интенсивности исходящего потока.
3. Синтезированную программную модель можно использовать на сетях СПРС для предсказания перегрузок и прогнозирования параметров оптимизации. Вычислительная часть полученного алгоритма может быть экспортирована в
практически любую программную среду и использована в качестве программного продукта на системе OSS оператора.
Литература
1. Шорнн О.А., Щучкин В.М. Использование интеллектуальных антенн в системах мобильной связи для снижения перегрузок//Труды МАИ, 2012. -№53.
2. Шорин О. А. Методы оптимального распределения частотно-временного ресурса в системах подвижной радиосвязи: Дис. ... д-ра. техн. наук / МТУСИ. - М., 2005.-351 с.
3. Демьянов А.И. Оценка параметров скачков нагрузки в сотовых сетях подвижной связи // Электросвязь, 2002. - №1. -С.33-36.
4. Джейкс У.К. Связь с подвижными объектами в диапазоне СВЧ. - М.: Связь, 1979. - 520 с.
5. Алехин Ю.Н. Дистанционное управление положением луча в панельных антеннах КАТНКЕ1Ы//Мобильные системы. — 2003. -№3. -С.62-66.
6. http://www.comba-telecom.com/Prd_Channel.aspx?ClasslD= 21 // 15.05.2012.
7. http://www.amphenol-jaybeam.com/remote-electrical-tilt-technology.php// 15.05.2012.
8. 3GPP TS 25.463 UTRAN Iuant interface: Remote Electrical Tilting (RET) antennas Application Part (RETAP) signaling.
Simulation modeling of inter-cell movement and tracking of local congestion in communication systems with mobile objects
Shcuchkin Vladimir M., postgraduate student of Moscow Technical University of Communications and Informatics.
Abstract
The simulation model of inter-cell movements of users of mobile radio communication network is developed in this paper. That allows one to simulate the behavior of mobile subscribers in the real-time. Presence of two types of subscribers is taken into account: being in the mode of handover (HO), and subscribers in the passive state (OR). A programming language GPSS (General Purpose Simulation System — general-purpose modeling system) is used as the simulation environment for the simulation of various systems, primarily queuing systems. Simultaneously with the modeling and the collection of information about intra-cell movements of subscribers in this model are implemented algorithms of calculation of obtained samples in order to predict local overloads in the selected cell in real time. As an example of the algorithm we present simulation results for macro-cell. Synthesized programming model can be used in mobile radio networks for the prediction of congestion and optimization parameters. Computation part of the resulting algorithm can be exported to virtually any programming environment and used as a software product in the operators OSS system. Intelligent systems of this class are part of the SON (Self Optimizing Network) concept developed by a consortium of 3GPP
References
1. Shorin O. A., Shchuchkin V M. Use of intellectual aerials in systems of mobile communication for decrease in overloads // Works MAI, 2012. — No. 53.
2. Shorin O. A. Methods of optimum distribution of a time-and-frequency resource in systems of a mobile radio communication: Dis. ... Drs. tehn. sciences / MTUCI. — M., 2005. — 351 p.
3. Demyanov A.I. An assessment of parameters of jumps of loading in cellular networks of mobile communication // Telecommunication, 2002. — No.1. — pp. 33-36.
4. Jeiks U.K. Communication with mobile objects in a range of the microwave oven. — M.: Svyaz, 1979. — 520 p.
5. Alekhin Yu.N. Remote control by the provision of a beam in panel KATHREIN aerials // Mobile systems, 2003. — No. 3. — pp. 62-66.
6. http://www.comba-telecom.com/Prd_Channel.aspx?ClassID=21 // 15.05.2012.
7. http://www.amphenol-jaybeam.com/remote-electrical-tilt-technology.php // 15.05.2012.
8. 3GPP TS 25.463 UTRAN Iuant interface: Remote Electrical Tilting (RET) antennas Application Part (RETAP) signaling.