(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210998113.X (22)申请日 2022.08.19 (71)申请人 武汉大学 地址 430072 湖北省武汉市武昌区珞珈山 街道八一路2 99号 (72)发明人 李昭　鲁杨　姜卫平　叶世榕　 (74)专利代理 机构 武汉科皓知识产权代理事务 所(特殊普通 合伙) 42222 专利代理师 郑勤振 (51)Int.Cl. G06F 30/20(2020.01) G06F 17/11(2006.01) G06F 17/16(2006.01) G06F 119/14(2020.01) (54)发明名称 组合卫星重力场的新方法及系统 (57)摘要 公开了一种组合卫星重力场 的新方法及系 统， 通过利用方差分量估 量(VCE)、 熵权法(EWM)、 变异系数法(CVW)等数学模型对不同机构重力场 数据的组合， 组合模型相较于原始模 型具有更高 的精度， 并且能够将原始模型中较为明显的误差 进行抵消， 组合前后模型间变化趋势由原始模型 间差异最大的0.011mm/month降低到0.001mm/ month、 水柱高的最大差异由0.945cm减小为 0.211cm。 组合卫星重力场是一种可行并且具有 较高精度的卫星重力场数据应用新模式， 该模式 能够为卫星重力场应用研究提供更加稳定与高 精度产品， 而且也是以后未来卫星重力场应用研 究的重要方式。 权利要求书2页 说明书7页 附图4页 CN 115455659 A 2022.12.09 CN 115455659 A 1.一种组合卫星重力场的方法， 其特征在于， 利用方差分量估量、 熵权法或变异系数法 确定n种卫星重力场产品的权重矩阵， 将 每种卫星重力场产品乘以相应的权重， 然后 将赋完 权重的卫星重力场进行累加得到组合重力场产品。 2.根据权利要求1所述的组合卫星重力场的方法， 其特征在于， 利用方差分量估量确定 n种卫星重力场产品的权 重矩阵的方法为： k种精度的观测值相应的误差方程： V1＝A1X‑L1 V2＝A2X‑L2 … Vk＝AkX‑Lk 式中， i＝1， …， k； Ai为ni×k设计矩阵， ni为第i类观测值的个数； Li为ni×1观测向量； Vi 为对应于Li的残差向量； Li的先验权阵为P1， P2，…， Pk， 并且假设Li和Lj(i≠j)是独立的， 令 及法方程数据 方差分量估计值由下式确定 式中， ni为第i类观测值的个数， Ni为对应的法矩阵， N‑1为总法矩阵的逆矩阵， 第i类观 测值对应 的单位权方差因子； Nj为第j类观测值对应 的法矩阵(j≠i)， 第j类观测值对应 的单位权方差因子； 利用下式计算单位权方差因子 。 式中， tr表示提取对角线元 素； 根据单位权方差因子得到n种卫星重力场产品的权 重矩阵。 3.根据权利要求1所述的组合卫星重力场的方法， 其特征在于， 利用熵权法确定n种卫 星重力场产品的权 重矩阵的方法为： 先对数据进行无量纲化处理， 得到归一化数值， 再根据信 息论确定信 息熵， 计算各个指 标的信息熵， 通过信息熵计算各指标的熵值权 重， 公式为： 权　利　要　求　书 1/2 页 2 CN 115455659 A 2式中： K1为原始数据标准化得到的值； xij(i＝1， 2， …， n； j＝1， 2， …， m)为各指标的数 值； Pij(i＝1， 2， …， n； j＝1， 2， …， m)为归一化数值； ej(j＝1， 2， …， m)为信息熵； wj为j指标 的熵值权 重； m为指标个数， n 为指标的样本数。 4.根据权利要求1所述的组合卫星重力场的方法， 其特征在于， 变异系数法确定n种卫 星重力场产品的权 重矩阵的方法为： 设有m个参评样本、 n个指标， Xij(i＝1...m， j＝1....n)表示第i个样本第j个指标用的 数值， 计算出 各指标的均值 和标准差Sj： 则各指标的变异系数为： 式中， Vj表示第j个指标的变异系数 对Vi作归一化处理， 得到各指标的权值 为： 式中， Wj表示第j个指标的权 重。 5.一种组合 卫星重力场的系统， 其特 征在于， 包括： 处理器； 用于存储所述处 理器可执行指令的存 储器； 其中， 所述处 理器被配置为执 行所述指令时实现权利要求1 ‑4任一项所述的方法。 6.一种非临时性计算机可读存储介质， 其上存储有计算机程序， 其特征在于， 所述计算 机程序被处 理器执行时实现如权利要求1至4任一项所述的方法。权　利　要　求　书 2/2 页 3 CN 115455659 A 3