子午线收敛角与长度比简捷算法.pdfVIP

第 1期 矿 山 测 量 NO．1 2010年 2月 MINE SURVEYING Feb．2010 doi：10．3969／j．issn．1001—358X．2010．01．023 子午线收敛角与长度比简捷算法 周复旦 ，赵长胜 ，倪 飞，王德高 (徐 州师范大学测绘 学院 江苏 徐州 22l116) 摘要：子午线收敛角和长度比的计算公式是大地测量中的常规计算，但是公式推 导复杂，初学者难以 掌握 。提 出适合于计算机计算的子午线收敛角与长度比的简捷算法，并编写了相应的数据处理程 序 ，其效果与用传统公式解算结果相同。本算法公式推导简单 ，便于理解与编程 。 ． 关键词：子午线收敛角；长度比；高斯投影 中图分类号 ：P22 文献标识码 ：B 文章编号：1001—358X(2010)01—0076—02 大地方位角与高斯平面方位角互相转换时必须 推导子午线 收敛角 的传统方法是：首先将 式 计算子午线收敛角和方向改正 ，大地线长度和高斯 (2)、式 (3)代人到式 (1)中，再将分母用泰勒级数展 平面上的直线长度互相转换时必须计算长度 比。子 开与分子相乘 ，简化后将反正切 函数按级数展开，最 午线收敛角和长度比的计算公式是大地测量 中的常 后获得子午线收敛角计算式。 规计算式 。但是这些公式推导复杂 ，教师难讲，学生 用大地坐标计算子午线收敛角的公式为： 难以掌握 。由于现代计算机广泛应用 ，没有必要像 y = sinB · + 1sinBcos2BI(1+3叼 +277)+ 手算时代那样推导公式。本文提出适合于计算机计 算的子午线收敛角和长度比的简捷算法，编写了相 sinBcosBI(2一t)+… (4) 应的数据处理程序 ，与用传统公式相 比，计算结果相 用高斯平面坐标计算子午线收敛角的公式为： 同，但算法简单 ，便于理解与编程 。 y 一 。(1+ 一 一2 )+ 。 1 -f-午 线 收 取 用 阳 ．问稳 舁 法 (2+5￡；+3t；) (5) 子午线收敛角 的理论值…为： 式中符号的意义请见文献 。 Ox／Oy (1) = tan ～ _ ／ 以上公式的推导繁琐 ，初学者不容易学习。在 计算机发达 的今天，没有必要将公式推导搞得如此 式中， 和 是高斯投影面坐标 、对经差 z的 复杂。建议采用下述算法 ： 偏导数 ，为： ①首先将式 (2)、式 (3)设计为两个函数 dx(o，