子午线弧长的计算方法及精度分析.pdfVIP

１０６ 测　 绘　 通　 报　 　 　 　 　 　 　 　 　 　 　 　 ２０１７年　 第８期 引文格式:刘学杰ꎬ杨丽坤.子午线弧长的计算方法及精度分析[Ｊ].测绘通报ꎬ２０１７(８):１０６￣１０９.ＤＯＩ:１０.１３４７４/ ｊ.ｃｎｋｉ.１１￣２２４６.２０１７.０２６４. 子午线弧长的计算方法及精度分析 １ ２ 刘学杰 ꎬ杨丽坤 (１. 河南省中纬测绘规划信息工程有限公司ꎬ河南 焦作４５４０００ꎻ２. 郑州工业贸易学校ꎬ河南郑州４５０００７) 摘要:计算子午线弧长除了采用经典的级数展开算法之外ꎬ还可通过数值积分与常微分方程数值解法进行求解ꎮ 为评价各种算法 的精度ꎬ本文选取大地纬度自０°—９０°、间隔距离为１°、１′、１″的３组样本数据ꎬ分别基于传统算法、数值积分算法和常微分方程数值 算法３大类 １１种算法计算得到各组样本所对应的子午线弧长结果ꎬ并从算法精度和运算速度两个方面对各种数值算法进行了分 析与评价ꎮ 实例表明三阶、四阶Ｒｕｎｇｅ￣Ｋｕｔｔａ算法不仅精度高ꎬ而且运算效率是其他算法的２倍多ꎬ研究结果为计算子午线弧长的 提供了有效的算法模型ꎮ 关键词:子午线弧长ꎻ数值积分ꎻ常微分方程ꎻ展开算法 中图分类号:Ｐ２２６　 　 　 　 　 　 文献标识码:Ａ　 　 　 　 　 　 文章编号:０４９４￣０９１１(２０１７)０８￣０１０６￣０４ Ｃａｌｃｕｌａｔｉｏｎ Ｍｅｔｈｏｄｓ ａｎｄ Ａｃｃｕｒａｃｙ Ａｎａｌｙｓｉｓ ｏｆ Ｍｅｒｉｄｉａｎ Ａｒｃ Ｌｅｎｇｔｈ １ ２ ＬＩＵ Ｘｕｅｊｉｅ ꎬＹＡＮＧ Ｌｉｋｕｎ (１. Ｚｈｏｎｇｗｅｉ Ｓｕｒｖｅｙｉｎｇ ａｎｄ Ｐｌａｎｎｉｎｇ Ｉｎｆｏｒｍａｔｉｏｎ Ｅｎｇｉｎｅｅｒｉｎｇ Ｃｏ.ꎬＬｔｄ. ｏｆ Ｈｅｎａｎ ＰｒｏｖｉｎｃｅꎬＪｉａｏｚｕｏ４５４０００ꎬＣｈｉｎａꎻ ２. Ｚｈｅｎｇｚｈｏｕ Ｉｎｄｕｓｔｒｙ ａｎｄ Ｔｒａｄｅ ＳｃｈｏｏｌꎬＺｈｅｎｇｚｈｏｕ４５０００７ꎬＣｈｉｎａ) Ａｂｓｔｒａｃｔ:Ｔｈｅｒｅ ａｒｅｓｅｖｅｒａｌｋｉｎｄｓｏｆ ａｌｇｏｒｉｔｈｍｓｆｏｒｃａｌｃｕｌａｔｉｎｇｔｈｅｍｅｒｉｄｉａｎａｒｃｌｅｎｇｔｈｅｘｃｅｐｔｔｈｅｃｌａｓｓｉｃａｌｅｘｐａｎｄｅｄａｌｇｏｒｉｔｈｍꎬｓｕｃｈａｓ ｎｕｍｅｒｉｃａｌ ｉｎｔｅｇｒａｔｉｏｎ ａｎｄｏｒｄｉｎａｒｙ ｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓｎｕｍｅｒｉｃａｌ ｓｏｌｕｔｉｏｎ.Ｉｎｏｒｄｅｒｔｏｓｔｕｄｙｔｈｅａｃｃｕｒａｃｙｏｆ ｅａｃｈａｌｇｏｒｉｔｈｍꎬｔｈｉｓｐａｐｅｒ ｓｅｌｅｃｔｅｄ ３ ｓｅｔｓ ｏｆ ｓａｍｐｌｅ ｄａｔａ ｗｉｔｈｉｎ ｔｈｅ ｇｅｏｄｅｔｉｃ ｌａｔｉｔｕｄｅ ｆｒｏｍ ０° ｔｏ ９０°ꎬｗｈｏｓｅ ｉｎｔｅｒｖａｌｓ ａｒｅ １°ꎬ１′ꎬ１″ꎬｒｅｓｐｅｃｔｉｖｅｌｙ. Ｂａｓｅｄ ｏｎ ｔｈｅ ｔｒａｄｉｔｉｏｎａｌ ｅｘｐａｎｄｅｄ ａｌｇｏｒｉｔｈｍｓꎬ ｎｕｍｅｒｉｃａｌ ｉｎｔｅｇｒａｌ ａｌｇｏｒｉｔｈｍｓ ａｎｄ ｎｕｍｅｒｉｃａｌ ｓｏｌｕｔｉｏｎ ｏｆ ｏｒｄｉｎａｒｙ ｄｉｆｆｅｒｅｎｔｉａｌ ｅｑｕａｔｉｏｎｓꎬ ｔｈｅ ｃｏｒｒｅｓｐｏｎｄｉｎｇ ｍｅｒｉｄｉａｎ ａｒｃｌｅｎｇｔｈｒｅｓｕｌｔｓｗｅｒｅ ｃａｌｃｕｌａｔｅｄ ａｎｄｔｈｅ ｑｕａｌｉｔｙ ｏｆ ｅａｃｈ ｎｕｍｅｒｉｃａｌ ａｌｇｏｒｉｔｈｍｗｉｔｈｒｅｇａｒｄｔｏ ａｌｇｏｒｉｔｈｍ ａｃｃｕｒａｃｙ ａｎｄ ｃｏｍｐｕｔａｔｉｏｎ ｓｐｅｅｄｗｅｒｅ ｅｖａｌｕａｔｅｄ.Ｔｈｅｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ３ａｎｄ４ ｏｒｄｅｒ Ｒｕｎｇｅ￣Ｋｕｔｔａ ａｌｇｏｒｉｔｈｍ ｎｏｔ ｏｎｌｙ ｈａｖｅｈｉｇｈｐｒｅｃｉｓｉｏｎｂｕｔ ｔｈｅ ｃｏｍｐｕｔｉｎｇ ｓｐｅｅｄ ｉｓ ｔｗｉｃｅ ｍｏｒｅ ｔｈａｎ ｏｔｈｅｒ