Избыточный измерительный модуль на микромеханических датчиках для малогабаритных систем персональной навигации

Избыточный измерительный модуль на микромеханических датчиках для малогабаритных систем персональной навигации

П.С. Маринушкин1, Т.Г. Нестеренко2

1 Сибирский федеральный университет, Красноярск 2Национальный исследовательский Томский политехнический университет

Аннотация: В статье рассматривается инерциальный измерительный блок с избыточным количеством комбинированных датчиков первичной информации, предназначенный для использования в составе малогабаритных систем персональной навигации. Делается вывод о целесообразности использования в качестве основания измерительного блока правильной четырехугольной пирамиды, на боковых гранях и основании которой расположены чувствительные элементы. В качестве чувствительных элементов в блоке используются комбинированные датчики ускорения, угловой скорости и напряженности магнитного поля, выполненные в едином корпусе. Обсуждается математическая модель и алгоритм обработки информации рассматриваемого избыточного измерительного блока. Ключевые слова: акселерометр, гироскоп, избыточность, инерциальная навигация, инерциальный измерительный модуль, комплексирование, магнитометр, микромеханический датчик, носимая электроника, персональная навигация.

Введение

Достижения последних двух десятилетий в области технологии микроэлектромеханических систем (МЭМС) сделали возможным создание малогабаритных инерциальных систем персональной (пешеходной) навигации, предназначенных как для военного, так и гражданского применения [1, 2]. Такие системы имеют ряд преимуществ перед радионавигационными системами, как то: автономность работы, высокая частота обработки информации; помехоустойчивость. В то же время они не лишены и недостатков, основным из которых является свойство накапливать с течением времени ошибки определения навигационных параметров, что определяет необходимость периодической коррекции измерений [3, 4]. Поэтому в настоящее время большой объем научных исследований направлен на создание технических и алгоритмических решений, обеспечивающих улучшение качества навигационной информации в инерциальных системах персональной навигации. В частности, одним из

разрабатываемых направлений, позволяющих снизить тенденцию к накоплению погрешностей, является комплексирование выходной информации избыточного числа неортогонально ориентированных друг относительно друга датчиков первичной информации (акселерометров, гироскопов и магнитометров).

Математическая модель избыточного измерительного модуля

Важным вопросом разработки малогабаритных систем персональной навигации, реализующих принцип функциональной избыточности, является проектирование блоков измерителей первичной информации, поскольку только определенные схемы пространственного расположения избыточных измерителей первичной информации помогают повысить точность и надежность навигационных систем [5 - 7].

Представляется, что, учитывая специфику малогабаритных систем персональной навигации, для их построения наиболее перспективными выглядят простейшие конфигурации измерителей кубической или пирамидальной формы, которые можно реализовать в миниатюрном исполнении, например на основе пластин кремний-на-изоляторе (SOI) с гибкими соединениями из полиимида [8]. В качестве возможного варианта выполнения такого измерителя рассмотрим избыточную схему (рис. 1), включающую в себя в качестве источников информации пять комбинированных датчиков, объединяющих акселерометр, гироскоп и магнитометр в едином корпусе. Таким образом, заданный набор датчиков образует три измерительных канала: кажущегося ускорения, угловой скорости и напряженности магнитного поля. Датчики расположены на боковых гранях и основании правильной четырехугольной пирамиды. Система координат блока представлена тремя взаимно перпендикулярными осями координат X, Y и Z, пересекающимися в точке O. Ось Z совпадает с осью симметрии пирамиды, а оси X и Y лежат в плоскости её основания. Ось

1

чувствительности первого датчика направлена вдоль оси 2, а оси чувствительности остальных датчиков - под углом а к плоскости основания перпендикулярно боковым граням пирамиды, и развернуты друг относительно друга на углы, кратные 20.

Рис. 1. - Конфигурация избыточного измерительного блока

Обозначив единичный вектор ориентации оси чувствительности 1-го комбинированного датчика через V, можно, в соответствии с рис. 1, записать:

¥{= 1 cosаcosP + j cosаsinP + k sinа ,

где 1, ] и и - взаимно перпендикулярные единичные вектора по осям ОХ, ОУ и О2 связанной с измерительным блоком системы координат.

Используя (1), матрицу измерений Н, преобразующую выходные сигналы из системы координат, связанной с датчиками, в систему координат измерительного блока, можно записать следующим образом:

Н

ж ж . ж

cos—cosB cos — smB sm —

2 2 2

cosacosP cosаsinP sinа

- cosаcosP cosаsinP sinа

- cosаcosP - cosаsinP sinа cosаcosP - cosаsinP sinа

