道路中边桩坐标放样正反算CASIO fx-5800P程序(全线贯通)

CASIO fx-4800P、fx-5800P型计算器用于线路施工曲线中线点坐标的计算程序中铁十局三建公司工程技术部摘要：本文介绍了CASIO fx-4800P 、fx-5800P型计算器程序编制用于铁路、公路曲线线路内任意中线点的坐标计算程序及使用方法。

本计算程序具有操作简便、计算快捷、应用广泛等特点、极大地减轻了测量工作者的内业工作量，对于测量工作者有较大的参考和指导作用。

关键词：曲线线路施工测量计算程序1．概述过去，线路中线施工放样基本依靠经纬仪和钢尺了来进行角度及距离测量。

对于曲线线路一般的测量方法是：经纬仪置于某一中线点上，采用偏角法拨角再用钢尺量距来定出中线点。

随着电子技术进步和经济发展，测量仪器和测量方法的不断改进，目前，全站仪已广泛地应用于工程施工测量中，极大的提高了测量工作效率。

但是，在进行铁路、公路工程的曲线线路施工测设时，需要在线路所在区域建立统一坐标系或独立坐标系，利用坐标变换的方法，将整个曲线的三个部分（第一缓和曲线、中间圆曲线、第二缓和曲线）统一到同一坐标系中。

根据坐标系的建立，计算出整个曲线内任意点的坐标，再采用全站仪利用极坐标方法进行施工放样。

前提是首先利用计算器计算出各中线点坐标，然后才能进行放样。

而普通型计算器不仅计算速度慢，且要求计算者必须正确地记忆很多计算公式，计算繁琐而且容易出错，满足不了现场测设工作的要求。

为了能够快速准确地为全站仪提供测设数据，发挥全站仪快速测设的特点，提高测量工作效率，应采用可编程的计算器，编制计算程序。

本文主要介绍应用CASIO fx-4800P型计算器的计算程序，供公司测量同行们参照使用。

2．计算程序QXZBJS（文件名：曲线坐标计算fx-4800P）Defm2:R：L：A：N“ZH：X=”：E“ZH：Y=”：F：“FWJ=”：K“ZH：LC=”： P=L2/(24R)-L4/(2688R3):M=L/2-L3/(240R2):T“T”=（R+P）tng(A/2)+M ◢G=RAπ/180：“S”S=G+L◢LbiA：｛C，V｝：C“CSDLC=”：V“HXPJ=”：D=C-K：D≤L＝＞I=D-D5/（40R2L2）：U=D3/(6RL)-D7/（336R3L3):J=√（I2+U2）：Goto1：≠＞D≤G＝＞O=90（2D-L）/( Rπ)：I=RsinO+M：U=R（1-cosO）+P：J=√（I2+U2）：Goto2：≠＞D=S-（C-K）： = D-D5/（40R2L2）：Z[2]=D3/6RL-D7/（336R3L3）：I=T+（T-Z[1]）cosA-Z[2]sinA：U=（T-Z[1]）sinA+Z[2]cosA：J=√（I2+U2）：Goto3：Lbi1:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-30D2/ （RLπ）：H=F-90D2/ （RLπ）：≠＞Q=F+30D2/ （RLπ）：H=F+90D2/（ RLπ）⊿ Goto4：Lbi2:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-tng-1(U/I)：H=F-O：≠＞Q= F+tng-1(U/I)：H=F+O⊿Goto4：Lbi3:{Q}:Q“Z=1；Y=2”：Q=1＝＞Q=F-tng-1(U/I):H=F-(A-90(S-(C-K))2/ （RLπ）): ≠＞Q= F+tng-1(U/I):H=F+(A-90(S-(C-K))2/ （RLπ）):⊿ Goto4：Lbi4:B=90+H:H＜0＝＞H“QXFWJ”=B+360◢≠＞H≥360＝＞H“QXFWJ”=H-360 ◢≠＞H“QXFWJ”=H◢⊿Goto5：Lbi5: X“CSD:X”=JcosQ+N+VcosB◢ Y“CSD:Y”=JsinQ+E+VsinB◢GotoA3．程序说明3.1 输入已知变量R—圆曲线半径，显示R？L—缓和曲线长，显示L？A—曲线转向角，显示A？E—直缓点纵坐标，显示ZH：X=？N—直缓点横坐标，显示ZH：Y=？F—第一切线方位角，即ZH至JD的方位角，显示FWJ=？K—直缓点里程，显示ZH：LC=？3.2 计算待求量T—切线长度，显示T= …S—曲线全长，显示S= …Z[3]—外矢距，即JD到QZ的距离，显示E0= …3.3 输入待求变量K—输入待求（测设）点的里程，显示LC=？V—横向偏距，即测设点左、右侧外移距的偏移量，若为中线点输入0；右侧输入“+”值，左侧输入“-”值。

