曲线任意里程中边桩坐标正反算(CASIO_fx-5800P计算器)程序.doc

《曲线任意里程中边桩坐标正反算(CASIO_fx-5800P计算器)程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、曲线任意里程中边桩坐标正反算(CASIOfx-5800P计算器)程序曲线任意里程中边桩坐标正反算(CASIOfx-5800P计算器)程序一、程序功能        本程序由一个主程序(XLJS)和两个子程——正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子

2、程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIOfx-4500P计算器及CASIOfx-4850P计算器上运行。  特别申明：(1).适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm；                    (2).在引用该核心计算部分时，请注明来源。（3）该程序可计算标准、非标准曲线，但不可以计算与线路成斜交角的坐标二、源程序1.主程序(XLJS)"1.KS+E=>X-Y"："2.X-Y=>KS+E"?N：Lbl1："KS"?Z："E

3、"?S：Prog"XXX"：1÷P→C：(P-R)÷(2HPR)→D：180÷π→E：IFN=1:ThenGoto2：ElseGoto3Lbl2：Abs(Z-O)→W：Prog"SUB1"："X=":X◢"Y=":Y◢"FWJ="：F-90◢Prog"GC"Goto1Lbl3："XS"?X："YS"?Y：X→I：Y→J：Prog"SUB2"："KS="：O+W→Z◢"E=":S◢Prog"GC"Goto32.正算子程序(SUB1)0.1739274226→A：0.3260725774→B：0.069431

4、8442→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+ScosF→X：Y+SsinF→Y3.反算子程序(FSZ

5、)G-90→T：Abs((Y-V)cosT-(X-U)sinT)→W：S=0：Lbl0：Prog"SUB1"：T+QEW(C+WD)→L：(J-Y)cosL-(I-X)sinL→S：IFAbsS<10^-6:ThenGoto1：ElseW+S→W：Goto0Lbl1：S=0：Prog"SUB1"：(J-Y)÷sinF→S数据库：IFZ≤下一起点里程:Then本线路起点里程→O：本线路X→U：本线路Y→V：本线路方位角→G：本线路长度→H：本线路起点半径→P：本线路止点半径→R：1、-1→Q(左偏Q=-1，

6、右偏Q=1，直线段Q=0）:ElseIFZ≤下一起点里程:Then本线路起点里程→O：本线路X→U：本线路Y→V：本线路方位角→G：本线路长度→H：本线路起点半径→P：本线路止点半径→R：1、-1→Q(左偏Q=-1，右偏Q=1，直线段Q=0）:IFEnd高程主程序GCLbI1：Prog"高程子程序名GCZL"：AbsR→R：2tan-1（T/R）→O：πRO/180→I：IFZ≤J-T:ThenZ-J→A：0→M：Q→V：Goto2：ElseIFZ≤J-T+I/2:Then-T→A：Z-J+T→M：Q→V

7、：Goto2：ElseIFZ≤J-T+I:Then-T→A:J-T+I-Z→M:V→-D:Goto2：ElseIFZ>J-T+I:Then-2T-Z+I+J→A:0→M:-D→V:Goto2:IFEndLbI2:IFW≤0:ThenGoto3:Else-R→R:Goto3LbI3:180M/(πR)→P:R(1/cosP-1)→C:Abs(R×tanP)→B:"GC=":E+(A+B)V+C→K◢竖曲线要素子程序“GCZL“IFZ≤下一起点里程:Then交点里程→J:交点高程→E：交点前坡度含正负（7.x

8、x/100）→Q：交点后坡度含正负（7.xx/100）→D：偏向判定（W＝0凹，W＝1凸）→W：切线长→T：竖曲线半径→R：ElseIFZ≤下一起点里程:Then交点里程→J：····································：R竖曲线半径"1.KS+E=>X-Y":"2.X-Y=>KS+E":?N:"K1"?O:"X0"?U:"Y0"?V:"FWJ"?G:"XLCD"?H:"R1"?P:"R2"?R:?