线路中边桩任意里程坐标正反算改.pdfVIP

线路中边桩任意里程坐标正反算改 含数据库、测站放样计算（第五次修改-正式版） 1.加编数据库及计算总调度程序,计算中不必逐项输入线元要素,提高 运算速度，避免现场忙中出错 2.可加入多条线路的数据库,内业输入数据,外业一目了然 3.计算直观,人性化 4.正算直接输入里程和边距,反算输入近似里程便可 5.增加了“计算点与测站点”的距离和方位角计算语句，方便直接放样 6.愿收获与大家共享,同时也希望大家提出心得和宝贵建议 7.核心计算程序摘自“yshf” 一.改动后的程序清单，增设数据库程序,可输入 n 条线路的数据库 1.(QXZDJS 计算总调度程序) M=1=Prog SJK1:Prog SUBSJK△←┘ M=n=Prog SJKn:Prog SUBSJK△←┘ 2.(SJK1 数据库程序) 1.SZ = XY ：2.XY = SZ ：｛N,S ｝：N ：SDKI←┘ S≤本线元终点里程＝＞U=本线元起点 X 坐标：V=本线元起点 Y 坐标： G=本线元起点正切线方位角：P=本线元起点曲率半径：R=本线元终点 曲率半径：O=本线元起点里程：H=本线元终点里程：Q=本线元左直右偏 向(0 或 1、-1):Goto0Δ←┘ S≤下一线元终点里程= ＞O=。。。。。。。。。。。。。。。。： Goto0Δ←┘ LB1 0 3.(SUBSJK 运算主程序) {UVGPROHQ}：UXA ：VYA ：GFA ：PRA ：RRB ：ODKA ： HDKB ：Q←┘ Deg:SH=Goto AΔD=(P-R)÷2PR(H-O) ：N=1=Goto 1 ：≠Goto 2Δ←┘ Lbl 1 ：{Z}：Z ：W=Abs(S-O) ：Prog SUB1 ：XS= ：X◢ YS= ：Y◢ QX-FWJ= ：F=F-90←DMS◢ I=0 ：J=0 ：Pol （X-CXC ，Y-EYC）：J0 ＝＞J=J+360 ： Δ“FWJ=”:J←DMS◢ “I=” ：I◢ Goto A←┘ Lbl 2 ：{XY}：I=X ：J=Y ：Prog SUB2 ：S=:S=O+W◢ Z=:Z◢ Lbl A 4. (SUB1 正算子程序) A=0.1739274226 ：B=0.3260725774 ：K=0.0694318442 ：L=0.3300094782 ： X=U+W(Acos(G+180QKW(1÷P+KWD)÷π)+Bcos(G+180QLW(1÷P+L WD)÷π)+Bcos(G+180(1-L)QW(1÷P+(1-L)WD)÷π)+Acos(G+180(1-K)Q W(1÷P+(1-K)WD)÷π)) ： Y=V+W(Asin(G+180QKW(1÷P+KWD)÷π)+Bsin(G+180QLW(1÷P+LW D)÷π)+Bsin(G+180(1-L)QW(1÷P+(1-L)WD)÷π)+Asin(G+180(1-K)QW( 1÷P+(1-K)WD)÷π)) ：F=G+180QW(1÷P+WD)÷π+90 ：X=X+ZcosF ： Y=Y+ZsinF 5.(SUB2 反算子程序) T=G-90 ：W=Abs((Y-V)cosT-(X-U)sinT) ：Z=0 ：Lbl 0 ：Prog SUB1 ： L=T+180QW(1÷P+WD)÷π ：Z=(J-Y)cosL-(I-X)sinL ： AbsZ1÷10^6=Goto1 ：≠W=W+Z ：Goto 0Δ←┘ Lbl 1 ：Z=0 ：Prog SUB1 ：Z=(J-Y)÷sinF 二. 使用说明： 1.把所有相关的”线元要素“依次输入”SJK“，如果有多条线路，可把 “SJK”分成 SJK1...SJKn 等 n 条线路数据库，然后对 “计算总调度程序” 中的 M1...Mn 备用 2.运算时直接调用“SJK”运行，程序提示输入里程“S” ？时，正算直接 输入待求点里程，反算输入所求点“近似“里程 3.如果在正算中途想进行反算，重新给”N“ ？赋值2便可 4.程序中“XC、YC”为测站坐标，“I ＝”、“FWJ ＝”为放样距离和方位角 5.程序中“QX-FWJ=”为计算点在中线上的正切线方位角，此方位角在 新线路输入数据库时非常重要，我们在每输入完一段曲线元要素后， 立即退出返回到“计算总调度程序(SUBSJK)”中，然后输入该段曲线元 的终点里程，其计算出来的 X，Y 中桩坐标应该与设计图纸(当然是经 过复查无误的图纸