公路测量卡西欧5800万能程序

一．使用说明：1.规定：（1）以道路中线的前进方向区分线路的左右偏转方向，线路左偏Q=-1，线路右偏Q=1，直线Q=0。

（2）当所求点位于线路中线时，Z=0；当位于左侧时，Z取负值；当位于右侧时，Z取正值。

（3）当线元为直线时，其起终点的曲率半径为无穷大，取10的45次方代替。

（4）当线元为圆曲线时，起终点曲率半径为圆曲线半径。

（5）当线元为缓和曲线时，起终点曲率半径与直线相接取10的45次方代替，与圆曲线相接取圆曲线半径。

2.输入与显示说明：（1）输入部分：N？：选择计算方式，输入1表示由里程边距计算坐标，输入2表示由坐标计算里程边距。

X0？线元起点X坐标。

Y0？线元起点Y坐标。

S0？线元起点里程。

F0？线元起点切线方位角。

LS？线元长度。

R0？线元起点曲率半径。

RN?线元终点曲率半径。

Q？线元左右偏标志（左偏=-1，右偏=1，直线=0）S?正算时所求点里程。

Z？正算时所求点距中线边距。

位于线路左侧输入负值，在右侧输入正值。

AR?斜交角度。

（为线路右角，在0-180度之间）X?反算时所求点X坐标。

Y？反算时所求点Y坐标。

（2）显示部分：XS=计算结果，所求点X坐标。

YS=计算结果，所求点Y坐标。

FS=计算结果，所求点切线方位角。

S=计算结果，反算时所求点里程Z=计算结果，反算时所求点边距。

二．程序清单1（用于任意曲线元的坐标计算）1.主程序(TYQXJS)″1.SZ=>XY″:″2.XY=>SZ″:?N:″X0″?U:″Y0″?V:″S0″?O: ″F0″?G:″LS″?H:″R0″?P:″RN″?R:?Q:1÷P→C:(P-R) ÷(2HPR) →D:180÷π→E:If N=1: THEN Goto 1:Else Goto 2:IfEnd:Lbl 1:?S:”D(Z-，Y+)=”?Z:90→T:”O=”?T:Abs(S-O) →W:Prog ″SUB1″:″XS″:X◢″YS″:Y◢″FS″:F◢Goto 1:Lbl 2:?X:?Y:I→X:J→Y:Prog ″SUB2″:″S″:O+W→S:◢″Z″:Z◢Goto 22.正算子程序(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))) →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) →F:X+Zcos(F+T) →X:Y+Zsin(F+T) →Y3.反算子程序(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 Abs(Z)<1E-6 :Then Goto 1:Else W+Z→W:Goto 0 : IfEnd:Lbl 1:0→Z:Prog ″SUB1″:(J-Y)÷sin(F)→Z 三．用于提前将曲线元要素输入子程序，计算时只输入里程桩号与边距。

卡西欧5800坐标正反算万能程序CASIO5800计算器程序曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序一、程序功能本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。

本程序也可以CASIO fx-5800P计算器计算器上运行。

二、源程序1.主程序(TYQXJS)"1.SZ => XY"："2.XY => SZ"：?N："X0"? U："Y0"? V："S0"? O："F0"? G："LS"? H："R0"? P："RN"? R：?Q←┘1÷P-> C：(P-R)÷(2HPR) -> D：180÷π->EIf N=1：Then Goto1：Else Goto 2：IfEndLbl 1：？S：？Z：Abs(S-O) -> W：Prog "SUB1"："XS"：X◢"YS"：Y◢Goto 1←┘Lbl 2：？X：？Y：X-> I：Y->J：Prog "SUB2"："S"：O+W-> S◢"Z"：Z◢Goto 22. 正算子程序(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))) -> 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. 反算子程序(SUB2)G-90-> T：Abs((Y-V)cosT-(X-U)sinT) -> W ：0-> Z：Lbl 0：Prog "SUB1"：T+QEW(C+WD) ->L：(J-Y)cosL-(I-X)sinL-> Z：If Abs(Z) <0.0000001：Then Goto1：Else W+Z ->W：Goto 0：IfEndLbl 1：0 ->Z：Prog "SUB1"：(J-Y)÷sinF-> Z三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

卡西欧FX5800---辛普森公式(万能公式)复化辛普森公式1.Lbl 0：“XA=”?A：“YA=”?B: “CA=”?C：“1/RA=”?D：“1/RB=”?E：“DKA=”?F:“DKB=”?G2.Lbl 1:“DKI=”?H：“DL=”?O：“DR=”?R：IF H＞G:THEN Goto0 IFEND3.(E-D)/Abs(G-F)→P:Abs(H-F) →Q: P×Q→I:D+I→T4.C+(I+2D)×Q×90/π→J5. C+(I/4+2D)Q×45/(2π) →M: C+(3I/4+2D)Q×135/(2π) →N6. C+(I/2+2D)Q×45/π→K7. A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12 →X8.B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12 →Y9.”FW=”:J▲DMS ▲10. “X=”:X▲11. “Y=”:Y▲12.“XL=”: X+Ocos(J-90) →U ▲13.“YL=”: Y+Osin(J-90) →V▲14.“XR=”:X+Rcos(J+90)→ W▲15.“YR=”:Y+Rsin(J+90)→ Z▲16.Goto 1程序结束程序说明：A- 曲线元起点A的坐标；B- 曲线元起点B的坐标；C- 曲线元起点A的切线坐标方位角；F- 曲线元起点A的里程；G- 曲线元起点B的里程；H- 曲线上待求点i的里程；D- 曲线元起点A的曲率；E- 曲线元终点B的曲率；XL-左边线点位X坐标；YL-左边线点位Y坐标；XR-右边线点位X坐标；YR-右边线点位Y坐标；X- 中线点位纵坐标；Y- 中线点位横坐标；DL-左边线距中线平距；DR-右边线距中线平距；该程序需要输入的数据为：（1）.曲线元起点A的坐标及切线坐标方位角，计算器上用“XA”，“YA”，“CA”显示；（2）.曲线元起点A和B的曲率，计算器上用I÷RA，I÷RB显示（曲线左偏时取“-”）；（3）.曲线元起点A和终点B的里程，计算器上用“DKA”，“DKB”显示；（4）.输入待求点里程和该点距左右的水平距离，计算器上用“DKI”，“DL”，“DR”显示；每算完一个待求点的中线及边线坐标，程序又让输入下一点的“DKI”，“DL”，“DR”当输入的“DKI”大于“DKB”时，此时输入下一个曲线元起点的曲率和里程，即可计算下一个曲线中线及边线点位坐标。

1、主程序：程序名：ZCX输入：”K”?K←If K >540980.584 (直线起点里程) And K<= （小于等于符号）541358.301:（直缓点里程）Then 2890632.365→A: 526217.878→B:79°34′47.8″→F:540980.584→J:Prog”ZX”(A,B,F,J分别为起算点坐标，方位角，里程) Prog”PJ”(运行求偏距的程序)IfEndIf K>541358.301 （直缓点里程）And K<=541488.301：（缓圆点里程）Then 2890700.680→A:526589.366→B: 79°34′47.8″→F: 541358.301→J:1→G:3500→R:130→L:Prog”HHQX”（A,B,F,J,G,L分别为直缓点坐标，方位角，里程，偏向，圆曲线半径，缓和曲线长度）Prog”PJ”IfEndIf K>541488.301（缓圆点里程）And K<=541694.275:（圆缓点里程）Then 2890723.4→A: 526717.363→B:80°38′38.4″→F:541488.301→J:1→G:3500→R:Prog”YQX”（A,B,F,J,G,L分别为缓圆点坐标，方位角，里程，偏向，圆曲线半径，缓和曲线长度）Prog”PJ”IfEndIf K>541694.275 （圆缓点里程）And K<=541824.275：（缓直点里程）Then 2890762.83→A:527050.912→B: 85°04′47.7″→F: 541824.275→J:-1→G:3500→R:130→L:Prog”HHQX”（A,B,F,J,G,L分别为缓直点坐标，方位角，里程，偏向，圆曲线半径，缓和曲线长度）Prog”PJ”IfEndIf K>541824.275 （缓直点里程）And K<=542167.602：（下一段直缓点里程）Then 2890762.83→A:527050.912→B: 85°04′47.7″→F: 541824.275→J: Prog”ZX”(A,B,F,J分别为缓直点坐标，方位角，里程)注：以上程序为一段标准曲线（包括直线，缓和曲线，圆曲线，缓和曲线，直线）的主程序，程序名为ZX，HHQX，YQX2、直线程序程序名：ZXA+(K – J )xcos(F) →X:”X0=”:X◢B+(K – J )xsin(F) →Y:”Y0=”:Y◢F→W:”F=”:W◢3、圆曲线程序程序名：YQXA+2xRxSin（90x(K- J)÷∏÷R）XCos(F+90xGx(K- J) ÷∏÷R) →X:”X0=”:X◢B+2xRxSin（90x(K- J)÷∏÷R）XSin (F+90xGx(K- J) ÷∏÷R) →Y:”Y0=”:Y◢F+180xGx(K- J) ÷∏÷R→W: F=”:W◢4、缓和曲线程序程序名：HHQXA+((K- J)-(K- J)5÷(90xR²xL²))xcos(F+30xGx(K- J) ²÷∏÷R÷L) →X: “X0=”:X◢B+((K- J)-(K- J)5÷(90xR²xL²))xsin(F+30xGx(K- J) ²÷∏÷R÷L) →Y: “Y0=”:Y◢F+90xGx(K- J) ²÷∏÷R÷L→W:”F=”:W◢5、偏距程序50→DimZ“Z=”?Z:(Z为偏离线路中线的距离)X+Zxcos(W-90) →U:”XL=”:U◢(左边X坐标)Y+Zxsin(W-90) →V:”YL=”:V◢(左边Y坐标)X+Zxcos(W+90) →D:”XR=”:D◢(右边X坐标)Y+Zxsin(W+90) →E:”YR=”:E◢(右边Y坐标)Return。

