曲线任意里程中边桩正反算5800P(修改版).doc

线路中边桩任意里程坐标正反算程序5800PMain（主程序名）Lbl 0：Cls1.SZ=NE 2.NE=SZ:Q=1 ?or 2?:?Q Q=1 = Goto1:Q=2 = Goto2Lbl 1:Cls:DK=?SS=0= Goto 0:ProgDATA （注：S=0时，程序会由Lbl 0开始重新运算）Lbl 1:Z=?Z:ProgQXJSSUB1(注：输入正算里程)Fix 4:Cls (注：小数点后保留4位有效数字)X=:Locate 3,1,N (注：输出正算X值)Y=:Locate 3,2,E(注：输出正算Y值)FS=:F DMS (注：DMS为一个符号,以度分秒格式显示，正算切线方位角)ProgFYJS:Goto 1Lbl 2: Cls:N=? B (注：输入反算点坐标X值)B=0= Goto 0:E=? C: Predict DK=?S: BN: CE:ProgQXJSSUB2: S=:SZ=:Z Goto 2（注：B=0时，程序会由Lbl 0开始重新运算，输入反算点坐标Y值，其中英文Predict/pridikt/中文含义“预计”，Predict DK=预计反算里程（呵呵，班门弄斧一下），作用：以便从程序库DATA提取相应的线元要素，否则计算器不停的计算，造成假死机）【主程序完，主程序经此改动后,即可加入测站点坐标进行距离与方位角的放样了】 FYJS （注：方位角和距离放样子程序名,与主程序ProgFYJS:Goto 1对应）ClrStatXJ=?G:YJ=?H （注：XJ,输入测站点X值，YJ，输入测站点Y值）Pol(N-G,E-H):ClsIf J0:ThenJ+360Y:Else JY:IfEndFY HD=:,I （注：字母I，正算值与测站点平距）FY HR=: YDMS (注：DMS度分秒格式，正算值至测站点方位角)【方位角和距离放样子程序完】DATA (数据库子程序名)Goto1 同时保存多个曲线时的指针，通过更改序号（如：1改成2）变换曲线Lbl 1 第一条曲线（如：左线数据库）IF S*（线元终点里程）:Then*A（线元起点方位角）:*O（线元起点里程）:*U（线元起点X）:*V（线元起点Y）:*P（线元起点曲率半径）:*R（线元终点曲率半径）: *L（线元起点至终点长度）: Return:IfEndbl 2 第二条曲线（如：右线数据库）IF S*（线元终点里程）:Then*A（线元起点方位角）:*O（线元起点里程）:*U（线元起点X）:*V（线元起点Y）:*P（线元起点曲率半径）:*R（线元终点曲率半径）: *L（线元起点至终点长度）: Return:IfEndIF S0.0001 = Goto1Lbl 2: 0Z：ProgQXJSSUB1:(C-E) sin(F+90) Z【反算子程序完】【使用说明】1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当曲线半径在左时，P、R取负值，当曲线半径在右时，P、R取正值，当曲线半径为无穷大（即直线）时，P、R以10的45次代替。(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。(3) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。 (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。 (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明 输入部分：1. SZ = NE2. NE = SZ Q ?选