卡西欧4800坐标正反算通用程序(终极篇)

作者：柠檬树     QQ：249910569 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 0 K——计算点的里程 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.0001=>Goto0△ 第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)×90÷π 4.数据库（命名为A） 第1行：W≤175.191=> GoTo2△ 第2行：A=175.191:D=428513.730:E=557954.037:F=92°26′40″:P=0:Q=1/240:L=70.417: W≤A+L=>GoTo 1△ 第3行：A=245.607: D=428507.298:E=558024.092:F=100°50′59.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=118°15′25.2″:P=1/240:Q=0:L= 55.104: W≤A+L=>Goto 1△ 第5行：A=373.627:D=428453.283:E=558138.912:F=124°50′4.5″0:P=0:Q=-1/180:L=67.222: W≤A+L=>Goto 1△Goto 2 第6行：Lbl 1:Prog “XYF”: Goto 3 第7行：Lbl 2:D=0:E=0:F=0 第8行：Lbl 3 A——曲线段起点的里程 D——曲线段起点的x坐标 E——曲线段起点的y坐标 F——曲线段起点的坐标方位角 P——曲线段起点的曲率(左负右正) Q——曲线段终点的曲率(左负右正) L——曲线段长度（尽量使用长度，为计算断链方便） 说明： （1）正算主程序可以计算一般边桩的坐标，如要计算类似涵洞端墙的坐标需增加两个变量，具体方法参考本程序集中的第1篇辛普生公式的坐标计算通用程序 （2）程序规定，左偏曲线曲率（半径倒数）输入负值，右偏曲线曲率输入正值，直线上点曲率输入0。 （3）本程序精度极高，不受曲线半径大小影响，即使极小半径的螺旋曲线等误差仅为万分之一（0.1mm），可以忽略。 （4）若是从大里程向小里程的反方向计算，则曲率取正方向时的负值，方位角减去(或加上)180度。 （5）有多个匝道的项目，可随时更改正反算主程序中的红色字体部分来调用其它线路的数据 （6）反算桩号偏差为1mm （7）可以计算任意线型的任意点坐标 5.坐标计算通用程序(命名为ZB) 本程序与上述中的几个程序无关，是在计算坐标时采用输入曲线元参数模式下的程序，为临时使用方便。 第1行：A”K0“:D”X0“:E”Y0“:F”F0“:P”1÷R1“:Q”1÷R2“: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