卡西欧5800程序(修改于湘桂)R（主程序，通过里程求中桩、边桩坐标，在输入测站坐标后并求极坐标放样要素）1→M :“NEWZ CIRCLE”？M : If M=-1: Then “JD”？U :“1-TR-1-TL” ? N : “AT ” ? O : “JX” ? S : “JY ” ? P : “L0 ” ? H : “R” ? R : “ZJ” ? A : IfEnd↙“OPP NAME ”？W : If W=-1 : Then “OPPX”？B : “OPPY”？C : Else Prog“CONTREL”:IfEnd↙Lbl 0: “C”？K : If Int(K÷1000)>0 : Then Int(K÷1000)→V :Else “V”？V ：K+1000V→K：IfEnd↙If M=1:Then Prog“CIRCLE” IfEnd↙H÷2-H^(3)÷(240R2)+(R+H2÷(24R))tan（A÷2）→T:πR(A-180H÷π÷R)÷180+2H→L↙U-T→G:K→Q↙If Q≤G+L÷2:Then Q-G→Q: Else G+L-Q→Q: IfEnd↙If Q＜0:Then Prog “L”:Goto1: IfEnd↙If Q≤H:Then Q-Q^(5)÷(40R2H2)→X:Q^(3)÷(6RH)-Q^(7)÷(336R^(3)H^(3)) →Y:90Q2÷(πRH) →F: Else 180(Q-H÷2)÷(πR) →F:Rsin(F)+H÷2-H^(3)÷(240R2) →X:R(1-Cos(F))+H2÷(24R) →Y: IfEnd↙Pol(T-X,Y) ↙If K≤G+L÷2:Then O-180-NJ→E:NF+O→F: “F=”:F▲DMS◢ ElseNA+NJ+ O→E: O+NA-NF→F: “F=”:F▲DMS◢IfEnd↙“X=”:S+Icos(E)→X◢“Y=”:P+Isin(E)→Y◢If F＜0:Then 360+F→F: IfEnd↙If F≥360:Then F-360→F: IfEnd↙Lbl1: “÷1”？Z: If Z=1:Then Prog “I”: IfEnd↙If Z=2:Then Prog “LD”: IfEnd↙If Z=-1:Then Prog “Q”: “X=”:X◢“Y=”: Y◢IfEnd↙Prog “P”↙Goto 0↙P（子程序，受主程序R调用，求极坐标放样要素:方位角及平距）Pol(X-B,Y-C): “D=”:I◢If J＜0:Then 360+J→J: Else 0+J→J: IfEnd↙“AQ=”:J▲DMS◢L（子程序，受主程序R调用，求直线上点位坐标及切线方位角）T-Q→Q:If K＜G:Then O+180→E: O→F: “X=”:S+Qcos(E)→X◢“Y=”:P+Qsin(E) →Y◢ Else O+NA→F: “X=” S+Qcos(F) →X◢“Y=”:P+Qsin(F) →Y◢IfEnd↙I（子程序，受主程序R调用，求边桩坐标）“RA”？E：“D”？D:E+F→E: If E≥360:Then E-360→E: IfEnd↙X+Dcos(E)→X: Y+Dsin(E)→Y: “X=”: X◢“Y=”: Y◢E→F↙Q（子程序，受主程序R调用，求副线边桩坐标）“1A”？E：“1D”？D：E+F→E↙X+Dcos(E)→X: Y+Dsin(E)→Y:E→F↙“2R”？E: “2D”？D: E+F→E↙X+Dcos(E)→X: Y+Dsin(E)→Y↙E→F↙LD（子程序，受主程序R调用，通过坐标反求里程和偏距）40→DimZ↙O→Z[1]:U→Z[2]:S→Z[3]: P→Z[4] ↙H÷2-H^(3)÷(240R2)→M: H2÷(24R)-H^(4)÷(2688R^(3))→V:90H2÷(πR)→D:(R+V)Tan(A÷2)+M→T: πRA÷180+H→K:O+NA→J:Prog“JD”:J→Z[21]:U-T+K→Z[22]: S-Tcos(O)→Z[5]:P-Tsin(O)→Z[6]:Z[5]+Mcos(O)+(R+V)cos(O+90N)→Z[7]:Z[6]+Msin(O)+ (R+V)sin(O+90N)→Z[8] ↙H-H^(5)÷(40R2H2)→Z:H^(3)÷(6RH)-H^(7)÷(336R^(3)H^(3))→C: Pol(C,Z):90-J→J:Prog“JD”: O+JN→J:Prog“JD”↙Z[5]+Icos(J)→Z[9]:Z[6]+Isin(J)→Z[10]:S+Tcos(Z[21])→Z[11]: P+Tsin(Z[21])→Z[12] ↙Pol(Z[10]-Z[8],Z[9]-Z[7]):90-J→J:Prog“JD”:J→Z[13] ↙J+(A-2D)N→J: Prog“JD”: J→Z[14] ↙“CEDIAN”◢Lbl 0: ↙“B”？B:If B=0:Then Prog “R”: IfEnd↙“XN=”？X: “YN=”？Y:Pol(Y-Z[8],X-Z[7]) ↙90-J→J: Prog“JD”:J→W:J-Z[13]→J↙Prog“JD”:J→Z[15] ↙W-Z[14]→J↙Prog“JD”:J→Z[16] ↙If N=1:Then Goto 1: IfEnd↙If N=-1:Then Goto 2: IfEnd↙Lbl 1: ↙If Z[15]＞180：Then Goto 3: IfEnd↙If Z[16]≤180：Then Goto 5:Else Goto 4 : IfEnd↙Lbl 2: ↙If Z[15]＜180：Then Goto 3: IfEnd↙If Z[16] ＞180：Then Goto 5:Else Goto 4 : IfEnd↙Lbl 3: ↙“H1”◢Pol(Y-Z[6],X-Z[5]):90-J→J: Prog“JD”↙J-Z[1]+90→J: Prog“JD”↙Isin(J)→G↙If G＞0:Then 0→L:Goto 6: IfEnd↙If G≤0:Then “DK=”:Z[2]-T+G→Z[19]◢If N=1:Then “LD=”:-NIcos(J)→Z[20]◢Else “LD=”:NIcos(J) →Z[20]◢IfEnd↙“SD=1,LD≠1”？D:If D=1:Then Prog “SDCQJC”: IfEnd↙Goto 0↙Lbl 4: ↙“Y”◢If N=-1:Then 360-Z[15]→Z[15]: IfEnd: “DK=”:Z[2]-T+H+πR(Z[15])÷180→Z[19]◢If N=1:Then “LD=”:-(I-R)→Z[20]◢Else “LD=” :I-R→Z[20]◢IfEnd↙“SD=1,LD≠1”？D:If D=1:Then Prog “SDCQJC”: IfEnd↙Goto 0↙Lbl 5: ↙“H2”◢Pol(Y-Z[12],X-Z[11]):90-J→J: Prog“JD”:J→W:Z[21]-90→J: Prog“JD”:W-J→J: Prog“JD”:Isin(J)→G:If G＜0:Then 0→L: -G→G↙Goto 7: IfEnd↙If G≥0:Then “DK=”:Z[22]+G→Z[19]◢IfEnd: If N=1:Then “LD=”:-NIcos(J)→Z[20]◢Else “LD=”:NIcos(J)→Z[20]◢IfEnd↙“SD=1,LD≠1”？D:If D=1:Then Prog “SDCQJC”: IfEnd↙Goto 0↙Lbl 6: ↙L+G→L:L-L^(5)÷(40R2H2)→Z:L^(3)÷(6RH)-L^(7)÷(336R^(3)H^(3)) →C:Pol(C,Z):90-J→J: Prog“JD”:J→W:Z[1]+WN→J: Prog“JD”:Z[5]+Icos(J)→Z[17]: Z[6]+Isin(J) →Z[18] ↙J+2WN→J: Prog“JD”:J-90→J: Prog“JD”:J→W↙Pol（Y-Z[18],X-Z[17]）:90-J→J: Prog“JD”↙J-W→J: Prog“JD”:Isin(J)→G:If Abs(G)＞0.001:Then Goto 6: IfEnd↙“DK=”:Z[2]-T+L+G→Z[19]◢If N=1:Then “LD=”:-Ncos(J) →Z[20]◢Else“LD=”:NIcos(J) →Z[20]◢IfEnd↙“SD=1,LD≠1”？D:If D=1:Then Prog “SDCQJC”: IfEnd↙Goto 0↙Lbl 7: ↙L+G→L:L-L^(5)÷(40R2H2)→Z:L^(3)÷(6RH)-L^(7)÷(336R^(3)H^(3)) →C:Pol(C,Z):90-J→J: Prog“JD”:J→W:Z[21]-180-WN→J: Prog“JD”↙Z[11]+Icos(J)→Z[17]:Z[12]+Isin(J)→Z[18]:J-2WN→J: Prog“JD”↙J-90→J: Prog“JD”:J→W↙Pol(Y-Z[18],X-Z[17]):90-J→J: Prog“JD”↙J-W→J:Prog “JD”↙Isin(J)→G: If Abs(G)＞0.001:Then Goto 7:IfEnd↙“DK=”:Z[22]-L-G→Z[19]◢If N=1:Then “LD=”:NIcos(J)→Z[20]◢Else“LD=”:-NIcos(J) →Z[20]◢IfEnd↙“SD=1,LD≠1”？D:If D=1:Then Prog “SDCQJC”: IfEnd↙Goto 0↙JD（子程序，受LD程序调用，计算角度归化值）If J＜0:Then J+360→J: IfEnd↙If J≥360:Then J-360→J: IfEnd↙CIRCLE（子程序，线路平面曲线要素，受主程序R调用）If K≤缓直点里程:Then 交点里程→U:交点X坐标→S:交点Y坐标→P:转角值→A:-1（曲线段左偏）或1（曲线段右偏）→N:圆曲线半径→R:缓和曲线长→H:直线段起始方位角→O:Goto 1: IfEnd↙If K≤……（同上输完全部曲线要素为止）Lbl 1↙CONTREL（控制点坐标数据库子程序，受主程序R和ZBFS调用）If W=控制点点号:Then 控制点X坐标→B:控制点Y坐标→C: IfEnd ↙If W=……（同上输完全部控制点数据为止）ZBFS（主程序，求两已知点方位角及平距,或已知两点求后方交会点坐标）Lbl 0:1→W: “NEW OPP”？W:If W=1:Then “OPPX”？X: “OPPY”？Y:Else Prog “CONTREL”:B→X:C→Y: IfEnd↙Lbl 1:-1→W: “NEW CPP”？W: If W=-1:Then “CPPX”？S: “CPPY”？T: Else Prog “CONTREL”:B→S:C→T: IfEnd↙Pol(S-X,T-Y): “D=”:I◢If J＜0:Then J+360→J: Else 0+J→J: IfEnd↙“AQ=”:J▲DMS◢Lbl 2:0→W: “ZBFS=0,SQS≠0”？W:If W=0:Then Goto 1: ElseProg “SQS”: IfEnd↙Goto 0↙SQS（子程序，受主程序ZBFS调用，通过角度边长求交会点坐标值）Lbl 0: “S1=”？A: “J1=？M: “J2=”？N: “S2=？D:N-M→Q: If Q ＜0:Then Q+360→Q: IfEnd↙If Q＜45°0′0″:Then “ERROR PLEASE REPEAT”◢Goto 0: IfEnd↙√(A2+D2-2ADos(Q))→E: “S3=”:E◢If Abs(A-D)≥60:Then “ERROR PLEASE REPEAT”◢Goto 0: IfEnd↙If Abs(I-E)＜0.01:Then cos-1((A2+E2-D2)÷(2AE))→R: cos-1((D2+E2-A2) ÷(2DE))→G: Else “ERROR PLEASE REPEAT”◢Goto 0: IfEnd↙If J＜0:Then J+360→J:IfEnd↙(X+Acos(J+R)+S+Dcos(J+180-G))÷2→F: “X=”:F◢(Y+Asin(J+R)+T+Dsin(J+180-G))÷2→P: “Y=”:P◢SQX（主程序，竖曲线求线路各里程点中桩设计高程）1→S：“NEWZ SQXK”？S: If S=-1:Then ？U:？R:？T:？N:？I:？F: “HB”？B: IfEnd↙Lbl 0: “C”？K : If Int(K÷1000)>0 : Then Int(K÷1000)→V :Else “V”？V ：K+1000V→K：IfEnd↙If S≠-1:Then Prog “SQXK”: IfEnd↙U-T→G:K→Q: If Q≤G:Then Q-U→Q: Else Q-G→W: IfEnd↙If Q≤G:Then “H=”:B+QN÷1000→H◢Goto 1↙Else If Q≤U+T:Then “H=”:B+(Q-U)N÷1000+FW2÷(2R)→H◢Goto 1↙Else If Q≥U+T:Then “H=”:B+(Q-U)I÷1000→H◢Goto 1:IfEnd↙Lbl 1: “P”？P: “O=”:H-P→O◢Goto 0↙SQXK（子程序，受主程序SQX调用，竖曲线要素）If K≤竖曲线圆直点里程：Then 变坡点里程→U:变坡点设计高程→B:竖曲线圆曲线半径→R:竖曲线切线长→T:竖曲线直线段坡度常数（上、下坡有正负之分）→N:变坡点常数（凸为-1，凹为1）→F: Goto 0: IfEnd↙If K≤……（同上输完全部竖曲线要素为止）Lbl 0↙SDCQJC（子程序，受LD程序调用，采集隧道支护位置坐标高程数据一次性计算出该支护位置的超欠挖量）40→DimZ↙Lbl 0:Z[19]→Z[23]:Prog “SQXZ”↙Z[24]-Z[27]→Z[30]: Z[23]→Z[31] ↙If Z[31]≤Z[30]:Then Z[31]-Z[24]→Z[31]:Else Z[31]-Z[30] →Z[32]: IfEnd↙If Z[31]≤Z[30]:Then Z[25]+(Z[31])×(Z[28])÷1000→Z[33]: Goto 1↙Else If Z[31]≤Z[24]+Z[27]:Then Z[25]+(Z[31]-Z[24])×(Z[28])÷1000+(Z[29])×(Z[32])2÷(2×Z[26])→Z[33]: Goto 1↙Else If Z[31]≥Z[24]+Z[27]:Then Z[25]+(Z[31]-Z[24])×I÷1000→Z[33]: Goto 1: IfEnd↙Lbl 1: “H0=”:Z[33]+内轨面至圆心高差常数→Z[33]◢“HN”？K: “Δh=”:K-Z[33]→Z[34]◢If Z[20]≥隧道中心线与设计中心线的偏移量ω:Then Z[20]-ω→Z[35]:Else ω- Z[20]→Z[35] : IfEnd↙“RN=”:√（（Z[34]）2+（Z[35]）2）→Z[36]◢“R0”？Z: “RCQ=”: Z[36]-Z→Z[37]◢SQXZ（子程序，受主程序SDCQJC调用，竖曲线要素）40→DimZ↙If Z[23]≤竖曲线圆直点里程：Then 变坡点里程→Z[24]:变坡点设计高程→Z[25]:竖曲线圆曲线半径→Z[26]:竖曲线切线长→Z[27]:竖曲线直线段坡度常数（上、下坡有正负之分）→Z[28]:变坡点常数（凸为-1，凹为1）→Z[29]:IfEnd↙文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持. If Z[23]≤……（同上输完全部竖曲线要素为止）Lbl 0↙。

