曲線任意里程中邊樁坐標(biāo)正反算(CASIO fx-4800P計(jì)算器)程序

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

2、(1). 適用于弧長(zhǎng)小于2倍半徑的各種線元坐標(biāo)正反算，精度優(yōu) 于1mm； (2). 在引用該核心計(jì)算部分時(shí)，請(qǐng)注明來(lái)源。二、源程序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=>Got

3、o 1：>Goto 2Lbl 1：SZ：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=XY"YS"=YF"FS"=F-90Goto 1Lbl 2：XY：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+WZ"Z"=ZGoto 22. 正算子程序(SUB1)A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：F=1-L：M=1-K：X=U+W(Acos(G+

4、QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Bcos(G+QEFW(C+FWD)+Acos(G+QEMW(C+MWD)：Y=V+W(Asin(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. 反算子程序(SUB2)T=G-90：W=Abs(Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sin

5、L：AbsZ<1E-6=>Goto1：>W=W+Z：Goto 0Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF三、使用說(shuō)明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)的曲率半徑為無(wú)窮大，以10的45次代替。(4) 當(dāng)線元為圓曲線時(shí)，無(wú)論其起點(diǎn)、止點(diǎn)與什么線元相接，其曲率半徑均等于圓弧的半徑。(5) 當(dāng)

6、線元為完整緩和曲線時(shí)，起點(diǎn)與直線相接時(shí)，曲率半徑為無(wú)窮大，以10的45次代替；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時(shí)，曲率半徑為無(wú)窮大，以10的45次代替；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。(6) 當(dāng)線元為非完整緩和曲線時(shí)，起點(diǎn)與直線相接時(shí)，曲率半徑等于設(shè)計(jì)規(guī)定的值；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。止點(diǎn)與直線相接時(shí)，曲率半徑等于設(shè)計(jì)規(guī)定的值；與圓曲線相接時(shí)，曲率半徑等于圓曲線的半徑。2、輸入與顯示說(shuō)明 輸入部分： 1. SZ => XY 2. XY = > SZ N ? 選擇計(jì)算方式，輸入1表示進(jìn)行由里程、邊距計(jì)算坐標(biāo) ；輸入2表示由坐標(biāo)反算里

7、程和邊距。 X0 ？線元起點(diǎn)的X坐標(biāo) Y0 ？線元起點(diǎn)的Y坐標(biāo) S0 ？線元起點(diǎn)里程 F0 ？線元起點(diǎn)切線方位角 LS ？線元長(zhǎng)度 R0 ？線元起點(diǎn)曲率半徑 RN ？線元止點(diǎn)曲率半徑 Q ？ 線 元左右偏標(biāo)志(左偏Q=-1，右偏Q=1，直線段Q=0） S ？ 正算時(shí)所求點(diǎn)的里程 Z ？正算時(shí)所求點(diǎn)距中線的邊距(左側(cè)取負(fù)，值右側(cè)取正值，在中線上取零) X ？反算時(shí)所求點(diǎn)的X坐標(biāo) Y ？反算時(shí)所求點(diǎn)的Y坐標(biāo) 顯示部分： XS=××× 正算時(shí)，計(jì)算得出的所求點(diǎn)的X坐標(biāo) YS=××× 正算時(shí)，計(jì)算得出的所求點(diǎn)的Y坐標(biāo) FS=×

8、15;× 正算時(shí)，所求點(diǎn)對(duì)應(yīng)的中線點(diǎn)的切線方位角 S=××× 反算時(shí)，計(jì)算得出的所求點(diǎn)的里程 Z=××× 反算時(shí)，計(jì)算得出的所求點(diǎn)的邊距四、算例 某匝道的由五段線元（直線+完整緩和曲線+圓曲線+非完整緩和曲線+直線）組成，各段線元的要素（起點(diǎn)里程S0、起點(diǎn)坐標(biāo)X0 Y0、起點(diǎn)切線方位角F0、線元長(zhǎng)度LS、起點(diǎn)曲率半徑R0、止點(diǎn)曲率半徑RN、線 元左右偏標(biāo)志Q）如下：S0 X0 Y0 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 9579.228 -1999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E45 0 1、正算 （注意：略去計(jì)算方式及線元要素輸入，請(qǐng)自行根據(jù)所