具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。

因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

下面结合附图，对本发明的一些实施方式作详细说明。在不冲突的情况下，下述的实施例及实施例中的特征可以相互组合。

某同步卫星每天24小时进行测轨，测轨数据包括地面测距主站单程测距、双程测距数据和副站四程测距数据以及各站气象数据，地面运控中心每天使用这些轨道观测数据完成该卫星的精密定轨。地面站设备改造与维护、设备切换、降雨等都会对同步卫星的外测数据造成影响，有可能会影响卫星的定轨精度，岗位值班人员每天进行同步卫星轨道确定任务时要根据外测数据质量的变化调整定轨策略。目前外测数据的总体质量情况没有有效的手段进行监视，基本还需要依靠专业技术人员进行复核，导致卫星的精密定轨需要过多的人工干预，既无法快速的获知同步卫星定轨处理全流程的运行状态，又导致人员工作强度大，不适应后续卫星平台在轨管理的现状。有鉴于此，本发明实施例提供了一种同步卫星外测数据的质量监控方法，用以缓解上文中所提出的技术问题。

实施例一

图1为本发明实施例提供的一种同步卫星外测数据的质量监控方法的流程图，如图1所示，该方法具体包括如下步骤：

步骤S102，接收地面测距主站发送的当前气象数据和实时测距数据。

具体的，本发明实施例提供的同步卫星外测数据的质量监控方法，主要是对实时测距数据的总体质量进行监控，实时测距数据包括：单程测距数据，双程测距数据，四程测距数据，本发明实施例提供的方法适用于对上述任意一种测距数据进行质量评估。为了能够降低气候对实时测距数据质量的影响，在对实时测距数据进行质量评估时，还需要将气象数据考虑进去，因此，除了需要接收地面测距主站发送的实时测距数据以外，还需要获取当前气象数据。

步骤S104，获取日常轨道精密定轨与预报输出的卫星星历文件和测距系统误差。

步骤S106，基于实时测距数据和卫星星历文件确定实时测距数据对应的理论测距数据。

测距系统误差为日常轨道精密定轨输出的参量，用于对实时测距数据进行修正。为了能够计算出实时测距数据对应的理论测距数据，方法运行时还需要获取日常轨道预报输出的卫星星历文件，基于卫星星历文件能够确定卫星在某一时刻下的地心位置矢量，且接收到的地面测距主站发送的实时测距数据中，每个测距数据均带有测距时标，因此，基于实时测距数据和卫星星历文件可以进行理论测距数据的计算。下文中将对理论测距数据的计算方法进行详细的描述。

步骤S108，基于当前气象数据、测距系统误差、实时测距数据和理论测距数据确定实时测距数据的随机误差。

在计算得到理论测距数据后，利用当前气象数据能够计算出大气折射修正量，测距系统误差也能对实时测距数据进行一定量的修正，对实时测距数据进行上述两种修正后，再与理论测距数据相减，即可确定出实时测距数据的随机误差；其中，随机误差的大小用于表征实时测距数据的数据质量，显然，随机误差越大表示实时测距数据与理论测距数据之间的差距较大，表明实时测距数据的数据质量较差或卫星轨道出现变化；反之，如果随机误差较小，表明实时测距误差的数据质量较好。对于随机误差的大小判断，本发明实施例不对其进行具体限定，用户可以根据实际需求设定相应的取值范围，从而快速判断出实时测距数据的质量是否合格。

本发明提供的同步卫星外测数据的质量监控方法，在获取到实时测距数据后，首先利用实时测距数据和日常轨道精密定轨与预报输出的卫星星历文件确定出实时测距数据对应的理论测距数据，然后再结合当前气象数据和日常精密定轨输出的测距系统误差确定出实时测距数据的随机误差。本发明方法将轨道动力学模型和观测模型等先验信息用于卫星外测数据的实时处理，提高了外测数据质量异常监测的敏感度，基于实时测距数据的随机误差变化可实现对卫星外测数据的实时监测评估，为及时掌握外测数据质量，调整定轨策略提供了数据支撑。

上文中对本发明实施例提供的同步卫星外测数据的质量监控方法进行了简要的描述，下面对其中涉及的相关方法步骤进行具体介绍。

