曲線任意里程中邊樁坐標(biāo)正反算及放樣fx-4850程序

1、曲線任意里程中邊樁坐標(biāo)正反算及放樣fx-4850程序第五次修改李正云本程序是在yshf及hangniu1973兩位師傅的成果上作的一次改動(dòng)，修改內(nèi)容主要有一下幾點(diǎn)：1、在變量符號上修改，目的是更能與工地實(shí)際結(jié)合；2、把原主程序分為兩個(gè)，更能靈活調(diào)用；3、參加放樣程序，做到坐標(biāo)計(jì)算到放樣一體化；4、使得整個(gè)測量放樣過程更加簡單。GAUSSL坐標(biāo)正反算fx 4850程序源程序1. 正算主程序GSZSI"X0" : S"Y0": OKO" : G"F0": H"KN": P"R0": R&q

2、uot;RN": Q' Q(-Z+Y) ：D=(P-R)十(2(H-OPR):KL L( -Z+Y) ：M ANG(YJJ) =90:注：此處假設(shè)不給 M賦值，那么可計(jì)算斜交點(diǎn)J=Abs(K-O : Prog"SUB1": FWJ= :F=F-M: X= ：X=X:Pause0: 丫二 ：丫二丫丄2. 反算主程序GSFSXY: Z2=X : Z3 =Y :I"X0" : S"Y0": OK0" : G"F0" : H"KN": P"R0" : R&

3、quot;RN": Q' Q(-Z+Y) :D=(P-R)十(2(H-OPR)：J=Abs(Y-S)cos(G-90) (X-I)si n(G-90)L=0: M M(YJJ) =90:Lbl 0 : Prog "SUB1":L=(Z3-Y)cos(G- 90+QJ(1 十 P+JD)X 180-n) -(Z2-X)sin(G- 90+QJ(1 十 P+JD) X 180十 n ):AbsL<1E-6=>Goto1:豐 >J=J+L: Goto 0 -1Lbl 1 : L=0 : Prog "SUB1": L=(Z3-

4、 Y) - si nF: K= :K= O+J:Pause0: L= :L=L A3. 正算子程序(SUB1)Defm 4:A=0.1184634425:B=0.2393143352:Z4=0.2844444444:C=0.0469100770:E=0.2307653449:Z1=0.5:X=l+J(Acos(G+QCJ(1 十 P+CJDX 180十 n )+Bcos(G+QEJ(1-P+EJDX 180十 n )+Z4 cos(G+QZ1J(1- P+Z1JD) X 180十 n )+Bcos(G+Q(1 - E)J(1 - P+(1-E)JD) X180十 n )+Acos(G +Q (

5、1- C)J(1 - P+(1-C)JD) X 180十 n ):Y=S+J(Asin(G+QCJ(1 十 P+CJDX 180十 n )+Bsin(G+QEJ(1 - P+EJD)X 180十 n )+Z 4sin(G+QZ1J(1 - P+Z1JD) X 180十 n )+Bsin(G+Q(1 - E)J(1 - P+(1-E)JD) X 180十 n )+Asin(G+Q (1 - C)J(1 - P+(1-C)JD) X 180十 n ):F=G+QJ(仟 P+JD) X 180十 n +M X=X+LcosF Y=Y+LsinF4. 曲線元要素?cái)?shù)據(jù)庫：DAT-MK> O=>

6、;K< H=> I=*:S=*:O=*:G=*:H=*:P=*:R=*:Q=*K> O=>K< H=> I=*:S=*:O=*:G=*:H=*:P=*:R=*:Q=*K> O=>K< H=> I=*:S=*:O=*:G=*:H=*:P=*:R=*:Q=*K> O=>K< H=> I=*:S=*:O=*:G=*:H=*:P=*:R=*:Q=*K> O=>KC H=> I=*:S=*:O=*:G=*:H=*:P=*:R=*:Q=*K> O=>KC H=>I=*：S=*:o=*:g

7、=*:h=*:p=*:R=*:Q=*“v注：如有多個(gè)曲線元要素繼續(xù)添參加數(shù)據(jù)庫DAT-M中5、M線坐標(biāo)正算組合程序MG-ZBProgDAT-M:Prog GSZS6 M線坐標(biāo)計(jì)算放樣組合程序MG-FYProgMG-ZB:Prog LTKZD: Prog FYM-GSFB7、M線坐標(biāo)反算組合程序ProgDAT-M:Prog GSFS 說明：一、程序功能及原理1 .功能說明：本程序由兩個(gè)主程序一一正算主程序GSZS反算主程序GSFS和兩個(gè)子程一一 正算子程序SUB1、線元數(shù)據(jù)庫DAT-M構(gòu)成，可以根據(jù)曲線段一一直線、圓曲 線、緩和曲線完整或非完整型的線元要素起點(diǎn)坐標(biāo)、起點(diǎn)里程、起點(diǎn)切線 方位角、終

8、點(diǎn)里程、起點(diǎn)曲率半徑、止點(diǎn)曲率半徑及里程邊距或坐標(biāo)，對該曲 線段范圍內(nèi)任意里程中邊樁坐標(biāo)進(jìn)行正反算。本程序可以在 CASIO fx-4800P 計(jì) 算器及CASIO fx-4850P計(jì)算器上運(yùn)行。由于參加了數(shù)據(jù)庫DAT-M,可實(shí)現(xiàn)坐 標(biāo)正反算的全線貫穿。組合程序5可實(shí)現(xiàn)M線的正算貫穿，組合程序7可實(shí)現(xiàn)M線的反算貫穿， 組合程序 6可實(shí)現(xiàn)坐標(biāo)計(jì)算到放樣一體化。2計(jì)算原理：利用 Gauss-Legendre 5 點(diǎn)通用公式計(jì)算線路中邊樁坐標(biāo)并計(jì)算放樣數(shù)據(jù)。利用待求點(diǎn)至線元起點(diǎn)切線作垂線，逐次迭代趨近原理反算里程及邊距。二、使用說明1、規(guī)定(1) 以道路中線的前進(jìn)方向 即里程增大的方向 區(qū)分左右；

9、當(dāng)線 元往左偏時(shí)，Q=-1 ;當(dāng)線元往右偏時(shí)，Q=1;當(dāng)線元為直線時(shí)，Q=0(2) 當(dāng)所求點(diǎn)位于中線時(shí)，L=0;當(dāng)位于中線左側(cè)時(shí)，L取負(fù)值； 當(dāng)位于中線右側(cè)時(shí)，L取正值。(3) 當(dāng)線元為直線時(shí)，其起點(diǎn)、止點(diǎn)的曲率半徑為無窮大，以 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)

10、線元為非完整緩和曲線時(shí)，起點(diǎn)與直線相接時(shí)，曲率半徑 等于設(shè)計(jì)規(guī)定的值; 與圓曲線相接時(shí)， 曲率半徑等于圓曲線的半徑。 止點(diǎn)與直線 相接時(shí)， 曲率半徑等于設(shè)計(jì)規(guī)定的值; 與圓曲線相接時(shí)， 曲率半徑等于圓曲線的 半徑。(7) 曲線元要素?cái)?shù)據(jù)庫DAT-M可根據(jù)線型不同分為各個(gè)線元段輸入到 DAT-M 中，即分為直線段、緩和曲線、圓曲線等。(8) 正算時(shí)可僅輸入里程和邊距及右交角可實(shí)現(xiàn)全線計(jì)算，但反算時(shí)只能通過首 先輸入里程K值讀取數(shù)據(jù)庫DAT-M計(jì)算器自動(dòng)將里程K所在線元數(shù)據(jù)賦給反算 主程序GSFSS行試算，試算出的里程和邊距須帶入正算主程序 GSZ中計(jì)算坐標(biāo)， 假設(shè)坐標(biāo)吻合那么反算正確。2、輸入與

