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

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

1、曲线任意里程中边桩坐标反正算（CASIOfx-5800P计算器）程序一，程序功能本程序由一个主程序（ZHUCHENXU）和几个子程序---正算子程序（SUB1）反算子程序（SBU2）数据库（SBU3）,算方位角程序（FA），使用方位角算坐标（NE），边坡放样程序（BIANPO）,隧道圆心放样程序（SDYX）构成，可以根据直线，圆曲线，缓和曲线（完整或非完整型）的线元要素（起点坐标，起点里程，起点切线方位角，线元长度，起点曲率半经，止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩进行正反算。二，源

2、程序增加边量22→DimZ1,主程序（ZHUCHENXU）’’1.SZ→XY,2.XY→SZ,3.SDYX,4BIANPO’’?N:Lb11:’’ZHUANGHAO=’’?S:Prog’’SUB3’’:N≥2=>Goto2:Abs(S-0)→W:’’BIANZHU=’’?→Z:Prog’’SUB1’’:’’XS=’’:X→X◢’’YS=’’:Y→Y◢’’FS=’’:(F-90)→F:F▼DMS◢Prog’’FA’’:Goto1:Lb12:’’CX=’’?X:X→I:’’CY=’’?Y:Y→J:’’DMG=’’

3、?→Z[7]:Prog’’SUB2’’:’’S=’’:(O+W)→S◢’’Z’’:Z→Z◢IfN=3:ThenProg’’SDYX’’:IfEnd:IfN=4:ThenProg’’BIANPO”:IfEnd:Goto22，正算子程序（SUB1）0.1739274226→Z[1]:0.3260725774→B:0.0694318442→K:0.3300094782→L:(1-L)→F:(1-K)→M:U+W(Z[1]cos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QE

4、FW(C+FWD))+Z[1]cos(G+QEMW(C+MWD)))→X:V+W(Z[1]sin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Z[1]sin(G+QEMW(C+MWD)))→Y:G+QEW(C+WD)+90→F:X+ZcosF→X:Y+ZsinF→Y3,反算子程序（SUB2）G-90→Z[2]:Abs((Y-V)cosZ[2]-(X-U)sinZ[2])→W:0→Z:Lb10:Prog’’SUB1’’:Z[2]+QEW(C+WD

5、)→L:(J-Y)cosL-(I-X)sinL→Z:IfAbs(z)＜1E-6:ThenGoto1:ElseW+Z→W:Goto0:IfEndLb11:0→Z:Prog’’SUB1’’:(J-Y)÷sinF→Z4,数据库（SUB3）IfS＜226100:ThenCls:Stop:IfEndIfS＜226255.833:Then226100→0:99037.9736→U:4282.3590→V:201.3266304→G:1045→P:1045→R:253.543→H:0→Q:Goto1:IfEnd…………..

6、……………IfS＞236005.651:ThenCls:IfEndLb11:1÷P→C:(P-R)÷(2HPR)→D:180÷∏→E:5．算方位角程序（FA）LbI0:‘‘XC=’’?A:‘‘YC=’’?T:X:Y0→J:POI(X-A,Y-T)→I:‘‘BC=’’:I◢IfJ＜0:Then360+J→J:ElseJ→J:IfEnd‘‘FW=’’:J▼DMS◢Goto06.使用方位角算坐标（NE）Lbl0：‘‘X=’’?X:‘‘Y=’’?Y:‘‘L=’’?L:‘‘F0’’?F:L×Cos(F)+X→N:‘‘N

7、=’’N◢L×Sin(F)+Y→E:‘‘E=’’E◢Goto07.边坡放样程序（BIANPO）Prog’’TCD’’:Prog’’BHD’’:IFS＜226255.731:then232.257→Z[17]:225844.269→Z[18]:-0.69→Z[19]:11789.745→Z[22]:Goto3:IfEndIFS＜226351.25:then230.837→Z[17]:226050→Z[18]:2.8→Z[19]:10^45→Z[22]:Goto3:IfEndIFS＜226688.75:then2

8、39.272→Z[17]:226351.25→Z[18]:2.8→Z[19]:-45000→Z[22]:Goto3:IfEndIFS＜227318:then243.997→Z[17]:226520→Z[18]:2.05→Z[19]:10^45→Z[22]:Goto3:IfEndIFS＜227982:then260.356→Z[17]:227318→Z[18]:2.05→Z[19]:-16000→Z[