可用曲線任意里程中邊樁坐標正反算

1、曲線任意里程中邊樁坐標正反算 (CASIO fx-4800P計算器 程序一、程序功能本程序由一個主程序 (TYQXjs和兩個子程正算子程序 (SUB1、反算子程序 (SUB2 序構(gòu)成,可以根據(jù)曲線段直線、圓曲線、緩和曲線(完整或非完整型的線元要素(起點坐標、起點里程、起點切線方位角、線元長度、起點曲率半徑、止點曲率半徑及里程邊距或坐標,對該曲線段范圍內(nèi)任意里程中邊樁坐標進行正反算。另外也可以將本程序中核心算法部分的兩個子程序移植到其它相關(guān)的程序中,用于對曲線任意里程中邊樁坐標進行正反算。 本程序也可以在 CASIO fx-4500P計算器及 CASIO fx-4850P計算 器上運行。二、源程

2、序1. 主程序 (TYQXjs"1.SZ => XY":"2.XY => SZ":N :U"X0" :V"Y0" :O"S0" :G"F0" :H"LS" :P"R0" :R"RN" :Q :C=1÷P :D=(P-R÷(2HPR:E=180÷:N=1=>Goto 1: >Goto 2Lbl 1:SZ:SZ :W=Abs(S-O:Prog "SUB1&q

3、uot;:X"XS"=XY"YS"=YGoto 1Lbl 2:XY:XY :I=X:J=Y:Prog "SUB2":S"S"=O+WZ"Z"=ZGoto 22. 正算子程序 (SUB1A=0.1739274226:B=0.3260725774:K=0.0694318442:L=0.3300094782:F=1-L:M=1-K:X=U+W(Acos(G+QEKW(C+KWD+Bcos(G+QELW(C+LWD+Bcos(G+QEFW(C+FWD+Acos(G+QEMW(C+MWD:Y=V+W(As

4、in(G+QEKW(C+KWD+Bsin(G+QELW(C+LWD+Bsin(G+QEFW(C+FWD+Asin(G+QEMW(C+MWD:F=G+QEW(C+WD+90:X=X+ZcosF:Y=Y+ZsinF3. 反算子程序 (SUB2T=G-90:W=Abs(Y-VcosT-(X-UsinT:Z=0:Lbl 0:Prog "SUB1":L=T+QEW(C+WD :Z=(J-YcosL-(I-XsinL:AbsZ<1E(E即 EXP-6=>Goto1: >W=W+Z:Goto 0Lbl 1:Z=0:Prog "SUB1":Z=(J-Y

5、÷sinF三、使用說明1、規(guī)定(1 以道路中線的前進方向(即里程增大的方向區(qū)分左右;當線元往左偏時,Q=-1;當線元往右偏時, Q=1;當線元為直線時, Q=0。(2 當所求點位于中線時, Z=0;當位于中線左鍘時, Z 取負值;當位于中線中線右側(cè)時, Z 取正值。(3 當線元為直線時,其起點、止點的曲率半徑為無窮大,以 10的 45次代替。(4 當線元為圓曲線時,無論其起點、止點與什么線元相接,其曲率半徑均等于圓弧的半徑。(5 當線元為完整緩和曲線時,起點與直線相接時,曲率半徑為無窮大,以 10的 45次代替;與圓曲線相接時,曲率半徑等于圓曲線的半徑。止點與直線相接時,曲率半徑為無

6、窮大,以 10的 45次代替;與圓曲線相接時,曲率半徑等于圓曲線的半徑。(6 當線元為非完整緩和曲線時,起點與直線相接時,曲率半徑等于設(shè)計規(guī)定的值;與圓曲線相接時,曲率半徑等于圓曲線的半徑。止點與直線相接時,曲率半徑等于設(shè)計規(guī)定的值;與圓曲線相接時,曲率半徑等于圓曲線的半徑。2、輸入與顯示說明輸入部分:1. SZ => XY2. XY = > SZN ? 選擇計算方式,輸入 1表示進行由里程、邊距計算坐標 ;輸入 2表示由坐標反算 里程和邊距。X0 ?線元起點的 X 坐標Y0 ?線元起點的 Y 坐標S0 ?線元起點里程F0 ?線元起點切線方位角LS ?線元長度R0 ?線元起點曲率半

7、徑RN ?線元止點曲率半徑Q ? 線 元左右偏標志 (左偏 Q=-1,右偏 Q=1,直線段 Q=0S ? 正算時所求點的里程Z ?正算時所求點距中線的邊距 (左側(cè)取負,值右側(cè)取正值,在中線上取零 X ?反算時所求點的 X 坐標Y ?反算時所求點的 Y 坐標顯示部分:XS=×××正算時,計算得出的所求點的 X 坐標YS=×××正算時,計算得出的所求點的 Y 坐標S=×××反算時,計算得出的所求點的里程Z=×××反算時,計算得出的所求點的邊距四、算例某匝道的由五段線元(直線 +完

8、整緩和曲線 +圓曲線 +非完整緩和曲線 +直線組成,各段線元的要素(起點里程 S0、起點坐標 X0 Y0、起點切線方位角 F0、線元長度LS 、起點曲率半徑 R0、止點曲率半徑 RN 、線 元左右偏標志 Q 如下:S0 X0 Y0 F0 LS R0 RN Q 500.000 19942.837 28343.561 125 16 31.00 269.256 1E(E即 EXP45 1E45 0 769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1 806.748 19766.566 28594.574 120 25 54.0

9、7 112.779 221.75 221.75 -1 919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.228 -1 999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0 1、正算(注意:略去計算方式及線元要素輸入,請自行根據(jù)所求點所在的線元輸入線元要素S=700 Z=-5 計算得 XS=19831.41785 YS=28509.72590S=700 Z=0 計算得 XS=19827.33592 YS=28506.83837S=700 Z= 5 計算得 XS=1

10、9823.25398 YS=28503.95084S=780 Z=-5 計算得 XS=19785.25749 YS=28575.02270S=780 Z=0 計算得 XS=19781.15561 YS=28572.16358S=780 Z= 5 計算得 XS=19777.05373 YS=28569.30446S=870 Z=-5 計算得 XS=19747.53609 YS=28654.13091S=870 Z=0 計算得 XS=19742.68648 YS=28652.91379S=870 Z= 5 計算得 XS=19737.83688 YS=28651.69668S=940 Z=-5.12

11、3 計算得 XS=19741. 59118 YS=28722.05802S=940 Z=0 計算得 XS=19736.47687 YS=28722.35642S=940 Z= 3.009 計算得 XS=19733.47298 YS=28722.531682、 反算X=19831.418 Y=28509.726 計算得 S=699.9999974 Z= -5 .00018164 X=19827.336 Y=28506.838 計算得 S=699.9996493 Z= 0.000145136 X=19823.25398 Y=28503.95084 計算得 S=699.9999985 Z= 5.00

12、0003137 X=19785.25749 Y=28575.02270 計算得 S=780.0000035 Z= -5 .000001663 X=19781.15561 Y=28572.16358 計算得 S=780.0000025 Z=- 0.000002979 X=19777.05373 Y=28569.30446 計算得 S=780.0000016 Z= 4.99999578X=19747.536 Y=28654.131 計算得 S=870.0001137 Z= -4.99941049X=19742.686 Y=28652.914 計算得 S=870.0003175 Z=- 0.00041814X=19737.837 Y=28651.697 計算得 S=870.0002748