在一个可选的实施方式中，如图2所示，上述步骤S106，基于实时测距数据和卫星星历文件确定实时测距数据对应的理论测距数据，具体包括如下步骤：

步骤S1061，获取实时测距数据的测距时标。

步骤S1062，基于测距时标和卫星星历文件确定实时测距数据对应的理论测距参考值。

具体的，通过上文中的描述可知，实时测距数据均携带测距时标，但实时测距数据所携带的测距时标并不是测距伪码信号的实际发射时刻，在本发明实施例中，计算理论测距参考值时暂不考虑测距伪码信号在空间中的传输时延，也即，认为测距伪码信号的发射时刻与接收测距伪码信号的时刻相同，对于四程测距数据，地面测距副站的接收时刻也与上述发射时刻相同，因此利用上述特征，可以计算出每一种实时测距数据(单程、双程、四程)对应的理论测距参考值。

步骤S1063，基于理论测距参考值和预设模糊周期对测距时标进行修正，得到测距伪码信号的实际发射时刻。

图3示出了一种伪码测距测量原理示意图，如图3所示，测量的测距数据R是通过计算最近的收“1”时间与最近的发“1”时间差得到的，同时测距时标打在最近的发“1”时标上，如果测距数据跨6个模糊度，当收端收到第1个脉冲时，发“1”端已经到了第7个脉冲，因此测距时标打在第7个脉冲上，此时测出的R值为未解模糊度的测距值，测距时标t也为未修正的时标，在这种测量体制下，利用外测数据进行卫星轨道确定时，必须先对获得的测距数据和对应的时标进行修正。

在本发明实施例中，计算测距伪码信号的实际发射时刻时需使用理论测距参考值与预设模糊周期，单程测距数据、双程测距数据和四程测距数据的预设模糊周期通过设备研制方获取。在计算单程、双程以及四程测距数据的测距伪码信号的实际发射时刻时，需要使用相应的理论测距参考值和预设模糊周期。下文中将对具体计算方法进行详细介绍。

步骤S1064，基于实际发射时刻和卫星星历文件确定实时测距数据对应的理论测距数据。

具体的，在确定了测距伪码信号的实际发射时刻之后，利用卫星星历文件即可确定出卫星在测距伪码信号的实际发射时刻的地心位置矢量，针对单程测距数据或双程测距数据，只需要再与地面测距主站的地心位置矢量相结合即可确定出单程测距数据和双程测距数据对应的理论测距数据；而对于四程测距数据，由于测距过程还涉及到地面测距副站，因此，在确定出卫星在测距伪码信号的实际发射时刻的地心位置矢量之后，还需要再与地面测距主站的地心位置矢量以及地面测距副站的地心位置矢量相结合才能确定出四程测距数据对应的理论测距数据。

上文中对如何基于实时测距数据和卫星星历文件确定实时测距数据对应的理论测距数据的过程进行了简要描述，下面对其中如何确定实时测距数据对应的理论测距参考值的方法进行具体介绍。

在一个可选的实施方式中，上述步骤S1062，基于测距时标和卫星星历文件确定实时测距数据对应的理论测距参考值，具体包括如下步骤：

步骤S10621，基于测距时标和卫星星历文件确定卫星在测距时标对应时刻的地心位置矢量。

具体的，卫星星历文件是记载卫星运行轨道数据的文件，利用卫星星历文件并给定测距时标，即可确定出测距时标对应时刻下卫星的地心位置矢量。

具体的，本发明实施例使用星历文件进行插值计算求解卫星的地心位置矢量，插值方法采用拉格朗日多项式插值，设函数在x₁,x₂,…,x_n(相当于本发明实施例中的卫星的地心位置矢量对应的时间)各点上的值分别为y₁,y₂,…,y_n(相当于本发明实施例中的卫星的地心位置矢量)，拉格朗日多项式插值的经典公式如下：利用上述插值公式，即可确定出卫星在测距时标对应时刻的地心位置矢量。一般地，卫星星历插值选择7-8阶多项式进行插值，可选的，n取值为8。

若实时测距数据为单程测距数据或双程测距数据，则执行步骤S10622；若实时测距数据为四程测距数据，则执行步骤S10623。

步骤S10622，基于卫星在测距时标对应时刻的地心位置矢量和地面测距主站的地心位置矢量确定单程测距数据对应的理论测距参考值或双程测距数据对应的理论测距参考值。

