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

程序源清单清单中子程序“4”和子程序“7”中的数值为示例高速公路全段的曲线元素(直曲表及逐桩坐标表由“海地2000公路优化设计系统”生成,逐桩坐标表用于检验程序输入的正确性,详见后附录一)ZBJSL1Defm35:Z[26]=0:Z[27]=1:Prog"1"ZBFSL1Defm34:Fixm:{XY}:Z[33]=X"Xp":Z[34]=Y"Yp":Z[26]=2π:Z[27]=1:Prog"1"1L1N=Z[16]:V=Z[17]:Z[26]=2π=>GotoJΔ{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:Lb1J:X=0:{X}:X"QXH":1L2Lb1Z:Z[20]=Ans:Z[19]=X:X≠0=>X=-1=>Lb1K:X=P:Y=Q:S=Z[18]+90:Z[15]=5:Prog"2":GotoM:≠>Prog"E":GotoDΔΔ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:1L4Lb10:{GF}:Z=Z+F"a0":X=X+GcosZ:Y=Y+GsinZ:G≠0=>Goto0ΔProg"D":Lb1D:Z[20]=π=>GotoWΔ1L5Lb19:Z[26]=2π:=>S=10^-6:GotoIΔ{D}:D"Lp"=-1=>Prog"3":Goto9Δ1L6Lb1L:Z[19]=0=>Goto8ΔL7Lb1W:X=Z[19]:Prog"7":π=Ans=>Z[27]=-π=>Z[26]=2π:Z[27]=1ΔGotoZΔX=0=>GotoZΔ1L8Lb18:Z[27]=-π=>"Lp=":D◢"Dp=":K◢Δ1L9S=D-C+T:Lb1I: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ΔΔΔΔΔΔΔ1L10S=W+90+Ans:X=A-TcosW+Rec(Pol(X,YN÷Z,J+W:Y=J+B-TsinW:Z[27]=-π=>GotoXΔZ[26]=2π=>Prog"DX"ΔZ[27]=-5=>Z[27]=-π:GotoLΔZ[27]≠1=>GotoIΔ{K}:K"Dp"=999=>Prog"2":GotoCΔ1L11Lb1X:"Xp=":X=X+KcosS◢"Yp=":Y=Y+KsinS◢L12Z[27]=-π=>GotoMΔLb1C:Pol(X-P,Y-Q:I=Ans:J<0=>J=J+360Δ"Fmp=":J◢L13J=J-Z[18]:J>360=>J=J-360ΔJ<0=>J=J+360Δ"Bp=":J◢"Sp=":I◢L14Lb1M:X=0:Y=1:{X}:X"DHDC(SJGCJS)YORN[N]":X=1=>Prog"H"ΔJ=0:{JXIO}:J"Hg"=0=>GotoA:≠>Z[11]=Z[11]+1ΔZ[11]≠1=>GotoEΔ{O}:Z[12]=O"BM":{O}:Z[13]=O"Hh":Lb1E:"Hc=":J=Z[12]-Z[13]+J◢1L15Z[15]=5=>GotoKΔO=0:{O}: