卡西欧5800计算器编程基础学习程序集.pdf

《卡西欧5800计算器编程基础学习程序集.pdf》由会员分享，可在线阅读，更多相关《卡西欧5800计算器编程基础学习程序集.pdf（42页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、 fx-5800P易学易用程序集 卡西欧(中国)贸易有限公司 序 内容 页码 1 目录 第1页 2 易学易用程序集 第3页 3 公路高铁程序 第8页 4 编程基础 第25页 5 实用用户公式集 第26页 6 基本操作 第28页 7 常见问题 第31页 8 编程技巧 第34页 9 数学功能 第37页 目录 1 版程序列表 备注 序 主程序 子程序 功能说明 1 极坐标放样程序（已知坐标计算边长、方位角计算程序） 2 已知边长方位角计算坐标程 3 建筑轴线偏移程序(或扩大基础程序) 4 测角前方交会坐标计算 5 测角后方交会坐标计算 6 坐标转换计算程序 7 直线中边桩坐标计算程序 8 圆曲线中边

2、桩坐标计算程序 9 完整缓和曲线中边桩坐标计算程序 10 ； 线元法匝道连续计算程序（万能程序） 11 ； 交点法曲线连续计算程序 12 竖曲线高程连续计算程序 序 参考书目作者 参考书目名称 1 王中伟教授 道路施工放样综合程序，王中伟教授在QQ（595077）空间把自己的程序电 子版贡献给广大用户，在此特别感谢！ 2 覃辉教授公路与铁路施工测量程序。 3 卡西欧fx-5800P 说明书。 4 其他可能参考的书籍。 5 特别感谢 特别感谢广西路桥的李虎（东华理工学院97届毕业），他无私奉献了他的公路 高铁测量程序。 6 免责申明 本书中程序仅供用户免费参考，卡西欧（中国）贸易有限公司不对用户

3、使用本 书中程序发生的任何问题负责。 7 联系方式用户对用本书疑问，请发邮件到： 。欢迎交流！ 8 客服热线卡西欧的客服热线：400-700-6655 。 9 程序传输本书的程序用户可以到全国省会城市的测绘店免费传输G版程序。 10 程序征集 如果您有非常好的实用程序，希望广大用户能免费得到分享，请和本书编者 联系:。 11 欢迎翻印 本书欢迎翻印；希望更多的用户能得到，欢迎流通。本书电子稿及其他的测量 程序源代码请到卡西欧网站下载: 下载中心计算器工程资料下载，或者到百度文库中搜索fx-5800P卡西 欧易学易用程序集下载。 G1 G2 G3 G4 G5 G6 G7 G8 G9 G10 G1

4、1 G11-1 G12 G12-1 G10-1G10-2 fx-5800P G2011-9-30 . 常用测量程序 1 极坐标放样程序（已知坐标计算边长、方位角计算程序） ，程序名 1 Deg：Fix 3 设置角度单位为十进制，3 位固定小数显示 2 XO=?A：YO=?B 提示输入起点的坐标（A,B） 3 XN=?C：YN=?D 提示输入端点的坐标（C,D） 4 Pol(CA),(DB) 调用极坐标函数计算距离和角度,自动保存到 I 和 J 5 J0J+360J 把小于 0 的角度换算成方位角 6 JL=：I 显示水平距离 I 7 FWJ=：JDMS 显示方位角 J 2 已知边长方位角计算坐

5、标程序，程序名 1 Deg：Fix 3 设置角度单位为十进制，3 位固定小数显示 2 XO=?A：YO=?B 提示输入起点的坐标（A,B） 3 L=?L：C=?C 提示输入所求点的距离 L 和方位角 C 4 A+Lcos(C) X 计算所求点的 X 坐标 5 B+Lsin(C) Y 计算所求点的 Y 坐标 6 X=：X 显示所求点的 X 坐标 7 Y=：Y 显示所求点的 Y 坐标 3 建筑轴线偏移程序(或扩大基础程序)。程序名 1 Deg：Fix 3 设置角度单位为十进制，3 位固定小数显示 2 ?A:?B:?C:?D: ?E 提示输入起点的坐标（A,B） 、起点方位角 C、轴线偏移距 D,E

