复化辛卜生公式可以算是公路线形坐标计算的万能公式。它不仅对直线、圆曲线、缓和曲线,不完整的缓和曲线(卵形曲线)都能用一套公式进行计算。当然它的精度可以人为判断,不过它使用起来还是比较繁琐,一不小心更容易出错,对直线、圆曲线的坐标计算更无必要用此公式。但是如果在Casio上编一套程序,对任何线形只要用这一套程序,只要输入起算点坐标、起始方位角、桩号几个参数,使用起来就很方便了。
笔者编程如下:
File1
FHXPS(复化辛卜生)
L1 Lb1 0
L2 {L}
L3 W=(1∕B-1∕A)∕(S-K)
L4 O=1∕A+W×(L-K)
L5 I=N+(O+1∕A)(L-K)×180∕π∕2◢
L6 V=1∕A+W×(L-K)/2
L7 U=N+(V+1∕A)(L-K)×180∕π∕4
L8 C=1∕A+W×(L-K)∕4
L9 D=N+(C+1∕A)(L-K)×180∕π∕8
L10 E=1∕A+W×(L-K)×3∕4
L11 F=N+(E+1∕A)×(L-K)×180×3∕π∕8
L12 R"YPPQX:1,ZPPQX:-1"
L13 R=1 =>Prog“YP”:=> R=-1=> Prog“ZP”⊿⊿
L14 Goto 0
File2
YP
L1 X=G+(L-K)∕6∕2×(cosN+4(cosD+cosF)+2cosU+cosI) ◢
L2 Y=H+(L-K)∕6∕2×(sinN+4(sinD+sinF)+2sinU+sinI) ◢
File3
ZP
L1 X=G+(K-L)∕6∕2×(cosN+4(cosD+cosF)+2cosU+cosI) ◢
L2 Y= H+(K-L)∕6∕2×(sinN+4(sinD+sinF)+2sinU+sinI) ◢
A ---- 起点半径(曲率)
B ---- 终点半径(曲率)
S ---- 终点桩号
K ---- 起点桩号
L ---- 待求点桩号
I ---- 待求点切线方位角
N ---- 起始切线方位角
V ---- 1∕2等分点曲率
U ---- 1∕2等分点切线方位角
C ---- 1∕4等分点曲率
D ---- 1∕4等分点切线方位角
E ---- 3∕4等分点曲率
F ---- 3∕4等分点的切线方位角
G---起算点X坐标
H---起算点Y坐标
当然你在程序里可以加上一个计算边桩的子程序。