由于单程测距数据和双程测距数据所涉及到的对象均为地面测距主站和卫星，在忽略测距伪码信号在空间的传输时延的前提下，可认为双程测距数据对应的理论测距参考值为单程测距数据对应的理论参考值的2倍，也即其中，ρ′₂表示双程测距数据对应的理论测距参考值，ρ′₁表示单程测距数据对应的理论参考值，表示卫星在测距时标对应时刻t₀的地心位置矢量，表示地面测距主站在测距时标对应时刻t₀的地心位置矢量。

步骤S10623，基于卫星在测距时标对应时刻的地心位置矢量、地面测距主站的地心位置矢量和地面测距副站的地心位置矢量确定四程测距数据对应的理论测距参考值。

由于四程测距数据所涉及到的对象为地面测距主站、卫星和地面测距副站，在忽略测距伪码信号在空间的传输时延的前提下，四程测距数据对应的理论测距参考值可利用算式来进行计算，其中，表示地面测距副站在测距时标对应时刻t₀的地心位置矢量(一般不随时间改变)。

上文中对如何计算实时测距数据对应的理论测距参考值的过程进行了详细的描述，下面对如何确定实时测距数据对应的理论测距数据进行具体介绍。

在一个可选的实施方式中，上述步骤S1064，基于实际发射时刻和卫星星历文件确定实时测距数据对应的理论测距数据，具体包括如下步骤：

步骤S10641，基于实际发射时刻和卫星星历文件确定卫星在实际发射时刻的地心位置矢量。

具体的，同上述步骤S10621，利用拉格朗日多项式插值公式即可求得卫星在测距伪码信号的实际发射时刻下的地心位置矢量。

若实时测距数据为单程测距数据或双程测距数据，则执行步骤S10642-步骤S10643；若实时测距数据为四程测距数据，则执行步骤S10644-步骤S10645。

步骤S10642，基于卫星在实际发射时刻的地心位置矢量和地面测距主站的地心位置矢量确定单程测距数据的测距时延或双程测距数据的测距时延。

步骤S10643，基于单程测距数据的测距时延和光速确定单程测距数据对应的理论测距数据；或者，基于双程测距数据的测距时延和光速确定双程测距数据对应的理论测距数据。

具体的，在计算任意一种测距数据(单程或双程或四程)对应的理论测距数据时，首先要计算该测距数据的测距时延，然后将该测距数据的测距时延与光速相乘即可求得该测距数据对应的理论测距数据。例如，在计算单程测距数据对应的理论测距数据时，首先要计算单程测距数据的测距时延，然后将单程测距数据的测距时延与光速相乘即可求得单程测距数据对应的理论测距数据。

单程测距数据的测距时延为测距伪码信号由地面测距主站发送至卫星的时延，单程测距数据的测距时延可使用算式进行表示，其中，τ₁表示单程测距数据的测距时延，t₂表示卫星接收到测距伪码信号的时刻，t₁表示地面测距主站发送测距伪码信号的时刻，表示t₂时刻卫星的地心位置矢量，表示t₁时刻地面测距主站的地心位置矢量。当前已求得了测距伪码信号的发射时刻t₁，以及卫星在测距伪码信号的实际发射时刻下的地心位置矢量地面站的精确坐标可知，可以通过坐标转换得到相应时刻的地心位置矢量，所以卫星接收到测距伪码信号的时刻t₂以及t₂时刻卫星的地心位置矢量为未知数。

本发明实施例使用以下方法来求取任意一段“测距信号的接收端位置随时间变化”的转发区间的测距时延，假设测距伪码信号的某一段转发区间为从目标A发送到目标B，且测距信号从目标A出发的时刻为t_s，那么从目标A转发到目标B的测距信号传输时延其中，目标B在t_s+τ时刻的地心位置矢量为未知数，目标A在t_s时刻的地心位置矢量为已知数，c表示光速，τ可利用如下算法来进行求解：

1)计算并判断Δτ＝τ″-τ′是否小于10^-12；其中，首次计算时，令τ′＝0；

2)若不小于，则将当前计算得到的τ″作为下一次迭代计算时的τ′，并利用步骤1)中的算式迭代计算下一代的τ″，直至Δτ收敛(Δτ<10^-12)，并将满足收敛条件的τ″作为目标A转发到目标B的测距信号传输时延τ。

