[发明专利]一种进行星间单差模糊度固定的GNSS网数据处理方法

权利要求书

1.一种进行星间单差模糊度固定的GNSS网数据处理方法，其特征在于，包括以下步骤：

步骤一，对GNSS观测网的双频伪距和载波观测数据进行处理，获得所有未知参数的实数解；

步骤二，进行宽巷星间单差模糊度固定；

步骤三，进行消电离层组合星间单差模糊度固定；

步骤四，调整未知参数的解；调整未知参数的解的方法为：

首先利用下式调整未知参数向量的估值：

其中，为非钟差参数实数解，和N_xb为总法矩阵N的分块矩阵，它们有如下关系：

则星间单差模糊度固定之后的非钟差参数的解向量为对应的协方差矩阵由下式计算：

其中和分别为星间单差模糊度固定之后的非差模糊度估值及其协方差矩阵；

最后，利用第j个历元钟差参数的法方程信息及其与的关系计算得到该历元的钟差解：

2.根据权利要求1所述的一种进行星间单差模糊度固定的GNSS网数据处理方法，其特征在于，步骤一中，以双频消电离层组合非差伪距和载波相位作为观测量，采用最小二乘方法进行参数估计，未知参数至少包括双频消电离层组合模糊度参数、接收机和卫星的钟差参数。

3.根据权利要求2所述的一种进行星间单差模糊度固定的GNSS网数据处理方法，其特征在于，采用最小二乘方法进行参数估计需要进行若干次迭代，以确保参数估值的质量，每次迭代计算包括以下步骤：

第一步，逐历元构造构建法方程，对应历元j＝1,2,3,…的法方程表示为：

其中，为未知参数向量；为历元j的钟差参数，只与历元j的观测量有关；为除钟差之外的其他未知参数，存在于多个历元的观测量当中；为以分块矩阵表示的法矩阵，为法方程的右项，为了避免法方程奇异，需要选取一个接收机或卫星的钟作为参考钟，将其钟差参数约束于给定的值；

第二步，对每个历元的法方程进行参数约化，消去钟差参数得到约化后的法方程：

保存对应于钟差参数的法方程信息：N_C,j、N_CX,j和Z_C,j；

第三步，将所有历元约化后的法方程叠加，得到总法方程：

其中N为总法矩阵；Z为总法方程右项，求解总法方程得到除钟差参数之外的所有未知参数的解向量和对应的协方差矩阵Σ＝N^-1；

第四步，利用第二步保存的钟差参数的法方程信息和第三步得到的非钟差参数的估值计算每个历元j的钟差参数：

4.根据权利要求1所述的一种进行星间单差模糊度固定的GNSS网数据处理方法，其特征在于，步骤二中，进行宽巷星间单差模糊度固定包括以下步骤：

第一步、采用MW组合计算宽巷非差模糊度：

利用MW组合计算与每个消电离层组合非差模糊度参数对应的宽巷非差模糊度及其标准差的估值：和

其中为测站r对卫星s的某一连续无周跳数据段中第i个历元的MW组合观测量；n为该数据段的历元总数；

第二步、利用宽巷非差模糊度作为观测量计算卫星端的宽巷FCB：

第三步、进行宽巷星间单差模糊度固定；

利用第一步得到的非差宽巷模糊度构造出所有可能的宽巷星间单差模糊度，并利用第二步得到的卫星端宽巷FCB对其进行改正，计算公式如下：

其中，和分别为接收机r到卫星s和l的宽巷非差模糊度，和分别为卫星s和l的宽巷FCB，为接收机r对卫星s和l的宽巷星间单差模糊度；

采用下式计算将固定到最近整数的成功概率p_w：

其中为对应于的标准差；[·]表示取最近的整数；

将满足固定条件的宽巷星间单差模糊度固定到最近的整数并将该单差模糊度作为候选星间单差模糊度；否则，在后续步骤中不考虑该单差模糊度。