《坐标正算主程序(命名为zbzs)》由会员分享,可在线阅读,更多相关《坐标正算主程序(命名为zbzs)(2页珍藏版)》请在金锄头文库上搜索。
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
2、-E):S=Icos(F-J):W=W+S 第 5 行:Abs(S)0.0001=Goto 0 第 6 行:T=Isin(J-F) 第 7 行:K=W 第 8 行:T”LEN”=TX1取样点的 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+
3、XI)90),0,S,4) 第 4 行:F=F+S(2P+SI)90 4. 道路线元参数道路线元参数(命名为命名为 A) 第 1 行:W175.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: WA+L =Goto 1 第 4 行:A=318.522: D=428482.988:E=5580
4、92.538:F=1181525.2:P=1/240:Q=0:L=55.104: WA+L =Goto 1 第 5 行:A=373.627:D=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曲线
5、段长度(尽量使用长度,为计算断链方便) 说明:(说明:(1)正算主程序可以计算一般边桩的坐标,如要计算类似涵洞端墙的坐标需增加两个变)正算主程序可以计算一般边桩的坐标,如要计算类似涵洞端墙的坐标需增加两个变量,具体方法参考本程序集中的第量,具体方法参考本程序集中的第 1 篇辛普生公式的坐标计算通用程序篇辛普生公式的坐标计算通用程序(2)程序规定,左偏曲线曲率(半径倒数)输入负值,右偏曲线曲率输入正值,直线上点曲率)程序规定,左偏曲线曲率(半径倒数)输入负值,右偏曲线曲率输入正值,直线上点曲率输入输入 0。(3)本程序精度极高,不受曲线半径大小影响,即使极小半径的螺旋曲线等误差仅为万分之一)本程
6、序精度极高,不受曲线半径大小影响,即使极小半径的螺旋曲线等误差仅为万分之一(0.1mm),可以忽略。),可以忽略。(4)若是从大里程向小里程的反方向计算,则曲率取正方向时的负值,方位角减去)若是从大里程向小里程的反方向计算,则曲率取正方向时的负值,方位角减去(或加上或加上)180 度。度。(5)有多个匝道的项目,可随时更改正反算主程序中的红色字体部分来调用其它线路的数据)有多个匝道的项目,可随时更改正反算主程序中的红色字体部分来调用其它线路的数据(6)反算桩号偏差为)反算桩号偏差为 1mm(7)可以计算任意线型的任意点坐标)可以计算任意线型的任意点坐标5. 坐标计算通用程序坐标计算通用程序(命
7、名为命名为 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 上测试成功,可以输出正确结果。