卡西欧5800公路坐标正反算程序

CASIO5800坐标正反算计算程序CASIO 9860/5800坐标正反算计算程序1. 坐标正算主程序(命名为ZBZS)第1行：Lbl 0:”K=”?K:”BIAN=”? Z:”α=”?B第2行：Prog “A”第3行：”X=”:N+Zcos(F+B)◢第4行：”Y=”:E+Zsin(F+B)◢第5行：”F=”:F◢第6行：Goto 0K——计算点的里程BIAN——计算点到中桩的距离（左负右正）α——取前右夹角为正2. 坐标反算桩号和偏距主程序(命名为ZBFS)第1行：”X1=”? C:”Y1=”?D: ”K1=”?K第2行：Lbl 0:Prog “A”第3行：Pol(C-N,D-E)第4行：List Ans[1]→I第5行：List Ans[2]→J第6行：Icos(F-J)→S:K+S→K第7行：Abs(S)>0.0001=>Goto 0第8行：”K1=”:K◢第9行：”BIAN=”:Isin(J-F)→Z◢X1——取样点的X坐标Y1——取样点的Y坐标K1——输入时为计算起始点(在线路内即可)，输出时为反算点的桩号Z——偏距（左负右正）3. 计算坐标子程序(命名为XYF)为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第1行：K-A→S：(Q-P)÷L→I第2行：N+∫(cos(F+X(2P+XI)×90÷π),0,S)→N第3行：E+∫(sin(F+X(2P+XI)×90÷π),0,S)→E第4行：F+S(2P+SI)×90÷π→F4. 数据库（命名为A）第1行：K≤175.191=>Stop第2行：175.191→A:428513.730→N:557954.037→E:92°26′40″→F:0→P:1/ 240→Q:70.417→L:K≤A+L =>GoTo 1（第一缓和曲线，圆半径为240）第3行：245.607→A: 428507.298→N:558024.092→E: 100°50′59.4″→F: 1/240→P:1/240→Q:72.915→L: K≤A+L =>Goto 1（第圆曲线，半径为240）第4行：318.522→A: 428482.988→N:558092.538→E: 118°15′25.2″→F: 1/240→P: 0→Q: 55.104→L: K≤A+L =>Goto 1（第二缓和曲线，圆半径为240）第5行：373.627→A:428453.283→N:558138.912→E:124°50′4.5″→F:0→P:-1/180→Q:67.222→L:K≤A+L=>Goto 1:Stop（第一缓和曲线，圆半径为180）第6行：Lbl 1:Prog “XYF”A——曲线段起点的里程N——曲线段起点的x坐标E——曲线段起点的y坐标F——曲线段起点的坐标方位角P——曲线段起点的曲率(左负右正)Q——曲线段终点的曲率(左负右正)L——曲线段长度（尽量使用长度，为计算断链方便）说明：（1）在9860中，程序中所有公式和部分函数结果均存储在ListAns列表数组中，要想多次调用最好随公式取出结果，并赋给变量。

1.主程序(1-ZB ZFS)Deg:Fix 6Lbl 9:”I”:”2.YOU JIN”:”MI MA =”?C: Goto 1↙”Lbl 4："3.ZH =>BG"：”1.ZH =>XY”2.XY =>ZH "：”3.ZH =>BG 4.TF =>JS ”:"N="?N↙If N =3;Then Prog “BG JS 20M”: IfEnd↙If N =4: Then Prog “4-TFJS”; IfEnd↙If N ≥5: Then Cls:Stop; IfEnd↙Lbl 9:”ZH=”?S:Prog “1.1-ZFS”:If N=1:Then Goto3:IfEnd↙If N=2:Then Goto2:IfEnd↙Lbl3:”Z=”?Z:”J=”?J:Abs(S-O)→W:Prog”1.2-ZFS”:X+Zcos(F +J)→X:Y+Zsin(F+J) →Y↙If Z=0(带下划线的0为数字0):Then “X=”:X◢"Y="：Y◢F-90→F:“F WJ=”：F▼DMS◢(5800里面的度分秒的)Prog “123”：“BG =>Z”:D◢Goto4: IfEnd↙If Z=1: Then”(Y)X=”: X◢“(Y)Y=”：Y◢F-90→F:“F WJ=”：F▼DMS◢(5800里面的度分秒的)Prog “123”：“BG =>Z”:D◢Goto4: IfEnd↙If Z=-1: Then”(Z)X=”: X◢“(Z)Y=”：Y◢F-90→F:“F WJ=”：F▼DMS◢(5800里面的度分秒的)Prog “123”：“BG =>Z”:D◢Goto4: IfEnd↙If Z＞1: Then”(Y)X=”: X◢“(Y)Y=”：Y◢F-90→F:“F WJ=”：F▼DMS◢(5800里面的度分秒的)Prog “123”：“BG =>Z”:D◢Goto4: IfEnd↙If Z＜1: Then”(Z)X=”: X◢“(Z)Y=”：Y◢F-90→F:“F WJ=”：F▼DMS◢(5800里面的度分秒的)Prog “123”：“BG =>Z”:D◢Goto4: IfEnd↙Lbl 2：“(Q)X =”？X:“(Q)Y=” ?Y:X →I:Y →J: Prog"1.3-ZFS"：O+W→S：“(SQ)ZH=”：S◢“(SQ)Z=”：Z◢Goto 4↙Lbl 1:If C=1989:Then Goto 4:0→C:Else Stop:IfEnd↙(密码符号段把原来的密码调到后面去了 )2. 正算子程序(1.2-ZFS )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→F3. 反算子程序(1.3-ZFS)G-90→T：Abs((Y-V)cos(T)-(X-U)sin(T)) →W:0→Z：Lbl 0：Prog "1.2-ZFS"：T+QEW(C+WD)→L：(J-Y)cos（L）-(I-X)sin（L）→Z：If Abs(Z)< -6：Then Goto1: Else W+Z→W:Goto 0↙Lb1 1:0→Z:Prog"1.2-ZFS":(J-Y)÷sin (F)→Z↙4. 数据库子程序(1.1-ZFS)If S＜***:Then Cls:Stop:IfEnd(全线路的起点里程桩号如：If S＜1497534.054:Then Cls:Stop:IfEnd小于1497534.054的全部失败)If S〉*****:Then Cls:Stop:IfEnd↙（同上）If S≦***（线元终点里程）:Then***→U（线元起点X）:***→V（线元起点Y）: ***→O（线元起点里程）:***→G（线元起点方位角）: ***→H（线元起点至终点长度）：***→P（线元起点曲率半径,在这里声明下，如果是直线段的话输1.×10+45就是输1. +45）:***→R（线元终点曲率半径同上）: 0或1、-1→Q：Goto 1 :IfEnd ↙If S≦……………………………………………：Goto 1 :IfEnd ↙………………[注意：输入方法：ZH-HY(ZH点注意P输入1.×10+45就是输1.+45，HY R按照曲线要素半径),HY-YH(P和R都按照曲线要素输入半径),YH-HZ(YH P按照曲线要素半径，R输1.×10+45),HZ-ZH(P和R 都输1.×10+45) Q为转角值看曲线上是左转角-1还是右转角1HZ-ZH Q为0]最后结尾记得一定要输入Lbl 1:1÷P→C:(P-R)÷(2HPR)→D:180÷π→E：Return↙5.BG JS 20M (子程序)Deg; Fix3↙Lbl1:”ZH=”?:Prog“BG SJK”:”Z=”?Z↙If Z=0:Then A→D:”(0)BG =>Z”:D◢Goto1: IfEnd↙If Z=1:Then A→D:”(Y)BG =>Z”:D◢Goto1: IfEnd↙If Z=-1:Then A→D:”(Z)BG =>Z”:D◢Goto1: IfEnd↙If Z＞1:Then A+(Z-1)(C/100)→D:”(Y)BG =>Z”:D◢Goto1: IfEnd↙If Z＜-1:Then A+(-Z-1)(-B/100)→D:”(Z)BG =>Z”:D◢Goto1: IfEnd↙6.BG SJK(标高数据库)If S =K: Then H→A:i(Z)→B: i(Y)→C: IfEnd↙注：K 为桩号里程，H为中桩高程，i(Z)为左横坡，i(Y)为右横坡。

