论文推荐｜陈成：辅助纬度与大地纬度间的无穷展开

《测绘学报》

构建与学术的桥梁 拉近与权威的距离

辅助纬度与大地纬度间的无穷展开

陈成¹

, 金立新^2,3, 边少锋¹, 李松林¹

1. 海军工程大学导航工程系, 湖北 武汉 430033;

2. 中铁第一勘察设计院集团有限公司, 陕西 西安 710043;

3. 甘肃铁道综合工程勘察院有限公司, 甘肃 兰州 730000

收稿日期：2018-01-05；修回日期：2018-10-27

基金项目：国家自然科学基金（41631072；41574009；41474061）

第一作者简介：陈成(1990-), 男, 博士生, 研究方向为地图投影与地球重力场。E-mail:chenchengzyq@163.com

摘要：利用无穷级数理论和拉格朗日反演定理，详细推导了大地测量和制图学中常用的辅助纬度与大地纬度间的无穷展开，主要表现为参考椭球第一偏心率的幂级数形式。通过建立一系列严格的系数递推公式，得到了等量纬度反解展开式和等角纬度反解展开式；同时，推导了古德曼函数的泰勒展开式，进而得到了等角纬度正解展开式；利用级数除法公式，得到了等距离纬度正解展开式系数的行列式表示。通过比较本文方法与计算机代数系统Mathematica直接推导求得的辅助纬度正反解展开式e⁰~e⁴⁰阶系数和相应的程序用时，表明本文算法是正确的、快速的。以CGCS2000参考椭球为例，对辅助纬度正反解进行了算例分析，也进一步验证了本文公式的正确性。

关键词：辅助纬度 大地纬度 偏心率 拉格朗日反演定理 古德曼函数

Infinite expansions of the auxiliary latitudes with respect to the geodetic latitude

CHEN Cheng¹, JIN Lixin^2,3, BIAN Shaofeng¹, LI Songlin¹

1. Department of Navigation, Naval University of Engineering, Wuhan 430033, China;

2. China Railway First Survey and Design Institute Group Co. Ltd., Xi’an 710043, China;

3. The General Engineering Survey Institute of Railways of Gansu Co. Ltd., Lanzhou 730000, China

Foundation support: The National Natural Science Foundation of China (Nos. 41631072; 41574009; 41474061)

First author: CHEN Cheng(1990—), male, PhD candidate, majors in map projection and earth gravity field. E-mail:chenchengzyq@163.com.

Abstract: By using the theory of infinite series and the Lagrange inversion theorem, the infinite expansions of the auxiliary latitudes in geodesy and cartography with respect to the geodetic latitude are given, which are the power series of the first eccentricity of the reference ellipsoid. The inverse solutions of the isometric and the conformal latitudes can be obtained from the recurrences of a number of known coefficients. The forward expansion of the conformal latitude is proposed on the basis of the Taylor expansion of the Gudermannian function, and the forward expansion of the rectifying latitude is given in which the coefficient can be expressed as a explicit determinant by the division of power series. Compared the coefficients of the auxiliary latitudes expansions at ordere⁰toe⁴⁰with the results by the computer algebra system Mathematica and the corresponding cost time, the algorithm of this paper is verified to be correct and fast. By taking CGCS2000 reference ellipsoid, the numerical experiment is made for the expansions and then the correctness is further proved.

Key words: auxiliary latitude geodetic latitude eccentricity Lagrange inversion theorem Gudermannian fuction

