匝道坐标计算程序

匝道坐标计算程序（辛甫生公式）
L1 LbI0:｛EG｝:A“XA”:B“YA”:C“CA”:D“1/RA”:E“1/RB”:F“KA”:G“KB”
L2 LbI1:｛HORS｝:H“KI”:O“DL”:R“DR”:S“PJ”:H＞G =＞ Goto2△
L3 P=（E-D）/Abs（G-F）:Q=Abs（H-F）:I=P×Q:T=D+I
L4 J=C+（I+2D）Q×90/π⊿M=C+（I/4+2D）Q×45/(2π):N=C+（3I/4+2D）Q×135/(2π)
L5 K=C+（I/2＋2D）Q×45/π
L6 X=A+Q(cosC+cosJ+4(cosM+cosN)+2cosK)/12⊿
L7 Y=B+Q(sinC+sinJ+4(sinM+sinN)+2sinK)/12⊿
L8 U“XL”=X+Ocos(J-180+S)⊿V“YL”=Y+Osin(J-180+S)⊿
L9 W“XR”=X+Rcos(J+S)⊿Z“YR”=Y+Rsin(J+S)⊿
L10 A=X:B=Y:C=J:D=T:F=H:Goto1
L11 LbI2:A=X:B=Y:D=E:F=G:C=J:Goto0
数据输入说明：
① XA、YA、CA分别为A点（起始点）的X、Y坐标及方位角。

② 1/RA、1/RB分别为A点和B点的曲率，以1/R表示。

曲线左偏时取“－”；曲线右偏时取“＋”；直线时为0；缓和曲线时： 1/RA为0或是〒1/RA、1/RB为0或是〒1/RB；圆曲线时均为〒1/R。

【左“－”右“＋”】
③ KA、KB分别为曲线原起点A和终点B的里程桩号。

④ XL、YL、XR、YR分别为左边桩和右边桩的X、Y坐标。

⑤ DL、DR分别为左、右边桩距中桩的水平距离。

⑥ 当一个曲线计算完后，即KI＞KB时，程序显示1/RB和KB，此时需输入下一个曲线原终点的曲率和桩号。

⑦ 计算无缓和曲线时，不能由直线段直接计算而应将KA、1/RA输入LR点桩号和〒1/R。

⑧ S为偏角，定义为路线前进方向顺时针旋转到的夹角。