6、. 3 A+Dcos(C) + Ecos(C-90) X 计算所求点的 X 坐标 4 B+D sin(C) + Esin (C-90) Y 计算所求点的 Y 坐标 X=：X 显示所求点的 X 坐标 6 Y=：Y 显示所求点的 Y 坐标 7 （） （） 易学易用程序易学易用程序集集 3 4 测角前方交会坐标计算，程序名 1 Deg：Fix 3 设置角度单位为十进制，3 位固定小数显示 2 ?A:?B:?C:?D: ?E : ?F 提示输入两个已知点的坐标（A,B）, （C,D） ，观测的两 个角度 E 和 F 3 (A(tanF)-1)+C(tanE)-1+(D-B)(tanE)-1+(tanF

7、)-1) X 计算所求点的 X 坐标 4 (B(tanF) -1)+D(tanE)-1+(A-C)(tanE)-1+(tanF)-1) Y 计算所求点的 Y 坐标 5 XP=：X 显示所求点的 X 坐标 6 YP=：Y 显示所求点的 Y 坐标 7 5 测角后方交会坐标计算，程序名 1 Deg：Fix 3：11DimZ 设置角度单位为十进制，3 位固定小数显示 2 ?A:?B:?C:?D: ?E : ?F:?G: ?H : ?K 提示输入两个已知点的坐标（A,B）, （C,D） ， （E,F） 观测的角度 G,H,K 3 Pol(CA),(DB) : J0J+360J : JZ4 计算点（）的方

8、位角 4 Pol(EA),(FB) : J0J+360J : JZ5 计算点（）的方位角 5 Z4- Z5Z1 方位角计算顶角 6 Z10 Z1+360Z1 计算点（）的顶角 并存贮到变量 7 Pol(EC),(FD): J0J+360J : JZ4 计算点（）的方位角 8 Pol(AC),(BD): J0J+360J : JZ5 计算点（）的方位角 9 Z4-Z5Z2 方位角计算顶角 10 Z20Z2+360Z2 计算点（）的顶角 并存贮到变量 11 Pol(AE),(BF) : J0J+360J : JZ4 计算点（）的方位角 12 Pol(CE),(DF) : J0J+360J : JZ