在大地测量和制图学中，除了大地纬度之外，还有6种常用的辅助纬度，它们都有着广泛的应用^[1-3]。如等量纬度常用于墨卡托投影，地心纬度和归化纬度常用于椭球的几何学，等角纬度、等距离纬度和等面积纬度常用于双重投影^[4-6]。在进行投影转换等分析计算时，经常用到它们与大地纬度间的解析展开式，众多国内外学者对此展开了研究^[7]。地心纬度、归化纬度和大地纬度间正反解函数关系式较简洁，文献[8]利用拉格朗日共轭级数法严格得到了相应的无穷级数展开。而其余4种辅助纬度与大地纬度间关系式较复杂，除等量纬度、等角纬度与大地纬度间正解具有明显的函数关系式外，一般表现为第一偏心率的幂级数形式，或者大地纬度的三角级数形式。由于直接推导的计算量较大，文献[9—10]也只分别得到了展开至e′⁶、e⁸阶等量纬度关于大地纬度的反解展开式。文献[11]首次系统研究了辅助纬度与大地纬度间的展开式，并给出了至e⁸阶的正反展开式。由于计算机代数系统可以进行符号运算，一些较复杂的人工推导可以交由计算机完成。因此，借助Mathematica、Maxima等计算机代数软件，文献[12—13]及其他文献(
http://geographiclib.sourceforge.net/html/auxlat.html;
https://www.academia.edu/7580468/Funciones_de_Latitud)得到了更高阶的等角纬度、等距离纬度和等面积纬度关于大地纬度正反解展开式。然而，除了地心纬度和归化纬度，其他几种辅助纬度与大地纬度间的无穷展开仍然没有给出，上述文献致力于直接推导求解有限阶的辅助纬度展开式，也没有给出展开式系数的一般公式。虽然这类无穷展开表现为幂级数和三角函数组合的形式，但利用傅里叶级数方法仍然很难得到具体的无穷展开式。另外，由于等面积纬度与大地纬度间的关系式存在多层函数嵌套，也很难得到等面积纬度与大地纬度间的无穷展开式。由此，本文拟通过泰勒展开定理和拉格朗日反演定理，推导等量纬度、等角纬度和等距离纬度与大地纬度间的无穷展开式，并取CGCS2000参考椭球，对辅助纬度正反解展开式进行算例分析。

1 等量纬度展开式

由文献[4, 14]，等量纬度q与大地纬度B的关系式为

(1)

式中，e是椭球第一偏心率，为一个小参数，gd^-1(x)=arctanh(sinx)，为古德曼函数的反函数^[15-16]。当大地纬度趋近于极点时，等量纬度趋向于无穷大。在进行数值计算时，为了避免舍入误差，gd^-1(x)通常修正为gd^-1(x)=arcsinh(tanx)。

根据反双曲函数的级数展开，等量纬度可展开成无穷级数

(2)

等量纬度与大地纬度的反解可通过等角纬度与大地纬度的反解展开式和等量纬度与等角纬度的闭合关系式而得到，但是等量纬度与大地纬度的反解也有直接关系式。为求得等量纬度与大地纬度的反解展开，应用拉格朗日反演定理^[17]，有

(3)

式中，gd(x)=arcsin(tanhx)，为古德曼函数^[15-16]。同样，为了避免舍入误差，取计算式gd(x)=arctan(sinhx)。

现在来求得一个有用的级数恒等式，利用组合函数和幂级数的性质^[16]，可得

(4)

式中，系数ζ_mk用递推形式表示为

(5)

部分系数为

(6)

令

(7)

将式(7)和式(4)代入式(3)得

(8)

通过逐次微分运算，可以建立递推式

(9)

进一步可得系数递推式

(10)

式中，初始值F₀⁰=1，0≤t≤m－1，r=－(t+1)+2r′，其中0≤r′≤t+1，－(t+1)≤r≤t+1。当|r|>t时，F_r^t=0，特别地

(11)

(12)

式(10)中的系数递推关系可用图解表示，如图 1所示。

图 1 F_r^t递推关系Fig. 1 Recursion relations for F_r^t

部分系数F_r^t为

(13)

因此，有

(14)

或者

(15)

式中，系数

(16)

相对于式(10)、式(13)中系数F_r^t的m、k，式(16)中对应的为m－k+s、k－s。式(14)和式(15)即为等量纬度与大地纬度的反解展开式，给出了一般项系数的计算式，是一种一般式。实际应用中，为了避免舍入误差，一般展开到e¹²或者n⁶(n为第三扁率^[14])，系数η_mk可用矩阵表示为(0≤k≤m－1)

(17)

特别地，有ηm0=F－m+1m－1ζm0=1，因此

(18)

式中，e′为椭球第二偏心率。

n与e的关系式为

n的幂级数，可借助

(19)

文献[9-10]分别给出了到e′⁶、e⁸阶等量纬度与大地纬度反解展开式的精度，文献[10]还讨论了反解展开式的精度。文献[18-20]也分别用数值方法、解析方法对大地纬度反解等量纬度作了一定的研究。利用第二偏心率与第一偏心率的关系式和恒等式sech²q=1－tanh²q，可验证文献[9—10]的结果分别与本文展开至e⁶、e⁸的结果一致。

2 等角纬度展开式2.1 古德曼函数的泰勒展开

根据泰勒展开定理，古德曼函数在有一增量Δ时，可表达为

(20)

式中，gd(m)(x)为古德曼函数对自变量的m阶导数。

利用导数递推公式和数学归纳法，容易得到

(21)

式中，系数G₁₁=－H₁₁=1，当k=0或k>m时G_mk=H_mk=0，其他情况可由下列递推公式计算

(22)

或者

(23)

特别地，有

(24)

以及一些低阶系数

(25)

(26)

反正弦函数也有类似的泰勒展开，可用于求解等面积纬度的展开式。但由于等面积纬度与大地纬度的关系式存在多层函数嵌套，很难求得一个简洁的无穷展开式，需借助Mathematica或者Maxima软件来求解一定阶的级数展开式。

2.2 等角纬度的正解展开式

根据文献[4, 14]，等角纬度φ与大地纬度的关系式为

(27)

将式(2)、式(4)代入式(27)，得

(28)

不难得到

(29)

(30)

(31)

式中，C_pr=(－1)pC2m+1m－pC2k+1k－r，0≤p≤m，0≤r≤k。

取当且仅当m为奇数且l=(m+1)/2时δ_ml=0，其余情况为δ_ml=1，以及

(32)

利用恒等式sinh(gd^-1B)=tanB和cosh(gd^-1B)=secB，连同式(20)、式(28)—式(31)，可得

(33)

式中，系数

(34)

式中，相对于式(31)系数C_pr中的整数m、k，此处应用m+s－l－1、l－s替代；展开到e⁸阶的系数参考文献[11]，展开到e¹⁰阶的系数和精度分析参考文献[13]，本文计算结果均与其一致，并补充e¹²阶的系数

(35)

2.3 等角纬度的反解展开式

将式(27)代入式(14)，并利用恒等式tanhq=sinφ和sechq=cosφ，有

(36)

因此

(37)

式中

(38)

或者

(39)

展开到e¹⁰阶的系数和精度分析仍然可以参考文献[13]，本文计算结果与其一致，并补充e¹²阶的系数

(40)

3 等距离纬度展开式

根据文献[4, 21]，等距离纬度ψ与大地纬度的关系式为

(41)

式中，X为子午线弧长，R为等距离半径，分别为(取椭球常半轴为单位长度)

(42)

(43)

根据幂级数展开或者傅里叶级数展开方法^[21-22]，可得

(44)

(45)

式中，系数x₀=B,r₀=1，

(46)

(47)

从式(44)、式(46)可以看出，子午线弧长X展开式关于大地纬度线性项的系数就是等距离半径R。因此，根据级数除法公式^[16]，有

(48)

式中，系数

(49)

(50)

或者

(51)

由于r_k=xk0(k≥1)，式(51)中行列式第一列r_kx₀－r₀x_k关于B的线性项r_kx₀－r₀xk0B为0，对b_m的计算没有影响，这也是显然的。因此，在计算b_mk时，可以去除所有关于B线性项和正弦项sin2lB(l≠k)，再通过行列式的逐次运算，可进一步得到

(52)

特别地

(53)

展开到e¹⁰阶的系数和精度分析参考文献[13]，本文计算结果与其一致，并补充e¹²阶的系数

(54)

对于等距离纬度关于大地纬度的反解展开式，可以根据正解公式，运用拉格朗日反演方法，或者直接建立常微分方程

(55)

再根据小参数展开的庞加莱方法^[23]

(56)

等距离纬度反解展开式的计算量很大，也很难得到比较简洁的递推公式，必须借助Mathematica或者Maxima软件计算得到一定阶的展开式。展开到e¹⁰阶的系数和精度分析可参考文献[13]，本文补充e¹²阶的系数

(57)

4 算例分析

为了进一步验证本文辅助纬度与大地纬度间无穷展开式的正确性，可将本文方法(递推法)得到的高阶展开式与计算机代数系统直接推导(直接法)得到的结果进行比较。其中，在利用计算机代数系统直接推导正解展开式时是对原函数进行泰勒展开，而在反解时直接应用拉格朗日反演定理。显然，本文求解展开式系数的递推方法也可以利用计算机代数系统编程实现。为了节省程序运行时间，应用拉格朗日反演定理时应先把三角级数乘积化简成倍角形式再进行求导，而在递推求解等距离纬度正解展开式时，应用如下递推公式替代式(46)、式(47)

声明：本站部分文章及图片转载于互联网，内容版权归原作者所有，如本站任何资料有侵权请您尽早请联系jinwei@zod.com.cn进行处理,非常感谢！