CASIO4850P 坐标正算、坐标反算程序.doc

程序源清单
清单中子程序“4”和子程序“7”中的数值为示例高速公路全段的曲线元素(直曲表及逐桩坐标表由“海地2000公路优化设计系统”生成,逐桩坐标表用于检验程序输入的正确性,详见后附录一)
ZBJS L1 Defm 35:Z[26]=0:Z[27]=1:Prog "1"
ZBFS L1 Defm 34:Fixm:{XY}:Z[33]=X"Xp":Z[34]=Y"Yp":Z[26]=2π:Z[27]=1:Prog "1"
1 L1 N=Z[16]: V=Z[17]:Z[26]=2π=>Goto JΔ{NV}:Z[15]=0:Z[11]=0:Z[16]=N"X0":Z[17]=V"Y0":P"Xc"=0=>K=Z[16]:L=Z[17]:{KL}: D"S1"E"S2"K"XH1"L"YH1"X"XH2"Y"YH2":Z[16]=K:Z[17]=L:T=Pol(X-K,Y-L:M=cos-1((T2+E2-D2)÷2ET:Fixm:U=J+180-M:P=X+EcosU:Q=Y+EsinUΔQ"Yc":Pol(N-P,V-Q:Fixm:Z[18]=J:Lb1 J:X=0:{X}: X"QXH":
1L2 Lb1 Z:Z[20]=Ans:Z[19]=X:X≠0=>X=-1=>Lb1K:X=P:Y=Q:S=Z[18]+90:Z[15]=5:Prog"2":Goto M:≠>Prog"E":Goto DΔΔ
1L3 {ABCRHWM}:W"Fm"C"JD"A"XJ"B"YJ":Z[10]=R:{R}:Z[8]=R"R2":R= Z[10]:Z[1]=H"LS1":{H}:Z[2]=H"LS2":R≠Z[8]=>M"LF"ΔH=Z[1]:Z=0:X=0:Y=0:
1L4 Lb1 0:{GF}:Z=Z+F"a0":X=X+GcosZ:Y=Y+GsinZ:G≠0=>Goto 0ΔProg"D":Lb1 D:Z[20]=π=>Goto WΔ
1L5 Lb1 9:Z[26]=2π:=>S=10^-6:Goto IΔ{D}:D"Lp"=-1=>Prog "3":Goto 9Δ
1L6 Lb1 L:Z[19]=0=>Goto 8Δ
L7 Lb1 W:X=Z[19]:Prog "7":π=Ans=>Z[27]=-π=>Z[26]=2π:Z[27]=1ΔGoto ZΔX=0=>Goto ZΔ
1L8 Lb1 8:Z[27]=-π=>"Lp=":D◢"Dp=":K◢Δ
1L9 S=D-C+T:Lb1 I:S>L=>X=T+ Z[23]+Rec(S-L+Z[14]+Z[3],N:Y=J:Z:≠>S<0=>X=S:Y=0:≠>R= Z[10]:H=Z[1]:Prog"C":H=0=>Prog "6":≠>S≤H=>Prog "5":≠>SS=S-H:Prog "6":≠>SZ[4]=M÷Abs(Z[8]-Z[10]:Z[10]H=Z[2]:R=Z[8]:O=Z[22]:S=L-S-H-O:I=Z[10]:Prog "A":U=Z-U+180:Prog "B":X=X+Rec(Pol(F,GN÷Z,(U-J)N÷Z:Y=Y+J:Z-O:≠>H=Z[1]:R=Z[10]:O=Z[21]:S=S-H-O:I= Z[8]:Prog "A":X=X+Rec(Pol(F,GN÷Z,(U+J)N÷Z:Y=Y+J:OΔ≠>R=Z[8]:H=Z[2]:Prog "C":SS=L-S-H:Prog "6":Prog
"B":Z-U:≠>S=L-S:Prog "5":Prog "B":Z-IΔΔΔΔΔΔΔ
1L10 S=W+90+Ans:X=A-TcosW+Rec(Pol(X,YN÷Z,J+W:Y=J+B-TsinW: Z[27]=-π=>Goto XΔZ[26]=2π=>Prog "DX"ΔZ[27]=-5=>Z[27]=-π:Goto LΔZ[27]≠1=>Goto IΔ{K}:K"Dp"=999=>Prog "2":Goto CΔ
1L11 Lb1 X: "Xp=":X=X+KcosS◢"Yp=":Y=Y+KsinS◢
L12 Z[27]= -π=>Goto MΔLb1 C: Pol(X-P,Y-Q:I=Ans:J<0=>J=J+360Δ"Fmp=":J◢
L13 J=J-Z[18]:J>360=>J=J-360ΔJ<0=>J=J+360Δ"Bp=":J◢"Sp=":I◢
L14 Lb1 M:X=0:Y=1:{X}:X"DHDC(SJGCJS) Y OR N [N]":X=1=> Prog "H"ΔJ=0:{JXIO}:J"Hg"=0=>Goto