1、高程异常高差改正数高程是重力位(数)在地固坐标参考系中的几何表达两点之间的高差是其重力位差在地球空间的几何表达。正(常)高在地固坐标参考系中客观唯一且精密可测都是满足唯一性和可测性要求的大地测量要素。地面点的大地高(椭球)H等于地面点正高h*与大地水准面(椭球)高N之和也等于地面点正常高h与该点高程异常 ζ之和式(11.7)是GNSS代替水准测量由N(或 ζ)测定地面点正高h*(或正常高h)的大地测量学依据。由式(11.7)得地面点的正高与正常高之差为可见任意地面点的正高h*与正常高h之差等于该点高程异常 ζ与大地水准面高N之差Δζ。Stokes 边值问题积分解是大地水准面及其外部整个地球空间的扰动位(Hofmann2006)可同时确定大地水准面高N及外部高程异常 ζ(广义Stokes公式)。特别地Stoke边值问题积分解约束了 ζ与N之间的解析函数关系(章传银2017)[功能]利用参考地球重力场位系数模型由计算点与目标点大地高计算地面及近地空间点高程异常径向梯度(cm/km)、高程异常高差改正数(m)的模型值。[输入]近地空间计算点文件参考重力场位系数模型文件。近地空间计算点文件记录格式:点号/点名经度(度小数)纬度(度小数)大地高(m)...目标大地高...位系数模型文件第1行约定为位系数尺度参数:地心引力常数GM与椭球长半轴a。尺度参数(GM、a)是位系数模型的表示基准。[参数设置]设置计算点文件格式输入计算点大地高和目标大地高在记录中的列序号输入位系数模型最大计算阶数。[输出]计算点模型高程异常高差改正数结果文件。在源点值文件记录的基础上增加一列该点的高程异常径向梯度(cm/km)和一列高程异常高差改正数(m)的模型值保留4位有效数字。参考重力场模型及计算阶数应保持不变。最大计算阶数一般不大于720。当计算点在地面、目标点位于大地水准面时程序计算正高与正常高之差即高程异常与大地水准面高之差。高程异常高差改正数m模型高程异常高差改正数或模型残差高程异常高差改正数或模型残差高差改正数实测重力校正值。高程异常径向梯度cn/km模型径向梯度或模型残差径向梯度或模型残差径向梯度实测重力校正值。2、程序编制程序首先参考了文献3和文献4的代码。但是用参考文献中的代码计算结果与上面有差异。因为参考文献使用的11.9式的近似公式实际上严格来说高程异常的径向梯度公式为用严格公式的计算结果如下图。计算结果差别在0.1mm以内。再看文件比较结果英文标题的文件是章老师程序计算结果中文标题的文件是本程序计算结果。计算结果差别均在0.1mm以内。参考文献过指定点模型重力等位面构造 - 章传银利用重力场模型计算等正(常)高面的重力位差和大地高差 - 章传银利用重力场模型计算正高与正常高差异 - 章传银利用重力场模型计算高程异常径向变化率 - 章传银Molodensky边值校正Fortran代码 - 章传银正高与正常高的定义与差别-CSDN博客PAGravf4.5参考说明书.pdf