任意里程中边桩坐标正反算专用程序.docVIP

曲线任意里程中边桩坐标正反算 CASIO-4850计算器 专用程序 一、程序功能及原理1.功能说明：本程序由一个主程序 TYQXjs 和四个子程——正算子程序 SUB1 、反算子程序 SUB2 等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。本程序也可以在CASIO fx-4800P计算器运行，但由于容量较小，不建议使用。2．计算原理：利用Gauss-Legendre 5点通用 可查阅 2006年第二期，由武汉大学花向红所著---道路平面测设计算的通用公式 公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb 0 该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。二、源程序1.主程序 TYQXjs Lbl 0：Defm 400:1.SZ XY,2.XY SZ：N： N 1 Goto 1：≠ Goto 2Δ←Lbl 1： SZ ：SZS≤Z[26] J 1：Prog “DAT1”：Goto 3ΔS≤Z[34] J 2：Prog “DAT1”：Goto 3ΔS≤Z[42] J 3：Prog “DAT1”：Goto 3ΔS≤Z[50] J 4：Prog “DAT1”：Goto 3Δ ……………S≤Z[8（N+1）+10] J n：Prog “DAT1”：Goto 3ΔLbl 3：W Abs S-O ：Prog SUB1：XS :X X◢YS :Y Y◢”FWI ”: F F-M :F→DMS◢Goto 4←┘Lbl 2： JXY ：JXY：Z[4] X：Z[5] Y：N 0Lbl??A：Isz N：N J: Prog DAT1: Z[11] U：Z[12] V: Z[9] G:A Z[9]-M:J J+1: Prog DAT1: Z[13] U: Z[14] V: Z[10] G:B Z[10]-M:Prog ZX1 ：Z[6]×Z[7] 0 J N：Prog DAT1 ：Goto BΔA Z[9]+M：B Z[10]+M：Prog ZX1 ：Z[6]×Z[7] 0 J N：Prog DAT1 ：Goto B：≠ Goto AΔLbl??B：Prog SUB2：S :S O+W◢Z :Z Z◢Goto 2Lbl 4:J 0:I Pol X-Z[1],Y-Z[2] :F J:F 0 F F+360Δ??”DIST ”:??I◢”FW ”:F→DMS◢ Goto 12. 正算子程序 SUB1 A 0.11座机电话号码:B 0.座机电话号码52:N 0.座机电话号码44:K 0.0座机电话号码0:L 0.座机电话号码49:Z[3] 0.5:X U+W Acos G+QEKW C+KWD +Bcos G+QELW C+LWD +Ncos G+QEZ[3]W C+Z[3]WD +Bcos G+QE 1-L W C+ 1-L WD +Acos G+QE 1-K W C+ 1-K WD : Y V+W Asin G+QEKW C+KWD +Bsin G+QELW C+LWD +Nsin G+QEZ[3]W C+Z[3]WD +Bsin G+QE 1-L W C+ 1-L WD +Asin G+QE 1-K W C+ 1-K WD ：M”ANG ”: F G+QEW C+WD +M：X X+ZcosF：Y Y+ZsinF3. 反算子程序 SUB2 M”ANG ”:T G-M：W Abs Y-V cosT- X-U sinT ：Z 0：Lbl 0：Prog SUB1：L T+QEW C+WD ：Z Z[5]-Y cosL- Z[4]-X sinL：AbsZ 1E-6 Goto1：≠ W W+Z：Goto 0Δ←┘Lbl 1：Z 0：Prog SUB1：Z Z[5]-Y ÷sinF4. 垂距计算子程序 ZX1 Z[6] Z[5]-Z[12] COS A- Z[4]-Z[11] SIN AZ[7] Z[5]-Z[14] COS B- Z[4]-Z[13] SIN B5.曲线元要素数据库：DAT1Lbl1:J 1 U Z[16]:V Z[17]:O Z[18]:G Z[19]:H Z[20]:P Z[21]:R Z[22]:Q Z[23]ΔJ 2 U Z[24]:V Z[25]:O Z[26]:G Z[27]:H Z[