文档介绍：两种大地子午线弧长计算方法的比较
胡洋
长安大学西安2604070210
摘要：在研究与大地椭球体有关的一些测量计算时，例如研究高斯投影计算和 弧度测量计算时，往往要用到大地子午线弧长，而大地子午线弧长的计算公式 涉及到椭圆积分，不能两种大地子午线弧长计算方法的比较
胡洋
长安大学西安2604070210
摘要：在研究与大地椭球体有关的一些测量计算时，例如研究高斯投影计算和 弧度测量计算时，往往要用到大地子午线弧长，而大地子午线弧长的计算公式 涉及到椭圆积分，不能用普通方法求出被积函数。采用变步长辛普森公式求它 的积分和按泰勒级数展开采用普通方法求出被积函数是两种截然不同的计算方 法，本文就两种不同的计算方法作出比较，用程序编程实现从而得到计算结果， 验证了两种方法的正确性和可靠性以及按泰勒级数展开的方法求解大地子午线弧长的优 之处。
关键字：大地子午线弧长，变步长辛普森公式，泰勒级数。
引言：
计算地球椭球子午线的弧长，是大地测量、天文测量、航空航天技术以及地理信息处理技术 中的一项基本内容。计算子午线的弧长涉及到椭圆积分，所以一般是将其展开成级数形式， 再用逐项积分的方法求出满足一定精度要求的计算公式。本文在引出子午线弧长的计算公式 之后，介绍了两种计算方法的原理，并用C语言程序实现了两种算法的电算，继而对两种 算出的结果进行了分析讨论。
一、子午线弧长的计算公式
我们知道，子午椭圆的一半，端点与极点相重合，而赤道又把子午线分成对称的两部分，因 此，推到从赤道开始到已知纬度B之间的子午线弧长的计算公式即可。
如下图，今取子午线上的某微分弧段PP'= d x，令P点纬度等于B, P'点纬度为B+ dB, P 点的子午圈曲率半径为M，于是有dj MdB，因此，为了计算从赤道开始到任意纬度B的
“ r B a（1 一e2）
平行圈之间的弧长，必须求出下面的积分值X = j °Md（式2）。又知M=（1 ..）切
（式3），故（式2）符合椭圆积分的定义（椭圆积分定义为可以表达为如下形式的任何函数
x
的积分：f （x） = j 0#[t,P（t）]dt，其中R是其两个参数的有理函数，P是一个无重根的3或 4阶多项式的平方根，而c是一个常数），因此（式1）为椭圆积分，无法用普通的积分方法 求出原函数。
、两种方法计算子午线弧长的原理
（一）变步长辛普森公式求子午线弧长的积分
a (1 一 e 2)
令M(B)= ，计算步骤如下:
V (1 - e2 sin2 B)3
1、
用梯形公式计算T = h[M（0） + M（B）]/2，其中n=1，h=B。且令Sn=Tn
n
2、
3、
— 1 _ h n-1 一
用变步长梯形法则计算T2广2 T2n + 2 M (B + h /2) k=0
S = (4T -T )/3
用辛普森公式计算2n ' 2 n n
4、
若IS - S l>8，则令2 n n n，h/2 n h，转2继续进行计算；否则结束即为所求的 2n n
积分值。8是事先给定的积分精度。
a (1 - e 2)
，
（二）泰勒级数展开后按普通方法积分
令 M(B)=-
V (1 - e2 sin2 B)3
1、按泰勒级数展开如下:
M - m + m sin2 B