公路曲线大地坐标施工放样程序（CASIO fx-58000）一、主程序28:""?:""?:""?:""?:.5^(3)/(240):""?:"":(/(24R )-G ^(4)/(2688R ^(3)))ta n (0.5)T "":/180+G L""?:"":[3]"1,2"?:""?:""?:""?:10:"Z B J SD im Z X C Z Y C C L S G R R G G R Q A A T Q R G L R A Z H O H Z O Z Y Z H J D N P J DE MF F L b K π→====-→==++→=→==+→=====▲▲▲"?:","?:""?:0:0:0:P r "":P r "":0:K V N W V B B K O S IfS T h e n S XY U o g C o g D G o to IfE n d ==-+=-→↵≤→→→↵[][]IfS L :T h e n T +C O S (A )(T +S -L )X :S i n (A )(T +S -L )Y :A U :P ro g "C ":P ro g "D ":G o to 0:IfE n d IfS L -G :T h e n L -S S :P ro g "B ":X Z1:YZ2:T -Z [2]S in (A )+C o s (A )(T -Z [1])X :Z [2]C o s (A )+S in (A )(T -Z [1])Y :L -S S :A -U U :P ro g "C ":P ro g "D ":G o ≥→→→↵≥→→→→→→→2222to 0:IfE n d IfS G :T h e n 90(2S -G )/(R )U :Q +R S in (U )X :G/(24R )-G ^(4)/(2688R ^(3))+R (1-C o s (U ))Y :P ro g "C ":P ro g "D ":G o to :If E n d If S G :T h e n K -0S :S -S ^(5)/(40RG)X :S ^(3)/(6R G )-S ^(7)/(336R ^(3)G ^(3))Y :90S /(R G )ππ↵≥→→→↵≤→→→U :P ro g "C ":P ro g "D ":G o to 0:I f E n d →↵二、子程序 BSC"X=Dta注：”/”表示÷三、程序变量说明以程序运行时计算器显示为准：XC?EC?提示输入架仪点的坐标LS?提示输入测设曲线缓和曲线长R?提示输入测设曲线半径A?提示输入测设曲线转角（以度分秒输入）T= L=显示计算的曲线元素（检验LS、R、A是否输入正确）ZH? 提示输入测设曲线直缓（圆）桩桩号（对比设计资料可检查直缓（直圆）桩桩号输入是否正确Y1,Z2?判断曲线转向，右转输入1，左转输入2JDN?JDE?提示输入交点的N和E坐标F?提示输入曲线计算起始边的坐标方位角（以度分秒格式输入）K?提示输入欲放样点的桩号V?提示输入欲放样点相对中点的法线方向的距离，位于曲线内侧输入负值，外侧为正值，测设中桩时输0B?提示输入欲放样点相对中点的切线方向的距离，位于曲线前进方向输入正值，后退方向输入负值。

