利用最小二乘直接法反演卫星重力场模型的MPI并行算法

针对海量卫星重力数据反演高阶次地球重力场模型的密集型计算任务与高内存耗用问题,基于MPI实现了最小二乘直接法恢复高阶次位系数的并行算法。引入并行读写、分块存储与分块计算等方式完成了设计矩阵的构建、法方程的形成与求解等密集型计算任务的并行算法,数值计算结果表明三者的并行相对效率峰值可分别达到95%、68%、63%。利用GOCE轨道跟踪和径向扰动重力梯度数据(共518 400个历元)分别反演了120、240阶次地球重力场模型,计算时间仅为40 min、7 h,内存耗用峰值仅为290 MB、1.57 GB;采用与GOCE同等噪声水平的观测数据恢复的重力场模型精度与GOCE已发布模型的解算精度相一致,联合GRACE和GOCE的解算模型能够实现二者独立信息的频谱互补,表明本文方法可高效稳定地恢复高阶次地球重力场模型。     

����GOCE��GRACE������ģ��λϵ��������Ϸ�

???GOCE??GRACE????????????????е???????????о?????GOCE??GRACE?????????λ??????????????????????λ???????????????????????????Э????????????????????GRACE?????????ITG??GRACE2010s??GOCE?????????GO_CONS_GCF_2_DIR_R1?????????????????????????????????????????????????????????????????GRACE??GOCE????????????????????????????徫?????????GOCE??GRACE????????GOCO01S??????????????????Э?????????????????????????????????m??λ?????????????