利用以上方法即可求得单程测距数据的测距时延τ₁。然后将τ₁与光速相乘即可求得单程测距数据对应的理论测距数据。根据可计算出卫星接收到测距伪码信号的时刻t₂，进一步结合卫星星历文件可计算出t₂时刻卫星的地心位置矢量

图4为本发明实施例提供的一种双程测距和四程测距的示意图，图4中，实线表示双程测距时测距伪码信号的流向，虚线表示四程测距时测距伪码信号的流向。双程测距数据的测距时延应表示为τ₁+τ₂，其中τ₁表示测距伪码信号由地面测距主站发送至卫星的时延，τ₂表示卫星将测距伪码信号返回至地面测距主站的时延，上文中已经计算出τ₁，t₂和那么利用上文中所描述的求取任意一段测距信号的接收端位置随时间变化的转发区间的测距时延的方法，可以求出τ₂和测距伪码信号返回至地面测距主站的时刻t₃，进而确定出双程测距数据的测距时延τ₁+τ₂。然后将τ₁+τ₂与光速相乘即可求得双程测距数据对应的理论测距数据。

步骤S10644，基于卫星在实际发射时刻的地心位置矢量、地面测距主站的地心位置矢量和地面测距副站的地心位置矢量确定四程测距数据的测距时延。

步骤S10645，基于四程测距数据的测距时延和光速确定四程测距数据对应的理论测距数据。

四程测距数据的测距时延应表示为τ₁+τ₃+τ₄+τ₅，其中τ₁表示测距伪码信号由地面测距主站发送至卫星的时延，τ₃表示卫星将测距伪码信号转发至地面测距副站的时延，τ₄表示地面测距副站将测距伪码信号转发至卫星的时延，τ₅表示卫星将测距伪码信号转发至地面测距主站的时延。

上文中已经计算出τ₁，t₂和τ₃的算式中未知参量为测距伪码信号转发至地面测距副站的时刻t₄，以及t₄时刻地面测距副站的地心位置矢量利用上文中所描述的求取任意一段测距信号的接收端位置随时间变化的转发区间的测距时延的方法，可以得出τ₃和t₄。

进一步的，τ₄的算式中未知参量为测距伪码信号转发至卫星的时刻t₅，以及t₅时刻卫星的地心位置矢量利用上文中所描述的求取任意一段测距信号的接收端位置随时间变化的转发区间的测距时延的方法，可以得出τ₄，t₅，进一步结合卫星星历文件可计算出

进一步的，τ₅的算式中未知参量为测距伪码信号转发至地面测距主站的时刻t₆，以及t₆时刻地面测距主站的地心位置矢量利用上文中所描述的求取任意一段测距信号的接收端位置随时间变化的转发区间的测距时延的方法，可以得出τ₅和t₆，进而确定出四程测距数据的测距时延τ₁+τ₃+τ₄+τ₅，然后将τ₁+τ₃+τ₄+τ₅与光速相乘即可求得四程测距数据对应的理论测距数据。

上文中对如何确定出任一种实时测距数据对应的理论测距数据的方法进行了详细的描述，下面对如何计算实时测距数据的随机误差的方法进行具体描述。

在一个可选的实施方式中，上述步骤S108，基于当前气象数据、测距系统误差、实时测距数据和理论测距数据确定实时测距数据的随机误差，具体包括如下步骤：

步骤S1081，对当前气象数据和实时测距数据进行野值辨识与剔除，得到目标气象数据和目标测距数据。

卫星外测数据为地面测距主站实时观测，要快速对数据予以处理，辨识和剔除数据包含的野值，必须建立既节省内存，又确保精度和实时性要求的外测数据野值剔除算法。本发明实施例中利用滑动数据拟合方法，不断改进拟合系数，提高拟合精度，再用函数通过预报以检验观测数据的合理性。

具体的，令y_j(j＝1～n)表示某时刻外测数据观测量，为了不失一般性，认为y_j是时间t的二阶函数：因为测量是等间隔的，即t_j＝j*h，其h为采样间隔(常值)，j为采样点点序。可以根据最小二乘原理，确定合适的拟合系数，以使将求得方差E最小时称为的最佳二乘拟合逼近。