都是我自己编的，用起来还可以。

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(0.5E) →FReC(Abs(F),D+0.5E)：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(0.5R(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) →HGoto 0Lbl 1＂END″说明：QZHAO：曲中点里程，H+B：曲中点高程，I1：第一坡度，I2：第二坡度，R：半径，T：切线长，DQ：待求点里程。

卡西欧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:反算出的偏距。

卡西欧5800计算程序（请仔细核对）（程序输完毕后记得先在普通模式扩展变量，输入：35→Dimz，按EXE键，出现Done,扩展成功）ZHYHZ（不运行，该程序为主线坐标的计算程序）E2÷24÷R→Z[1]: E÷2-E^（3）÷(240R2) →Z[2]: F2÷24÷R→Z[3]: F÷2-F^3÷(240R2) →Z[4]: (R+Z[1])tan(H÷2)+Z[2]-（Z[1]-Z[3])÷sin （H）→Z[5]： (R+Z[3])tan (H÷2)+Z[4]+(Z[1]-Z[3])÷sin (H)→Z[6]: tan-1 ((R+Z[1])÷(Z[5]-Z[2]))→ Z[7]: tan-1 ((R+Z[3])÷(Z[6]-Z[4])) →Z[8]: E÷R×90°÷π→Z[9]: F÷R×90°÷π→Z[10]: G+DH→Z[11]: A-Z[5]→Z[12]: Z[12]+E→Z[13]: Z[13]+( 90°-Z[7]-Z[9])Rπ÷180°→Z[14]: Z[13]+(H-Z[9]-Z[10]) Rπ÷180°→Z[15]: Z[15]+F→Z[16]: B+Z[5]cos (G-180°)→Z[17]: C+Z[5]sin (G-180°)→Z[18]: B+Z[6]cos （Z[11]）→Z[19]: C+Z[6]sin （Z[11]）→Z[20] ←┚If K≤Z[14]:Then E→Q: K-Z[12] →L: Else F→Q: Z[16]-K→L:IfEnd ←┚If K≤Z[12]: Then Goto 1：Else If K＞Z[12] And K≤Z[13]: Then Goto 2: Else If K＞Z[13] And K≤Z[14]: Then Goto 3: Else If K ＞Z[14] And K≤Z[15]: Then Goto 4: Else If K＞Z[15] And K≤Z[16]: Then Goto 5: Else If K＞Z[16]: Then Goto 6:IfEnd:IfEnd:IfEnd:IfEnd:IfEnd:IfEnd←┚LbI 1:G→Z[26]: Z[17]-Lcos (G-180°)→Z[21]: Z[18]-Lsin(G-180°) →Z[22]: Goto 9←┚LbI 2: L-L^5÷(40R2Q2) →Z[23]: L^（3）÷(6RQ)- L^（7）÷(336R^（3）Q^（3）) →Z[24]: G+DL2×90°÷(QRπ) →Z[26]: Goto 7←┚LbI 3: (L-Q)×180°÷(Rπ)+Z[9]→Z[25]: Rsin （Z[25]）+Z[2]→Z[23]: R(1-cos （Z[25])）+Z[1]→Z[24]: G+DZ[25]→Z[26]: Goto 7←┚LbI 4: (L-Q)×180°÷(Rπ)+Z[10]→Z[25]: Rsin（ Z[25]）+Z[4]→Z[23]: R(1-cos （Z[25]）)+Z[3]→Z[24]: Z[11]-DZ[25]→Z[26]: Goto 8←┚LbI 5: L-L^5÷(40R2Q2)→Z[23]: L^（3）÷(6RQ)- L^7÷(336R^（3）Q^（3）)→Z[24]: Z[11]-DL2×90°÷(QRπ)→Z[26]: Goto 8←┚LbI 6: Z[11]→Z[26]: Z[19]-L×cosZ[11]→Z[21]: Z[20]-Lsin （Z[11]）→Z[22]: Goto 9←┚LbI 7: Z[17]+Z[23]cos（G）+Z[24]cos(G+D×90°)→Z[21]: Z[18]+ Z[23]sin（G）+Z[24]sin（G+D×90°)→Z[22]: Goto 9←┚LbI 8: Z[19]+Z[23]cos(Z[11]+180°)+Z[24]cos(Z[11]+ D×90°)→Z[21]: Z[20]+Z[23]sin(Z[11]+180°)+Z[24]sin（Z[11]+D×90°)→Z[22]: Goto 9←┚LbI 9: ＂JJ＂?N:＂PJ＂?P: ＂X=＂：Z[21]+Pcos(z[26]+N)→X◢＂Y=＂: Z[22]+Psin(Z[26]+N)→Y◢＂QXFWJ=＂DMS◢←┚ZXZB（该程序为主线程序，只运行该程序）＂CZX＂?S:＂CZY＂?W ←┚LbI 0: ＂K＂?K: If K＜54074.263（第1条曲线的终点桩号带大桩号，也就是是缓直点桩号，下面同上）: Then Goto A：Else If K＜56247.409 （第2条曲线的终点桩号带大桩号）And K≥54074.263（第1条曲线的终点桩号带大桩号）: Then Goto B:Else If K＜58007.266 And K ≥56247.409: Then Goto C:IfEnd:IfEnd:IfEnd←┚（有几条曲线就是几个IfEnd）LbI A: 53203.87266（第1条曲线的交点桩号）→ A: 3560815.660（第1条曲线交点的X坐标）→B:484158.7230（第1条曲线交点的Y 坐标）→ C:76°10＇54.56＂（计算方位角）→G: 16°49＇2.2＂（转角）→ H:-1（路线左右偏，左-，右+）→D:4500（半径）→R:430→（第1缓和曲线的长度）E:430（第2缓和曲线的长度）→F:Prog＂ZHYHZ＂: GotoW←┚（下面同上）LbI B: 55278.9644→ A: 3561878.191→B: 485952.830→ C: 59°21＇52.34＂→G: 19°22＇48.7＂→ H:+1→D:4500→R:430→E:430→F:Prog＂ZHYHZ＂: GotoW←┚LbI C: 57236.8579→ A: 3562263.305→B: 487888.000→ C: 78°44＇41＂→G:14°13＇19.1＂→ H:-1→D:4500→R:430→E:430→F:Prog＂ZHYHZ＂: GotoW←┚LbI W:＂PJ1＂?T←┚If T=0：Then X→Z: Y→M: GotoP: Else GotoV: IfEnd←┚LbI V:＂JJ1＂?O←┚＂X1=＂: X+Tcos（Z[26]+N+O）→U◢＂Y1=＂：Y+Tsin（Z[26]+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=＂DMS◢Goto0显示：1、CZX？（输入仪站X坐标）2、CZY？(输入仪站Y坐标)3、K？（输入所求桩号，带大桩号）4、JJ？（输入夹角，正交输90°，斜交输夹角）5、PJ？（输入偏距，左负右正，中桩输0）6、X=（所求X坐标）7、Y=（所求Y坐标）8、QXFWJ（所求曲线方位角）9、PJ1?（输入左右边桩再偏距离，前负后正，不需要再偏输0）10、JJ1？（输入左右边桩再偏夹角，上面输0，该步不会显示）11、X1= (所求左右边桩再偏X坐标)12、Y1= (所求左右边桩再偏Y坐标)13、JL= （仪站点至所求点的距离）14、FWJ=（仪站点至所求点的方位角）再往下继续计算Z-ZBFSLbI 0＂DQX＂?W:＂DQY＂?S←┚＂JDX＂?A:＂JDY＂?B:＂A0＂?G:＂ZJ＂?H:＂Z(-1)Y(1)＂?D:＂R＂?R:＂LS1＂?E:＂LS2＂?F:＂K＂?K←┚LbI 1←┚Norm 1: Prog＂ZHYHZ＂: Pol(W-Z[21],S-Z[22]): If J＞0：Then J →J:I→I: Else J+360°→J: I=I:IfEnd←┚J-Z[26]→Z[27]:Icos(Z[27])→Z[28]:K+Z[28]→K: If Abs（Z[28]≥0.001：Then Goto1：Else Fix 3：IfEnd：Isin（Z[27]）→Z[29]: Fix 3:＂FSJL=＂:Z[29] ◢＂FSK=＂:K◢Goto 0反算时要在反算区域正算一次，获取曲线要素显示1、DQX?(输入反算时的X坐标)2、DQY?(输入反算时的Y坐标)3、然后一直按＂EXE＂出现下面的字符4、FSJL= (所求反算左右偏距，左负右正)5、FSK= (所求反算的桩号)6、再按继续运算WN（不运行）＂ZS=1,FS=2＂？N: 1÷P→C:(P-R)÷(2HPH)→D:180÷π→E: If N=1: Then Goto 1: Else Goto 2: IfEnd←┚LbI 1:1→I:0→Z：Abs(K-Ｏ)→W: Prog＂ZS＂←┚F-90→F←┚＂JJ＂?K:＂PJ＂?L: ＂ZX=＂：Z[7]+Lcos(F+K)→Z[3]◢＂ZY=＂: Z[8]+Lsin(F+K)→Z[4]◢＂QXFWJ=＂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÷π→E0.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）））→Z[7]：V+W（Asin（G+QEKW（C+KWD））+Bsin（G+QELW（C+LWD））+Bsin（G+QEFW（C+FWD））+Asin（G+QEMW（C+MWD）））→Z[8]：G+QEW（C+WD）+90→F：Z[7]+ZcosF→Z[7]：Z[8]+ZsinF→Z[8]3反算子程序（FS）（不运行）G-90→T：Abs（（Z[8]-V）cosT-（Z[7]-U）sinT）→W：0→Z←┚Lbl 0←┚prog＂ZS＂：T+QEW（C+WD）→L：（J-Z[8]）cosL-（I-Z[7]）sinL→Z：If Abs(Z)＜0.001: Then Goto1：Else 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≤218.50（第1条单一曲线的终点下面同上）: Then Goto A：Else If K≤373.307 And K＞218.50: Then Goto B:Else If K≤526.691 And K＞373.307: Then Goto C：Else If K≤946.255 And K＞526.691: Then Goto D IfEnd:IfEnd:IfEnd: IfEnd←┚LbI A:3803701.508（起点X坐标）→U:494281.813（起点Y坐标）→V:0起点桩号→O:293°47＇49＂（方位角）→G:218.50（曲线长度）→H:3990.75 （起点半径）→P:3990.75 （止点半径）→R:+1（路线左右偏，左负右正）→Q:Prog ＂WN＂:Goto WLbI B:3803795.099→U:494084.403→V:218.50→O:296°56＇2.4＂→G:154.807→H:10^45 →P:10^45 →R:0→Q:Prog ＂WN＂:Goto W LbI C:3803865.221→U:493946.388→V:373.307→O: 296°56＇2.4＂→G:153.384→H:10^45 →P:496.5 →R:+1→Q:Prog ＂WN＂:Goto W LbI D:3803941.562→U:493813.538→V:526.691→O: 305°47＇3.15＂→G:419.564→H:496.5 →P:496.5 →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+K+O）→U◢＂Y1=＂ :Z[4]+Tsin（F+K+ 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=＂DMS◢（仪站至测点的方位角）Goto01．规定以道路中线前进方向（即里程增大的方向）区分左右；当线元往左偏时Q=-1；当元线往右偏时，Q=1：当元线为直线时，Q=0°2．当所求点位于中线时，Z=0；当位于中线左侧时，Z取负值：当位于中线右侧时，Z取正值°当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次方代替°3．当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径°4．当线元为完整的缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替°与圆曲线相接时，曲率半径等于圆曲线的半径，止点跟直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径°5．当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线衔接时，曲率半径等于圆曲线的半径°止点与直线相接时，曲率半径等于设计规定的值；与圆曲线衔接时，曲率半径等于圆曲线的半径°二、输入与显示说明1.输入部分：2.. 输入3.ZS=1，FS=2，选择计算方式，输入1正算；输入2反算4.数据库的输入5.U？线元起点的X坐标6.V？线元起点的Y坐标7.O？线元起点的里程8.G？线元的起点切点方位角9.H？线元的长度10.P？线元起点的曲率的半径11.R？线元止点的曲率的半径.12.ZY？线元左右偏标志（左偏Q=-1，右偏Q=1，直线段Q=0）13.计算时：14.K？正算时所求的里程15.JJ？夹角，（正交90°，斜交直接输角度）16.PJ？偏距（左负右正）17..结果18.ZX=正算时，计算得出的测点的X坐标19.ZY=正算时，计算得出的测点的Y坐标20.QXFWJ=正算时，测点对应中线点的切线方位角21.PJ1?：左右桩号再偏距离，不需要时输022.JJ1? :左右边桩再偏角度23.JL=仪站与测点的距离24.FWJ=仪站与测点的方位角25.反算时输入：26.X？反算时所求点的X坐标27.Y？反算时所求点的Y坐标28.结果29.FSZH=反算时，计算得出的所求点的里程30.FSJL=反算时，计算得出的所求点的边距（负为左偏，正为右偏）GCJSAbs(C-D)×R÷2→Z[1]: A-Z[1]→Z[2]: A+Z[1]→Z[3]:If D-C＞0：Then 1→E:Else -1→E:IfEnd←┚If K≤Z[2] ：Then Goto 1：Else If K＞Z[2] And K≤A: Then Goto 2: Else If K＞A And K≤Z[3]: Then Goto 3: Else IF K>Z[3]:then Goto 4:IfEnd:IfEnd:IfEnd:IfEnd←┚LbI 1： Z[2]-K→L: B-(Z[1]+L)×C→H: Goto 5←┚LbI 2: K-Z[2]→L: B-(Z[1]-L)×C+E×L2÷R÷2→H: Goto 5←┚LbI 3: Z[3]-K→L: B+(Z[1]-L)×D+E×L2÷R÷2→H: Goto 5←┚LbI 4: K-Z[3]→L: B+(Z[1]+L)×D→H:Goto 5←┚LbI 5:＂EL=＂：H◢←┚主线高程（ZXGC）LbI 0: ＂K＂?K: If K≤552.156（第1条竖曲线的终点下面同上）: Then Goto A：Else If K≤1352.805 And K＞552.156: Then Goto B:Else If K≤2045.311 And K＞1352.805: Then Goto C:Else If K≤2493.942 And K＞2045.311: Then Goto D：Else If K≤3265.244 And K＞2493.942: Then Goto E：IfEnd:IfEnd:IfEnd:IfEnd: IfEnd←┚LbI A: 501.436（边坡点桩号）→ A: 226.395（边坡点高程）→B: -0.0257（边坡点前坡比）→ C: -0.005（边坡点后坡比）→D:4900（半径）→R: Prog＂GCJS＂: Goto F←┚（下面同上）LbI B: 1303.305→A: 222.386→B: -0.005→ C: 0.04→D:2200→R: Prog＂GCJS＂: Goto F←┚LbI C: 1985.311→A: 249.666→B: 0.04→ C: -0.04→D:1500→R: Prog＂GCJS＂: Goto F←┚LbI D: 2425.578→A: 232.055→B: -0.04→ C: -0.0218→D:7500→R: Prog＂GCJS＂: Goto F←┚LbI E: 2849.694→A: 222.823→B: -0.0218→ C: -0.03→D:6000→R: Prog＂GCJS＂: Goto F←┚LbI F：If K≤309.562（第1条曲线的终点）: Then Goto 1：Else If K≤1053.022 （第2条曲线终点）And K＞309.562: Then Goto 2:Else If K≤1354.954 And K＞1053.022: Then Goto 3:Else If K≤1613.025 And K＞1354.954: Then Goto 4：Else If K≤1863.124 And K＞1613.025: Then Goto 5：Else If K≤2937.433 And K＞2312.560: Then Goto 6：Else If K≤3145.553 And K＞2937.433: Then Goto 7：IfEnd:IfEnd:IfEnd:IfEnd:IfEnd:IfEnd:IfEnd:←┚LbI 1: -0.015（直线段横坡）→ F: -0.02（圆曲线线段横坡）→G: 126.9（第1缓和曲线起点）→ I: 176.9（第1缓和曲线终点）→J:259.562（原曲线终点）→L: 309.562(第2缓和曲线终点)→M: 1(方向,左负右正)→N :Goto G←┚（下面同上）LbI 2: -0.015→ F: -0.02→G: 860.749→ I: 920.749→J:993.022→L: 1053.022→M:-1→N :Goto G←┚LbI 3: -0.015→F: -0.02→G: 1248.441→I: 1283.441→J:1319.954→L: 1354.954→M:+1→N :Goto G←┚LbI 4: -0.015→F: -0.02→G: 1439.731→I: 1489.731→J:1563.025→L: 1613.025→M:-1→N :Goto G←┚LbI 5: -0.015→ F: -0.02→G: 1712.165→ I: 1757.165→J:818.124→L: 1863.124→M: 1→N :Goto G←┚LbI 6: -0.015→F: -0.02→G: 2121.544→I: 2121.544→J:2312.560→L: 2312.560→M: 1→N :Goto G←┚LbI 7: -0.015→F: -0.02→G: 2752.223→I: 2812.223→J:2877.433→L: 2937.433→M: 1→N :Goto G←┚LbI 8: -0.015→F: -0.02→G: 3060.515→I: 3060.515→J:3145.553→L: 3145.553→M:-1→N :Goto G←┚LbI G:IfK≤I:Then Goto H :Else If K≤J And K＞I: Then Goto I: Else If K≤L And K＞J: Then Goto J: Else If K≤M And K＞L: Then Goto K: Else If K＞M: Then Goto L: IfEnd:IfEnd:IfEnd:IfEnd: IfEnd←┚LbI H:F→Z:F→Y: Goto 9LbI I:J=I=＞Goto J:（(K-I)× (Abs(F)+ Abs(G))）÷(J-I) -Abs(F)→Z[1]:If N=-1: Then Goto M: Else Goto N: IfEnd←┚LbI M:Z[1]→Y:If Z[1]≤-F:Then F→Z:Else –Z[1]→Z: IfEnd :Goto 9LbI N:Z[1]→Z:If Z[1]≤-F:Then F→Y:Else –Z[1]→Y: IfEnd :Goto 9LbI J:If N=-1:Then G→Z:-G→Y:Else -G→Z: G→Y: IfEnd: Goto 9←┚LbI K: M=L=＞Goto L: (M-K)×(Abs(F)+ Abs(G))÷(M-L)- Abs(F)→Z[2]: If N=-1: Then Goto P: Else Goto Q: IfEnd←┚LbI P:Z[2]→Y:If Z[2]≤-F:Then F→Z:Else –Z[2]→Z: IfEnd :Goto 9LbI Q:Z[2]→Z:If Z[2]≤-F:Then F→Y:Else –Z[2]→Y: IfEnd :Goto 9LbI L:F→Z:F→Y: Goto 9Lbl 9:＂ZC=1,YC=2＂?X:If X=1:Then Goto V:EIse Goto W: IfEnd Ibl V:” BJZ”?J:”BHZ=”:H+J×Z◢Goto 0Ibl W:” BJY”?S:”BHY=”:H+S×Y◢Goto 0显示K？输入桩号显示结果：EL= 所求点中桩高程显示“ZC=1,YC=2”？左侧边桩输1，右侧边桩输2显示“BJZ”或“BJY”输入左右边桩的距离显示结果：“BHZ=”或“BHY=”所求点左右边桩的高程。