Угол в в (2) для рассматриваемой схемы избыточного блока целесообразно принять равным 45°, чтобы гарантировать равномерное расположение комбинированных датчиков в пространстве и обеспечить

1К1 Инженерный вестник Дона. №3 (2016) Н| ivdon.ru/ru/magazine/arcliive/n3y2016/3709

одинаковую точность определения всех компонент каждого из измеряемых векторов вдоль разных осей. Что касается критериев выбора угла а, то величина угла а является оптимальной, если ковариационная матрица Р ошибок оценки компонент измеряемого вектора имеет наименьший след (сумму диагональных элементов) [9, 10]:

Р = и2 (Нг Н ) ,

где о - среднеквадратические отклонение ошибок измерений (датчики в каждом измерительном канале считаем равноточными).

Тогда, минимизируя след матрицы Р, получим, что оптимальная величина угла а для рассматриваемого избыточного блока составляет 24,095 градуса. Это означает, что при данной ориентации осей чувствительности датчиков дисперсия ошибок по осям ОХ, ОУ, О2 уменьшится в 2cos(а) , т.е. в д/1,667 раз (по сравнению с блоком, в котором используется минимально

необходимое число датчиков с ортогонально расположенными осями чувствительности) и соответствующая матрица измерений Н примет следующий вид:

0 0 1

а результирующая матрица измерений всех измерительных каналов блока Н2 запишется как:

Н 05x3 05x3 "

Н г = 05x3 Н 05x3

Ахз 05x3 Н

где 0иХи - нулевая матрица размера (и*и).

Далее для рассматриваемой избыточной схемы связь вектора состояния х с измерениями г можно выразить с помощью следующей модели измерений:

г = Н Е х + V ,

где г - вектор выходных сигналов комбинированных датчиков; х -компоненты входного вектора вдоль осей ОХ, ОУ и О2; V - вектор случайных помех, сопровождающих измерения.

