修正后5800坐标正反算分段曲线程序(1)

5800 程序( 曲线任意里程中边桩坐标正反算) 一、程序功能 本程序由一个主程序(QXZBJS)和两个子程 正算子程序 (SUB1)、反算子程序(SUB2)序构成，可以根据曲线段直线、 圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起 点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半 径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进 行正反算。另外也可以将本程序中核心算法部分的两个子程序移植 到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。 特别申明：(1). 适用于弧长小于 2 倍半径的各种线元坐标正反算， 精度优于 1mm； QXZBJS(主程序名) LbI 0:4Dimz:Cls ?S Prog“SUB3“ “1.SZ XY“ “2.XYSZ“ ? N:Cls 1PC：(P-R)(2HPR) D ：180 E N=2=Goto 2 LbI 1:Cls:?S S=0=Goto 0 (注: 当 S 输入 0 时,程序会由 LbI 0 开始重新运行) ?Z:Abs(S-O)W Prog“SUB1“:Cls “XS=“:X “YS=“:Y “FS=“:F-90 DMS Goto 0 LbI 2:Cls:?X X=0=Goto 0 (注: 当 X 输入 0 时,程序由会 LbI 0 开始重新运行) ?Y:X I:YJ Prog“SUB2“:O+WS:Cls “S=“:S “Z=“:Z+0.0000（此处为距线路中心偏距离左负右正） Goto 3 LbI 3 “A=”：(S-104822) 0.008+378.51A ?B If BA+6.445:Then“K=”: （Z2+(B-(A+4.43)2）-2.65K Else If BA+4.474:Then “K=”: (Abs(Z)+0.7)2+(B-(A+3.73)2)-3.64K Else If BA+1.79:Then “K=”: ( Abs(Z)+6.72)2+(B-(A+2.54)2)-9.78K Else If BA-1.92:Then “K=”: (Abs(Z)-3.06+Abs(B-A-1.79) tan(22415)K IfEnd: IfEnd: IfEnd: IfEnd: Goto 0 SUB1(子程序名 ) 0.1739274226A:0.3260725774B:0.0694318442K:0.330009482L 1-LF:1-KM G+QEKW(C+KWD)Z1 G+QELW(C+LWD)Z2 G+QEFW(C+FWD)Z3 G+QEMW(C+MWD)Z4 Acos(Z1)X X+Bcos(Z2)X X+Bcos(Z3)X X+Acos(Z4)X U+WXX Asin(Z1) Y Y+Bsin(Z2)Y Y+Bsin(Z3)Y Y+Asin(Z4)Y V+WYY G+QEW(C+WD)+90F X+Zcos(F)X Y+Zsin(F)Y 此处把 A B K 字母替换成数字！ 0.1739274226A:0.3260725774B:0.0694318442K:删除，便于隧道程序提出 字母！ SUB2(子程序名 ) G-90 T (Y-V)cos(T)-(X-U)sin(T)W Abs(W)W:0Z LbI 4:Prog“SUB1“ T+QEW(C+WD)L (J-Y)cos(L)-(I-X)sin(L)Z If Abs(Z)0.000001:Then 0Z: Prog“SUB1“ (J-Y)sin(F)Z Else W+ZW:Goto 4:IfEnd SUB3(子程序名 ) Goto 1 Lb1 1：If S96101.247：Then 578608.630U(线元起点的 X 坐标):492589.754 V( 线 元起点的 Y 坐标):95951.247O(线元起点里程):390901G （线元起点切线方位角） :150H（线元长度）:1045P（线元起点曲率半径）:1200R（线元终点曲率半径） :-1Q（线 元左右偏标志(左偏 Q=-1，右偏 Q=1，直线段 Q=0) ）:Return:Ifend Lb1 1：If S97063.338：Then 578726.881U:492681.998V:96101.247O:35 3409 G:962.091H:1200P:1200R:-1 Q:Return:Ifend 为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个 Lbl，每增加一个工程增加一个文件 注:1E-6 即是 10 的负 6 次方, 亦即是等于 0.000001 输入时按 SHIFT-log-(-)-6 亦可直接将 If Abs(Z)1E-6 写成 If Abs(Z)0.000001 1. SZ = XY 2. XY = SZ N ? 选择计算方式，输入 1 表示进行由里程、边距计算坐标 ；输入 2 表示由坐标反算里程和边距。 X0 ？线元起点的 X 坐标 Y0 ？线元起点的 Y 坐标 S0 ？线元起点里程 F0 ？线元起点切线方位角 LS ？线元长度 R0 ？线元起点曲率半径 RN ？线元止点曲率半径 Q ？ 线 元左右偏标志(左偏 Q=-1，右偏 Q=1，直线段 Q=0) S ？ 正算时所求点的里程 Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零) X ？