卡西欧5800道路施工逐桩坐标计算程序说明：本程序包括了一个完整复合曲线的所有计算，其中需要输入的要素也是对应的完整复合曲线的要素，对不完整的曲线需要将要素进行必要的计算，转换后方能使用该程序，本程序已经包括了边桩几斜交边桩坐标的计算，参见要素说明.本程序是通过4800，4850程序的改编而来，详细程序如下：?D:?A:?R:?C:?N:?E:?F （输入计算要素，见后）C2÷6÷R－C4÷336÷R3+C6÷42240÷R5－R(1-COS（90C÷R÷∏）)→XC－C3÷40÷R2＋C5÷3456÷R4－RSin（90C÷R÷∏）－C7÷599040÷R6→Y90C÷∏÷R→B(R＋X)tan(Abs(A))÷2)＋Y→T(Abs(A)－2B) ∏R÷180＋2C→LD－T→GG＋C→HG＋L÷2→IG＋L－C→KG＋L→MA÷Abs(A) →SF＋A÷2＋90S→U(R＋X)÷cos(A÷2)→VN＋V cos(U)→BE＋V sin(U)→OLbl 1?Z 输入待求点桩号If Z≤GThenT＋G－Z→LF＋180→VF→UGoto 2IfEndIf Z≤HThenZ－G→LGoto 7IfEendThenF＋A÷2＋180＋90S＋(Z－I)×180S÷R÷∏→LL＋90S→UGoto 5IfEendIf Z≤MThenM－Z→LGoto 7IfEndIf Z＞MThenZ－M＋T→LF＋A→VV→UGoto 2IfEndLbl 2“X=”:N＋Lcos(V)→X ▲“Y=”:E＋Lsin(V)→Y ▲Goto 6Lbl 3“X=”:X＋Abs(J)cos(V)→X ▲“Y=”:Y＋Abs(J)sin(V)→Y ▲Goto 1Lbl 4“X=”:N＋Tcos(J)＋Vcos(Y)→X ▲“Y=”:E＋Tsin(J)＋Vsin(Y)→Y ▲Goto 6Lbl 5“X=”:B＋Rcos(L)→X ▲“Y=”:O＋Rsin(L)→Y ▲Goto 6Lbl 7L－L5÷(40 R2 C2)＋L9÷3456÷R4÷C4－L13÷599040÷（C6R6）→VL3÷（6RC）－L7÷336÷（C3R3）＋L11÷42240÷（C5R5）→JStan-1(J÷V) →U√（V2+J2）→VIf Z≤G+CThenGoto 8If Z≤MThenGoto 9IfEendLbl 8F+180→JF+U→YF+90SL2÷(RC∏)→UGoto 4Lbl 9F+A→JJ+180－U→YJ－90SL2÷(RC∏) →UGoto 4Lbl 6“BZ=”:?→J 开始计算边桩坐标，J为斜距，左侧输负，右侧输正90→V 斜交角度，如果不是正交改为：？→VIf J=0 如果不计算边桩坐标J输0直接返回到桩号输入ThenGoto 1IfEndIf J＜0ThenU－180＋V→VGoto 3IfEndIf J＞0ThenU+V→VGoto 3IfEnd要素说明：D—交点桩号，该交点为完整复合曲线的大交点，如果不满足条件需要转换A—转角值，路线右偏输正，左偏输负。

