4800计算程序(精)

《4800计算程序(精)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、曲线任意里程中边桩坐标正反算(CASIOfx-4850P计算器)程序一、程序功能及原理1.功能说明：本程序由一个主程序(TYQXjs)和四个子程：正算子程序(SUB1)、反算子程序(SUB2)等构成，可以根据曲线段、直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。2．计算原理：利用Gauss-Legendre5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（

2、即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。二、源程序1.主程序(TYQXjs)Lbl0：Defm50:"1.SZ=>XY,2.XY=>SZ"：N：N=1=>Goto1：≠>Goto2Δ←Lbl1：{SZ}：SZS≤Z[18]=>J=1：Prog“DAT1”：Goto3ΔS≤Z[26]=>J=2：Prog“DAT1”：Goto3Δ            ……………S≤Z[8（N+1）+2]=>J=n：Prog“DAT1”：Goto3ΔLbl3：W=Abs(S-O)：Prog"SUB1"："XS=":X=X◢"YS=":Y

3、=Y◢”FWI=”:F=F-M:F→DMS◢Goto4←┘Lbl2：{XY}：XY：Z[4]=X：Z[5]=Y：N=0Lbl  A：IszN：A=Z[8N+3]-M：B=Z[8（N+1）+3]-M：Prog"ZX1"：Z[6]×Z[7]<0=>J=N：ProgDAT1"：GotoBΔA=Z[8N+3]+M：B=Z[8（N+1）+3]+M：Prog"ZX1"：Z[6]×Z[7]<0=>J=N：ProgDAT1"：GotoB：≠>GotoAΔLbl  B：Prog"SUB2"："S=":S=O+W◢"Z=":Z=Z◢Goto2Lbl4:J=0:I=Pol(X-Z[1],Y-Z[2]):

4、F=J:F<0=>F=F+360Δ  ”DIST=”:  I◢”FW=”:F→DMS◢Goto12.正算子程序(SUB1)A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:Z[3]=0.5:X=U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(G+QEZ[3]W(C+Z[3]WD))+Bcos(G+QE(1-L)W(C+(1-L)WD))+Acos(G+QE(1-K)W(C+(1-K)WD))):Y=V+W(Asin(G+QEKW(C+

5、KWD))+Bsin(G+QELW(C+LWD))+Nsin(G+QEZ[3]W(C+Z[3]WD))+Bsin(G+QE(1-L)W(C+(1-L)WD))+Asin(G+QE(1-K)W(C+(1-K)WD)))：M”ANG=”:F=G+QEW(C+WD)+M：X=X+ZcosF：Y=Y+ZsinF3.反算子程序(SUB2)M”ANG=”:T=G-M：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl0：Prog"SUB1"：L=T+QEW(C+WD)：Z=(Z[5]-Y)cosL-(Z[4]-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：G

6、oto0Δ←┘Lbl1：Z=0：Prog"SUB1"：Z=(Z[5]-Y)÷sinF4.垂距计算子程序(ZX1)Z[6]=(Z[5]-Z[8N+1])COSA-(Z[4]-Z[8N])SINAZ[7]=(Z[5]-Z[8(N+1)+1])COSB-(Z[4]-Z[8(N+1)])SINB5.曲线元要素数据库：DAT15Lbl1:J=1=>U=Z[8]:V=Z[9]:O=Z[10]:G=Z[11]:H=Z[12]:P=Z[13]:R=Z[14]:Q=Z[15]ΔJ=2=>U=Z[16]:V=Z[17]:O=Z[18]:G=Z[19]:H=Z[20]:P=Z[21]:R=Z[22]:Q

7、=Z[23]Δ                      ……………………………J=N=>U=Z[8N]:V=Z[8N+1]:O=Z[8N+2]:G=Z[8N+3]:H=Z[8N+4]:P=Z[8N+5]:R=Z[8N+6]:Q=Z[8N+7]Δ(注：如有多个曲线元要素继续添加入数据库DAT1中)Lbl2：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180÷π三、使用说明1、规定