CC BY
161
DOI: 10.15593/RZhBiomeh/2015.1.02 УДК 531.391:612.76
Российский
¡Курная
Биомеханики
www.biomech.ru
НАВИГАЦИЯ ПЕШЕХОДА С ИСПОЛЬЗОВАНИЕМ БЕСПЛАТФОРМЕННОЙ ИНЕРЦИАЛЬНОЙ НАВИГАЦИОННОЙ СИСТЕМЫ (БИНС), УСТАНОВЛЕННОЙ НА СТОПЕ
Механико-математический факультет Московского государственного университета им. М.В. Ломоносова, Россия, 119991, Москва, Ленинские горы, 1, е-mail: [email protected], [email protected]
Аннотация. Целью работы является исследование метода определения местоположения пешехода в пространстве с использованием бесплатформенной инерциальной навигационной системы (БИНС) на микроэлектронных механических сенсорах. Данный метод эффективен для применения в помещениях и в других местах, где отсутствует сигнал GPS. Координаты вычисляются интегрированием модельных уравнений БИНС. Основная проблема метода - «уход» координат из-за низкой точности используемых в БИНС микромеханических акселерометров и датчиков угловой скорости. Чтобы компенсировать накопление ошибок, метод интегрирования дополняется технологией коррекции по нулевой скорости (zero velocity update). Для этого БИНС крепится на стопе. Движение стопы можно разделить на две фазы: фаза опоры и фаза переноса. На фазе опоры, где БИНС считается неподвижной, можно уточнить координаты и ориентацию, используя информацию о равенстве нулю линейной скорости чувствительной массы БИНС. Особенностью данной работы является применение при коррекции «сигма-точечного» («unscented») варианта фильтра Калмана, что способствует повышению точности при выходе ошибок навигации из линейной зоны. Приведены результаты моделирования и обработки экспериментов. Траектория стопы в имитационной модели и в экспериментах представляет собой квадрат. В экспериментах используется блок инерциальных датчиков ADIS16405 производства компании Analog Devices. Результаты можно оценить как обнадеживающие. Возможно, использование описанной технологии позволит строить кинограммы движения стопы человека без использования системы видеоанализа.
Ключевые слова: навигация пешехода, бесплатформенные инерциальные навигационные системы, коррекция по нулевой скорости, Unscented Kalman Filter.
Задача навигации пешехода с помощью бесплатформенных инерциальных навигационных систем (БИНС) в настоящее время активно исследуется [5, 6]. Данный метод эффективен для применения в помещениях и в местах, где отсутствует сигнал GPS. Основная проблема метода - «уход» координат из-за низкой точности используемых микроэлектронных механических сенсоров. Поэтому при навигации необходима дополнительная коррекция, например, с применением магнитометров, сигналов Wi-Fi, и т.п. Целью данной работы является исследование метода определения местоположения пешехода в пространстве с использованием БИНС, установленной на
© Болотин Ю.В., Фатехрад М., 2015
Болотин Юрий Владимирович, д.ф.-м.н., профессор кафедры прикладной механики и управления, Москва
Ю.В. Болотин, М. Фатехрад
Введение
стопе. Координаты вычисляются интегрированием модельных уравнений. Чтобы компенсировать накопление ошибок, в опоре проводится коррекция с использованием информации о нулевой скорости стопы в фазе опоры. Определение фазы опоры -самостоятельная задача, требующая внимания. Особенностью данной работы в сравнении с [5, 6] является применение коррекции для коррекции сигма-точечного алгоритма фильтра Калмана - в английском варианте Unscented Kalman Filter (UKF) [7], что способствует повышению точности при выходе ошибок навигации из линейной зоны [2].
Постановка задачи
Пусть БИНС крепится на теле человека и с ее помощью требуется определить его координаты. Метод инерциальной навигации [2] используется там, где GPS не работает, например, в помещении. Для низкоточной БИНС свойственны большие ошибки в интегрировании модельных уравнений [2]; чтобы уменьшить ошибки, предлагается использовать специальный алгоритм коррекции. Для этого БИНС крепится на стопе. Движение стопы можно разделить на две фазы: фаза опоры и фаза переноса. На фазе опоры, где БИНС считается неподвижной, можно уточнить координаты и ориентацию, используя технологию коррекции по нулевой скорости (zero velocity update) [6].
Введем две системы координат: Cx - неподвижная относительно Земли и Mz -жестко связанная с чувствительным элементом БИНС, установленной на стопе. Ось x3
вертикальна; оси xl, x2 привязаны к ориентирам на местности (рис. 1).
Относительную ориентацию Cx, Mz определим углами Эйлера-Крылова:
Cx У > 9 > Y > Mz. (1)
3 2 1 w
Пусть ю - угловая скорость БИНС; f - кажущееся ускорение чувствительной
массы M БИНС. Под кажущимся ускорением здесь и ниже понимается удельная сила, создаваемая акселерометром и действующая со стороны корпуса БИНС на чувствительную массу. Обозначим Ш2, f проекции Ш, f на оси системы Mz; Шх,
f - проекции Ш, f на оси системы Cx .
Рис. 1. Примерное положение систем координат и блока БИНС. Чувствительная масса БИНС смещена относительно каблука ботинка
Введем модель погрешности датчиков угловой скорости и акселерометров в
виде
ш; =ш2 +у г +6ш2, /' = / + 8/. (2)
Здесь ш' - измеренная угловая скорость; - измеренное кажущееся
ускорение; - дрейф датчиков угловой скорости; 5шг - шум датчиков угловой скорости; 5/; - шум акселерометров.
Требуется определить координаты х чувствительной массы навигационной системы в системе координат Сх на некотором отрезке ходьбы [?0, ^ ], используя
измерения ш'(/), /' (?) и начальные условия. Длина отрезка - от долей минут до нескольких минут.
Алгоритм навигации
При решении поставленной задачи примем следующие упрощающие предположения, допустимые при использовании низкоточной инерциальной навигационной системы: угловой скоростью вращения Земли можно пренебречь; вектор силы тяжести постоянен в пространстве.
Алгоритм состоит из следующих блоков: решение модельных уравнений, выставка, определение фазы опоры, коррекция.
Модельное уравнение. Модельное уравнение описывает движение модельной точки и модельную ориентацию приборного трехгранника. Запишем модельные уравнения в виде [2, 3]:
1
д = -д х = ух, *х=В'агГ,+8х, В1х=В1(д). (3)
Здесь е К4 - кватернион поворота от Сх к модельному трехграннику Мг'; ух - модельная линейная скорость чувствительной массы относительно Земли; В\х -
матрица преобразования от системы Сх к Мг'; О обозначает умножение кватернионов [4];
Г (П „
Q' =
Vй, у
В дискретном времени tm=m А/, т = 0,1, 2, • • • уравнение Пуассона записывается в виде д^т+1) = Ад(1т) О д(?т), где введены обозначения
Ад = [Ap0,Ap]T , Ар о = cos
л к ■
Ар = т—г sin
к = ш 'А.
С использованием известных измерений /.'(?); [?0 < ? < ^] и известных
начальных условий х'(?0), V'(?0), 9(0), у'(?0), ) можно решить модельное уравнение и получить модельные координаты, скорости и углы х' (?), V' (?), 6'(?), у'(?),
относительно системы координат Сх как функции времени.
Алгоритм выставки. Рассмотрим задачу определения начальных условий х' (?0), V' (?0), 9(0), у'(?0), ) для модельных уравнений, называемую задачей
начальной выставки. Во время начальной выставки инерциальная система предполагается неподвижной. Без ограничения общности предполагаем, что во время выставки истинные значения координат х(/0) и скоростей v(t0) равны нулю, и, так как начальное положение системы считается известным, можно утверждать, что х '(г0) = 0, V '(¿0) = 0. Поскольку для грубых датчиков угловой скорости, не чувствующих
угловую скорость Земли, угол курса у не определяется, считаем у' (г0) = 0. Ищем начальные значения углов крена и тангажа 0' (г0 ), у' (г0 ), используя соотношения
f =" g = J x о x
Г 0 ^ 0
V g0 J
BZx (0, ву') =
cos в' 0 - sin в'
sin в 'sin у' cos у' cos в 'sin у' sin в 'cos у' - sin у' cos в 'cos у'
, Л = Bl ■ f,
где - нормальное значение ускорения силы тяжести [2]. Отсюда
Л =-£о ■0 Ю, Л2 = £о ■ с°8 0 Ю ■у Ю.
Разрешая последние уравнения относительно 0(0), у'(О, получим
в '(t0) = arcsin
fz
• л
g<
у '(t0) = arcsin
0 J
f'
J z-
V g0cos в '(t0)
(4)
0) J
Нахождение фазы опоры. Предположим, что определен момент гк1 начала (к -1 )-й опорной фазы и найдены модельные переменные х (к-1), V (к-1), 0 (к-1), У ХЧ-1), У ХЧ-1) в этот момент. Рассмотрим алгоритм поиска следующего момента начала опоры ^. Решая модельное уравнение на отрезке [гк_г < г < и ], где и взято заведомо больше, чем ^, найдем матрицу В^ (г) = В'ат (г) преобразования от трехгранника Мг' к Сх и проекции кажущегося ускорения на оси Сх:
ЛХ (г) = в; (г) ■ л: (г), [^ < г < г*].
Введем некоторые пороговые значения юшах и длшах. Определим фазу опоры как интервал [гк, г, на котором выполнены неравенства
(t)
<шп
fx, '(t) - gx
<Afm
(5)
Алгоритм коррекции в форме иКГ. Сигма-точечная форма фильтра Калмана -«Umcented Kalman FiUer»> (UKF) - является распространением линейной фильтрации Калмана на случай, когда уравнения динамики и наблюдения системы нелинейны и не могут быть адекватно линеаризованы. Этот алгоритм использует подход «сигма-точечного преобразования» - «unscented traщform»> (иТ), который подразумевает выбор некоего минимального набора точек (сигма-точек) для задания математического ожидания и ковариации искомого вектора с наложенными на него случайными помехами. По сигма-точкам строятся нелинейные функции прогноза, которые затем используются для вычисления матрицы ковариации. Основная идея иТ в предположении, что легче аппроксимировать функцию распределения вероятности случайной величины с ненулевым математическим ожиданием, чем аппроксимировать произвольную нелинейную функцию преобразования. Данная методика подразумевает выбор специальным методом набора точек (сигма-точек) в пространстве параметров,
который с достаточной точностью описывает статистические характеристики (математическое ожидание, дисперсию и отчасти моменты более высоких порядков) искомого случайного вектора [7].
Перейдем к описанию алгоритма Unscented Kalman Filter. Рассматривается следующая нелинейная система:
Ук+1 = f(Ук) + , % = h( Ук) + 8zk.
Здесь f - закон движения; ук - вектор состояния; Wk - шум в динамике; Zk - измерение; h - функция измерений; 8zk - шум в измерениях; к - номер фазы опоры.
Пусть у - оценка у; Р - ковариационная матрица ошибок оценки; % - сигма-точки вокруг y',L- размерность вектора состояния у; X - масштабный коэффициент; W - вес; ук - априорная оценка состояния; К - коэффициент обратной связи; ук -апостериорная оценка состояния. Инициализация фильтра имеет вид
Уа=Е[уа], Р0= Е[(у0 -у0)О0-%)т].
Сигма-точки в некоторый момент к — 1 определены формулами
ЗС*-1 = [&-1, Я-1 + y/(L + X)Pk_! , Я-1 - J(L + X)Pk_J, к е {1,..со}. Априорные оценки в следующий момент к определены формулами прогноза
2 L 2 L
Цк-1 =f(Xk-1), , Рк
1=0 1=0
Здесь веса задаются формулами
wW = W(c) = + (1 _ a2) w(C) = w(m) =-1-,
0 X + L ° X + L V 1 1 2(X + L)
-3 7
где a ~ 10 , X = a L _ L - настроечные параметры [7].
Априорные оценки измерения и ковариации их ошибок определены формулами
1=0
{Z^-z-J, РУкЧ •
1=0 1=0
Апостериорные оценки в момент к вычисляются по формулам коррекции:
К = Py^At Ук = yk+K{zk-zk) , Рк = Рк -КР2к2кКт .
Применительно к рассматриваемой задаче навигации f - результат интегрирования уравнений (3) на интервале переноса \t(k),или опоры \tk,t(k)];
размерность вектора состояния L = 9; y = [х, v, 0, у, т ; h = [х3, vl3 v2, v3]т, что соответствует предположению о движении по ровной местности. В качестве моментов коррекции рассматриваются как моменты t начала опорной фазы, так и моменты t(1) конца опорной фазы.
Имитационное моделирование
Целью имитационного моделирования является проверка потенциальной точности метода. Для изучения особенностей показаний БИНС, укрепленной на стопе в процессе ходьбы, необходимо сымитировать движения стопы. Здесь строится весьма упрощенная имитационная модель, не претендующая на биомеханическую достоверность. Рассмотрим один шаг пешехода, предполагая все шаги правой и левой ноги одинаковыми и рассматривая для определенности движение правой стопы. Шаг состоит из четырех участков:
Г0 - состояние статической опоры. В этом положении скорость и угловая скорость стопы равны нулю.
Г - процесс отрыва стопы: пятка отрывается от поверхности земли, а носок продолжает касаться земной поверхности.
Г - процесс переноса: стопа полностью оторвана от земли.
Г - процесс прислоения стопы: пятка касается земной поверхности, и носок постепенно опускается.
Обозначим Т0, Т, Т, Т длительности участков Г, Г, Г, Г . Моменты
начала соответствующих участков обозначим 1к, ^, , 43) (рис. 2). В данной модели предполагается, что, когда правая нога находится в состоянии Г, левая нога находится в состоянии Г. Таким образом, длительности Г и Г равны. Отсюда следует, что, когда левая нога находится в состоянии Г, правая нога находится в состоянии Г и Т и Т равны. Описанные основные участки движения стопы циклично повторяются при ходьбе.
В имитационной модели задаются значения х, V, 9, у, у в точках опоры ^ и точках максимальной высоты подъема стопы (£2) + /) / 2, данные точки соединяются сплайном третьего порядка, с помощью которого вычисляются значения координат, углов, скорости и ускорения в промежуточных точках (рис. 3-5). Примерное положение блока БИНС на стопе показано на рис. 1, где О - пятка. Координаты чувствительного элемента М БИНС в системе координат Ог предполагаются постоянными.
Пусть готг - радиус-вектор ОМ в системе отсчета Ог; х, vx, ах - координаты, абсолютная скорость и абсолютное ускорение точки М относительно системы координат Сх; гсох - радиус-вектор СО в системе отсчета Сх; Ог - система координат жестко связана со стопой.
Рис. 2. Участки и длительность участков одного шага
1 1,5 Время, с
Рис. 3. Графики координат и скоростей точки О при имитации:
■ х^, ¿/х^ / , - - - - Хт, — ¿/х^ / dt
Рис. 4. Графики зависимости угла и угловой скорости тангажа от времени:
-- 0: —
1 1,5 Время, с
Рис. 5. Графики зависимости углов крена и тангажа (пунктир) и соответствующих угловых скоростей (сплошная линия) от времени при имитации:
-- у; — с/у/сЛ; —¡Зху/Ш
1 м
2 з
Введем матрицы элементарных поворотов (1)
"1 0 0 " cos 9 0 - sin 9 cos у sin у 0
В(у) = 0 cos у sin у ) N 0 1 0 , ад= - sin у cos у 0
0 - sin у cos у sin 9 0 cos 9 0 0 1
Тогда имеют место формулы
~0~ V
B„ = В(у)Вф)В(у), <D= = 5(y)5(0) v(0 0 + B( y)9(0 1 + y(0 0
1 0 0
Координаты, скорость и ускорение точки М в системе отсчета Сх определены формулами
* = К Гот2 + ГсоХ , V* = К! К X Гош) + Кш, = К К X Юг Х,Гош, + Юг Х Гош.г) + ^ •
Ускорение точки М в проекции на оси Ог определено формулами а 2 = Вхах; 8 г = ^тх8х • Кажущееся ускорение /2 от удельной силы, действующей на чувствительную массу БИНС, определяется формулой / = аг - 8 .
При вычислении показаний датчиков (2) 5/, 5шг моделируются как белые шумы со средними квадратами (СКО) о[5/ ], о[5юг ], а - как случайные величины, постоянные в течение шага, со средними квадратами о[5уг ].
Модельная траектория стопы представляет собой квадрат. Приняты следующие значения параметров:
М = 0,01 с, о[5шг] = 0,05, о[5^] = 0,02, о[5/] = 0,1^ .
с с с
В табл. 1 приведена зависимость среднеквадратического отклонения ошибок навигации от длины шага. При ходьбе короткими шагами погрешность навигации уменьшается, что может объясняться большим числом точек коррекции. В табл. 2 показана зависимость ошибок навигации от параметров погрешностей акселерометров (м/с2). При увеличении погрешностей датчиков, как и следовало ожидать, точность уменьшается. На рис. 7 представлены графики истинной траектории и ее оценки.
Таблица 1
Зависимость СКО погрешностей определения координат и углов тангажа, крена и курса стопы от длины шага
Параметр Длина шага
0,5 м 1 м
а[ 5х] 0,064 0,071
а[ 59] 0,047 0,051
а[ 5у] 0,048 0,052
а[ 5у] 0,011 0,012
Таблица 2
Зависимость СКО погрешностей определения координат и углов тангажа, крена и курса стопы от величины СКО шума акселерометров о[5/ ]
Параметр Величина СКО шума акселерометров
а[5/ ] = 0,1 м/с2 а[5/ ] = 0,2 м/с2 а[5/ ] = 1,0 м/с2
а[ 5х] 0,071 0,075 0,240
а[ 59] 0,051 0,053 0,056
а[ 5у] 0,052 0,054 0,061
а[ 5у] 0,012 0,013 0,031
Навигация при реальных данных
В этом эксперименте используется блок инерциальных датчиков с микроэлектронными механическими сенсорами ADIS16405 производства компании Analog Devices. Прибор включает три магнетометра, три акселерометра и три гироскопа. Частота съема переменная, порядка 270 Гц. Данные записываются в компьютер через USB-кабель. Прибор был закреплен на деревянной дощечке, примотанной скотчем к стопе испытуемого.
В эксперименте испытуемый прошел по квадрату размером 4,7*4,7 м и вернулся в исходную точку. Угол курса, координаты и скорости в начальный момент t0 примем, как указывалось выше, равными нулю, а углы крена и тангажа определяются начальной выставкой.
*Ю = 0, V(io) = 0, е '(О = 0,10, у '(t0) = 0,13, V'(to) = 0.
Для того чтобы найти фазы опоры [tk, t(k)], примем следующие пороговые значения параметров на интервале [tk, t„ ]:
wmax < 0,039рад, fmsx < 0,15 м •
Рис. 7. График траектории движения стопы при имитации: истинная; — модельная
6 8 10 12 14 Время, с
Рис. 8. Зависимость продольной координаты стопы от времени: — истинная; — модельная; —скорректированная
Рис. 9. Зависимость угла тангажа стопы
от времени: — истинная; — модельная; —-скорректированная
с
с
Используем корректирующие измерения в начале и в конце фазы опоры:
x3(tk) = 0, vi(tk) = 0, v2(tk) = 0, = 0,
Хз«?) = 0, v1(tk1)) = 0, V2(tf) = 0, V3(tk1}) = 0.
Для навигации необходима предварительная калибровка БИНС, т.е. определение масштабов, перекосов и смещений показаний датчиков. С учетом ошибок масштабов и перекосов уравнения измерений имеют вид
/[ = К/'+ b,, ш' = Кш'+ b .
J z f j / ' 2 ш ш
Определение матриц Кf, Kffl можно проводить один раз за время жизни БИНС. Для определения К^, Kffl применен метод [1] с использованием грубого поворотного стенда (калибровка проведена А.В. Козловым). После введения масштабных матриц К, К уравнения измерений можно переписать в виде
/' = /"+ bf, ш2 =ш"+ b®,
где /", ш'' - масштабированные показания датчиков; b^ , Ьж - постоянные векторы.
Калибровка смещений нулей ввиду их нестабильности необходима перед каждым экспериментом. Калибровка смещений нулей b^ акселерометров проведена
следующим образом: блок инерциальных датчиков разворачивается в разные положения и фиксируется в течение некоторого времени At; (несколько секунд). Далее
вычисляется < /'' >. - средние показания акселерометров /'' в течение времени At .
Для определения b решается нелинейная задача минимизации:
N 9 II ||2
J(bf ) = BII4 " < /" > +bf ]2 ^ min .
После минимизации значение J уменьшилось в тринадцать раз, т.е. смещения уточнены по крайней мере в три раза. Значения дрейфов датчиков угловой скорости
вычисляются как среднее < ю" > перед началом движения.
На рис. 10 приведена траектория чувствительной массы на плоскости. Уход координат за 30 с составил порядка метра. На рис. 11 показана зависимость углов от времени. Углы крена и тангажа имеют практически нулевой дрейф. Дрейф угла курса значительно меньше, чем дрейф соответствующего датчика угловой скорости.
Следует отметить, что динамический диапазон используемых датчиков угловой скорости был 6 рад/с, в то время как угловая скорость по тангажу достигала 10 рад/с. Чтобы исправить ситуацию, показания датчиков угловой скорости на участках насыщения экстраполировались посредством полинома пятого порядка, со вторым порядком гладкости в точках склейки. Это позволило значительно улучшить точность навигации, хотя неучтенные угловые скорости, без сомнения, остались. Проблема насыщения, в совокупности с низкой точностью датчиков, привела к ошибке навигации порядка 1 м. Отметим при этом, что без коррекции модельные уравнения «разваливались» в течение двух шагов. Ожидается, что переход к новым датчикам позволит качественно улучшить ситуацию.
г=1
i f
.............
i i
О
1
хи м
Рис. 10. Траектория стопы на плоскости в эксперименте
0 5 10 15 20 25 30 40
Рис. 11. Зависимость углов курса, крена и тангажа стопы от времени в эксперименте:
Выводы
Разработана методика навигации пешехода с использованием установленной на стопе бесплатформенной инерциальной навигационной системы. Для коррекции по нулевой скорости применен алгоритм Unscented Kalman Filter. Методика проверена на модельных данных, проведены эксперименты с БИНС «ADIS16405». Алгоритмы дали ошибку порядка 1 м при движении в течение 30 с. Результаты показали, что угловая скорость стопы по тангажу даже при медленной ходьбе достигает 8 рад/с, что приводило к насыщению показаний датчиков угловой скорости. Поэтому, чтобы добиться хороших результатов, следует использовать датчики угловой скорости с динамическим диапазоном не менее 10 рад/с. По-видимому, для повышения точности навигации в течение длительного времени требуется дополнительная коррекция с использованием, например, информации о прямолинейности движения на каких-то участках.
Заметим в заключение, что достоинством инерциальных датчиков в сравнении с системой видеоанализа является высокая точность при измерениях в течение коротких промежутков времени, однако при изменениях в течение длительного времени наблюдаются значительные «уходы» координат. Возможно, использование описанной технологии позволит строить кинограммы движения стопы и других звеньев человека без использования системы видеоанализа.
Благодарности
Работа выполнена при поддержке РФФИ (грант 13-01-00604А). При экспериментах использовалось оборудование системы панорамной виртуальной реальности лаборатории математического обеспечения имитационных динамических систем (МОИДС).
Список литературы
1. Вавилова Н.Б., Парусников Н.А., Сазонов И.Ю. Калибровка бескарданных навигационных систем при помощи грубых одностепенных стендов // Современные проблемы математики и механики. -2009. - Т. 1. - С. 212-223.
2. Голован А.А., Парусников Н.А. Математические основы навигационных систем. Ч. I. - М.: Изд-во МГУ им. М.В. Ломоносова, 2010. - 132 с.
3. Голован А.А., Парусников Н.А. Математические основы навигационных систем. Ч. II. - М.: Изд-во МГУ им. М.В. Ломоносова, 2011. - 170 с.
4. Ишлинский А.Ю. Ориентация, гироскопы и инерциальная навигация. - М.: Наука, 1976. - 672 с.
5. Bebek O., Suster M., Rajgopal S., Fu M. Personal navigation via shoe mounted inertial measurement units // IEEE/RAS-EMBS International Conference on Intelligent, Robots and Systems. - Taipei, 2010. - P. 10521058.
6. Feliz R., Zalama E. Pedestrian tracking using inertial sensors // Journal of Physical Agents. - 2009. - Vol. 3, № 1. - P. 35-42.
7. Wan E., Merwe R. The unscented kalman filter for nonlinear estimation // IEEE/AS-SPCC Adaptive Systems for Signal Processing, Communications, and Control Symposium. - Alta, 2000. - P. 153-158.
PEDESTRIAN NAVIGATION WITH A STRAPDOWN INERTIAL NAVIGATION SYSTEM (SINS) MOUNTED ON THE FOOT
Y.V. Bolotin, M. Fatehrad (Moscow, Russia)
The paper investigates an approach to pedestrian positioning with a strapdown inertial navigation system (SINS) based on microelectronic mechanical sensors. This approach is particularly urgent in the indoor applications where the GPS signal is unavailable. The position is calculated by integrating the navigation equations of SINS. The main problem with the approach is coordinates error drift caused by low accuracy of micromechanical accelerometers and gyros used in SINS. To compensate the cumulative error, integrating the navigation equations is complemented with the zero velocity update technology. To implement this technology, SINS is mounted on the foot of the pedestrian. The foot motion can be divided into two phases, i.e. the support phase and the transport phase. In the support phase where SINS is assumed to be immobile, the coordinates and orientation of SINS can be corrected using the information on its zero velocity. A novel technique for the discussed problem tried in this paper is application of the unscented Kalman filter. This filter provides better results than the extended Kalman filter when the system errors are far into the zone of nonlinearity. Results of both simulation and experiments are presented and discussed. The pedestrian trajectory in both cases is a rectangle. For the experiment, the inertial measurement unit ADIS16405 by the Analog Devices Company is used. The results can be estimated as encouraging. A possibility of applying this method to restore the sequence of frames for biomechanical analysis without a videoanalysis system can be seen.
Key words: pedestrian navigation, SINS, zero velocity update, Unscented Kalman Filter.
Получено 17 октября 2014