对于点序j观测数据y_j，分别取其前n个数据：y_j-n,…y_j-1，算式两侧分别对多项式拟合系数a₀,a₁,a₂求导，并令其为零，可得系数的拟合值为： 其中，α＝2ⁿ⁺¹，

由此便可确定出观测数据的拟合函数

这样就可以利用拟合值检验原观测值y_j是否失真，若失真，则利用拟合值予以修复，以上是以第j点数据为例说明利用拟合法修复数据的具体思想。

由于数据观测是连续的，即由于目标在不断运动，拟合函数需根据运动变化予以更新，才能提高野值辨识和剔除的效率。因此，利用新观测的数据，辅助以下滑动方法，不断改进系数拟合值提高拟合精度，再用该函数通过预报以检验观测数据的合理性。

通过上文中的描述可知，外测数据野值辨识与剔除的流程如下：

(1)构造n个外测数据缓冲区；

(2)利用缓冲区累计的n个数据y₁,…y_n构造拟合函数

(3)对于利用步骤(2)求得新观测数据y_n+1的预报值

(4)判断是否成立，其中，σ₂表示野值判门限值，一般取17倍的设备观测精度指标值；

(5)若不等式成立，则说明y_n+1为正常数据点，且缓冲区利用先进先出原则去除第一点y₁，补充y_n+1到缓冲区以重复以上过程；

(6)若不等式不成立，则利用代替y_n+1，重复上述步骤(1)-(4)。

在判断过程中，如果连续若干个(例如5个)新的观测数据均不能让步骤(4)中的不等式成立，则将缓冲区中的数据丢弃，重新选取并构造新的拟合函数。

因此，在计算随机误差之前，首先需要利用上述方法对当前气象数据和实时测距数据进行野值辨识与剔除，得到目标气象数据和目标测距数据。

步骤S1082，基于目标气象数据确定测距数据对应的大气折射修正量。

气象数据是修正观测数据折射误差的基本依据，因此，可以根据地面测距主站发送的气象数据实时对测距数据进行实时修正，在本发明实施例中，具体为基于目标气象数据计算出测距数据对应的大气折射修正量，以供后续步骤中计算实时测距数据的随机误差。

步骤S1083，基于理论测距参考值和预设模糊周期对目标测距数据进行修正，得到修正后的测距数据。

上文中介绍了对测距时标进行修正的过程，与时标修正的原理类似，对于目标测距数据也需要进行修正，具体的，本发明实施例使用算式来计算修正后的测距数据，其中，ρ″_c表示修正后的测距数据，ρ_c表示实时测距数据，ρ′_c表示实时测距数据对应的理论测距参考值，c表示光速，ΔT表示实时测距数据对应的预设模糊周期。在对单程、双程以及四程测距数据进行修正时，需要使用相应的理论测距参考值和相应的预设模糊周期。

步骤S1084，基于大气折射修正量、测距系统误差、修正后的测距数据和理论测距数据确定实时测距数据的随机误差。

在计算出大气折射修正量、测距系统误差、修正后的测距数据和理论测距数据之后，本发明实施例采用算式来计算实时测距数据的随机误差，其中，Δσ表示实时测距数据的随机误差，ρ″_c表示修正后的测距数据，表示实时测距数据对应的理论测距数据，Δρ_trop表示大气折射修正量，Δρ_s表示测距系统误差。

对于一个跟踪弧段内的数据，还可以计算随机误差均方差来进一步评估外测数据质量，随机误差均方差的估计值为其中，M表示跟踪弧段内测量数据的总数。

在一个可选的实施方式中，上述步骤S1063，基于理论测距参考值和预设模糊周期对测距时标进行修正，得到测距伪码信号的实际发射时刻，具体包括如下内容：

基于算式t′＝t-N·ΔT对测距时标进行修正，得到测距伪码信号的实际发射时刻；其中，t′表示测距伪码信号的实际发射时刻，t表示测距时标对应时刻，N表示模糊数，N＝|ρ′_c/Δρ|，ρ′_c表示理论测距参考值，Δρ＝ΔT*c，ΔT表示预设模糊周期，c表示光速。

