Matlab gravitysphericalharmonic函数示例坐标参数取值疑问

核心前提说明 #

gravitysphericalharmonic 函数默认采用地心地球固连（ECEF）坐标系：原点在地球质心，x轴指向本初子午线与赤道的交点，z轴指向地球北极，符合右手坐标系规则。

问题1：南极上空25000米位置z坐标的取值逻辑 #

• EGM96重力模型匹配的参考椭球极半径（地心到南北极点的直线距离）为6356.751e3米，和你提到的赤道半径6378.137e3米的差值，是地球椭球扁率带来的固有差异。
• 南极处于ECEF坐标系的z轴负方向，上空25000米位置的地心距离为极半径加高度：6356751 + 25000 = 6381751 米，因此z坐标取为-6381.751e3。

问题2：赤道表面位置x坐标为负的原因 #

• 球谐重力模型的重力加速度始终指向地心，坐标的正负仅会改变加速度分量的符号，不会影响模值的正确性：如果位置取x轴正方向[6378.137e3, 0, 0]，输出的gx分量为负；取x轴负方向，输出的gx分量为正。
• 官方示例选择负x坐标只是为了让输出的重力分量数值为正，更符合常规使用习惯，没有特殊物理含义。

内容的提问来源于stack exchange，提问作者Mei Win