5800计算器全线坐标计算放样程序（修改版）“XLZBJSCX” ◢ （第一个程序）Lb1 0 ↙Cls : Fix 4 : 30→Dimz ↙“XHS="?G ( 后视点 X) :"YHS="?L ( 后视点 Y) :"XZJ="?M ( 置镜点 X) :"YZJ="?N ( 置镜点 Y) :Pol(G-M,L-N):"DH=":I ( 后视距) ◢ J<0=>J+360→J:"FH=":J►DMS ◢ ( 后视方位角 ) Lbl 1 ↙ ( If （如果的意思） And （和字的意思） Then （然后的意思） )“K=”?K ◢ （计算里程） ( 下面有色的是数据库 )If （如果） K< 51760.052 本曲线缓直点桩号 And （和）K≥ 51048.785 上一个曲线缓直点桩号： Then （然后）本曲线缓直点桩号51760.052 →Z[1] : 上一个曲线缓直点桩号51048.785 →Z[2] ： 1 -1 →O（注：左偏曲线输入 - 1→O, 右偏曲线输入1→O） : 偏角12 ’ 23 ’ 19.5 ’ →A ：半径3289.486 →R : 第一缓和曲线长度0 →Z[6] : 第二缓和曲线长度0 →Z[7] : 交点 X 坐标→B : 交点 Y 坐标→C : 小里程向交点方位角→E : 交点向大里程方位角→F : Goto 2 : IfEnd ↙ ( 重兰字是输入的数字 ) ………… （曲线段分段输入）补充直线段输入如下If （如果） K< 本段直线终点里程 And （和）K≥ 本段直线起点里程 :Then （然后）1→O: 本段直线终点里程→Z[3]: 终点坐标X→Z[16]: 终点坐标Y→Z[17]: 方位角→E:Goto 4:IfEnd ↙Lb1 2 ↙ （曲线要素计算）Z[6] ÷ 2- Z[6]^ 3 ÷ (240R^2)+ Z[6]^ 5 ÷ (34560*R^4) →Z[8] ↙ （ M1 不输）Z[7 ] ÷ 2- Z[7]^ 3 ÷ (240R^2)+ Z[7]^ 5 ÷ (34560 R^4) →Z[9] ↙ （ M2 不输）Z[6]^ 2 ÷ (24R)- Z[6]^ 4 ÷ (2688R^3) →Z[10] ↙ （ P1 不输）Z[7]^ 2 ÷ (24R)- Z[7]^ 4 ÷ (2688R^3) →Z[11] ↙ （ P2 不输）(πAR ) ÷ 180+0. 5 × ( Z[6]+ Z[7])→S ↙ （曲线总长）90 × Z[6 ] ÷ ( R × π) →Z[14] ↙ （第一缓和曲线总偏角）( ×÷以后自己改 )9 0 × Z[7 ] ÷ ( R × π) →Z[15] ↙ （第二缓和曲线总偏角 , 可以省略）Z[8] ＋（ R+Z[10])TAN(A/2)-(Z[10]-Z[11] )/SIN A→Z[12] ↙ ( 切线 T1)Z[9] ＋（ R+Z[11])TAN(A/2)+(Z[10]-Z[11] )/SIN A→Z[13] ↙ ( 切线 T2)B+ Z[12]*COS (E+180)→ Z[16] ↙ （ ZH 点 X ）C+ Z[12]*SIN(E+180)→ Z[17] ↙ （ ZH 点 Y ）Z[1]-S→Z[3] ↙ (ZH 点里程 )Z[3]+ Z[6]→Z[4] ↙ (HY 点里程 )Z[1]- Z[7]→Z[5] ↙ (YH 点里程 )GOTO 3 ↙LB1 3 ↙ ( 判断里程点与曲线关系 )If K≤Z[3] And K> Z[2] : Then Goto 4 : IfEnd ↙If K≤Z[4] A nd K> Z[3] : T hen Goto 5 : IfEnd ↙If K≤Z[5] A nd K> Z[4] : T hen Goto 6 : IfEnd ↙If K≤Z[1] A nd K> Z[5] : T hen Goto 7 : IfEnd ↙LB1 4 ↙ （里程小于直缓点直线独立坐标）K- Z[3] →X : 0→Y : E→T : Prog“TYZBCX” ：Goto 1 ↙Lb1 5 ↙ （第一缓和曲线独立坐标）K- Z[3] →H ↙H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →X ↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Y ↙90*H^2/( R*π* Z[6]) →T ↙IF O >0 ：Then T +E→T : Else E-T →T : T<0=>360+T→T : IfEnd ↙PRrog“TYZBCX” ：Goto 1 ↙Lb1 6 ↙ （圆曲线独立坐标）K- Z[4] →H ↙H*180/( R*π)+ Z[14]→T ↙R*SIN( T)+ Z[8]→X ↙R*(1-COS (T))+ Z[10]→Y ↙IF O >0 ：Then T +E→T : Else E-T →T : T<0=>360+T→T : IfEnd ↙Prog“TYZBCX” ：Goto 1 ↙Lb1 7 ↙ （第二缓和曲线独立坐标）Z[1] -K →H↙H-H^5/(40*R^2* Z[7]^2)+H^9/(3456*R^4* Z[7 ]^4) →U ↙H^3/(6*R* Z[7])-H^7/(336*R^3* Z[7]^3) →V ↙90*H^2/( R*π* Z[7]) →T ↙Z[13]COS (A)+ Z[12]-U*COS( A)-V*SIN (A)→X ↙Z[13]*SIN( A)-U*SIN( A)+V*COS (A)→Y ↙IF O >0 ： Then F-T→T : T<0=>360+T→T : Else F+T →T : IfEnd ↙Prog“TYZBCX” ：Goto 1 ↙子程序：“TYZBCX” ↙ （统一坐标计算）（第 2 个程序）IF O<0 : Then -Y→Y : IfEnd ↙“QXJ=” :T ◢ （计算里程点切线方位角，可以不显示）Z[16]+X*COS (E)-Y*SIN( E)→Z[18] ↙Z[17]+X*SIN (E ） +Y*COS (E ）→Z[19] ↙“XI=” : Z[18] ◢ （ XI 中线 X ）“ YI =” : Z[19] ◢ （ YI 中线 Y ）Pol(Z[18]-M,Z[19]-N):"DI=":I ◢ （中桩放样距）J<0=>J+360→J:"FI=": J ►DMS ◢ （中桩放样方位角）“ PJ =”?P ◢ ( 输入边桩与线路夹角 PJ ，左偏– 90 右 +90 )“ PD =”?D ◢ （输入边桩距 PD ）Z[18]+D*COS(T+P) →Z[20] ↙Z[19]+D*SIN(T+P) →Z[21] ↙“XP=”: Z[20] ◢ （ XP 边桩 X ）“YP=”: Z[21] ◢ （ YP 边桩 Y ）Pol(Z[20]-M,Z[21]-N):"DP=":I ◢ （边桩放样距）J<0=>J+360→J:"FP=":J ►DMS ◢ （边桩放样方位角）Return ↙ （以上都要输进计算器）注解不输K 里程 XI 中线 X YI 中线 Y PD 输入边桩距 XP 边桩 X 坐标YP 边桩 Y 坐标 I 边桩放样距 PJ 输入边桩与线路夹角，左偏– 90 右 +90 ) I 中桩放样距 T 计算里程点切线方位角，可以不显示卡西欧FX5800全线贯通万能正、反算程序FX5800计算器的积分程序（正反算、全线贯通、新线路）终极版ZHUCHENGXU 主程序"1.ZS,2.FS" ?→Q输入1正算，输入2反算“NEW=0,OLD≠0”?ZIf Z=0:Then “X0=”?A:“Y0=”?B:“C0=”?C:“1/R0=”?D:“1/RI=”?E:“SP=”?F:“EP=”?G:Ifend:Q=2=>Goto 2Lbl 1 :“KM=,<0Stop”?H:H<0=>Stop:“PJ=”?O:“PY=”?LLbl Z:Z=1=> Prog“01”：Z=2=> Prog“02”选择数据库文件，可增加H- F→X:0.5(E-D)÷(G-F)→NC+(XD+NX2)*180÷π→P：P<0=>P+360→P：P>360=>P-360→PA+∫(cos(C+(XD+NX2)*180÷π),0,X)+Lcos(P+O)→UB+∫(sin(C+(XD+NX2)*180÷π),0,X)+Lsin(P+O)→VQ=2=>Goto 4：Cls：Fix 3"Xn=":Locate 4，1，U："Yn=": Locate 5，2，V：“FWJ=”：P▶DMS◢Norm 2：Cls：Goto 1Lbl 2:“XD=,<0,STOP”？R：R<0=>Stop:“YD=”？S“KMDG=”？H ：90→O:0→L：Goto Z （H线路范围内的任意桩号）Lbl 4:Pol（R-U，S-V）：J<0 => J+360→JWhile abs（Icos（J-P））≤0.001:P-J>180=> J+360→J： P-J<-180=> P+360→P:IF P-J>0:then -I→L:else I→L ifendGoto 3: Whileend：H+Icos（J-P）→H:Goto ZLbl 3:Cls：Fix 3“KM=”: Locate4，1，H：“PY=”: Locate4，2，L◢Norm 2：Cls：Goto 201（数据库子程序）If H<=第一曲线终点桩号：then 第一曲线起点X→A：第一曲线起点Y→B：第一曲线起点方位角→C：起点曲率→D：终点曲率→E：起点桩号→F：终点桩号→G：return：ifend……………程序说明：1、该程序可以计算任意线形（直线、圆曲线、缓和曲线、不完整曲线）任意桩号的坐标（正算，输入1），也可根据坐标计算该点到线路的距离及垂足桩号（反算，输入2）；2、（NEW=0，OLD≠0）？如果要计算的点为数据库线路中的点，则输入数据库编号（以整数1、2、3…代替输入）；如果在数据库中没有要计算线路的数据，则输入曲线要素X0：曲线起点X坐标；Y0：曲线起点Y坐标；C0:曲线起点方位角；R0-1、 RI-1:曲线起点、终点曲率，直线为0，曲线左偏输入负值，右偏输入正值；SP、 EP：曲线起点桩号，终点桩号；KM：待求点桩号；PJ：正斜交的设定；PY：偏中距离，线路上的点输入0，右偏输入+值，左偏输入-值；3、正算显示坐标及切线方位角；反算输入线路的任意桩号（此桩号越接近真实值计算速度越快）、待求点坐标，显示待求点桩号及偏中距离；4、正算子程序为积分公式编写而成；反算子程序为角度趋近的方法编写，计算速度有点慢。

道路中边桩坐标放样正反算CASIO fx-5800P程序（全线贯通）
程序的特点：
1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写；
2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐；3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。

.二、程序代码
.
.
.
.
.
.
注：路线数据库子程序ROAD－DATA1是根据计算示例项目的直曲表编写，大家使用时应按各自项目的直曲表改写或新建。

