任意曲线任意里程中边桩坐标正反算程序(适用于casio5800计算器)

《任意曲线任意里程中边桩坐标正反算程序(适用于casio5800计算器)》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、CASIO5800计算器程序20080610版一、程序功能：本程序由一个主程序（MAIN）及两个子程序：坐标正算程序（SUB1）、坐标反算程序（SUB2）构成，可以根据曲线段的组成部分：直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、止点坐标、起点切线方位角、线元长度、起点曲率半径、止点曲率半径、里程、边距或坐标），对该曲线段范围A任意里程屮边桩坐标进行正反算。修改部分说明：这次修改将主程序中的注释语句全部取消。二、源程呼：1、主程序（MAIN）3->DimZ：n（l.SZ-XY,2.XY-SZ）n：nN=n?-*N：IfN=l：ThenGoto1:ElseGo

2、to4：IfEnd说明：DimZ是为程序增加变量，l.SZ—XY是利用里程（S）和边距（Z）计算坐标（XY）;2.XY-SZ是利用坐标（XY）计算里程（S〉和边距（Z）。Lbl1：nS=n?-*S："Z=n?-*Z："ANGLE"?-*T：IfT=90：ThenGoto2：ElseGoto3:IfEnd-Lbl2：ProgMDATln：Abs（S-O）->W：ProgHSUBln：nXS=n：X遑nYS=":YJGoto5—Lbl3:Prog^DATl"：Abs（S-O）-*W：ProgMSUB1M：UXS="："YS=":AGoto5说明：正算时调用语句，提示输入里程（S

3、）、边距（Z）,进入线元要素数据库，进入正算子程序进行计算，并显示坐标（X、Y）。进入第三段程序语句。F:”FWJ"=F-90—DMSJ句在进行曲线位置构造物放样计算时，可以直接显示切线方位角，这样的话就可以进行任意构造物的放样测量了。T"ANGLE"句及后面的判断语句构成斜交构造物中轴线的测量计算，其中的T如果是计算左右边桩的话输入90度即可，如果是斜交的屮轴线坐标计算的话，直接输入线路前进方向与构造物屮轴线的右交角即可。另外计算切线方位角的话必须按中桩计算。Lbl4："X="?-*X：nY="?-Y：X-*I：Y-J：ProgHDATln：PrognSUB2n：Fix3：

4、"S=":O+W^uZ=n：Z2Goto5*-说明：反算吋调用语句，提示输入坐标（X、Y）,并赋值，进入线元要素数据库，进行反算子程序，取小数点后三位数，并显示里程（S）、边距（Z）。转入第二段程序语句，重复执行程序。Lbl5：0-J：Pol（X-Z[l],Y-Z[2]）：J->F：IfF<0：ThenF+360-F：ElseIfF>360：ThenF-360-*F:IfEnd：IfEnd-—说明：进行极坐标放样的角度及距离计算，角度以度分秒形式显示。Z[l]、Z[2]为测站点坐标，自动变量。这句还需要改正。其实这句也没什么大用，有了坐标了，直接输入仪器进行放样就可以了，不必

5、要非得计算偏角和斜距吧。2、正算子程序（SUB1）0.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)+T-*F

6、：X+ZcosF—X:Y+ZsinF—Y:Return说明：利用复化辛普森公式进行坐标正算。输入即可。已经进行过复核。3、反算子程序(SUB2)G-90-*Z[3]：Abs((Y-V)cosZ[3]-(X-U)sinZ[3])-*W：0->Z：Lbl0：Prog"SUBl":Z[3]+QEW(C+WD)-L：(J-Y)cosL-(I-X)sinL-Z：IfAbsZZ：PrognSUBlu：(J-Y)4-sinE->Z：Return说明：利用坐标正算的公式进行反算，并给出计算点的屮桩里程

7、及坐标。4、曲线元要素数据库(DAT1)IfS<止点里程桩号AndS&起点里程桩号：ThenXXX.XXX-U：XXX.XXX-V：XXX.XXX-O：XXX.XXX-G：XXX.XXX-H：XXX.XXX一P:XXX.XXX-Q：IfEnd-IfS＜止点里程桩号And8＞起点里程桩号：ThenXXX.XXX-U：XXX.XXX-V：XXX.XXX-O：XXX.XXX一G:XXX.XXX-H：XXX.XXX—P:XXX.XXX-Q：IfEnd*-14-P-*C：（P-R）+（2HPR）-*D：1804-