Casio 5800P 主线坐标 高程 匝道坐标程序 2011-3-24修正.doc

Casio 5800P 计算程序 对程序进行部分更新，匝道计算比以前的公式精度高，程序公式简洁，缺点是计算速度稍慢。 2010/5/24 加入极坐标放样程序，对边桩坐标计算进行了改进，不在同时计算左右边桩 2010/6/30 匝道坐标计算程序 A03ZDZB 主程序名 “N=”?O:If O=0:Then Goto 0:Else “E=”?Q:IfEnd Lbl 0:”A=1,B=2,C=3,D=4,E=5,H=6,J=7”?A:”K=”?K If A=1 Then prog “A”:Goto 1:Ifend A匝道数据库 If A=2 Then prog “B”:Goto 1:Ifend B匝道数据库 If A=3 Then prog “C”:Goto 1:Ifend C匝道数据库 If A=4 Then prog “D”:Goto 1:Ifend D匝道数据库 If A=5 Then prog “E”:Goto 1:Ifend E匝道数据库 If A=6 Then prog “H”:Goto 1:Ifend H匝道数据库 If A=7 Then prog “J”:Goto 1:Ifend J匝道数据库 Lbl 1:prog “ZDJS”:Goto 0 主程序完。ZDJS为计算过程子程序 ZDJS 计算子程序 30→DimZ:K-F→I:(D-C)/(Z-F)*I→L:V+(2C)*I*90/π→P:If P0:Then 360+P→P:Else If P360:Then P-360→P:IfEnd:IfEnd:“HR=”:PDMS Norm 1: “HD=”?M:If M=0: Then Goto 7:Else “HJ=(Z=1,Y=+)”?G:IfEnd Lbl 7: Fix 3:“X=”:N+I*(cos(V)+4∑(cos(V+7.5(2H+1)*I*((2H+1)*L/12+2C)/ π),H,0,5)+2∑(cos(V+15HI(HL/6+2C)/ π),H,1,5)+cos(P))/36+Mcos (P+G)→X “Y=”:E+I*(sin(V)+4∑(sin(V+7.5(2H+1)*I*((2H+1)*L/12+2C)/ π),H,0,5)+2∑(sin(V+15HI(HL/6+2C)/ π),H,1,5)+sin(P))/36+Msin (P+G)→Y X-Z[29]:Y-Z[30]:If O=0:Then Goto 8:Else Prog “FS”:IfEnd Lbl 8:Norm 1: “10x”?S:K+10S→K 程序结束 FS 极坐标计算子程序 30-DimZ:Pol(Z[29]-O,Z[30]-Q):If J0:Then J+360-J:IfEnd: “H=”:JDMS Fix 3:“HD=”:I 程序结束 匝道数据库输入范例： A 直线、圆曲线、缓和曲线（包括卵形曲线）一种线形为一组要素 If K≤167.43:Then 0:→F:167.43→Z:4218460.504→N: 415641.464→E:39°48′47.5″→V:-1/62.75→C:0→D:IfEnd If K＞167.43: Then 167.43→F:334.684→Z:4218611.934→N: 415610.526→E:323°22′28.4″→V:0→C:1/300→D:IfEnd If K＞334.684: Then .........其他匝道数据输入同此 程序说明： O=测站点X坐标，输0为不计算测站点到待求点方位角与距离 Q=测站点Y坐标 A= 计算A匝道输入1，B输入2，C输入3......... K=待求点里程 HR=待求点的中线前进方位角 HD=中线到边桩的距离，输入正值即可，计算中桩坐标输入0 HJ=中线与边桩夹角，左边-值 右边+值 X=待求点X坐标 Y=待求点Y坐标 (H)J=测站点至待求点方位角 (HD)I=测站点至待求点距离 10x=桩号自动递增或递减，如果是计算20m递增的里程，可以输入2，计算循环后，K自动递增，输入-2自动递减。 数据库字母含义： F=线元起点里程 Z=线元终点里程 N=起点坐标X E=起点坐标Y V=起算点方位角 C=起点曲线半径的倒数，直线为0，曲线左转为负，右转为正 D=终点曲线半径的倒数，直线为0，曲线左转为负，右转为正 主线坐标程序 A01ZXZB主程序名 30-DimZ :Norm 1: “N=”?O:If O=0:Then Goto 9:Else “E=”?Q:IfEnd Lbl 9: “LI CHENG=”?K:prog “B01Z