4850坐标正反算(交点法与线元法的联合应用)全线
通
f(x)—4850全线通坐标正、反算 Deg:Defm 11:
M≠0＝＞GOTO X:≠＞GOTO 0⊿ ---当M=0,为反算；当M≠0时，为正算。 LbI 0: \"OX,OY,OS”：｛XYS｝: ---中线外任意点X、Y坐标，S为一个曲线附近任意一个桩号。 Z〔8〕=X: Z〔9〕=Y+0.0001:
S”S=”:Z=1＝＞GOTO Y:≠＞Z=2＝＞GOTO Z⊿ LbI X: ｛Z｝:
Z=-1＝＞GNEAFRC: ｛S｝:S:GOTO 1:≠＞Z=-2＝＞NEFPRQ○: ｛S｝:S:GOTO 6｛S｝:≠＞Z=1＝＞GOTO Y:≠＞Z=2＝＞GOTO Z⊿⊿⊿⊿⊿ LbI Z:
S≤等缓曲线终点桩号＝＞G=等缓曲线交点桩号：N=交点（X）坐标：
E=交点（Y）坐标：A=转角(左转为负,右转为正)：F=起点切线方位角：R=半径：C=缓和曲线长度：GOTO 1⊿
S≤线元终点桩号＝＞Q=线元起点桩号：○=线元终点桩号；N=线元起点（X）坐标：E=线元起点（Y）坐标：F=线元起点切线方位角：P=线元起点曲率半径(左转为负,右转为正)：R=线元终点曲率半径(左转为负,右转为正)：GOTO 6⊿ S≤„„„„„„„„. LbI Y:
S≤等缓曲线终点桩号＝＞G=等缓曲线交点桩号：N=交点（X）坐标：
E=交点（Y）坐标：A=转角(左转为负,右转为正)：F=起点切线方位角：R=半径：C=缓和曲线长度：GOTO 1⊿
S≤线元终点桩号＝＞Q=线元起点桩号：○=线元终点桩号；N=线元起点（X）坐标：E=线元起点（Y）坐标：F=线元起点切线方位角：P=线元起点曲率半径(左转为负,右转为正)：R=线元终点曲率半径(左转为负,右转为正)：GOTO 6⊿ S≤„„„„„„„„. 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: Z〔1〕=G-T: Z〔2〕=Z〔1〕+C: Z〔3〕=Z〔1〕+L/2: Z〔4〕=Z〔1〕+L-C: Z〔5〕=Z〔1〕+L: A＜0＝＞I=-1:≠＞I=1: ⊿ U=F+A/2+90I: V=W+R: B=N+VcosU: ○=E+VsinU: LbI 2：
S≤Z〔1〕＝＞L=T+Z 〔1〕-S: V=F+180: U=F: GOTO 3⊿
S≤Z〔2〕＝＞L=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 4⊿
S≤Z〔4〕＝＞L=F+A/2+180+90I+180(S-Z〔3〕)I/πR :
U=L+90I: GOTO 5⊿
S≤Z〔5〕＝＞L= Z〔5〕-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 4⊿
S＞Z〔5〕＝＞L=S-Z〔5〕+T: V=F+A: U=V: GOTO 3⊿ LbI 3： X=N+LcosV: Y=E+LsinV: U=U:
M=0＝＞GOTO 8: ≠＞GOTO 7⊿ LbI 4:
X=N+TcosP+VcosQ: Y=E+TsinP+VsinQ: U=U:
M=0＝＞GOTO 8: ≠＞GOTO 7⊿ LbI 5: X=B+RcosL: Y=○+RsinL: U=U:
M=0＝＞GOTO 8: ≠＞GOTO 7⊿ LbI 6:
Z〔6〕=(R-P)/Abs(○-Q) Z〔7〕= Abs(S-Q): I= Z〔6〕* Z〔7〕: U=F+90(I+2P)* Z〔7〕/π: A=F+45(I/4+2P) Z〔7〕/2π: B=F+135(3I/4+2P)* Z〔7〕/2π: C=F+45(I/2+2P) Z〔7〕/π:
X=N+ Z〔7〕/12*(cosF+4*(cosA+cosB)+2cosC+cosU): Y=E+ Z〔7〕/12*(sinF+4*(sinA+sinB)+2sinC+sinU): M=0＝＞GOTO 8: ≠＞GOTO 7⊿ LbI 7:｛DV｝: G=X+Dcos(U+V): H=Y+Dsin(U+V): U<0＝＞U=U+360⊿ U>360＝＞U=U-360⊿ S:Pause 0:G: Pause 0:H▲
｛D｝:D:D≠0＝＞GOTO 7: ≠＞GOTO 9⊿ LbI 8:J=0:
Z〔10〕=(X-Z〔8〕)sin(U-90)-(Y-Z〔9〕)cos(U-90):
Abs Z〔10〕≤0.001＝＞S=S+ Z〔10〕:POI(Z〔8〕-X, Z〔9〕-Y): J=J-U:
J＜0＝＞J=J+360: ≠＞J=J ⊿
J＜180＝＞Z〔11〕=1: ≠＞Z〔11〕=-1⊿ “OK”:
Z〔11〕*I:Pause 0:S”S=”: Pause 0:J-＞DMS▲
GOTO 0: ≠＞S=S+ Z〔10〕:Z=1＝＞GOTO Y: ≠＞Z=2＝＞GOTO Z⊿GOTO 0 ⊿ LbI 9: K:S=S+K:GOTO X
