卡西欧5800计算器编程入门程序集

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圆曲线中边桩坐标计算程序9完整缓和曲线中边桩坐标计算程序10；线元法匝道连续计算程序（万能程序）11；交点法曲线连续计算程序12竖曲线高程连续计算程序序参考书目作者参考书目名称1王中伟教授《道路施工放样综合程序》，王中伟教授在QQ（595077）空间把自己的程序电子版贡献给广大用户，在此特别感谢！2覃辉教授《公路与铁路施工测量程序》。3卡西欧《fx-5800P说明书》。4其他可能参考的书籍。5特别感谢特别感谢广西路桥的李虎（东华理工学院97届毕业），他无私奉献了他的公路高铁测量程序。6免责申明本书中程序仅供用户免费参考，卡西欧（中国）贸易有限公司不对用户使用本书中程序发生的任何问题负责。7联系方式用户对用本书疑问，请发邮件到：gongchengji@casio.com.cn。欢迎交流！8客服热线卡西欧的客服热线：400-700-6655。9程序传输本书的程序用户可以到全国省会城市的测绘店免费传输G版程序。10程序征集如果您有非常好的实用程序，希望广大用户能免费得到分享，请和本书编者联系:gongchengji@casio.com.cn。11欢迎翻印本书欢迎翻印；希望更多的用户能得到，欢迎流通。本书电子稿及其他的测量程序源代码请到卡西欧网站下载:http://www.casio.com.cn/客户支持﹥下载中心﹥计算器﹥工程资料下载，或者到百度文库中搜索《fx-5800P卡西欧易学易用程序集》下载。G1G2G3G4G5G6G7G8G9G10G11G11-1G12G12-1G10-1G10-2fx-5800PG2011-9-30...常用测量程序1极坐标放样程序（已知坐标计算边长、方位角计算程序），程序名1Deg：Fix3设置角度单位为十进制，3位固定小数显示2"XO="?A："YO="?B提示输入起点的坐标（A,B）3"XN="?C："YN="?D提示输入端点的坐标（C,D）4Pol((C–A),(D–B))调用极坐标 函 关于工期滞后的函关于工程严重滞后的函关于工程进度滞后的回复函关于征求同志党风廉政意见的函关于征求廉洁自律情况的复函 数计算距离和角度,自动保存到I和J5J<0J+360→J把小于0的角度换算成方位角6"JL="：I◢显示水平距离I7"FWJ="：JDMS◢显示方位角J2已知边长方位角计算坐标程序，程序名1Deg：Fix3设置角度单位为十进制，3位固定小数显示2"XO="?A："YO="?B提示输入起点的坐标（A,B）3"L="?L："C="?C提示输入所求点的距离L和方位角C4A+Lcos(C)→X计算所求点的X坐标5B+Lsin(C)→Y计算所求点的Y坐标6"X="：X◢显示所求点的X坐标7"Y="：Y◢显示所求点的Y坐标3建筑轴线偏移程序(或扩大基础程序)。程序名1Deg：Fix3设置角度单位为十进制，3位固定小数显示2?A:?B:?C:?D:?E提示输入起点的坐标（A,B）、起点方位角C、轴线偏移距D,E.3A+Dcos(C)+Ecos(C-90)→X计算所求点的X坐标4B+Dsin(C)+Esin(C-90)→Y计算所求点的Y坐标"X="：X◢显示所求点的X坐标6"Y="：Y◢显示所求点的Y坐标7（）（）易学易用程序集34测角前方交会坐标计算，程序名1Deg：Fix3设置角度单位为十进制，3位固定小数显示2?A:?B:?C:?D:?E:?F提示输入两个已知点的坐标（A,B）,（C,D），观测的两个角度E和F3(A(tanF)-1)+C(tanE)-1+(D-B))÷((tanE)-1+(tanF)-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坐标75测角后方交会坐标计算，程序名1Deg：Fix3：11→DimZ设置角度单位为十进制，3位固定小数显示2?A:?B:?C:?D:?E:?F:?G:?H:?K提示输入两个已知点的坐标（A,B）,（C,D），（E,F）观测的角度G,H,K3Pol((C–A),(D–B)):J<0J+360→J:J→Z[4]计算点（）的方位角4Pol((E–A),(F–B)):J<0J+360→J:J→Z[5]计算点（）的方位角5Z[4]-Z[5]→Z[1]方位角计算顶角6Z[1]<0Z[1]+360→Z[1]计算点（）的顶角并存贮到变量7Pol((E–C),(F–D)):J<0J+360→J:J→Z[4]计算点（）的方位角8Pol((A–C),(B–D)):J<0J+360→J:J→Z[5]计算点（）的方位角9Z[4]-Z[5]→Z[2]方位角计算顶角10Z[2]<0Z[2]+360→Z[2]计算点（）的顶角并存贮到变量11Pol((A–E),(B–F)):J<0J+360→J:J→Z[4]计算点（）的方位角12Pol((C–E),(D–F)):J<0J+360→J:J→Z[5]计算点（）的方位角13Z[4]-Z[5]→Z[3]方位角计算顶角14Z[3]<0Z[3]+360→Z[3]计算点（）的顶角并存贮到变量15tan(Z[1])tan(G)÷(tan(Z[1])-tan(G))→Z[7]计算点（）的参数16tan(Z[2])tan(H)÷(tan(Z[2])-tan(H))→Z[8]计算点（）的参数17tan(Z[3])tan(K)÷(tan(Z[3])-tan(K))→Z[9]计算点（）的参数18(Z[7]A+Z[8]C+Z[9]E)÷(Z[7]+Z[8]+Z[9])→Z[10]计算所求点的X坐标19(Z[7]B+Z[8]D+Z[9]F)÷(Z[7]+Z[8]+Z[9])→Z[11]计算所求点的Y坐标20"XP="：Z[10]◢显示所求点的X坐标21"YP="：Z[11]显示所求点的Y坐标（）（）（）4226坐标转换计算程序，程序名1Deg:Fix3设置角度单位为十进制，3位固定小数显示2"X0"?A:"Y0"?B输入施工坐标系原点在统一坐标系中的X、Y坐标（A,B）3"ANGLE"?E输入统一坐标系的X轴顺时针旋转至施工坐标系的X轴的角度值4Cls:"LEFT(1),RIGHT(ELSE)"?F判断施工坐标系是左手系还是右手系5"XPP"?C:"YPP"?D待换算平面点在施工坐标系中的X、Y坐标（C,D）6[[A][B]]→MatA向矩阵A赋值7[[cos(E),-sin(E)][sin(E),cos(E)]]→MatB向矩阵B赋值8F≠1-D→D根据施工坐标系的类型决定变量D的取值9[[C][D]]→MatC向矩阵C赋值10MatA+MatB×MatC→MatD◢计算转换后的坐标值，并将结果显示在矩阵D中11两平面坐标系转换计算示意图：（）（）（）57直线中边桩坐标计算程序,程序名1Deg：Fix3设置角度单位为十进制，3位固定小数显示2"X0="?X:"Y0="?Y:"K0="?E:"FWJ="?J输入直线起点坐标(X,Y)，起点桩号E，起点方位角J3Lbl1设置程序标签4?K:?L:?A提示输入所求点桩号K，边距L，和中线夹角A5X+(K-E)cos(J)+Lcos(J+A)→C计算所求边桩点的X坐标6"XP=":C◢显示所求边桩点的X坐标7Y+(K-E)Sin(J)+LSin(J+A)→D计算所求边桩点的Y坐标8"YP=":D◢显示所求边桩点的Y坐标9Goto1转向程序标签1继续执行程序10备注：该程序请从起点往终点计算。8圆曲线中边桩坐标计算程序,程序名1Deg：Fix3设置角度单位为十进制，3位固定小数显示2"X="?A:"Y="?B:"FWJ="?C:"KO="?K:"R="?R输入圆曲线起点坐标（A,B），起点桩号K，起点切线方位角C,半径R（左负右正）3Lbl1设置程序标签4"K="?L提示输入所求点桩号L，5"B="?E:"N="?N提示输入所求点边距E和以及夹角N（左负右正）6(((L-K)÷R)÷2)r→M计算圆心角的一半，并把弧度转换成角度值。r的输入[FUNCTION][5](ANGLE)[2](r)7Abs(2Rsin(M))→F计算所求点的弦长，并取绝对值。8A+Fcos(C+M)+Ecos(C+2M+N)→G计算所求边桩点的X坐标9B+Fsin(C+M)+Esin(C+2M+N)→H计算所求边桩点的Y坐标10"XP=":G◢显示所求边桩点的X坐标11"YP=":H◢显示所求边桩点的Y坐标12Goto1转向程序标签1循环执行程序13备注：该程序请从起点往终点计算。69完整缓和曲线中边桩坐标计算程序,程序名1Deg：Fix3设置角度单位为十进制，3位固定小数显示2"X0="?X:"Y0="?Y:"FWJ="?C:"KO="?K:"A="?A:"R="?R输入缓和曲线起点坐标(X,Y)，起点切线方位角C，起点桩号K，缓和曲线参数A,小半径R（左负右正）3Lbl1设置程序标签4"K="?L提示输入所求点桩号L，5"B="?E:"D="?D提示输入所求点边距E和以及夹角D（左负右正）6Abs(R)÷R→S计算偏转系数并存储到S7L-K→P计算所求点的曲线长并存储到P8P-P∧（3）÷(40R2)+P∧（5）÷(3456R∧（4）)-P∧（7）÷(599040R∧（6）)+P∧（9）(÷175472640R∧（8）)→M计算切线支距X坐标并存储到M9P2÷(6R)-P∧（4）÷(336R∧（3）)+P∧（6）÷(42240R∧（5）)-P∧（8）÷(9676800R∧（7）)+P∧（10）÷(3530096640R∧（9）)→N计算切线支距Y坐标并存储到N10M2+N2→Q计算弦长并存储到Q11cos-1（M÷Q)→T计算弦长对应的偏角值。并存储到T12(P2÷A2÷2)r→U计算切线夹角并存储到U13X+Qcos(C+ST)+Ecos(C+SU+D)→G计算所求边桩点的X坐标14Y+Qsin(C+ST)+Esin(C+SU+D)→H计算所求边桩点的Y坐标15"XP=":G◢显示所求边桩点的X坐标16"YP=":H◢显示所求边桩点的Y坐标17Goto1转向程序标签1循环执行程序18备注：该程序只适合完整缓和曲线。7fx-5800P程序使用流程：本程序数据和主程序是分开的，编程时将不同的工程数据存放到不同的数据文件里，如A匝道，文件名为A，将匝道A所有的曲线线元参数输入A文件里。运行时只要运行文件名A的程序就可以了，具体运行流程见下图：输出所计算位置与测站坐标X,Y距离和方位角DIS,FW输出所计算位置坐标X,Y输入距离所求点到所求桩号中线距离K，左－右+输出P点所在线路的桩号L和距离K,左－右+输入中桩切向右侧夹角M输入所求P点的坐标XP,YP输入计算桩号L输入计算区域附近桩号输入测站点坐标X0,Y0运行数据程序文件如ZX、B、C等公路高铁程序计算下一点输入桩号L计算下一点第四位小数点后为0.0001，如1100.0001第四位小数点不能是0.0001，如1100.0038（一）交点法数据文件编辑交点法编写数据文件必须是对称型的，即直线段→缓和曲线段→圆曲线段→缓和曲线段→直线段，（如果任意一端没有直线段，则把直线段长度看做是0），另外圆曲线两侧缓和曲线的旋转常数必须相等，并且和直线段连接处的半径必须是无穷大。交点法数据文件编写一般是根据设计图纸提供的平面曲线参数一览表提供的参数来编写，每个弯道包括：弯道起点方位角（C），交点X坐标（D），交点Y坐标（E），缓和曲线长度（F，当没有设缓和曲线时，F=0），交点转交（G，向左转弯，G为负值，向右转弯，G取正值），交点桩号（H），弯道圆曲线半径（R）。下图是一段市政道路设计参数数据。根据上图提供的数据，可以编辑成如下的数据文件：文件名（O表示变量字母O，其余的0均为零）在程式中，有两个条件转移语句即IfL>0：Then98°39°35.12°→C:4774.384→D:2415.861→E:140→F:31°17°23°→G:410.007→H:600→R:IfEndIfL>1060：Then129°56°58.19°→C:4206.421→D:3093.946→E:70→F:-33°50°48°→G:1285.437→H:600→R:IfEnd„„3→DimZ：“X0”?A：”Y0”?B：“Ln”?L：Abs(L)Int(Abs(1000L))÷1000→Z*3+：Lbl0：IfZ*3+≠0.0001：Then?L：Else“Xp”?X：”Yp”?Y：X→Z*1+：Y→Z*2+：IfEnd：Lbl1:IfL>0：Then98°39°35.12°→C：4774.384→D：2415.861→E：140→F：31°17°23°→G：410.007→H：600→R：IfEnd：IfL>1060Then129°56°58.19°→C：4206.421→D：3093.946→E：70→F：-33°50°48°→G：1285.437→H:600→R：IfEnd：Prog”XLJS”：IfZ[3+≠0.0001：ThenGoto0：Else(Z[1]－X)cos(O)＋(Z[2]－Y)sin(O)→N：L＋N→L：Int(Abs(1000L))÷1000+0.0001→L：－(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：IfAbs（N）≥0.001：ThenGoto1：Else“L=“：L◢“K=”：K◢IfEnd：Goto0:IfEnd数据文件的编写9如果还有其他弯道，可以继续后加。在这些存放设计参数的语句程序表达式是固定的。编辑曲线参数时，每个曲线参数放在一个IfL>***.***（两个弯道中间直线段上的任意桩号）:Then***°**°**.**°→C（弯道起点方位角）：****.***→D（交点X坐标）：****.****→E（交点Y坐标）：***.***→F（缓和曲线长度：**°**°**°→G（转角，向左转为负值，向右转为正值）：***.***→H（交点里程桩号）:***→R（圆曲线半径）：IfEnd条件式语句里。如果有多个弯道，一直按上述形式编下去。变量说明：序语句说明1L>***.***弯道参数的起点范围，可以是弯道起点，也可以是弯道前直线段里的任意点桩号。2***°**°**.**°→C弯道起点方位角3****.***→D、****.****→E交点X坐标、Y坐标4***.***→F缓和曲线长度，如果没设臵缓和曲线，则缓和曲线长度看做05**°**°**°→G转角，向左转为负值，向右转为正值6***.***→H交点里程桩号7***→R圆曲线半径【注意】运行程序时，第一步是输入测站坐标X0，Y0，然后输入近似桩号，主要是为了加快收敛，输入时注意当根据桩号和距中距离计算所求点坐标时，请不要把小数点第四位输入0.0001，如K21+369.0921，请把桩号输入21369.092即可，相反，当根据坐标反算桩号和距中距离时，请在输入时把第四位小数数位0.0001，如K21+200，输入时请输入21200.0001。10（二）线元法数据文件编辑线元法数据文件编辑适合于各种各样的线型组合，特别是在大型立交区或者是不规则的线型组合地方，如小区道路、重丘区山路等有特别明显的有点，本程序为了减少使用者设计参数的输入，特别把直线与缓和曲线以及圆曲线连接的地方只要输入缓和曲线参数就可以了，这样可以减少输入很多设计参数，同时节约大量的内存。看下图：该图是一个半喇叭型的互通立交，一共是5条匝道，我们以D匝道为例，来看看线元法数据文件的建立。113→DimZ：“X0”?A:“Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))÷→Z[3]:Lbl0:IfZ[3]≠0.0001:Then?L:Else“Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl1:IfL≥0:Then302°18°26°→W：494347.926→X：477989.873→Y：55→P:40→R：39.216→I:39.289→J：IfEnd:IfL≥78.504：Then494380.677→X：477971.271→Y：358°28°45.9°→W：40→R：-65→P:39.216→I:78.504→J：IfEnd:Prog””:Prog””:IfZ[3]≠0.0001:ThenGoto0:Else(Z[1]－X)cos(O)＋(Z[2]－Y)sin(O)→N:L＋N→L:-(Z[1]－X)sin(O)+(Z[2]－Y)cos（O）→K：IfAbs（N）≥0.001：ThenGoto1：Else“L=“:L◢“K=”:K◢IfEnd:Goto0:IfEnd下表是D匝道的线元表：D匝道线元数据 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 ：从上表中可以看出，D匝道由三个线元组成，第一段是缓和曲线段，根据上表判断，该缓和曲线属于第一段缓和曲线，其最小半径处桩号为J=K0+039.289，该处坐标X=494347.926,Y=477989.873，半径为R=40米，方位角W=302°18′26.0″，缓和曲线旋转常数P=55，在缓和曲线终点K0+039.289处，有一圆曲线与其相连接，长度I=39.216。第二段是圆曲线，第三段又是一段缓和曲线，根据上表判断，该曲线属于第二缓和曲线，曲线最小半径位于与园曲线连接处，桩号J=K0+78.504，坐标为X=494380.677，Y=477971.271,方位角W=358°28′45.9″，最小半径R=40，缓和曲线旋转常数P=-65（为了判断是第一段还是第二缓和曲线，我们通过P的正负来判断，第一缓和段P统一用正数表示，第二缓和段统一用负号表示。起点处相连的圆曲线长度I可以看做0，也可以看做是39.216.由此我们整理下。从整个线型全部是向右转向，所以R全部取正数。具体如下：第一段缓和曲线线元参数说明：第二段缓和曲线线元参数最小半径处坐标X=494347.926Y=477989.873最小半径R=40最小半径处方位角：W=302°18′26.0″缓和曲线旋转常数P=55与之相连的圆曲线长度I=39.216最小半径处桩号J=39.289最小半径处坐标X=494380.677Y=477971.271最小半径R=40最小半径处方位角：W=358°28′45.9″缓和曲线旋转常数P=65与之相连的圆曲线长度I=39.216最小半径处桩号J=W=358°28′45.9″此段数据计算范围可以从K0+000~K0+78.504此段数据计算范围可以从K0+78.504到匝道结束。当桩号L位于K0+000~K0+78.504时，曲线参数按第一段编写，当L>78.504时，按第二段缓和曲线编写，当然第二段缓和曲线计算范围也可以从圆曲线起点K0+039.289到K0+78.504中间任意桩号开始。下面是D匝道的数据程序文件（文件名“”）：（O表示变量字母O，其余的0均为零）线元名称起点桩号偏向半径或缓和曲线旋转常数A线元长度(m)线元起点坐标、方位角X(m)Y(m)方位角缓和曲线DK0+000.00右偏5539.289494341.902478027.825260°38′55.7″圆曲线DK0+039.289右偏4039.216494347.926477989.873302°18′26.0″缓和曲线DK0+078.504右偏65105.625494380.677477971.271358°28′45.9″12在上面程式中，每个缓和曲线的线元参数放在一个If„Then„„IfEnd循环语句中，If后面紧跟计算的起点范围，终点范围从下段的起点范围结束。变量定义：序变量说明1X、Y输入缓和曲线半径最小处的X、Y坐标2R缓和曲线最小半径，当线路左转时，半径数据前面加“-”号，以判断线路转向。3W输入缓和曲线半径最小处的方位角4P缓和曲线旋转常数，第一缓和曲线长度为正数，第二缓和曲线常数为负数。当不设缓和曲线时，P=05I与缓和曲线最小半径端连接的圆曲线长度。6J缓和曲线小半径段的桩号。再看看E匝道数据文件编辑E匝道线元设计数据表线元名称起点桩号偏向半径或缓和曲线旋转常数A线元长度(m)线元起点坐标、方位角X(m)Y(m)方位角直线EK0+000.0071.296494360.619477706.29379°49′57.0″缓和曲线EK0+071.296右偏12072494373.205477776.47079°49′57.0″圆曲线EK0+143.296右偏20060.070494381.631477847.87090°08′44.6″缓和曲线EK0+203.366右偏9545.125494372.527477907.018107°21′16.3″缓和曲线EK0+248.491左偏7564.841494355.878477948.932113°49′05.6″圆曲线EK0+313.333左偏86.7517.800494337.373478010.65892°24′19.0″133→DimZ：“X0”?A:”Y0”?B:“Ln”?L:Abs(L)-Int(Abs(1000L))/1000→Z[3]:Lbl0:IfZ[3]≠0.0001:Then?L:Else“Xp”?X:”Yp”?Y:X→Z[1]：Y→Z[2]：IfEnd:Lbl1:IfL≥0:Then494381.6316→X：477847.870→Y：90°08°44.6°→W：200→R：120→P:60.07→I:143.296→J：IfEnd:IfL≥203.366：Then494372.527→X：477907.018→Y：107°21°16.3°→W：200→R：-95→P:60.07→I:203.366→J：IfEnd:IfL≥248.491：Then494337.373→X：478010.658→Y：92°24°19°→W：－86.75→R：75→P:17.8→I:313.333→J：IfEnd:Prog””:Prog””：IfZ[3]≠0.0001:ThenGoto0:Else(Z[1]－X)cos(Ｏ)＋(Z[2]－Y)sin(Ｏ)→N:L＋N→L:-(Z[1]－X)sin(Ｏ)+(Z[2]－Y)cos（Ｏ）→K：IfAbs（N）≥0.001：ThenGoto1：Else“L=“:L◢”K=”:K◢IfEnd:Goto0:IfEnd在匝道线元表中，我们可以看出，E匝道式由一段直线段→缓和曲线段→圆曲线段→缓和曲线段→缓和曲线段→圆曲线段连接起来的。值得注意的是，在位于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匝道数据线元名称起点桩号偏向半径或缓和曲线旋转常数A线元长度(m)线元起点坐标、方位角X(m)Y(m)方位角圆曲线CK0+000.00右偏3256.876495005.991478056.67756°19′42.3″直线CK0+056.87621.277494991.262478104.123158°09′54.4″缓和曲线CK0+078.153右偏6530.179494971.511478112.036158°09′54.4″圆曲线CK0+108.332右偏14030.431494943.127478122.242164°20′25.7″缓和曲线CK0+138.763右偏6530.179494913.167478127.219176°47′40.6″直线CK0+168.941右偏85422.023494883.008478126.740182.°58′11.8″缓和曲线CK0+590.964右偏8585494461.552478104.874182.°58′11.8″圆曲线CK0+675.964右偏8572.739494379.485478086.681211°37′04.2″直线FK0+000.00111.525494591.311478420.033252°32′00.2″缓和曲线FK0+111.505右偏10090.909494557.835478313.650252°32′00.2″圆曲线FK0+202.434右偏110101.114494542.812478224.689276°12′33.6″缓和曲线FK0+303.548右偏120130.909494595.310478142.42032852′35.7″缓和曲线GK0+000.000右偏175133.152494718.802478112.963182°58′11.8″圆曲线GK0+133.152右偏230149.124494587.600478093.368199°33′17.5″缓和曲线GK0+282.286右偏190149.72494472.335478002.890236°42′21.1″15计算程序是核心程序，但是道路计算程序是固定的，使用者不需要做任何修改，也不直接运行道路计算程序分为两部分，第一部分是常规的对成型线型组合计算的程序（文件名：），交点法编辑的数据文件直接调用；第二部分是适用于非对称的线型组合计算文件（文件名：），在数据文件调用完“”文件后，还要调用“”文件。程序代码：（为了输入方便代码分行显示，使用者不必分行输入。）Deg:Fix4G÷Abs(G)→V:F÷2-F^(3)÷（240R2）+F^（5）÷(34560R^（4）)-F^（7）÷(8386560R^（6）)+F^（9）÷(3158507520R^（8）)→Q:F2÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:Q+(R+P)Vtan(G÷2)→T:πRGV÷180+F→I:IFL≤H-T:ThenL-H→J:Goto1:IfEnd:IfL<H-T+F:ThenL-H+T→J:C→Z:1→S:Goto2:IfEnd:IfL<H-T+I-F:ThenL-H+T→J:C→Z:1→S:Goto3:IfEnd:IfL<H-T+I:ThenH-T+I-L→J:C+G→C:-1→S:Goto2:IfEnd:IfL≥H-T+I:ThenL-H-I+2T→J:C+G→C:Goto1:IfEnd:Lbl1:D+Jcos(C)→X:E+Jsin(C)→Y:C→O:GotoC:Lbl2:(J-J^(5)÷(40(RF)2)+J^(9)÷(3456(RF)^(4))-J^(13)÷(599040(RF)^(6))+J^(17)÷(175472640(RF)^(8))-T)S→M:(J^(3)÷(6RF)-J^(7)÷(336(RF)^(3))+J^(11)÷(42240(RF)^(5))-J^(15)÷(9676800(RF)^(7))+J^(19)÷(3530097000(RF)^(9)))V→N:C+90J2SV÷(πRF)→O:GotoBLbl3:180(J-0.5F)÷(πR)→O:Q+Rsin(O)-T→M:(P+R(1-cos(O)))V→N:C+VO→O：GotoBLblB:D+Mcos(C)-Nsin(C)→X:E+Msin(C)+Ncos(C)→Y:GotoC:LblC:IfZ*3+≠0.0001:Then?K:"M"?U:K<0=>180-U→U:O+U×Abs(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:IfX-A<0:ThenW+180→W:IfEnd:IfW<0:ThenW+360→W:IfEnd:√((X-A)2+(Y-B)2)→S:"DIS":S◢"FW":WDMS◢IfEnd道路计算程序16非对称线型组合计算程序（）（O表示变量字母O，其余的0均为零）最小半径处坐标X=494347.926Y=477989.873最小半径R=40最小半径处方位角：W=302°18′26.0″缓和曲线旋转常数P=55与之相连的圆曲线长度I=39.216最小半径处桩号J=39.289Abs(P+10^(-10)）÷(P+10^(-10))→S：Abs(R)÷R→V:Abs(R)→R:P2÷R→F:180(I+F)V÷(πR)→G:90FSV÷(πR)→O:IfS=1:ThenW-O→C:IfEnd:IfS=-1:ThenW-O-G→C:IfEnd:C<0C+360→C:C≥360C-360→C:F2÷(24R)-F^(4)÷(2688R^(3))+F^(6)÷(506880R^(5))-F^(8)÷(154828800R^(7))→P:F÷2-F^(3)÷(240R2)+F^5÷(34560R^(4))-F^(7)÷(8386560R^6)+F^9÷(3158507520R^(8))→Q:Q+(R+P)Vtan(G÷2)→T:IfS=1:ThenJ-F+T→H:C→Z:ElseJ-I-F+T→H:C+G→Z:IfEnd:(F-F^(3)÷(40R2)+F^(5)÷(3456R^(4))-F^(7)÷(599040R^(6))+F^(9)÷(175472640R^(8))-T)S→M:(F2÷(6R)-F^(4)÷(336R^(3))+F^(6)÷(42240R^(5))-F^(8)÷(9676800R^(7))+F^(10)÷(3530097000R^(9)))V→N:X-Mcos(Z)+Nsin(Z)→D:Y-Msin(Z)-Ncos(Z)→E17程序运行分为两种，第一种是根据根据使用者输入的桩号和距中距离计算出所求点的坐标和测站点到该点的方位角和水平距离，可以直接利用计算器提供的距离和方位角进行极坐标放样。另一种是根据使用者输入的实测坐标计算出该点的桩号和距中距离，可以根据桩号和距中距离进行高边坡、隧道断面点的设计高程、以及进行路面高程控制计算方面很方便。（一）据桩号和距中距离计算所求点的坐标、及方位角和水平距离。要进行程序的线路计算，直接运行数程序文件就可以了，不需要理会计算程序。如果计算某桩号L距离中线K处的坐标X，Y和该点到测站点的距离和方位角。运行提示如下：序变量说明1X0?提示输入测站点X坐标2Y0？提示输入测站点Y坐标（测站点的位数和数据程序文件里坐标取位要一致）3Ln?计算点近似桩号，近似桩号输入第四位小数后不能等于0.0001，4L？输入计算点准确桩号，近似桩号输入第四位小数后不能等于0.0001。5K？提示输入距离中桩的桩号如果计算点位于线路的左幅，则在数据前加入“-”号，如果是计算中桩坐标，则输入06M？输入该点与线路的前进方向的按顺时针转过的角度，常于斜交的涵洞，桥梁等构造物7X显示所求点的X坐标8Y显示所求点的Y坐标9FW显示所求点至测站点的方位角10DIS显示所求点至测站点的距离要继续计算则从L?开始，计算新点的坐标，具体流程见程序流程图的左边部分。例：要计算D匝道的DK0+020中桩坐标、DK0+080左边5.3米的坐标、Dk0+240斜交75°58′11.8″、距中右侧5米处的坐标，并计算出该点至坐标为X=495005.991，Y=478056.677控制点的水平距离和方位角。程序运行18运行文件“”提示如下：序变量说明1X0?495005.991测站点X坐标2Y0？478056.677测站点Y坐标3Ln?0输入小输点第四位起不等于0.0001的桩号4L？20输入所求点桩号5K？0输入距中距离6M？90输入斜交角度7X494341.482显示所求点X坐标8Y478007.903显示所求点Y坐标9DIS666.2959显示所求点至测站平距10FW184°11′52.42″显示所求点至测站方位角11L？80输入下点的桩号12K？-5.313M？9014X=494382.229015Y=477965.959416DIS630.32417FW188°16′29.5″18L？24019K？5.020M？75°58°11.8°21X=494455.311022Y=478101.549123DIS552.50524FW175°20′29.6″19（二）据坐标反算桩号和距中，根据某点的坐标，计算出该点在线路中的桩号和距离，程序运行时提示如下：序变量说明1X0?提示输入测站点X坐标2Y0？提示输入测站点Y坐标（测站点的位数和数据程序文件里坐标取位要一致）3Ln?计算点近似桩号，近似桩号输入第四位小数后不能等于0.0001，4Xp？所求点P的X坐标5Yp？所求点P的Y坐标6L显示所求点的桩号7K显示所求点的距中距离，如果位于线路的左边，则显示结果为负值，如果K为0时，则该点在线路中线上。计算新点的桩号和距中距离，则从?Xp开始。如此循环计算。具体流程见流程图的右边部分例：求点X=494382.2290，Y=477965.9594在D匝道中的位置，运行文件“”如下：序变量说明1X0?495005.991测站X坐标(可输入0)2Y0？478056.677测站Y坐标(可输入0)3Ln?60.0001输入小输点第四位起等于0.0001的桩号4Xp？494382.2290输入所求点X坐标5Yp？477965.9594输入所求点Y坐标6L79.999999显示所求点桩号7K-5.3000显示所求点距中距离,“-“表示该点位于中桩左侧20（一）本程序在高边坡放样中的应用高边坡放样是目前公路道路测量工作普遍要做的一项工作，边坡坡顶线或者填方边坡坡脚线是边坡放样最繁琐的工作，利用本程序进行边坡放样，工作步骤如下：①司镜员走到目的地后，先测出一个点的坐标和高程；②利用本程序求出该点在线路中的桩号和距中距离，然后根据计算出来的桩号和设计图纸，内插出道路边桩高程，根据道路边桩高程和道路设计宽度、坡度、设计碎落平台的宽度和间隔高度、实测点的高程计算出按实测高程，边坡距中距离。③根据高程算出的距中距离和程序计算出立镜点至距中距离差值dK，指挥司镜员靠近道路中线或者往外偏移多少米。④重新做1~3工作，直到dK小于边坡宽度控制精度要求为止。如上图，司镜员在P点立镜，测出P点坐标和高程后，用本程序计算出P点到的中桩桩号和距离中桩距离K，根据中桩桩号和有关设计数据内插出边坡坡脚点的设计高程Hs和距离中桩的距离K1，根据Hc、和Hs的高差计算出K2，按照P点的高程，边坡顶距中桩的理论宽度K理论，具体计算如下：K理论=K1+h*n1+m+h*n2+m+（Hc-Hs-2h）*n3则dK=K-K理论dK求出来后，可以指挥司镜员靠近道路dK米，然后重复以上步骤，直到dK值能满足边坡控制精度。K1hhmmdKK2程序应用21R2R2'（二）本程序在隧道断面欠挖、超挖放样中的应用。利用本程序进行隧道断面开挖的欠挖和超挖计算，看下图：用全站仪测出P点的坐标和标高后，根据坐标算出P点的中桩桩号和距中距离K，根据中桩桩号查看P点的路面高程，算出P点至设计路面的高差h，根据P点所在圆弧的圆心O2在断面中的位置（K2，h2）算出P点到O2点距离R2’，用R2’-R2的差距就是P点的超挖或者欠挖的数据。（三）本程序在路面施工控制放样中的应用路面施工时，路面的轴线偏位、路面宽度、路面高程是测量的重要控制指标，也是道路施工控制最严格的一道测量控制。利用本程序输入桩号和距中距离，可以直接计算出所求点至测站点的水平距离和方位角，从而直接通过极坐标法直接放出该点在实地的位置。（四）处理线路设计中的“断链”问题断链是线路设计中经常出现的问题，断链分为长链与短链，如果是短链，那么只要弯道的计算范围从短链处分界就可以了，如果是长链，则需要把线路分成两个文件，从开始出现断链处分，这样可以避免桩号相同，但位置不同的情况。（五）本程序在道路竣工验收过程中的应用道路竣工验收，其主要方面就是轴线偏位和宽度，可以通过现场测量道路中线的坐标，通过坐标反算出该点和设计中线的偏移值K（轴线偏位值）。22公路纵断面计算也是采用计算程序和数据程序分开的原则，数据文件是专门存放道路纵断面设计参数，对于多个工程，可以建立多个数据文件，这样每个数据文件通过调用计算程序就可以计算出该点的高程。①数据文件：每个纵断面由纵坡和竖曲线组成，本程序每一个变坡点计算范围是从竖曲线前的纵坡任意桩号开始到竖曲线后纵坡上的任一桩号结束。下表是某线路路基纵断面设计的参数一览表：路基纵断面设计桩号变坡点高程半径切线长竖曲线起点桩号竖曲线终点桩号纵坡K037+950.0080.56K038+650.0075.933000.0070.40K038+579.599K038+720.401-0.6614%K038+900.0086.014500.0079.372K038+820.628K038+979.3724.0320%K039+245.0087.756000.0059.597K039+185.403K039+304.5970.5043%K039+410.0091.866000.0052.227K039+357.773K039+462.2272.4909%K039+590.0093.2110000.0052.214K039+532.786K039+647.2140.7500%99.841.8943%根据上表编程数据文件如下：文件名从上面程序看，"Dh"?E:Lbl0:"L="?L:"I="?I:"K"?K语句是固定的，“Dh”提醒使用者输入所求点距离路面的厚度，比如路面底基层距离路面厚度为37cm，则输入0.37，表示所求点是路面底基层顶面高程，“L=”，输入所求点桩号,“I”=表示输入路面横坡值，“K=”表示输入路面某点至中桩距离，均为正值。以上数据均在程序运行时输入。在每一个变坡点处，要求输入的参数有：序变量说明1C；输入变坡点高程2D：输入边坡点桩号3M：输入变坡点小里程段纵坡坡度4N：输入变坡点大里程侧纵坡坡度5R：竖曲线半径L表示竖曲线前纵坡任意里程的桩号，本程序每段竖曲线参数计算