CASIO5800道路放样程序主程序文件名“XYNE”“XC”？→O：“YC”？→P：Lbl 0：“K”？→K：“W”？→W：Prog“AD”：Prog“AB”：（R+L2÷(24R)-L^(4)÷(2688R^(3))）tan(Abs(V)÷2)+L÷2-L^(3) ÷(240R2) →T:E-Tcos(A) →E : F-Tsin (A) →F : E+Tcos(A)+Tcos(A+V) →M: F+Tsin(A)+Tsin(A+V) →N: 兀R（Abs（V））÷180→D：K-U→Z：If K≤U Then E+Zcos（A）-Wsin（A）→X：F+Z sin（A）+W cos（A）→Y：“X”：X⊿“Y”:Y⊿Prog“12”：Goto0:ElseIf K＜U+L: Then Prog“1”： E+Gcos(A)-SHsin(A)-Wsin(A+SQ)→X : F+Gsin(A)+SHcos(A)+Wcos (A+SQ)→Y : “X”：X⊿“Y”:Y⊿Prog “12”：Goto0:ElseIf Z＜D: Then 90(2Z+L) ÷兀÷R→Q：Rsin(Q)+L÷2-L^(3) ÷240÷R2→G：R(1-cos(Q))+L2÷24÷R→H: E+Gcos(A)-SHsin(A)-Wsin(A+SQ) →X: F+Gsin(A)+SHcos(A)+Wcos (A+SQ)→Y : “X”：X⊿“Y”:Y⊿Prog“12”：Goto0:ElseIf Z＜D+L: Then D+L-Z→Z：Prog“1”：M-Gcos(A+V)-SHsin(A+V)-Wsin(A-SQ+V) →X: N-Gsin(A+V)+SHcos(A+V)+Wcos (A-SQ+V)→Y : “X”：X⊿“Y”:Y ⊿Prog“12”：Goto0:ElseIf Z＞D+L: Then K-(U+D+L)→Z： M+Zcos(A+V)-Wsin(A+V) →X: N+Zsin(A+V) +Wcos (A+V)→Y : “X”：X⊿“Y”:Y⊿Prog“12”：Goto0:Else子程序“1”Z-Z^(5)÷40÷R 2÷L 2→X: Z^(3) ÷6÷R÷L-Z^(7) ÷366÷R^(3) ÷L^(3) →Y: 90Z 2÷R÷L÷兀→Q: X→G:Y→H:子程序“AB”If V≤0 ：Then -1→S：Else If V＞0 ：Then 1→S：Else显示程序“12”PoL（X-O，Y-P）→I：I：J：“I”：I⊿If J≤0：Then 360+J→J：“J”：J⊿矩阵程序“2”Mat B[1，1]→E: Mat B[1，2]→F:Mat B[1，3]→V:Mat B[1，4]→R: Mat B[1，5]→L: Mat B[1，6]→U: Mat B[1，7]→A:数据库嵌入程序“AD”If K≤ZH桩号：Then Prog“A”：Goto 0：Else If K≤ZH桩号：Then Prog“B”：Goto 0：Else If K≤ZH桩号：Then Prog “C”：Goto 0：Else If K≤ZH桩号：Then Prog“?”：Goto 0：………… Else Lblo：数据库“A”If K≤ZH桩号：Then [[交点X坐标，交点Y坐标，转角V （左-,右+）,半径R,缓和曲线L,直缓桩号U,方位角A]] →Mat B:Porg “2”:Else If K≤ZH桩号：Then [[交点X坐标，交点Y坐标，转角V（左-,右+）,半径R,缓和曲线L,直缓桩号U,方位角A]] →Mat B:Porg “2”:Else。

