简洁一个交点计算坐标与反算桩号程序.docVIP

Deg：Defm4 Lb0:｛Q｝：Q≠0=｛L｝：L： GOTO1⊿｛VH｝：L=0：V“XA”：H“YA” Lb1：N=1 Lb2：Fixm：Prog“JD”:U=E：M=F N=N+2:Prog“JD”：P=E：T=F N=N-1:Prog“JD” Pol(P-E,T-F:G=J：Pol(E-U,F-M U=G-J:AbsU＞180=＞U=U-360AbsU÷U⊿ D=U÷AbsU：M=S÷2-S∧3÷240R2 P=S2÷24R-S∧4÷2688R∧3：@=A2÷24R-A4÷2688R3 第二缓曲内移量 Abs（S-A）＞0=〉T=（@-P）÷（2tg0.5U）+0.5(@＋P+2R)×tg0.5U+M Z[3]=（P－@）÷（2tg0.5U）＋0.5(@＋P＋2R)×tg0.5U+A÷2-A3÷240R2: ≠=T=(R+P)Abstg0.5U+M：Z[3]=T⊿ @=AbsU÷(180÷π)×R+(S+A)÷2 曲线全长 L＞Z+@-T=＞C=0:K=G:X=L－（Z+@-T）＋Z[3]：Y=Q：≠=＞ LZ-T+@-A=＞W=Z-T-L+@:C=Z[3]：S=A:K=G-180:Prog“LU”：Y=-Y： ≠=＞LZ-T+S=＞C=T：K=J：W=L-Z＋T: Y=D（P+R-Rec（R-DQ，（（W-0.5S）÷R×(180÷π：X=M+J： ≠=＞L＞Z-T=＞W=T＋L－Z:C=T：K=J:Prog“LU”： ≠=＞C=T:K=J:X= L－Z＋T ：Y=Q⊿⊿⊿⊿ Pol（X-C，Y：J=J＋K X=E+Rec（I，J：Y=J+F Pol（Z[1]-X，Z[2]-Y Q≠0=X“XK”◢Y“YK”◢I“JL” ◢ J“FWJ”=180+J◢GOTO 0⊿ Pol(V-X,H-Y L＞Z+@-T=＞K=K：≠=＞LZ-T+@-A=＞K=K-180-DZ[4]：≠=＞LZ-T+S=＞K=K＋180D（W-0.5S）÷πR：≠=＞L＞Z-T=＞K=K＋DZ[4]：≠=＞K=K⊿⊿⊿⊿ @=Icos（J-K Abs@≤0.005=＞L“ZH HAO”=L+@◢ @“B K”=Isin(J-K ◢ ≠=＞L=L+@：GOTO 1⊿ 1 此程序约780字节。如果要进行极坐标放样，应先把测站坐标X、Y输入计算器寄存器Z[1] 、Z[2]。 2 如要计算线路高程可在J“FWJ”=180+J◢和@“B K”=Isin(J-K ◢后面运行子程序Prog“SU” 3 当半宽Q输入数字“0”时，程序将根据输入的坐标XA、YA求算桩号L+@及 半宽@“B K”；当半宽Q输入数字非“0”时，则以输入的桩号L、Q计算线路坐标XK◢ YK◢与测站坐标的距离I“JL”◢及方位角 J“FWJ”=180+J◢ 4 U=G-J 转角左“-”，右“+”。 T 第一切线长度 Z[3] 第二切线长度 5 @ 区别阿拉伯数字0与英文字母O； Z[1] Z[2] 测站坐标X、Y； 6 V“XA”:H“YA”需要求算桩号及半宽的坐标；L 所求桩号； Q 所求半宽 7 D=U÷AbsU 运算符号及自行判断左右拐左“-1”右“+1”； T第一切线长度 8 L＞Z+@-T 缓直点桩号； LZ-T+@-A圆缓点桩号； LZ-T+S 9 Z-T直缓点桩号； K＋180D（W-0.5S）÷πR 对应曲线长度切线方位角 LU Z[4]=W2÷2RS×180÷π：X=W-W∧5÷40R2S2-W9÷3456R4S4-DQsin Z[4] Y=DW∧3÷6RS-DW∧7÷336R∧3S∧3+Qcos Z[4] JD N=1=＞E=23535.019：F=11333.96⊿ N=2=＞Z=121112.164：E=19114.638：F=13055.709：R=4500：S=320 A=320 N=3=＞E=12357.419：F=10728.347⊿ N=1只输入线路第一交点“前直线段”任意中心桩号点坐标E、F即可 N=2输入交点参数 N=3终点只输入线路终交点的“后直线段”任意中心桩号坐标E、F即可 Z 交点桩号； EF交点坐标； R 半径； S 第一缓和曲线长度LS1； A 第二缓和曲线长度LS2； 如果无缓和曲线则S=0、A=0；如果等缓和曲线则S=LS1、A=LS1、 SU LQZ：Fixm：W=1：U=0 Lb2: Prog“SD” C=I:D=J W=W+2:Prog“SD”E=I:F=J W=W-1: Prog“SD” C=(J-D)÷(I-C：D=(F-J)÷(E-I：T=0.5RAbs(D-C LI+T=GoTo2⊿ LI-T=U=(T-Ab