4850计算器坐标正反计算程序.docVIP

曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序（修改版） 程序功能及原理1.功能说明：本程序由一个主程序(XYZ)和两个子程——正算子程序(A)、反算子程序(B)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。2．计算原理：利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。 源程序 主程序(XY-ZD) Defm 2:F=1：(F=1正算，F=2反算)Z[1]=90（右边线与路中线右夹角） Prog”THB”:Fix 3:”X=”:X▲”Y=”:Y▲（F=2时，为“KM=“：Z▲”D=“：D▲） 2次程序：THB Z“KM“:A=第一段线终点桩号：B=第二段线终点桩号：C=第三段线终点桩号：E=第四段线终点桩号：F=第五段线终点桩号：……（有多少段就加多少段;超过十五段，要另加子程序转过来） ZA=Prog”C1”:≠ZB=Prog”C2”: ≠ ZC=Prog”C3”: ≠ ZE=Prog”C4”: ≠ ZF=Prog”C5” ΔΔΔΔΔ……（有几段线，则打几个三角） C=1÷P：S=(P-R)÷(2HPR)：E=180÷π：F=1=Goto 1：≠Goto 2Δ←Lbl 1：D：W=Z-O：Prog A：Goto 3 Lbl 2：{XY}：XY：I=X：J=Y：Prog B：Z=O+W：D=D：Goto 3 Lbl 3：F=1=X:Y :≠Z:D:3. 正算子程序(A)A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:M=0.5:X=U+W(Acos(G+QEKW(C+KWS))+Bcos(G+QELW(C+LWS))+Ncos(G+QEMW(C+MWS))+Bcos(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))): Y=V+W(Asin(G+QEKW(C+KWS))+Bsin(G+QELW(C+LWS))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)WS)))：Z[2]=G+QEW(C+WS)+Z[1]：X=X+Dcos Z[2]：Y=Y+Dsin Z[2]4. 反算子程序(B)T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：D=0：Lbl 0：Prog A：L=T+QEW(C+WS)：Z=(J-Y)cosL-(I-X)sinL：AbsZ1E-6=Goto1：≠W=W+Z：Goto 0Δ←┘Lbl 1：Z=0：Prog A：Z=(J-Y)÷sinZ[2]5.曲线元要素数据库：C1~C5…….O=xxxx.xxx:U=xxxx.xxx:V=xxxx.xxx:G=xxxx.xxx:H=xxxx.xxx:P=xxxx.xxx:R=xxxx.xxx:Q=x (注：一段线一个子程序，分别提供这几个要素数据。)三、使用说明1、规定????????(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。????????(2) 当所求点位于中线时，D=0；当位于中线左铡时，D取负值；当位于中线中线右侧时，D取正值。????????(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。????????(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。???????? (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。???????? (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径(7)曲线元要素数据库（C1~Cn）可根据线型不同分为各个