发明名称 一种基于微惯性传感器的室内定位方法
摘要 本发明公开了一种基于微惯性传感器的室内定位方法,其具体步骤为:搭建包括微惯性传感器(1)、磁阻电子罗盘(2)、信息采集处理装置(3)、导航解算处理器(4)、显控装置(5)、通信总线A(6)、通信总线B(7)和通信总线C(8)的微惯性室内定位系统;微惯性室内定位系统开机后,误差补偿模块(11)补偿测量数据误差,静动态测量数据辨识模块(12)辨识载体静态与动态数据,初始对准模块(13)解算载体初始姿态角,捷联惯性导航解算模块(14)解算载体定位信息,显控装置(5)输出载体定位信息。本方法解决了卫星定位在室内存在盲区、其它定位技术信息采集工作量大、算法复杂、数据处理困难、硬件设备成本高的问题。
申请公布号 CN105136142A 申请公布日期 2015.12.09
申请号 CN201510666248.6 申请日期 2015.10.16
申请人 北京机械设备研究所 发明人 卞伟伟;葛晓飞;高文劭
分类号 G01C21/18(2006.01)I;G01C21/20(2006.01)I 主分类号 G01C21/18(2006.01)I
代理机构 中国航天科工集团公司专利中心 11024 代理人 岳洁菱;姜中英
主权项 一种基于微惯性传感器的室内定位方法,其特征在于具体步骤为:第一步 搭建微惯性室内定位系统微惯性室内定位系统,包括:微惯性传感器(1)、磁阻电子罗盘(2)、信息采集处理装置(3)、导航解算处理器(4)、显控装置(5)、通信总线A(6)、通信总线B(7)和通信总线C(8);其中微惯性传感器(1)包括:三轴微惯性加速度计(9)和三轴微惯性陀螺仪(10);导航解算处理器(4)包括:误差补偿模块(11)、静动态测量数据辨识模块(12)、初始对准模块(13)和捷联惯性导航解算模块(14);微惯性传感器(1)与磁阻电子罗盘(2)封装于一体,均通过通信总线A(6)与信息采集处理装置(3)连接;信息采集处理装置(3)通过通信总线B(7)与导航解算处理器(4)连接;导航解算处理器(4)通过通信总线C(8)与显控装置(5)连接;显控装置(5)用于人机交互;第二步 微惯性室内定位系统开机在需要定位的室内载体上安装好微惯性室内定位系统后,控制显控装置(5)给微惯性室内定位系统加电,微惯性传感器(1)与磁阻电子罗盘(2)开始采集载体数据并经数据采集处理装置(3)处理后传输给导航解算处理器(4);第三步 误差补偿模块(11)补偿测量数据误差导航解算处理器(4)接收到微惯性传感器(1)和磁阻电子罗盘(2)的测量数据后,由误差补偿模块(11)进行误差补偿;考虑三轴微惯性加速度计(9)的零偏、安装误差、随机漂移误差项,忽略二阶以上动态小量误差,建立三轴微惯性加速度计(9)的误差模型为:<img file="dest_path_image002.GIF" wi="152" he="25" />(1)式中,<img file="dest_path_image004.GIF" wi="17" he="18" />为三轴微惯性加速度计(9)输出值;<img file="dest_path_image006.GIF" wi="21" he="25" />为三轴微惯性加速度计(9)零偏;<img file="dest_path_image008.GIF" wi="25" he="24" />为三轴微惯性加速度计(9)标度系数;<img file="dest_path_image010.GIF" wi="21" he="26" /><img file="dest_path_image012.GIF" wi="88" he="24" />为三轴微惯性加速度计(9)<img file="dest_path_image014.GIF" wi="10" he="18" />轴对<img file="dest_path_image016.GIF" wi="14" he="21" />轴的正交误差系数;<img file="dest_path_image018.GIF" wi="14" he="16" />为载体运动输入加速度;<img file="dest_path_image020.GIF" wi="16" he="20" />为三轴微惯性加速度计(9)随机误差;考虑三轴微陀螺仪(10)的零偏、安装误差、正交误差和随机漂移误差,忽略二阶以上动态小量误差,建立三轴微陀螺仪(10)的误差模型为:<img file="dest_path_image022.GIF" wi="193" he="25" />(2)式中,<img file="dest_path_image024.GIF" wi="18" he="20" />为三轴微陀螺仪(10)输出值;<img file="dest_path_image026.GIF" wi="22" he="25" />为三轴微陀螺仪(10)的零偏;<img file="dest_path_image028.GIF" wi="26" he="25" />为三轴微陀螺仪(10)的标度系数;<img file="dest_path_image030.GIF" wi="17" he="16" />为载体运动输入角速度;<img file="dest_path_image032.GIF" wi="14" he="16" />为三轴微陀螺仪(10)随机误差;<img file="dest_path_image034.GIF" wi="22" he="26" /><img file="338482dest_path_image012.GIF" wi="88" he="24" />为三轴微陀螺仪(10)<img file="193306dest_path_image014.GIF" wi="10" he="18" />轴对<img file="159993dest_path_image016.GIF" wi="14" he="21" />轴的安装误差系数;<img file="dest_path_image036.GIF" wi="24" he="26" /><img file="866787dest_path_image012.GIF" wi="88" he="24" />为与加速度有关的一次项误差系数;建立磁阻电子罗盘(2)的误差模型为:<img file="dest_path_image038.GIF" wi="362" he="25" />(3)式中,<img file="dest_path_image040.GIF" wi="17" he="18" />为磁阻电子罗盘(2)的输出值;<img file="dest_path_image042.GIF" wi="22" he="25" />为预先输入的载体航向角;<img file="491279dest_path_image004.GIF" wi="17" he="18" />、<img file="dest_path_image044.GIF" wi="17" he="18" />、<img file="dest_path_image046.GIF" wi="17" he="20" />、<img file="dest_path_image048.GIF" wi="18" he="18" />、<img file="dest_path_image050.GIF" wi="17" he="18" />为磁阻电子罗盘(2)的补偿系数;误差补偿模块(11)对测量数据进行误差补偿后传输给静动态测量数据辨识模块(12);第四步 静动态测量数据辨识模块(12)辨识载体静态与动态数据静动态测量数据辨识模块(12)接收到经误差补偿的测量数据后,根据微惯性传感器(1)测量数据的变化情况辨识出载体的静态与动态,并将测量数据分解为静态数据与动态数据;其中,静态数据用于载体姿态角的初始对准,动态数据用于载体的姿态更新、速度更新与位置更新;第五步 初始对准模块(13)解算载体初始姿态角初始对准模块(13)根据静动态测量数据辨识模块(12)的静态数据进行载体姿态角的初始对准解算;由三轴微加速度计(9)输出的加速度静态信息确定载体的俯仰角<img file="dest_path_image052.GIF" wi="14" he="20" />和横滚角<img file="dest_path_image054.GIF" wi="14" he="18" />,公式为:<img file="dest_path_image056.GIF" wi="214" he="32" />(4)<img file="dest_path_image058.GIF" wi="138" he="26" />(5)式中,<img file="dest_path_image060.GIF" wi="24" he="26" />、<img file="dest_path_image062.GIF" wi="24" he="28" />、<img file="dest_path_image064.GIF" wi="24" he="25" />为加速度计的输出值,<img file="dest_path_image066.GIF" wi="16" he="18" />为重力加速度;由磁阻电子罗盘(2)确定载体的航向角<img file="954489dest_path_image040.GIF" wi="17" he="18" />:<img file="dest_path_image068.GIF" wi="50" he="25" />(6)通过式(4)、式(5)和式(6)得到载体初始静止状态下的初始姿态角;初始对准模块(13)将载体的初始姿态角信息传送给导航解算处理器(4)内的捷联惯性导航解算模块(14);第六步 捷联惯性导航解算模块(14)解算载体定位信息捷联惯性导航解算模块(14)根据静动态辨识模块(12)辨识出的动态数据与初始对准模块(13)获得的载体初始姿态角,采用捷联惯性导航算法对结果进行连续积分处理,并转换到导航坐标系中,先解算出载体的姿态角信息:俯仰角<img file="533106dest_path_image052.GIF" wi="14" he="20" />、航向角<img file="918957dest_path_image040.GIF" wi="17" he="18" />和横滚角<img file="24054dest_path_image054.GIF" wi="14" he="18" />;再解算出载体的速度信息:横向速度<img file="dest_path_image070.GIF" wi="20" he="25" />、前向速度<img file="dest_path_image072.GIF" wi="20" he="26" />和纵向速度<img file="dest_path_image074.GIF" wi="18" he="24" />;最后解算出载体的位置信息横向位移<img file="dest_path_image076.GIF" wi="14" he="16" />、前向位移<img file="dest_path_image078.GIF" wi="16" he="18" />和纵向位移<img file="dest_path_image080.GIF" wi="14" he="14" />;载体的姿态角信息解算:载体坐标系<img file="dest_path_image082.GIF" wi="14" he="20" />至导航坐标系<img file="dest_path_image084.GIF" wi="14" he="16" />的姿态矩阵<img file="dest_path_image086.GIF" wi="17" he="18" />为:<img file="dest_path_image088.GIF" wi="586" he="76" />(7)由初始对准模块(13)得到的初始姿态角通过公式(7)得到初始状态矩阵<img file="dest_path_image090.GIF" wi="20" he="25" />,为载体的姿态更新提供了初值,由初始状态矩阵求出初始四元数的值;<img file="dest_path_image092.GIF" wi="297" he="102" />(8)<img file="dest_path_image094.GIF" wi="152" he="50" />(9)用经过误差补偿的三轴微陀螺仪(10)测量的角速度<img file="dest_path_image096.GIF" wi="30" he="26" />、<img file="dest_path_image098.GIF" wi="32" he="28" />、<img file="dest_path_image100.GIF" wi="30" he="26" />,对式(8)进行四阶龙格—库塔法解算,并依据式(9)做归一化处理,实现四元数的实时更新;<img file="dest_path_image102.GIF" wi="410" he="78" />(10)经过式(8)和式(9)实时计算出四元数后,由式(10)完成姿态矩阵的更新,并根据式(7)和式(10)获得载体的姿态角信息;载体的速度信息解算:经过补偿的三轴微加速度计(9)测量的比力信息<img file="dest_path_image104.GIF" wi="22" he="26" />、<img file="dest_path_image106.GIF" wi="22" he="28" />、<img file="dest_path_image108.GIF" wi="22" he="25" />,通过姿态矩阵<img file="28587dest_path_image086.GIF" wi="17" he="18" />与室内坐标系的关系转化到室内坐标系,通过一次积分进行室内坐标系内载体的速度更新;<img file="dest_path_image110.GIF" wi="157" he="78" />(11)结合前面解算出的载体的姿态角信息,对式(8)进行二阶龙格‑库塔法解算,得到载体的速度信息;载体的位置信息解算:<img file="dest_path_image112.GIF" wi="88" he="78" />(13)对式(13)再进行一次积分解算,得到载体的位置信息;导航解算处理器(4)输出载体的姿态角信息,包括:俯仰角<img file="969868dest_path_image052.GIF" wi="14" he="20" />、航向角<img file="104218dest_path_image040.GIF" wi="17" he="18" />和横滚角<img file="814554dest_path_image054.GIF" wi="14" he="18" />,速度信息,包括:横向速度<img file="119765dest_path_image070.GIF" wi="20" he="25" />、前向速度<img file="735292dest_path_image072.GIF" wi="20" he="26" />和纵向速度<img file="948098dest_path_image074.GIF" wi="18" he="24" />,位置信息,包括:横向位移<img file="965470dest_path_image076.GIF" wi="14" he="16" />、前向位移<img file="487588dest_path_image078.GIF" wi="16" he="18" />和纵向位移<img file="91875dest_path_image080.GIF" wi="14" he="14" />;解算得到的载体定位信息由通信总线C(8)传输给显控装置(5);第七步 显控装置(5)输出载体定位信息显控装置(5)将接收到的载体的姿态角信息、速度信息、位置信息输出显示,从而完成了基于微惯性传感器的室内定位。
地址 100854 北京市海淀区北京142信箱208分箱