一种区域定位导航增强信息提取算法
【技术领域】
[0001] 本发明设及导航信号增强信息提取算法，具体为一种区域定位导航增强信息提取 算法，广泛适用于精密定位导航场景中。
【背景技术】
[0002] 全球卫星导航系统（GNSS，Global化vigationSatelliteSystem)快速发展的 近20年间，各个系统发射卫星数量不断增加，合理利用运些资源有助于提高GNSS的定位精 度。同时，需要实时精密定位的用户数量快速增长，且动态环境导致的多普勒频偏、城市建 筑物密集带来的多径效应和峡谷森林等复杂环境造成的遮挡作用大大影响了接收机定位 导航的精度和实时性，单纯依赖卫星信号已无法满足定位导航需求，因此为用户提供辅助 信息来增强定位导航性能具有重要的研究意义。
[0003] 目前，国内外提供的GNSS定位导航增强信息主要是从静态观测站获取数据，尚未 从动态用户接收机出发去研究，而动态用户却是GNSS提供服务的主要对象，也是拥有各类 原始信息的主要来源。基于海量动态历史GNSS数据的区域定位导航增强信息提取算法可 W解决静态数据分析的局限性。
【发明内容】
[0004] 本发明为了解决基于静态数据分析的定位导航增强信息提取算法具有局限性的 问题，提供了一种区域定位导航增强信息提取算法。
[0005] 本发明是采用如下的技术方案实现的：一种区域定位导航增强信息提取算法，包 括W下步骤：
[0006] 四个卫星系统历史中频数据经过卫星信号接收机捕获、跟踪、位同步和帖同步处 理，由卫星信号接收机输出各路卫星信号的伪距测量值和导航电文；
[0007] 四个卫星系统各自的权矩阵分别为
4,分别代表GPS、抓S、化0ANSS、Galileo卫星系统，目i为第j颗卫星的高度角；口！为单点 定位时观测卫星在天顶时的伪距观测单位权方差，则初始权矩阵P= (Pi，P2，P3，P4);
[0008] 不同系统n颗卫星的伪距测量值观测方程组化矩阵AP=AAx+e;式中Ax为 位置偏差，A表示系统矩阵，e表示系统噪声；
[0009] 对n颗卫星的伪距测量值观测方程组化矩阵利用初始权矩阵P进行加权，得到 PAP=PAAx+Pe，利用序贯最小二乘原理求解得到位置偏差AX= (ATpA)lATpAP， 并结合伪距测量值观测方程组化矩阵得到误差V=AAX-AP ;
[0010] 选择基于M残差的方差分量估计验后定权法来对误差进行迭代，使得误差V最 小并得到最佳权矩阵P。:基于M残差的方差分量估计验后定权法的公式为：=W，式中
[0011]
[001引S。=k1付(N(i)NiN( "Ni)-化2付(N(i)Ni)+rii，i=1，2, 3, 4，
[001引Sij=kitr(N(i)NiN(i)Ni)，i=l，2,3,4，j= 1，2,3,4,A^, =4^户乂，其中A康 示第i个卫星系统单独定位时的系统矩阵，ni表示第i个卫星系统的可见卫星数目， N=N1+N2+N3+N4,其中
第i个系统的第K次更新后的权矩阵 式中C为任一常数，迭代计算至<=<=< =< 时，误差V最小，此时得到
， 的权矩阵P/W就是最佳权矩阵P。；
[0014] 卫星信号的伪距测量值和导航电文通过最佳权矩阵P。进行加权融合，得到定位导 航增强信息；
[0015] 定位导航增强信息与相应的=维位置坐标进行对应，得到区域定位导航增强信 息，完成区域定位导航增强信息的提取。
[0016] 本发明在普通的方差分量估计法中加入了精密调控估计性能的M估计，提高了方 差分量估计法的精确度和抗粗差性能，从而可W有效的解决多系统组合定位时，对多系统 伪距进行有效加权的问题；同时，通过将定位导航增强信息与位置信息进行对应，为增强信 息增加区域属性，从而提高增强信息的精度。通过使用该种区域定位导航增强信息提取算 法，可W有效的提高卫星系统定位的精度和实时性，特别是在遮挡严重、多径效应严重、可 见卫星数目不足的恶劣环境中，在本算法提取出来的区域增强信息的辅助下，可W实现伪 距精密单点定位。
【附图说明】
[0017] 图1为卫星信号处理流程图。
[001引图2为区域定位导航增强信息提取流程图。
【具体实施方式】
[0019] 一种区域定位导航增强信息提取算法，包括W下步骤：
[0020] 四个卫星系统历史中频数据经过卫星信号接收机捕获、跟踪、位同步和帖同步处 理，由卫星信号接收机输出各路卫星信号的伪距测量值和导航电文；
[0021]四个卫星系统各自的权矩阵分别为
，i=l，2,3， 4,分别代表GPS、抓S、化0ANSS、Galileo卫星系统，目i为第j颗卫星的高度角；口>5为单点 定位时观测卫星在天顶时的伪距观测单位权方差，则初始权矩阵P= (Pi，P2，P3，P4);
[0022] 不同系统n颗卫星的伪距测量值观测方程组化矩阵AP=AAx+e;式中Ax为 位置偏差，A表示系统矩阵，e表示系统噪声；
[0023] 对n颗卫星的伪距测量值观测方程组化矩阵利用初始权矩阵P进行加权，得到 PAP=PAAX+Pe，利用序贯最小二乘原理求解得到位置偏差AX= (ATpA)lATpAP， 并结合伪距测量值观测方程组化矩阵得到误差V=AAX-AP ;
[0024] 选择基于M残差的方差分量估计验后定权法来对误差进行迭代，使得误差V最 小并得到最佳权矩阵P。:基于M残差的方差分量估计验后定权法的公式为：没=r.，式中
[0027]
，其中A康 示第i个卫星系统单独定位时的系统矩阵，rii表示第i个卫星系统的可见卫星数目， N=N1+N2+N3+N4,其中：
第i个系统的第K次更新后的权矩阵
式中C为任一常数，迭代计算至
时，误差V最小，此时得到 的权矩阵P/W就是最佳权矩阵P。；
[0028] 卫星信号的伪距测量值和导航电文通过最佳权矩阵P。进行加权融合，得到定位导 航增强信息；
[0029] 定位导航增强信息与相应的=维位置坐标进行对应，得到区域定位导航增强信 息，完成区域定位导航增强信息的提取。
[0030] 具体实施时，基于M残差的方差分量估计验后定权法中的参数
现对ki与k2中的参数A1,，A,,，入3,做如下推导：在使用基 于M残差的方差分量估计时，设有m类彼此独立的观测，第i类观测的个数为ni，相应 的单位权方差为<，权矩阵为Pi，函数模型与随机模型分别为：Vi=AlAX-Li(i=
正i. = (A.i，左，巧='塔(气，C2，...，C, ).i= 1，2,. . .，m，根据Huber给出的关于M估计的定义，不同类观测的解方程为：
，根据 M估计理论，可引入前文中的=个参数：
./ = ，分别是第i类观测的第一、二S多余参数。
[0031] 多系统历史中频数据需经历捕获、跟踪、位同步和帖同步四个处理过程才能得到 各路卫星信号的伪距和导航电文，运是由GNSS软件接收机完成的，在捕获阶段，接收机分 别从数字中频数据中捜索、捕获各颗可见卫星信号，并且从中获得当前运个卫星信号的载 波频率和C/A码相位的粗略估计值，在跟踪阶段，跟踪环路将通过码环和载波环来逐步精 细估计运两个信号参量，再经过位同步和帖同步之后输出各路卫星信号伪距和导航电文数 据比特。
[0032] 因为每一颗卫星对应着不相同的信号的多普勒频移、C/A码序列及其相位等一些 信号参量，所W必须对每颗卫星的信号分别进行独立的跟踪与处理。在典型接收机的基带 数字信号处理部分，常有N个图1所示的处理通道，中频数据需要从全部N个接收机信号 通道的入口输入到处理通道，然后各个通道各自输出其所跟踪的那颗卫星的伪距和导航电 文。
[0033] 得到海量伪距测量值和导航电文后，需要对海量历史GNSS伪距和导航电文进行 加权融合，再将加权融合后的定位导航增强信息与相应的=维位置坐标进行对应，来完成 区域定位导航增强信息的提取，因此海量历史GNSS伪距和导航电文的最佳权矩阵的确定 是区域定位导航增强信息提取的关键。本发明使用高度角先验定权和基于M残差的方差 分量估计验后定权法来确定最佳权矩阵，区域定位导航增强信息提取算法流程图如图2所 /J、- 〇
[0034] 由于导航电文中包含了钟差、轨道误差等重要信息，加权后的导航电文即为精度 较高的钟差、轨道误差校正值等等，也即本发明中所指的增强信息。在导航电文加权的过 程中，已经根据伪距观测方程对其进行了定位解算，也即得到了各个定位点的位置坐标，因 此将各个加权后的导航电文与相应的位置坐标对应起来，形成一个数值对（locationID， ANSI)，完成区域定位导航增强信息的提取。
【主权项】
1. 一种区域定位导航增强信息提取算法，其特征在于包括以下步骤： 四个卫星系统历史中频数据经过卫星信号接收机捕获、跟踪、位同步和帧同步处理，由 卫星信号接收机输出各路卫星信号的伪距测量值和导航电文； 四个卫星系统各自的权矩阵分别为：，1 = 1，2,3,4，分 别代表GPS、BDS、GLOANSS、Galileo卫星系统，Θ j为第j颗卫星的高度角；为单点定位 时观测卫星在天顶时的伪距观测单位权方差，则初始权矩阵P= (PubP3J4); 不同系统η颗卫星的伪距测量值观测方程组化矩阵Λ P =ΑΛχ+ε ;式中Λ X为位 置偏差，A表不系统矩阵，ε表不系统噪声； 对η颗卫星的伪距测量值观测方程组化矩阵利用初始权矩阵P进行加权，得到P Λ P =PA Λ χ+Ρ ε，利用序贯最小二乘原理求解得到位置偏差Λ X = (AtPA) 1AtP Λ P，并结合 伪距测量值观测方程组化矩阵得到误差V = A Λ X- Λ P ; 选择基于M残差的方差分量估计验后定权法来对误差进行迭代，使得误差V最小 并得到最佳权矩阵Pc:基于M残差的方差分量估计验后定权法的公式为：，式中Slj= k!tr(N( 1}ΝΧΝ( liN1), i = I, 2, 3, 4, j = I, 2, 3, 4,，其中 A1 表示第 i 个 卫星系统单独定位时的系统矩阵，Ii1表示第i个卫星系统的可见卫星数目，N = Ni+N2+N3+N4，其中： ，第i个系统的第K次更新后的权矩阵 式中c为 ， 任一常数，迭代计算至，误差V最小，此时得到的权矩阵P1(k)就是最佳 权矩阵P。； 卫星信号的伪距测量值和导航电文通过最佳权矩阵Pc进行加权融合，得到定位导航增 强信息； 定位导航增强信息与相应的三维位置坐标进行对应，得到区域定位导航增强信息，完 成区域定位导航增强信息的提取。
【专利摘要】本发明涉及导航信号增强信息提取算法，具体为一种区域定位导航增强信息提取算法，包括以下步骤：多系统历史中频数据经过卫星信号接收机处理，由卫星信号接收机输出各路卫星信号的伪距和导航电文；通过高度角先验定权确定初始权矩阵；初始权矩阵通过基于M残差的方差分量估计验后定权法来确定最佳权矩阵；各路卫星信号的伪距和导航电文与最佳权矩阵进行加权融合，得到定位导航增强信息；定位导航增强信息与三维坐标进行对应，完成区域定位导航增强信息的提取。本发明有效的提高卫星系统定位的精度和实时性，特别是在遮挡严重、多径效应严重、可见卫星数目不足的恶劣环境中，在本算法提取出来的区域增强信息的辅助下，可以实现伪距精密单点定位。
【IPC分类】G01S19/37
【公开号】CN105182379
【申请号】CN201510651325
【发明人】李灯熬, 吴刚, 赵菊敏
【申请人】太原理工大学
【公开日】2015年12月23日
【申请日】2015年10月10日