曲线任意里程中边桩坐标正反算
(CASIO fx-4850P
计算?/p>
)
程序
(
修改
版)
默认分类
2009-07-23 10:57:59
阅读
724
评论
1
字号:大中小
一、程序功能及原理
1.
功能说明:本程序由一个主程序
(TYQX
js
)
和两个子?/p>
—?/p>
正算子程?/p>
(SUB1)
、反算子程序
(SUB2)
序构成,可以根据曲线?/p>
—?/p>
直线、圆曲线、缓和曲线(完整或非完整型)的线元要素(起点坐标、起?/p>
里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径)及里程边距或坐标,对该曲线段范围
内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关?/p>
程序中,用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在
CASIO fx-4500P
计算器及
CASIO
fx-4850P
计算器上运行?/p>
2
.计算原理:利用
Gauss-Legendre 5
点通用公式计算线路中边桩坐标并计算放样数据?/p>
二、源程序
1.
主程?/p>
(TYQX
js
)
Defm 2:
"
1.SZ => XY
"
?/p>
"
2.XY => SZ
"
?/p>
N
?/p>
Prog”DAT1?U
"
X0
"
?/p>
V
"
Y0
"
?/p>
O
"
S0
"
?/p>
G
"
F0
"
?/p>
H
"
LS
"
?/p>
P
"
R0
"
?/p>
R
"
RN
"
?/p>
Q
?/p>
C=1÷
P
?/p>
D=(P-R)÷
(2HPR)
?/p>
E=180÷π
?/p>
N=1=>Goto 1
?/p>
?gt;Goto 2Δ?/p>
Lbl 1
?/p>
{SZT}
?/p>
SZ:T”ANG?/p>
?/p>
W=Abs(S-O)
?/p>
Prog
"
SUB1
"
?/p>
"
XS=
"
:X=X
?/p>
"
YS=
"
:Y=Y
?/p>
Goto 3←┘
Lbl 2
?/p>
{XY}
?/p>
XY
?/p>
I=X
?/p>
J=Y
?/p>
Prog
"
SUB2
"
?/p>
"
S=
"
:S=O+W
?/p>
"
Z=
"
:Z=Z
?/p>
Goto 2
Lbl 3:J=0:I=Pol(X-Z[1],Y-
Z[2]):F=J:F<0=>F=F+360:?gt;F>360=>F=F
-
36ΔΔ
”DIST=?
I
?/p>
”FW=?F→DMS
?/p>
Goto 1
2.
正算子程?/p>
(SUB1)
A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.230
7653449:M=0.5:F
=1-L:M=1-K:X=U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(G+QEMW(C+MWD))
+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+QEMW(C+MWD))+Bsin
(G+QE(1-L)W(C+(1-L)WD))+Asin (G+QE(1-K)W(C+(1-K)WD)))
?/p>
F=G+QEW(C+WD)+T
?/p>
X=X+ZcosF
?
Y=Y+ZsinF
3.
反算子程?/p>
(SUB2)
T=G-90
?/p>
W=Abs((Y-V)cosT-(X-U)sinT)
?/p>
Z=0
?/p>
Lbl 0
?/p>
Prog
"
SUB1
"
?/p>
L=T+QEW(C+