具体的，本发明实施例利用算式计算测距伪码信号的实际发射时刻，Δρ＝ΔT*c可理解为模糊距离，其中，单程测距数据、双程测距数据和四程测距数据对应的预设模糊周期跟设备有关。在计算单程、双程以及四程测距数据的测距伪码信号的实际发射时刻时，需要使用相应的理论测距参考值和相应的预设模糊周期。

在一个可选的实施方式中，上述步骤S1082，基于目标气象数据确定测距数据对应的大气折射修正量，具体包括如下内容：

利用算式确定测距数据对应的大气折射修正量，其中，ΔR表示大气折射修正量，C表示大气折射率指数，E_N表示理论测站仰角，表示地面折射率，T表示绝对温度，P表示地面大气压，P_e表示地面水气压。

为简化气象观测，综合考虑可将大气折射率指数C取统计平均值1.4142E-4(1/米)，绝对温度T＝T₀+T′，其中，T₀＝273.15(℃)，T′表示地面温度，地面水气压U表示地面相对湿度。

综上所述，本发明实施例提供了一种高精度的同步卫星外测数据的质量监控方法，该方法将轨道动力学模型(使用了卫星星历)和观测模型(结合了测距系统误差)等先验信息用于同步卫星外测数据的实时处理，提高了外测数据质量异常监测的敏感度，基于实时测距数据的随机误差变化可实现对卫星外测数据的实时监测评估，为及时掌握外测数据质量，调整定轨策略提供了数据支撑。

图5示出了多站双程和四程测距数据的质量监控界面，图5中左侧视图表示地面测距主站发送的经模糊距离和时标修正、以及野值辨识与剔除后的原始观测数据，显示的数字为当前观测时刻的测量值。图5中右侧视图为随机误差，显示的数字为该弧段统计的随机误差均值。从图5中可以看出，副站2失锁，外测数据出现了异常。

实施例二

本发明实施例还提供了一种同步卫星外测数据的质量监控装置，该同步卫星外测数据的质量监控装置主要用于执行上述实施例一所提供的卫星外测数据的质量监控方法，以下对本发明实施例提供的同步卫星外测数据的质量监控装置做具体介绍。

图6是本发明实施例提供的一种同步卫星外测数据的质量监控装置的功能模块图，如图6所示，该装置主要包括：接收模块10，获取模块20，第一确定模块30，第二确定模块40，其中：

接收模块10，用于接收地面测距主站发送的当前气象数据和实时测距数据。

获取模块20，用于获取日常轨道精密定轨与预报输出的卫星星历文件和测距系统误差。

第一确定模块30，用于基于实时测距数据和卫星星历文件确定实时测距数据对应的理论测距数据。

第二确定模块40，用于基于当前气象数据、测距系统误差、实时测距数据和理论测距数据确定实时测距数据的随机误差；其中，随机误差的大小用于表征实时测距数据的数据质量。

本发明提供的同步卫星外测数据的质量监控装置，包括：接收模块10，用于接收地面测距主站发送的当前气象数据和实时测距数据。获取模块20，用于获取日常轨道精密定轨与预报输出的卫星星历文件和测距系统误差。第一确定模块30，用于基于实时测距数据和卫星星历文件确定实时测距数据对应的理论测距数据。第二确定模块40，用于基于当前气象数据、测距系统误差、实时测距数据和理论测距数据确定实时测距数据的随机误差；其中，随机误差的大小用于表征实时测距数据的数据质量。本发明装置将轨道动力学模型和观测模型等先验信息用于同步卫星外测数据的实时处理，提高了外测数据质量异常监测的敏感度，基于实时测距数据的随机误差变化可实现对同步卫星外测数据的实时监测评估，为及时掌握外测数据质量，调整定轨策略提供了数据支撑。

可选的，第一确定模块30包括：

获取单元，用于获取实时测距数据的测距时标。

第一确定单元，用于基于测距时标和卫星星历文件确定实时测距数据对应的理论测距参考值。

第一修正单元，用于基于理论测距参考值和预设模糊周期对测距时标进行修正，得到测距伪码信号的实际发射时刻。

第二确定单元，用于基于实际发射时刻和卫星星历文件确定实时测距数据对应的理论测距数据。

可选的，第一确定单元具体用于：

基于测距时标和卫星星历文件确定卫星在测距时标对应时刻的地心位置矢量。