Вектора модели (3) имеют вид: х = [[ л Л И иу иг Г,

г = [ ... г 15 Г ,

£ = [  115 Г ,

где/х,/у,/2, юх, Юу, ю2, Их, Иу, Иг, - ускорения (кажущиеся), угловые скорости и напряженности магнитного поля соответственно вдоль осей ОХ, ОУ и 02; zг -выходной сигнал 1-го измерителя ( = 1, ..., 15); - ошибка 1-го измерителя.

Алгоритм обработки информации избыточного измерительного блока

В целом, алгоритм функционирования избыточных навигационных несколько усложняется по сравнению с неизбыточными [1, 2]. В частности, возникает необходимость пересчета составляющих измеряемых векторов с осей чувствительности датчиков в систему координат измерительного блока.

Данная задача сводится к решению системы линейных алгебраических уравнений (3). Но, поскольку число уравнений в системе (3) больше, чем число неизвестных, то она является переопределенной и её точного решения не существует. Поэтому для решения данной системы целесообразно воспользоваться методом наименьших квадратов, схема которого применительно к решению системы (3) состоит из нахождения матрицы, обратной матрице измерений НЕ и умножения полученного результата на

1

вектор измерений г. Получившийся результат дает вектор-столбец значений неизвестных х:

х = Н * г

где Н* - псевдообратная матрица, вычисляемая по формуле:

Н * = (Н * Н е) Н * •

Сказанное можно обобщить в виде схемы, представленной на рис. 2.

Комбинированный датчик I

Ш/ь

Комбинированны {Л 0)2 ^

датчик 2

Комбинированный^

датчик 3

Комбинирован ный^ [_4' С°4'

датчик 4

Комбинированный^) 5'

датчик 5

Фильтрация

Л

/у

/: щ

щ

сог Нх Ну И-

Рис. 2. - Схема обработки информации избыточного измерительного блока Полученный в результате вектор х, содержащий оценки компонент векторов кажущегося ускорения, угловой скорости и напряженности магнитного поля, уже может быть обработан процессорным модулем навигационной системы, в результате чего и будут определены параметры местонахождения пользователя.

Заключение

Рассмотрены принципы построения инерциального измерительного блока с избыточным количеством комбинированных датчиков первичной информации. При условии его технической реализации в микроэлектронном исполнении избыточный измерительный блок может быть использован в качестве чувствительного элемента малогабаритных навигационных систем.

В частности, на его основе в ходе дальнейшей работы предполагается построить персональную навигационную систему, характеризующуюся повышенной точностью и надежностью и, при этом, лишенную основных недостатков радионавигационных систем (необходимость наличия дополнительной инфраструктуры, низкая помехоустойчивость).

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 16-38-50039.

Литература

1. Fischer C., Sukumar P. T., Hazas M. Tutorial: implementing a pedestrian tracker using inertial sensors // Pervasive Computing, IEEE. 2013. Vol. 12. No. 2. pp.17-27. DOI: 10.1109/MPRV.2012.16

2. J. O. Nilsson, I. Skog, P. Handel, K. Hari, "Foot-mounted INS for everybody - an open-source embedded implementation", Proceedings of the 2012 IEEE/ION Position, Location and Navigation Symposium (April 2012), pp. 140-145.

3. Ломакин М.А. Особенности построения модели погрешности МЭМС-датчиков при решении навигационной задачи // Инженерный вестник Дона, 2014, №2 URL: ivdon.ru/ru/magazine/archive/n2y2014/2398.

4. Хацько Н.Е. О возможности использования инерциальных датчиков низкого и среднего класса точности в системах автоматического управления полетом летательного аппарата // Инженерный вестник Дона, 2013, №3 URL: ivdon.ru/magazine/archive/n3y2013/1756.

5. J. Cheng, J. Dong, R. J. Landry, D. A. Chen, "Novel optimal configuration form redundant MEMS inertial sensors based on the orthogonal rotation method," Sensors, vol. 14(8), pp. 13661-13678, 2014.

6. A. Waegli, S. Guerrier, J. Skaloud, "Redundant MEMS-IMU integrated with GPS for performance assessment in sports, " Position, Location and Navigation Symposium, 2008 IEEE/ION, pp. 1260-1268, 2008.

7. Seong Yun Cho, Chan Gook Park, "A calibration technique for a redundant IMU containing low-grade inertial sensors," ETRI Journal, vol. 27, no. 4, pp. 418— 426, 2005.

8. S. A. Zotov, M. C. Rivers, A. A. Trusov, A. M. Shkel, "Folded MEMS pyramid inertial measurement unit," Sensors Journal, IEEE, vol. 11, no. 11, pp. 2780-2789, 2011.

9. Неортогональная БИНС для малых БПЛА. URL: special.habrahabr.ru/kyocera/p/114513/ (дата обращения 01.04.2016).

10. Алешкин М. В. Математические модели, методы и алгоритмы обработки избыточной информации измерительного блока / В. В. Алешкин, А. С. Матвеев, М. В. Алешкин // Интернет и инновации: сб. трудов Междунар. конф. Саратов: СГТУ, 2008. С. 377-380.

References

1. Fischer C., Sukumar P. T., Hazas M. Tutorial: implementing a pedestrian tracker using inertial sensors. Pervasive Computing, IEEE. 2013. Vol. 12. No. 2. pp.17-27. DOI: 10.1109/MPRV.2012.16

2. Nilsson J. O., Skog I., Handel P., Hari K. Foot-mounted INS for everybody - an open-source embedded implementation. Proceedings of the 2012 IEEE/ION Position, Location and Navigation Symposium (April 2012), pp. 140-145.

3. Lomakin M.A. Inzenernyj vestnik Dona (Rus), 2014, №2 URL: ivdon.ru/ru/magazine/archive/n2y2014/2398

4. Khats'ko N.E. Inzenernyj vestnik Dona (Rus), 2013, №3 URL: ivdon.ru/magazine/archive/n3y2013/1756

5. Cheng J., Dong J., Landry R. J., Chen D. A. Novel optimal configuration form redundant MEMS inertial sensors based on the orthogonal rotation method. Sensors, 2014, vol. 14(8), pp. 13661-13678.

6. Waegli A., Guerrier S., Skaloud J. Redundant MEMS-IMU integrated with GPS for performance assessment in sports. Position, Location and Navigation Symposium, 2008 IEEE/ION, 2008, pp. 1260-1268.

7. Cho S. Y., Park C. G., A calibration technique for a redundant IMU containing low-grade inertial sensors. ETRI Journal, 2005, vol. 27, no. 4, pp. 418-426.

8. Zotov S. A., Rivers M. C, Trusov A. A., Shkel A. M. Folded MEMS pyramid inertial measurement unit. Sensors Journal, IEEE, 2011, vol. 11, no. 11, pp. 2780-2789.

9. Neortogonal'naya BINS dlya malyh BPLA [Non-orthogonal Strapdown Inertial Navigation System for Small Unmanned Aerial Vehicles]. URL: special.habrahabr.ru/kyocera/p/114513/ (accessed 01.04.2016).

10. Aleshkin M. Matematicheskie modeli, metody i algoritmy obrabotki izbytochnoj informatsii izmeritel'nogo bloka [Mathematical models, methods and algorithms of processing the redundant information of the measurement unit]. Internet i innovatsii: sbornik trudov Mezhdunarodnoj konferentsii [Internet and Innovations: Proceedings of the International Conference]. Saratov, 2008. pp. 377-380.