常用超高阶次缔合勒让德函数计算方法对比分析

超高阶重力场模型的构制与应用的数值稳定性取决于超高阶次缔合勒让德函数的计算方法。介绍了4种常用的缔合勒让德函数的递推计算方法:标准前向列推法、标准前向行推法、Belikov列推法和跨阶次递推法,并对计算程序做了改进,从计算速度、计算精度两方面分析比较了阶次高至3 060阶的各种方法的优劣。数值试验证明,跨阶次递推法是计算超高阶次缔合勒让德函数较优的方法。

超高阶重力场模型计算高程异常算法与实现

详细阐述了利用地球重力场模型位系数计算高程异常的数学模型。选取标准向前列递推法和跨阶次递推法作为缔合勒让德函数的计算方法,分析了两种算法的适用范围。以高精度的超高阶EIGEN-6C4重力场模型为例,选择全球范围内5°×5°的格网组成的2701个点和具有代表性的12个点计算高程异常,将计算结果和ICGEM网站计算结果比较验证了算法的可靠性,并绘制了全球范围内5m等高距的高程异常等值线图,实现了全球高程异常的可视化。

不同Legendre函数递推公式对计算球谐函数定积分的影响

针对球谐函数定积分计算中Legendre函数递推问题展开研究,分析了标准向前列推法、Belikov法、跨阶次法、X数法以及顾及麦克劳林级数展开式对球谐函数定积分计算的影响。利用Eigen6c-4地球重力场模型计算扰动引力梯度径向分量,分析不同方法之间的差异。实验表明,不考虑麦克劳林级数展开式时4种方法的相对精度在高纬度地区较差,但计算模型扰动引力径向分量的精度一致,结合麦克劳林级数式可提高高纬度地区定积分计算的相对精度,但会降低中低纬度地区定积分计算的精度,并且对高纬度地区扰动引力径向分量的影响极小,但会严重降低低纬度地区扰动引力梯度计算的精度。