11、顯示說明 1 輸入局部：X0 ？線元起點(diǎn)的 X 坐標(biāo)Y0 ？線元起點(diǎn)的 Y 坐標(biāo)K0 ？線元起點(diǎn)里程F0 ？線元起點(diǎn)切線方位角KN ？線元終點(diǎn)里程R0 ？線元起點(diǎn)曲率半徑RN ？線元止點(diǎn)曲率半徑Q ？ 線元左右偏標(biāo)志(左偏Q=-1，右偏Q=1,直線段Q=0)K ？ 正算時(shí)所求點(diǎn)的里程L ？ 正算時(shí)所求點(diǎn)距中線的邊距 (左側(cè)取負(fù)值，右側(cè)取正值， 在中線上取零 )ANG?正算邊樁時(shí)左右邊樁連線與線路中線的右交角X ？ 反算時(shí)所求點(diǎn)的 X 坐標(biāo)Y ? 反算時(shí)所求點(diǎn)的 Y 坐標(biāo)M ? 斜交右角線元要素?cái)?shù)據(jù)庫中K> O=>K< H=>中的0和H分別為該段線元起點(diǎn)里程 和終點(diǎn)里程A

12、、 B、Z4 是 Gauss-Legendre 求積公式中的插值系數(shù)C 、E、Z1 是 Gauss-Legendre 求積公式中的求積節(jié)點(diǎn) 2顯示局部：X=xxx正算時(shí)，計(jì)算得出的所求點(diǎn)的 X坐標(biāo)Y=xxx正算時(shí)，計(jì)算得出的所求點(diǎn)的 Y坐標(biāo)K=xxx反算時(shí)，計(jì)算得出的所求點(diǎn)的里程L=xxx反算時(shí)，計(jì)算得出的所求點(diǎn)的邊距三、算例 某匝道的由五段線元直線+完整緩和曲線 +圓曲線 +非完整緩和曲線 +直線組成， 各段線元的要素起點(diǎn)里程 SO、起點(diǎn)坐標(biāo)X0 丫0、起點(diǎn)切線方位角F0、線元長度LS起點(diǎn)曲率半徑R0止點(diǎn)曲率半徑RN線元左右偏標(biāo)志 Q如下：S0X0Y0F0LSR0RNQ500.0001994

