论道路坐标计算积分计算的问题01——为什么有积分计算的问题.docVIP

论道路坐标计算中积分计算的问题01——为什么有积分计算的问题 2009年05月29日 00:31 事情先要从一个人发来的电子邮件和一本书说起。 5月21日，我收电子邮件，连续两封来自测量工程师同行的措辞犀利的信给了我极大的震撼，之所以震撼，主要还是他比较极端的语气，他对我去年出版的《CASIO fx-5800P计算器与道路坐标放样计算》一书给予了几乎全部的否定，不过，我还是怀着虚心的心境与他进行了若干次的电子邮件交流，前前后后有十来封电子邮件。具体内容这里不表，有兴趣了可查看我的日志：一位测量技术“批评家”的信。 我这里要说的是他在信中提到的一个测量工程事故，为说明清楚，我整理如下： ……让我勃然大怒的是您书中坐标正反算采用积分公式，这是很多从事测量工作不深的人摔跟头的地方，此公式让他们看到计算曲线太简单了，从数学上讲积分是严密的，数值分析是可靠的，问题出在卡西欧计算器上，计算器采用复化辛普森公式，对积分区间划分1~9段，这是随机的，当计算曲线长，半径大，划分段数小时误差大得吓死人，潭邵高速有两座桥曲线编程采用积分公式计算测量的，误差达1.2米，全桥报废。因此看到您的编程采用113页4-64公式时骂您误人子弟，因为您的书发行了，上千万人在学习它，两句语句就能解决核心程序，让没有深入研究的人欣喜若狂，认为它是可靠的也是可以采用的，结果也是危险的。…… ……那是多年以前的事，只记得是潭邵***标，那两座桥处在缓和曲线和圆曲线上，半径很大很长，至少有6公里，测量是吴***，总工是邹***，吴当时刚从学校毕业三年，计算器是CASIO4800，编程用的是计算器积分公式，这对于一个刚毕业的学生来讲选择它编程很简单，放线后河流改道围堰施工，把柱子盖梁全做好了，架桥机往上架梁，发觉梁不是长了就是短了。查原因，加工厂没问题。查测量资料与程序，与会的几个测量工程师算的数据用积分公式与复化Simpson算的数据有多个结果，而用Gauss-Legndre5结点，缓圆曲线参数方程算的结果一样，同行护短，都想把责任推给卡西欧。致电卡西欧中国代理，回复：“说明书说得很清楚，采用高斯近似积分，计算段数取1~9段，计算结果有误差，你们要调整tol值提高精度，你们对计算结果要复核，我们的计算器是近似计算，不是精确计算，6公里长的曲线分几段积分外行，近似懂么”。 编程没错，卡西欧没错，错在那里？错在没有复核，结果总工、测量工程师开除留司察看，禁止直接使用计算器上积分公式编程去放高精度的桥梁。 以上两段，第一段是前一封电子邮件的，后一段是在我看了他前一封电子邮件，请他把当时的数据、具体细节提供给我研究之后，他发给我的，他没有具体参数，但在信中还是提到了几个关键数据，而且经过几轮交流，他也明显对我有所好感，信中的公司名称、项目名称、当事人姓名都用了真实名字（这里为了避免麻烦，都隐去或者使用化名）。 从描述的情节来看，这个测量事故的核心在于计算，而计算的核心又在于积分计算。本日志就说道路坐标的积分计算。 说到道路坐标的计算，是道路施工放样计算的一个最基本、最常见的计算类型，不论是道路主线还是立交匝道，中线无非是三种线型：直线、圆曲线、缓和曲线，前二者有精确的坐标计算公式，第三者麻烦一些，计算公式是无穷多项式，在保证精度的前提下可取前面几项即可，缓和曲线坐标计算之麻烦还在于实际工程中还存在大量不完整缓和曲线（如卵形曲线），这给缓和曲线坐标计算带来更大的难度，曾几何时，卵形曲线的参数计算及坐标计算思路、方法多多，不断见于各种期刊杂志。 其中，有一种采用积分计算道路坐标的数学方法使得以上问题好象迎刃而解：不管哪种线型均使用同一计算公式、局部缓和曲线（卵形曲线）也不例外，具体见本人整理的连载日志： 道路中线坐标计算的通用数学模型01——曲线上任意点的切线方位角和坐标的计算 道路中线坐标计算的通用数学模型02——路线坐标积分计算的数值算法 道路中线坐标计算的通用数学模型03——道路中线外一点的坐标计算 道路中线坐标计算的通用数学模型04——路线外一点反求对应桩号及距离 坐标计算积分的关键数学公式如下： . . 关于此公式的积分计算，对于卡西欧计算器，由于有内部的积分函数，可在程序中使用积分函数计算，而计算机编程语言则没有积分函数，只能使用数值方法完成，同样，卡西欧计算器编程中也可使用数值方法来进行。说句实话，我之前一直钟情于使用卡西欧计算器的积分函数来完成上面公式的计算，因为代码编写简洁，在我编写出版的《CASIO fx-5800P计算器与道路坐标放样计算》一书中也是如此。 让我感到奇怪的是，在我接触到的其他人编写的卡西欧计算器程序，关于这一公式的程序实现，不管是在正式出版的技术论文、书籍教材，还是流传于网上的各种计算程序，