4800计算程序(精).docVIP

曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序 一、程序功能及原理1.功能说明：本程序由一个主程序(TYQXjs)和四个子程正算子程序(SUB1)、反算子程序(SUB2)等构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb0=该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。二、源程序1.主程序(TYQXjs)Lbl 0：Defm 50:1.SZ = XY,2.XY = SZ：N： N=1=Goto 1：≠Goto 2Δ←Lbl 1：{SZ}：SZS≤Z[18]=J=1：Prog “DAT1”：Goto 3ΔS≤Z[26]=J=2：Prog “DAT1”：Goto 3Δ???????????? ……………S≤Z[8（N+1）+2]=J=n：Prog “DAT1”：Goto 3ΔLbl 3：W=Abs(S-O)：Prog SUB1：XS=:X=XYS=:Y=Y◢”FWI=”: F=F-M :F→DMS◢ Goto 4←┘Lbl 2：{XY}：XY：Z[4]=X：Z[5] =Y：N=0Lbl??A：Isz N：A=Z[8N+3]-M：B=Z[8（N+1）+3]-M：Prog ZX1 ：Z[6]×Z[7]0=J=N：Prog DAT1 ：Goto BΔA=Z[8N+3]+M：B=Z[8（N+1）+3]+M：Prog ZX1 ：Z[6]×Z[7]0=J=N：Prog DAT1 ：Goto B：≠ Goto AΔLbl??B：Prog SUB2：S=:S=O+WZ=:Z=Z◢ Goto 2Lbl4:J=0:I=Pol(X-Z[1],Y-Z[2]):F=J:F0=F=F+360Δ??”DIST=”:??I◢”FW=”:F→DMS◢ Goto 12. 正算子程序(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+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：Lbl 0：Prog SUB1：L=T+QEW(C+WD)：Z=(Z[5]-Y)cosL-(Z[4]-X)sinL：AbsZ1E-6=Goto1：≠W=W+Z：Goto 0Δ←┘Lbl 1：Z=0：Prog SUB1：Z=(Z[5]-Y)÷sinF4. 垂距计算子程序(ZX1)Z[6]=(Z[5]-Z[8N+1])COS A-(Z[4]-Z[8N])SIN AZ[7]=(Z[5]-Z[8(N+1)+1])COS B-(Z[4]-Z[8(N+1)])SIN B5.曲线元要素数据库：DAT1Lbl1: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=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中)Lbl 2 ：UX0：VY0：OS0：GF0：HLS：PR0：RRN：Q：C=1÷P：D=(