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}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△ {U}：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=〃◢ （结果显示该桩号方位角） {BO}：B〃S〃O〃⊿〃 （输入边桩距离，交角） P=P+O
L〃XB〃=D+BcosP◢ （结果显示边桩X坐标） M〃YB〃=F+BsinP◢ （结果显示边桩Y坐标）
以上是坐标计算程序，括号内是程序计算的大致原理及说明，中间部分为直线、圆曲线、缓和曲线计算的各种公式，大家也知道，书上也有。
该程序是最基本的，如再加几条语句就可以处理断链。
该程序还需变量输入，把整条线路曲线要素一次输入到4800中，就可以计算线路任意中、边桩坐标。