子午线弧长的数值计算方法探讨.pdfVIP

2007年第4期 河南测绘 第86期 子午线弧长的数值计算方法探讨 李学民 田耀永 陈 卫 、J ； 河南省测绘工程院 驻马店 463000 = + 7 — 4 ， 6 一 、 引言 e 一 予午线弧长计算足椭球人地测量的一个基本数学问题，它的计算公式中含仃椭圆积分，椭圆积分无 分析解。人地测量中一般的做法足按二项式定理展开，展开后为积分方便，又需要按有关公式，化正弦 幂函数为余弦的倍角函数，合并同类项后逐项积分。子午线弧长在大地测量和地图制图中有着广泛的用 途，如用于推算地球形状大小的弧度测量，地图投影中的gauss投影计算等。 二、子午线弧长计算的一般求法 从相关的参考书或教科书中可以知道，子午线弧长可表示为如下椭网积分： 3 =口(1一e ) (1一e sin )一jdB 由于这足个椭圆积分，不可能用一般的积分方法求出其解，必须按二项式定理将(1一e sin ) 展 开为级数形式，然后逐步地进行积分。于是有 r1一e sin B 2一三e 2 cos 2B+4_25 e 4一—1—5 e 4 cos 2B+—1—5 e 4 cos 4 + 4 64 16 64 一 一 525e 6 cos2B+ os4 一旦 e 6 cos6B+．．．． 512 256 512 令常系数： ： 1十三e z+竺e + e +．．… 4 64 256 ： 三e z+堕e +堕 e +．．… 4 16 512 ： e + e一+ ．．… 64 256 旦 e +．．… ： 512 3 贝IJ：(1一P sin B) = 一flCOS2B+ycos6B一2cos6B+ =a(1一e 2) ( 一flcos2B+ycos4B一2cos6B+……)dB 积分解X=ac 一鲁sin2B+ sin4B- 6 … ㈩ ． 12． 2007年第4期 }If『南测绘 第86期 (1)式即为赤道 了午线卜某点的了午线弧K公式。代入相应的椭球元素，即可求得该椭球元素下 的某点的了午线弧长。(本义所有公式中的 取弧度值。) 将1975年同际椭球元素a：6378140．00000，b：6356755．28815代入 (1)式，得： 6367452．1242 160