[发明专利]一种多分析中心精密站坐标产品综合方法及系统

权利要求书

1.一种多分析中心精密站坐标产品综合方法，其特征在于，所述方法包括：

S1、对多个分析中心提供的GNSS站坐标产品文件进行预处理，利用分析中心提供的先验信息对预处理得到的原始法方程系统进行初步先验约束消除，利用相似变换矩阵消除残余约束；

S2、采用非潮汐海洋、大气和水文负载产品对目标框架进行非线性改正，将消除先验约束之后的分析中心站坐标产品对齐到目标框架，并计算目标框架下的分析中心站坐标解；所述步骤S2具体包括如下分步骤：

计算站点在当前历元t的非潮汐海洋、大气和水文负载形变，非潮汐海洋、大气和水文负载形变记为；

在目标框架中引入非潮汐海洋、大气和水文负载，计算目标框架在当前历元的坐标，计算公式为：

；

t和t₀分别表示当前历元和参考历元；表示目标框架解在当前历元的位置；表示目标框架解在参考历元的位置；表示目标框架解在参考历元的速度；

计算分析中心站坐标产品中核心站点目标框架解的坐标值；

对分析中心无约束的法方程系统添加统一约束：

；

其中，是由核心站坐标构成的约束矩阵，为核心站坐标构成的相似变换矩阵，为相应的转置矩阵，为目标框架下的分析中心站坐标解，s表示分析中心编号；

计算目标框架下的分析中心站坐标解；

S3、对于所有的分析中心站坐标产品，根据分析中心站坐标解建立站坐标综合模型，得到总的综合法方程系统，利用最小二乘原理得到综合解估计值和相似变换参数以及站坐标残差，实现多分析中心站坐标产品综合；

S4、利用站点间的相关性信息识别分析中心站坐标产品中的异常站点，并利用方差分量估计方法进行各个分析中心产品的权重分配；

所述利用站点间的相关性信息识别分析中心站坐标产品中的异常站点具体包括：

将每个分析中心的站坐标产品划分为一类观测数据，采用下式计算每个分析中心的后验方差：

；

其中，表示后验方差；nsit表示分析中心s的站坐标产品中的站点总数；设表示分析中心s中的跟踪站isit的坐标后验残差，表示分析中心s中的跟踪站k的坐标后验残差，表示站点isit和k之间的相关信息；表示综合法方程系统的信息矩阵逆；表示分析中心站坐标观测量的权阵；表示分析中心站坐标观测方程的设计矩阵；tr(·)为迹运算函数；

逐一计算去除每个站点后的方差改变量，并对进行从大到小排序处理，满足如下关系式时剔除站点观测量：

；

其中，ndel表示待剔除的异常站坐标数量；

同时，若下式成立：

；

则将第idel个站点标记为异常站点，在站坐标综合中剔除；

S5、重复以上步骤S3~S4，直至各个分析中心的权重趋于稳定或达到预设的最大迭代次数，得到权重分配结果；

S6、基于权重分配结果求解总的综合法方程系统，输出综合站坐标。

2.根据权利要求1所述的多分析中心精密站坐标产品综合方法，其特征在于，所述预处理具体包括：初步判定各精密站坐标产品的格式，并统计精密站坐标产品中站点数、核心站数以及存储原始法方程系统信息。