1、FHXPSzhu（主程序）1.ZS,2.FS ?Q 输入1正算，输入2反算 “NEW=0,OLD0”?Z:If Z=0:Then “X0=”?A:“Y0=”?B:“C0=”?C:“1/R0=”?D:“1/RI=”?E:“SP=”?F:“EP=”?G:Ifend:Q=2=Goto2Lbl 1 :“KM=,0 STOP”?H:HSTOP:“PJ=”?O:“PY=”?LLbl Z:Z=1= Prog“01”：Z=2= Prog“02” :Fix 3H- FM：（E-D）/（ G-F）*MN:C+（N+2 D）M*90/P：PP+360PA+M(cos(C)+4cos(C+7.5*(2K+1)M(2K
2、+1)N/12 +2D)/),K,0,5)+2(cos(C+15KM(KN/6+2D)/),K,1,5)+cos(P)/36+Lcos(P+O)XB+M(sin(C)+4sin(C+7.5*(2K+1)M(2K+1)N/12 +2D)/),K,0,5)+2(sin(C+15KM(KN/6+2D)/),K,1,5)+sin(P)/36+Lsin(P+O)YQ=2=Goto 4X=:XY=:Y“FWJ=”：PDMSNorm 2：Goto 1Lbl 2:“XD=,0,STOP”？R：RSTOP:“YD=”？S“KMDG=”？H ：90O:0L：Goto Z （H线路范围内的任意桩号）Lbl 4:Po
3、l（R-X，S-Y）：J J+360JWhile abs（Icos（J-P）0.001:P-J180= J+360J： P-J P+360P:IF P-J0:then -IL:else IL ifendGoto 3: Whileend：H+Icos（J-P）H:Goto ZLbl 3:Fix 3“ZHUANGHAO”:H“JULI”:LNorm 2：Goto 201（数据库子程序）If H=第一曲线终点桩号：then 第一曲线起点XA：第一曲线起点YB：第一曲线起点方位角C：起点曲率D：终点曲率E：起点桩号F：终点桩号G：return：ifend程序说明：1、该程序可以计算任意线形（直线、圆曲
4、线、缓和曲线、不完整曲线）任意桩号的坐标（正算，输入1），也可根据坐标计算该点到线路的距离及垂足桩号（反算，输入2）；2、（NEW=0，OLD0）？如果要计算的点为数据库线路中的点，则输入数据库编号（以整数1、2、3代替输入）；如果在数据库中没有要计算线路的数据，则输入0，并按提示输入曲线要素（X0：曲线起点X坐标，Y0：曲线起点Y坐标，C0:曲线起点方位角，1/R0:曲线起点曲率，直线为0， 1/RI:曲线终点曲率，直线为0，曲率左偏输入负值，右偏输入正值，SP:曲线起点桩号，EP：曲线终点桩号），KM为待求点桩号，PJ为待求点与线路的偏角（正交输入90，斜交输入实际角度），“PY”为偏中距离，线路上的点输入0，右偏输入+值，左偏输入-值；3、 正算显示坐标及切线方位角；反算输入线路的任意桩号（此桩号越接近真实值计算速度越快）、待求点坐标，显示待求点桩号及偏中距离；正算子程序为复化辛普森公式编写而成，此时线元段数N取6，已经可以满足精度要求；反算子程序为角度趋近的方法编写，计算速度有点慢。
安全达人 