.
.
四、程序变量清单
.
.
五、计算流程示例
1．中桩坐标计算示例
计算任务：计算高速公路K6+100～K6+700段的中桩坐标及切线方位角（桩距20m），并在导线点（2807118.026，474113.687）上架设全站仪，计算各中桩的极坐标放样数据。

使用道路中边桩坐标放样正反算程序RAOD－2的操作流程见下表。

.
.
.
.
2．边桩坐标计算示例
计算任务：计算HY高速公路主线K6+100～K6+700段的边桩（左右各12.25米）坐标（桩距20m）。

使用道路中边桩坐标放样正反算程序RAOD－2的操作流程见下表。

.
.
3．坐标反算计算示例
计算任务：根据前面计算的K6+100的中桩、左右边桩的坐标，反算对应桩号及偏距，并进行验证。

使用道路中边桩坐标放样正反算程序RAOD－2的操作流程见下表。

.。

本程序非本人原创，结合上传者长处，进行改编，愿收获和大家共享，不足之处多多指教。

希望过路神仙把竖曲线加进去带算隧道超欠挖，谢谢曲線任意里程中邊樁坐標正反算5800p本程序由内蒙@文哥和师父阿玖哥共同改编完成主程序 ZBJSLbl 6："1.SZ→ XY"："2.XY→ SZ"：？S：？N：Prog“SUB0”1÷P→C：(P-R)÷(2HPR) →D：180÷π→E:N=1 =>Goto 1: Goto 2Lbl 1:?Z: Abs(S-O)→W: Prog"SUB1":"XS=":X◢"YS=":Y◢"FS=":F-90 ►DMS◢“XC” ?C：“YC” ? E：0→I：0→J：Pol（X-C，Y-E）:“I=”：I◢If J∠0: Then “J=”: J+360►DMS◢IfEndIf J >0：Then “J=” ：J►DMS◢IfEnd：Goto 6Lbl 2?X:?Y:X→I:Y→J:Prog"SUB2":O+W→S:Cls:"S=":S◢"Z=":Z◢Goto 6正算子程序SUB10.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F：1-K→M：U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bc os(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD)))→X：V+W(Asin(G+QEKW(C+KWD))+Bsi n(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QE MW(C+MWD))) →Y：G+QEW(C+WD)+90→F：X+Zcos（F）→X：Y+Zsin（F）→Y反算子程序 SUB2G-90→T(Y-V)cosT-(X-U)sin(T) →WAbs(W)→W:0→ZLbl6：Prog "SUB1"T+QEW(C+WD) →L：(J-Y)cos（L）-(I-X)sin（L）→ZIF Abs（Z）<1E-6：Then0→Z:Prog "SUB1":(J-Y)÷sin（F）→Z:Else W+Z→W：G oto6:IfEnd数据库子程序SUB0Goto 1(线元可输入多条，分离式可在前多加一位，匝道一样。

线元法曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序（附带高程）一、功能及原理说明1. 功能说明：本程序由一个主程序(1-MAIN)和七个子程序——正算子程序(1-SUB-ZS)、反算子程序(1-SUB-FS)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。

本程序在CASIO FX-5800P计算器运行。

2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<=0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出反算结果。

3. 程序输入计算器后，请根据统计串列List X，Y，freq[N]在程序中实际可能被使用的维数，将其统计矩阵串列改为包含相应维数的任意数据的矩阵，即在list列表中随机输入相应维数的数据，保证矩阵大小符合维数要求。

4. 本程序正算速度在1-2秒左右，反算比正算慢点，可根据需要调整精度加快速度。

本程序可建立曲线要素数据库及高程变坡点数据库，一次输入整条线路数据，计算时自动调用数据库，实现全线贯通，也可临时手动输入线元计算要素进行计算。

5. 本程序由小骆在前人的基础上改进而成，经个人测试可用。

学识浅薄，不足之处，在所难免，欢迎提出改进意见。

二、源程序1.主程序(1-MAIN)Deg:fix 320→DimZLbl 0:cls：“INPUT(0),ZX(1),YX(2)”?I (选择手动输入计算参数0，还是调用线路数据1,2）If I=0:Then Prog “1-DAT1”:IfEnd“SZ=>XY(1),XY=>SZ(2)”?N（正算，反算）If N=1 :Then Goto 1Else If N=2 :Then Goto 2Else Goto 3IfEnd:IfEndLbl 1:“KP=”?S（输入待求桩号）If S<0:Then Goto 0:IfEnd“JL(m)=”?Z （输入偏距）If Z≠0:Then “ANGLE→R(Deg)=”?M:IfEnd （输入斜交右角）If I≠0:Then Prog “1-DAT2”:IfEndS-O→W:If W<0:Then Goto 0:Else If W>H:Then cls:locate 6,2,"KP OUT"◢Goto 0:IfEnd:IfEnd(前半条针对“DAT”情况，后半条针对“INPUT”情况。

C A S I O f x-5800P计算与道路坐标放样计算
一、前言
本程序是《CASIO fx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。

原道路坐标放样计算程序只基于道路的单个基本型曲线，有效计算范围仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。

改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是：
1．使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦；
2．坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）；
3．坐标反算方面，输入平面坐标，即可计算对应的桩号和距中距离（含左右信息）；
4．对于存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现
程序的特点：
1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写；
2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐；
3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。

程序代码
计算器操作正算说明记住要输入正算1
C匝道计算结果。

目录一、坐标正算基本公式………………………………………………………０２二、坐标反算原理 (04)三、高程数据库录入变换 (05)四、计算器程序………………………………………………………０7 ０1、ZBZS(坐标正算)………………………………………………………０702、ZBFS(坐标反算)………………………………………………………０803、ＧＣJF(高程积分） (09)04、PＪFＹ(坡脚放样)………………………………………………………1０05、ＪFCＸ（积分程序) (11)0６、ZＢFY(坐标放样)………………………………………………………１１０７、ＤT(递推） (12)08、H P(横坡） (13)０9、LK(路宽)………………………………………………………１410、SJK1(平面数据库) (14)11、SＪK２(纵面数据库)………………………………………………1４12、SJK3(左路宽度数据库） (15)１3、SＪＫ4(右路宽度数据库) (15)1４、SJＫ5(横坡数据库)......................................................1６１5、SＪK６(下边坡数据库) (16)16、SJK7(左上边坡数据库)………………………………………………1７17、ＳＪK８（右上边坡数据库）………………………………………………1８五、后记 (19)CASIO 58０0计算器公路工程测量程序一、正算所涉及得计算公式 X R i d XαβBd Y d l d βI图表 错误!不能识别的开关参数。

在图1中,A 点为回旋曲线起点,B 点为回旋曲线止点,I 点为所求坐标点。

设: Ａ点得X 坐标为X A ,Y 坐标为Y Ａ ,Ａ点得切线方位角为α,A 点得曲率为ρＡ,A 点得里程为ＬＡ，B 点得曲率为ρＢ,Ｂ点得里程为L Ｂ,I 点得曲率为ρI ,Ｉ点得里程为ＬI 。

fx-5800P道路中边桩坐标放样正反算程序ROAD-2（修订版）道路中边桩坐标放样正反算程序主程序 ROAD-2主程序：ROAD-2Deg:Fix 3:20→DimZ↓"ZS[1],FS[2]"?J↓If J=1:Then Goto 0:Else Goto 4:IfEnd↓Lb10↓"Z -1,Y +1"? I↓"LC"?P:"PJ"?D↓If I<0:Then Prog“ROAD-DATAZ” :Else Prog“ROAD-DATAY” :IfEnd↓M+O→N↓Prog"ROAD-SUB1"↓Prog"ROAD-SUB2"↓"X=":F+Dcos(Z+90)→Z[11]◢"Y=":G+Dsin(Z+90)→Z[12]◢"α=":Z DMS◢Goto 0↓Lbl 4↓"Z -1,Y +1"? I↓"X"?U:"Y"?V:"GSLC"?P↓Lbl 5↓If I<0:Then Prop“ROAD-DATAZ” :Else Prop“ROAD-DATAY” :IfEnd↓M+O→N↓Prog"ROAD-SUB1"↓Prog"ROAD-SUB2"↓Z-90→A↓(V-G)cos(A)-(U-F)sin(A)→H↓If Abs(H)>0.0001:Then P+H→P:Goto 5:IfEnd↓"LC=":P◢"PJ=":(G-V)÷sin(A)→D◢Goto 4↓平曲线要素及主点桩号计算子程序： ROAD-SUB1If O<0:Then -1→W:Else 1→W:IfEnd:WO→A↓B^(2)÷24÷R-B^(4)÷2688÷R^(3)→Z[6]↓C^(2)÷24÷R-C^(4)÷2688÷R^(3)→Z[7]↓B÷2-B^(3)÷240÷R^(2)+B^(5)÷34560÷R^(4)→Z[8]↓C÷2-C^(3)÷240÷R^(2)+C^(5)÷34560÷R^(4)→Z[9]↓Z[8]+(Z[7]-Z[6])÷(2×tg(A÷2))+(Z[7]+Z[6]+2×R)×tg(A÷2)÷2→S↓Z[9]+(Z[6]-Z[7])÷(2×tg(A÷2))+(Z[7]+Z[6]+2×R)×tg(A÷2)÷2→T↓RAπ÷180+( B+C)÷2→L↓L-B-C→Q↓(R+(Z[6]+Z[7])÷2)÷cos(A÷2)-R→E↓K-S→Z[1]↓Z[1]+B→Z[2]↓Z[2]+Q÷2→Z[3]↓Z[1]+L-C→Z[4]↓Z[4]+C→Z[5]↓中线坐标计算子程序： ROAD-SUB2X-Scos(M)→F:Y-Ssin(M)→G↓If P>Z[1]:Then Goto 1:IfEnd↓X→F:Y→G:P-K→Z[16]:0→Z[17]:M→A:M→Z:Goto 5↓Lbl 1↓If p>Z[2]:Then Goto 2:IfEnd↓P-Z[1]→L:L→Z[14]:B→Z[15]:Prog"ROAD-SUB3"↓WZ[17]→Z[17]:M→A:M+90WL^(2)÷(BRπ)→Z:Goto 5↓Lbl 2↓If P>Z[4]:Then Goto 3:IfEnd↓P-Z[1]→L:90(2L-B)÷R÷π→A↓Rsin(A)+Z[8]→Z[16]:W(R(1-cos(A))+Z[6]→Z[17]↓M+WA→Z:M→A:Goto 5↓Lbl 3↓If P>Z[5]:Then Goto 4:IfEnd↓Z[5]-P→L:L→Z[14]:C→Z[15]:Prog"ROAD-SUB3"↓X+Tcos(N)→F:Y+Tsin(N)→G↓—Z[16]→Z[16]:WZ[17]→Z[17]:N→A:N-90WL^(2)÷(CRπ→Z:Goto 5↓Lbl4↓X→F:Y→G:P-Z[5]+T→Z[16]:0→Z[17]:N→A:N→Z:Goto 5↓Lbl 5↓F+Z[16]cos(A)-Z[17]sin(A)→F↓G+Z[16]sin(A)+Z[17]cos(A)→G↓If Z<0:Then Z+360→Z:IfEnd↓缓和曲线局部坐标计算子程序 ROAD-SUB3If Z[14]=0:Then 0→Z[16]:0→Z[17]:Else↓Z[14]-Z[14]^(5)÷40÷(RZ[15])^(2)+Z[14]^(9)÷3456÷(RZ[15])^(4)→Z[1 6]↓Z[14]^(3)÷6÷(RZ[15])-Z[14]^(7)÷336÷(RZ[15])^(3)+Z[14]^(11)÷42240÷(RZ[15])^(5)→Z[17]↓IfEnd↓路线数据库子程序 ROAD-DATAZ12200.836→Z[18]:13481.073→Z[19]:If PZ[19]→Z[18]:14682.846→Z[19]:If PZ[19]→Z[18]:16353.635→Z[19]:If PZ[19]→Z[18]:17328.058→Z[19]:If PZ[19]→Z[18]:18555.915→Z[19]:If PZ[19]→Z[18]:19690.214→Z[19]:If PLbl 1:IfP<12200.836 Or P>19690.214:Then Cls:Locate 5,2,"LC ERROR"◢Stop:Ifend:Return路线数据库子程序 ROAD-DATAY17157.391→Z[18]:18529.883→Z[19]:If PZ[19]→Z[18]:19680.759→Z[19]:If PLbl 1:IfP<17157.391 Or P>19680.759:Then Cls:Locate 5,2,"LC ERROR"◢Stop:Ifend:Return。

FX-5800p计算器公路测量常用程序一、程序功能本程序由6个主程序、5个次子程序及5个参数子程序组成。

主要用于公路测量中坐标正反算，设计任意点高程及横坡计算,桥涵放样，路基开挖口及填方坡脚线放样。

程序坐标计算适应于任何线型.二、源程序1.主程序1：一般放样反算程序(①正算坐标、放样点至置仪点方位角及距离；②反算桩号及距中距离)程序名:1ZD-XYLb1 0:Norm 2F=1：(正反算判别，F=1正算，F=2反算,也可以改F前加？，改F为变量)Z[1]=90（与路线右边夹角）Prog＂THB＂：F=1=>Goto 1:F=2=>Goto 2Lb1 1:Ｆix 3:＂Ｘ=＂：Locate 6,4,Ｘ◢＂Ｙ＝＂：Locate 6,4,Ｙ◢Ｐrog＂3JS”:Goto 0:Lb1 2:Fix 3:＂ＫＭ＝＂：Locate 6,4,Ｚ◢Ｇoto 0２．主程序２：高程序横坡程序(设计任意点高程及横坡)程序名：2GCLbI 0:Norm 2“KM”?Z:?D:Prog”H”:Fix3:”H=”:Locate 6,4,H◢“I=”: Locate 6,4,I◢Goto 03.主程序3：极坐放样计算程序(计算放样点至置仪点方位角及距离)程序名：3JSX：Y：1268．123→K(置仪点X坐标)2243．545→L（置仪点Y坐标，都是手工输入,也可以建导线点数据库子程序,个人认为太麻烦）Y-L→E：X-K→F：Pol(F,E):IF J<0:ThenJ+360→J:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J:(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示) Fix 4:” FWJ=”: Locate 6,4,J◢(不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢来直接显示)Fix 3:”S=”:Locate 6,4,I◢4．主程序4：涵洞放样程序（由涵中心桩号计算出各涵角坐标、在主程序3中输入置仪点坐标后计算放样点至置仪点方位角及距离)程序名：4JH-XYLbI 0:Norm 290→Z[1](涵洞中心桩与右边夹角，手工输入，也可以修改成前面加？后变为变量)1→F:Prog”THB”:?L:Z[2]-Z[1] →E:X+Lcos(E) →X:Y+Lsin(E) →Y:Fix 3: ＂Ｘ=＂：Locate 6,4,Ｘ◢＂Ｙ＝＂：Locate 6,4,Ｙ◢Ｐrog＂3JS”:Goto 0:5．主程序5:路基开挖边线及填方坡脚线放样程序（输入大概桩号及测量坐标、地面标高计算出偏移距离、桩号、距中距离、填挖高度）程序名：5FBXLbI 0:Norm 2: 18→DimZ:2→F:90→Z[1]:Prog“THB”:Z:D:”M0”?M:M→Z[4]:D→Z[3]:Prog”6GD”:L→Z[6]:If D<0:Then 0.75-L→D:GotoH:Else L-0.75→D:Goto H:IfEndLbI H:Prog”H”:H-0.03-Z[4] →Z[5]:Z[6] →L:If Z[5]<0:Then –z[5] →G:Goto W:Else Z[5] →G:Goto T:LbI W:Prog “W0”:Z[10]+Z[11] →A: If G>A:Then Goto 1:Else If G>Z[10]:Then Goto 2:Else Goto 3:IfEnd:LbI 1:L+Z[12]+Z[13]+Z[14]+(G-A)×Z[9]+Z[11]×Z[8]+Z[10]×Z[7]:Goto Z:LbI 2:L+Z[12]+Z[13]+(G-Z[10])×Z[8]+Z[10]×Z[7]:Goto Z:LbI 3:L+Z[12]+G×Z[7]:Goto z:LbI T:L+0.5→N:If G>Z[17]:Then (N+Z[18]+(G-Z[17])×Z[16]+Z[17]×Z[15])→S:Goto Z:Else (N+G×Z[15])→S:Goto z:LbI Z:Z[3]→D:Fix 2:Abs(D)-S→T:”L0=”:L Locate 6,4,T◢＂Ｄ＝＂：Locate 6,4,Ｄ◢“TW=”: Locate 6,4,Z[5]◢Goto 06．主程序6：路基标准半幅宽度计算程序(对于设计有加宽渐变的有用，如路基宽度无变化，则把此程序直接输入半幅宽度值至L)程序名：6GDProg “G0”Z-C→E:(B-A)×E/S+A→L:L:7．坐标计算次程序（ＴＨＢ）程序名：THB18→DimZ：＂ＫＭ＂?Z：Prog ＂Ｘ０＂1÷P→Ｃ： (P-R)÷(2HPR) →Ｓ：180÷π→Ｅ：Ｆ=1=>Goto 1：Ｆ=2=>Goto 2←┘Lbl 1：?Ｄ： Abs(Ｚ－O) →W：Prog "Ａ"：X：Ｙ：Ｇoto 3LbI 2:Ｘ：Ｙ：Ｘ→I：Ｙ→J：Prog ＂Ｂ＂：Ｏ＋Ｗ→Z：Ｄ→Ｄ：Ｇoto 3LbI 3:IF Ｆ＝１Then Ｘ：Ｙ：Else Ｚ：Ｄ8. 正算子程序(Ａ)程序名：A0.1184634425→A： 0.2393143352→B：0.2844444444→N 0.046910077→K：0.2307653449→L：0.5→M： U+W(Acos(G+QEKW(C+KWＳ))+Bcos(G+QELW(C+LWＳ))+Ncos(G+QEMW(C+MWＳ))+Bcos(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))) →X：V+W(Asin(G+QEKW(C+KWＳ))+Bsin(G+QELW(C+LWＳ))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)WＳ))) →Y：G+QEW(C+WＳ)+Ｚ［１］→Ｚ［２］：X+Ｄcos（Ｚ［２］）→Ｘ： Y+Ｄsin（Ｚ［２］）→Ｙ9. 反算子程序(Ｂ)程序名：BG-90→Ｔ： Abs((Y-V)cos(T)-(X-U)sin(T)) →Ｗ：0→Ｄ：Lbl 0：Prog "Ａ"： T+QEW(C+WＳ) →Ｌ： (J-Y)cos(L)-(I-X)sin(L)→Ｄ：IF Abs(Ｄ)<0.01:Then Goto1：Else W+Ｄ→W:Goto 0←┘Lbl 1:0→D：Prog "Ａ"：(J-Y)÷sin(Ｚ［２］) →D:10．高程计算子程序（H）程序名：HProg “S0”:R:T:C:G:I:C-T→F:Z-F→L:C+T→E:G-TI→Q:If T=O:Then Q+LI→H:Goto 0:Else IfZ<F:Then Q+LI→H:Goto 0:Else If Z≤E:Then Q+LI+L2÷2÷R→H:Goto 0:LbI 0:H:If D=0:Then Goto I:Else Prog “I”:H+V→H:Goto I:LbI I:H:I:11.高程超高计算程序（I）程序名：IProg”I0”:W=1=> Goto 0:W=2=>Goto 1:LbI 0:If L=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×((N-M)×(Z-C)÷L+M)→V:Goto 2:IfEnd: LbI 1:If L=0:Then Abs(D)×M→V:Goto 2:ElseAbs(D)×(((3((Z-C)÷L)2-2((Z-C)÷L)∧(3))×(N-M))+M)→V:Goto 2:IfEnd:LbI 2:Abs(D)→E:V÷E→I:I(E-K)→V:12．数据子程序(附后示例)①程序名：X0（坐标计算要素程序）If Z≥25900 And Z≤26615.555:Then 25900→O:11587.421→U:1847.983→V:101。

道路直线段中边桩坐标放样计算CASIO fx-5800P
程序
发表于2009年04月05日23:21阅读（200）评论（14）分类:
举报
王中伟 技术交流
三、数学模型
如图所示，已和輕甲元的起点坐标0（%必）、起点切絃方位角旳、直统甲元伙產z・
.V
亘妹孕元的终点（或直裟上呆点〉坐标ZU..X；）计算如下］
£=兀■厶COS偽I y.m + Zfink f
•四、程序代码
五、程序变量清单
道苗
五、计算流程示例需要对设计文件中的逐桩坐标进行验算，将断链前后两段路线分别进行计算。

1.计算断链前路段的中桩坐标使用道路直线段中边桩坐标放样计算程序LINE 的操作流程见下表。

2-计算断链后路段的中桩坐标使用道路直线段中边桩坐标放样计算程序LINE的操作流程见下表。

日志转到我空间分享给好友复制网址隐藏签名档小字体选用该信纸上一篇下一篇返回日志列表转载自钱¥霸转载丁2010年02月01日10:09阅读（2）评论（0）分类：工程技术丛书举报道路中边桩坐标放样正反算CASIO fx-5800P程序（全线贯通）一、前言本程序是《CASIOx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。

原道路坐标放样计算程—仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。

改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是：1. 使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦；2. 坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）;3. 坐标反算方面，输入平■面坐标，即可计算对应的桩号和距中距离（含左右信息）；4. 对丁存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现（可参照立交匝道程序中匝道的选择）。

程序的特点：1. 可进行中桩坐标的正、反算，程序代码简洁，便丁阅读和改写；2. 主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐；3. 使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。

二、道路示例项目基本资料基本资料同《CASIO fx-5800P计算与道路坐标放样计算〉〉第6章HY高速公路第2合同段（合同段起止桩号：K4+800~K9+600。

这里摘取直线、曲线及转角表资料如下（若图片不活晰，请参见参见教材P161附录1）:三、程序代码注：路线数据库子程序ROAD DATA促根据计算示例项目的直曲表编写，大家使用时应按各自项目的直曲表改写或新建。

卡西欧5800公路坐标正反算程序卡西欧5800公路坐标正反算程序是一种用于测量和计算公路坐标的工具，可以帮助测量员准确地确定公路上的位置和坐标。

本文将详细介绍卡西欧5800公路坐标正反算程序的功能、使用方法以及相关注意事项。

一、功能介绍卡西欧5800公路坐标正反算程序具有以下主要功能：1. 坐标正算：根据已知的起点坐标、方位角和距离，计算出目标点的坐标。

2. 坐标反算：根据已知的起点坐标和目标点坐标，计算出方位角和距离。

3. 坐标转换：支持不同坐标系之间的转换，如高斯投影坐标系、经纬度坐标系等。

4. 坐标纠正：根据已知的控制点坐标和实测点坐标，进行坐标纠正和校正。

二、使用方法1. 坐标正算：步骤1：进入卡西欧5800公路坐标正反算程序。

步骤2：输入起点坐标（经度、纬度）、方位角和距离。

步骤3：点击计算按钮，程序将自动计算并显示目标点的坐标。

2. 坐标反算：步骤1：进入卡西欧5800公路坐标正反算程序。

步骤2：输入起点坐标（经度、纬度）和目标点坐标（经度、纬度）。

步骤3：点击计算按钮，程序将自动计算并显示方位角和距离。

3. 坐标转换：步骤1：进入卡西欧5800公路坐标正反算程序。

步骤2：选择需要转换的坐标系类型。

步骤3：输入待转换的坐标信息。

步骤4：点击转换按钮，程序将自动进行坐标转换并显示结果。

4. 坐标纠正：步骤1：进入卡西欧5800公路坐标正反算程序。

步骤2：输入已知的控制点坐标和实测点坐标。

步骤3：点击纠正按钮，程序将自动进行坐标纠正和校正，并显示结果。

三、注意事项1. 在使用卡西欧5800公路坐标正反算程序之前，确保已正确安装和配置软件，并且设备连接正常。

2. 在进行坐标计算和转换时，输入的坐标信息应准确无误，以确保计算结果的准确性。

3. 在进行坐标纠正时，应选择合适的控制点和实测点，并确保其坐标信息准确无误。

4. 在使用过程中，如遇到问题或需要进一步了解功能细节，可以参考使用手册或联系卡西欧客服进行咨询。

卡西欧fx-5800P坐标正反算、高程及边坡放样程序《ZS》正算主程序Lb1 0：“K=”?S(所求里程)：“PJ=”?Z(偏距，左偏负右偏正)：Prog “PM-SJ”： Abs(S-O) → W：Prog "SUB1"↙Fix 3： "XS="：X◢（计算得X坐标）"YS="：Y◢（计算得Y坐标）Norm 2："FW=":F-90→F：F►DMS◢（所求里程设计线处切线方位角）Prog“SQX-SJ”： Prog“SQX”： Fix 3: “H=”：H◢（设计线处该里程的设计高程）Norm 2：Goto 0《FS》反算主程序Lb1 0：“GK=”?S（所求点的估计里程，目的调取数据库）：?X：?Y：Prog“PM-SJ”： X →I： Y→J：Prog "SUB2"：Fix 3:"K="： O+W→S◢（所求点的里程）"PJ="： Z◢(所求点的偏距，左偏负右偏正)Prog“SQX-SJ”：Prog“SQX”：“H=”：H◢（反算求的里程设计线处设计高程）Norm 2↙Goto 0（此处的Goto 0做《LJWF-2JBP》路基挖方2级边坡放样的子程序时要去掉）《SUB1》正算子程序1÷P→C：(P-R)÷(2HPR)→D：180÷π→E：0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F：1-K→M：U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD))) →X：V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) → Y：G+QEW(C+WD)+90→F：X+Zcos(F)→ X：Y+Zsin(F)→Y《SUB2》反算子程序G-90→T： Abs((Y-V)cos(T)-(X-U)sin(T)) → W： 0→ Z：Lbl 0：Prog "SUB1"：T+QEW(C+WD) →L：(J-Y)cos(L)-(I-X)sin(L) → Z：If AbsZ<0.000001：then Goto1：Else W+Z→W：Goto 0：If End←┘Lbl 1：0→Z：Prog "SUB1"：(J-Y)÷sin(F)→ Z《PM-SJ》子程序（平面线形数据库）If S≥34600（线元起点里程）And S＜35613.616（线元止点里程）：Then 2467949.131 →U（线元起点X坐标）:536385.574 →V（线元起点Y坐标）：34600 →O（线元起点里程）：202°10°11.1°→G（线元起点切线方位角）：1013.616 →H（线元长度）：10^（45）→P （线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正，直线为0）：If End ↙If S ≥35613.616And S＜35723.616：Then 2467010.452→U：536003.084→V：35613.616→O：202°10°11.1°→G：110→H：10^（45）→P：2800→R：-1→Q：If End←┘If S≥35723.616And S＜36520.648：Then 2466908.316→U：535962.244→V：35723.616→O：201°02°39.47°→G：797.032→H：2800→P：2800→R：-1→Q：If End←┘If S≥36520.648 And S＜36630.648:Then 2466133.988 →U:535785.048→V：36520.648 →O:184°44°5.31°→G：110 →H：2800→P：10^（45）→R：-1 →Q：If End ↙If S ≥36630.648And S＜36687.250：Then2466024.256→U：535777.405→V：36630.648→O：183°36°33.68°→G：56.602→H：10^（45）→P：10^（45）→R：0→Q：If End ←┘（可以一直往下加数据库）《SQX-SJ》子程序（竖曲线数据库）If S≤36687.250（竖曲线终点里程）：Then +0.0053（前坡度：上坡正下坡负）→E：-0.004（后坡度：上坡正下坡负）→F：15000→R（竖曲线半径）：35500→G（变坡点里程）：98.920→C（变坡点高程）：If End←┘（可以一直往下加数据库）《SQX》子程序（竖曲线计算公式）E-F→J：If J＞0：Then-R→R：Else R→R：Abs（R×J÷2）→T：If End←┘If S≤G-T：Then G-S→L：C-LE→H：If End←┘If S＞G-T And S≤G： Then G-S→L：（S+T-G）^2÷2÷R→N：C-LE+N → H：If End←┘If S＞G And S≤G+T： Then S-G→L：（G +T- S）^2÷2÷R→N：C+LF+N → H：If End←┘If S＞G+T ：Then S-G→L： C+LF → H：If End←┘三、使用说明(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

5800计算器公路全线坐标正、反算计算程序FX5800全线贯通万能正、反算程序（一体化、超好用、短小、易懂）FX5800计算器的积分程序（正反算、全线贯通、新线路）终极版ZHUCHENGXU主程序"1.ZS,2.FS" ?→Q←┘输入1正算，输入2反算“NEW=0,OLD≠0”?Z←┘IfZ=0:Then“X0=”?A:“Y0=”?B:“C0=”?C:“1/R0=”?D:“1/RI=”?E:“SP=”?F:“EP=”?G:Ifend:Q=2=>Goto 2←┘Lbl1 :“KM=,<0 Stop”?H:H<0=>Stop:“PJ=”?O:“PY=”?L←┘LblZ:Z=1=> Prog“01”：Z=2=> Prog“02”←┘选择数据库文件，可增加H- F→X:0.5(E-D)÷(G-F)→N←┘C+(XD+NX2)*180÷π→P：P<0=>P+360→P：P>360=>P-360→P←┘-A+∫(cos(C+(XD+NX2)*180÷π),0,X)+Lcos(P+O)→U←┘B+∫(sin(C+(XD+NX2)*180÷π),0,X)+Lsin(P+O)→V←┘Q=2=>Goto4：Cls：Fix 3←┘"Xn=":Locate4，1，U："Yn=": Locate5，2，V：“FWJ=”：PDMS◢Norm 2：Cls：Goto1←┘Lbl2:“XD=,<0,STOP”？R：R<0=>Stop:“YD=”？S←┘“KMDG=”？H ：90→O:0→L：GotoZ←┘（H线路范围内的任意桩号）Lbl4:Pol（R-U，S-V）：J<0 => J+360→J←┘Whileabs（Icos（J-P））≤0.001:P-J>180=> J+360→J： P-J<-180=> P+360→P:If P-J>0:then -I→L:else I→L :Ifend:Goto3: Whileend：H+Icos（J-P）→H:GotoZ←┘Lbl3:Cls：Fix 3←┘“KM=”: Locate4，1，H：“PY=”: Locate4，2，L◢Norm 2：Cls：Goto2←┘01（数据库子程序）If H<=第一曲线终点桩号：then第一曲线起点X→A：第一曲线起点Y→B：第一曲线起点方位角→C：起点曲率→D：终点曲率→E：起点桩号→F：终点桩号→G：return：ifend ……………程序说明：1、该程序可以计算任意线形（直线、圆曲线、缓和曲线、不完整曲线）任意桩号的坐标（正算，输入1），也可根据坐标计算该点到线路的距离及垂足桩号（反算，输入2）；2、（NEW=0，OLD≠0）？如果要计算的点为数据库线路中的点，则输入数据库编号（以整数1、2、3…代替输入）；如果在数据库中没有要计算线路的数据，则输入曲线要素X0：曲线起点X坐标；Y0：曲线起点Y坐标；C0:曲线起点方位角；R0-1、 RI-1:曲线起点、终点曲率，直线为0，曲线左偏输入负值，右偏输入正值；SP、 EP：曲线起点桩号，终点桩号；KM：待求点桩号；PJ：正斜交的设定；PY：偏中距离，线路上的点输入0，右偏输入+值，左偏输入-值；3、正算显示坐标及切线方位角；反算输入线路的任意桩号（此桩号越接近真实值计算速度越快）、待求点坐标，显示待求点桩号及偏中距离；4、正算子程序为积分公式编写而成；反算子程序为角度趋近的方法编写，计算速度有点慢。

CASIO fx-5800P线元法坐标正反算程序说明：本程序适用于卡西欧计算器 CASIO fx-5800P,可对全线贯通坐标正反算、竖曲线高程计算。

该程序可计算任意线型，包含（直线、圆曲线、缓和曲线、卵形曲线）等,还可以能通过坐标反推该点里程和距中线距离,适用测量员专用。

主程序名：ABCYT第1步Deg：Fix 3：10→DimZ第2步Lbl 3：＂1.DK=>XY＂：＂2.XY=>DK＂：＂Q＂?W：＂DK＂?S：Prog＂ABCYTSJ＂：If P=0：Then 10^(45)→P：IfEnd：If R=0：Then 10^(45)→R：IfEnd第3步1÷P→C：(P-R)÷(2HPR)→D：180÷π→E：If W=1：Then Goto 1：Else Goto2：IfEnd第4步Lbl 1：＂W＂?Z：＂α＂?N：Abs(S-O)→W：Prog ＂ABCYTZ＂第5步Cls：＂F=＂：Locate 3,1,F°：＂X=＂：Locate 3,2,X：＂Y=＂：Locate 3,3,Y◢第6步Prog＂ABCYTSQX＂：Cls：＂H=＂：Locate 3,1,H◢第7步1→W：90→N：Goto 3第8步Lbl 2：?X：?Y：X→I：Y→J：Prog＂ABCYTF＂：O+W→S第9步Cls：＂K=＂：Locate 3,1,S：＂S=＂：Locate 3,2,Z◢第10步2→W：Goto 3正算子程序名：ABCYTZ第1步0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L第2步1-L→F：1-K→M第3步U+W×(A×cos(G+Q×E×K×W×(C+K×W×D))+B×cos(G+Q×E×L×W×(C+L×W×D))+B×cos(G+Q×E×F×W×(C+F×W×D))+A×cos(G+Q×E×M×W×(C+M×W×D)))→X第4步V+W×(A×sin(G+Q×E×K×W×(C+K×W×D))+B×sin(G+Q×E×L×W×(C+L×W×D))+B×sin(G+Q×E×F×W×(C+F×W×D))+A×sin(G+Q×E×M×W×(C+M×W×D)))→Y第5步G+Q×E×W×(C+W×D)→F：F+ N→Z[1]第6步X+Z×cos（Z[1]）→X：Y+Z×sin（Z[1]）→Y反算子程序名：ABCYTF第1步Lbl 2：（S-O）→W：0→Z：Prog ＂ABCYTZ＂：F-90→Z[9]：(J-Y)×cos(Z[9])-(I-X)×sin(Z[9])→Z[10]第2步If Abs(Z[10])>0.001：Then S+Z[10]→S：Goto 2：Else Goto 1：IfEnd第3步Lbl 1：(Y-J)÷sin(Z[9])→Z第4步Pol(X-I，Y-J)：If Z<0：Then -1×I→Z：Else 1×I→Z：IfEnd数据库名：ABCYTSJ第1步Goto1第2步Lbl 1：If S<7586.707 Or S>13346.96：Then Cls：Locate 2,2,＂PQX＂：Locate 4,3,＂CHAOXIAN＂：Locate 10,4,＂→Stop＂◢第3步Stop：IfEnd第4步Lbl 1：IF S<7946.707：Then 98°56′56″→G：7586.707→o：3378605.445→U：453648.704→V：0→P：4500→R：360→H：1→Q：Return：IfEnd第5步Lbl 1：IF S<11766.03：Then 101°14′26″→G：7946.707→o：3378544.714→U：454003.518→V：4500→P：4500→R：3819.323→H：1→Q：Return：IfEnd第6步Lbl 1：IF S<12126.03：Then 149°52′11″→G：11766.03→o：3376389.890→U：457018.324→V：4500→P：0→R：360→H：1→Q：Return：IfEnd第7步Lbl 1：IF S<13346.96：Then 152°09′41.68″→G：12126.03→o：3376073.846→U：457190.654→V：0→P：0→R：1220.93→H：0→Q：Return：IfEnd第n步……………………………………………………数据输入说明：第1步Goto1第2步Lbl 1：If S<本条线路起点里程S>本条线路止点里程Cls：Locate 2,2,＂PQX＂：Locate 4,3,＂CHAOXIAN＂：Locate 10,4,＂→Stop＂◢第3步Stop：IfEnd第4步Lbl 1：If S<本线元止点里程：Then线元起点切线方位角→G：线元起点桩号→O：线元起点坐标X→U：线元起点坐标Y→V：线元起点半径（直线为0、曲线为半径）→P：线元止点半径（直线为0、曲线为半径）→R：线元长度→H：线元转向（左转为-1、右转为1、直线为0）→Q：Return：IfEnd第n步Lbl 1：每增加一行则为增加一个线元要素。

CASIO fx—5800P计算器路线坐标计算程序（单个交点）主程序名CALXY （计算中桩、边桩坐标）Lbl 1：“ZZ”? →Z[27] ：“D(-+m)”? →D：If D≠0 ：Then “∠(DMS)”? →V：Ifend』If Z[27]≤Z[17]：Then Z[18]+ (Z[27] —Z[17]) Cos( F)+D Cos( F+V)：“X=”：Ans▲Z[19]+ (Z[27] —Z[17]) Sin( F )+D Sin( F+V)：“Y=”：Ans▲：Ifend』If Z[27]＞Z[17] And Z[27]≤Z[20]：Then Z[27] —Z[17] →L：Prog “HHXY”：Z[18] +I Cos( F+ JH）+D Cos( F+WH+V)：“X=”：Ans▲：Z[19] +I Sin(F+ JH）+D Sin( F+WH+V)：“Y=”：Ans▲：Ifend』If Z[27]＞Z[20] And Z[27]≤Z[24] —Z[6]：Then Z[27] —Z[20] →L：Prog “YUXY”：Z[21] +I Cos( Z[23]+ JH）+D Cos( Z[23]+WH+V)：“X=”：Ans▲：Z[22] +I Sin(Z[23]+ JH）+D Sin( Z[23]+WH+V)：“Y=”：Ans▲：Ifend』If Z[27]＞Z[24] —Z[6] And Z[27]≤Z[24]：Then Z[24] —Z[27] →L：Prog “HHXY”：Z[25] +I Cos( C—JH）—D Cos( C—WH+V)：“X=”：Ans▲：Z[26] +I Sin(C—JH）—D Sin( C—WH+V)：“Y=”：Ans▲：Ifend』If Z[27] ＞Z[24]：Then Z[25] —(Z[27] —Z[24]) Cos( C )—D Cos( C+V)：“X=”：Ans▲Z[26] —(Z[27] —Z[24]) Sin( C —D Sin( C+V)：“Y=”：Ans▲：Ifend』Goto 1子程序名JD （输入曲线参数，计算曲线要素）30→Dim Z：Fix 5：“JD”? →Z[1] ：“XJD”? →Z[2] ：“YJD”? →Z[3]：“F0(DMS)”? →F：“A（-+ DMS）”？→Z[4]：“R”? →R：“LH1”? →Z[5]：“LH2”? →Z[6]：√￣(Z[4]2)→A：Int(Z[4] ÷A)→H』Z[5] →S：Prog “PQ”：Q →Z[7]：P→Z[8]：B →Z[9]：R →S：Prog “PQ”：Q →Z[10]：P→Z[11] ：B →Z[12] 』（R+ Z[8]）Tan(A÷2)+ Z[7] —(Z[8] —Z[11]) ÷Sin(A)→Z[13] ：“T1=”：Ans ▲（R+ Z[11]）Tan(A÷2)+ Z[10] —(Z[8] —Z[11])÷Sin(A)→Z[14] ：“T2=”：Ans▲』√￣((Z[13] —Z[7])2+( R+ Z[8])2) —R →Z[15]：“E=”：Ans▲R( A —Z[9] —Z[12] )÷57.2958+ Z[5] + Z[6] →Z[16]：“L=”：Ans▲』Z[1] —Z[13]→Z[17]：Z[2] —Z[13]Cos( F) →Z[18]：Z[3] —Z[13] Sin( F) →Z[19]：If Z[5]＞0 ：then Z[17] +Z[5]→Z[20]：Z[5] →L：Prog “HHXY”：Z[18] +I Cos( （F+ JH）)→Z[21]：Z[19] +I Sin( （F+ JH）)→Z[22]：（F+ WH）→Z[23] ：Else Z[17]→Z[20]：Z[18]→Z[21]：Z[19]→Z[22]：F→Z[23]：Ifend』Z[17] + Z[16]→Z[24]：Z[2] +Z[14] Cos( （F+ Z[4]）)→Z[25]：Z[3] +Z[14] Sin( （F+ Z[4]）)→Z[26]：F+ Z[4] —180→C』Prog “CALXY”』子程序名PQS÷2—S3÷（240 R2）→Q：S2÷（24 R）—S4÷（2384 R3）→P：28.6479S ÷R→B』子程序名HHXYL—L5÷(40S2 R2)+L9÷(3456 R4S4) →X：L3÷（6S R）—L7÷（336S3 R3）+L11÷（42240 R5S5）→Y：28.6479L2÷(S R) →W：Pol（X,Y）』子程序名YUXY57.2958L÷R →W：RSin(W)→X：R（1—Cos(W）)→Y：Pol（X,Y）』CASIO fx—5800P计算器路线坐标计算程序（单个交点）使用说明首先运行子程序“JD”，输入曲线参数，计算曲线要素。

CASIO fx-5800P ZXZFS直线坐标正反算程序V1.5说明：该程序适用于计算器 CASIO fx-5800P直线坐标正算、反算，可连续计算中桩及边桩坐标其斜交等，还可以反算线路外一点距线路中心的里程和垂距。

增加高程计算时错误主程序名：ZXZFS第1步 Deg：Fix 3“1.JL”：“2.2HJ”?→VV=1=>PROG “JL”:V=2=>PROG “2HJ”。

第2步 Lbl C：＂1.DK=>XY＂：＂2.XY=>DK＂：？W：Lbl 1：＂K0+＂？E第3步 If W=1：Then Goto A：Else If W=2：Then Goto B：Else Goto C：IfEnd第4步 Lbl A：＂E＂？G90=H (90为与中线夹角)第5步 B+(E-A)×cos(F)+G×cos(F+H)→X第6步 C+(E-A)×sin(F)+G×sin(F+H)→Y第7步 Cls：＂X=＂：Locate 3,1,X：＂Y=＂：Locate 3,2,Y◢第8步＂ZJD(X)＂？S：＂ZJD(Y)＂？T第9步 Pol(X-S,Y-T)：I→M第10步 If J<0：Then J+360→N：Else J→N：IfEnd第11步 Cls：＂JL=＂：Locate 4,1,M：＂FY=＂：Locate 4,2,N°◢第12步 Goto C第13步 Lb1 B：＂X=＂？K：＂Y=＂？L第14步 Pol(K-B，L-C)：I→D：If J<0：Then J+360→P：Else J→P：IfEnd第15步 If F<P：Then P-F→Ｏ：Else F-P→Ｏ：IfEnd第16步 Sin(Ｏ)×D→Q第17步 A+Q÷tan(Ｏ)→R第18步 sin(Ｏ)×D→Q：If F<P：Then 1×Q→Q：Else -1×Q→Q：IfEnd第19步 Cls：＂DK=＂：Locate 4,1,R：＂PJ=＂：Locate 4,2,Q◢第20步 Goto C数据库名：ZXZFSSJ第1步 Goto 1第2步 Lb1 1:If E<259105:Then 256100→A:2575647.999→B:495400.3447→C:221°47ˊ20.24″→F:Return:IfEnd第n步……………………………………………………数据库说明：第1步 Goto 1第2步 Lb1 1:If E<7586.706:Then 7152.5561→A:3378672.9776→B:453219.8377→C:98°56ˊ55.65″→F:Return:IfEnd第2步 Lb1 1:If E<止点桩号:Then 起点桩号→A:起点坐标X→B:起点坐标Y→C:起点方位角F→F:Return:IfEnd第n步 Lbl 1：每增加一行则为增加一个竖曲线要素。

CASIO fx-5800P坐标正反算程序[已核对]一、程序内容1、主程序名:ABCYTDeg:Fix 3Lbl 4:＂1.DK=>XY＂:＂2.XY=>DK＂:＂J＂?N:＂DK＂?S:Prog＂CYT1＂↙1÷P→C:(P-R)÷(2HPR)→D:180÷π→E:N=1=>Goto1:Goto2:↙Lbl 1:＂W＂?Z:Abs(S-o)→W:Prog＂CYT2＂:＂X=＂:X◢＂Y=＂:Y◢F-90→F:＂FWJ=＂:F▲DMS◢Goto4↙Lbl 2:?X:?Y:X→I:Y→J:Prog＂CYT3＂:o+W→S:＂DK=＂:S◢＂W=＂:Z◢Goto4↙2、正算子程序名:CYT20.1739274226→A:0.3260725774→B:0.0694318442→K:0.3300094782→L:1-L→F:1-K→M:U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C +MWD)))→X:V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C +MWD)))→Y:G+QEW(C+WD)+90→F:X+Zcos（F）→X:Y+Zsin（F）→Y↙3、反算子程序名:CYT3G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T))→W:0→Z:Lbl 0:Prog ＂CYT2＂:T+QEW(C+WD)→L:(J-Y)cos（L）-(I-X)sin（L）→Z:IF Abs（Z）<10^(-6):Then 0→Z:Prog ＂CYT2＂:(J-Y)÷sin（F）→Z↙Else W+Z→W :Goto 0:IfEnd↙4、数据库名:CYT1Goto1↙Lbl 1:IF S<7946.707:Then 98°56′56″→G:7586.707→o:3378605.445→U:453648.704→V:10^(45)→P:4500→R: 360→H: 1→Q:Return:IfEnd↙Lbl 1:IF S<11766.03:Then 101°14′26″→G:7946.707→o:3378544.714→U:454003.518→V:4500→P:4500→R: 3819.323→H: 1→Q:Return:IfEnd↙Lbl 1:IF S<12126.03:Then 149°52′11″→G:11766.03→o:3376389.890→U:457018.324→V: 4500→P: 10^(45)→R: 360→H: 1→Q:Return:IfEnd↙Lbl 1:IF S<13346.96:Then 152°09′41.68″→G:12126.03→o:3376073.846→U:457190.654→V:10^(45)→P: 10^(45)→R:1220.93→H:0→Q:Return:IfEnd↙CYT1数据库说明:Goto1↙表示同时保存多个曲线要素Lbl 1:IF S<***线元终点里程:Then***线元起点方位角→G:***线元起点里程→O:***线元起点X→U:***线元起点Y→V:***线元起点曲率半径→P:***线元终点曲率半径→R: ***线元起点至终点长度→H:0→Q:Return:IfEnd↙（H=0或1、-1，详情附后）。