卡西欧5800公路坐标正反算程序卡西欧5800公路坐标正反算程序是一种用于测量和计算公路坐标的工具。

它基于卡西欧5800全站仪的技术和功能，能够准确地测量和计算公路的坐标信息。

下面我将详细介绍卡西欧5800公路坐标正反算程序的工作原理和使用方法。

一、工作原理：卡西欧5800全站仪是一种高精度的测量仪器，它能够通过测量角度和距离来确定地点的坐标。

在公路测量中，我们通常会在公路上设置一系列的控制点，然后使用全站仪测量这些控制点的坐标。

通过这些测量数据，我们可以计算出公路上其他点的坐标。

卡西欧5800公路坐标正反算程序基于这些测量数据，通过一系列的计算和算法，能够准确地计算出公路上其他点的坐标。

具体的计算方法包括正算和反算。

正算是指已知控制点的坐标，通过测量数据和算法计算出其他点的坐标。

反算是指已知某个点的坐标，通过测量数据和算法计算出该点在公路上的位置。

二、使用方法：1. 设置控制点：首先，在公路上选择一些具有代表性的控制点，并使用全站仪进行测量。

在测量过程中，需要确保全站仪的水平和垂直仪器准确，以获得准确的测量数据。

2. 导入测量数据：将测量数据导入卡西欧5800公路坐标正反算程序。

该程序通常会提供一个界面，用于导入和管理测量数据。

3. 正算：如果需要计算公路上其他点的坐标，可以选择正算功能。

在正算功能中，需要输入已知控制点的坐标和测量数据。

程序将根据这些数据进行计算，并给出其他点的坐标。

4. 反算：如果需要计算某个点在公路上的位置，可以选择反算功能。

在反算功能中，需要输入已知点的坐标和测量数据。

程序将根据这些数据进行计算，并给出该点在公路上的位置。

5. 结果输出：计算完成后，卡西欧5800公路坐标正反算程序将给出计算结果。

结果通常以表格或图形的形式呈现，以便用户查看和使用。

三、注意事项：1. 测量准确性：为了获得准确的测量结果，使用卡西欧5800全站仪时需要注意测量准确性。

在测量过程中，需要保持仪器的稳定和准确，并避免测量误差。

卡西欧fx-5800 坐标正算程序“SO”算方位Lbl 0“D”：Pol(X-A，Y-B) →D◢If J≥0:Then J◢Else “J”:J+360◢IfEnd“ZZ”正算直线坐标A+Rec(Abs(P-K),Q) →XB+J→Y“X”：X+Rec(D,Q+N×90) →X◢“Y”:Y+J→Y◢P?所求桩号 A，B：起点x,y K:起点桩号 Q：方位角D：边桩距离 N？左偏“-1”右偏：“+1”X= Y= 所求点坐标“YY”正算圆曲线坐标Q+C×（P-K）÷S×180÷∏→Z“X”：A+C×S×（Sin(Z)-Sin(Q)）+Rec(D,Z+N×90) →X◢“Y”:B-C×S×（Cos(Z)-Cos(Q)）+J→Y◢P:所求点桩号 A，B：起点x,y K：起点桩号Q：方位角 S：曲线半径 D？边桩距离N？左偏“-1”，右偏“+1”X= y= 所算点坐标 C：线路左转“-1”右转“+1”“H1”正算缓和曲线坐标Abs(P-K) →LT2÷L→HL2×180÷∏÷2÷T2→VQ+C×V→ZL-L∧（3）÷40÷H2+L∧（5）÷3456÷H∧（4）-L∧（7）÷H∧（6）÷599040→EL2÷6÷H- L∧（4）÷336÷H∧（3）+ L∧（6）÷42240÷H∧（5）- L∧（8）÷H∧（7）÷（6×1612800）→FIf C<0：Then-F→F：IfEndA+E×Cos(Q)-F×Sin(Q) →X“X”：X+Rec(D,Z+M×N×90) →X◢B+E×Sin(Q)+F×Cos(Q) →Y“Y”：Y+J◢P？所求点桩号 A，B：起点x,y K:起算点里程桩号Q：前切线方位角 T（是圆曲线要素的A1）=√（根号）缓曲总长×半径C=切线线路方向左转“-1”，右转“+1“M=第一缓曲“+1”，第二缓曲”-1“N：算边桩左偏“-1”，右偏“+1”D：边桩距离坐标正算数据库例子“肇兴隧道右线数据库”Lbl 0“P”?P:“N”？N：“D”？DIf P≥16836.5 And P≤17009.624:Then 66881.803→A：14635.782→B：16836.5→K:280032030.50→Q:902.044→T:1→C: 1→M:Prog“H1”：IfEnd (第一缓曲线)If P≥17009.624 And P≤17529.681:Then 66914.52→A：14465.781→B：17009.624→K:281035049.370→Q:4700→S：Prog“YY”:IfEnd （圆曲线）If P≥17539.681 And 17702.805≥P:Then 67102.343→A：13799.117→B：17702.805→K:108059031.50→Q:902.044→T:-1→C: -1→M:Prog“H1”：IfEnd (第二缓曲线)If P≥17702.805 And P≤19550:Then 67102.343→A：13799.17→B：17702.805→K:288059031.50→Q: Prog“ZZ”：IfEnd （直线）。

