高速公路测量计算CASIO程序全套
坐标主程序
XY(坐标程序)
LbI 2
｛K｝
PROG “SUB3”
X=X◢
Y=Y◢
｛v｝
PROG “SUB 2”：⊿èv=√2
GOTO 2
FANS
｛MNWT｝
M“X0=” ： N“Y0=”
W“X-=” ： T“Y-=”
Z=tan-1((T-N)/(W-M+0.19)
LbI 7
｛UV｝
U“JIAO=”
V“JU=”
U=U+Z：GOTO 3：⊿èW＞M
U=U+Z-180
LbI 3
P=M+VcosU◢
Q=N+VsinU◢
｛K｝
K“DGLC=”
LbI 6
Prog “SUB 3”
O=（Q-Y+（P-X）/（tanE+0.19））×sinE
K=K+O
O=Abs O
GOTOèO≥0.17 6：⊿
K=K◢
V=（Q-Y+（X-P）×tanE）×cosE ◢
GOTO 7
子程序
ZHIX 1
H= K-80100 （80+100为起点里程，直线段）
A=3026441.522 （A、B为K80+100坐标）
B=515965.236
E=270.265 为K80+100方位角
X=A+HcosE
Y=B+HsinE
QUX 1 （K80+100~ K80+948.114曲线右偏）
R=900：L=170
H=K-80310.927 （80+310.927为ZH点里程）
I=K-80480.927 （80+480.927为HY点里程）
A=3026442.499èI≤0
B=515754.311 （A、B为ZH点坐标）
PROG “SUB 1”
E=270.265+3×G （270.265为ZH点方位角）
X=A+Fcos（E-2G）
Y=B+Fsin（E-2G）
Goto 2 ⊿
J=K-80778.114
A=3026448.633 ：B=515584.489èJ≤0 （A、B为HY点坐标）
C=I×90÷π÷R
D=2×R×sinC
E=275.6763+2C （275.6763为HY点方位角）
X=A+D×cos（E-C）
Y=B+D×sin（E-C）
Goto 2 ⊿
S=K-80948.114
A=3026606.193 ：B=515149.157èS≤0 （A、B为HZ点坐标）
H=-S
PROG “SUB 1”
E=300.0072-3×G （300.0072为HZ点方位角）
X=A+Fcos（E+180+2G）
Y=B+Fsin（E+180+2G）
Goto 2 ⊿
LbI 2
备注：曲线左偏与右偏，只有中间几个“＋－”号的区别。以下不再备注）
QUX 2 （K80+948.114~ K81+475.631曲线左偏）
R=1000：L=174.535
H=K-80948.114
I=K-81122.649
A=3026606.193èI≤0
B=515149.157
PROG “SUB 1”
E=300.0072-3×G
X=A+Fcos（E+2G）
Y=B+Fsin（E+2G）
Goto 2 ⊿
J=K-81301.097
A=3026689.019èJ≤0 ：B=514995.593
C=I×90÷π÷R
D=2×R×sinC
E=295.007-2C
X=A+D×cos（E+C）
Y=B+D×sin（E+C）
Goto 2 ⊿
S=K-81475.631
A=3026784.297 ：B=514657.014èS≤0
H=-S
PROG “SUB 1”
E=279.783+3×G
X=A+Fcos（E+180-2G）
Y=B+Fsin（E+180-2G）
Goto 2 ⊿
LbI 2
ZHIX 2
H= K-81475.632
A=3026784.297
B=514657.014
E=279.7828
X=A+HcosE
Y=B+HsinE
QUX 3 （K82+420.711~ K83+354.430曲线左偏）
R=1200：L=180
H=K-82420.711
I=K-82600.711
A=3026944.880èI≤0
B=513725.679
PROG “SUB 1”
E=279.7828-3×G
X=A+Fcos（E+2G）
Y=B+Fsin（E+2G）
Goto 2 ⊿
J=K-83174.43
A=3026971.015èJ≤0 ：B=513547.632
C=I×90÷π÷R
D=2×R×sinC
E=275.4856-2C
X=A+D×cos（E+C）
Y=B+D×sin（E+C）
Goto 2 ⊿
S=K-83354.43
A=3026814.4519 ：B=512821.7913èS≤0
H=-S
PROG “SUB 1”
E=243.7953+3×G
X=A+Fcos（E+180-2G）
Y=B+Fsin（E+180-2G）
Goto 2 ⊿
LbI 2
SUB 1
C=H-H^5÷40÷R2÷L2
D= H^3÷6÷R÷L- H^7÷336÷R^3÷L^3
F=√(C2+D2)
G=H2×30÷π÷R÷L
SUB 2
｛U“L=”｝
M“XL”=X+Ucos（E-90°）◢
N“YL”=Y+Usin（E-90°）◢
｛V“R=”｝
P“XR”=X+Vcos（E+90°）◢
Q“YR”=Y+Vsin（E+90°）◢
SUB 2（也可改为如下，但在输入外矢距时注意，线左为负数，线右为正数）
｛U｝
U“[L-，R+]=”
P“EX”=X+Ucos（E+90°）◢
Q“EY”=Y+Usin（E+90°）◢
SUB3
PROGèK≤80310.927 “ZX 1”：GOTO 1：⊿
PROG “QX 1”： GOTO 1：⊿èK≤80948.114
PROG “QX 2”：èK≤81475.632 GOTO 1：⊿
PROG “ZX 2”：GOTO 1：⊿èK≤82420.711
PROG “QX 3”： GOTOèK≤83354.43 1：⊿
LbI 1
根据坐标反算线路桩号及左右偏值
File ZBFS
1: M“X0” : N“Y0” : A“PJ” : R : L“LS” : T : O“L” : H“ZH” : F“FWJ”
2: Z[1]=P“XJD”+T cos ( F+180) : Z[2]=K“YJD”+T sin ( F+180)
3: U“R+1,L-1”
4: Z[3]=H+L : Z[4]=H+O－L : Z[5]=H+O
5: Lbl 0 : C“LC” : Z[7]=Abs(C－H)
6: Prog“KL”
7: S=－((X－M)sin(Q+90)－(Y－N) cos (Q+90) : Abs S<0.0001 C“ LC ” ◢
W=“LP ”= ( X－M ) sin Q－(Y－N) cos Q◢ C = C+S : Goto 0 △
File KL
1: C≤H Z[8]=－Z[7] : Z[9] = 0 : Z[10] = 0 : Goto 2 △
2: C≤Z[3] Z[8]= Z[7]－Z[7] ^5÷40÷R2÷L2 : Z[9] = Z[7]^3÷6÷R÷L－Z[7] ^7÷336÷R^3÷L^3 : Z[10] = 90 Z[7] 2÷R÷L÷Л : Goto 2 △
3: C≤Z[4] Z[10] = 90 L÷R÷Л +180( Z[7] －L)÷R÷Л : Rsin (Z[10]+0.5L－L^3÷240÷R2 : Z[9] = R(1－cosZ[10] )+L2÷24÷R : Goto 2 △
4: C≥Z[5] X=－( Z[7] －O) : Y = 0 : Z[10] = A : Goto 1 △
5: C≥Z[4] Z[10]= A－90 (O－ Z[7] ) 2 ÷R÷L÷Л :X = (O－Z[7] )－(O－Z[7] )^5÷40÷R2÷L2 : Y = (O－Z[7] )^3÷6÷R÷L－(O－Z[7] )^7÷336÷R^3÷L^3 △
6: Lbl 1 : Z[8] = T+(T－X) cos A－Y sin A : Z[9] =(T－X) sin A+ Y cos A : Goto 2 △
7: Lbl 2 : Pol (Z[8],Z[9]) : Z[7] =F+JU : Q=F+Z[10] U : X=Z[1] + IcosZ[7] : Y=Z[2] + I sin Z[7]
注：
M、N：已知坐标
A：偏角
R：半径
L：缓和曲线长
T：切线长
O：曲线长
H：直缓点桩号
F：第一切线方位角
U：路线前进方向左“-”右“+”
高速公路的一些线路计算
一、缓和曲线上的点坐标计算
已知：①缓和曲线上任一点离ZH点的长度：l
②圆曲线的半径：R
③缓和曲线的长度：l0
④转向角系数：K(1或－1)
⑤过ZH点的切线方位角：α
⑥点ZH的坐标：xZ，yZ
计算过程：
说明：当曲线为左转向时，K=1，为右转向时，K=-1，
公式中n的取值如下：
当计算第二缓和曲线上的点坐标时，则：
l为到点HZ的长度
α为过点HZ的切线方位角再加上180°
K值与计算第一缓和曲线时相反
xZ，yZ为点HZ的坐标
切线角计算公式：
二、圆曲线上的点坐标计算
已知：①圆曲线上任一点离ZH点的长度：l
②圆曲线的半径：R
③缓和曲线的长度：l0
④转向角系数：K(1或－1)
⑤过ZH点的切线方位角：α
⑥点ZH的坐标：xZ，yZ
计算过程：
说明：当曲线为左转向时，K=1，为右转向时，K=-1，
公式中n的取值如下：
当只知道HZ点的坐标时，则：
l为到点HZ的长度
α为过点HZ的切线方位角再加上180°
K值与知道ZH点坐标时相反
xZ，yZ为点HZ的坐标
三、曲线要素计算公式
公式中各符号说明：
l——任意点到起点的曲线长度（或缓曲上任意点到缓曲起点的长度）
l1——第一缓和曲线长度
l2——第二缓和曲线长度
l0——对应的缓和曲线长度
R——圆曲线半径
R1——曲线起点处的半径
R2——曲线终点处的半径
P1——曲线起点处的曲率
P2——曲线终点处的曲率
α——曲线转角值
四、竖曲线上高程计算
已知：①第一坡度：i1(上坡为“＋”，下坡为“－”)
②第二坡度：i2(上坡为“＋”，下坡为“－”)
③变坡点桩号：SZ
④变坡点高程：HZ
⑤竖曲线的切线长度：T
⑥待求点桩号：S
计算过程：
五、超高缓和过渡段的横坡计算
已知：如图，
第一横坡：i1
第二横坡：i2
过渡段长度：L
待求处离第二横坡点（过渡段终点）的距离：x
求：待求处的横坡：i
解：d=x/L
i=(i2-i1)(1-3d2+2d3)+i1
六、匝道坐标计算
已知：①待求点桩号：K
②曲线起点桩号：K0
③曲线终点桩号：K1
④曲线起点坐标：x0，y0
⑤曲线起点切线方位角：α0
⑥曲线起点处曲率：P0(左转为“－”，右转为“＋”)
⑦曲线终点处曲率：P1(左转为“－”，右转为“＋”)
求：①线路匝道上点的坐标：x,y
②待求点的切线方位角：αT
计算过程：
CASIO fx-4800P )
F1 XLCS （主程序）
L1 Norm:Deg:U=O″A0″:Prog ″1″:Q=U:C″X-JD″:D″Y-JD″:U=A″A:R+L-″:
Prog ″1″:B=Abs U:R:S″L0″:E″K-ZH″
L2 Fix 3:M=.5S-S＾3÷240R2:P=S2÷24R:T″T″=(R+P)tan.5B+M◢F″L″=πRB÷
180+S◢F=F+E:I=0:J=0:Norm
L3 K″ZJ:XY＝＞1″:K≠1＝＞I=L″K″:U=0:Prog ″4″: G=X:H=Y:≠＞G″X″:H″Y″◣{K}:
K″HS:XY＝＞1″:K=1＝＞X=V″X″:Y=W″Y″:≠＞I=W″K″:U=0:Prog ″4″◣Fix 3
L4 Pol(X-G,Y-H:I″D0″=I◢Fix 4:N=J:J<0＝＞J=J+360◣Prog ″2″:J″A0″=J◢Norm:
U=0:{U}:U″AB″:Prog ″1″:K=U
L5 Lb1 0:U=0:{U}:U″CS: XY＝＞1″:U=0＝＞{UZ}:I=Z″K″:U″BZ:R+L-″:A<0＝＞U
=-U◣Prog ″4″: ≠＞{XY}:X:Y◣Prog ″5″: Prog ″6″:Goto 0
F2 1 （十进制→六十进制子程序）
L1 U=Int U+Frac U÷.6+Frac 100U÷90
F3 2 （六十进制→十进制子程序）
L1 60Frac J:J=Int J+.01Int Ans+.006Frac Ans
F4 3 （缓和曲线上任意点坐标计算子程序）
L1 Y=RS:X=I-I＾5÷40Y2:Y=I＾3(1-I＾4÷56Y2)÷6Y+URec(1,90I2÷πY:X=X-UJ
F5 4 （中桩、边桩坐标计算子程序）
L1 I＞F＝＞X=T-UsinB+Rec(I-F+T,B:Y=J+UcosB:Goto 0◣
L2 I＞F-S＝＞I=F-I:Prog ″3″: Rec(1,B:U=X:X=T+TI-XI-YJ:Y=TJ-UJ+YI:Goto 0◣
L3 I＞E+S＝＞Y=P+R-Rec(R-U,180(I-E-.5S)÷πR:X=M+J:Goto 0◣
L4 I＞E＝＞I=I-E:Prog ″3″: ≠＞X=I-E:Y=U◣
L5 Lbl 0: A<0＝＞Y=-Y◣Rec(1,Q:U=X:X=C-TI+XI-YJ:Y=D-TJ+UJ+YI
F6 5 （测设数据输出子程序）
L1 Pol(X-G,Y-H:J=J-N:J＜0＝＞J=J+360◣J=K+J:J≥360＝＞J= J-360◣
Fix 4:Prog ″2″:J″AC″=J◢Fix 3:I″DC″=I◢Norm
F7 6 （测设时移桩数据计算子程序）
L1 Lbl 0:Norm:L=0:J=0
L2 Lbl 1:U=0:{U}:U″SC″:U≠0＝＞J=J+U:L=L+1:Goto 1◣J≠0＝＞{U}:U″V″:
Prog ″1″:Fix 3:J″SD″=J÷L×sin U◢J″MOVE″=I-J◢Goto 0◣