一种rt-ppp星历改正数据优选算法
【技术领域】
[0001]本发明涉及卫星导航技术领域，具体涉及到一种RT-ppp星历改正数据优选算法。
【背景技术】
[0002] 在RT-ppp定位过程中，IGU星历特别是其卫星钟差精度往往成为影响精度和收敛 时间的关键性因素。鉴于目前尚不存在比较好的星历质量评估定权方法，提出此解决方案。 IGU的预报星历和观测星历时间间隔仅仅为3小时，但是其轨道精度尤其是卫星钟种差精度 却相差甚远，不同型号GPS卫星的质量、形状及光压模型都存在差异，导致其预报轨道精度 不同，从而造成IGU的预报轨道短时间内的质量具有非常强的相关性。
【发明内容】
[0003] 针对现有技术中RT-PPP定位存在上述的不足，本发明提供了一种RT-PPP星历改正 数据优选算法，能够实时为ppp解算提供经过自动轨道质量控制的星历，提高了实时ppp的 定位精度。
[0004] 本发明的技术方案为，一种RT-PPP星历改正数据优选算法，包括以下步骤：
[0005] 连续采集IGU精密星历的数据；
[0006] 根据采集的数据对IGU观测星历和其对应此时刻的IGU预报星历进行求差；
[0007]根据IGU观测星历和其对应此时刻的IGU预报星历的差值进行第一次定权并得到 权阵；
[0008] 根据上述权阵通过轨道松弛法进行参数约束并进行伪距定位解算。
[0009] 上述的优选算法，其中，所述根据采集的数据对IGU观测星历和其对应此时刻的 IGU预报星历进行求差的步骤具体为，将IGU观测轨道与对应的时刻之前储存下来的预报轨 道进行比较，并且求差，得到轨道参数差值 Vl，其中包括归算到单位为米的卫星钟种差。
[0010] 上述的优选算法，其中，所述根据ICTJ观测星历和其对应此时刻的IGU预报星历的 差值进行第一次定权并得到权阵的步骤具体为，取定标准差〇〇，则初始权阵P = σο/V;伪距观 测方程为:L=A*X+1，Ρ，其中L为伪距观测值向量，Α为设计矩阵，X是坐标位置，1为L-L〇。
[0011] 上述的优选算法，其中，所述根据上述权阵通过轨道松弛法进行参数约束并进行 伪距定位解算的步骤包括：
[0012] 将权阵应用于实时定位的最小二乘定位方程，并最终求得方差阵；
[0013] 将方差阵通过轨道松弛法进行参数约束并对有约束的轨道进行伪距定位解算。
[0014] 上述的优选算法，其中，所述根据上述权阵通过轨道松弛法进行参数约束并进行 伪距定位解算的步骤中，轨道松弛约束条件：
其中〇sxl后为验标准差，通过解伪距 观测方程得到，轨道参数Sxi和k对应t分布，取定k = 2.5达到98.8%的置信度。
[0015] 上述的优选算法，其中，若所述轨道松弛约束条件满足：
取定k = 2.5达到 98.8 %的置信度，则更新标准差为：
[0016]上述的优选算法，其中，将求得的新标准差归算到权阵中重新进行伪距定位最小 二乘计算直到轨道参数均满足轨道松弛约束条件为止，其中所述的轨道松弛约束条件为
，其中σδχ1后为验标准差，通过解伪距观测方程得到，轨道参数δ χ1和k对应t分布，取 定k = 2.5达到98.8 %的置信度。
[0017]本发明提供的一种RT-PPP星历改正数据优选算法，包括以下步骤:连续采集I⑶精 密星历的数据;根据采集的数据对IGU观测星历和其对应此时刻的IGU预报星历进行求差； 根据IGU观测星历和其对应此时刻的IGU预报星历的差值进行第一次定权并得到权阵;本发 明适用于RT-PPP定位，能够实时为PPP解算提供经过自动轨道质量控制的星历，提高了实时 PPP的定位精度，以及利用短时间内预报星历质量的强相关性进行星历权重判定，为轨道松 弛提供权重初始值，并且无需为接收机增加新硬件或者服务器，仅仅需要在软件层面上可 以更改现有的RT-PPP接收机的自动化轨道质量控制。
【附图说明】
[0018] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明及其特征、外 形和优点将会变得更明显。在全部附图中相同的标记指示相同的部分。并未刻意按照比例 绘制附图，重点在于示出本发明的主旨。
[0019] 图1为本发明提供的一种RT-PPP星历改正数据优选算法的示意图。
【具体实施方式】
[0020] 在下文的描述中，给出了大量具体的细节以便提供对本发明更为彻底的理解。然 而，对于本领域技术人员而言显而易见的是，本发明可以无需一个或多个这些细节而得以 实施。在其他的例子中，为了避免与本发明发生混淆，对于本领域公知的一些技术特征未进 行描述。
[0021] 为了彻底理解本发明，将在下列的描述中提出详细的步骤以及详细的结构，以便 阐释本发明的技术方案。本发明的较佳实施例详细描述如下，然而除了这些详细描述外，本 发明还可以具有其他实施方式。
[0022] 参照图1所示，本发明提供了一种RT-PPP星历改正数据优选算法，包括以下步骤： [0023]步骤S1:连续采集IGU精密星历的数据，不断采集IGU超快精密星历，并将数据用于 后续的计算。
[0024]步骤S2:根据采集的数据对IGU观测星历和其对应此时刻的IGU预报星历进行求 差，将IGU观测轨道与对应的时刻之前储存下来的预报轨道进行比较，并且求差，得到轨道 参数差值^，其中包括归算到单位为米的卫星钟种差。
[0025]步骤S3:根据IGU观测星历和其对应此时刻的IGU预报星历的差值进行第一次定权 并得到权阵，具体为，取定标准差〇〇,则初始权阵Ρ = σ〇/ν;伪距观测方程为:L = A*X+1，P，其 中L为伪距观测值向量，A为设计矩阵，X是坐标位置，1为L-Lo,也就是说对最新的ICTJ观测星 历和其对应此时刻的IGU预报星历进行求差，并根据此差值进行第一次定权，得到权阵Po。
[0026] 步骤S4:根据上述权阵通过轨道松弛法进行参数约束并进行伪距定位解算