9、5 计算点（）的方位角 13 Z4-Z5Z3 方位角计算顶角 14 Z30：Then 983935.12C:4774.384D: 2415.861E:140F:311723G:410.007H:600R:IfEnd If L1060：Then 1295658.19C:4206.421D: 3093.946E:70F:-335048G:1285.437H:600R:IfEnd 3DimZ： “X0”?A： ”Y0”?B： “Ln”?L：Abs(L) Int(Abs(1000L)1000Z*3+ ：Lbl 0：If Z*3+0.0001：Then ?L：E lse “Xp”?X： ”Yp”?Y：X

10、Z*1+：YZ*2+：IfEnd：Lbl 1:If L0 ：Then 983935.12C：4774.384D：2415 .861E：140F：311723G：410.007H：600R：IfEnd：If L1060 Then 1295658.19C：4206. 421D：3093.946E：70F：-335048G：1285.437H:600R：IfEnd：Prog”XLJS” ：If Z3+0.0001： Then Goto 0：Else (Z1X)cos(O)(Z2Y)sin(O)N：LNL：Int(Abs(1000L)1000+0.0001L：( Z1X)sin(O)+(Z2Y)cos

11、（O）K：If Abs（N）0.001：Then Goto 1：Else “L=“：L “K=” ：K IfEnd：Goto 0: IfEnd 数据文件的编写数据文件的编写 9 如果还有其他弯道，可以继续后加。在这些存放设计参数的语句程序表达式是固定的。编辑曲线参数时，每 个曲线参数放在一个 If L*.*（两个弯道中间直线段上的任意桩号） :Then *.*C（弯道 起点方位角） ：*.*D（交点 X 坐标） ：*.*E（交点 Y 坐标） ：*.*F（缓和曲线长度：* *G （转角， 向左转为负值， 向右转为正值） ： *.*H （交点里程桩号） :* R （圆曲线半径） ： IfEnd 条

12、件式语句里。如果有多个弯道，一直按上述形式编下去。 变量说明： 序 语句 说明 1 L*.* 弯道参数的起点范围，可以是弯道起点，也可以是弯 道前直线段里的任意点桩号。 2 *.*C 弯道起点方位角 3 *.*D、*.*E 交点X坐标、Y坐标 4 *.*F 缓和曲线长度，如果没设臵缓和曲线，则缓和曲线长度看做 0 5 *G 转角，向左转为负值，向右转为正值 6 *.*H 交点里程桩号 7 *R 圆曲线半径 【注意】运行程序时，第一步是输入测站坐标 X0，Y0，然后输入近似桩号，主要是为了加快收敛， 输入时注意当根据桩号和距中距离计算所求点坐标时，请不要把小数点第四位输入 0.0001，如 K2

13、1+369.0921，请把桩号输入 21369.092 即可，相反，当根据坐标反算桩号和距中距离时，请在输入 时把第四位小数数位 0.0001，如 K21+200，输入时请输入 21200.0001。 10 （二）线元法数据文件编辑 线元法数据文件编辑适合于各种各样的线型组合，特别是在大型立交区或者是不规则的线型组合地方， 如小区道路、重丘区山路等有特别明显的有点，本程序为了减少使用者设计参数的输入，特别把直线与缓和 曲线以及圆曲线连接的地方只要输入缓和曲线参数就可以了，这样可以减少输入很多设计参数，同时节约大 量的内存。看下图： 该图是一个半喇叭型的互通立交，一共是 5 条匝道，我们以 D

14、匝道为例，来看看线元法数据文件的建立。 11 3DimZ： “X0”?A: “Y0”?B: “Ln”?L: Abs(L)-Int(Abs(1000L)Z3:Lbl 0:If Z30.0001: Then ?L:Else “Xp”?X:”Yp”?Y:XZ1：YZ2：IfEnd:Lbl 1:If L0 :Then 3021826W：494347.926X：477989. 873Y：55P: 40R：39.216I: 39.289J：IfEnd: If L78.504：Then 494380.677X：477971.271 Y：3582845.9W：40R： -65P:39.216I: 78.504

15、J：IfEnd:Prog”: Prog”: If Z30.0001:Then Goto 0:Else (Z1X)cos(O)(Z2Y)sin(O)N:LNL:-(Z1X)sin(O)+(Z2Y)c os（O）K：If Abs（N）0.001：Then Goto 1：Else “L=“:L “K=”:K IfEnd:Goto 0: IfEnd 下表是 D 匝道的线元表： D 匝道线元数据分析：从上表中可以看出，D 匝道由三个线元组成，第一段是缓和曲线段，根据上表判 断，该缓和曲线属于第一段缓和曲线，其最小半径处桩号为 J=K0+039.289，该处坐标 X=494347.926,Y=477989

16、. 873， 半径为 R=40 米， 方位角 W=3021826.0， 缓和曲线旋转常数 P=55， 在缓和曲线终点 K0+039.289 处， 有一圆曲线与其相连接，长度 I=39.216。第二段是圆曲线，第三段又是一段缓和曲线，根据上表判断，该曲线 属于第二缓和曲线， 曲线最小半径位于与园曲线连接处， 桩号J=K0+78.504， 坐标为X=494380.677， Y=477971.271, 方位角 W=3582845.9， 最小半径 R=40， 缓和曲线旋转常数 P=-65 （为了判断是第一段还是第二缓和曲线， 我们通过 P 的正负来判断，第一缓和段 P 统一用正数表示，第二缓和段统一用

17、负号表示。起点处相连的圆曲线 长度 I 可以看做 0，也可以看做是 39.216.由此我们整理下。从整个线型全部是向右转向，所以 R 全部取正数。 具体如下： 第一段缓和曲线线元参数说明： 第二段缓和曲线线元参数 最小半径处坐标 X=494347.926 Y=477989. 873 最小半径 R=40 最小半径处方位角：W=3021826.0 缓和曲线旋转常数 P=55 与之相连的圆曲线长度 I=39.216 最小半径处桩号 J=39.289 最小半径处坐标 X=494380.677 Y=477971.271 最小半径 R=40 最小半径处方位角：W=3582845.9 缓和曲线旋转常数 P=

18、65 与之相连的圆曲线长度 I=39.216 最小半径处桩号 J= W=3582845.9 此段数据计算范围可以从K0+000 K0+78.504 此段数据计算范围可以从 K0+78.504 到匝道结束。 当桩号 L 位于 K0+000 K0+78.504 时，曲线参数按第一段编写，当 L78.504 时，按第二段缓和曲线编 写，当然第二段缓和曲线计算范围也可以从圆曲线起点 K0+039.289 到 K0+78.504 中间任意桩号开始。下面 是 D 匝道的数据程序文件（文件名“”）：（O 表示变量字母 O，其余的 0 均为零） 线元 名称 起点桩号 偏向 半径或缓和曲 线旋转常数 A 线元长

19、度 (m) 线元起点坐标、方位角 X(m) Y(m) 方位角 缓和曲线 DK0+000.00 右偏 55 39.289 494341.902 478027.825 2603855.7 圆曲线 DK0+039.289 右偏 40 39.216 494347.926 477989. 873 3021826.0 缓和曲线 DK0+078.504 右偏 65 105.625 494380.677 477971.271 3582845.9 12 在上面程式中，每个缓和曲线的线元参数放在一个 If Then IfEnd 循环语句中，If 后面紧跟计 算的起点范围，终点范围从下段的起点范围结束。 变量定义：

20、 序 变量 说明 1 X、Y 输入缓和曲线半径最小处的 X、Y 坐标 2 R 缓和曲线最小半径，当线路左转时，半径数据前面加“-”号，以判断线路转向。 3 W 输入缓和曲线半径最小处的方位角 4 P 缓和曲线旋转常数，第一缓和曲线长度为正数，第二缓和曲线常数为负数。当不 设缓和曲线时，P=0 5 I 与缓和曲线最小半径端连接的圆曲线长度。 6 J 缓和曲线小半径段的桩号。 再看看 E 匝道数据文件编辑 E 匝道线元设计数据表 线元名称 起点桩号 偏向 半径或 缓和曲 线旋转 常数 A 线元长 度(m) 线元起点坐标、方位角 X(m) Y(m) 方位角 直线 EK0+000.00 71.296

21、494360.619 477706.293 794957.0 缓和曲线 EK0+071.296 右偏 120 72 494373.205 477776.470 794957.0 圆曲线 EK0+143.296 右偏 200 60.070 494381.631 477847.870 900844.6 缓和曲线 EK0+203.366 右偏 95 45.125 494372.527 477907.018 1072116.3 缓和曲线 EK0+248.491 左偏 75 64.841 494355.878 477948.932 1134905.6 圆曲线 EK0+313.333 左偏 86.75 1

22、7.800 494337.373 478010.658 922419.0 13 3DimZ： “X0”?A:”Y0”?B: “Ln”?L:Abs(L)-Int(Abs(1000L)/1000Z3:Lbl 0:If Z30.0001: Then ?L: Else “Xp ”?X:”Yp”?Y:XZ1：YZ2：IfEnd:Lbl 1:If L0 :Then 494381.6316X：477847.870Y：900844.6 W：200R：120P:60.07I: 143.296J：IfEnd: If L203.366：Then 494372.527X：477907.018Y： 1072116.3W

23、：200R：-95P:60.07I: 203.366J：IfEnd: If L248.491：Then 494337.373X： 478010.658Y：922419W：86.75R：75P:17.8I: 313.333J：IfEnd: Prog”: Prog” ” ：If Z30.0001:Then Goto 0:Else (Z1X)cos()(Z2Y)sin() N:LNL:-(Z1X)sin( )+(Z2Y)cos（）K：If Abs（N）0.001：Then Goto 1：Else “L=“:L ”K=”:K IfEnd:Goto 0: IfEnd 在匝道线元表中，我们可以看出，E 匝

24、道式由一段直线段缓和曲线段圆曲线段缓和曲线段缓 和曲线段圆曲线段连接起来的。值得注意的是，在位于 K0+248.491 处，是两段反向的缓和曲线连接点， 在编写数据程序文件时， 分三段编写， 第一段从K0+000到K0+203.366， 第二段从K0+203.366到K0+248.919， 第三段从 K0+248.919 到线路终点 K0+331.13。具体文件如下：（O 表示变量字母 O，其余的 0 均为零） 14 为了让使用者熟悉掌握数据程序文件的编写线元法的数据程序文件，现在将 C、F、G 等匝道的线元数据 列出来，供使用者练习。C 匝道数据；F 匝道数据；G 匝道数据 线元名 称 起点

25、桩号 偏向 半径或 缓和曲 线旋转 常数 A 线元长 度(m) 线元起点坐标、方位角 X(m) Y(m) 方位角 圆曲线 CK0+000.00 右偏 32 56.876 495005.991 478056.677 561942.3 直线 CK0+056.876 21.277 494991. 262 478104. 123 1580954.4 缓和曲线 CK0+078.153 右偏 65 30.179 494971. 511 478112.036 1580954.4 圆曲线 CK0+108.332 右偏 140 30.431 494943.127 478122.242 1642025.7 缓和曲

26、线 CK0+138.763 右偏 65 30.179 494913.167 478127.219 1764740.6 直线 CK0+168.941 右偏 85 422.023 494883.008 478126.740 182.5811.8 缓和曲线 CK0+590.964 右偏 85 85 494461.552 478104.874 182.5811.8 圆曲线 CK0+675.964 右偏 85 72.739 494379.485 478086.681 2113704.2 直线 FK0+000.00 111.525 494591.311 478420.033 2523200.2 缓和曲线

27、FK0+111.505 右偏 100 90.909 494557.835 478313.650 2523200.2 圆曲线 FK0+202.434 右偏 110 101.114 494542.812 478224.689 2761233.6 缓和曲线 FK0+303.548 右偏 120 130.909 494595.310 478142.420 3285235.7 缓和曲线 GK0+000.000 右偏 175 133.152 494718.802 478112.963 1825811.8 圆曲线 GK0+133.152 右偏 230 149.124 494587.600 478093.36

28、8 1993317.5 缓和曲线 GK0+282.286 右偏 190 149.72 494472.335 478002.890 2364221.1 15 计算程序是核心程序，但是道路计算程序是固定的，使用者不需要做任何修改，也不直接运行道路计算 程序分为两部分，第一部分是常规的对成型线型组合计算的程序（文件名：），交点法编辑的数据文 件直接调用； 第二部分是适用于非对称的线型组合计算文件 （文件名：） ， 在数据文件调用完 “” 文件后，还要调用“”文件。程序代码：（为了输入方便代码分行显示，使用者不必分行输入。） Deg:Fix 4 GAbs(G)V: F2-F(3)（240R2）+F（5

29、）(34560R（4）)-F（7）(8386560R（6）)+F（9）(3158507 520R（8）)Q:F2(24R)-F(4)(2688R(3)+F(6)(506880R(5)-F(8)(154828800R(7)P:Q+(R+P)Vtan(G2) T:RGV180+FI: IF LH-T:Then L-HJ:Goto 1:IfEnd: If LH-T+F:Then L-H+TJ:CZ:1S:Goto 2:IfEnd: If LH-T+I-F:Then L-H+TJ:CZ:1S: Goto 3:IfEnd: If LH-T+I:Then H-T+I-LJ:C+GC:-1S:Goto 2:

30、 IfEnd: If LH-T+I:Then L-H-I+2TJ:C+GC:Goto 1:IfEnd: Lbl 1: D+Jcos (C)X: E+Jsin(C)Y: CO: Goto C:Lbl 2: (J-J(5)(40(RF)2)+J(9)(3456(RF)(4)-J(13)(599040(RF)(6)+J(17)(175472640(RF)(8)-T)SM: (J(3)(6RF)-J(7)(336(RF)(3)+J(11)(42240(RF)(5)-J(15)(9676800(RF)(7)+J(19)(3530097000(RF)(9) VN: C+90J2SV(RF)O: Goto

31、B Lbl 3: 180(J-0.5F) (R)O: Q+Rsin (O)-TM:(P+R(1-cos (O)VN:C+VOO： Goto B Lbl B: D+Mcos(C)-Nsin(C)X: E+Msin(C)+Ncos(C)Y: Goto C: Lbl C:If Z*3+0.0001:Then ?K: M?U:K180-U U:O+UAbs(K)(K+10(-10) Z:X+Abs(K)cos(Z)X: :Y+Abs(K)sin(Z)Y: X:X Y:Y tan-1(Y-B)(X-A+10(-10)W: If X-A0:Then W+180W: IfEnd: If W0:Then W+

32、360W:IfEnd:(X-A)2+(Y-B )2)S:DIS: S FW: W DMS IfEnd 道路计算程序道路计算程序 16 非对称线型组合计算程序（）（O 表示变量字母 O，其余的 0 均为零） 最小半径处坐标 X=494347.926 Y=477989. 873 最小半径 R=40 最小半径处方位角：W=30218 26.0 缓和曲线旋转常数 P=55 与之相连的圆曲线长度 I=39.216 最小半径处桩号 J=39.289 Abs(P+10(-10)）(P+10(-10) S：Abs (R)RV:Abs(R) R:P2RF: 180(I+F)V(R) G:90FSV(R) O:

33、If S=1:Then W-OC: IfEnd: If S=-1:Then W-O-GC :IfEnd:C0 C+360C:C360 C-360C: F2(24R)-F(4)(2688R(3)+F(6)(506880R(5)-F(8)(154828800R(7)P:F2-F(3)(240R2)+F5( 34560R(4)-F(7)(8386560R6)+F9(3158507520R(8)Q: Q+(R+P)Vtan(G2)T:If S=1:Then J-F+TH:CZ: Else J-I-F+TH:C+GZ:IfEnd: (F-F(3)(40R2)+F(5)(3456R(4)-F(7)(599

34、040R(6)+F(9)(175472640R(8)-T)SM:(F2(6R)-F(4) (336R(3)+F(6)(42240R(5)-F(8)(9676800R(7)+F(10)(3530097000R(9)VN: X-Mcos(Z)+Nsin(Z)D:Y-Msin(Z)-Ncos(Z)E 17 程序运行分为两种，第一种是根据根据使用者输入的桩号和距中距离计算出所求点的坐标和测站点到该 点的方位角和水平距离，可以直接利用计算器提供的距离和方位角进行极坐标放样。另一种是根据使用者输 入的实测坐标计算出该点的桩号和距中距离， 可以根据桩号和距中距离进行高边坡、 隧道断面点的设计高程、 以及进行

35、路面高程控制计算方面很方便。 （一）据桩号和距中距离计算所求点的坐标、及方位角和水平距离。要进行程序的线路计算，直接运行数 程序文件就可以了，不需要理会计算程序。如果计算某桩号 L 距离中线 K 处的坐标 X，Y 和该点到测站点的 距离和方位角。运行提示如下： 序 变量 说明 1 X0? 提示输入测站点 X 坐标 2 Y0？ 提示输入测站点 Y 坐标 （测站点的位数和数据程序文件里坐标取位要一致） 3 Ln? 计算点近似桩号，近似桩号输入第四位小数后不能等于0.0001， 4 L？ 输入计算点准确桩号，近似桩号输入第四位小数后不能等于 0.0001。 5 K？ 提示输入距离中桩的桩号如果计算点

36、位于线路的左幅，则在数据前加入“-”号，如 果是计算中桩坐标，则输入 0 6 M？ 输入该点与线路的前进方向的按顺时针转过的角度，常于斜交的涵洞，桥梁等构造 物 7 X 显示所求点的 X 坐标 8 Y 显示所求点的 Y 坐标 9 FW 显示所求点至测站点的方位角 10 DIS 显示所求点至测站点的距离 要继续计算则从 L?开始，计算新点的坐标，具体流程见程序流程图的左边部分。例：要计算 D 匝道的 DK0+020 中桩坐标、DK0+080 左边 5.3 米的坐标、Dk0+240 斜交 755811.8、距中右侧 5 米处的坐标， 并计算出该点至坐标为 X=495005.991，Y=478056

37、.677 控制点的水平距离和方位角。 程序运行程序运行 18 运行文件“”提示如下： 序 变量 说明 1 X0? 495005.991 测站点X坐标 2 Y0？ 478056.677 测站点Y坐标 3 Ln? 0 输入小输点第四位起不等于0.0001的桩号 4 L？ 20 输入所求点桩号 5 K？ 0 输入距中距离 6 M？ 90 输入斜交角度 7 X 494341.482 显示所求点X坐标 8 Y 478007.903 显示所求点 Y 坐标 9 DIS 666.2959 显示所求点至测站平距 10 FW 1841152.42显示所求点至测站方位角 11 L？ 80 输入下点的桩号 12 K？

38、 -5.3 13 M？ 90 14 X= 494382.2290 15 Y= 477965.9594 16 DIS 630.324 17 FW 1881629.5 18 L？ 240 19 K？5.0 20 M？ 755811.8 21 X= 494455.3110 22 Y= 478101.5491 23 DIS 552.505 24 FW 1752029.6 19 （二）据坐标反算桩号和距中，根据某点的坐标，计算出该点在线路中的桩号和距离，程序运行时提示如下： 序 变量 说明 1 X0? 提示输入测站点X坐标 2 Y0？ 提示输入测站点Y坐标（测站点的位数和数据程序文件里坐标取位要一致）

39、3 Ln? 计算点近似桩号，近似桩号输入第四位小数后不能等于0.0001， 4 Xp？ 所求点P的X坐标 5 Yp？ 所求点P的Y坐标 6 L 显示所求点的桩号 7 K 显示所求点的距中距离，如果位于线路的左边，则显示结果为负值，如果K为0时，则该点在线路 中线上。 计算新点的桩号和距中距离，则从?Xp 开始。如此循环计算。具体流程见流程图的右边部分 例：求点 X=494382.2290，Y=477965.9594 在 D 匝道中的位置，运行文件“”如下： 序 变量 说明 1 X0? 495005.991 测站X坐标(可输入0) 2 Y0？ 478056.677 测站Y坐标(可输入0) 3 L

40、n? 60.0001 输入小输点第四位起等于0.0001的桩号 4 Xp？ 494382.2290 输入所求点X坐标 5 Yp？ 477965.9594 输入所求点Y坐标 6 L 79.999999 显示所求点桩号 7 K -5.3000 显示所求点距中距离, “-“ 表示该点位于中桩左侧 20 （一）本程序在高边坡放样中的应用 高边坡放样是目前公路道路测量工作普遍要做的一项工作，边坡坡顶线或者填方边坡坡脚线是边坡放样 最繁琐的工作，利用本程序进行边坡放样，工作步骤如下： 司镜员走到目的地后，先测出一个点的坐标和高程； 利用本程序求出该点在线路中的桩号和距中距离，然后根据计算出来的桩号和设计图

41、纸，内插出道 路边桩高程，根据道路边桩高程和道路设计宽度、坡度、设计碎落平台的宽度和间隔高度、实测点 的高程计算出按实测高程，边坡距中距离。 根据高程算出的距中距离和程序计算出立镜点至距中距离差值 dK，指挥司镜员靠近道路中线或者 往外偏移多少米。 重新做 13 工作，直到 dK 小于边坡宽度控制精度要求为止。 如上图，司镜员在 P 点立镜，测出 P 点坐标和高程后，用本程序计算出 P 点到的中桩桩号和距离 中桩距离 K，根据中桩桩号和有关设计数据内插出边坡坡脚点的设计高程 Hs 和距离中桩的距离 K1， 根据 Hc、和 Hs 的高差计算出 K2，按照 P 点的高程，边坡顶距中桩的理论宽度 K

42、理论，具体计算如下： K理论=K1+h*n1+m+h*n2+m+（Hc-Hs-2h）*n3 则 dK=K- K理论 dK 求出来后，可以指挥司镜员靠近道路 dK 米，然后重复以上步骤，直到 dK 值能满足边坡控制精度。 K1 h h m m dK K2 程序应用程序应用 21 R2 R2 （二） 本程序在隧道断面欠挖、 超挖放样中的应用。 利用本程序进行隧道断面开挖的欠挖和超挖计算， 看下图： 用全站仪测出 P 点的坐标和标高后，根据坐标算出 P 点的中桩桩号和距中距离 K，根据中桩桩号查看 P 点的路面高程，算出 P 点至设计路面的高差 h，根据 P 点所在圆弧的圆心 O2 在断面中的位置（

43、K2，h2）算 出 P 点到 O2 点距离 R2，用 R2-R2 的差距就是 P 点的超挖或者欠挖的数据。 （三）本程序在路面施工控制放样中的应用 路面施工时，路面的轴线偏位、路面宽度、路面高程是测量的重要控制指标，也是道路施工控制最严格 的一道测量控制。利用本程序输入桩号和距中距离，可以直接计算出所求点至测站点的水平距离和方位角， 从而直接通过极坐标法直接放出该点在实地的位置。 （四）处理线路设计中的“断链”问题 断链是线路设计中经常出现的问题，断链分为长链与短链，如果是短链，那么只要弯道的计算范围从短 链处分界就可以了，如果是长链，则需要把线路分成两个文件，从开始出现断链处分，这样可以避免

44、桩号相 同，但位置不同的情况。 （五）本程序在道路竣工验收过程中的应用 道路竣工验收，其主要方面就是轴线偏位和宽度，可以通过现场测量道路中线的坐标，通过坐标反算出 该点和设计中线的偏移值 K（轴线偏位值）。 22 公路纵断面计算也是采用计算程序和数据程序分开的原则，数据文件是专门存放道路纵断面设计参数， 对于多个工程，可以建立多个数据文件，这样每个数据文件通过调用计算程序就可以计算出该点的高程。 数据文件：每个纵断面由纵坡和竖曲线组成，本程序每一个变坡点计算范围是从竖曲线前的纵坡任意桩 号开始到竖曲线后纵坡上的任一桩号结束。下表是某线路路基纵断面设计的参数一览表： 路基纵断面设计 桩号 变坡点

45、高程 半径 切线长 竖曲线起点桩号 竖曲线终点桩号 纵坡 K037+950.00 80.56 K038+650.00 75.93 3000.00 70.40 K038+579.599 K038+720.401 -0.6614% K038+900.00 86.01 4500.00 79.372 K038+820.628 K038+979.372 4.0320% K039+245.00 87.75 6000.00 59.597 K039+185.403 K039+304.597 0.5043% K039+410.00 91.86 6000.00 52.227 K039+357.773 K039+4

46、62.227 2.4909% K039+590.00 93.21 10000.00 52.214 K039+532.786 K039+647.214 0.7500% 99.84 1.8943% 根据上表编程数据文件如下：文件名 从上面程序看，Dh?E:Lbl 0:L=?L:I=?I:K?K 语句是固定的，“Dh”提醒使用者输入所求点距离路面 的厚度，比如路面底基层距离路面厚度为 37cm，则输入 0.37，表示所求点是路面底基层顶面高程， “L=”， 输入所求点桩号,“I”=表示输入路面横坡值，“K=”表示输入路面某点至中桩距离，均为正值。以上数据均 在程序运行时输入。在每一个变坡点处，要求输入的参数有： 序 变量 说明 1 C； 输入变坡点高程 2 D： 输入边坡点桩号 3 M： 输入变坡点小里程段纵坡坡度 4 N： 输入变坡点大里程侧纵坡坡度 5 R： 竖曲线半径 L 表示竖曲线前纵坡任意里程的桩号，本程序每段竖曲线参数计算范围是从竖曲线前的纵坡上任意点到竖曲 线后纵坡上任意点的桩号。 Dh?E:Lbl 0:L=?L:I=?I:K=?K If L37950:Then 75.93C:38650D:-0.006614M:0.04032N:3000R:IfEnd If L38820.628:Then 86.01C:38900D:0.040