CC BY
76
Труды Карельского научного центра РАН № 4. 2014. С. 29-35
УДК 519.25
АВТОМАТИЧЕСКАЯ КАЛИБРОВКА ЛОКАЛЬНЫХ СИСТЕМ ПОЗИЦИОНИРОВАНИЯ НА ОСНОВЕ ПОСТРОЕНИЯ КАРТЫ СИЛ СИГНАЛОВ
Р. В. Воронов1, О. В. Лукашенко2, А. П. Мощевикин1,3
1 Петрозаводский государственный университет
2Институт прикладных математических исследований Карельского научного центра РАН 3ЗАО «РТЛ-Сервис»
Для определения местоположения объекта в закрытых помещениях часто используют локальные системы позиционирования, основанные на методе построения шаблонов карт слышимости от нескольких базовых радиостанций (fingerprinting technique). Для создания этих шаблонов необходимо провести ряд замеров силы сигнала во многих точках здания. Изменение окружающей обстановки (перестановка мебели, движущиеся объекты, изменение положения антенн и пр.) значительно ухудшает точность локации. В статье предложен способ автоматической рекалибрации системы расчетов локации при изменении внешних условий. Таким образом, для учета изменений нет необходимости проводить повторную процедуру обучения системы.
Ключевые слова: локальные системы позиционирования, локация, базовая станция, сила сигнала.
R. V. Voronov, О. V. Lukashenko, А. P. Moschevikin. AUTOMATIC RECALIBRATION OF A LOCAL POSITIONING SYSTEM BASED ON THE FINGERPRINTING TECHNIQUE
Local positioning systems based on the fingerprinting technique are widely used for localization indoors. This technique requires preliminary measurements of the received radio signal strength in many points inside a building. The changes in surroundings (rearrangement of furniture, changing antenna orientation, blocking of line-of-sight by large moving objects, etc.) substantially decrease the localization accuracy. This paper discusses the proposed method of the automatic recalibration of the position calculation system in case of changes in the environment. Thus it is not necessary to carry out additional measurements in order to consider these changes.
Key words: local positioning systems, location, base station, signal strength.
Введение
Известно, что спутниковые системы навигации невозможно использовать в закрытых пространствах: в зданиях, шахтах и т. д. В таких случаях разрабатываются специальные
системы локации на основе беспроводных сетей датчиков. При этом для определения местоположения мобильных объектов используют различные методы, в том числе основанные на определении мощности входного сигна-
0
л а [4, 5]. В помещении устанавливаются стационарные базовые станции сети, каждая из которых имеет свою зону покрытия. Для создания подобных систем используют беспроводные технологии локальных сетей WiFi, ZigBee, nanoLOC, UWB, Bluetooth и др.
Под локацией понимают определение местоположения каких-либо объектов. Необходимым условием для расчета локации является предварительная настройка устанавливаемой системы — разбиение пространства помещения на зоны и определение для каждой из них сил принимаемых сигналов от базовых станций. Обычно зона не превышает размеров комнаты (иногда в комнатах можно выделить несколько зон). Для каждой зоны заранее формируется набор средних значений сил сигнала от базовых станций (шаблоны слышимости). Когда мобильный объект в очередной раз выходит в эфир, его сигнал регистрируется всеми базовыми станциями в пределах слышимости. Этот набор сил сигналов сравнивается с заранее полученными шаблонами. В качестве искомой зоны локации выбирается зона с тем шаблоном, который наилучшим образом совпадает с зарегистрированным набором сил сигналов [4-6].
Необходимо отметить, что сила сигнала на приемной стороне является дискретной величиной, поскольку радиопередача осуществляется в цифровой форме, и входной сигнал регистрируется, например, с использованием аналого-цифровых преобразователей или управляемых автоматических регуляторов усиления.
Пусть N — число точек (или зон) пространства, в которых может находиться мобильный объект, D — число базовых станций. Множество сил сигналов от базовых станций является дискретным. Тогда в ходе настройки системы необходимо для каждой базовой станции d определить вероятности Pd(x\j) того, что будет зафиксирована сила сигнала х при условии нахождения мобильного объекта в зоне j = 1,... ,N.
В ходе эксплуатации системы местоположение мобильного объекта определяется следующим образом. Пусть pj — априорная вероятность нахождения объекта в точке j (j =
1,..., N). Предположим, что зарегистрирован вектор сил сигнала х = (х\,... ,хв), где -сила сигнала базовой станции d. Тогда апостериорная вероятность qj нахождения объекта в точке j равна
в
Рз П РлЫи)
* = 1Г^Б------------‘ Ш
X) Рк П Ра{ха\к) к=1 <1=1
Поскольку сила сигнала на приемной стороне сильно зависит от числа отражений, величины и формы препятствий по пути распространения электромагнитной волны, любые изменения в помещении (перестановка мебели, изменение ориентации антенн базовых станций и т. д.) приводят к необходимости заново настраивать систему.
Кроме этого, в ходе инсталляции системы позиционирования всегда стараются уменьшить количество используемых базовых станций с целью минимизации стоимости оборудования. Поэтому зачастую используют следующий итерационный сценарий. Сначала устанавливают «узловые» базовые станции в расчетных зонах здания, обучают систему, строят шаблоны и карты слышимости, а затем по мере необходимости в разных частях здания добавляют еще несколько станций, улучшая точность локации в этих зонах. Поскольку процесс обучения системы занимает достаточно много времени, в некоторых случаях его можно заменить автоматической процедурой формирования новых шаблонов сил сигнала в новых зонах или от ранее не использованных базовых станций. Понятно, что эффективность автоматической конфигурации зависит от плотности расположения базовых станций, т. е. от фактора пересечения зон радиопокрытия разных базовых станций.
В данной статье предлагается подход, позволяющий произвести автоматическую перенастройку системы. Подход основан на оценке неизвестных параметров распределения мощности сигнала в ходе эксплуатации системы. Такая перенастройка может быть осуществлена в случаях:
• добавления новых базовых станций;
• изменения конфигурации помещения;
• добавления новых зон, которые ранее не были нанесены на карту;
• изменения условий работы базовой станции (наклон антенны и прочее).
0
Модель потерь мощности сигнала в эфире
В большинстве работ [7-9], в которых рассматривается процесс распространения электромагнитных волн для диапазона частот выше 1 ГГц, приводится следующая модель затухания:
PL(s)
PL0 + IO7 lg — s0
+ C) s > so, (2)
где PL(s) - потеря мощности сигнала между передатчиком и приемником (path loss), выраженная в децибелах. Слагаемое в квадратных скобках описывает детерминированную составляющую потерь, зависящую от расстояния s между радиоустройствами, а С - случайная компонента. PLq — это потеря мощности на заранее известном расстоянии sо, например, в 1 м. Характер зависимости уменьшения мощности сигнала с увеличением расстояния зависит от значения показателя 7, связанного с особенностями окружения. Например, для офисных помещений 7 может изменяться от значения 1,7 для условий прямой видимости (LOS, Line-of-Sight) между передатчиком и приемником до значений 3,5 и выше для условий отсутствия прямой видимости (NLOS, Non-Line-of-Sight) [7]. Некоторые исследователи вводят в модель параметры, отвечающие за число стен, встречающихся на пути распространения сигнала, их толщину и материал [11].
Вторая компонента в формуле (2) описывает изменчивость локальных условий на стороне приемника радиосигнала. Она включает в себя учет множественных путей распространения электромагнитных волн, огибания ими препятствия, интерференционного сложения фаз и амплитуд. £ зависит от конфигурации помещений и расположения объектов и помех на пути распространения волны, ориентации и диаграммы направленности антенны, размера антенны, степени поляризации излучения и др. Компонента £ может принимать как положительные, так и отрицательные значения. Часто предполагается, что £ - случайная величина, имеющая нормальное распределение N(0, а) с нулевым математическим ожиданием и некоторым стандартным отклонением а [7-9].
Содержательно формула (2) показывает, что потеря мощности сигнала (а значит, и сама мощность сигнала на принимающей стороне) может моделироваться как нормальная случайная величина N(a,a), где параметры
а, а > 0.
На рисунке в качестве иллюстрации величины потерь мощности сигнала внутри помещений приведена карта слышимости одной точки доступа в технологии RealTrac [2, 3, 10]. Числа в квадратах обозначают потери мощности сигнала, выраженные в децибелах. Местоположение стационарного устройства А590 обозначено темным квадратом с перекрестием внутри. Мощность выходного сигнала точки доступа равнялась 20 дБм (100 мВт). Числа внутри ячеек метровой длины соответствуют средним потерям сигнала в децибелах. Например, для числа 64 (желтая ячейка), соответствующего уровню потерь в эфире PL = 64 дБм, зарегистрированная сила входного сигнала на мобильном радиоузле, находящемся в данном квадрате, равнялась —44 дБм.
Приведенная карта радиопокрытия также подтверждает возможность расчета локации на основе силы входного сигнала. Если зарегистрированное значение потерь в эфире от точки доступа А590 оказалось, скажем, около 80 дБм, то набор вероятных мест нахождения объекта резко уменьшается (это либо небольшая часть длинного коридора, либо одна из комнат). Если на одном этаже установлено 3-5 таких точек доступа, то достаточно легко добиться точности локации на уровне комнаты (room level).
Метод автоматической калибровки
СИСТЕМЫ ЛОКАЦИИ
Рассмотрим случай, когда имеется некоторое множество базовых станций, для которых известны распределения вероятностей сил сигнала в каждой зоне пространства. Тем самым мы имеем возможность определять вероятность нахождения объекта в определенной зоне.
Пусть добавляется еще одна базовая станция, для которой нам не известны распределения вероятностей сил сигнала в зонах. Периодически мобильный узел выходит в эфир и точки доступа регистрируют силу входного сигнала. Для простоты предположим, что сила входного сигнала от новой базовой станции во время ьго по счету измерения является непрерывной случайной величиной X\ с плотностью Рг(х). Таким образом, предполагается, что мобильный объект меняет свое местоположение от измерения к измерению, причем в момент регистрации измерения местоположение объекта неизвестно.
Проводится М независимых измерений, в результате которых сила сигнала приняла значения х\,..., хм ■ Таким образом, значение
Карта слышимости устройства А590. Числа внутри ячеек метровой длины соответствуют средним потерям сигнала в децибелах
- это реализация случайной величины г = 1
В соответствии с моделью (2) предполагается, что условная плотность распределения силы сигнала в зоне ] равна
р{х\г,а^) =
<7\/27Г
ехр -
2сг2
, 3 = 1,■
(3)
где параметр и является известным и одинаковым для всех зон ] = 1,N. Параметры -условные математические ожидания значений сил сигнала в зонах ] = 1,..., N неизвестны. Требуется найти их точечные оценки.
Для удобства введем обозначения х = (ж1, • • • 5 ям') ПО (&1, • • • 5 Одг)-
Пусть дц - известные вероятности нахождения объекта в точке ] в момент регистрации г-го измерения. Вероятности могут быть рассчитаны по формулам (1) теми базовыми станциями, для которых имеются готовые карты сил сигналов, либо иным способом. По формуле полной вероятности для плотности распределения Рг(х) силы сигнала в момент измерения г справедливо следующее выражение:
N
Рг(.х) =
3=1
Будем искать оценки параметров оу с помощью метода максимального правдоподобия.
Запишем функцию правдоподобия для измерений XI, ..., хм-
М N
Ь(х; а) = П ^2р(х^; щ) ду. (4)
%=1 з=1
Строго говоря, функция (4) не является классической функцией правдоподобия, так как содержит плотности вероятностей различных распределений. Но если дц не будут зависеть от г, то есть <7^ = для всех г = 1,М, то (4) является функцией правдоподобия.
Для получения точечных оценок параметров а найдем максимум этой функции. Удобнее это сделать при помощи поиска максимума логарифмической функции правдоподобия:
М / N
1пЬ(х\а) = I ) .
1=1 \л=1
Частные производные по логарифмической функции правдоподобия имеют следующий вид:
дак
1=1
(5)
©
где
Сг/г(®) —
ехР (-(Жі2д-°"'> ) Яік Е ехР (-(%~°з)2) ду
7 = 1 4 7
(6)
м
г=1
~~М
X)
г=1
а = /(а).
(8)
Для численного решения системы уравнений (8) будем использовать метод простых итераций [1]. Опишем итерационную процедуру:
а(°) = (0,..., 0),
а<*+1> =/(аЮ), і = 1, 2,...
Последовательность вычисляется до
тех пор, пока 11 а№ — а^+1^ 11 > є, где є - заданная точность.
Если в качестве начального приближения взять а= (0,..., 0), то
м
X/ Чік%і
(1) _ і=1_ ак ~ м
(9)
Величины Сгк{о) в дискретном случае имеют смысл вероятностей нахождения объекта в точке к при условии регистрации измерения Х\.
Приравняв к нулю производные, получим систему уравнений:
X/ Яік
г=1
Обозначим
Ък (а) = --------
X) с1к(а)
1=1
При фиксированных значение гик величина 7может быть проинтерпретирована как вес г-го измерения по отношению к зоне к. Тогда
м
ак = ^2ък(а)хи к = 1,...,ЛГ. (7)
г=1
То есть математическое ожидание силы сигнала в точке к равно взвешенной сумме сил сигналов. Введем обозначения:
м
1к{а) = '^2'Угк{а)хг, к =
г=1
/(а) = (/1(а),...,/лг(а)).
В этих обозначениях система (7) запишется следующим образом:
что уже будет близко к решению уравнения (8).
Эксперименты
Для проверки адекватности предложенного подхода было проведено два вычислительных эксперимента. Первый заключался в имитации перемещения мобильного объекта по зонам помещения и программной генерации регистрируемых сил сигналов от базовых станций. Второй эксперимент проводился на реальных данных, собранных в ходе эксплуатации действующей системы локации.
Опишем имитационный эксперимент. Вначале необходимо было определить число зон N, число базовых станций И, число измерений М, выбрать значение среднеквадратического отклонения а и условные математические ожидания сил сигнала базовых станций (1 = 1,..., .О в зонах ] = 1,... Все, кроме одной, базовые станции использовались для определения локации мобильного объекта, после чего для последней базовой станции строилась карта зон видимости. В ходе эксперимента для каждого из измерений г = 1,М случайным образом выбиралась зона ] и для всех базовых станций с1 = 1,..., О программно генерировались силы сигнала хм как реализации нормально распределенных случайных величин с математическими ожиданиями ajc^ и среднеквадратическим отклонением а. На основе значений Хга от всех, кроме одной, базовых станций (1 = 1,..., -О — 1 по формуле (1) вычислялись вероятности дц нахождения мобильного объекта в зонах в момент г-го измерения. После имитации всех измерений с помощью предложенного в статье алгоритма вычислялись оценки математических ожиданий
сил сигнала для последней базовой (с? = И) станции во всех зонах ] = 1,..., N. Найденные значения сравнивались с при помощи
метрики А = Х^=1 |а? - а?т>|-
При реализации имитационного эксперимента были использованы следующие значения параметров: число зон N = 10, число базовых станций И = 5. Для каждой из пяти базовых станций (б. с.) для каждой зоны бы-
ли назначены значения математических ожиданий сил сигнала (см. табл. 1).
ния aj сравнивались с при помощи метрики А = ^2f= і \dj-ajD\- Итоги эксперимента
таковы. Значение А = 1,2, число зон с различным значением — а^о\ представлено в табл. 3.
Таблица 3. Число зон с различным значением
К — а]о\
ац^І число зон
[0-1) 56
[1-2) 28
[2-3) 9
[3-4) 4
>4 3
Вычисления проводились для различных значений М и а. В табл. 2 представлена зависимость А от М и а.
Таблица 2. Вычисленные значения Д
а 25 50 75 100 125 150
0,5 0,29 0,55 0,48 0,41 0,46 0,53
1 0,87 0,66 0,59 0,5 0,54 0,48
1,5 1,25 0,7 0,72 0,65 0,6 0,53
2 1,24 1 0,87 0,77 0,7 0,5
2,5 1,94 0,87 1,08 1,34 0,91 0,91
3 2Д 1,14 1,17 1,3 0,78 1,02
Была отмечена быстрая сходимость алгоритма, что дает основание использовать в качестве приближенного значения для формулу (9).
Теперь опишем эксперимент, проведенный на реальных данных. В помещении 1Т-парка Петрозаводского государственного университета было выделено 100 зон (ЛА = 100) и задействовано б базовых станций (И = 6). Для определения оценок математических ожиданий сил сигнала базовых станций й =
1,..., -О было проведено по 80 замеров в каждой зоне 2 = 1,..., N.
После этого в каждой зоне было проведено по 10 измерений (М = 1000). Используя данные от первых пяти базовых станций, по формуле (1) в момент каждого измерения г вычислялись вероятности дц нахождения объекта в зонах j = 1,..., N. Затем с использованием предложенного алгоритма вычислялись оценки математических ожиданий в зонах ] = 1,..., N для шестой базовой станции. Параметр а был равен трем. Найденные значе-
Заключение
В статье рассмотрена задача автоматической рекалибрации системы расчетов локации при изменении внешних условий. Предложен способ оценки неизвестных параметров распределения мощности сигнала в ходе эксплуатации системы. По мнению авторов, результаты экспериментов позволяют считать, что предложенный метод может быть использован на практике.
Исследования проводились в рамках деятельности МИП ООО «Наносети» и ЗАО «РТЛ-Сервис», а также поддерживались Петрозаводским государственным университетом (Программа стратегического развития ПетрГУ на 2012-2016 гг.), Министерством образования и науки РФ (гос. контракт 14-ВВВ.21.0162), Фондом СР МФП в НТС, Американским фондом гражданских исследований и развития (С1ЮР) и Министерством экономического развития Республики Карелия.
Литература
1. Калиткин Н. Н. Численные методы. М.: Наг ука, 1978. 512 с.
2. Мощевикин А. П., Галов А. С., Волков А. С. Локация в беспроводных сетях датчиков стандарта папоЬОС // Информационные технологии. 2011. № 8. С. 43-47.
3. Мощевикин А. П., Галов А. С., Волков А. С. Точность расчета локации в беспроводных сетях датчиков стандарта папоЬОС // Информационные технологии. 2012. № 9. С. 37-41.
0
4. Chandrasekaran G., Ergin M., Yang J., Liu S., Chen Y., Gruteser M., Martin R. Empirical Evaluation of the Limits on Localization Using Signal Strength: Beyond Cramer-Rao Bounds // Proceedings of IEEE SECON. 2009. P. 1-9.
5. Elnahrawy E., Li X., Martin R. The Limits of Localization Using Signal Strength: A Comparative Study // Proceedings of IEEE SECON. 2004. P. 406-414.
6. Galov A., Moschevikin A., Voronov R. Combination of RSS localization and ToF ranging for increasing positioning accuracy indoors // Proceedings of the 11th International Conference on ITS Telecommunications (ITST). 2011. P. 299304.
7. Ghassemzadeh S. S., Jana R., Rice C. IV., Turin W., Tarokh V. Measurement and modeling of an ultra-wide bandwidth indoor channel // IEEE Tran. Commun. 2004. Vol. 52. P. 1786-1796.
8. Molisch A.F., Balakrishnan K., Chong C. C., Emami S., Fort A., Karedal J., Kunisch J., Schantz H., Schuster U., Siwiak K. IEEE 802.15.4a ChannelModel - FinalReport, Tech. Rep., Document IEEE 802.1504-0062-02-004a, 2005.
9. Molisch A. F. Ultrawideband Propagation Channels - Theory, Measurement, and Modeling // IEEE Transactions on Vehicular Technology. 2005. Vol. 54, No. 5. P. 1528-1545.
10. Moschevikin A., Galov A., Soloviev A., Mikov A., Volkov A., Reginya S. Realtrac technology overview // EvAAL 2013, Communications in Computer and Information Science series CCIS. 2013. Vol. 386, P. 60—71.
11. Osama W. Ata, Ala’Eddin M. Shahateet, Mutaz I. Jawadeh, Adnan I. Amro. An Indoor Propagation Model Based on a Novel Multi Wall Attenuation Loss Formula at Frequencies 900 MHz and 2.4 GHz // Wireless Personal Communications. 2013. Vol. 69, Issue 1. P. 23-36.
СВЕДЕНИЯ ОБ АВТОРАХ:
Воронов Роман Владимирович
доцент, к. т. н.
Петрозаводский государственный университет, математический факультет пр. Ленина, 33, Петрозаводск,
Республика Карелия, Россия, 185910 эл. почта: [email protected] тел.: (8142) 711068
Лукашенко Олег Викторович
младший научный сотрудник, к. ф.-м. н. Институт прикладных математических исследований Карельского научного центра РАН ул. Пушкинская, 11, Петрозаводск,
Республика Карелия, Россия, 185910 эл. почта: [email protected] тел.: (8142) 763370
Мощевикин Алексей Петрович
доцент, к. ф.-м. н.
Петрозаводский государственный университет,
физико-технический факультет
ул. Университетская, 10а, Петрозаводск,
Республика Карелия, Россия, 185910
эл. почта: [email protected]
тел.: (8142) 719671
Voronov, Roman
Petrozavodsk State University, Faculty of Mathematics 33 Lenina St., 185910 Petrozavodsk, Karelia, Russia e-mail: [email protected] tel.: (8142) 711068
Lukashenko, Oleg
Institute of Applied Mathematical Research,
Karelian Research Centre, Russian Academy of Sciences 11 Pushkinskaya St., 185910 Petrozavodsk,
Karelia, Russia
e-mail: [email protected] tel.: (8142) 763370
Moschevikin, Alexey
Petrozavodsk State University,
Faculty of Physical Engineering
10A Universitetskaya St., 185910 Petrozavodsk,
Karelia, Russia
e-mail: [email protected] tel.: (8142) 713237