卡西欧5800计算器坐标正反算程序主程序坐标正反算程序名称：ZBZFSLb1 A↙30→Dim Z :"0=ZS,1=FS"?Z ↙If Z=0:Then Goto B:IfEnd↙ (Z=0进入里程点坐标正算）If Z=1:Then Prog"FSLCZ "↙ (Z=1进入反算里程边桩）Lb1 B ↙Prog”ZBZS” ↙子程序反算里程桩名称：FSLCZLb1 1↙"XK="?X:"YK="?Y↙（输入任意测点的X、Y坐标）“K=”?K ↙（试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）Lb1 2↙Prog"SJK ":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos(W)-(X-Z[18])*Sin(W))→S↙If S<0.0001:Then Goto 4:Else Goto 3:Ifend↙Lb1 3↙K+S→K: Prog"SJK":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z[18])*Sin (W)→Q↙If Q<0.0001 :Then Q→S: Goto 4:Else if Q<S: Then K+Q→K:Goto 2 :Else if Q>S :Then K-Q→K:Goto 2:Ifend:Ifend:Ifend↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:"DP(-Z+Y)=":I◢(偏距）"K=":K+S→K◢（里程）Goto 1↙子程序坐标正算名称：ZBZSLb1 0 ↙“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=":JDMS◢(后视方位角)Lb1 1↙“K=”?K ：（输入所需计算里程）Prog"SJK "↙XI : Z[18]↙YI: Z[19] ↙Pol(Z[18]-M,Z[19]-N): J<0=>J+360→J↙“PJ=”?P↙ (输入桩与线路夹角)“PD=”?D↙ （输入桩距中线的距离）Z[18]+D*Cos(T+P) →Z[20] ↙Z[19]+D*Sin(T+P) →Z[21] ↙“X=”: Z[20] ◢ （放样坐标X）“Y=”: Z[21] ◢ （放样坐标Y）Pol(Z[20]-M,Z[21]-N):"D=":I◢（放样距）J<0=>J+360→J:"F=":JDMS◢（放样方位角）Goto 1↙子程序数据库名称：SJKif K<本段曲线终点里程And K≥上段曲线终点里程：Then 本段曲线终点里程→Z[1] : 上段曲线终点里程（第一段曲线输起点的里程）→Z[2] ：1→O （注：左偏曲线输入-1→O,右偏曲线输入1→O）: 半径→R :曲线偏角→A:第一缓和曲线→Z[6] : 第二缓和曲线→Z[7] : 交点X→B :交点Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : Pr og”JSPB”:Return: Ifend↙if…………Prog”JSPB”:Retur n:Ifend（曲线段分段输入）↙补充直线段输入如下 (只需输线路的最后一段直线数据）if K<本段直线终点里程And K≥本段直线起点里程:Then 本段直线终点里程→Z[3]:终点X→Z[16]:终点Y→Z[17]:方位角→E: Z[16]+ (K- Z [3])*Cos(E)→Z[18]:Z[17]+ (K- Z[3])*Sin(E)→Z[19] : Return:Ifend ↙子程序计算判别名称：JSPBLb1 2 ↙（曲线要素计算）Z[6]/2- Z[6]^3/(240*R^2)+ Z[6]^5/(34560*R^4) →Z[8] ↙ （M1）Z[7]/2- Z[7]^3/(240*R^2)+ Z[7]^5/(34560*R^4) →Z[9] ↙ （M2）Z[6]^2/(24*R)- Z[6]^4/(2688*R^3) →Z[10] ↙ （P1）Z[7]^2/(24*R)- Z[7]^4/(2688*R^3) →Z[11] ↙ （P2）π*A*R/180+0.5*( Z[6]+ Z[7])→Z[25] ↙ （曲线总长）90* Z[6]/(R*π) →Z[14] ↙ （第一缓和曲线总偏角）90* 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]- Z[25]→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] And K> Z[3] : Then Goto 5 : Ifend ↙if K≤Z[5] And K> Z[4] : Then Goto 6 : Ifend ↙if K≤Z[1] And K> Z[5] : Then Goto 7 : Ifend ↙Lb1 4 ↙（里程小于直缓点直线独立坐标）K- Z[3] →Z[23] : 0→Z[24] : E→T : Goto 8↙Lb1 5 ↙（第一缓和曲线独立坐标）K- Z[3] →H ↙H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →Z[23] ↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Z[24] ↙90*H^2/( R*π* Z[6]) →T ↙if O>0 ：Then T +E→T : Eles E-T →T : T<0=>360+T→T : Ifend ↙Got o 8 ↙Lb1 6 ↙（圆曲线独立坐标）K- Z[4] →H ↙H*180/( R*π)+ Z[14]→T ↙R*Sin(T)+ Z[8]→Z[23] ↙R*(1-Cos(T))+ Z[10]→Z[24] ↙if O>0 ：Then T +E→T : Eles E-T →T : T<0=>360+T→T : Ifend ↙Goto 8 ↙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)→Z[23] ↙Z[13]*Sin(A)-U*Sin(A)+V*Cos(A)→Z[24] ↙if O>0 ：Then F-T→T : T<0=>360+T→T : Else F+T →T : Ifend ↙Goto 8 ↙Lb1 8if O<0 : Then - Z[24]→Z[24] : Ifend ↙Z[16]+Z[23]*Cos(E) -Z[24]*Sin（E）→Z[18] ↙Z[17]+Z[23]*Sin(E)+Z[24]*Cos(E)→Z[19] ↙Return↙此程序大约占计算器2900字节！此程序本人已亲自验算，无误！！子程序反算里程桩名称：FSLCZLb1 1↙"XK="?X:"YK="?Y↙（输入任意测点的X、Y坐标）“K=”?K ↙（试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）Lb1 2↙Prog"SJK ":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos(W)-(X-Z[18])* Sin(W))→S↙If S<0.0001:Then Goto 4:Else Goto 3:Ifend↙Lb1 3↙K+S→K: Prog"SJK":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z [18])*Si n (W)→Q↙If Q<0.0001 :Then Q→S: Goto 4:Else if Q<S: Then K+Q →K: Goto 2 :Else if Q>S :Then K-Q→K:Goto 2:Ifend:Ifen d:Ifend↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:"DP(-Z+Y)=":I◢(偏距）"K=":K+S→K◢（里程）Goto 1↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:------------差个撒样"DP(-Z+Y)=":I◢(偏距）"K=":K+S→K◢（里程）Goto 1↙。

目录一、坐标正算基本公式………………………………………………………０２二、坐标反算原理 (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 。

卡西欧5800计算器程序坐标正反算（含高程）功能:坐标正反算(含高程),把要素内置化(无需改程序文件,即可更换路线,同时存两条线要素)另含两个边仰(一个横向坡,一个纵向坡)坡放样模块.帮助文件未完成程序浏览:FileName:RESET 初始化程序Norm 1:50→C:12345→J"RESET PW"I:I=J=>500→DimZ 为数据库增加额外变量500个，在SET、SETPFDYS程序根据实际再增减变量50→Z[C+22]"PASSWORDS"I:I→Z[C+39] 重设要素保护密码Cls:StopFileName:DATLOCK 要素保护密码确定认程序Cls:Norm 1:50→C"PASSWORDS"I:Cls:I≠Z[C+39]=>StopFileName:SHELL(外壳程序)50→C ;在扩充变量预留前50个给别的程序用。

如不够就适量加大。

RESET，SHELL，SET，SETPFDYS，这几个程序中C值必需一致Z[C+35]→I:"STATION-N"I:I→Z[C+35] 设置测站N坐标Z[C+36]→I:"STATION-E"I:I→Z[C+36] 设置测站E坐标Z[C+44]→I:"STATION-Z"I:I→Z[C+44] 设置测站Z坐标Z[C+45]→I:"STATION-HI"I:I→Z[C+45] 设置仪高Z[C+41]→I:"GC-DH"I:I→Z[C+41] 放样点高差常数Norm 1Z[C+23]→N 当前分段要素N坐标Z[C+24]→E 当前分段要素E坐标Z[C+25]→M 当前分段要素起点桩号Z[C+26]→H 当前分段要素起点方位角（正北）单位：弧度Z[C+27]→A 当前分段要素起点曲率有左偏负右偏正（注意不半径）Z[C+28]→R 当前分段要素终点曲率有左偏负右偏正（注意不半径）Z[C+29]→L 当前分段要素长度Z[C+3]→WZ[C+32]→I:"DAT1 2 3"I:I→Z[C+32] 平曲线要素数库类型选择1为内置式，2文件式，3实时输入(查看当前要素值)0→IZ[C+4]→I:"1 2 3"I:I→Z[C+4] 选择本程序模式默认为坐标正算，1坐标反算，2横向边仰坡放样，3，纵向边仰坡放样(隧道进出口用到) Lbi 1Z[C+4]→II=1=>Prog"TURNZH"I=2=>Prog"HXBYP"I=3=>Prog"ZXBYP"Lbi ADeg:Norm 1:ClsZ[C+1]→G:"ZH"G:G→Z[C+1] ;G桩求桩号Z[C+2]→B:"JL"B:B→Z[C+2] ;B横向距离，左正右负Lbi BProg "ZBJS" ;坐标正算Prog "GCJS" ;高程计算Fix 3:Cls ;设置三位小数"ZH＝":Locate 4,1,G ;第一行显示桩号"X＝":Locate 3,2,X ;第二行显示X(N)坐标"Y＝":Locate 3,3,Y ；第三行显示Y(E)坐标Prog "GCJS" ;高程计算"Z":Locate 2,4,Z+Z[C+41] ；第四行显示Z坐标Locate 10,4,Z[C+2] ；第四行显示横向距离Do:I+1→I:I=1000=>Goto 1:LpWhile Not(Getkey=57 Or Getkey=27) ;锁定键盘，并在几十秒后自动返回Getkey=57=> Goto 1Z+Z[C+41]+Z[C+49]-Z[C+44]-Z[C+45]→Z[C+48]Pol(Z[C+35]-X,Z[C+36]-Y)Cls:"DL":Locate 3,1,I:Locate 12,1,Z[C+45]"DH":Locate 3,2,Z[C+48]:Locate 12,2,Z[C+49]"FWJ"J<0=J+360→J:J◤DMS◢J→Z[C+46]I→Z[C+47]Goto 1FileName：ZBJS 坐标计算程序Prog"READDAT"RadG-M→QIF AR=0 :Then If A=R :Then 1→J:Else 3→J :IfEnd:Else If A=R :Then 2→J:Else 3→J:IfEnd:IfEndIf J=1 :Then H→F:N+QCos(H)→X:E+QSin(H)→Y:IfEnd 直线段直接计算If J=2 :Then H+QR→F:Rec(Abs(2Sin(Abs)÷R),H+QR÷2):N+I→X:E+J→Y:IfEnd 圆弧段直接计算If J=3 :Then Goto 5 :IfEndX+BCos(F+W)→XY+BSin(F+W)→YDeg:ReturnLbi 5 用五点通用坐标计算计算缓和段(R-A)÷L→KKQ2→J→P:H+IP+JP2→U→P:H+IP+JP2→V→P :H+IP+JP2→D→P:H+IP+JP2→F→P:H+IP+JP2→T→I→J→ON+Q(ICos(U)+JCos(V)+OCos(D)+JCos(F)+ICos(T))→XE+Q(ISin(U)+JSin(V)+OSin(D)+JSin(F)+ISi n(T))→YH+AQ+KQ2→FX+BCos(F+W)→XY+BSin(F+W)→YDeg:ReturnFileName:TURNZH(坐标反算)Norm 1Z[C+37]→I:"DQD-N"I:I→Z[C+37] 输入待求点N坐标I=-1=>Prog"INFWJDLDH"Z[C+38]→I:"DQD-E"I:I→Z[C+38] 输入待求点E坐标0→B：M+L÷2→G：Prog"ZBJS"Lbi S:RadZ[C+37]-X→I:Z[C+38]-Y→JIf I=0 And J=0 :Then Goto A:IfEndPol(I,J)Lbi A:Rec(I,J-F):J→Z[C+2]G+I→G:IF Abs(I)> :Then Prog"ZBJS":Goto S↙G→Z[C+1]FileName：INFWJDLDH 实测坐标（用方位角，距离，高差）输入程序Deg:Norm 1:ClsZ[C+46]→J:"DQD-FWJ"J:J→Z[C+46]Z[C+47]→I:"DQD-DL"I:I→Z[C+47]Z[C+48]→K:"DQD-DZ"K:K→Z[C+48]Z[C+49]→P:"RHT"P:P→Z[C+49]Rec(I,J)Z[C+35]+I→Z[C+37]Z[C+36]+J→Z[C+38]Z[C+44]+Z[C+45]+K-P→Z[C+43]FileName:HXBYP 横向边仰坡放样程序Cls:Norm 1Z[C+40]→I:"QPD-DL"I:I→Z[C+40] 输入起坡点与中桩距离常数Z[C+41]→I:"QPD-DZ"I:I→Z[C+41] 输入起坡点与中桩高差常数Z[C+42]→I:"i"I:I→Z[C+42] 输入边仰坡坡度，左仰坡(路堑)为正，右仰坡(路堑)为负，左边坡为负，右边坡为正Cls:Z[C+43]→I:"DQD-Z"I:I→Z[C+43] 输入实测高程Prog"TURNZH"Prog "GCJS"(Z[C+43]-Z-Z[C+41])Z[C+42]+Z[C+40]→BFix 3:B-Z[C+2]◢显示与上一次测点的偏差B→Z[C+2]FileName:ZXBYP 纵向仰坡放样程序(隧道进口使用)Cls:Norm 1:Z[C+33]→I:"QPD-ZH"I:I→Z[C+33] 输入起坡点桩号Z[C+34]→I:"QPD-Z"I:I→Z[C+34] 输入起坡点高程Z[C+42]→I:"i"I:I→Z[C+42] 仰坡时(隧道进出口仰坡)，进口为正，出口为负)Cls:Z[C+43]→I:"DQD-Z"I:I→Z[C+43] 输入实测高程Prog"TURNZH"Z[C+42](Z[C+43]-Z[C+34])+Z[C+33]→GFix 3:G-Z[C+1]◢显示与上一次测点的偏差G→Z[C+1]FileName:SET 设置程序Lbi SNorm 150→C50→Z[C+22]0→KZ[C+5]→I:"1 Or 2"I:I→Z[C+5] ;选择1线或者2线180Z[C+3]÷π→J:"XZJJ"J:Jπ÷180→Z[C+3] ;斜桩夹角If I=1:Then Z[C+8]→Z[C+6]:Z[C+9]→Z[C+7]Z[C+12]→Z[C+16]:Z[C+13]→Z[C+17]:Z[C+18]→Z[C+20]Else Z[C+10]→Z[C+6]:Z[C+11]→Z[C+7]Z[C+14]→Z[C+16]:Z[C+15]→Z[C+17]::Z[C+19]→Z[C+20]Z[C+7]+3→Z[C+31]0→I:"1PQX,2SQX,3YSDS"I ;1输入平曲线要素，2输入竖曲线要素，3输入1、2线夹平竖曲线交点(变坡点)数量以让程序规划内存建立要素数据为，其它值退出本程序I=1 => Goto 1:I=2 => Goto 2:I=3 => Goto 3I≠0=>Prog"DATLOCK"Prog"SETPFDYS"StopLbi 1 ;平曲线要素输入部分Z[C+6]→D-1→Z[D+3]For 1→J T o Z[C+16]Cls:"PQX":Locate 9,1,JZ[D+1]→I:"N"I:I→Z[D+1] ;请输入N坐标Z[D+2]→I:"E"I:I→Z[D+2] ;请输入E坐标Z[D+3]→I:"R"I:I→Z[D+3] ;请输入半径RZ[D+4]→I:"LS1"I:I→Z[D+4] ;请输入缓和曲线1长度Z[D+5]→I:"LS2"I:I→Z[D+5] ;请输入缓和曲线2长度Next1→Z[D-2]Goto SLbi 2 ;竖曲线要素输入部分Z[C+7]→D-1→Z[D+3]For 1→J T o Z[C+17]Cls:"SQX":Locate 9,1,JZ[D+1]→I:"BPDZH"I:I→Z[C+1] ;请输入变坡点桩号Z[D+2]→I:"BPDGC"I:I→Z[C+2] ;请输入变坡点高程Z[D+3]→I:"R"I:I→Z[C+3] ;请输入半径RD+3→DNext1→Z[D]Goto SLbi 3 ;要素点数输入部分Z[C+12]→I:"1PQXDS"I:I→Z[C+12] ;输入1线平曲线点数,最小值为3Z[C+13]→I:"1SQXDS"I:I→Z[C+13] ;输入1线竖曲线点数,最小值为3Z[C+18]→I:"1JDZH"I:I→Z[C+18] ;输入1线平曲线第一个交点的交点桩号Z[C+14]→I:"2PQXDS"I:I→Z[C+14] ;输入2线平曲线点数,最小值为3Z[C+15]→I:"2SQXDS"I:I→Z[C+15] ;输入2线竖曲线点数,最小值为3Z[C+19]→I:"2JDZH"I:I→Z[C+19] ;输入2线平曲线第一个交点的交点桩号C+Z[C+22]→Z[C+8] 1线路平曲线要素指针偏移基数Z[C+8]+5Z[C+12]→Z[C+9] 1线路竖曲线要素指针偏移基数Z[C+9]+3Z[C+13]→Z[C+10] 2线路平曲线要素指针偏移基数Z[C+10]+5Z[C+14]→Z[C+11] 2线路竖曲线要素指针偏移基数Z[C+11]+3Z[C+15]→Z[C+21] 要素数据库结束指针Z[C+21]+1→Z[C+30]Z[C+21]→DimZ 增加额外变量-1→KGoto SFileName:SETPFDYS 交点要素=>分段要素Norm 150→CRadZ[C+6]+5→DZ[C+21]+1→FF+6→DimZ-1→Z[F] ; 设置分段要素开始标志Z[C+20]→V ; 从数据库读入当前线第一个平曲线交点桩号0→GPol(Z[D+1]-Z[D-4],Z[D+2]-Z[D-3]) ；计算第一条直线的距离I，计算方位角JJ→HZ[D-4]→Z[F+1] 直线段Z[D-3]→Z[F+2]V-I→Z[F+3]For 1→G T o Z[C+16]-2Pol(Z[D+1]-Z[D-4],Z[D+2]-Z[D-3]) ；计算第本交点与下一交点的距离I，计算方位角JI→W:J→HPol(Z[D+6]-Z[D+1],Z[D+7]-Z[D+2])J-H→K ；计算转角K，负值是左转角，正值是右转角IF Abs(K)>π:Then If K>0:Then K-2π→K:Else K+2π→K:IfEnd:IfEndZ[D+3]→R ;从数据库读入半径RZ[D+4]→L ;从数据库读入缓和曲线1长度Z[D+5]→M ;从数据库读入缓和曲线2长度L÷2÷R→B ;计算缓和曲线1的B0角M÷2÷R→O ;计算缓和曲线2的B0角(Abs(K)-B-O)R→P ;计算圆曲线长度L2÷24÷R-L^(4) ÷2384÷R^(3) →Q ;内移值计算方法2L÷2-L^(3) ÷240÷R2→S ;切线增长值计算方法2M2÷24÷R-M^(4) ÷2384÷R^(3) →J ;内移值计算方法2M÷2-M^(3) ÷240÷R2→I ;切线增长值计算方法2Abs(π÷2- Abs(K)) →ARtan(Abs(K÷2))+S+J÷Cos(A)-QTan(A)→T ;计算切线1Rtan(Abs(K÷2))+I+Q÷Cos(A)-JTan(A)→U ;计算切线2H<0=>H+2π→HK<0 =>-1→AH→Z[F+4]0→Z[F+5]0→Z[F+6]F+6→FF+6→DimZRec(T,H) 缓和段1 Z[D+1]-I→NZ[D+2]-J→EN→Z[F+1]E→Z[F+2]V-T→Z[F+3]L=0=>Goto AH→Z[F+4]0→Z[F+5]A÷R→Z[F+6]F+6→DimZLbi A 圆弧段0→I:0→JL=0=>Goto BPol(L-L^(3)÷40÷R2,L2÷6÷R-L^(4)÷336÷R^(4)) Lbi B Rec(I,H+AJ)N+I→NE+J→EN→Z[F+1]E→Z[F+2]V-T+L→Z[F+3]H+AB→Z[F+4]A÷R→Z[F+5]A÷R→Z[F+6]F+6→FF+6→DimZM=0=>Goto C 缓和段2 Rec(2RSin÷R),H+AB+AP÷R÷2) N+I→NE+J→EN→Z[F+1]E→Z[F+2]V-T+L+P→Z[F+3]H+AB+AP÷R→Z[F+4]A÷R→Z[F+5]0→Z[F+6]Lbi CF+6→FF+6→DimZRec(U,H+K)Z[D+1]+J→Z[F+1]Z[D+2]+I→Z[F+2]V-T+L+P+M→Z[F+3]V-T+L+P+M-U+I→VD+5→DNextH+K→Z[F+4]0→Z[F+5]0→Z[F+6]F+9→DimZ1→Z[F+7]V→Z[F+9] 要素终点桩号DegFileName：READDATZ[C+32]→II=3=>Prog"READDAT3"If GM+L :Then ;验证要素是否可用，否则重设要素I=1=>Prog"READDAT1"I=2=>Prog"READDAT2"Prog"REALDAT3"IfEndFileName:READDAT1 平曲线分段要素读取程序1Z[C+30]→F0→JLbi 1 搜索要素If G<z[f+3]:then< p="">-1=Z[F]=>Goto E:F-6→F:-1→J:Goto 1Else 1=Z[F+1]=>Goto E:-1=J=>Goto 2:1→J:F+6→F:Goto 1 IfEndLbi 2F→Z[C+30]Z[F+1]→Z[C+23] 当前分段要素N坐标Z[F+2]→Z[C+24] 当前分段要素E坐标Z[F+3]→Z[C+25] 当前分段要素起点桩号Z[F+4]→[C+26] 当前分段要素起点方位角（正北）Z[F+5]→Z[C+27] 当前分段要素起点曲率有左偏负右偏正（注意不半径）Z[F+6]→Z[C+28] 当前分段要素终点曲率有左偏负右偏正（注意不半径）</z[f+3]:then<>。

