4850计算器坐标正反算交点法程序.doc

4850坐标正反算(交点法与线元法的联合应用)全线通 f(x)4850全线通坐标正、反算Deg:Defm 11:M0GOTO X:GOTO 0 -当M=0,为反算；当M0时，为正算。LbI 0: OX,OY,OS”：XYS: -中线外任意点X、Y坐标，S为一个曲线附近任意一个桩号。Z8=X:Z9=Y+0.0001:S”S=”:Z=1GOTO Y:Z=2GOTO ZLbI X: Z: Z=-1GNEAFRC: S:S:GOTO 1:Z=-2NEFPRQ: S:S:GOTO 6S:Z=1GOTO Y:Z=2GOTO ZLbI Z:S等缓曲线终点桩号G=等缓曲线交点桩号：N=交点（X）坐标：E=交点（Y）坐标：A=转角(左转为负,右转为正)：F=起点切线方位角：R=半径：C=缓和曲线长度：GOTO 1S线元终点桩号Q=线元起点桩号：=线元终点桩号；N=线元起点（X）坐标：E=线元起点（Y）坐标：F=线元起点切线方位角：P=线元起点曲率半径(左转为负,右转为正)：R=线元终点曲率半径(左转为负,右转为正)：GOTO 6S.LbI Y:S等缓曲线终点桩号G=等缓曲线交点桩号：N=交点（X）坐标：E=交点（Y）坐标：A=转角(左转为负,右转为正)：F=起点切线方位角：R=半径：C=缓和曲线长度：GOTO 1S线元终点桩号Q=线元起点桩号：=线元终点桩号；N=线元起点（X）坐标：E=线元起点（Y）坐标：F=线元起点切线方位角：P=线元起点曲率半径(左转为负,右转为正)：R=线元终点曲率半径(左转为负,右转为正)：GOTO 6S.LbI 1：P=C2/24R-C4/2688R3+C6/678450R5：Q=C/2-C3/240R2+C5/41472R4：B=90C/ R：T=（R+P）tan(AbsA/2)+Q:W=(R+P)/cos(A/2)-R:L=(AbsA-2B) R/180+2C:Z1=G-T:Z2=Z1+C:Z3=Z1+L/2:Z4=Z1+L-C:Z5=Z1+L:A0I=-1:I=1: U=F+A/2+90I:V=W+R:B=N+VcosU:=E+VsinU:LbI 2：SZ1L=T+Z 1-S:V=F+180:U=F:GOTO 3SZ2L=S-Z 1:Q=L-L5/40R2 C2+L9/3456(RC)4:P=L3/6RC-L7/336(RC)3+L11/42240(RC)5:V=(Q2+P2):L=30IL2/R C:Q=F+Itan-1(P/Q):P=F+180:U=F+3L:GOTO 4SZ4L=F+A/2+180+90I+180(S-Z3)I/R :U=L+90I:GOTO 5SZ5L= Z5-S:Q=L-L5/40R2 C2+L9/3456(RC)4:P=L3/6RC-L7/336(RC)3+L11/42240(RC)5:V=(Q2+P2):L=30IL2/R C:Q=F+A+180-I tan-1(P/Q):P=F+A:U=F+A-3L:GOTO 4SZ5L=S-Z5+T:V=F+A:U=V:GOTO 3LbI 3：X=N+LcosV:Y=E+LsinV:U=U:M=0GOTO 8: GOTO 7LbI 4:X=N+TcosP+VcosQ:Y=E+TsinP+VsinQ:U=U:M=0GOTO 8: GOTO 7LbI 5:X=B+RcosL:Y=+RsinL:U=U:M=0GOTO 8: GOTO 7LbI 6:Z6=(R-P)/Abs(-Q)Z7= Abs(S-Q):I= Z6* Z7:U=F+90(I+2P)* Z7/:A=F+45(I/4+2P) Z7/2:B=F+135(3I/4+2P)* Z7/2:C=F+45(I/2+2P) Z7/:X=N+ Z7/12*(cosF+4*(cosA+cosB)+2cosC+cosU):Y=E+ Z7/12*(sinF+4*(sinA+sinB)+2sinC+sinU):M=0GOTO 8: GOTO 7LbI 7:DV:G=X+Dcos(U+V):H=Y+Dsin(U+V):U360U=U-360S:Pause 0:G: Pause 0:HD:D:D0GOTO 7: GOTO 9LbI 8:J=0:Z10=(X-Z8)sin(U-90)-(Y-Z9)cos(U-90):Abs Z100.001S=S+ Z10:POI(Z8-X, Z9-Y):J=J-U:J0J=J+360: J=J J180Z11=1: Z11=-1“OK”: