《CASIO fx-4850P线路计算9.0版·匝道版》程序设计 作者：徐光辉 2004年9月26日
10．匝道高程计算子程序CG
V=Z[43]:N=Z[66+Z[1]]:H=Z[102+Z[1]]:W=H+2? Lbl 0?
W>H+N+1?W=H+N+1:Goto 1? P=Z[W+3+2N]?
K>Z[W]-P?W=W+1:Goto 0? Lbl 1?
V=W-1?Goto 2? V=W-1:Z[43]=V? J=Z[V+2+N]?
R=(J-Z[V+1+N])÷(Z[V]-Z[V-1])? S=(Z[V+3+N]-J)÷(Z[V+1]-Z[V])? T=Z[V+3+2N]? Lbl 2? L=K-Z[V]?
L≤-T?Y=J+RL:Goto 3? L≥T?Y=J+SL:Goto 3? A=2T÷(S-R)? X=L+A(R+S)÷2? Y=J+XX÷2÷A-ARS÷2? Lbl 3? Y=Y-Z[1+Z[1]]? Q=0?C=0:Goto 7? A=Abs Q?
A≤Z[46]?C=0:Goto 7?
I=-1?
Q>0?I=1?
Z[38]=Z[46+I]+Z[46]?
N=Z[72+Z[1]]:H=Z[108+Z[1]]:W=H+1?
21
《CASIO fx-4850P线路计算9.0版·匝道版》程序设计 作者：徐光辉 2004年9月26日 Lbl 5?
W>H+N-1?W=H+N-1:Goto 6? K>Z[W]?W=W+1:Goto 5? Lbl 6?
@=Z[W-1+N]:X=Z[W+N]:@=I@:X=IX? Z[W-1+N]=0?@=2? Z[W+N]=0?X=2?
Z=(Z[W]-K)÷(Z[W]-Z[W-1])?
M=(X-@)(1-3ZZ+2ZZZ)+@? （当超高过渡段按直线性计算时，此行换为 M=(@-X)Z+X?） A≤Z[38]?C=(Z[46]-A)M÷100:Goto 7? C=Z[31](Z[38]-A)÷100-M(Z[38]-Z[46])÷100? M>Z[31]?C=(Z[46]-A)M÷100? Lbl 7? C=Y+C
11．反求匝道桩号子程序FZ Z[1]=Int Z[1]? Z[1]>6?Z[1]=6? Z[1]<1?Z[1]=1?
N=Z[60+Z[1]]:H=Z[96+Z[1]]? Z[41]=X:Z[42]=Y? X=Z[H-3]? Y=Z[H-2]? U=Z[H-1]? W=1:@=0
Z[41]=X?Z[42]=Y?@=100@+W?
Lbl 0?
W>N-1?Goto 1?
Z[31]=X:Z[32]=Y:Z[33]=U?
Z[34]=Z[H+W]:Z[35]=Z[H+N+W-1]:Z[36]=Z[H+N+W]:Z[37]=Z[H+W-1]?
22
《CASIO fx-4850P线路计算9.0版·匝道版》程序设计 作者：徐光辉 2004年9月26日 Z=Z[34]:Prog “ZD”? Z[38]=X:Z[39]=Y:Z[40]=U?
Z[31]=Z[41]?Z[32]=Z[42]?C=π:Goto 6
?
Pol(Z[31]-Z[41],Z[32]-Z[42]):C=J+π-Z[33]? Lbl 6?
Z[38]=Z[41]?Z[39]=Z[42]?D=π:Goto 7
?
Pol(Z[38]-Z[41],Z[39]-Z[42]):D=J+π-Z[40]? Lbl 7?
Cos C>0?cos D≤0?@=100@+W?
W=W+1? Goto 0? Lbl 1?
@=0?Q=0:”NO”:K=0◢ Goto 5?
@>100?{@}:@”SEL QXD:”? W=1:X=Z[H-3]:Y=Z[H-2]:U=Z[H-1]? Lbl 2?
Z[31]=X:Z[32]=Y:Z[33]=U?
Z[34]=Z[H+W]:Z[35]=Z[H+N+W-1]:Z[36]=Z[H+N+W]:Z[37]=Z[H+W-1]? W≥@?Goto 3?
Z=Z[34]:Prog “ZD”:W=W+1? Goto 2? LbL 3?
X=Z[41]:Y=Z[42]?
X=Z[31]?Y=Z[32]?K=Z[37]:Q=0:Goto 4
?
Q=Z[31]:K=Z[32]:G=Z[33]-π÷2:Prog “ZH”? Z[35]≠Z[36]?Prog “HZ”:Goto 4? Z[35]=0?K=Z[37]+W:Q=-H:Goto 4? R=1÷Z[35]?
Pol(H+R,W):K=Z[37]+RJ?
23
《CASIO fx-4850P线路计算9.0版·匝道版》程序设计 作者：徐光辉 2004年9月26日 Z[35]>0?Q=R-I:≠?Q=R+I:K=K-Rπ?
Lbl 4?
“KD=”:Int (1000K+.5)÷1000◢ “WB=”:Int (1000Q+.5)÷1000◢ Lbl 5
12．坐标转换计算子程序ZH H=X-Q:W=Y-K? Pol(H,W)? G=J-G?
H=Icos G:W=Isin G
13．反求缓曲长度子程序HZ R=Z[35]：P=Z[36]:S=Z[34]-Z[37]? B=RS÷(R-P)? P>R?Z=1:≠?Z=-1? D=-ZSB?
C=Abs (S÷(P-R))? C=0?C=1? A=BD÷2÷C? T=Z[33]-ZDD÷2÷C?
I=-D+D^5÷40÷C÷C-D^9÷3456÷C^4+D^13÷599040÷C^6?
J=-D^3÷6÷C+D^7÷336÷C^3-D^11÷42240÷C^5+D^15÷9676800÷C^7? Z[38]=Z[31]+ZIcos T-Jsin T? Z[39]=Z[32]+ZIsin T+Jcos T? A=Z[33]+A? B>0?A=A+π
? A>2π?A=A-2π? A<0?A=A+2π
?
I=Abs B:J=Abs (PS÷(R-P))?
24
《CASIO fx-4850P线路计算9.0版·匝道版》程序设计 作者：徐光辉 2004年9月26日 R=Abs R:P=Abs P?
R>P?R=1÷R:≠?R=1÷P? I>J?U=I:≠?U=J?
Q=Z[38]:K=Z[39]:G=A-π÷2:Prog “ZH”? Z=-Z:U=RU?
A=Z÷46448640÷U^8? C=Z÷241920÷U^6? D=H÷3840÷U^5? E=Z÷945÷U^4? F=ZW÷384÷U^4? G=H÷48÷U^3? I=Z÷(15UU)? J=ZW÷(8UU)? L=H÷2÷U? M=Z? P=ZW? S=W? Lbl 0?
N=AS^17-CS^13+DS^10-ES^9+FS^8-GS^6+IS^5-JS^4+LSS-MS+P? Q=17AS^16-13CS^12+10DS^9-9ES^8+8FS^7-6GS^5+5IS^4-4JS^3+2LS-M? @=S-N÷Q?
Abs (@-S)<10-6?Goto 1? S=@? Goto 0? Lbl 1?
T=Z@@@÷6÷U-Z@^7÷336÷U^3+Z@^11÷42240÷U^5-Z@^15÷9676800÷U^7? P=@-@^5÷40÷U÷U+@^9÷3456÷U^4-@^13÷599040÷U^6? T=H?P=W?Q=0:Goto 2?
Pol(H-T,W-P)?
Q=Sin (J+Z@@÷2÷U-π÷2)?
25