公路路线座标正反算(500计算器)添加隧道超欠挖.doc

公路路线座标正反算（积分公式）通用程序 公路路线座标正反算 由于现在计算机普及，计算机功能日益强大，宜采用较简单的积分公式，便于计算机处理。 单线元通用积分公式如下 M = (1.0/Re-1.0/Rs)/Ls; x=∫｛cos(Ta + L/Rs + 0.5*M *L*L)，0，L｝; y=∫｛sin(Ta +L/Rs + 0.5*M *L*L)，0，L｝; a(i)= Ta +L/Rs + 0.5*M *L*L Rs:缓和曲线起点半径 Re:缓和曲线止点半径 Rs，Re （NE坐标系下，右偏为正，左偏为负） Ta:缓和曲线起点的真北方位角 Ls:不完整缓和曲线长度。 此公式为缓和曲线在坐标系下任意位置的通用积分公式，能完全适应缓和曲线左偏、右偏、Rs Re 、Rs Re等各种情况，不必先凑成完整缓和曲线，降低算法的复杂程度。虽然此公式是由缓和曲线推导出来，也可和于直线与圆曲线，可降低计算机编程的复杂程度。 Fx-5800计算机程序 QXJS-000 主程序 Lbl 4:“1.SZ=NE”:“2.NE=SZ”:?Q:Z[12]→S：?S：Prog“QXJS-SUB0”↙ Lbl 0:Q=1 = Goto1:Q=2 = Goto2:↙ Lbl 1:?Z:Prog“QXJS-SUB1”:“N=”:N◢:“E=”:E◢:“F=”:F DMS◢: Goto3↙ Lbl 2: “N=”:?B: “E=”:?C:B→N: C→E:Prog“QXJS-SUB2”: “S=”:S◢: “Z=”:Z◢: Goto4↙ Lbl 3： J＜0＝＞J＝J+360： “FW＝”：J DMS◢ Goto4↙ S→Z[12]:Z→Z[13] Prong “SD” QXJS-SUB0 数据库子程序 IF S***（线元终点里程）:Then***→A（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R（线元终点曲率半径）: ***→L（线元起点至终点长度）: Return:IfEnd↙ IF S35693.064:Then 257°13′1″→A:35643.064→O: 3140794.382→U: 514646.536→V:10^(45)→P:965→R:50→L: Return:IfEnd↙ QXJS-SUB1 正算子程序 0.5（1÷R-1÷P）÷L→D:S-O→X↙ U+∫(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 QXJS-SUB2 反算子程序 Lbl 1:0→Z：1→Q：Prog“QXJS-SUB0”: Prog“QXJS-SUB1”↙ 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“QXJS-SUB1”:(C-E) ÷sin(F+90) →Z 隧道超欠挖程序 程序名：SD Prog“SD1”: 0.84→A:”CQHD”?A:A→Z[15]:0→M:”Hc1=”?M:M→Z[8]:Prog“SQX”:Cos(Z[10]) ×Z[4] →E: H+Z[5] →Z[5]: H+Z[7] →Z[7]: E+Z[5] →Z[10]:If Z[8]Z[10]:Then Goto R:Else Goto S:IfEnd:LbI R: √((Z[13]-Z[9]) 2+(Z[8]-Z[5]) 2 )- Z[4]-Z[15] →L:Goto L: LbI S: Z[4]-Z[6] →S: √(S 2-(Z[7]-Z[5]) 2 ) ×S÷Abs(S)→T:Abs(Z[13]-Z[9])-T→T:√(T 2+(Z[8]-Z[7]) 2 )- Z[6]-Z[15] →L:Goto L: LbI L: Fix 3: ＂=＂: Locate 6,4, ＂=＂:Z[8]-H ＂Ｋ＂]： ＂＂Prong “S1” 0.5R Abs(I-J)→T If I-J﹤0: Then -1→V:Else 1→V：IfEnd K-Z→B: If K≤Z-T:Then BI+A →H:Return:IfEnd: If K≤Z:Then A+BI-V(T+K-Z)2÷（2R）→H:Return:IfEnd If K≤Z+T:Then A+BJ-V(T-K+Z)2÷（2R