5800计算程序本程序使用高斯四节点法为计算内核,可以计算任意里程偏距坐标,反算里程和边距,计算任意点高程,分正常计算和设计计算,支持多条线路.支持两种以中线为旋转轴的超高方式,支持加宽,支持带绿化带,支持开口线和坡脚线计算.本程序仅提供大家参考,计算器主程序：ZHU-CHENG-XULbi0:“1，ZS=FS，2ZS，3FS，4XY＝＞SG，5。

”？U:U=1＝＞Prog”ZS=FS”:U=2＝＞Porg”ZS”:U=3＝＞Porg”FS”: U=4＝＞Porg”XY＝＞SG”:…………Goto0://计算器总的主程序，进入选择各种分支计算程序。

1为公路三维计算，2为普通正算，3为普通反算，4为大地坐标转施工坐标。

选择错误重新选择。

公路三维部分正反算选择程序：ZF=FS“1LZ＝＞XY，2XY＝＞LZ,3BIANPO-FY”？U： //正反算选择，正算选1，反算选2If U=1: ThenProg”ZS-XH”:IfEnd://进入正算循环主程序IfU=2:Then Prog”FS-XH”:IfEnd: //进入反算循环主程序IfU=3:ThenProg”BIANPO-FY”:IfEnd://进入边坡开挖主程序正算循环主体程序ZS-XH“1PT-2SJ”？W：//普通计算和设计边距计算选择“XL-XZ“？U //选择线路1～NIfW=1: Then //运行普通计算模式边踞自由输入Lbi0:”L=”?L:”Z=”?Z: //输入桩号和边距Porg”ZS-XH-1“：Goto0:IfEnd:运行正算循环子程序IfW=2:Then //运行设计计算模式边距按设计输入Lbi1:”L=”?//输入桩号Porg”0。

5LUKUAN”://录入公路加宽后，左右侧设计边距Porg”ZS-XH-1“：Goto1:IfEnd:运行正算循环子程序正算循环子程序ZS-XH-1Prog”PM-XL-XZ”: //进入平面线路选择程序对线路线元进行选择Prog”CHAOXIAN-CL” //里程超限处理Prog”GL-ZS”: //进入坐标正算程序Prog”GL-BZ-ZS” //进入公路边桩正算程序Prog”GC-XL-XZ”: //进入高程线路选择程序对线路线元进行选择Prog”GC”://进入高程部分提取中桩高程IfW=1:ThenProg”B-H”://普通计算时进入边桩高程部分提取边桩高程ElseZ[3] →Z: //按设计计算时候代入减去半边绿化带后的路基边距Prog”ZS-XS”: //进入正算显示程序，提取正算三维结果反算循环主体程序FS-XH“XL-XZ”？U：//选择线路 1～NLbio：“L”？L：“X=”？T：“Y=”？S： //输入起算假定桩号，实测X YProg”GL-FS”: /进入公路反算程序Prog”FS-XS”:Goto0: //进入反算显示程序坐标正算程序GL-ZS.1739274226→I:.3260725774→J:.0694318442→K:.3300094782→E //四节点法计算系数录入1-E→F：1-K→M：Q+KrO(C+KD) →K: 数据库要用X、Y、Q作为起点坐标方位角Q+ErO(C+ED) →E:Q+FrO(C+FD) →F:Q+MrO(C+MD) →M:X+O(ICos(K)+JCos(E)+JCos(F)+ICos(M)) →X: //计算出X坐标Y+O(ISin(K)+JSin(E)+JSin(F)+ISin(M)) →Y: //计算出Y坐标Q+Or(C+D) →Q： //计算出算点方位角边桩坐标计算程序ZS-BZX+Zcos(Q+90) →G //计算出边桩X坐标Y+ZSin(Q+90) →M //计算出边桩Y坐标公路反算程序GL-FS20→V: //V值设置为20，为第一次反算控制Do:Porg”PM-XL-XZ”://进入平面线路选择程序Porg:”CHAOXIAN-CL”://里程超限处理程序IfV＞10: // V越小反算速度可能越慢，但是半径过小有可能算不出来，10为建议值Then Porg”GL-JY-ZS”://进入公路简易正算二节点程序Else Porg”GL-ZS”: //进入公路正算4节点程序IfEnd:T-X→I:S-Y→J: //试算点与测点X，Y距离差Pol(I,J):J-Q→J: //方位角转换为接近试算点到测点的值：Rec(I,J): //计算出X差和Y差（X里程，J边）I+(1-CJ-2DJ) →V //计算出里程大致差距（此公式为加速反算）L+V→L：//修改后的里程Abs(V) →V: //修改值的绝对值LpWhile V＞1m: //小于一毫米通过，大于一毫米循环计算中桩高程计算程序GCIf I-J＜0:Then1→H:IfEnd: //开口向上部分赋值IfI-J=0:ThenGoto1: IfEnd: //无竖曲线部分赋值IfI-J＞0:Then-1→IfEnd: //开口向下时赋值(L-C)I+K+(L-C)^2/(2E)H→I：//计算得出竖曲线部分中桩高程Goto2:Lbi1:K+(L-C)I→I: Lbi2: //计算得出无竖曲线部分中桩高程边桩高程计算程序B-H（L-C）/（D-C）→P: //计算点占本线元长比例3P^(2^(S))-2P^(3^(S)) →P: //计算得线性和高次抛物线系数F+P（T-F）→E： //初步计算计算点横坡If V=1:Then If N= -1:Then //直线和直缓，缓直,曲线左偏时Z＜0：ThenIfE＞Z[4]:ThenZ[4]→E: IfEnd : //计算左侧时，横坡大于一般横坡取一般横坡I+Abs(Z)E→O: IfEnd : //左侧高程IfZ＞0:Then I+ZF→O:IfEnd: //右侧高程得数If N=1:Then IfZ＜0:Then I+Abs(Z)E:IfEnd: //曲线右偏时，左侧高程得得IfZ＞0:Then If E＞Z[4]:Then Z[4]→E:IfEnd //右侧高程横坡大于一般横坡，取一般横坡I+ZE→O:IfEnd: IfEnd: IfEnd: //右侧高程得数IfV=2:ThenIf Z＞0:Then I+ZE→O:IfEnd://圆曲线全超高和S曲线超高渐变部分右侧高程得数IfZ＜0:ThenI-Abs(Z)E→O:IfEnd:IfEnd: // 圆曲线全超高和S曲线超高渐变部分左侧高程得数设计标准路面宽度子程序 CHANGSHU12→Z[1]：//选择12米为路面半宽度，根据设计修改1.5→Z[2]：//选择1.5米为绿化带半宽-0.02 →Z[4] //一般横坡录入路面宽度子程序 0。

