发明名称 一种低成本AHRS辅助GPS确定整周模糊度装置及方法
摘要 本发明公开一种低成本AHRS辅助GPS确定整周模糊度装置及方法,本装置由AHRS系统和GPS构成,其中AHRS由MEMS三轴陀螺仪、MEMS三轴加速度计和MEMS三轴磁力计构成,GPS模块由两个标准GPS天线和单频GPS接收机构成,所述MEMS三轴陀螺仪、MEMS三轴加速度计、MEMS三轴磁力计和单频GPS接收机与一个数据处理与控制单元相连,所述两个标准GPS天线与单频GPS接收机相连。同时本装置采用的卡尔曼滤波方法,对数据进行融合处理,在GPS出现周跳或其它异常情况时提供AHRS的姿态角输出,保证系统有效无间断运行,不会因为GPS系统失常而影响装备的可用性。
申请公布号 CN103744101B 申请公布日期 2016.03.30
申请号 CN201410000540.X 申请日期 2014.01.02
申请人 上海大学 发明人 高同跃;费涛;葛海浪;王冬冬;史何定;沈高炜
分类号 G01S19/55(2010.01)I 主分类号 G01S19/55(2010.01)I
代理机构 上海上大专利事务所(普通合伙) 31205 代理人 何文欣
主权项 一种低成本AHRS辅助GPS确定整周模糊度的方法,采用的装置由AHRS系统和GPS模块构成,其中AHRS由MEMS三轴陀螺仪(12)MEMS三轴加速度计(11)和MEMS三轴磁力计(13)构成,GPS模块由两个标准GPS天线(14、15)和单频GPS接收机(16)构成,所述MEMS三轴陀螺仪(12)、MEMS三轴加速度计(11)、MEMS三轴磁力计(13)和单频GPS接收机(16)与一个数据处理与控制单元(17)相连,所述两个标准GPS天线(14、15)与单频GPS接收机(16)相连;该方法的操作步骤如下:(1).GPS数据采集:单频GPS接收机(16)通过天线(14、15)采集卫星的L1载波信号,信号内容包括卫星星历、载波相位,然后该接收机对(16)星历进行处理,得出接收机(16)的在WGS‑84坐标系下的坐标;(2).卡尔曼滤波:卡尔曼滤波初始的基线向量由AHRS确定的姿态角反算出来作为初始值,这样基线向量的初始值越接近真值,则可以缩短滤波时间,滤波方程如下:<img file="339927dest_path_image001.GIF" wi="324" he="199" />其中,<img file="479166dest_path_image006.GIF" wi="31" he="26" />为状态向量,<img file="dest_path_image007.GIF" wi="341" he="29" />,<img file="951735dest_path_image008.GIF" wi="41" he="26" />是<img file="65185dest_path_image006.GIF" wi="31" he="26" />预测的K+1时刻的状态向量,<img file="dest_path_image009.GIF" wi="50" he="26" />是k+1时刻的状态向量,系统噪声的方差矩阵为<img file="52733dest_path_image010.GIF" wi="22" he="25" />;<img file="dest_path_image011.GIF" wi="38" he="26" />为一步转移矩阵,<img file="277040dest_path_image012.GIF" wi="129" he="51" />,观测噪声的方差矩阵为<img file="dest_path_image013.GIF" wi="30" he="25" />,<img file="406714dest_path_image014.GIF" wi="28" he="26" />是K时刻系统的方差,<img file="dest_path_image015.GIF" wi="37" he="26" />是K时刻预测的K+1时刻的方差,<img file="691064dest_path_image016.GIF" wi="45" he="26" />是K+1时刻系统的方差,<img file="dest_path_image017.GIF" wi="14" he="18" />为单位矩阵,<img file="103591dest_path_image018.GIF" wi="31" he="25" />为双差载波相位观测量;<img file="555883dest_path_image002.GIF" wi="227" he="39" />其中,<img file="193907dest_path_image020.GIF" wi="25" he="25" />是k卫星和n卫星的载波相位的双差值;<img file="dest_path_image021.GIF" wi="265" he="102" />(7)其中,<img file="375490dest_path_image022.GIF" wi="21" he="26" />、<img file="dest_path_image023.GIF" wi="21" he="28" />、<img file="627479dest_path_image024.GIF" wi="21" he="25" />是接收机(16)到<img file="dest_path_image025.GIF" wi="14" he="20" />、<img file="dest_path_image027.GIF" wi="9" he="18" />卫星的单位向量双差值的三维分量;而<img file="527302dest_path_image028.GIF" wi="78" he="26" />,<img file="dest_path_image029.GIF" wi="18" he="26" />为接收机(16)到<img file="155730dest_path_image025.GIF" wi="14" he="20" />卫星的单位向量:<img file="457398dest_path_image030.GIF" wi="273" he="53" />(8)其中,<img file="dest_path_image031.GIF" wi="20" he="22" />、<img file="817972dest_path_image032.GIF" wi="21" he="25" />、<img file="dest_path_image033.GIF" wi="20" he="21" />为<img file="267408dest_path_image025.GIF" wi="14" he="20" />卫星在WGS‑84坐标系下的坐标;<img file="637210dest_path_image034.GIF" wi="14" he="16" />、<img file="dest_path_image035.GIF" wi="16" he="18" />、<img file="527805dest_path_image036.GIF" wi="14" he="14" />为接收机(16)在WGS‑84坐标系下的坐标;<img file="dest_path_image037.GIF" wi="16" he="20" />为<img file="324860dest_path_image038.GIF" wi="18" he="24" />载波的波长,长度为0.190293673;(3).整周模糊度去相关:采用整数高斯变换法去相关,把整周模糊度的误差协方差矩阵进行上三角变换,然后确定整周模糊度的搜索空间;(4).基线约束:通过增加约束找到正确的整周模糊度<img file="dest_path_image039.GIF" wi="30" he="20" />,由载波相位的双差观测原理得到:<img file="dest_path_image041.GIF" wi="269" he="41" />(9)式(9)中:<img file="996013dest_path_image042.GIF" wi="22" he="21" />为双差载波相位观测量,<img file="169505dest_path_image044.GIF" wi="48" he="31" />是<img file="dest_path_image045.GIF" wi="21" he="25" />中与<img file="dest_path_image047.GIF" wi="56" he="32" />对应的部分,<img file="dest_path_image049.GIF" wi="63" he="36" />为基线向量;<img file="242503dest_path_image039.GIF" wi="30" he="20" />为整周模糊度组合;<img file="210459dest_path_image050.GIF" wi="19" he="25" />为观测噪声;<img file="dest_path_image051.GIF" wi="154" he="102" />(10)<img file="368908dest_path_image022.GIF" wi="21" he="26" />、<img file="80512dest_path_image023.GIF" wi="21" he="28" />、<img file="945700dest_path_image024.GIF" wi="21" he="25" />是接收机到<img file="84557dest_path_image025.GIF" wi="14" he="20" />、<img file="667985dest_path_image027.GIF" wi="9" he="18" />卫星的单位向量双差值的三维分量;由式(9)可得基线矢量:<img file="dest_path_image053.GIF" wi="313" he="37" />(11)式(11)中:<img file="980018dest_path_image054.GIF" wi="34" he="26" />是<img file="965292dest_path_image045.GIF" wi="21" he="25" />中与<img file="dest_path_image055.GIF" wi="37" he="26" />对应的部分,<img file="275050dest_path_image056.GIF" wi="34" he="28" />是<img file="408091dest_path_image054.GIF" wi="34" he="26" />的转置矩阵,<img file="461498dest_path_image055.GIF" wi="37" he="26" />为基线向量;<img file="301278dest_path_image039.GIF" wi="30" he="20" />为整周模糊度组合,<img file="781938dest_path_image042.GIF" wi="22" he="25" />是载波相位的双差观测量,<img file="dest_path_image057.GIF" wi="148" he="26" />;可求出基线的长度<img file="136696dest_path_image058.GIF" wi="33" he="26" />:<img file="dest_path_image059.GIF" wi="128" he="33" />(12)式(12)中:<img file="993793dest_path_image055.GIF" wi="37" he="26" />为基线向量,<img file="688080dest_path_image060.GIF" wi="37" he="28" />是<img file="339641dest_path_image055.GIF" wi="37" he="26" />的转置矩阵;求出这个基线长度以后,可进行基线约束检验:<img file="dest_path_image061.GIF" wi="89" he="30" />(13)式(13)中:<img file="447275dest_path_image058.GIF" wi="33" he="26" />是解算出的基线长度,<img file="842484dest_path_image062.GIF" wi="16" he="18" />为事先测好的基线长度,<img file="dest_path_image063.GIF" wi="17" he="18" />为基线约束的门限值,可根据具体要求选取;在基线约束的过程中,剔除不满足基线约束的整周模糊度组合;(5).姿态角约束: AHRS系统在GPS 对姿态角没有精确解出之前已经给出了载体的航向角<img file="391277dest_path_image064.GIF" wi="41" he="25" />及俯仰角<img file="dest_path_image065.GIF" wi="41" he="25" />,因此可继续对整周模糊度组合进行校验,从满足基线约束的多组备选值中,逐个求解载体航向及俯仰角,把不符合姿态角约束的剔除,假设基线的安装月载体的纵轴重合,那么可求出航向角<img file="541635dest_path_image066.GIF" wi="17" he="16" />及俯仰角<img file="dest_path_image067.GIF" wi="17" he="22" />:<img file="808669dest_path_image068.GIF" wi="94" he="45" />(14)<img file="dest_path_image069.GIF" wi="157" he="30" />(15)式(14、15)中:<img file="7569dest_path_image070.GIF" wi="14" he="15" />、<img file="dest_path_image071.GIF" wi="15" he="18" />、<img file="473185dest_path_image036.GIF" wi="14" he="14" />为基线在地理坐标系下的东北天三个方向的坐标;<img file="732128dest_path_image066.GIF" wi="17" he="16" />为航向角;<img file="486458dest_path_image067.GIF" wi="17" he="22" />为俯仰角;姿态约束检验:<img file="223469dest_path_image072.GIF" wi="97" he="28" />(16)<img file="dest_path_image073.GIF" wi="97" he="28" />(17)式(16、17)中:<img file="809172dest_path_image064.GIF" wi="41" he="25" />为AHRS解算出的载体航向角,<img file="239016dest_path_image065.GIF" wi="41" he="25" />为AHRS解算出的载体俯仰角,<img file="215062dest_path_image074.GIF" wi="16" he="20" />为姿态角约束的门限值,可根据具体要求选取;(6).基线矢量结算,根据约束条件筛选出的整周模糊度求解基线向量,进而求出精确的载体姿态;    在求解出整周模糊度以后就可带入公式(11)得到载体在WGS‑84坐标系下的坐标,经过坐标转换矩阵,假设本系统基线的安装与载体的纵轴重合,则载体坐标系和地理坐标系重合:<img file="dest_path_image075.GIF" wi="317" he="76" />(19)得出基线在载体坐标系下的坐标;(7).姿态解算,由上面根据式(19)解算出基线矢量,然后根据式(14)、(15)可求出基线在载体坐标系下的航向角和俯仰角。
地址 200444 上海市宝山区上大路99号
您可能感兴趣的专利