[发明专利]iGMAS多分析中心多卫星系统精密轨道产品综合方法

权利要求书

1.iGMAS多分析中心多卫星系统精密轨道产品综合方法，其特征是，包括：

S100多个分析中心获取各卫星的轨道数据，根据轨道数据对各分析中心的精密轨道产品进行预处理；

S200对各卫星，取该卫星所对应的所有分析中心的轨道数据的中位数，作为该卫星初始的综合轨道；

S300根据历史轨道产品精度将各分析中心提供的不同卫星的轨道数据分成p类观测数据，分别建立各分析中心提供的各类观测数据和综合轨道间的相似变换方程，并根据相似变换方程构建对应的误差方程；

S400采用最小二乘法解算误差方程，分别获得各分析中心的各不同类观测数据相对于综合轨道的变换参数、综合轨道、以及各卫星的轨道观测值残差；

S500采用等价权函数获得各卫星的权重，基于各分析中心的各类观测数据中各卫星的权重，采用方差分量估计法估计各分析中心的各类观测数据的方差分量，根据方差分量计算并更新各分析中心的各类观测数据的权重；

S600根据各类观测数据的权重更新综合轨道，以更新后的综合轨道作为后一次迭代的输入；

S700重复步骤S300～S600进行迭代，直至收敛，最后一次迭代更新后的综合轨道即最终的综合轨道；

步骤S500具体为：

S510采用等价权函数获得各卫星的权重P_l^i,j：

其中，P_l^i,j表示第i个分析中心的第j类观测数据中卫星l的权重；表示第i个分析中心的第j类观测数据中卫星l的中误差；表示第i个分析中心的第j类观测数据的所有卫星的抗差中误差；c₀和c₁表示淘汰域的临界值，为经验数值；

S520基于各分析中心的各类观测数据中各卫星的权重P_l^i,j，采用方差分量估计法计算各分析中心的各类观测数据的方差分量，根据方差分量更新各类观测数据的权重；

所述方差分量估计方程为其中，H表示方差分量参数的系数矩阵，DIAG{}表示对角阵函数，表示克罗内克积，I_n表示由n个数值1组成的列向量，n为分析中心数，m^i,j表示分析中心i的第j类观测数据中的卫星数，q^i,j表示分析中心i的第j类观测数据的卫星历元数；表示由各分析中心各类观测数据的方差分量组成的向量；表示所有分析中心第j类观测数据的方差分量序列的中位数；W也表示方差分量参数的系数矩阵，V^i,j表示分析中心i的第j类观测数据的残差向量；表示分析中心i的第j类观测数据的先验权重，初始值设为1，后续迭代计算中采用进行更新；表示根据验后方差分量估计的分析中心i的第j类观测数据的权重；i＝1,2...n，j＝1,2...p。

2.如权利要求1所述的iGMAS多分析中心多卫星系统精密轨道产品综合方法，其特征是：

S100中所述的根据轨道数据对各分析中心的精密轨道产品进行预处理，具体为：

统计精密轨道产品中卫星系统数、卫星数、轨道数据时长、轨道数据间隔以及各精密轨道产品所在的时间系统；

对需要进行产品一致性改正的各分析中心的精密轨道产品进行产品一致性改正。

3.如权利要求1所述的iGMAS多分析中心多卫星系统精密轨道产品综合方法，其特征是：

步骤S300中，所构建单历元的误差方程如下：

该误差方程为第i个分析中心的第j类观测数据的误差方程，其中，表示第i个分析中心的第j类观测数据；表示第i个分析中心的第j类观测数据对应的观测值误差；表示第j类观测数据的综合轨道，记为第j类综合轨道；A^i,j表示第i个分析中心的第j类观测数据的观测矩阵；表示第i个分析中心的第j类观测数据对应的第b类变换参数；E^i,j为大小为3×3的单位矩阵，表示第i个分析中心的第j类观测数据相对应的第b类变换参数；表示第i个分析中心的第j类观测数据对应的第b类旋转参数的系数矩阵；表示第i个分析中心的第j类观测数据对应的第b类尺度参数的系数矩阵；表示第i个分析中心的第j类观测数据对应的第b类平移参数；表示第i个分析中心的第j类观测数据对应的第b类旋转参数；表示第i个分析中心的第j类观测数据对应的第b类尺度参数。

4.如权利要求1所述的iGMAS多分析中心多卫星系统精密轨道产品综合方法，其特征是：

步骤S400进一步包括子步骤：

S410对变换参数附加相应的约束方程其中，表示第i个分析中心的第j类观测数据相对于综合轨道的变换参数，表示变换参数的权重，其初始值为1；

S420采用步骤S300方法建立所有分析中心所有观测历元的第j类观测数据的误差方程，采用最小二乘原则分别解算各分析中心的第j类观测数据相对综合轨道的变换参数、相应的综合轨道、以及观测值残差，其中，j依次取1、2、……p；l表示卫星编号，l＝1,2，……m；P_l^i,j表示第i个分析中心的第j类观测数据中卫星l的权重，初始的权重设为1；表示第i个分析中心第j类观测数据中卫星l的观测值误差。