全線直曲計算程序通用

1、經(jīng)實踐在5800上可行，反算運行時間就是慢點有4秒左右才出結果。5800測量程序加入了數(shù)據(jù)庫，綜合了網(wǎng)上的一些而成。非本人原創(chuàng)，希望能與廣大朋友交流。5800坐標正反程序(加入了數(shù)據(jù)庫) 1. 主程序(TYQXJS)Lbl 4："1.SZ => XY"："2.XY => SZ"：？N：？S：Prog“SUB0”1÷PC: (P-R)÷(2HPR) D:180÷E:N=1 => Goto1： Goto2：Lbl 1:? Z： Abs(S-O) W：Prog "SUB1"： "XS

2、="：X"YS="：YF-90F： “FS=”：FDMSGoto4Lbl 2：？X： ？Y：XI：YJ： Prog“SUB2”：O+WS：“S=”：S“Z=”：ZGoto42. 正算子程序(SUB1)0.1739274226A： 0.3260725774B： 0.0694318442K： 0.3300094782L： 1-LF： 1-KMU+W(Acos(G+QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Bcos(G+QEFW(C+FWD)+Acos(G+QEMW(C+MWD)X： V+W(Asin(G+QEKW(C+KWD)+Bsin(G+QEL

3、W(C+LWD)+Bsin(G+QEFW(C+FWD)+Asin(G+QEMW(C+MWD) Y：G+QEW(C+WD)+90F： X+Zcos（F）X： Y+Zsin（F）Y2. 反算子程序(SUB2)G-90T： Abs(Y-V)cos(T)-(X-U)sin(T) W： 0Z：Lbl 0：Prog "SUB1"：T+QEW(C+WD) L： (J-Y)cos（L）-(I-X)sin（L）Z：IF Abs（Z）<0.000001：Then0Z:Prog "SUB1" :(J-Y)÷sin（F）ZElse W+ZW:Goto4:IfEn

4、dSUB0 數(shù)據(jù)庫子程序Goto1 同時保存多個曲線時的指針Lbl 1：IF S<*（線元終點里程）:Then*G（線元起點方位角）:*O（線元起點里程）:*U（線元起點X）:*V（線元起點Y）:*P（線元起點曲率半徑）:*R（線元終點曲率半徑）: *H（線元起點至終點長度）：0或1、-1Q：Return:IfEndLbl 1：IF S<*（線元終點里程）:Then*G（線元起點方位角）:*O（線元起點里程）:*U（線元起點X）:*V（線元起點Y）:*P（線元起點曲率半徑）:*R（線元終點曲率半徑）: *H（線元起點至終點長度）：0或1、-1Q：Return:IfEnd.為了便于解

5、讀，每增加一個線元增加一行語句，每增加一條曲線增加一個Lbl，每增加一個工程增加一個文件。三、使用說明1、規(guī)定(1) 以道路中線的前進方向（即里程增大的方向）區(qū)分左右；當線元往左偏時，Q=-1；當線元往右偏時，Q=1；當線元為直線時，Q=0。(2) 當所求點位于中線時，Z=0；當位于中線左鍘時，Z取負值；當位于中線中線右側時，Z取正值。(3) 當線元為直線時，其起點、止點的曲率半徑為無窮大，以10的45次代替。(4) 當線元為圓曲線時，無論其起點、止點與什么線元相接，其曲率半徑均等于圓弧的半徑。(5) 當線元為完整緩和曲線時，起點與直線相接時，曲率半徑為無窮大，以10的45次方代替；與圓曲線相

6、接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半徑為無窮大，以10的45次方代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。(6) 當線元為非完整緩和曲線時，起點與直線相接時，曲率半徑等于設計規(guī)定的值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半徑等于設計規(guī)定的值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。2、輸入與顯示說明輸入部分：1. SZ => XY2. XY = > SZN ? 選擇計算方式，輸入1表示進行由里程、邊距計算坐標 ；輸入2表示由坐標反算里程和邊距。S ？正算時所求點的里程，反算時為所求點的近似里程，反算時輸入的近似里程不能夸線元，

7、若夸線元的話，需進行第二次反算，并以第一次的結果做為近似里程，否則反算的里程及邊距不對， Z ？正算時所求點距中線的邊距(左側取負，值右側取正值，在中線上取零)X ？反算時所求點的X坐標Y ？反算時所求點的Y坐標顯示部分：XS=××× 正算時，計算得出的所求點的X坐標YS=××× 正算時，計算得出的所求點的Y坐標FS=××× 正算時，所求點對應的中線點的切線方位角S=××× 反算時，計算得出的所求點的里程Z=××× 反算時，計算得出的所求點的邊距

8、要注意區(qū)別字母O和數(shù)字0，為回車符號。ClrMemory:1N: 5800計算器由坐標計算面積小程序LbI 0: E:F:"XO"?X:"Y0"?Y:(XF-EY)/2G:Q+GQ:XE:YF:Abs(Q)S:Fix 1:"   S=":Locate 6,4,S Norm 2:"   N=":Locate 6,4,N N+1N:Goto 0:(說明：至少輸入三點坐標，才能計算正確的面積，N為了計輸了的點數(shù)）這段時間下雨，閑就沒事，就搞點小程序，還有需要計算器程序的，留言說聲，

9、盡能力幫你編，同時也是跟大家一起學習交流。CASIO fx5800計算器 工程測量與公路測量程序1、坐標正算”X0=?X:”Y0=?Y:”I=?I:”J=?J X+ICos(J)U:Y+ISin(J)V ”X=:U （待求點的X坐標）”Y=:V （待求點的Y坐標）說明：X0 Y0：已知點坐標I：兩點的距離J：方位角2、坐標反算Lbl 0 ”X1=?X:”Y1=?Y:”X2=?U:”Y2=?V Pol(U-X,V-Y): J<0 J+360J ”I=:I ”J=:J DMS GOTO 0說明：X1Y1：第一點的坐標，X2Y2第二點的坐標，I：兩點的距離，J：方位角3、圓曲線”X0=?X:

10、”Y0=?Y:”FWJ=?D:”QDHAO=?G: ”ZDHAO=?C: ”R=?R Lbl 1 ”DQHAO=?L L<G Or LC GOTO 2 ”PY=?K 180(L-G)÷(R) E 2RSin(0.5E) F ReC(Abs(F),D+0.5E)：Cls ”X=:X+I+KCos(D+E+90) A ”Y=:Y+J+KSin(D+E+90) B D+Eo:o<0 360+o0o ”FWJ=:o DMS Goto 1 Lbl 2 ”END說明：X0Y0：起始點的坐標，F(xiàn)WJ：起始點的方位角，QDHAO：起點里程，ZDHAO終點里程，R：半徑，DQHAO：待求點

11、里程，PY：偏移量4、豎曲線”QZHAO =?J: ”H+B=?B:”I1=?C:”I2=?D:”R=?R:”T=?T: (D-C)÷100) ÷Abs(D-C) ÷100)A:Abs(0.5R(D-C) ÷100)S ”T:S J-TX JTY Lbl 0 ”DQ=?L If L<X:Then Goto 1:Else If LY:Then Goto 1:If End:If End ”GC=:B+C(L-J) ÷100+A(L-J+T) ÷(2R) H Goto 0 Lbl 1 ”END說明：QZHAO：曲中點里程，H+B：曲中

12、點高程，I1：第一坡度，I2：第二坡度，R：半徑，T：切線長，DQ：待求點里程。5正算偏移”X0=?A: ”Y0=?B:”FWJ=?D: Lbl 0 ”LS=?L: ”PY=?C ReC(L,D) ”X=:A+I+CCos(D+ 90) X ”Y=:B+J+CSin(D+ 90) Y Goto 0說明：X0 Y0：已知點坐標，LS：兩點距離，F(xiàn)WJ：方位角6.萬能公式”X0=?A:”Y0=?B:”FWJ=?C:”1÷R-QD?D:”1÷R-ZD?E:”QDHAO?F: ”ZDHAO?G Lbl 1 ”DQHAO=?H:”PY=?o If H<F:Then Goto 2

13、:Else If HG:Then Goto 2:If End:IfEnd (E-D)÷Abs(G-F) P:Abs (H-F) Q:PQI:C+90Q(I+2D) ÷J:C+45Q(I÷8+D) ÷M:C+135Q(3I÷8+D) ÷N:C+45Q(I÷2+2D) ÷K A+Q(Cos(C) +4(Cos(M)+Cos(N)+2Cos(K)+Cos(J)÷12U:B+Q(Sin(C)+4(Sin(M)+Sin(N)+2Sin(K)+Sin(J)÷12V ”X=:U+oCos(J+90) X ”Y

14、=:V+oSin(J+90) Y ”FWJ=: J DMS Goto 1 Lbl 2 ”END說明：X0Y0：起開點坐標，F(xiàn)WJ：起開點方位角，1÷R-QD：起開點曲率，1÷R-ZD終點曲率，QDHAO：起點里程，ZDHAO終點里程，DQHAO：待求點里程，XY：待求點坐標6、標高計算”QZHAO=?J:”QDHAO=?G:”ZDHAO=?K:”BHB=?B:”SX=?Q:”I1=?C:”I2=?D:”R=?R:”T=?T (D-C)÷100W:W÷Abs(W)A:”T:0.5RAbs(W) S ”o-YIGAO=?o J-TX:J+TY Lbl 1 ”LDQHAO=?L:”PY=?U:”XIELV=?V If L<X:Then Goto 2:Else If LY:Then Goto 3:If End:I End ”HL=:B+0.01C(L-J)+A(L-J+T) ÷(2R)+0.01UV+QQ ”HO=:o-HZ Goto 1 Lbl 2 L<G Goto 4 ”HL=:B+0.01 (L-J)C+0.01UV+Q H ”H0=:o-HH Goto 1 Lab 3 L>K Got