Casiofx-4500PA计算器的道路中线测设程序.doc

Casio fx-4500PA计算器的道路中线测设程序 黄卫华1 张英2 （1 江西省高等级公路管理局 南昌 330013） （2 江西省南昌市公路局 南昌 30006） 摘 要：本文介绍Casio fx-4500PA计算器配合全站仪测设道路中心线的方法，和常用到的计算程序，供大家参考使用，特别是需要完整代码的测设人员，免去从头学习编写的烦恼。 关键词：信息工程；坐标法;中线放样;计算程序代码 0 前 言 本文就坐标法测设道路中心线时用到的计算方法，利用Casio fx-4500PA计算器配合全站仪（或测距仪）在道路中线测设的计算程序作详细介绍。 基本测设原理 坐标法放样就是置全站仪于导线点或某一已知坐标的点上，用极坐标的方法测设中线。如图1所示，将仪器置于导线点D0，以导线点D1作为定向后视点，欲标定中线上一点M，只要知道θ角和距离S即可。当D0（x0，y0）、D1（x后视，y后视）和M（x1，y1）坐标已知，则可以求出θ和S值。 点M的坐标需要根据平面线形理论进行计算，一般由设计文件图纸提供，但是实地的情况复杂，有时图纸并不详细，要另外自行计算，而且野外工作气候恶劣，操作不便，不宜携带太多图纸，如果只带一张导线点数据表和直线曲线转角表以及一台计算器，通过计算器计算中线上任一点M的坐标，工作就更加便捷，失误也会减少很多。 图1 程序设计 Casio fx-4500PA计算器的编程空间容量只有1103 steps，变量从A～Z共26个（当然，可以用Defm来扩充，但每扩充一个变量，编程空间容量就会减小8steps，不合算）。为了能在有限的空间内编写出功能较完整，使用起来又方便的代码，就必须优化压缩程序，又不失可读性。程序分为三个模块，模块1是用于计算上面提到的θ和S值，流程图见图2；模块2用于计算中线上任一点M的坐标，流程示意图见图3；模块3用于输入保存单个平曲线的要素。整个程序占用容量812steps。 数据变量对应表 表1 变量名 表 示 内 容 J 交点的桩号 C 交点的X坐标 D 交点的Y坐标 F 交点前进方向的方位角 A 平曲线转角 R 平曲线半径 H 缓和曲线长 O 平曲线是右转时为0，左转时为1 G 起点桩号（一般是前一交点的HZ或YZ点） E 终点桩号（一般是后一交点的ZH或ZY点） M 立仪点D1的X坐标 N 立仪点D1的Y坐标 Z 需测设的中线上一点桩号 上表中的变量值除M、N、Z在模块1——PM FY中输入外，其余在模块3——INDATA中输入，都是曲线要素数据，在使用时不得通过其他方式改变其值，否则会影响计算结果。 结束程序输入立仪器点D1的坐标（X0，Y0） 结束程序 输入立仪器点D1的坐标（X0，Y0） Y是否超限ZBEGIN或ZEND?N输入桩号Z Y 是否超限 ZBEGIN 或ZEND? N 输入桩号Z WHAT ZH=? 桩号Z0? 桩号Z0? N N 调用文件名CO-XY计算桩号Z的坐标输入已知点M的坐标(X1,Y1)Y 调用文件名CO-XY计算桩号Z的坐标 输入已知点M的坐标(X1,Y1) Y 计算S和θ(FUJ) 计算S和θ(FUJ)值 图2 模块1——PM FY流程图显示结果 显示结果 模块3 INDATA保存的数据（曲线交点桩号、坐标，曲线前进方向直线方位角，曲线转角，曲线半径，缓和曲线长） 模块3 INDATA保存的数据（曲线交点桩号、坐标，曲线前进方向直线方位角，曲线转角，曲线半径，缓和曲线长） 模块1 PM FY 得到的桩号Z 计算B、P、Q、T、L和ZH、HY、QZ、YH、HZ点的桩号 判断桩号Z是在哪一段范围内 小于 ZH点 大于 HZ点 在HY- YH之间 在第二缓和 曲线上 在第一缓和 曲线上 结束模块 Lbl1 计算Z的坐标（X，Y） Lbl6 计算Z的坐标（X，Y） Lbl25 算Q、P Lbl25 算Q、P Lbl34 算L、Q、P 计算Z的坐标（X，Y） Lbl CS计算S、U 计算V、W、K 图3 模块2 图3 模块2——CO-XY流程示意图 3 Casio fx-4500PA计算源程序 现将代码列出： Filename：PM FY L1 {MNZ}：M″X0=″：N″ Y0=″：Z″WHAT ZH″ L2 Z＜0={XY}：X″ X1=″：Y″ Y1=″：≠Fixm：ZG=″Begin″▲Goto E：≠ ZE =″End″▲Goto E：≠Prog CO-XY△△△ L3 Pol((X-M),(Y-N))：V″ S″=V▲W0=W=360+W：≠W△W″ FUJ″=W▲ L4 Lbl E Filename：CO-XY L1 Fixm：B=180H/(2πR)：P=HH/(24R)：Q=H/2-HH