曲线任意里程中边桩坐标正反算程序.docVIP

曲线任意里程中边桩坐标正反算程序 程序功能 本程序由一个主程序（TYQXjs）和两个子程序——正算子程序（SUB1）、反算子程序 SUB2 构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及CASIO fx-4850P计算器上运行。 源程序 主程序（TYQXjs） “1.SZ XY” : “2.XY SZ”:N:U“X0”:V“Y0”:O“S0”:G“F0”：H“LS”：P“R0”：R“RN”：Q：C 1÷P:D P-R ÷ 2HPR :E 180÷π：N 1 Goto 1:≠ Goto 2△ Lbl 1: SZ :SZ:W Abs S-O :Prog“SUB1”:X”XS” X▲ Y”YS” Y▲ Goto 1 Lbl 2: XY :XY:I X:J Y:Prog”SUB2”:S”S” O+W▲ Z”Z” Z▲ Goto 2 2．正算子程序（SUB1） A 0.1座机电话号码6:B 0.座机电话号码74:K 0.0座机电话号码2：L 0.座机电话号码82:F 1-L:M 1-K: X U+W Acos G+QEKW C+KWD +Bcos G+QELW C+LWD +Bcos G+QEFW C+FWD +Acos G+QEMW C+MWD : Y V+W Asin（G+QEKW CKWD ）+Bsin G+QELW C+LWD +Bsin G+QEFW C+FWD + Asin G+QEMW C+MWD : F G+QEW C+WD +90:X X+ZcosF:Y Y+ZsinF 3．反算子程序 SUB2 T G-90:W Abs（（Y-V）cosT- X-U sinT）:Z 0:Lbl 0：Prog“SUB1”：L T+QEW C+WD : Z J-Y cosL- I-X sinL：AbsZ 1E-6 Goto 1: ≠ W W+Z:Goto 0△ Lbl 1:Z 0:Prog“SUB1”：Z J-Y ÷sinF 使用说明 规定 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q -1；当线元往右偏时，Q 1；当线元为直线时，Q 0。 当所求点位于中线时，Z 0；当位于中线左侧时，Z取负值；当位于中线右侧时，Z取正值。 当线元为直线时，其起点、终点的曲率半径为无穷大，以10的45次代替。 当线元为圆曲线时，无论其起点、终点与什么线元相接，其曲率半径均等于圆弧的半径。 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于全曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值，与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于全曲线的半径。 输入与显示说明 输入部分： 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 ? 反算时所求点的X坐标 Y ? 反算时所求点的Y坐标 显示部分 XS ×××× 正算时，计算得出的所求点的X坐标 YS ×××× 正算时，计算得出的所求点的Y坐标 S ×××× 反算时，计算得出的所求点的里程 Z ×××× 反算时，计算得出的所求点的边距 算例 某匝道是由五段线元（直线+完整缓和曲线+全曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如下： S0 X0 Y0 F0 LS R0 RN Q 500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0 769.256 19787