数值计算方法在高斯投影反算中的应用.pdfVIP

第3期 李学民．数值计算方法在高斯投影反算中的应用 75 文章编号：1672—8262(2008)03—75—03 中图分类号：P226 文献标识码：B 数值计算方法在高斯投影反算中的应用 李学民 (河南省测绘工程院，河南郑州 450003) 摘 要：以子午线弧长的数学模型为基础，构造了求底点纬度的牛顿迭代算法，并对牛顿迭代法的适用条件进行了说 明，经过实例计算，该方法具有计算简便、精度高等特点。 关键词：高斯投影；底点纬度；牛顿迭代法 l 引 言 这就是牛顿迭代法。牛顿迭代法在非线性方程的 在大地测量中，高斯投影反算问题是一个重要的 解算中应用广泛，具有显著的优越性。牛顿迭代法具 测量学问题。通过高斯投影反算公式可以看出底点纬 有局部收敛性，在实际问题中，当实际问题本身能提供 度B 是一个非常重要的参数。当B 的值确定以后， 接近于根的初始近似值时，就可保证迭代序列收敛；当 高斯投影反算问题就迎刃而解。在各种文献中，一般 初值难以确定时，由(2)式确定的就不一定收敛。为 是构造直接计算B 的公式。在实际工作中，还可以以 此，需考虑牛顿迭代法的非局部收敛性。下面给出牛 顿迭代法非局部收敛的条件： 子午线弧长的数学模型为基础来确定B 的值。这样 设_厂( )在有根区间[。，b]上二阶导数存在，且满 求B 的问题就转化为求一元非线性方程的解。非线 足： 性方程的求解是一个常见的问题，通常是将非线性方 。)八b)0； 程线性化为级数形式后求解。但这种方法对于精度要 ②厂 ( ) 0， ∈[。，b]； 求较高时，线性化过程既费时费力，又非常繁琐、冗长， ③厂 ( )不变号， ∈[。，b]； 很容易出错。现在随着计算机技术的迅猛发展和计算 ④初值 。∈[。，b]且使 ( 。) 。)0。 机编程语言的普及，我们可以借助于比较成熟的数值 则牛顿迭代序列 }收敛于f( )=0在[。，b]内 计算方法，把大部分计算工作交给计算机去完成，把更 的唯一根。条件①保证了根的存在；条件②表明函数 多的精力放到问题的解决方法上。下面就介绍利用非 单调，根唯一；条件③表明_厂( )的图形在[。，b]上凹向 线性方程的数值解法来计算高斯投影反算中B 的详 不变；条件④规定了初值 。的选取条件。 细过程。 2．2 求B 的牛顿迭代法 曰 r 以子午线弧长公式 =。(1一e )I (1— 2 底点纬度B 的牛顿迭代法