1、1. 坐标正算主程序(命名为 ZBZS)第 1 行：Lbl 0:W,T,B:W”K=”:T”LEN=”:B第 2 行：Prog “A”第 3 行：X=D+Tcos(F+B)第 4 行：Y=E+Tsin(F+B)第 5 行：F=F 第 6 行：Goto 0K计算点的里程LEN计算点到中桩的距离（左负右正）B取前右夹角为正2. 坐标反算桩号和偏距主程序(命名为 ZBFS)第 1 行：U,V,K:U”X1”:V”Y1”: W”K1”第 2 行：I=0:J=0第 3 行：Lbl 0:Prog “A”第 4 行：Pol(U-D,V-E):S=Icos(F-J):W=W+S第 5 行：Abs(S)0.00
2、01=Goto 0第 6 行：T=Isin(J-F)第 7 行：K=W第 8 行：T”LEN”=T X1取样点的 X 坐标Y1取样点的 Y 坐标K1输入时为计算起始点 (在线路内即可)，输出时为反算点的桩号T偏距（左负右正）3. 计算坐标子程序(命名为 XYF)为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第 1 行：S=W-A：I=(Q-P)L第 2 行：D=D+(cos(F+X(2P+XI)90),0,S,4)第 3 行：E=E+(sin(F+X(2P+XI)90),0,S,4)第 4 行：F=F+S(2P+SI)904. 道路线元参数(命名为 A)第 1 行：W1
3、75.191= GoTo 2第 2 行：A=175.191:D=428513.730:E=557954.037:F=922640:P=0:Q=1/240:L=70.417: WA+L =GoTo 1第 3 行：A=245.607: D=428507.298:E=558024.092:F=1005059.4:P=1/240:Q=1/240:L=72.915: W A+L =Goto 1第 4 行：A=318.522: D=428482.988:E=558092.538:F=1181525.2:P=1/240:Q=0:L=55.104: WA+L =Goto 1第 5 行：A=373.627:D=
4、428453.283:E=558138.912:F=124504.50:P=0:Q=-1/180:L=67.222: WA+L=Goto 1Goto 2第 6 行：Lbl 1:Prog “XYF”: Goto 3第 7 行：Lbl 2:D=0:E=0:F=0第 8 行：Lbl 3A曲线段起点的里程D曲线段起点的 x 坐标E曲线段起点的 y 坐标F曲线段起点的坐标方位角P曲线段起点的曲率(左负右正)Q曲线段终点的曲率(左负右正)L曲线段长度（尽量使用长度，为计算断链方便）说明：（1）正算主程序可以计算一般边桩的坐标，如要计算类似涵洞端墙的坐标需增加两个变量，具体方法参考本程序集中的第 1 篇辛普
5、生公式的坐标计算通用程序（2）程序规定，左偏曲线曲率（半径倒数）输入负值，右偏曲线曲率输入正值，直线上点曲率输入 0。（3）本程序精度极高，不受曲线半径大小影响，即使极小半径的螺旋曲线等误差仅为万分之一（0.1mm），可以忽略。（4）若是从大里程向小里程的反方向计算，则曲率取正方向时的负值，方位角减去(或加上)180 度。（5）有多个匝道的项目，可随时更改正反算主程序中的红色字体部分来调用其它线路的数据（6）反算桩号偏差为 1mm（7）可以计算任意线型的任意点坐标5. 坐标计算通用程序(命名为 ZB)本程序与上述中的几个程序无关，是在计算坐标时采用输入曲线元参数模式下的程序，为临时使用方便。第 1 行：A”K0“:D”X0“:E”Y0“:F”F0“:P”1R1“:Q”1R2“:L“L0“第 2 行：Lbl 0:K,T,B:K”K”:T”LEN”:B第 3 行：S=Abs(K-A)：I=(Q-P)L第 4 行：D=D+(cos(F+X(2P+XI)90),0,S,4)第 5 行：E=E+(sin(F+X(2P+XI)90),0,S,4)第 6 行：F=F+S(2P+SI)90第 7 行：X=D+Tcos(Z+B)第 8 行：Y=E+Tsin(Z+B)第 9 行：F=F （可省掉）第 10 行：Goto 0本程序已在 4800/4850 上测试成功，可以输出正确结果。
