卡西欧5800缓和曲线计算程序

Lbl 0:{EG}:A “ X0 ” ：B “ Y0 ” ：C “ C0 ” ：D “ 1/R1 ” ：E “ 1/R2 ” ：F “ DK1 ” ：G “ DK2 ”Lbl 1:{HOR}:H “ DKi”:O “ DL”:R”DR”:H>G ＝＞Goto 2 ◣

P= （E-D ）/Abs(G-F):Q= Abs(H-F):I=P*Q

J=C+(I+2*D)*Q*90/ π◢

M=C+(I/4+2*D)*Q*45/2 π :N=C+(3*I/4+2*D)*Q*135/2/ π :K=C+(I/2+2*D)*Q*45/ π

X=A+Q*(cos C+4*( cos M + cos N)+2* cos K+ cos J)/12 ◢

Y=B+Q*(sinC+4*(sin M + sin N)+2* sin K+ sin J)/12 ◢

Goto 1 ◣

说明：A 为缓和曲线起点处的纵坐标（X 坐标）；

B 为缓和曲线起点处的横坐标（Y 坐标）；

C 为缓和曲线起点处的切线方位角；

D 为缓和曲线起点处的曲率半径；

E 为缓和曲线终点处的曲率半径；

F 为缓和曲线起点处的里程；

G 为缓和曲线终点处的里程；

H 为缓和曲线中所要放样点处的里程；

注意事项：D 和E 值分别为该缓和曲线前后两段曲线元半径的倒数（即1 / 半径），特别需要强调的是，当曲线右转的时候曲率半径输入时为正值，左转的时候输入时为负值，ZH 点的半径输“ 0 ”。