坐标反算（源程序）
File ZBFS
1: M“X0” : N“Y0” : A“PJ” : R : L“LS” : T : O“L” : H“ZH” : F“FWJ”
2: Z[1]=P“XJD”+T cos ( F+180) : Z[2]=K“YJD”+T sin ( F+180)
3: U“R+1,L-1”
4: Z[3]=H+L : Z[4]=H+O－L : Z[5]=H+O
5: Lbl 0 : C“LC” : Z[7]=Abs(C－H)
6: Prog“KL”
7: S=－((X－M)sin(Q+90)－(Y－N) cos (Q+90) : Abs S<0.0001 C“ LC ” ◢ W=“LP ”= ( X－M ) sin Q－(Y－N) cos Q◢ C = C+S : Goto 0 △

File KL
1: C≤H Z[8]=－Z[7] : Z[9] = 0 : Z[10] = 0 : Goto 2 △
2: C≤Z[3] Z[8]= Z[7]－Z[7] ^5÷40÷R2÷L2 : Z[9] = Z[7]^3÷6÷R÷L－Z[7] ^7÷336÷R^3÷L^3 : Z[10] = 90 Z[7] 2÷R÷L÷Л : Goto 2 △
3: C≤Z[4] Z[10] = 90 L÷R÷Л +180( Z[7] －L)÷R÷Л : Rsin (Z[10]+0.5L－L^3÷240÷R2 : Z[9] = R(1－cosZ[10] )+L2÷24÷R : Goto 2 △
4: C≥Z[5] X=－( Z[7] －O) : Y = 0 : Z[10] = A : Goto 1 △
5: C≥Z[4] Z[10]= A－90 (O－ Z[7] ) 2 ÷R÷L÷Л :X = (O－Z[7] )－(O－Z[7] )^5÷40÷R2÷L2 : Y = (O－Z[7] )^3÷6÷R÷L－(O－Z[7] )^7÷336÷R^3÷L^3 △
6: Lbl 1 : Z[8] = T+(T－X) cos A－Y sin A : Z[9] =(T－X) sin A+ Y cos A : Goto 2 △
7: Lbl 2 : Pol (Z[8],Z[9]) : Z[7] =F+JU : Q=F+Z[10] U : X=Z[1] + IcosZ[7] : Y=Z[2] + I sin Z[7]

注： M、N：已知坐标 A：偏角 R：半径 L：缓和曲线长 T：切线长 O：曲线长 H：直缓点桩号 F：第一切线方位角 U：路线前进方向左“-”右“+”



整理自本站论坛　谢谢　火溪河　提供