测绘学报 ›› 2018, Vol. 47 ›› Issue (11): 1437-1445.doi: 10.11947/j.AGCS.2018.20170659

• 大地测量学与导航 • 上一篇    下一篇

超高阶重力场模型最小二乘快速实现

田家磊1,2, 李新星1,3, 吴晓平1, 邢志斌1   

  1. 1. 信息工程大学地理空间信息学院, 河南 郑州 450001;
    2. 63850部队, 吉林 白城 137000;
    3. 武汉大学测绘学院, 湖北 武汉 430079
  • 收稿日期:2017-11-18 修回日期:2018-05-14 出版日期:2018-11-20 发布日期:2018-11-29
  • 通讯作者: 李新星 E-mail:minibad@126.com
  • 作者简介:田家磊(1988-),男,博士生,研究方向为物理大地测量。E-mail:tianjialei3037@163.com
  • 基金资助:
    国家自然科学基金(41274029;41404020;41774018;41504018;41674082);信息工程大学自立课题(2017503902;2016601002)

Fast Realization of Ultra-high-degree Geopotential Model by Improved Least-squares Method

TIAN Jialei1,2, LI Xinxing1,3, WU Xiaoping1, XING Zhibin1   

  1. 1. Institute of Geospatial Information, Information Engineering University, Zhengzhou 450001, China;
    2. 63850 troops, Baicheng 137000, China;
    3. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China
  • Received:2017-11-18 Revised:2018-05-14 Online:2018-11-20 Published:2018-11-29
  • Supported by:
    The National Natural Science Foundation of China (Nos. 41274029;41404020;41774018;41504018;41674082);The Self-reliance Project of Information Engineering University (Nos. 2017503902;2016601002)

摘要: 采用最小二乘方法解算超高阶重力场模型,不可避免会遇到大型矩阵的计算,直接求解是难以实现的。本文从重力场模型的基本观测方程出发,利用正余弦函数和面球谐因子的正交性,分析系数矩阵及法矩阵的特点,在法矩阵块对角化的基础上,利用系数矩阵求解法矩阵时“次m”递增的特点,对法矩阵求解方程进行约化、对Legendre函数的计算和存储方式进行了设计,结合缔合Legendre函数关于赤道的对称性,解决了大型矩阵存储及计算效率低下的难题,实现了超高阶重力场模型最小二乘方法的小存储、高效率的解算。通过试验模拟,改进后的方法相比传统块对角方法效率提高300倍,利用此方法可以在普通PC机上快速、高精度地解算2160阶次超高阶重力场模型,算法精度相比数值积分方法至少提高了5个数量级,并且在一定程度上可以评估原始观测数据的精度。

关键词: 块对角矩阵, 最小二乘, 超高阶地球重力场模型, Legendre函数

Abstract: Using the least-squares method to solve the ultra-high-degree geopotential model is very difficult because this method needs to calculate massive data.In this paper,we first analyzed the characteristics of coefficient matrix and normal matrix based on the basic observation equation of geopotential model using the orthogonality of trigonometric function and surface coherence factor.Then,we reduced the matrix solution equation and designed the calculation and storage method of Legendre function based on normal matrix block diagonalization using the increment characteristic of "secondary M" during the process of solving normal matrix by coefficient matrix solution method.At last,the problems of storage and low computational efficiency of large matrix were solved by combining the equator symmetry characteristic of Legendre function,reduction of the matrix solution equation and design of the calculation and storage method of Legendre function.Through the experimental test,the improved method compared to traditional block diagonal method canimprove the efficiency by 300 times,by using this method we can solve ultra-high-degree geopotential model in ordinary PC,and the precision of the model has been improved by 5 orders of magnitude at least when compared with the numerical integral method.At the same time,the accuracy of the original data can be evaluated.

Key words: block-diagonal, least-squares, ultra-high-degree geopotential model, Legendre function

中图分类号: