曲线任意里程中边桩正反算5800P(修改版).docVIP

线路中边桩任意里程坐标正反算程序5800P ※Main???（主程序Lbl 0：Cls ＂1.SZ=2.NE=SZ＂:?Q Q=1 = Goto1:Q=2 = Goto2↙ Lbl 1:Cls:＂＂?S= Goto 0:Prog＂＂Lbl 1:＂＂?Z:Prog＂＂Fix 4:Cls↙ (注：小数点后保留4位有效数字) ＂X=＂:Locate 3,1,N＂Y=＂:Locate 3,2,E (注：输出正算Y值) ＂＂ ?DMS◢ (注：?DMS为一个符号,以度分秒格式显示，正算切线方位角) Prog＂＂Goto 1↙Lbl 2: Cls:＂＂? B＂= Goto 0:＂＂? C: ＂?S: B→N: C→E:Prog＂＂: ＂=＂:S◢＂=＂:ZGoto 2↙ （注：B=0时，程序会由Lbl 0开始重新运算，输入反算点坐标Y值，其中英文Predict/pri’dikt/中文含义“预计”，Predict DK=预计反算里程（呵呵，班门弄斧一下），作用：以便从程序库DATA提取相应的线元要素，否则计算器不停的计算，造成假死机） 【主程序完，主程序经此改动后,即可加入测站点坐标进行距离与方位角的放样了※FYJS （注：方位角和距离放样子程序名,与主程序Prog＂＂Goto 1对应）ClrStat↙XJ=?G:YJ=?H↙ （注：XJ,输入测站点X值，YJ，输入测站点Y值） Pol(N-G,E-H):Cls↙If J0:ThenJ+360→Y:Else J→Y:IfEnd↙＂FY =＂:,I＂FY =＂: Y?DMS?DMS度分秒格式，正算值至测站点方位角) 【方位角和距离放样子程序完】 ※DATA????????? (数据库子程序Goto1↙?????????????? 同时保存多个曲线时的指针Lbl 1↙ 第一条曲线（如：左线数据库）IF S***（线元终点里程）:Then***→A（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→L（线元起点至终点长度）: Return:IfEndbl 2↙ 第二条曲线（如：右线数据库） IF S***（线元终点里程）:Then***→A（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→L（线元起点至终点长度）: Return:IfEndIF S***:Then***→A:***→O:***→U:***→V:***→P:***→R: ***→L: Return:IfEnd↙ 【为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件。※QXJSSUB1???????????? (正算子程序0.5（1÷R-1÷P）÷L→D:S-O→XU+∫(cos(A+(X÷P+DX2)×180÷π,0,X)→N↙V+∫(sin(A+(X÷P+DX2)×180÷π),0,X)→E↙A+(X÷P+DX2)×180÷π→F↙N+Zcos(F+90) →N:E+Zsin(F+90) →E↙ 【正算子程序完，其中字母＂O＂用大写带下划线表示，数字＂0＂为小写不带下划线。】※QXJSSUB2??????????? (反算子程序Lbl 1:0→Z：1→Q：Prog＂＂Prog＂＂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＂＂:(C-E) ÷sin(F+90) →Z使用说明1、规定????????(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当曲线半径在左时，P、R取负值，当曲线半径在右时，P、R取正值，当曲线半径为无穷大（即直线）时，P、R以10的45次代替。????????(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。????????(3) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。???????? (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。???????? (6) 当线元为非完整缓和曲线时