曲线任意里程中边桩坐标正反1doc4850.docVIP

曲线任意里程中边桩坐标正反算 CASIO?fx-4850P计算器 程序 作者：未知????文章来源：lq52搜集????点击数： 319????更新时间：2005-12-22 一、程序功能及原理1.功能说明：本程序由一个主程序 TYQXjs 和两个子程——正算子程序 SUB1 、反算子程序 SUB2 序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。2．计算原理：利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。二、源程序1.主程序 TYQXjs Defm 2:1.SZ XY：2.XY SZ：N：Prog”DAT1”:UX0：VY0：OS0：GF0：HLS：PR0：RRN：Q：C 1÷P：D P-R ÷ 2HPR ：E 180÷π：N 1 Goto 1：≠ Goto 2Δ←Lbl 1： SZT ：SZ:T”ANG”：W Abs S-O ：Prog SUB1：XS :X XYS :Y Y◢ “FS ”：F F◢Goto 3←┘Lbl 2： XY ：XY：I X：J Y：Prog SUB2：S :S O+WZ :Z Z◢Goto 2Lbl 3:J 0:I Pol X-Z[1],Y-Z[2] :F J:F 0 F F+360:≠ F 360 F F-36ΔΔ??”DIST ”:??I◢”FW ”:F→DMS◢ Goto 12. 正算子程序 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 C+KWD +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+ZsinF3. 反算子程序 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 Goto1：≠ W W+Z：Goto 0Δ←┘Lbl 1：Z 0：Prog SUB1：Z J-Y ÷sinF4.曲线元要素数据库：DAT1J 1 U xxxx.xxx:V xxxx.xxx:O xxxx.xxx:G xxxx.xxx:H xxxx.xxx:P xxxx.xxx:R xxxx.xxx:Q xxΔJ 2 U xxxx.xxx:V xxxx.xxx:O xxxx.xxx:G xxxx.xxx:H xxxx.xxx:P xxxx.xxx:R xxxx.xxx:Q xxΔ……….J n U………. Δ 注：如有多个曲线元要素继续添加入数据库DAT1中 三、使用说明1、规定 1 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q -1；当线元往右偏时，Q 1；当线元为直线时，Q 0。 2 当所求点位于中线时，Z 0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。 3 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。 4 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 5 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 6 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。 7 曲线元要素数据库（DAT1）可根据线型不同分为各个线元段输入到DAT1中，即分为直线段、缓和曲线、圆曲线等。2、输入与显示说明 输入部