13、2.83728343.561125 1631.00269.2561E451E450769.25619787.34028563.378125 1631.0037.4921E45221.75-1806.74819766.56628594.574120 2554.07112.779221.75221.75-1919.52719736.07228701.893911730.6380.285221.759579.228-1999.81219744.03828781.65980 4050.00100.0001E451E450注：該算例中線元要素 Ls 為程序修改前須輸入的線元長度，程序修改后改為 輸入線元

14、終點(diǎn)里程 KN放樣程序 FYUXJ:VYJ:WXH:Z YH:XY:I=W -U:J=Z-V:Pol(I,J):J<0=>J=J+360A“ HSJL= :1 丄O' HF' =J:X=0=>U=U+LCos(J+P/'):V=V+LSin(J+P'/'): X= ：U:PauseO: 丫二 :V:工>I=X-U:J=Y- V: Pol(l,J):J<0=>J=J+360:AH' QF' =J:J=H-O:J<0=>J=J+360: A“ANG= :J f DMS該處不輸Pause。語句亦可

15、實(shí)現(xiàn)“ ANG=和“ QSJL=的同時(shí) 出現(xiàn)“QSJL=:l凌鐵大橋控制點(diǎn)數(shù)據(jù)庫 LTKZDA=21:B=22:C=23:D=24:E=25:F=26:G=27:H=28:I=29T “ZJDH :T=11=>U=*:V=*T=12=>U=*:V=*T=13=>U=*:V=*T=15=>U=*:V=*T= 4 =>U=*:V=*T= A=>U=*:V=*T= B=>U=*:V=*T= C=>U=*:V=*注：如有多個(gè)控制點(diǎn)繼續(xù)添參加數(shù)據(jù)庫 LTKZD中N “HSDH :N=11=>W=*:Z=*N=12=>W=*:Z=*N=13=&g

16、t;W=*:Z=*N=15=>W=*:Z=*N= 4 =>W=*:Z=*N= A=>W=*:Z=*N= B=>W=*:Z=*N= C=>W=*:Z=*注：如有多個(gè)控制點(diǎn)繼續(xù)添參加數(shù)據(jù)庫 LTKZD中組合程式 LTFYProg LTKZD: Prog FY 說明：一、程序功能及原理1、功能說明：FY程序有兩個(gè)功能，當(dāng)前視點(diǎn)坐標(biāo) X輸入為零時(shí)，程式執(zhí)行轉(zhuǎn)點(diǎn)功能，即首先 通過全站儀測定任意轉(zhuǎn)點(diǎn)距置鏡點(diǎn)距離 L 及置鏡點(diǎn)與轉(zhuǎn)點(diǎn)連線和置鏡點(diǎn)與后視 點(diǎn)連線的夾角P/'并輸入計(jì)算器，即可計(jì)算出任意轉(zhuǎn)點(diǎn)坐標(biāo)；另一功能是坐 標(biāo)放樣，即輸入放樣點(diǎn)X、丫坐標(biāo)，從而計(jì)算出與置鏡點(diǎn)

17、距離“QSJL=以及全站儀歸零放樣撥角“ ANG= 前視方位角后視方位角。凌鐵大橋控制點(diǎn)數(shù)據(jù)庫LTKZD功能是與FY程序配合后，可實(shí)現(xiàn)坐標(biāo)放 樣時(shí)僅輸入置鏡點(diǎn)、后視點(diǎn)的編號，從而減少其數(shù)據(jù)的大量輸入。2、計(jì)算原理： 利用坐標(biāo)計(jì)算和方位角的根本知識以及計(jì)算器極坐標(biāo)轉(zhuǎn)換功能鍵 POL二、使用說明1、規(guī)定當(dāng)前視點(diǎn)坐標(biāo) X 輸入為零時(shí)，程式執(zhí)行轉(zhuǎn)點(diǎn)功能2、輸入與顯示說明1輸入局部：XJ ？置鏡點(diǎn) X 坐標(biāo)YJ ？置鏡點(diǎn)丫坐標(biāo)XH ？后視點(diǎn) X 坐標(biāo)YH ？后視點(diǎn) 丫坐標(biāo)X？放樣點(diǎn) X 坐標(biāo)?放樣點(diǎn)丫坐標(biāo)YL ？ 實(shí)測轉(zhuǎn)點(diǎn)距置鏡點(diǎn)距離P'u ?實(shí)測置鏡點(diǎn)與轉(zhuǎn)點(diǎn)連線和置鏡點(diǎn)與后視點(diǎn)連線的夾角T ？ 置鏡點(diǎn)點(diǎn)號 ZJDHN ？ 后視點(diǎn)點(diǎn)號 HSDH11、12、13、15、4、A、B、