路桥施工测量常用计算资料大全：5800计算器程序1.doc

FHXPSzhu（主程序） "1.ZS,2.FS" ?→Q 输入1正算，输入2反算 “NEW=0,OLD ≠0”?Z:If Z=0:Then “X0=”?A:“Y0=”?B:“C0=”?C:“1/R0=”?D:“1/RI=”?E:“SP=”?F:“EP=”?G:Ifend:Q=2=>Goto2 Lbl 1 :“KM=,<0 STOP”?H:H<0=>STOP:“PJ=”?O:“PY=”?L Lbl Z:Z=1=> Prog“01”：Z=2=> Prog“02” :Fix 3 H- F→M：（E-D）/（ G-F）*M→N:C+（N+2 D）M*90/π→P：P<0=>P+360→P A+M(cos(C)+4∑cos(C+7.5*(2K+1)M((2K+1)N/12 +2D)/π),K,0,5)+2∑(cos(C+15KM(KN/6+2D)/π),K,1,5)+cos(P))/36+Lcos(P+O)→X B+M(sin(C)+4∑sin(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)→Y Q=2=>Goto 4 "X=":X◢ "Y=":Y◢ “FWJ=”：P▶DMS◢ Norm 2：Goto 1 Lbl 2:“XD=,<0,STOP”？R：R<0=>STOP:“YD=”？S “KMDG=”？H ：90→O:0→L：Goto Z （H线路范围内的任意桩号） Lbl 4:Pol（R-X，S-Y）：J<0 => J+360→J While abs（Icos（J-P））≤0.001:P-J>180=> J+360→J： P-J<-180=> P+360→P:IF P-J>0:then -I→L:else I→L ifend Goto 3: Whileend：H+Icos（J-P）→H:Goto Z Lbl 3:Fix 3 “ZHUANGHAO”:H◢ “JULI”:L◢ Norm 2：Goto 2 01（数据库子程序） If H<=第一曲线终点桩号：then 第一曲线起点X→A：第一曲线起点Y→B：第一曲线起点方位角→C：起点曲率→D：终点曲率→E：起点桩号→F：终点桩号→G：return：ifend …………… 程序说明： 1、该程序可以计算任意线形（直线、圆曲线、缓和曲线、不完整曲线）任意桩号的坐标（正算，输入1），也可根据坐标计算该点到线路的距离及垂足桩号（反算，输入2）； 2、（NEW=0，OLD≠0）？如果要计算的点为数据库线路中的点，则输入数据库编号（以整数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，已经可以满足精度要求；反算子程序为角度趋近的方法编写，计算速度有点慢。