Нелинейная задача расчета зон обслуживания базовых станций

Тип работы:
Реферат
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 621. 396. 41
Л. Л. Егоров, В.А. Кологривов
Нелинейная задача расчета зон обслуживания базовых станций
Сформулирована нелинейная задача расчета зон покрытия базовых станций сотовой сети стандарта GSM и на основе метода линеаризации нелинейных систем алгебраических уравнений Ньютона-Рафсона предложен нелинейный алгоритм метода наименьших квадратов для ее решения. Приведена оценка эффективности предложенного метода при решении задач планирования и оптимизации сети.
Ключевые слова: метод наименьших квадратов, нелинейная задача, итерационный алгоритм, базовая станция, зона покрытия.
Введение
Ранее [1, 2] были предложены и исследованы алгоритмы и методики расчетов оптимальных зон покрытия кластеров базовых станций сотовых сетей связи. Суть этих методов сводится к решению переопределенной системы линейных алгебраических уравнений вида A • r = d методом наименьших квадратов (МНК) (1) и методом взвешенных наименьших квадратов (МВНК) (2)
(-1
r =
AT • A
• AT • d, (1)
г = А • А]& quot-1 • АТ • d, (2)
где, А = [ау ^ - матрица эластичности зон покрытия, определяемая интенсивностями нагрузок соседних БС уг и yj — г = |г-| - вектор-столбец неизвестных, соответствующих радиусам круговых зон покрытия- d = - вектор известных расстояний пролетов между
соседними базовыми станциями (БС).
В обоих случаях искомыми величинами данных расчетов являются радиусы зон обслуживания г (у1, yj, ^), функционально зависящие от предполагаемых нагрузок БС
У1, yj и расстояний пролетов между БС. Важно отметить, что при решении данных задач средние интенсивности нагрузок на БС уг, yj были величинами постоянными и независимыми от размеров зон обслуживания. Полученные при этом результаты не противоречат физическим и геометрическим представлениям.
Однако более реалистично предполагать зависимость нагрузок базовых станций, а следовательно, и коэффициентов матрицы эластичности от размеров зон обслуживания. Такая задача расчета оптимальных размеров зон покрытия является нелинейной. При этом предполагается задание двумерной плотности абонентов на обслуживаемой территории, однако здесь для простоты изложения алгоритма решения нелинейной задачи плотность абонентов р предполагается постоянной.
Интенсивности нагрузок БС при этом пропорциональны площадям или квадратам радиусов зон покрытия:
2
Уг =р-^ =Р-я-г, (3)
где р — плотность абонентов- Si — площадь зоны покрытия г-й БС- г — радиус г-й зоны покрытия. Соответственно, элементы матрицы эластичности при этом приобретают вид
^ = k = k (4)
Уг + Уj Г2 + Г
где k — коэффициент, учитывающий влияние мощности передатчика БС, высоты подвеса антенны и т. д.
Решение нелинейной задачи предполагает задание начальных значений, в данном случае — интенсивностей нагрузок БС. В качестве начальных значений интенсивностей нагрузок БС можно использовать данные, полученные на этапе планирования сети, либо данные с контроллеров БС, принадлежащих интересующему нас кластеру (см. рис. 1), полученные в процессе эксплуатации сети предшествующей этапу оптимизации.
БС1
сП2
БС2
сПА
tfl3
о23
/БСЗ
о24
& lt-з35
(ЛЬ
БС4
045
БС5
Рис. 1. Симметричный кластер (группа БС) размерностью д = 5
В линейном случае данные о нагрузках БС уь однозначно определяют элементы матрицы эластичности, А = [ау ], что при известных расстояниях
пролетов между БС dij с использованием соотношений (1) и (2) позволяет рассчитать оптимальные размеры г = |г-| зон обслуживания БС.
В нелинейном случае задача решается итерационными методами. Здесь на каждом шаге итерации вычисляются размеры зон обслуживания и с использованием заданной плотности абонентов пере-считываются интенсивности нагрузок БС. Условием завершения итерационного процесса является уравновешивание размеров зон обслуживания и нагрузок БС, при котором значимое изменение размеров зон от итерации к итерации прекращается.
Целями данной статьи являются изложение нелинейного алгоритма МНК (НМНК), нахождение зон покрытий БС путем решения переопределенной системы нелинейных алгебраических уравнений, анализ полученных вариантов решений.
Решение переопределенной нелинейной системы
Решение переопределенных систем нелинейных уравнений может быть реализовано различными методами. Наиболее известным является алгоритм Ньютона-Рафсона, основанный на линеаризации исходной нелинейной системы уравнений [3] путем разложения нелинейной вектор-функции в окрестности предполагаемого корня в ряд Тейлора и удержании линейных членов ряда. В результате приходим к линеаризованной системе алгебраических уравнений
J (X (k)) • ДХ (к) = ^(Х (к)), (5)
где F'-(X (k)) = J (X (k)) — матрица первых производных- ДХ (к) = Х (кь1) — Х (к) — вектор приращений неизвестных или разностный вектор решений на (к+1)-й и к-й итерациях.
Для переопределенной исходной нелинейной системы переопределенная линеаризованная система алгебраических уравнений может быть решена методом наименьших квадратов (МНК) или (МВНК)
-1
ДХ (к) = - JT (Х (к)) •Х (к)) • JT (Х (к)) •Х (к)).
(6)
В результате итерационная запись метода Ньютона-Рафсона или нелинейного алгоритма МНК приобретает вид
х (й+1) = Х (к) — JT (Х (к)) •Х (к)) • JT (Х (к)) •Х (к)).
-1
(7)
Сходимость метода Ньютона-Рафсона существенно зависит от начальных значений
Х (0), и для обеспечения сходимости и снижения требуемого числа итераций обычно модифицируют выражение (7)
х (й+1) = Х (к) — ц • ^ (Х (к)) •Х (к)) • JT (Х (к)) •Х (к)),
-1
(8)
где ti — дополнительный коэффициент, подбираемый на каждом шаге итерации.
В нашем случае система нелинейных алгебраических уравнений и нелинейная вектор-функция имеют вид
АЩ) • R = D, (9)
F® = АЩ) • R — D = 0, (10)
где R — вектор неизвестных (радиусов зон обслуживания) — АЩ) — матрица эластичности, элементы которой являются функциями радиусов- В — вектор расстояний пролетов БС. В покомпонентном виде нелинейное уравнение т-го пролета между Ь-й и у-й БС как компонент нелинейной вектор-функции имеет вид
т (г г)=к •
г2 + г2
• Г + к •
Г2 + Г2
• г- ^
(11)
2
г
г
соответственно, т-я строка матрицы первых производных содержит два отличных от нуля компонента:
J[m,•••] =
k • (г? • Г + 3 • г2)& quot- 2* п * г3) k •(г2 • (г2 + 3 • п2) — 2* п, * г3)
Г + г, 2)2
(П2 + п, 2)2
(12)
БС1
?/12
БС2
Таким образом, нелинейная задача расчета зон покрытия сети сотовой связи может быть решена нелинейным методом наименьших квадратов, основанным на линеаризации нелинейных систем уравнений методом Ньютона-Рафсона.
Постановка задачи исследования
Решение нелинейной задачи будем находить, используя Ма^аЬ [4]. Расчет зон покрытий произведем на примере симметричного (рис. 1) и несимметричного (рис. 2) кластеров с равномерной нагрузкой на БС.
Для симметричного кластера расстояния между БС равны ?12 = ?25 = ?45 = ?14 = 10 км, ?13 = ?23 = = ?35 = ?34 = 7,071км- для несимметричного — ?12 = =?14 =10км- ?13 = ?23 = ?34 = 7,07 км- ?25 = ?45 = 12,61км- ?35 = 9,89 км. Далее сравним результаты решения нелинейной задачи нахождения зон покрытий БС с аналогичными вычислениями, полученными ранее в [2, 3] по методу наименьших квадратов (МНК) и модифицированному методу взвешенных наименьших квадратов (ММВНК).
Количественную и качественную оценку произведем по значению величины «нормы» Д [1], которая показывает суммарную среднеквадратическую ошибку решения. Слагаемые нормы, то есть разности суммы радиусов зон обслуживания пары БС и расстояния между этой парой, пропорциональны площадям перекрытия/недопокрытия:
с!14
БС5
Д = л/& gt-((п + п,) — ?г,)2, (13)
где Д — суммарная среднеквадратическая ошибка решения- п и г, — радиусы зон обслуживания соответствующих БСг и БС, — ?г, — пролет между БСг и БС,.
Варианты решений
При равномерной плотности абонентов на местности для симметричного кластера исследованы различные варианты решений. Наилучший результат вычислений по величине наименьшей среднеквадратической ошибки приведен в сравнительной табл. 1.
Рис. 2. Несимметричный кластер (группа БС) размерностью д=5
Оценка зон покрытий БС для симметричного кластера
Таблица 1
Метод вычисления Величина Д по МНК Величина Д по ММВНК Величина Д по методу Ньютона-Рафсона
Д Радиусы Б С, км Д Радиусы Б С, км Д Радиусы Б С, км
Нагрузка на БС У1−5 = 1 Эрл 5,4 • 10−15 Г1=5 Г2=5 г3 =2,071 г4 =5 Г5=5 5,4 • 10−15 Г1=5 Г2=5 г3 =2,071 г4 =5 Г5=5 0,0259 г1 =4,9523 г2 =4,9523 г3 =2,3518 г4 =4,9523 г5 =4,9523
Результаты вычислений для несимметричного кластера при тех же исходных данных приведены в табл. 2.
Таблица 2
Оценка зон покрытий БС для несимметричного кластера
Метод вычисления МНК ММВНК Метод Ньютона-Рафсона
Д Радиусы Б С, км Д Радиусы Б С, км Д Радиусы Б С, км
Нагрузка на БС У1−5 = 1 Эрл 0,645 7 г1 =5,0665 г2 =4,9311 г3 =2,2706 г4 =4,6699 г5 =8,0963 0,6457 г1 =5,0665 г2 =4,9311 г3 =2,2706 г4 =4,6699 г5 =8,0963 0,074 2 г1 = 4,935 г2 = 4,8945 г3 = 2,3433 г4 = 4,8814 г5 = 7,9385
Графическая интерпретация наилучших вариантов, полученных по методу Ньютона-Рафсона, приведена на рис. 3 и 4.
размерностью д = 5 при У15 = 1 Эрл размерностью д = 5 при У15 = 1 Эрл
Анализ полученных результатов
Результаты вычислений радиусов зон обслуживания БС по нелинейному алгоритму метода наименьших квадратов не противоречат физическим представлениям и адекватно учитывают геометрическую конфигурацию кластеров при заданной плотности нагрузки. В рассмотренных вариантах получены результаты, близкие к ранее достигнутым методами МВНК и МНК.
Решение нелинейной системы находится малым количеством итераций. Различным конфигурациям кластеров БС соответствуют свои неповторяющиеся решения.
Наблюдается зависимость решения от начального приближения, что приводит к различным вариантам покрытия кластера с одной и той же конфигурацией.
Результаты решения нелинейной задачи предложенным методом характеризуются большей среднеквадратической ошибкой по сравнению с линейной задачей в случае симметричного кластера, что, вероятно, обусловлено предварительной линеаризацией исходной нелинейной системы уравнений, и меньшей среднеквадратической ошибкой в случае несимметричного кластера, что обусловлено возможностью варьирования и подбора начальных условий в исходной системе уравнений.
Целесообразным является рассмотрение альтернативного метода решения переопределенной нелинейной системы уравнений, первый этап которого включал бы МНК, а на втором этапе линеаризацию.
Адекватное использование предложенного алгоритма предполагает задание двумерной плотности абонентов на обслуживаемой территории.
Литература
1. Егоров Л. Л. Алгоритм расчета зон покрытия базовых станций сотовой связи / Л. Л. Егоров, В. А. Кологривов // Доклады ТУСУРа (Томск). — 2007. — 2(16). -С. 157−162.
2. Егоров Л. Л. Алгоритм расчета зон покрытия базовых станций сотовой связи / Л. Л. Егоров, В. А. Кологривов, С. В. Мелихов // Доклады ТУСУРа (Томск). — 2009. -5(19). — С. 15−20.
3. Влах И. Машинные методы анализа и проектирования электронных схем / И. Влах, К. Сингхал / Пер. с англ. — М.: Радио и связь, 1988. — 560 с.
4. Потемкин В. П. Справочник по применению системы РС Ма^аЬ. — М.: Физматлит, 1993. — 112 с.
Егоров Леонид Леонидович
Аспирант каф. средств радиосвязи ТУСУРа Тел.: 41−37−09
Эл. почта: Yegoroff@sibmail. com Кологривов Василий Андреевич
Канд. техн. наук, доцент каф. средств радиосвязи ТУСУРа Тел.: 41−37−09
Эл. почта: mrc@main. tusur. ru
L.L. Yegorov, V.A. Kologrivov
Nonlinear problem of the calculation of a base station'-s coverage area
The nonlinear problem of calculation the GSM-standard base stations'- coverage area based on the linearization of Newton-Raphson nonlinear systems of algebraic equations has been stated. The nonlinear least squares method algorithm of its solving has been offered. The estimation of the efficiency of an offered method in solving problems with network planning and optimization has been made.
Keywords: least squares method, nonlinear problem, iterative algorithm, base station, coverage zone.

ПоказатьСвернуть
Заполнить форму текущей работой