若实时测距数据为单程测距数据或双程测距数据，则基于卫星在测距时标对应时刻的地心位置矢量和地面测距主站的地心位置矢量确定单程测距数据对应的理论测距参考值或双程测距数据对应的理论测距参考值。

若实时测距数据为四程测距数据，则基于卫星在测距时标对应时刻的地心位置矢量、地面测距主站的地心位置矢量和地面测距副站的地心位置矢量确定四程测距数据对应的理论测距参考值。

可选的，第二确定单元具体用于：

基于实际发射时刻和卫星星历文件确定卫星在实际发射时刻的地心位置矢量。

若实时测距数据为单程测距数据或双程测距数据，则基于卫星在实际发射时刻的地心位置矢量和地面测距主站的地心位置矢量确定单程测距数据的测距时延或双程测距数据的测距时延。

基于单程测距数据的测距时延和光速确定单程测距数据对应的理论测距数据；或者，基于双程测距数据的测距时延和光速确定双程测距数据对应的理论测距数据。

若实时测距数据为四程测距数据，则基于卫星在实际发射时刻的地心位置矢量、地面测距主站的地心位置矢量和地面测距副站的地心位置矢量确定四程测距数据的测距时延。

基于四程测距数据的测距时延和光速确定四程测距数据对应的理论测距数据。

可选的，第二确定模块40包括：

野值辨识与剔除单元，用于对当前气象数据和实时测距数据进行野值辨识与剔除，得到目标气象数据和目标测距数据。

第三确定单元，用于基于目标气象数据确定测距数据对应的大气折射修正量。

第二修正单元，用于基于理论测距参考值和预设模糊周期对目标测距数据进行修正，得到修正后的测距数据。

第三确定单元，用于基于大气折射修正量、测距系统误差、修正后的测距数据和理论测距数据确定实时测距数据的随机误差。

可选的，第一修正单元具体用于：

基于算式t′＝t-N·ΔT对测距时标进行修正，得到测距伪码信号的实际发射时刻；其中，t′表示测距伪码信号的实际发射时刻，t表示测距时标对应时刻，N表示模糊数，N＝|ρ′_c/Δρ|，ρ′_c表示理论测距参考值，Δρ＝ΔT*c，ΔT表示预设模糊周期，c表示光速。

可选的，第三确定单元具体用于：

利用算式确定测距数据对应的大气折射修正量，其中，ΔR表示大气折射修正量，C表示大气折射率指数，E_N表示理论测站仰角，表示地面折射率，T表示绝对温度，P表示地面大气压，P_e表示地面水气压。

实施例三

参见图7，本发明实施例提供了一种电子设备，该电子设备包括：处理器60，存储器61，总线62和通信接口63，所述处理器60、通信接口63和存储器61通过总线62连接；处理器60用于执行存储器61中存储的可执行模块，例如计算机程序。

其中，存储器61可能包含高速随机存取存储器(RAM，RandomAccessMemory)，也可能还包括非不稳定的存储器(non-volatile memory)，例如至少一个磁盘存储器。通过至少一个通信接口63(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接，可以使用互联网，广域网，本地网，城域网等。

总线62可以是ISA总线、PCI总线或EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图7中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。

其中，存储器61用于存储程序，所述处理器60在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的流过程定义的装置所执行的方法可以应用于处理器60中，或者由处理器60实现。

处理器60可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器60中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器60可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(Digital SignalProcessing，简称DSP)、专用集成电路(Application Specific Integrated Circuit，简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器61，处理器60读取存储器61中的信息，结合其硬件完成上述方法的步骤。

本发明实施例所提供的一种同步卫星外测数据的质量监控方法、装置和电子设备的计算机程序产品，包括存储了处理器可执行的非易失的程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行前面方法实施例中所述的方法，具体实现可参见方法实施例，在此不再赘述。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。

所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。

在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，或者是该发明产品使用时惯常摆放的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。

此外，术语“水平”、“竖直”、“悬垂”等术语并不表示要求部件绝对水平或悬垂，而是可以稍微倾斜。如“水平”仅仅是指其方向相对“竖直”而言更加水平，并不是表示该结构一定要完全水平，而是可以稍微倾斜。

在本发明的描述中，还需要说明的是，除非另有明确的规定和限定，术语“设置”、“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。