卡西欧5800坐标正反算程序卡西欧5800坐标正反算程序1.主程序(TYQXJS)Lb1 4:“1.SZ=>XY”：“2.XY=>SZ”:？N:？S:Prog“SUB0”↙1÷P→C：（P-R）÷（2HPR）→D：180÷π→E:N=1=>Goto 1:Goto 2: ↙Lb1 1:?Z:Abs(S-O) →W: Prog“SUB1”:“XS=”：X◢“YS=”:Y◢F-90→F:“FS=”F▲DMS◢ Goto 4↙Lb1 2:?X:?Y:X→I:Y→J:Prog“SUB2”：O+W→S:“S=”：S◢“Z=”:Z◢Goto 4↙2.正算子程序（SUB1）0.1739274226→A:0.3260725774→B:0.0694318442→K:0.330 0094782→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+MW D)))→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. 反算子程序（SUB2）G-90→T:Abs((Y-V)cos(T)-(X-U)sin(T)) →W:0→Z:Lb1 0:Prog“SUB1”:T+QEW(C+WD) →L:(J-Ycos(L)-(I-X)s in(L) →Z:If Abs(Z)〈0.000001：Then 0→Z:Prog“SUB1”:(J-S)÷sin(F) →Z:Else W+Z→W:Goto 0:IfEnd↙4. 数据库子程序（SUB0）Goto1↙同时保存多个曲线时的指针Lb1 1:If S〈 000.000（线元终点里程）：Then 000.000→G(线元起点方位角)：000.000→O （线元起点里程）：000.000→U(线元起点X坐标)：000.000→Y（线元起点Y坐标）：000.000→P(线元起点曲率半径)：000.000→R（线元终点曲率半径）：000.000→H（线元起点至终点长度）：0或1、-1→Q：Eeturn:IfEnd↙Lb1 1:If S〈 000.000（线元终点里程）：Then 000.000→G(线元起点方位角)：000.000→O （线元起点里程）：000.000→U(线元起点X坐标)：000.000→Y（线元起点Y坐标）：000.000→P(线元起点曲率半径)：000.000→R（线元终点曲率半径）：000.000→H（线元起点至终点长度）：0或1、-1→Q：Eeturn:IfEnd↙………………….为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lb1，每增加一个工程增加一个文件。

CASIO 5800计算器程序一、程序功能：本程序由一个主程序（AK）及两个子程序：坐标正算程序（SUB1）、坐标反算程序（SUB2）平曲线数据库（PQX）构成，可以根据曲线段的组成部分：直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、止点坐标、起点切线方位角、线元长度、起点曲率半径、止点曲率半径、里程、边距或坐标），对该曲线段范围内任意里程中边桩坐标进行正反算。

二、源程序：1、主程序（AK）A排在前面，方便运行程序，不用往下翻页了。

Deg：“0→ZX,1→A,2→B”?Z ：“1→JDF,2→XYF”?E：“X(CZ)”?A：“Y”?B 说明：设置角度模式，选择线路0为主线，1为A匝道。

可以自己再增加。

选择计算方法，交点法还是线元法计算坐标，1为交点法，2为线元法。

输入测站坐标，如不需要计算极坐标，即放样的时候需要输入坐标到仪器里面去，则X(CZ)输0，Y不用管可以为任何值，建议把测站坐标输进去，省了放一个点输一次坐标数据到仪器里面去。

只有一条线路可以不需要黄色部分。

“1.Ki→XY,2.XY→Ki”?N：Lbl 0：“Ki”?K：Prog“PQX”IfE=1：Then Prog“JD”：IfEndN=1＝＞Goto 1：N=2＝＞Goto 2Lbl 1?LK-O→W：Prog“SUB1”说明：1为正算，即知道桩号和偏距（注意只能计算边桩，不能计算斜交的点，构造物可以通过后面的ZHIXIAN程序计算），2为反算，即知道坐标，求桩号和偏距。

Ki为桩号，正算直接输，反算需要输入附近桩号，第一次计算最好输好桩号之后，把反算出来的桩号再输进去，重新反算下，如果和第一次反算的桩号一样，那就说明反算的桩号正确。

调用平曲线数据库，输入偏距L（偏距左输负数，右输正数），如果选择的是交点法计算，那需要输入是数据库是交点数据，调用正算子程序SUB1。

If A≠0：Then Pol（X-A,Y-B）：“JULI=”：I◢J<0=>J+360→J：“FWJ=”：J ▶DMS◢Else “Xi=”：X◢“Yi=”：Y◢“FWJ=”：F ▶DMS◢If End ：Goto0说明：正算，如果前面X(CZ)输0，则直接显示X、Y坐标和方位角，否则显示仪器到放样点的距离和方位角。

卡西欧5800道路坐标测量程序（通俗易懂版）主程序（自己给取个名吧）2→Dimz“1.SZ＝﹥XY2.XY＝﹥SZ”◢“N=”?NIf N=1:Then Goto 1:Else Goto 2:IfEndLbl 1“S=”?S“Z=”?Z“J=”?TIf S﹤0:Then Prog “SYL”:Else Prog “SYR”:IfEnd(可拿掉)1÷P→C:(P-R) ÷(2HPR)→D:180÷∏→E:Abs(Abs(S)-O)→W:Prog “SUB1”“QJ=”:Z[1]◣DMS◢“X=”:X◢“Y=”:Y◢Goto 1Lbl 2“S=”?SIf S<0:Then Prog “SYL”:Else Prog “SYR”:IfEnd(可拿掉) 90→T1÷P→C:(P-R) ÷(2HPR)→D:180÷∏→E“X=”?X“Y=”?YX→I:Y→JProg “SUB2”O+W→S（红色加粗加下画线为字母，常规为零）“S=”:S◢Z→Z“Z=”:Z◢Goto 2子程序（SUB1）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)))→XV+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+F WD))+Asin(G+QEMW(C+MWD))) →YG+QEW(C+WD)→Z[1]Z[1]+T→FX+Zcos(F)→XY+Zsin(F)→YReturn子程序(SUB2)G-90→Z[2]:Abs((Y-V)cos(Z[2])-(X-U)sin(L)→W:0→ZLbl 0Prog “SUB1”Z[2]+QEW(C+WD)→L(J-Y)cos(L)-(I-X)sin(L)→ZIf Abs(Z)<10(-6):Then Goto 1:EIse W+Z→W:Goto 0:IfEndLb1 10→Z:Prog “SUB1”(J-Y)÷sin(F)→ZReturn数据库（SYR）If S>29599.031(本段线元计算范围起点里程) And S≤29759.031(本段线元计算范围终点里程):Then 3311072.206(本段线元起点X坐标)→U:502909.1014(本段线元起点Y坐标)→V: 29599.031(本段线元起点里程)→O:256-25-59.11(本段线元起点方位角)→G:160(本段线元长度)→H:1400(本段线元起点半径)→P:10(45)(本段线元终点半径)(直线段为10的45次方)→R:-1(线元转向标识)(左偏为-1,右偏为1,直线为0)→Q:IfEnd输入说明正算 1 输入 N:输入1进入计算坐标模式S:里程桩号Z:偏距J:与线路角度2 显示 QJ:切线方位角X:X坐标Y:Y坐标反算 1 输入 N:输入2进入坐标反算里程偏距模式S:输入相近里程X:输入X坐标Y:输入Y坐标2 显示 S:反算出的里程桩号Z:反算出的偏距。

CASIO fx—5800计算器工程测量与公路测量程序1、坐标正算”X0=〃X:”Y0=〃Y:”I=〃I:”J=〃JX+ICos(J)→U:Y+ISin(J)→V”X=〃:U （待求点的X坐标）”Y=〃:V （待求点的Y坐标）说明：X0 Y0：已知点坐标I：两点的距离J：方位角2、坐标反算Lbl 0”X1=〃X:”Y1=〃Y:”X2=〃U:”Y2=〃VPol(U-X,V-Y):J<0 J+360→J”I=〃:I”J=〃:J DMSGOTO 0说明：X1Y1：第一点的坐标，X2Y2第二点的坐标，I：两点的距离，J：方位角3、圆曲线”X0=〃X: ”Y0=〃Y:”FWJ=〃D:”QDHAO=〃G: ”ZDHAO=〃C: ”R=〃RLbl 1”DQHAO=〃LL<G Or L＞C GOTO 2”PY=〃K180(L-G)÷(πR) →E2RSin →FReC(Abs(F),D+：Cls”X=〃:X+I+KCos(D+E+90) →A”Y=〃:Y+J+KSin(D+E+90) →BD+E→o:o<0 360+o→0o”FWJ=〃:o DMSGoto 1Lbl 2”END〃说明：X0Y0：起始点的坐标，FWJ：起始点的方位角，QDHAO：起点里程，ZDHAO 终点里程，R：半径，DQHAO：待求点里程，PY：偏移量4、竖曲线”QZHAO =〃J: ”H+B〃=B:”I1=〃C:”I2=〃D:”R=〃R:”T=〃T:((D-C)÷100) ÷Abs((D-C) ÷100)→A:Abs(D-C) ÷100)→S”T〃:SJ-T→XJ＋T→YLbl 0”DQ=〃LIf L<X:Then Goto 1:Else If L＞Y:Then Goto 1:If End:If End”GC=〃:B+C(L-J) ÷100+A(L-J+T) ÷(2R) →HLbl 1”END〃说明：QZHAO：曲中点里程，H+B：曲中点高程，I1：第一坡度，I2：第二坡度，R：半径，T：切线长，DQ：待求点里程。

卡西欧5800P坐标正算、反算程序
坐标正算ZB ZS
已知一个起算的坐标点，根据测得的方位角和平距推算另一点的坐标值。

程序如下:
“X0=”?B: 仪器点X坐标
“Y0=”?C: 仪器点Y坐标
Lbl 0:
“FWJ=”?F: 观测方位角
“PJ=”?D: 观测平距
Rec(D,F):
B+I→X:
C+J→Y:
Cls:
Fix 3:
“Xp=”: 待求点X坐标
“Yp=”: 待求点Y坐标
Locate 7,1,X: 显示在第一行
Locate 7,2,Y◢显示在第二行
Goto 0
注解：
（B，C）是仪器点的坐标，F为方位角，D为平距。

在本程序中，采用了循环（多次）运算多点与仪器点的关系计算法。

如果不必要，只算一个，那么把Lbl 0和Goto 0删除就行了。

坐标反算ZBFS
根据两已知坐标点求算其方位角及距离。

程序如下:
“X0=”?B: 仪器点X坐标
“Y0=”?C: 仪器点Y坐标
“Xp=”?X: 待求点X坐标
“Yp=”?Y: 待求点Y坐标
ol(X-B,Y-C):
J<0 => J+360→J:
I→D:
J→F:
Cls:
Fix 3:
“FWJ=”: 计算的方位角
F►DMS◢按度分秒显示
“PJ=”: 计算的平距
Locate 10,3,D◢平距显示在第三行
注解:
执行程序计算时,首先询问四变量两个点的坐标（B，C）、（X，Y）。

为了明了显示结果具体是何意思，特地使用字符串“PJ=”表示平距，字符串“FWJ=”表示方位角，但要注意是字符串而非变量名。

CASIO-5800积分法线路坐标正反算—雨中听风施工测量工作中较高档的计算器已很普及，施工放样的数据处理也应该采用程序化计算，以避免出错和提高测量工作的效率。

下面介绍本人炎汝(炎陵至汝城湘粤界)高速公路施工测量工作中应用的CASIO-5800计算器线元积分法坐标正反算程序（程序应用了测量空间里李艳阳老师提供的通用积分公式，在此表示感谢）。

该程序适用于任意线元的坐标、里程及偏距计算，只需更新数据库或者手工输入线元要素便可计之所需。

1.积分坐标正反算主程序(JF.ZBZFS)Lb1 0: Cls↓〝1.SZ => XY〞↓〝2.XY =>SZ〞:?Q↓〝灵活LH.K=0〞↓〝左线ZX.K=1〞↓〝右线YX.K=2〞:?K↓IF K=1:Then Prog〝YS1〞：Else IF K=2:Then Prog〝YS2〞: Else IF K=0:Then 〝QS.FW〞?A: 〝QS.LC〞?O:〝QS.X〞?U:〝QS.Y〞?V: 〝R0〞?P: 〝RN〞?R: 〝XY.CD〞?L↓Ifend：Ifend：Ifend:？S：S -（字母）O →X：IF X <（数字）0 Or X > L:Then 〝=CUOWU="◢Goto 0 ：Ifend :Q=2=>Goto 2↓Lbl 1：？Z：Prog 〝ZS〞↓Cls：Fix 3：〝N=〞：Locate 3,1,N:〝E=〞：Locate 3,2,E:〝FW=〞：F▲DMS◢Norm 2 :Goto 0↓Lbl 2：〝CX〞？B：〝CY〞？C：B→N：C→E： Prog〝FS〞↓Cls: Fix 3：〝S=〞：Locate 4,1,S:〝Z=〞：Locate 4,2,Z◢Norm 2 :Goto 02. 正算子程序(ZS)0.5(1÷R-1÷P)÷L→D:S-(字母)O→X↓U+∫（cos(A+(X÷P+DX²)×180÷（圆周率）∏)，（数字）0，X）→N↓V+∫（sim(A+(X÷P+DX²)×180÷（圆周率）∏)，（数字）0，X）→E↓A+(X÷P+DX²)×180÷（圆周率）∏→FN+Zcos(F+90)→N:E+Zsim(F+90)→E↓3. 反算子程序(FS)Lb1 1：(数字)0→Z: 1→Q: Prog〝ZS〞↓Pol(N-B+10^(-46) ,E-C+10^(-46)) ↓Isim(F-90-J)→W:S+W→S↓Abs(W)>0.0001=> Goto 1↓Lb1 2: (数字)0→Z : Prog〝ZS〞↓（C-E）÷sin(F+90)→Z数据库子程序(炎汝11标左线)（YS1）Goto 1 ↓Lb1 1 ↓IF S≥45871.346 And S<47089.759:Then194°17°39.6°→A:45871.346→O: 3966.911→U: 733.518→V:10＾45→P: 10＾45→R: 1218.413→L：Return:IfEnd↓IF S≥47089.759 And S<47254.759:Then194°17°39.6°→A:47089.759→O: 2786.220→U: 432.688→V:10＾45→P: -1480→R: 165→L： Return:IfEnd↓IF S≥47254.759 And S<47905.918:Then191°06°01.73°→A:47254.759→O: 2625.621→U: 394.932→V:-1480→P: -1480→R: 651.159→L：Return:IfEnd↓IF S≥47905.918 And S<48070.918:Then165°53°30.93°→A:47905.918→O: 1979.924→U:411.883→V:-1480→P:10＾45→R: 165→L： Return:IfEnd↓IF S≥48070.918 And S<48332.968:Then162°41°53°→A:48070.918→O: 1821.528→U: 458.013→V: 10＾45→P: 10＾45→R: 262.050→L:Return:IfEnd↓IF S≥48332.968 And S<48477.968:Then162°41°53°→A:48332.968→O: 1571.335→U: 535.949→V: 10＾45→P: 800→R: 145→L: Return:IfEnd↓IF S≥48477.968 And S<48915.797:Then167°53°25.75°→A: 48477.968→O: 1431.708→U: 574.858→V: 800→P: 800→R: 437.829→L： Return:IfEnd↓IF S≥48915.797 And S<49060.797:Then199°14°51.64°→A: 48915.797→O: 1000.161→U:547.941→V: 800→P: 10＾45→R: 145→L：Return:IfEnd↓IF S≥49060.797 And S<49556.014:Then204°26°24.40°→A: 49060.797→O: 866.451→U:491.983→V: 10＾45→P: 10＾45→R: 495.217→L：Return:IfEnd↓IF S≥49556.014 And S≤50236.938:Then204°26°24.40°→A: 49556.014→O: 415.608→U:287.091→V: -2500→P: -2500→R:680.924→L：Return:IfEnd↓Lb1 1: A：O：U：V：P：R:L数据库子程序(炎汝11标右线)（YS2）Goto 1 ↓Lb1 1 ↓IF S≥45700.089 And S<47055.633:Then194°17°39.6°→A:45700.089→O: 4117.556→U:500744.038→V:10＾45→P: 10＾45→R:1355.544→L:Return:IfEnd↓IF S≥47055.633 And S<47220.633:Then194°17°39.6°→A:47055.633→O: 2803.979→U: 409.351→V:10＾45→P: -1480→R: 165→L：Return:IfEnd↓IF S≥47220.633 And S<47871.792:Then191°06°01.73°→A:47220.633→O: 2643.381→U: 371.595→V:-1480→P:-1480→R: 651.159→L：Return:IfEnd↓IF S≥47871.792 And S<48036.792:Then165°53°30.93°→A: 47871.792→O: 1997.684→U:388.546→V:-1480→P: 10＾45→R: 165→L：Return:IfEnd↓IF S≥48036.792 And S<48291.984:Then162°41°53°→A: 48036.792→O: 1839.287→U: 434.676→V: 10＾45→P: 10＾45→R: 255.192→L：Return:IfEnd↓IF S≥48291.984 And S<48436.984:Then162°41°53°→A: 48291.984→O: 1595.643→U: 510.571→V: 10＾45→P: 800→R: 145→L：Return:IfEnd↓IF S≥48436.984 And S<48957.197:Then167°53°25.75°→A: 48436.984→O: 1456.015→U:549.480→V: 800→P: 800→R:520.213→L：Return:IfEnd↓IF S≥48957.197 And S<49102.197:Then205°08°52.79°→A: 48957.197→O: 948.224→U:491.452→V: 800→P: 10＾45→R: 145→L：Return:IfEnd↓IF S≥49102.197 And S<49510.315:Then210°20°25.50°→A: 49102.197→O: 820.975→U:422.046→V: 10＾45→P: 10＾45→R: 408.118→L：Return:IfEnd↓IF S≥49510.315 And S<49650.315:Then210°20°25.50°→A: 49510.315→O: 468.753→U:215.891→V: 10＾45→P: -1000→R:140→L：Return:IfEnd↓IF S≥49650.315 And S<49909.924:Then206°19°46.96°→A: 49650.315→O: 346.337→U:148.025→V: -1000→P: -1000→R:259.609→L：Return:IfEnd↓IF S≥49909.924 And S≤50049.924:Then191°27°18.76°→A: 49909.924→O: 101.402→U:64.201→V: -1000→P: 10＾45→R:140→L：Return:IfEnd↓Lb1 1: A：O：U：V：P：R：L数据库变量说明：S≥（线元的起点里程）［可要亦可不要］。

CASIO-5800计算器公路测量常用程序一、程序功能本程序由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,Ｚ◢＂Ｄ＝＂：Locate 6,4,Ｄ◢Ｇoto 0２．主程序２：高程序横坡程序(设计任意点高程及横坡)程序名：2GCLbI 0:Norm 2“KM”?Z:?D:Prog”H”:Fix 3:”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:Goto H: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 IfG>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,Ｚ◢＂Ｄ＝＂：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 3 LbI 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)÷si n(Ｚ［２］) →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:ThenQ+LI→H:Goto 0:Else If Z<F:Then Q+LI→H:Goto 0:Else If Z≤E:ThenQ+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:ElseAbs(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:Then25900→O:11587.421→U:1847.983→V:101。

卡西欧5800程序完整版＂ZS=1,FS=2＂？N: 1÷P →C:(P-R)÷(2HPH)→D:180÷π→E: If N=1: Then Goto 1: Else Goto 2: IfE nd←┚LbI 1:1 →I:0 →Z：Abs(K-Ｏ) →W: Prog＂ZS＂←┚F-90→F←┚＂JJ＂?K:＂PJ＂?L: ＂ZX=＂：X+Lcos(F+K) →Z[3]◢＂ZY=＂: Y+Lsin(F+K) →Z[4]◢＂QXFWJ=＂:F DMS◢Abs(Z) →Z:Goto 8←┚LbI 2:＂X＂?I:＂Y＂?J: Prog＂FS＂:＂FSZH=＂:Ｏ+W→K◢＂FSJL=＂:Z◢Lbl 8正算子程序（ZS）（不运行）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+ZcosF→X：Y+ZsinF→Y3反算子程序（FS）（不运行）G-90→T：Abs（（Y-V））cos（T）-（X-U）sin（T）→W：0→Z←┚Lbl 0←┚prog＂ZS＂：T+QEW（C+WD）→L：（J-Y）cosL-（I-X）sinL→Z：If Abs(Z)＜0.001: Then Goto1：E lse W+Z→W：IfEnd :Goto 0←┚Lbl 1：0→Z：Prog＂ZS＂：（J-Y）÷sinF→ZZX-WN（运行主程序，该程序适用于匝道，对全段完整非对称曲线计算太过麻烦）＂CZX＂?S(仪站X坐标):＂CZY＂?W（仪站Y坐标）LbI 0: ＂K＂?K: If K≤126.9（第1条单一曲线的终点下面同上）: Then Goto A：Else If K≤176.9 And K＞126.9: Then Goto B:Else If K≤259.562 And K＞176.9: Then Goto C:Else If K≤309.562 And K ＞259.562: Then Goto D：IfEnd:IfEnd:IfEnd:IfEnd←┚LbI A:92094.336（起点X坐标）→U:8287.962（起点Y坐标）→V:0起点桩号→O:145°11＇54＂（方位角）→G:126.9（曲线长度）→H:10^(45) （起点半径）→P:10^(45) （止点半径）→R:+1（路线左右偏，左负右正，直线输0）→Q:Prog ＂WN＂:Goto WLbI B:91990.134→U:8360.388→V:126.9→O:145°11＇54＂→G:50→H:10^(45) →P:150 →R:+1→Q:Pr og ＂WN＂:Goto WLbI C:91947.609→U:8386.570→V:176.9→O:154°44＇51＂→G:82.662→H:150 →P:150 →R:+1→Q:Pr og ＂WN＂:Goto WLbI D:91867.101→U:8399.992→V:259.562→O:186°19＇19＂→G:50→H:150 →P: 10^(45) →R:+1→Q: Prog ＂WN＂:Goto WLbI W:＂PJ1＂?T←┚（该项是只左右边桩再左右偏距，输0为不再左右偏，可以直接出方位角、距离，适用于桥梁桥台、盖梁、涵洞八字墙等放样，以路线前进方向，前-，后+）If T=0：Then Z[3]→Z: Z[4]→M: Goto P: Else Goto V: IfEnd←┚LbI V:＂JJ1＂? O←┚（左右边桩再偏夹角，用180°减去路线夹角）＂X1=＂: Z[3]+Tcos（F+N+ O）→U◢＂Y1=＂:Z[4]+Tsin（F+N+ O）→V◢U→Z: V→M: GotoP←┚LbI P:Pol(Z-S,M-W): If J≤0：Then J+360°→J: IfEnd←┚＂JL=＂:I◢（仪站至测点的距离）＂FWJ=＂:J DMS◢（仪站至测点的方位角）Goto0备注：1．规定以道路中线前进方向（即里程增大的方向）区分左右；当线元往左偏时Q=-1；当元线往右偏时，Q=1：当元线为直线时，Q=0。