基于向量化的扰动引力矢量快速并行算法  

Fast Parallel Algorithm of Disturbing Gravity Vector Based on Vectorization

作　　者：黄炎 王庆宾[1] 李国强[2] 冯进凯 谭勖立 Huang Yan;Wang Qingbin;Li Guoqiang;Feng Jinkai;Tan Xuli(Institute of Geospatial Information,Information Engineering University,Zhengzhou 450001,China;The Command Department,Information Engineering University,Zhengzhou 450001,China)

机构地区：[1]信息工程大学地理空间信息学院,郑州450001 [2]信息工程大学指挥系,郑州450001

出　　处：《吉林大学学报（地球科学版）》2021年第6期1863-1871,共9页Journal of Jilin University：Earth Science Edition

基　　金：国家自然科学基金项目(41574020);信息工程大学校自立课题(2105070232)。

摘　　要：无人飞行器在飞行过程中受地球外部扰动引力影响,对其精确控制需要计算飞行轨迹点扰动引力。为有效恢复地球外部引力场,常采用地球重力场位系数模型进行计算,但其计算耗时随模型阶数的升高而呈指数增加。本文提出一种扰动引力矢量计算的向量化方法,推导出利用地球重力场模型计算扰动引力三分量的向量化表达公式,并利用CUDA(compute unified device architecture)异构并行算法进行并行化以达到快速计算单点扰动引力的目的。通过仿真实验可知,本文所提方法可有效减少单点扰动引力计算耗时,利用跨阶次递推法进行单点计算加速比可达8以上,最高可达13.20;利用Belikov递推法进行单点计算加速比可达6以上,最高可达8.99。Unmanned aerial vehicle(UAV) is affected by the disturbing gravity of the Earth during its flight. To control the UAV accurately, it is necessary to calculate the disturbing gravity of the flight trajectory points. In order to effectively restore the external gravitational field of the Earth, the Earth gravitational model is often used, but its calculation time increases exponentially with the increase of the order of the model. In this paper, an improved vectorization method for the calculation of perturbation gravity is proposed, and a CUDA(compute unified device architecture)heterogeneous parallel algorithm is used for parallelization to achieve the purpose of rapid calculation of single point perturbation gravity. The simulation results show that the proposed method can effectively reduce the time-consuming of single-point calculation of perturbation gravity: The acceleration ratio of single-point calculation can be more than eight times and the highest can be up to 13.20 times by using the step-by-step recursion method;By using Belikov recursion method, the acceleration ratio of single point calculation can be more than six times and the maximum can be 8.99 times.

关 键 词：向量化 并行计算 扰动引力 地球重力场位系数模型 加速比 

分 类 号：P221[天文地球—大地测量学与测量工程]