一种适用任意距离的大地主题解算方法的研究和试验.pdfVIP

第31卷第2期 测绘科学 voI．31No．2 2006年 3月 ScienceofSurveyingandMapping Mar．2006 一 种适用任意距离的大地主题解算方法的研究和试验 何 菊 ，胡 鹏 ，胡 海 (武汉大学资源与环境学院，武汉 430079) 【摘 要】介绍了嵌套系数法解算大地主题问题的原理和公式；并以详尽的算例对该方法进行 了实验和分析，与 贝塞尔方法以及高斯平均约数算法相 比较，指出该算法是一种适用于任意大地距离的高精度算法，并且不存在奇 异 问题。 【关键词】大地问题正反算；大地距离；方位角 【中图分类号】P22；P226 【文献标识码】A 【文章编号】1009—2307(2006)02—0029—03 1 引 言 我们通常遇到计算大地椭球上距离和位置的问题 ，即 大地测量主题的正反算问题。百余年来，许多数学家和测 量学家研究出几十种方法来解决这类问题。其 中大部分适 用于短距离 (120krn)，一部分适用于中距离 (400krn)，只有 几种适用于长达 20000km 的距离。比较典型的有勒让尔级 数式、高斯平均引数公式 、贝塞尔解算方法 、陈俊勇长公 式等。 图1 椭球面 图2 球面 目前很多大地主题解算方法虽然比较完善，但经过许 多人的实践发现这些正反算法 (特别是反算)仍然存在一些 这就是贝塞尔两个非常著名的边长公式和经差投影公式。 明显的缺陷 J，主要表现在： 贝塞尔方法是在这两个公式的基础上引入系数，逐步展开 1)应用范围受到一定的限制。尽管 目前很多算法在理 积分 ，再查表得到其正反算解。 论上适用于任意大地距离 ，但 实际上受数值精度的影响， 2．2 由贝塞尔方法改进的嵌套系数法 它们几乎都只适用于一定范围而不能解算任意大地距离的 嵌套系数法_3J还是以这两个基本公式为出发点，采用不 主题 。 同的积分路线和各级数展开式中函数幂及其系数本身所具有 2)奇异问题。当两大地点位于同一子午圈或同在赤道 的叠乘性质，通过嵌套约化，推导出计算椭球改正项的三个 时，此算法不能计算，必须采用别的方法解算。 嵌套系数Kl，K2，K3和两个叠乘性质改正数 、△ 。 这些问题的存在，使得大地主题计算问题变得很复杂。 2．2．1 正算公式和步骤 大地上的坐标点有远有近，如果所使用的解算公式不能计 1)根据公式(1)式有ds：(1+e,2s “)¨2da将其展开， 算任意距离的正反算，则对于某些计算结果会产生很大误 设 t：1／4er2sinUn，得到：(1+esin2“)／2= (1+t-(3／4)t0+ 差，或者选择不同的算法公式而使问题 (特别是对于大批量 (5／4)t一(175／64)t)一(t—t+(15／8)t一(35／8)t4) d 的计算)变得复杂化。下面介绍一种嵌套系数法，此算法的 一 ((1／4)t 一 (3／4)t + (35／16)t)耐 d一 ((1／8)t 一 基础是贝塞尔的两个微分公式，但算法不 同，我们也可以 (5／1B)t)a~6a一(5／6)t4(z6：8o。 把它看作是贝塞尔的一种改进算法。 2)设 K1：1+t{1一t／4[3一t(5—11t)]}，K2：t{1