任意里程中邊樁坐標(biāo)正反算專用程序

1、曲線任意里程中邊樁坐標(biāo)正反算 (CASIO-4850計(jì)算器 專用程序一、程序功能及原理1. 功能說明:本程序由一個(gè)主程序 (TYQXjs和四個(gè)子程正算子程序 (SUB1、 反算子程序 (SUB2等構(gòu)成, 可以根據(jù)曲線段直線、圓曲線、緩和曲線(完整或非完整型的線元要素(起點(diǎn)坐標(biāo)、起點(diǎn)里程、起 點(diǎn)切線方位角、線元長度、起點(diǎn)曲率半徑、止點(diǎn)曲率半徑及里程邊距或坐標(biāo),對(duì)該曲線段范圍內(nèi)任意里 程中邊樁坐標(biāo)進(jìn)行正反算。本修改版程序既可實(shí)現(xiàn)正算全線貫通,亦可實(shí)現(xiàn)反算全線貫通。本程序也可以 在 CASIO fx-4800P 計(jì)算器運(yùn)行,但由于容量較小,不建議使用。2.計(jì)算原理:利用 Gauss-Legendre

2、 5點(diǎn)通用 (可查閱 <測(cè)繪工程 >2006年第二期,由武漢大學(xué)花向 紅所著 -道路平面測(cè)設(shè)計(jì)算的通用公式 公式正算線路中邊樁坐標(biāo)、 線外測(cè)點(diǎn)至曲線元起點(diǎn)和終點(diǎn)的 垂距的符號(hào)是否相異(即 Dca ×Dcb<0=>該測(cè)點(diǎn)在其線元內(nèi)進(jìn)行判斷并利用該線元要素反算中樁里程、 支距,最后計(jì)算出放樣數(shù)據(jù)。二、源程序1. 主程序 (TYQXjsLbl 0:Defm 400:"1.SZ => XY ,2.XY => SZ" :N :N=1=>Goto 1: >Goto 2Lbl 1:SZ:SZS Z26=>J=1:Prog “

3、 DA T1” :Goto 3S Z34=>J=2:Prog “ DA T1” :Goto 3S Z42=>J=3:Prog “ DA T1” :Goto 3S Z50=>J=4:Prog “ DA T1” :Goto 3S Z8(N+1 +10=>J=n:Prog “ DA T1” :Goto 3Lbl 3:W=Abs(S-O:Prog "SUB1" :"XS=":X=X"YS=":Y=Y” FWI=” : F=F-M :F DMS Goto 4Lbl 2:JXY:JXY :Z4=X:Z5 =Y:N=0Lbl

4、 A:Isz N :N=J: Prog" DA T1": Z11=U:Z12=V: Z9=G:A= Z9-M:J=J+1: Prog "DA T1": Z13=U: Z14=V: Z10=G:B= Z10-M:Prog "ZX1" :Z6×Z7<0=>J=N:Prog DA T1" :Goto B A=Z9+M: B=Z10+M:Prog "ZX1" :Z6×Z7<0=>J=N:Prog "DA T1" :Goto B : > Goto

5、 A Lbl B:Prog "SUB2" :"S=":S=O+W"Z=":Z=ZGoto 2Lbl 4:J=0:I=Pol(X-Z1,Y-Z2:F=J:F<0=>F=F+360” DIST=” : I” FW=” :F DMS Goto 12. 正算子程序 (SUB1A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:Z3=0.5:X=U+W(Acos (G+QEKW(C+KWD+Bcos(G+QELW(C+LWD+Ncos(

6、G+QEZ3W(C+Z3WD+Bcos(G+QE(1-LW(C+(1-LWD+Acos(G+QE(1-KW(C+(1-KWD: Y=V+W(Asin(G+QEKW(C+KWD+Bsin(G+QELW(C+LWD+ Nsin(G+QEZ3W(C+Z3WD+Bsin (G+QE(1-LW(C+(1-LWD+Asin (G+QE(1-KW(C+(1-KWD:M ” ANG=” : F=G+QEW(C+WD+M:X=X+ZcosF:Y=Y+ZsinF3. 反算子程序 (SUB2M ” ANG=” :T=G-M:W=Abs(Y-VcosT-(X-UsinT:Z=0:Lbl 0:Prog "SU

7、B1" :L=T+QEW(C+WD :Z=(Z5-YcosL-(Z4-XsinL:AbsZ<1E-6=>Goto1: >W=W+Z:Goto 0Lbl 1:Z=0:Prog "SUB1" :Z=(Z5-Y÷sinF4. 垂距計(jì)算子程序 (ZX1Z6=(Z5-Z12COS A-(Z4-Z11SIN AZ7=(Z5-Z14COS B-(Z4-Z13SIN B5. 曲線元要素?cái)?shù)據(jù)庫:DA T1Lbl1:J=1=>U=Z16:V=Z17:O=Z18:G=Z19:H=Z20:P=Z21:R=Z22:Q=Z23J=2=>U=Z24:V=

8、Z25:O=Z26:G=Z27:H=Z28:P=Z29:R=Z30:Q=Z31J=3=> U=Z32:V=Z33:O=Z34:G=Z35:H=Z36:P=Z37:R=Z38:Q=Z39J=4=> U=Z40:V=Z41:O=Z42:G=Z43:H=Z44:P=Z45:R=Z46:Q=Z47J=N=> U=Z8(N+1 :V=Z8(N+1 +1:O=Z8(N+1 +2:G=Z8(N+1 +3:H=Z8(N+1 +4:P=Z8 (N+1 +5:R=Z8(N+1 +6:Q=Z8(N+1 +7(注:如有多個(gè)曲線元要素繼續(xù)添加入數(shù)據(jù)庫 DAT1中 Lbl 2 :U"X0&qu

9、ot; :V"Y0" :O"S0" :G"F0" :H"LS" :P"R0" :R"RN" :Q :C=1÷P :D=(P-R÷(2HPR:E=1 80÷三、使用說明1、規(guī)定(1 以道路中線的前進(jìn)方向(即里程增大的方向區(qū)分左右;當(dāng)線路往左偏時(shí), Q=-1;當(dāng)線路往右偏 時(shí), Q=1;當(dāng)線路為直線時(shí), Q=0。(2 當(dāng)所求點(diǎn)位于中線時(shí), Z=0;當(dāng)位于中線左鍘時(shí), Z 取負(fù)值;當(dāng)位于中線中線右側(cè)時(shí), Z 取正值。(3 當(dāng)線元為直線時(shí),其起點(diǎn)、止點(diǎn)的

10、曲率半徑為無窮大,以 10的 45次代替。(4 當(dāng)線元為圓曲線時(shí),無論其起點(diǎn)、止點(diǎn)與什么線元相接,其曲率半徑均等于圓弧的半徑。(5 當(dāng)線元為完整緩和曲線時(shí),起點(diǎn)與直線相接時(shí),曲率半徑為無窮大,以 10的 45次代替;與圓曲 線相接時(shí),曲率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時(shí),曲率半徑為無窮大,以 10的 45次代替;與 圓曲線相接時(shí),曲率半徑等于圓曲線的半徑。(6 當(dāng)線元為非完整緩和曲線時(shí), 起點(diǎn)與直線相接時(shí), 曲率半徑等于設(shè)計(jì)規(guī)定的值; 與圓曲線相接時(shí), 曲率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時(shí),曲率半徑等于設(shè)計(jì)規(guī)定的值;與圓曲線相接時(shí),曲率半 徑等于圓曲線的半徑。(7曲線元要素?cái)?shù)據(jù)庫

11、(DA T1 可根據(jù)線型不同分為各個(gè)線元段輸入到 DAT1中, 即分為直線段、 緩和曲線、 圓曲線等。(8正算時(shí)可僅輸入里程 Lp 和邊距 Dp 及右交角 ANG(普通情況下,此角為 90度,特殊情況 下當(dāng)斷面與主線不正交時(shí),亦可能此角不為 90度 全線計(jì)算,反算時(shí)通過輸入測(cè)點(diǎn)的 X 、 Y 坐標(biāo)和右交角 ANG 后計(jì)算器自動(dòng)判斷該點(diǎn)所屬曲線元并利用該線元的曲線要素執(zhí)行反算中樁里程 Lp 及支距 Dp (缺點(diǎn)是使用 CACIO4850計(jì)算器反算時(shí)較慢, PC-E500或 PDA 計(jì)算速度較快 。 2、輸入與顯示說明(一 、 輸入部分:1. SZ => XY2. XY = > SZ1

12、、 N ? 選擇計(jì)算方式,輸入 1表示進(jìn)行由里程、邊距計(jì)算坐標(biāo) ;輸入 2表示由坐標(biāo)反算 里程和邊距 2、 X0 ?線元起點(diǎn)的 X 坐標(biāo)3、 Y0 ?線元起點(diǎn)的 Y 坐標(biāo)4、 S0 ?線元起點(diǎn)里程5、 F0 ?線元起點(diǎn)切線方位角6、 LS ?線元長度7、 R0 ?線元起點(diǎn)曲率半徑8、 RN ?線元止點(diǎn)曲率半徑9、 Q ? 線 元左右偏標(biāo)志 (左偏 Q=-1,右偏 Q=1,直線段 Q=010、 S ? 正算時(shí)所求點(diǎn)的里程11、 Z ?正算時(shí)所求點(diǎn)距中線的邊距 (左側(cè)取負(fù),值右側(cè)取正值,在中線上取零 12、 ANG ?正算邊樁時(shí)左右邊樁連線與線路中線的右交角13、 J ?曲線元數(shù)據(jù)庫曲線段判斷系數(shù)

13、(J=1、 2 .n 14、 X ?反算時(shí)所求點(diǎn)的 X 坐標(biāo)15、 Y ?反算時(shí)所求點(diǎn)的 Y 坐標(biāo)16、 M? 斜交右角 (普通情況下,此角一般為 90017、 Z1-測(cè)站點(diǎn) X 坐標(biāo)(運(yùn)行程序前先預(yù)置于計(jì)算器內(nèi)18、 Z2- 測(cè)站點(diǎn) Y 坐標(biāo)(運(yùn)行程序前先預(yù)置于計(jì)算器內(nèi)19、 A 、 B 、 N 是 Gauss-Legendre 求積公式中的插值系數(shù) 20、 K 、 L 、 Z3 是 Gauss-Legendre 求積公式中 的求積節(jié)點(diǎn)21、 N=1時(shí):Z8、 Z9、 Z10、 Z11、 Z12、 Z13、 Z14、 Z15分別是各曲線元的 X 、 Y 坐標(biāo)、起始 里程 S0、初始方位角 F

14、0、線元長度 Ls 、起點(diǎn)半徑 R0、終點(diǎn)半徑 Rn 、曲線偏向系數(shù) QN=n時(shí):Z8N、 Z8N+1、 Z8N+2、 Z8N+3、 Z8N+4、 Z8N+5、 Z8N+6、 Z8N+7分別是各 22、曲 線元的 X 、 Y 坐標(biāo)、起始里程 S0、初始方位角 F0、線元長度 Ls 、起點(diǎn)半徑 R0、終點(diǎn)半徑 Rn 、曲線偏向系 數(shù) Q23、正算中的 Z18、 Z26、 Z34、 Z42 Z8(N+1 +2分別為各線元中的終點(diǎn)里程24、主程序中的“ Defm 50”應(yīng)視具體線元的數(shù)量,適當(dāng)增加擴(kuò)展變量。25、本次修改版不考慮 CASIO fx-4500PA,因其容量太小,擴(kuò)展變量實(shí)在太少而不能勝

15、任。(二 、顯示部分:XS=×××正算時(shí),計(jì)算得出的所求點(diǎn)的 X 坐標(biāo)YS=×××正算時(shí),計(jì)算得出的所求點(diǎn)的 Y 坐標(biāo)S=×××反算時(shí),計(jì)算得出的所求點(diǎn)的里程Z=×××反算時(shí),計(jì)算得出的所求點(diǎn)的邊距四、算例某匝道的由五段線元 (直線 +完整緩和曲線 +圓曲線 +非完整緩和曲線 +直線 組成, 各段線元的要素 (起 點(diǎn)里程 S0、起點(diǎn)坐標(biāo) X0 Y0、起點(diǎn)切線方位角 F0、線元長度 LS 、起點(diǎn)曲率半徑 R0、止點(diǎn)曲率半徑 RN 、 線 元左右偏標(biāo)志 Q 如下:S0 X0 Y0

16、 F0 LS R0 RN Q500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E45 0769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.75 -1806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.75 -1919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9 579.228 -1999.812 19744.038 28781.659 80 40

17、 50.00 100.000 1E45 1E45 01、正算(注意:略去計(jì)算方式及線元要素輸入,請(qǐng)自行根據(jù)所求點(diǎn)所在的線元輸入線元要素S=700 Z=-5 計(jì)算得 XS=19831.41785 YS=28509.72590S=700 Z=0 計(jì)算得 XS=19827.33592 YS=28506.83837S=700 Z= 5 計(jì)算得 XS=19823.25398 YS=28503.95084S=780 Z=-5 計(jì)算得 XS=19785.25749 YS=28575.02270S=780 Z=0 計(jì)算得 XS=19781.15561 YS=28572.16358S=780 Z= 5 計(jì)算得

18、XS=19777.05373 YS=28569.30446S=870 Z=-5 計(jì)算得 XS=19747.53609 YS=28654.13091S=870 Z=0 計(jì)算得 XS=19742.68648 YS=28652.91379S=870 Z= 5 計(jì)算得 XS=19737.83688 YS=28651.69668S=940 Z=-5.123 計(jì)算得 XS=19741. 59118 YS=28722.05802S=940 Z=0 計(jì)算得 XS=19736.47687 YS=28722.35642S=940 Z= 3.009 計(jì)算得 XS=19733.47298 YS=28722.531682、 反算X=19831.418 Y=28509.726 計(jì)算得 S=699.9999974 Z= -5 .00018164X=19827.336 Y=28506.838 計(jì)算得 S=699.9996493 Z= 0.000145136X=19823.25398 Y=28503.95084 計(jì)算得 S=6