[发明专利]卫星导航系统GNSS定姿测量值的处理方法及GNSS定姿测量仪

权利要求书

1.一种卫星导航系统GNSS定姿测量值的处理方法，其特征在于，包括：

读取目标载体的载体姿态信息的历元初始值；

对所述历元初始值进行姿态量测处理，以获取所述载体姿态信息的历元量测值；

对前一历元最终值进行估值处理，以获取所述载体姿态信息的历元预报值；

对所述历元量测值和所述历元预报值进行抗差自适应融合滤波处理，以获取所述载体姿态信息的当前历元最终值。

2.根据权利要求1所述的方法，其特征在于，在读取目标载体的载体姿态信息的历元初始值之前，所述方法还包括：

实时读取所述目标载体的载体姿态信息的历元观测值；

对所述历元观测值进行估计处理，以获取载体姿态信息的所述历元初始值。

3.根据权利要求1所述的方法，其特征在于，对所述历元初始值进行姿态量测处理，以获取所述载体姿态信息的历元量测值包括：

通过GNSS静态相对定位算法来获取所述目标载体的多条基线在载体坐标系下的第一基线向量组k为历元号，BB为载体坐标系；

通过载波相对定位算法来获取所述目标载体的多条基线在地心坐标系下的第二基线向量组EE表示地心坐标系；

通过旋转矩阵将所述第二基线向量组转换成在当地水平坐标系下的第三基线向量组LL表示当地水平坐标系，所述第三基线向量组

根据旋转矩阵来获取所述第一基线向量组与所述第三基线向量组的函数关系：dXLLk=R312kdXBBk;]]>

通过泰勒级数展开至一阶来将非线性函数转换为线性函数AkδX~k-Lk=Vk,Pk,]]>其中，

Ak=∂R312∂αdXB,2k∂R312∂βdXB,2k∂R312∂γdXB,2k.........∂R312∂αdXB,ik∂R312∂βdXB,ik∂R312∂γdXB,ik.........∂R312∂αdXB,nk∂R312∂βdXB,nk∂R312∂γdXB,nk,]]>Lk=dXL,2k-B2kdXB,2k...dXL,ik-BikdXB,ik...dXL,nk-Bn-1kdXB,nk,]]>

Pk=ΣdXL,2k00.........0ΣdXL,ik0.........00ΣdXL,nk,]]>表示第i条基线k历元的解算精度，V_k为残差，δX~k=[AkTPkAk]-1[AkTP‾kLk],]]>Bik=R312k(α0,β0,γ0),]]>(α₀，β₀，γ₀)为所述历元初始值，(α，β，γ)表征载体姿态信息，k为历元号，i为基线编号；

通过最小二乘解算法对所述P_k进行求解，以获取所述历元量测值X~k=Xk0+δX~k.]]>