5800全线坐标正算.docVIP

GAUSSLE坐标正反算fx－4850程序源程序正算主程序??ZS“K”？K：“L(-Z??+Y)”？L：90→M：（注：此处若不给M赋值，则可计算斜交点）Prog”DAT” ： (P-R)÷(2(H-O)PR) →D Abs(K-O) →J：ProgSUB1” “X=”：X Y=”： Y正算子程序 SUB1Dimz 0.1184634425→A： 0.2393143352→B： 0.2844444444Z[4]： 0.0469100770： 0.2307653449： 0.5Z[1] I+J(Acos(G+QCJ(1÷P+CJD)×180÷π)+Bcos(G+QEJ(1÷P+EJD)×180÷π)+Z[4]cos(G+QZ[1]J(1÷P+Z[1]JD)×180÷π)+Bcos(G+Q(1-E)J(1÷P+(1-E)JD)×180÷π)+Acos(G+Q (1-C)J(1÷P+(1-C)JD) ×180÷π)) →X S+J(Asin(G+QCJ(1÷P+CJD)×180÷π)+Bsin(G+QEJ(1÷P+EJD)×180÷π)+Z[4]sin(G+QZ[1]J(1÷P+Z[1]JD)×180÷π)+Bsin(G+Q(1-E)J(1÷P+(1-E)JD)×180÷π)+Asin(G+Q (1-C)J(1÷P+(1-C)JD) ×180÷π)) →Y G+QJ(1÷P+JD) ×180÷π+M→F： X+LcosF： Y+LsinF. 曲线元要素数据库：DAT线元起点的X坐标线元起点的坐标线元起点线元起点线元点线元起点曲率半径线元点曲率半径线 元左右偏标志线元起点的X坐标线元起点的坐标线元起点线元起点线元点线元起点曲率半径线元点曲率半径线 元左右偏标志………… (注：如有多个曲线元要素继续添加入数据库DAT-M中) I=线元起点的X坐标：S=线元起点的Y坐标：O=线元起点里程：G=线元起点切线方位角：H=线元终点里程：P=线元起点曲率半径：R=线元止点曲率半径：Q=线 元左右偏标志(注：如有多个曲线元要素继续添加入数据库DAT中)?说明：??K ？ 正算时所求点的里程L ？??正算时所求点距中线的边距(左侧取负值，右侧取正值，在中线上取零)??X ？ 反算时所求点的X坐标?????? Y ？ 反算时所求点的Y坐标????M ???斜交右角说明：一、程序功能及原理1.功能说明：本程序由两个主程序——正算主程序(GSZS)、反算主程序(GSFS)和两个子程——正算子程序(SUB1)、线元数据库(DAT-M)构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本程序可以在CASIO fx-4800P计算器及 CASIO fx-4850P计算器上运行。由于加入了数据库(DAT-M)，可实现坐标正反算的全线贯通。????组合程序5可实现M线的正算贯通，组合程序7可实现M线的反算贯通，组合程序6可实现坐标计算到放样一体化。2．计算原理：利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。????利用待求点至线元起点切线作垂线，逐次迭代趋近原理反算里程及边距。二、使用说明1、规定??????(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时， Q= -1；当线元往右偏时，Q=1；当线元为直线时，Q=0。??????(2) 当所求点位于中线时，L=0；当位于中线左侧时，L取负值；当位于中线右侧时，L取正值。??????(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。????(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。?????? (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。(7)曲线元要素数据库（DAT-M）可根据线型不同分为各个线元段输入到DAT-M中，即分为直线段、缓和曲线、圆曲线等。 (8)正算时可仅输入里程和边距及右交角可实现全线计算，但反算时只能通过首先输入里程K值读取数据库DAT-M，计算器自动将里程K所在线元数据赋