卡西欧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。

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：待求点里程。

卡西欧5800测量实用程序(含直线、圆曲线、大地转施工程序)一、直线计算程序0→I:0→J:”X0”?D:”Y0”?E:”X1”?B:”Y1”?C:Pol(B-D,C-E):J →A:If A<0:Then A+360→A:Else A→A:IfEnd:Lbl 0:?O:?S:If O≠0:Then Goto 1:IfEnd:D+Scos(A) →X:”X=”:X ◢E+Ssin(A) →Y:”Y=”:Y ◢Goto 0:Lbl 1:D+Scos(A)+Ocos(A+90) →X:”X=”:X ◢E+Ssin(A)+Osin(A+90) →Y:”Y=”:Y ◢Goto 0说明：输入程序时注意区别字母O 与数字0程序运行时符号说明X0? Y0? 分别输入直线起点的XY 坐标值X1? Y1? 分别输入直线终点的XY 坐标值O? 输入边桩与中桩的距离( 左边桩为负值, 右边桩为正值), 如计算中桩坐标输入0S? 输入所求点到直线起点的距离二、圆曲线计算程序0→I:0→J: ”X0”?C:”Y0”?D:”X1”?E:”Y1”?F:?R:”L:-1 R:1”?N:”ZY”?W:Lbl 0:”LN”?T:T-W→O:Pol(E-C,F-D):I→S:J→A: If A<0:Then A+360→A:Else A →A:IfEnd:sin-1(S÷(2R)) →K:2∏RK÷180→L:180O÷(2∏R) →G:(2R)sin(G) →H:C+Hcos(A-KN+GN) →X:”X=”:X ◢D+Hsin(A-KN+GN) →Y:”Y=”:Y ◢R-0.5√(4R2-H2) →Q:”Q=”:Q ◢“S=”?V:If V=0:Then Goto 0:IfEnd:”L:-90 R:90”?U:X+Vcos(A-KN+2GN+U) →X:”X=”:X ◢Y+Vsin(A-KN+2GN+U) →Y:”Y=”:Y ◢Goto 0说明：输入程序时注意区别字母O 与数字0程序运行时符号说明X0? Y0? 分别输入直线起点的XY 坐标值X1? Y1? 分别输入直线终点的XY 坐标值R? 输入圆曲线半径L:-1 R:1? 圆曲线向左转弯时输入-1, 向右转弯时输入1ZY? 输入起点桩号LN? 输入所求点桩号( 按下exe 键后显示中桩坐标)S? 输入所要求的边桩与中桩的距离L:-90 R:90 要求左边边桩输入-90, 右边边桩输入90( 按下exe 键后显示边桩坐标)三、大地转施工“QK=”？Z:“QX=”？A:“QY=”?B:“ZX=”？O :“ZY=”？Q ↙Pol( O -A,Q-B) ↙J<0=>J+360→J ↙J>360=>J-360→J ↙J→F:“JF=”：F ◆DMS ◢Lb1 0 ↙？X:?Y ↙Pol （X-A,Y-B):I→D ↙J<0=>J+360→J ↙J>360=>J-360→J ↙J-F→T:T→G ↙T<0=>T+360→G ↙T>360=>T-360→G ↙Z+Dcos （G)→K:“K=”：K ◢Dsin （G)→C:“C=”：C ◢GoTo 0 ↙说明：QK-- 输入线路起始点里程QX-- 输入线路起始点中桩X 坐标QY-- 输入线路起始点中桩Y 坐标ZX-- 输入线路前进方向中桩X 坐标ZY-- 输入线路前进方向中桩Y 坐标JF=-- 显示线路测量坐标方位角X-- 输入线路附近任意点测量X 坐标Y-- 输入线路附近任意点测量Y 坐标K-- 显示计算结果：里程桩号，即工程X 坐标C-- 显示计算结果：左右偏距，即工程Y 坐标, 中桩=0 ，左偏= 负值，右偏= 正值。

卡西欧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。

一．使用说明：1.规定：（1）以道路中线的前进方向区分线路的左右偏转方向，线路左偏Q=-1，线路右偏Q=1，直线Q=0。

（2）当所求点位于线路中线时，Z=0；当位于左侧时，Z取负值；当位于右侧时，Z取正值。

（3）当线元为直线时，其起终点的曲率半径为无穷大，取10的45次方代替。

（4）当线元为圆曲线时，起终点曲率半径为圆曲线半径。

（5）当线元为缓和曲线时，起终点曲率半径与直线相接取10的45次方代替，与圆曲线相接取圆曲线半径。

2.输入与显示说明：（1）输入部分：N？：选择计算方式，输入1表示由里程边距计算坐标，输入2表示由坐标计算里程边距。

X0？线元起点X坐标。

Y0？线元起点Y坐标。

S0？线元起点里程。

F0？线元起点切线方位角。

LS？线元长度。

R0？线元起点曲率半径。

RN?线元终点曲率半径。

Q？线元左右偏标志（左偏=-1，右偏=1，直线=0）S?正算时所求点里程。

Z？正算时所求点距中线边距。

位于线路左侧输入负值，在右侧输入正值。

AR?斜交角度。

（为线路右角，在0-180度之间）X?反算时所求点X坐标。

Y？反算时所求点Y坐标。

（2）显示部分：XS=计算结果，所求点X坐标。

YS=计算结果，所求点Y坐标。

FS=计算结果，所求点切线方位角。

S=计算结果，反算时所求点里程Z=计算结果，反算时所求点边距。

二．程序清单1（用于任意曲线元的坐标计算）1.主程序(TYQXJS)″1.SZ=>XY″:″2.XY=>SZ″:?N:″X0″?U:″Y0″?V:″S0″?O: ″F0″?G:″LS″?H:″R0″?P:″RN″?R:?Q:1÷P→C:(P-R) ÷(2HPR) →D:180÷π→E:If N=1: THEN Goto 1:Else Goto 2:IfEnd:Lbl 1:?S:”D(Z-，Y+)=”?Z:90→T:”O=”?T:Abs(S-O) →W:Prog ″SUB1″:″XS″:X◢″YS″:Y◢″FS″:F◢Goto 1:Lbl 2:?X:?Y:I→X:J→Y:Prog ″SUB2″:″S″:O+W→S:◢″Z″:Z◢Goto 22.正算子程序(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))) →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) →F:X+Zcos(F+T) →X:Y+Zsin(F+T) →Y3.反算子程序(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 Abs(Z)<1E-6 :Then Goto 1:Else W+Z→W:Goto 0 : IfEnd:Lbl 1:0→Z:Prog ″SUB1″:(J-Y)÷sin(F)→Z 三．用于提前将曲线元要素输入子程序，计算时只输入里程桩号与边距。