发明名称 一种基于陀螺测量信息约束的GNSS载波相位模糊度求解方法
摘要 一种基于陀螺测量信息约束的GNSS载波相位模糊度求解方法,其步骤为:(1)以当前卫星接收机观测到卫星的几何精度因子最小为优化目标,得到四颗可见卫星作为主卫星;(2)根据多接收机天线对四颗主卫星的载波相位观测值,计算初始的双差载波相位模糊度搜索范围;(3)根据陀螺测量信息,以动态过程载体方位角变化值为约束条件,压缩模糊度搜索范围;(4)基于压缩后的模糊度搜索范围,求解得到正确的模糊度。本发明具有原理简单、解算速度快、精度高、稳定性好等优点。
申请公布号 CN103454664B 申请公布日期 2016.02.24
申请号 CN201310364318.3 申请日期 2013.08.20
申请人 中国人民解放军国防科学技术大学 发明人 李涛;何晓峰;祝建成;唐康华;潘献飞;胡小平;罗兵
分类号 G01S19/49(2010.01)I;G01S19/04(2010.01)I 主分类号 G01S19/49(2010.01)I
代理机构 湖南兆弘专利事务所 43008 代理人 周长清
主权项 一种基于陀螺测量信息约束的GNSS载波相位模糊度求解方法,其特征在于,步骤为:(1)以当前卫星接收机观测到卫星的几何精度因子最小为优化目标,得到四颗可见卫星作为主卫星;(2)根据多接收机天线对四颗主卫星的载波相位观测值,计算初始的双差载波相位模糊度搜索范围;所述步骤(2)的具体流程为:(2.1)得到双差载波相位方程;在同一时刻,同时对接收机和卫星进行差分,得到双差载波相位方程:<img file="FDA0000791824500000011.GIF" wi="1381" he="126" />其中,<img file="FDA0000791824500000012.GIF" wi="117" he="71" />为分别对接收机1、2和卫星i、j进行双差得到的载波相位观测,<img file="FDA0000791824500000013.GIF" wi="110" he="79" />为分别对接收机1、2和卫星i、j进行双差得到的双差整周模糊度,i和j为卫星编号,λ为载波波长,<img file="FDA0000791824500000014.GIF" wi="44" he="55" />为基线矢量,ε为观测噪声;(2.2)、计算初始的双差载波相位模糊度搜索范围;对于两个接收机即单基线的情况,假设两个天线同时观测到n+1颗卫星,观测一个历元可以组成n个双差载波相位方程,其中包含3个基线矢量参数,n个双差整周模糊度参数,组成如下矩阵形式的误差方程:ΔV=ΔA·ΔX‑ΔL               (5) 其中,残差向量<img file="FDA0000791824500000015.GIF" wi="798" he="95" />设计矩阵<img file="FDA0000791824500000016.GIF" wi="1061" he="543" />其中<img file="FDA0000791824500000017.GIF" wi="822" he="79" />代表历元t星站间单位矢量之差的各分量;ΔX=[a<sub>x</sub> a<sub>y</sub> a<sub>z</sub> N<sup>1</sup> … N<sup>n</sup>]<sup>T</sup>包括基线向量[a<sub>x</sub> a<sub>y</sub> a<sub>z</sub>]和n个双差整周模糊度待估参数[N<sup>1</sup> … N<sup>n</sup>],当载体在静态情况时,待估参数为n+3维;载波相位观测向量<img file="FDA0000791824500000018.GIF" wi="1698" he="101" />为t时刻第j 个双差载波相位观测值;组成法方程为:(ΔA<sup>T</sup>ΔPΔA)ΔX‑ΔA<sup>T</sup>ΔPΔL=0            (6) 其中,<img file="FDA0000791824500000021.GIF" wi="677" he="303" /><img file="FDA0000791824500000022.GIF" wi="502" he="255" />解法方程,得到未知参数矢量:ΔX=(ΔA<sup>T</sup>ΔPΔA)<sup>‑1</sup>ΔA<sup>T</sup>ΔPΔL       (7) 得到的双差整周模糊度解称为浮点解;对于模糊度矢量<img file="FDA0000791824500000023.GIF" wi="78" he="77" />的N个元素分别进行历元t检验:<img file="FDA0000791824500000024.GIF" wi="1436" he="151" />以便将实数解<img file="FDA0000791824500000025.GIF" wi="87" he="70" />附近的满足上述条件的整数<img file="FDA0000791824500000026.GIF" wi="100" he="71" />都挑选出来;式中1‑α为置信水平;<img file="FDA0000791824500000027.GIF" wi="426" he="93" />为实数解<img file="FDA0000791824500000028.GIF" wi="86" he="74" />的验后均方差;(Q<sub>XX</sub>)<sub>ii</sub>为权逆阵中第i行第i列的元素;<img file="FDA0000791824500000029.GIF" wi="113" he="101" />则可根据自由度f=n‑u和置信水平(1‑α),从历元t分布的数值表中查取;采用上述方式,分别将r个元素中每个元素的整数候选值都挑出来,进行排列组合,共得到<img file="FDA00007918245000000210.GIF" wi="278" he="133" />个不同的整数组合,N<sub>i</sub>为<img file="FDA00007918245000000211.GIF" wi="86" he="71" />置信区间中候选整数值的个数,则N即为整周模糊度搜索空间;(3)根据陀螺测量信息,以动态过程载体方位角变化值为约束条件,压缩模糊度搜索范围;(4)基于压缩后的模糊度搜索范围,求解得到正确的模糊度。
地址 410073 湖南省长沙市砚瓦池正街47号中国人民解放军国防科学技术大学机电工程与自动化学院