4850计算器坐标正反计算程序.doc

曲线任意里程中边桩坐标正反算(CASIOfx-4850P计算器)程序(修改版):本程序由一个主程序(XYZ)与两个子程——正算子程序(A)、反算子程序(B)序构成,可以根据曲线段——直线、圆曲线、缓与曲线(完整或非完整型)线元要素(起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径)及里程边距或坐标,对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分两个子程序移植到其它相关程序中,用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIOfx-4500P计算器及CASIOfx-4850P计算器上运行。:利用Gauss-Legendre5点通用公式计算线路中边桩坐标并计算放样数据。源程序主程序(XY-ZD)Defm2:F=1:(F=1正算,F=2反算)Z[1]=90(右边线与路中线右夹角)Prog”THB”:Fix3:”X=”:X▲”Y=”:Y▲(F=2时,为“KM=“:Z▲”D=“:D▲)2次程序:THBZ“KM“:A=第一段线终点桩号:B=第二段线终点桩号:C=第三段线终点桩号:E=第四段线终点桩号:F=第五段线终点桩号:……(有多少段就加多少段;超过十五段,要另加子程序转过来)ZProg”C1”:≠>ZProg”C2”:≠>ZProg”C3”:≠>ZProg”C4”:≠>ZProg”C5”ΔΔΔΔΔ……(有几段线,则打几个三角)C=1÷P:S=(P-R)÷(2HPR):E=180÷π:F=1=>Goto1:≠>Goto2Δ←Lbl1:D:W=Z-O:Prog"A":Goto3Lbl2:{XY}:XY:I=X:J=Y:Prog"B":Z=O+W:D=D:Goto3Lbl3:F=1=>X:Y:≠>Z:D:(A)A=:B=:N=:K=:L=:M=:X=U+W(Acos(G+QEKW(C+KWS))+Bcos(G+QELW(C+LWS))+Ncos(G+QEMW(C+MWS))+Bcos(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))):Y=V+W(Asin(G+QEKW(C+KWS))+Bsin(G+QELW(C+LWS))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)WS))):Z[2]=G+QEW(C+WS)+Z[1]:X=X+DcosZ[2]:Y=Y+DsinZ[2](B)T=G-90:W=Abs((Y-V)cosT-(X-U)sinT):D=0:Lbl0:Prog"A":L=T+QEW(C+WS):Z=(J-Y)cosL-(I-X)sinL:AbsZ<1E-6=>Goto1:≠>W=W+Z:Goto0Δ←┘Lbl1:Z=0:Prog"A":Z=(J-Y)÷sinZ[2]:C1~C5…….O=:U=:V=:G=:H=