任意精度的子午线弧长递归计算.pdf

第２７卷第５期 大地 测 量 与地 球 动 力 学 Ｖｏｌ．２７Ｎｏ．５ ２００７年１０月 ＪＯＵＲＮＡＬＯＦＧＥＯＤＥＳＹＡＮＤＧＥＯＤＹＮＡＭＩＣＳ 　 Ｏｃｔ．，２００７ 　　文章编号：１６７１５９４２（２００７）０５００５９０４ 任意精度的子午线弧长递归计算 刘仁钊　伍吉仓 （同济大学测量与国土信息工程系，上海　２０００９２） 摘　要　根据子午线弧长计算公式，推导出基于递归关系的任意精度的子午线弧长的一种新的计算方法，以及基 于Ｃ／Ｃ＋＋语言的子午线弧长具体方法的实现，最后在ＶＣ６．０中进行了验证。结果表明，该方法可满足不同精度 要求的弧长计算。 关键词　子午线弧长　公式推导　递归计算　Ｃ／Ｃ＋＋程序　精度分析 中图分类号：Ｐ２２６　　　　文献标识码：Ａ ＲＥＣＵＲＳＩＶＥＣＯＭＰＵＴＡＴＩＯＮＯＦＭＥＲＩＤＩＡＮＡＲＣ ＬＥＮＧＴＨＷＩＴＨＤＩＳＣＲＥＴＩＯＮＡＲＹＰＲＥＣＩＳＩＯＮ ＬｉｕＲｅｎｚｈａｏａｎｄＷｕＪｉｃａｎｇ （ＤｅｐａｒｔｍｅｎｔｏｆＳｕｒｖｅｙｉｎｇａｎｄＧｅｏｉｎｆｏｒｍａｔｉｃｓ，ＴｏｎｇｊｉＵｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ　２０００９２） Ａｂｓｔｒａｃｔ　Ａｎｅｗｍｅｔｈｏｄｆｏｒｒｅｃｕｒｓｉｖｅｃｏｍｐｕｔａｔｉｏｎｏｆｍｅｒｉｄｉａｎａｒｃｌｅｎｇｔｈｗｉｔｈｄｉｓｃｒｅｔｉｏｎａｒｙｐｒｅｃｉｓｉｏｎｉｓｄｅ ｄｕｃｅｄｏｎｔｈｅｂａｓｉｓｏｆｔｈｅａｒｃｆｏｒｍｕｌａ，ａｎｄｉｔｓｄｅｔａｉｌｅｄａｌｇｏｒｉｔｈｍａｓｗｅｌｌａｓｐｒｏｇｒａｍｉｎＣ／Ｃ＋＋ａｒｅｇｉｖｅｎ．Ｔｈｉｓ ｍｅｔｈｏｄｆｏｒａｒｃｌｅｎｇｔｈｃｏｍｐｕｔａｔｉｏｎｉｓｐｒｏｖｅｄｔｈｒｏｕｇｈＭｉｃｒｏｓｏｆｔＶＣ＋＋６．０，ｔｈｅｒｅｓｕｌｔｉｎｄｉｃａｔｅｓｔｈａｔｄｉｆｆｅｒｅｎｔｐｒｅｃｉ ｓｉｏｎｓｒｅｑｕｉｒｅｍｅｎｔｉｎｖａｒｉｏｕｓｃｏｎｄｉｔｉｏｎｓｃａｎｂｅｍｅｔｗｉｔｈ． Ｋｅｙｗｏｒｄｓ：ｍｅｒｉｄｉａｎａｒｃｌｅｎｇｔｈ，ｆｏｒｍｕｌａｄｅｄｕｃｔｉｏｎ，ｒｅｃｕｒｓｉｖｅｃｏｍｐｕｔａｔｉｏｎ，Ｃ／Ｃ＋＋ｐｒｏｇｒａｍ，ｐｒｅｃｉｓｉｏｎａｎａｌｙｓｉｓ 法，研究了子午弧长的正反解，文献［８］在现有弧长 １　引言 公式上作了有益的探索。本文根据文献［９］的研究 众所周知，子午线弧长是大地测量学中的一个 思路推导出通过递归方法计算子午线弧长的算法， 基本量，在大地坐标计算、ＧＰＳ定位和导航、ＧＩＳ地 该算法过程简洁，逻辑性强，无需针对具体的椭球， 理坐标数据表达、高斯平面直角坐标计算以及经典 适用于任意弧长精度。 的弧度测量中，都需要计算两点之间的子午线弧长。 由于子午线弧长的计算涉及到椭圆积分，所以在现 ２　子午线弧长的计算 ［１，２］ 有的教材 中都是将其展开成级数形式，再用逐 文献［２］将椭圆积分化成正弦倍角函数，根据 项积分的方法求出满足一定精度要求的计算公式。 －１０ 精度系数项展开至ｅ ，计算从赤道纬度０°～纬度 为了使公式在应用中表达更简洁或精度更高，不少 φ的一段子午线弧长的公式如下： 文献对此进行过探讨［３～５］。文献［６］甚至将级数展 φ ３ ２ ２ ２ －