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

曲线任意里程中边桩坐标正反算(CASIO fx-4850P计算器)程序(修改版)
程序功能及原理
:本程序由一个主程序(XYZ)和两个子程——正算子程序(A)、反算子程序(B)序构成,可以根据曲线段——直线、圆曲线、缓和曲线(完整或非完整型)的线元要素(起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径)及里程边距或坐标,对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中,用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。
:利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据。
源程序
主程序(XY-ZD)
Defm 2:F=1:(F=1正算,F=2反算)Z[1]=90(右边线与路中线右夹角)
Prog”THB”:Fix 3:”X=”:X▲”Y=”:Y▲(F=2时,为“KM=“:Z▲”D=“:D▲)
2次程序:THB
Z“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=>Goto 1:≠>Goto 2Δ←
Lbl 1:D:W=Z-O:Prog "A":Goto 3
Lbl 2:{XY}:XY:I=X:J=Y:Prog "B":Z=O+W:
D=D:Goto 3
Lbl 3:F=1=>X:Y :≠>Z:D:
3. 正算子程序(A)
A=:B=:N=:K=:L=:M=:X=U+W(Acos(G+QEKW(C+KW
S))+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+Dcos Z[2]:Y=Y+Dsin Z[2]
4. 反算子程序(B)
T=G-90:W=Abs((Y-V)cosT-(X-U)sinT):D=0:Lbl 0:Prog "A":L=T+QEW(C+WS):Z=(J-Y)cosL-(I-X)sinL:AbsZ<1E-6=>Goto1:≠>W=W+Z:Goto 0Δ←┘
Lbl 1:Z=0:Prog "A":Z=(J-Y)÷s