大地主题常微分方程组解算的数值方法 ——以MathCAD为工具

利用数学软件MathCAD的常微分方程求解函数Rkadapt()直接解算大地主题正反算常微分方程组,可一次性求解n个节点的大地元素值;在此基础上建立线性插值函数,可计算任意点的大地元素值。方法简洁、通用,求解精度高,适用于长短距离的大地主题解算。大地主题反算的起点方位角A1由Bessel函数方法求取。在大地主题解算的基础上,选择深圳地区最西与最东的二等GPS控制点Ⅱ3及Ⅱ54,估算了该区高斯投影6°带及3°带直角坐标系以及深圳独立坐标系的投影变形值。