FX5800P全线任意里程中边桩坐标正反算 曲线任意里程中边桩坐标正反算.doc

FX5800P全线任意里程中边桩坐标正反算（加了数据库）本程序非本人原创，结合上传者长处，进行改编，愿收获和大家共享，不足之处多多指教。希望过路神仙把竖曲线加进去带算隧道超欠挖，谢谢曲線任意里程中邊樁坐標正反算5800p本程序由内蒙@文哥和师父阿玖哥共同改编完成主程序ZBJSLbl6："1.SZ→XY"："2.XY→SZ"：？S：？N：Prog“SUB0”1÷P→C：(P-R)÷(2HPR)→D：180÷π→E:N=1=>Goto1:Goto2Lbl1:?Z:Abs(S-O)→W:Prog"SUB1":"XS=":X◢"YS=":Y◢"FS=":F-90►DMS◢“XC”?C：“YC”?E：0→I：0→J：Pol（X-C，Y-E）:“I=”：I◢IfJ∠0:Then“J=”:J+360►DMS◢IfEndIfJ>0：Then“J=”：J►DMS◢IfEnd：Goto6Lbl2?X:?Y:X→I:Y→J:Prog"SUB2":O+W→S:Cls:"S=":S◢"Z=":Z◢Goto6正算子程序SUB10.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L：1-L→F：1-K→M：U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD)))→X：V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD)))→Y：G+QEW(C+WD)+90→F：X+Zcos（F）→X：Y+Zsin（F）→Y反算子程序SUB2G-90→T(Y-V)cosT-(X-U)sin(T)→WAbs(W)→W:0→ZLbl6：Prog"SUB1"T+QEW(C+WD)→L：(J-Y)cos（L）-(I-X)sin（L）→ZIFAbs（Z）<1E-6：Then0→Z:Prog"SUB1":(J-Y)÷sin（F）→Z:ElseW+Z→W：Goto6:IfEnd数据库子程序SUB0Goto1(线元可输入多条，分离式可在前多加一位，匝道一样。例：左幅为K129+500,右幅输线元参数里程为1129+500，其他不变，前面1为任意数字，计算机便于区分)Lbl1：IFS<线元终点里程：Then@@@→O(线元起点里程)：@@@→U(线元起点X坐标)：@@@→V(线元起点Y坐标)：@@@→G(线元起点计算方位角)：@@@→P(线元起点半径)：@@@→R(线元止点半径)：@@@→H(线元长度)：@@@→Q(线元左、右偏标志，左偏-1，右偏1，直线为0)：Return：IfEndIFS<线元终点里程：Then@@@→O(线元起点里程)：@@@→U(线元起点X坐标)：@@@→V(线元起点Y坐标)：@@@→G线元(起点计算方位角)：@@@→P(线元起点半径)：@@@→R(线元止点半径)：@@@→H(线元长度)：@@@→Q(线元左、右偏标志，左偏-1，右偏1，直线为0)：Return：IfEnd。。。。。。。。。。。。。。。。。。。。。。。。。。。一程序功能本程序由一个主程序(ZBJS)和3个子程——正算子程序(SUB1)、反算子程序(SUB2)、数据库子程序（SUB0）构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIOfx-4500P计算器及CASIOfx-4850P计算器上运行。二、使用说明1、规定(1)以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。(2)当所求点位于中线时，Z=0；当位于中线左侧时，Z取负值；当位于中线中线右侧时，Z取正值。(3)当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。(4)当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。(5)当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。(6)当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止