外空扰动重力计算中的截断误差

给出了由地面重力数据计算外空扰动重力矢量的公式，并根据Wong and Gore截断理论给出了外区重力异常对计算点影响的公式。根据外空扰动重力与地面数据分辨率及其覆盖范围之间的关系，将重力异常分成不同频段、不同分辨率，分别计算了截断误差。利用全球重力位模型，计算出不同频段的截断误差，并给出了各频段相应的积分半径。对于模拟和检验航空矢量重力数据有一定的参考价值。     

扰动引力快速确定的替代算法

扰动引力的点质量表达方法使得扰动引力的计算模式具有最为简单的结构 ,但它本身也具有计算量大的缺点 ,与实时计算的目标相矛盾。为解决这一矛盾 ,文中用三次等距B样条函数进行试验。结果表明 ,在已知预设弹道坐标的条件下 ,三次等距B样条函数方法是适于进行快速精确逼近弹道扰动引力的数值算法     

扰动引力快速确定的替代算法

扰动引力的点质量表达方法使得扰动引力的计算模式具有最为简单的结构,但它本身也具有计算量大的缺点,与实时计算的目标相矛盾.为解决这一矛盾,文中用三次等距B样条函数进行试验.结果表明,在已知预设弹道坐标的条件下,三次等距B样条函数方法是适于进行快速精确逼近弹道扰动引力的数值算法.     

扰动重力水平分量对惯导系统的位置误差影响

研究了不同运动状态下扰动重力水平分量（HDG）对高精度惯导系统（inertial navigation system,INS）的位置误差影响。首先推导了HDG对INS误差影响的状态空间方程,进而推导出3种运动条件下INS位置误差与HDG之间的解析关系式,设计了基于惯导解算求解上述影响的方法。在匀速运动条件下,分别通过解析式与惯导解算两种方法计算了相同HDG引起的INS位置误差。解析式计算结果表明,±80 mGal（1 mGal=10-5 m/s2）范围内变化的HDG约可引起最大约3 000 m的INS位置误差;对两种方法计算结果的比较显示,所得INS位置误差的量级与变化情况基本一致,两组结果验证了各自方法的有效性。     

重力异常阶方差模型的构建及在扰动场元频谱特征计算中的应用

利用最新的全球引力位模型-EGM2008对经典的重力异常阶方差模型进行了分析比较,分析表明,经典的阶方差模型由于限于当时的观测条件,已经不能准确地描述扰动场元在各个频段的频谱分布。在Moritz阶方差模型基础上,利用EGM2008位模型获得的2160阶阶方差重新构建了新的分段重力异常阶方差模型-TSD模型,该模型与EGM2008位模型计算的阶方差比较其标准差和均值分别为0.25mgal2 、0.0 。利用TSD模型计算了不同频段内大地水准面高、重力异常、扰动重力、垂线偏差四个重力场扰动场元的频谱特征,计算结果表明：扰动场元频谱分布较之传统分析结果有较大的变化,其中重力异常、扰动重力及垂线偏差在中、低频部分的能量有明显的增加而高频及甚高频部分的比重有明显的减少。     

全球及局部海洋扰动重力反演的快速解析方法

从经典边值问题理论及球谐函数理论出发,在空域推导获得了由大地水准面高以及垂线偏差计算扰动重力的解析计算公式,为利用卫星测高数据反演海洋扰动重力提供了理论基础。针对全球海洋区域和局部海洋区域的扰动重力反演,在前人已有工作基础上,提出了改进的基于一维FFT的精确快速算法,保证了计算结果与原解析方法完全一致,且计算速度提高约20倍。该算法在提高计算效率的同时避免了由于引入FFT而产生的混叠、边缘效应问题,而且对观测数据的序列长度没有硬性要求,使得应用更加灵活。利用EGM2008地球重力场模型分别生成了2.5'分辨率大地水准面高数据和垂线偏差数据,按照本文提出的改进方法(采用全球积分计算)分别反演获得了全球及局部海洋区域的扰动重力。经比较分析,由大地水准面和垂线偏差分别反演获得的扰动重力其差异在0.8×10^-5 m/s²以内,这说明两种反演方法是基本一致的,但在数据包含系统误差的情况下,由垂线偏差反演扰动重力具有一定优势。     

扰动重力梯度的非奇异表示

在局部指北坐标系中用地心球坐标来表示扰动重力梯度张量,当计算点趋近于两极时,由于Legendre函数的一阶和二阶导数以及分母上所含余纬的正弦函数,将导致扰动重力梯度张量的计算出现无穷大。因此,本文引入了Legendre函数的一阶和二阶导数以及 无奇异性的计算公式,并且进一步推导了 无奇异性的计算公式。在将Legendre函数的一阶和二阶导数以及 、 无奇异性的计算公式代入到扰动重力梯度张量各分量的求解中时,又充分考虑了m等于0,1,2以及其它量时的复杂情况,建立了扰动重力梯度张量各分量无奇异性的详细计算模型。通过模拟实验表明,本文所建立的详细计算模型不仅能够完全满足当前卫星重力梯度张量计算的精度要求,而且模型稳定、可靠、易于编程实现。     

扰动引力多核快速并行算法

远程飞行器和弹道导弹的精确控制与导航会受到地球外部扰动引力场的影响,实践中常利用地球重力场位系数模型计算扰动引力矢量,但其计算效率随着所用模型阶次的升高而显著降低。针对这一问题,文中提出利用OpenMP并行算法快速计算轨迹点扰动引力矢量,在不损失精度的前提下有效提高了计算效率。通过数值实验可知,所提方法可以提高扰动引力矢量的计算效率,加速比达到6倍。这也为重力场元快速计算提供了参考方案。     

扰动重力场精密确定与逼近效果分析

介绍了局部重力场逼近的点质量方法;详细论述了组合点质量模型的构制方法;利用实测重力异常数据建立了某区域组合点质量模型;分析了各层点质量的散射误差、逼近效果、随高度的变化特征。结果表明,组合点质量模型能以较高精度逼近外部扰动重力场。     

航空和地面重力测量数据的一体化处理方法

在三维经典扰动位协方差模型基础上,推导获得局部重力异常协方差模型、局部扰动重力协方差模型,以及扰动重力和重力异常的互协方差模型,这3个模型在忽略小项后是一致的,即航空扰动重力数据与地面重力异常数据的联合处理采用一个模型即可,大大简化了处理流程.以最小二乘配置理论及本文推导的协方差模型为基础,构建了航空扰动重力测线数据与地面重力异常数据一体化处理的具体方法,利用该方法可推估地面待求点的重力异常数据,且处理过程自动完成了航空扰动重力数据的向下延拓,同时也大大削弱了航空数据中存在的系统偏差影响.在某试验区开展了一体化处理试验,结果表明,采用一体化处理得到的重力异常数据比单独使用5'分辨率航空重力数据推估获得的数据精度提高2.4 mGal,比单独使用2.5'分辨率地面重力数据推估获得的数据精度提高1.2 mGal.