缓和曲线任一点坐标计算程序

4800公路测量放样计算程序

CASIO fx-4800P计算器程序
DA DI——N E (公路弯道坐标计算)

K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”
M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（240R2）▲
N“HY”=M＋S▲
O“QZ”=M＋（πRAbs F÷180＋S）÷2▲
Z[1]“YH”=M＋πRAbs F÷180▲
Q“HZ”= Z[1]＋S▲
T=K－M▲
L=Q－M▲
E=（R＋S2÷（24R））÷cos（F÷2）－R▲
Lb1 0：：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△
：U“I——B OUT＋IN－”
Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋ S2÷（24R）－Ucos H：△
F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Tsin Abs F－Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsin Abs F－Tsin Abs F－Wcos F：△△
C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△
I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲
J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲
G=G＋20
Goto 0
程序运行
输入：JD？ ——交点桩号
R？ ——圆曲线半径
LS？ ——缓和曲线长度
PJ——L－R＋？ ——偏角 右偏为正
ZH——N？ ——ZH点N坐标
ZH——E？ ——ZH点E坐标
JD——N？ ——JD点N坐标
JD——E？ ——JD点E坐标
输出：ZH= ——ZH点桩号
HY= ——HY点桩号
QZ= ——QZ点桩号
YH= ——YH点桩号
HZ= ——HZ点桩号
T= ——切线长
L= ——曲线长
E= ——外矢距
再输入：ZHUANG HAO？ ——输入任一点桩号
I——B OUT＋IN－？ ——距中线距离（外侧为正）
再输出：COM——N= ——计算出N坐标
COM——E= ——计算出E坐标
以此循环——再输入——再输出
程序中各字母所代表的含义：

A——ZH点N坐标
B——ZH点E坐标
C——JD点N坐标
D——JD点E坐标
E——外矢距
F——偏角（右偏为正）
G——弯道内任一点桩号
H——弯道内任一点切线角
I——计算出的N坐标
J——计算出的E坐标
K——交点桩号
L——曲线长
M——ZH点桩号
N——HY点桩号
O——QZ点桩号
P——ZH--JD方位角
Q——HZ点桩号
R——圆曲线半径
S——缓和曲线长度
T——切线长
U——距中线距离（外侧为正）
V——支距
W——支距
X——弯道内坐标
Y——弯道内坐标
Z——任一点曲线长度
Z[1]——YH点桩号

另个程序
公路逐桩坐标计算程序
(可以计算对称、不对称缓和曲线)
Lb1 0
Z=？
V=？
W=V+2：Fixm
｛K｝
Lb1 1
K>Z[W+5Z+4]=>W=W+1:Goto 1⊿ (判断桩号在哪个交点范围，就是该交点曲线起点至下一交点曲线起点)
S=K-Z[W+5Z+3] （计算该桩号与曲线起点的距离）
R=Z[W+2Z+2]：L=Z[W+3Z+2]:E=Z[W+4Z+2] （读取该交点曲线要素R、Ls1 、Ls2）
Pol（Z[W]-Z[W-1],Z[W+Z+2]-Z[W+Z+1]） （计算该交点与下一交点直线方位角）
J<0=>J=J+360⊿
A=J
Pol（Z[W-1]-Z[W-2],Z[W+Z+1]-Z[W+Z]） （计算该交点与上一交点直线方位角）
J<0=>J=J+360⊿
C=A-J:A=J （计算偏角）
W=V+2=>Goto2⊿ （如果桩号在起点与第一交点曲线起点之间，则转 Lb1 2 ）
I=Abs（tan（c÷2））
M=L÷2-L^3÷240R^2:N=E÷2-E^3÷240R^2
P=L^2÷6R-L^4÷336R^3-R（1-cos（90L÷πR））
Q=E^2÷6R-E^4÷336R^3-R（1-cos（90E÷πR））
D=（P-Q）I÷2 : F=（P+Q+2R）I÷2
M=F+M-D:Q=F+N+D
N=πRAbsC÷180+（L+E）÷2
X=Z[W-1]-McosA
Y=Z[W+Z+1]-MsinA
M=Z[W-1]+Qcos（A+C）
V=Z[W+Z+1]+Qsin（A+C）
Q=AbsC÷C
S≤L=>P=0:Goto3⊿ （如果桩号在第一缓和曲线内，则转 Lb1 3）
S≤N-E=>S=S-L:Goto4⊿（如果桩号在圆曲线内，则转 Lb1 4）
S≤N=>S=N-S
Q=-Q:A=A+C-180:X=M:Y=V:L=E：P=180:Goto3 ⊿ （如果桩号在第二缓和曲线内，则转 Lb1 3）
P=A+C:S=S-N:D=M+ScosP:F=V+SsinP
Goto6 （如果桩号在直线内，则转 Lb1 6）
Lb1 2
P=A+C
D=Z[W-1]+ScosP
F=Z[W+Z+1]+SsinP:Goto6
Lb1 3
I=S-S^5÷40R^2÷L^2+S^9÷3456R^4÷L^4
J=Q（S^3÷6RL-S^7÷336R^3÷L^3）
P=P+A+90QS^2÷πRL:Goto5
Lb1 4
M=90（2S+L）÷πR
I=RsinM+L÷2-L^3÷240R^2
J=Q（L^2÷24R+R（1-cosM））
P=A+QM
Lb1 5
D=X+IcosA-jsinA:F=Y+JcosA+IsinA
Lb1 6
D〃X=〃◢ （结果显示X坐标）
F〃Y=〃◢ （结果显示Y坐标）
P〃AT=〃◢ （结果显示该桩号方位角）
：B〃S〃O〃⊿〃 （输入边桩距离，交角）
P=P+O
L〃XB〃=D+BcosP◢ （结果显示边桩X坐标）
M〃YB〃=F+BsinP◢ （结果显示边桩Y坐标）

以上是坐标计算程序，括号内是程序计算的大致原理及说明，中间部分为直线、圆曲线、缓和曲线计算的各种公式，大家也知道，书上也有。
该程序是最基本的，如再加几条语句就可以处理断链。
该程序还需变量输入，把整条线路曲线要素一次输入到4800中，就可以计算线路任意中、边桩坐标。