CASIO5800坐标计算程序
CASIO5800坐标正反算计算程序
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列表数组中，要想多次调用最好随公式取出结果，并赋给变量。
CASIO5800计算器中边桩坐标计算程序
CASIO5800计算器中边桩坐标计算程序一．GSGL（主程序）1.“K0”? W:“X0”? U:“Y0”? V:“R”?R: “LH”?L:“A0”?A:“B”?B:“R=1,L=-1”?N:“S-BACK"?E:“T-BACK"?F:2.Lb1 0:3.“KI”? K: E-F S:4.(R+L2÷(24R)-L^(4)÷(2688R^(3)))×tan(B÷2)+L÷2-L^(3)÷(240R2)T:5.S-T H:K-W J:J＜0 =〉Goto 0:6.If H>J: Then U+Jcos(A)X:V+Jsin(A)Y:7.“XI=”:X “YI=”:Y8.A G:Prog “SUB”: Goto 0: IfEnd:9.L=0 =〉Goto 1:10.K-H-W Z: Z-Z^(5)÷(90R2L2)C:11.30 Z2÷(πLR)D: Z＞L =〉Goto 1:12.U+Hcos(A)+Ccos(A+ND)X:V+Hsin(A)+Csin(A+ND)Y:13.“XI=”:X“YI=”:Y14.A+180NZ2÷(2πLR)G:Prog “SUB”: Goto 0:15.Lb1 1:16.πRB÷180-L C: K-H-W-L M: 90M÷(πR)I: 2Rsin(I)O:17.L-L^(3)÷(90R2)P:30LN÷(πR)Q:A+3Q+NI D: M>C =〉Goto 2:18.U+Hcos(A)+Pcos(A+Q)+Ocos(D)X: V+Hsin(A)+Psin(A+Q)+Osin(D)Y:19.“XI=”:X “YI=”:Y20.A+N(180L÷2÷π÷R+180M÷(πR))G: Prog “SUB”: Goto 0:21.Lb1 2:22.L=0 =〉Goto 0:23.A+NB G:24.H-J+L+πRB÷180P:P<-.003 =〉Goto 0:25.P-P^(5)÷(90R2L2)C: 30P2÷(πRL)D:A+NB-180-ND Q:26.U+Scos(A)+Tcos(G)+Ccos(Q)X: V+Ssin(A)+Tsin(G)+Csin(Q)Y:27.“XI=”:X “YI=”:Y28.A+NB-180NP2÷(2πRL)G: Prog “SUB” :Goto 0:二．SUB（子程序）1．Lb1 0:2. 0I：0O：0M:3. “Y=1，N=0”? I：I≠1 =〉Goto 2：4. G+90C：G-90D：“R=1，L=-1”? O：“S-LONG”? M：5. If O=1：Then X+Mcos(C)P : Y+Msin(C)Q：6. “XR=”：P “YR=”：Q7. Else X+Mcos(D)P: Y+Msin(D)Q：8. “XL=”：P “YL=”：Q IfEnd:9. Goto 0：10. Lb1 2说明1）K0：起点的桩号2）X0：起点的X坐标3）Y0：起点的Y坐标4）KI：需计算的桩号5）R：圆曲线半径6）LH：缓和曲线段长度7）A0：起点方位角（角度）8）B：偏角（角度）9）R=1，L=-1：表示路线左右偏，右1，左-110）S-BACK：起点至交点长度11）T-BACK：一般取012）Y=1，N=0：是否计算边桩，是1，否013）R=1，L=-1：左右边桩，右1，左-114）S-LONG：表示中桩至边桩距离。
CASIO-5800测量程序
CASIO-5800测量程序CASIO-5800测量程序5800中边桩任意点坐标计算程序一、文件名：HZ ZH 输入文件名Deg:Fix 4 基本设置"FWJ∠α="？Q8起点切线方位角"ZX(m)="？O: "ZY(m)= "？P8 置镜点坐标"HX(m)= "？M: "HY(m)= "？N 8 后视点坐标Pol(A－O,B－P) 8If J<0:Then J＋360"H:Else J"H:IfEnd8"QX(m)="?U: "QY(m)="?V: "QDK(m)="?G8 起算点坐标和里程Ldl 08"GDK(m)="?Z 8 观测点里程G－Z"C8"GX(m)=":U＋cos(Q) ×C"Ay 置镜点至中桩X坐标"GY(m)=":V＋sin(Q) ×C"By 置镜点至中桩Y坐标Pol(A－O,B－P) 8 置镜点到中桩方位角J－H"J8If J－H<O:Then J＋360"J:Else J"J:IfEnd8"SPJ∠α=":J}DMSy显示顺拔夹角"I(m)=":Iy 置镜点至中桩距离Ldl 18"BJ(m)= "?S 8 边距If S=0:Then Goto 0:IfEnd8"JIA JIAO L(-),R(+)="?E8Q＋E"W8"GX(m)=":A＋cos(W) ×S"Xy 显示观测点X 坐标"GY(m)=":B＋sin(W) ×S"Yy 显示观测点Y 坐标Pol(X－O,Y－P) 8If J<0:Then J＋360"J:IfEnd8J－H"J8 置镜点到观测点方位角If J<0:Then J＋360"J:Else J"J:IfEnd8"SPJ∠α=":J}DMSy显示顺拔夹角"I(m)=":Iy 显示置镜点到观测点距离Goto 1 回到前面重复下一观测点二、文件名：ZH HY 缓和曲线段内任意点坐标计算程序“Deg:Fix 4“FWJ∠a =”？Q“ZX(m)=”?O:”ZY(m)=”?P“HX(m)=”?M:”HY(m)=”?N Pol(M-O,N-P)If J<0:Then J+360→H：Else J→H:IfEnd “HWJ∠a =”:H?DMS◢“ZHUAN XIANG L（-1），R（1）=”？T“R（m）=”?R“LS(m)=”?D“KZH=”?G“ZH X(m)=”?U“ZH Y(m)=”?VLbI 0“QDK=”?ZAbs(Z-G)→CC-C^5/(40R^2D^2)+C^9/(3456 R^4D^4) →WC^3/(6RD)-C^7/(336R^3D^3)+C^11/(42240R^5D^5)→KQ+(C^2/(2RD)×(180/Π)) →FIf J<0:Then F+360→F：Else F→F:IfEnd “QFWJ ∠a =”:F?DMS◢“QS X（m）=”：U+W×Cos(Q)+K×Cos(Q+90) →A: A ◢“QS Y（m）=”：V+W×Sin(Q)+K×Sin(Q+90) →B: B ◢Pol(A-O,B-P)If J<0:Then J+360→J：Else J→J:IfEnd J-H→JIf J<0:Then J+360→J：Else J→J:IfEnd“SPJ∠a=”:J?DMS◢“I（m）=”:I◢LbI 1“BJ(m)=”?SIf S=0:Then Goto 0: IfEnd“JIA JIAO L(-),R(+) =”?E“BZ X=”：A+SCos(F+E) →X: X ◢“BZ Y=”：B+S Sin(F+E) →Y：Y ◢Pol(X-O,Y-P)If J<0:Then J+360→J：Else J→J:IfEnd J-H→JIf J<0:Then J+360→J：Else J→J:IfEnd “∠BIAN ZHUANG PJ=”:J?DMS◢“I（m）=”:I◢Goto 0三、文件名：HY YH 园和曲线段内任意点坐标计算程序“Deg:Fix 4“FWJ∠a=”？Q起算点方位角“ZX(m)=”?O:”ZY(m)=”?P置镜点坐标“HX(m)=”?M:”HY(m)=”?N后视点坐标Pol(M-O,N-P)If J<0:Then J+360→H：Else J→H:IfEnd “HFWJ ∠a=”:H?DMS◢后视方位角“ZHUAN XIANG L（-1），R（1）=”？T曲线转向左负右正“R（m）=”?R曲线半径“LS(m)=”?D缓和曲线长度“KZH(m)=”?G曲线起点ZH里程“XZH(m)=”?U曲线起点X“YZH(m)=”?V曲线起点YLbI 0“QDK(m)=”?Z观测点里程Abs(Z-G)→CR×Sin(180×(C- D）/（πR）+90D/（πR）)+D/2-D^3/(240R ^2)→WR×（1-Cos(180(C-D）)/（πR）+90D/（πR））+ D^2/(24R)- D^4/(2688R ^3)→KQ+(180（C-D）/（πR）+90D/（πR）)→F If F<0:Then F+360→J：Else F→J:IfEnd “Q FWJ=”：F?DMS◢观测点切线方位角“QS X(m)=”：U+√(W²+K²)×Cos((Q)+tan?¹(K÷W)) →A: A ◢观测点X坐标“QS Y(m)=”：V+√(W²+K²)×Sin((Q)+tan?¹(K÷W)) →B: B ◢观测点Y坐标Pol(A-O,B-P)If J<0:Then J+360→J：Else J→J:IfEndJ-H→JIf J<0:Then J+360→J：Else J→J:IfEnd“SPJ∠a=”:J?DMS◢观测点所拨偏角“I（m）=”:I◢置镜点到观测点距离“BIAN JU(m)=”?S边距，左负右正If S=0:Then Goto 0: IfEnd“JIA JIAO L(-),R(+) =”?E边桩至中桩切线方向的夹角左负右正“BZ X(m)=”：A+SCos(F+E) →X: X ◢“BZ Y(m)=”：B+S Sin(F+E) →Y：Y ◢Pol(X-O,Y-P)If J<0:Then J+360→J：Else J→J:IfEnd J-H→JIf J<0:Then J+360→J：Else J→J:IfEnd “BZ ∠a =”:J?DMS◢“I（m）=”:I◢Goto 0四、文件名：YH HZ 另一端缓和曲线段内任意点坐标计算程序“Deg:Fix 4“FWJ∠a=”？Q起算点方位角“ZX(m)=”?O:”ZY(m)=”?P置镜点坐标“HX(m)=”?M:”HY(m)=”?N后视点坐标Pol(M-O,N-P)If J<0:Then J+360→H：Else J→H:IfEnd “HFWJ ∠a=”:H?DMS◢后视方位角“ZHUAN XIANG L（-1），R（1）=”？T曲线转向左负右正“R（m）=”?R曲线半径“LS(m)=”?D缓和曲线长度“KZH(m)=”?G曲线起点ZH里程“XZH(m)=”?U曲线起点X“YZH(m)=”?V曲线起点YLbI 0“QDK(m)=”?Z观测点里程Abs(G-Z)→CC-C^5/(40R^2D^2)+C^9/(3456 R^4D^4) →WC^3/(6RD)-C^7/(336R^3D^3)+C^11/(42240R^5D^5)→KQ-(C^2/(2RD)×(180/Π)) →FIf F<0:Then F+360→F：Else F→F:IfEnd If F>360:Then F-360→F：Else F→F:IfEnd “QFWJ∠a =”：F?DMS◢“QS X(m)=”：U-W×Cos(Q)+K×Cos(Q+90) →A: A ◢“QS Y(m)=”：V-W×Sin(Q)+K×Sin(Q+90) →B: B ◢Pol(A-O,B-P)If J<0:Then J+360→J：Else J→J:IfEndJ-H→JIf J<0:Then J+360→J：Else J→J:IfEnd“SPJ∠a=”:J?DMS◢“I（m）=”:I◢LbI 1“BJ(m)=”?SIf S=0:Then Goto 0: IfEnd“JIA JIAO L(-),R(+) =”?E“BZ X=”：A+SCos(F+E) →X: X ◢“BZ Y=”：B+S Sin(F+E) →Y：Y ◢Pol(X-O,Y-P)If J<0:Then J+360→J：Else J→J:IfEnd J-H→JIf J<0:Then J+360→J：Else J→J:IfEnd “BZ∠a =”:J?DMS◢“I（m）=”:I◢Goto 0。
卡西欧fx-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 （直线）。
CASIO fx-5800计算器测量坐标程序
主程序坐标正反算程序名称：001XYLB1 A↙Cls : Fix 4 : 30→Dim Z :"v=1,v=0"?Z ↙If Z=1:Then Goto B:IfEnd↙(Z=1进入里程点坐标正算）If Z=0:Then Prog"FSLCBZ ": Goto A:IfEnd↙(Z=0进入反算里程边桩）Lb1 B ↙“K=”?K ：（计算里程）Prog"SJK1 ":Prog”ZBFY” :Goto A↙子程序数据库 SJK1IF 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 : Prog”XLZBJSCX”:Return: IFEND↙IF…………Prog”XLZBJSCX”:Return: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] :E →T:Return:IFEND子程序坐标正算，名称：XLZBJSCXLB1 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 : ELSE E-T →T : T<0=>360+T→T : IFEND ↙GOTO 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 : ELSE 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*S IN 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 ↙LLb1 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↙子程序反算里程边桩，名称：FSLCBZLb1 1↙"POINT=X"?X:" POINT=Y"?Y↙（输入任意测点的XY坐标）:“K=”?K （试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）↙Lb1 2↙Prog"SJK1 ":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"SJK1":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]:"W(-Z,Y+)=":I◢(偏距）（由于该程序不能准确判断边桩左右方向，暂作修改，取消左右边判定）"K=":K+S→K◢（里程）Return↙子程序坐标放样：ZBFYLB1 0 ↙"P(N)="?M:"P(E)="?N ↙“a=”?P◢(输入边桩与线路夹角，左-右+)“W=”?D◢（输入边桩距）Z[18]+D*COS(T+P) →Z[20] ↙Z[19]+D*SIN(T+P) →Z[21] ↙“XP=”: Z[20] ◢（边桩X）“YP=”: Z[21] ◢（边桩Y）Pol(Z[20]-M,Z[21]-N):"HD=":I◢（边桩放样距）J<0=>J+360→J:"HAR=":J→DMS◢（边桩放样方位角）Return↙直线短链要将该直线分成二条直线来对待编制数据库文件即可，曲线长链要单独建判定语句子程序并然后编制二个数据库（按来向里程和去向里程分别推出曲线起终点里程，以此编制曲线段的数据库）。
卡西欧计算器5800程序
1、ZXKW（直线开挖）“D=”?D:“PD=”?F:“ZXJ=”?L:Lbl 0:“Y=”?Y:“Z=”?Z“T=”:Abs(Y)-((Z-D)*F+L)→T◢“ZL=”:Z-(D+0.3)◢GoTo 0本程序适用于直线段渠道开挖，坐标为工程坐标，D为坡脚处高程，PD为坡度，ZXJ为待放样边坡坡脚距中心线的水平距离，T表示棱镜所应水平移动距离，T为正值时向中心线方向移动，负值时向中心线外移动。
“ZL=”表示为垂直向下开挖深度。
Lbl 0前各数值为定量，后各数值为变量，Lbl 0后的0为01234的0.2、ZXDM(直线段面)“M=”?M:“N=”?N:“Q=”?Q:“O=”?O:L bl 0:“X=”?X:“Y=”?Y:“Z=”?Z“F=”(X-M)COS(Q)+(Y-N)SIN(Q)+O→F◢“G=”-(X-M)SIN(Q)+(Y-N)COS(Q) →G◢“T=”:Abs(G)-((Z-96)*2+8.88-0.67) ◢GoTo 0本程序适用于直线段渠道开挖，坐标为大地坐标，M、N为某段直线起点在大地坐标中的大地坐标值，O为直线起点处桩号值，F为桩号，G为中心距，T表示棱镜所应水平移动距离，T为正值时向中心线方向移动，负值时向中心线外移动。
3、ZJZB(直角坐标)“A=”?A:“B=”?B:Lbl 0:“C=”?C:“D=”?DPol(C-A,D-B)“L=”:I→L◢J≤0 ＝＞ J+360°→J“Q=”:J→Q◢GoTo 0本程序适用于坐标反算，A、B为起算点坐标值，C、D为待定点坐标值，L为起算点与待定点间水平距离，Q为起算点与待定点间的方位角。
4、JZB(极坐标)“A=”?A:“B=”?B:Lbl 0:“L=”?L: “Q=”?QRec(L,Q)“X=”:I+A→X◢“Y=”:J+B→Y◢GoTo 0本程序适用于坐标正算，A、B为起算点坐标值，L为起算点与待定点间水平距离，Q为起算点与待定点间的方位角，X、Y为待定点坐标。
卡西欧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 三．用于提前将曲线元要素输入子程序，计算时只输入里程桩号与边距。
casio-fx5800程序
卡西欧fx5800程序此程序只需输入桩号、中边夹角、中边距离、即可计算全线的中边桩坐标、高程。
（线形为对称形，非对称行因为没时间研究一直搁浅，其实只是加上Ls1、Ls2、T1、T2等要素，也可自行研究）1、主线坐标计算主程序A01ZB (子程序名B0ZB)8→DimZ （K待求点桩号）Lbl 9:Norm 1:”K”?K:Prog”B0ZB”:F-T→A:A+S→B:A+W→D: D-S→C If K<C:Then U→H:J→J:-1→O:Else V→H:-J→J:1→O:Ifend:N+OTcos(H)→Z[1]:E+OTsin(H)→Z[2]:Goto 0:Lbl 0 （K待求点桩号）K<A=>Goto 1:K<B=>Goto 2:K<C=>Goto 3:K<D=>Goto 4:K≥D=>Goto 5:Lbl 1:0→Z[3]:H+JZ[3]→Z[4]:(K-A)cos(H)→Z[7]:(K-A)sin(H)→Z[8]: Goto 6Lbl 2:30(K-A)^2/R/S/→Z[3]:H+3JZ[3]→Z[4](K-A)-(K-A)^(5)/90/R^2/S^2→ZZcos(H+JZ[3])→Z[7]:Zsin(H+JZ[3])→Z[8]:Goto 6Lbl 3:90S/R/+180(K-B)/R/→Z[3]:H+JZ[3]→Z[4]Rsin(Z[3])+S/2-S^(3)/240/R^2→IR(1-cos(Z[3]))+S^2/24/R→ZIcos(H)-JZsin(H)→Z[7]:Isin(H)+JZcos(H)→Z[8]:Goto 6Lbl 430(K-D)^(2)/R/S/→Z[3]:H+3JZ[3]→Z[4](D-K)-(D-K)^(5)/90/R^(2)/S^(2)→Z-Zcos(H+JZ[3])→Z[7]:-Zsin(H+JZ[3])→Z[8]:Goto 6Lbl 50→z[3]:H-JZ[3]→Z[4]:(K-D)cos(H)→Z[7]:(K-D)sin(H)→Z[8]:Goto 6 Lbl 6 边桩坐标Norm 1:”HR=”显示中线的方位角(度分秒)Fix 3:”X=”:Z[1]+Z[7]→“Y=”:Z[2]+Z[8]→Norm 1:”HJ”?Q:”DL”?P:Fix 3:”XL=”:X+Pcos(Z[4]-Q)→Z[5]“YL=”:Y+Psin(Z[4]-Q)→Z[6]Norm 1:”DR”?LFix 3:”XR=”:X+Lcos(Z[4]-Q+180)→”YR=”:Y+Lsin(Z[4]-Q+180)→Norm 1:”10x K(+,-)”?G:K+10G→K:Goto 9主程序完子程序(数据库，这是以哈大客运专线长春的一段为例，标段里程DK751000-DK775+036.5) B0ZBIf K≤754068.16:Then 750863.625→F:4912830.447→N:504790.687→E:-1→J:11000→R:370→S:911.477→T:1820.783→W:26°03′28.19″→U:18°30′4.05″→V:IfendIf K>754068.16:Then 755230.815→F:4916973.988→N:506177.187→E: 1→J:10000→R:430→S:1162.655→T:2319.528→W:18°30′4.05″→U:29°19′38.38″→V:IfendIf K>********:Then ********→F: ********→N: ********→E: ********→J: ********→R: ********→S: ********→T: ********→W:**°**′**″→U:**°**′**″→V:Ifend（If K>本交点的缓直点:Then 下一交点的桩号→F: 下一交点X坐标→N: 下一交点Y坐标→E: 转向→J: 半径→R: 缓和曲线全长→S: 切线→T: 曲线全长→W:前方位角→U:后方位角→V:Ifend）按此，全线要素可全部输入算例要素表与坐标表见附表0123456789 数字ABCDEFGHIJKLMNOPQRSTUVWXYZ 字母⊿转移结束符号DimZ、=>、→、DMS 、Lbl、If、Then、Else、Ifend、Goto、Norm、Fix均为一个整体2、高程主程序A02GC (子程序名B0GC 竖曲线要素，HP1、HP2横坡数据，V=-1为左半幅横坡，其他为右半幅)2→DimZLbl 0:?K:Prog”B0GC”:.005RAbs(I-L)→T:Z-T→A:Z+T→B:K-A→C:B-K →D:K-Z→EIf L>I:Then 1→J:Else L<I=>-1→J:IfendIf K<A:Then 0→C:Else K>B=>0→D:IfendIf K≦Z:Then I→Y:C→G:Else K>Z=>L→Y:D→G:IfendLbl 1:Fix 3:”H”:H+EY/100+JG^2/2/R→XNorm 1:”HD”?U:”Z=-1,Y=1”?V:If V<0:Then Goto 2:Else Goto 3:IfendLbl 2:Prog“HP1”:”Hi”:OGoto 4Lbl 3:Prog”HP2”:”Hi”:OGoto 4Lbl 4:Norm 1:”D”?P:Fix 3:”H1”:X-U+.01OP→“10x K(+,-)”?W:K+10W→K:Goto 0主程序完子程序数据库B0GCIf K≦753225:Then 749700→Z:181.776→H:30000→R:0→I:0.3→L:IfendIf K>753225:Then 753300→Z:192.576→H:30000→R:0.3→I:0.8→L:IfendIf K>下一竖曲线的前变坡点的桩号:Then下一竖曲线交点里程→Z: 下一竖曲线交点高程→H: 下一竖曲线半径→R:前纵坡→I:后纵坡→L:Ifend可输入全线竖曲线要素子程序（横坡数据库）HP1If K≦753225:Then -2→O: IfendIf K>753225:Then -2+(-1/100)*(K-753225)→O: IfendIf K>753325:Then -3→O: IfendIf K>******:Then *****→O: IfendIf K>横坡变坡点桩号:Then 横坡或者公式→O: Ifend高程程序就不举例子了3、卵形曲线（辛普生法则，匝道及主线均可算）此程序为廊坊至涿州段高速路面标的算例。
卡西欧5800完整程序
卡西欧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=”所求点左右边桩的高程。
CASIO5800计算器程序(JXT)
CASIO5800坐标正反算（线元法）LXJS(主程序)“1.SZ=>NE,2.NE=>SZ”?Q↙LbI 0:Q=1=>Goto 1: Q=2=>Goto 2↙LbI 1:“K”?S:“Z”?Z:“ANG”?M:Prog“SUB1”:“X=”:N◢“Y=”:E◢“FWJ=”:F▶DMS◢Goto 1↙LbI 2:“X=”?B:“Y=”?C:Prog“SUB2”:“K=”:S◢“Z=”:Z◢Goto 2↙SJK(数据库)Goto 1↙LbI 1↙If S＜线元终点桩号:Then线元起点方位角→A:线元起点桩号→O:线元起点X坐标→U: 线元起点Y坐标→V:线元起点半径→P: 线元终点半径→R:线元长度→L:Returu:IfEnd↙注意:更换曲线时，需加时LbI 1。
SUB1(正算子程序)0.5（1÷R-1÷P）÷L→D:S-O→X↙U+∫(cos(A+(X÷P+DX2)×180÷π,O,X)→N↙V+∫(sin(A+(X÷P+DX2)×180÷π,0,X)→E↙A+(X÷P+DX2)×180÷π→F↙N+Zcos(F+M)→N:E+Zsin(F+M)→E↙SUB2(反算子程序)Lbl 1:0→Z:1→Q:Prog“S JK”:Prog“SUB1”↙Pol(N-B+10^(-46),E-C+10^(-46)):Isin(F-90-J)→W:S+W→S↙Abs(W)>0.0001=>Goto1↙Lbl 2: 0→Z:Prog“SUB1”:(C-E)÷sin(F+90)→Z↙CASIO5800高程计算GCJS（主程序）LbI 1:Fix 3:?K:Prog“SJK”:R Abs(D-C)÷2→T:If D-C＞0:Then 1→W:Else -1→W: IfEnd↙If K＜A:Then -C→I:Else D→I:IfEnd↙Abs(A-K)→L↙If L-T＞0:Then Z+IL→H:“HS=”:H◢Else Z+IL+W(T-L)2÷2÷R→H:“HS=”:H◢IfEnd↙Goto 1↙SJK(数据库)曲线段If K＜终点里程:Then前坡度→C: 后坡度→D:半径→R:竖曲线交点桩号前→A:竖曲线交点高程→Z:Goto 1: IfEnd↙直线段If K＜终点里程:Then坡度→C: 坡度→D:1×10^(45)→R:起点桩号→A:起点高程→Z:Goto 1: IfEnd↙。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
CASIO5800坐标计算程序
转载自陈福龙转载于2010年05月13日 20:34 阅读(50) 评论(0) 分类：原创文学
举报
一、前言
本程序是《CASIO fx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。
原道路坐标放样计算程序只基于道路的单个基本型曲线，有效计算范围仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。
改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是：
1．使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦；
2．坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）；
3．坐标反算方面，输入平面坐标，即可计算对应的桩号和距中距离（含左右信息）；
4．对于存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现（可参照立交匝道程序中匝道的选择）。
程序的特点：
1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写；
2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐；
3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。
二、道路示例项目基本资料
基本资料同《CASIO fx-5800P计算与道路坐标放样计算》第6章HY高速公路第2合同段（合同段起止桩号：
K4+800~K9+600）。
这里摘取直线、曲线及转角表资料如下（若图片不清晰，请参见参见教材P161附录1）：
.
.
.
三、程序代码
.
.
.
.
.
.
注：路线数据库子程序ROAD－DATA1是根据计算示例项目的直曲表编写，大家使用时应按各自项目的直曲表改写或新建。
.
.
四、程序变量清单
.
.
五、计算流程示例
1．中桩坐标计算示例
计算任务：计算HY高速公路K6+100～K6+700段的中桩坐标及切线方位角（桩距20m），并在导线点（2807118.026，474113.687）上架设全站仪，计算各中桩的极坐标放样数据。
使用道路中边桩坐标放样正反算程序RAOD－2的操作流程见下表。
.
.
下表给出了由道路路线CAD软件计算的相关结果，供大家进行计算验证。
.
.
2．边桩坐标计算示例
计算任务：计算HY高速公路主线K6+100～K6+700段的边桩（左右各12.25米）坐标（桩距20m）。
使用道路中边桩坐标放样正反算程序RAOD－2的操作流程见下表。
.
.
下表给出了由道路路线CAD软件计算的相关结果，供大家进行计算验证。
.
.
3．坐标反算计算示例
计算任务：根据前面计算的K6+100的中桩、左右边桩的坐标，反算对应桩号及偏距，并进行验证。
使用道路中边桩坐标放样正反算程序RAOD－2的操作流程见下表。
.
.。