卡西歐4800P曲線(xiàn)任意里程中邊樁坐標(biāo)正反算

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

2、曲 線(xiàn)任意里程中邊樁坐標(biāo)進(jìn)行正反算。本程序也可以在CASIO fx-4500P 計(jì)算器及 CASIO fx-4850P 計(jì)算器上運(yùn)行。特別申明：(1). 適用于弧長(zhǎng)小于 2 倍半徑的各種線(xiàn)元坐標(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" ：

3、P"R0" ： R"RN" ： Q ： C=1÷ P ： D=(P-R) ÷ (2HPR) ：E=180 ÷：N=1=>Goto 1 ： >Goto 2 Lbl 1 ：SZ ：SZ ：W=Abs(S-O) ：Prog "SUB1" ：X"XS"=XY"YS"=Y F"FS"=F-90 Goto 1 Lbl 2 ：XY ：XY ：I=X ：J=Y ：Prog "SUB2" ：S"S"=O+W Z&

4、quot;Z"=Z Goto 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+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+Zcos

5、F ：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)sinL ： AbsZ<1E-6=>Goto1 ： >W=W+Z ： Goto 0 Lbl 1 ：Z=0 ：Prog "SUB1" ：Z=(J-Y) ÷sinF三、使用說(shuō)明1 、規(guī)定(1) 以道路中線(xiàn)的前進(jìn)方向(即里程增大的方向)區(qū)分左右；當(dāng)線(xiàn) 元往左偏時(shí)，Q=-1 ；當(dāng)線(xiàn)元往右偏時(shí)，Q=1 ；當(dāng)

6、線(xiàn)元為直線(xiàn)時(shí)，Q=0 。(2) 當(dāng)所求點(diǎn)位于中線(xiàn)時(shí)， Z=0 ；當(dāng)位于中線(xiàn)左鍘時(shí)， Z 取負(fù)值； 當(dāng)位于中線(xiàn)中線(xiàn)右側(cè)時(shí)， Z 取正值。(3) 當(dāng)線(xiàn)元為直線(xiàn)時(shí)，其起點(diǎn)、止點(diǎn)的曲率半徑為無(wú)窮大，以 10 的 45 次代替。(4) 當(dāng)線(xiàn)元為圓曲線(xiàn)時(shí)，無(wú)論其起點(diǎn)、止點(diǎn)與什么線(xiàn)元相接，其曲 率半徑均等于圓弧的半徑。(5) 當(dāng)線(xiàn)元為完整緩和曲線(xiàn)時(shí)，起點(diǎn)與直線(xiàn)相接時(shí)，曲率半徑為無(wú) 窮大，以 10 的 45 次代替；與圓曲線(xiàn)相接時(shí)，曲率半徑等于圓曲線(xiàn)的半 徑。止點(diǎn)與直線(xiàn)相接時(shí)，曲率半徑為無(wú)窮大，以 10 的 45 次代替；與圓 曲線(xiàn)相接時(shí)，曲率半徑等于圓曲線(xiàn)的半徑。(6) 當(dāng)線(xiàn)元為非完整緩和曲線(xiàn)時(shí)，起點(diǎn)與直線(xiàn)

7、相接時(shí)，曲率半徑等 于設(shè)計(jì)規(guī)定的值；與圓曲線(xiàn)相接時(shí)，曲率半徑等于圓曲線(xiàn)的半徑。止點(diǎn)與 直線(xiàn)相接時(shí)，曲率半徑等于設(shè)計(jì)規(guī)定的值；與圓曲線(xiàn)相接時(shí)，曲率半徑等 于圓曲線(xiàn)的半徑。2 、輸入與顯示說(shuō)明輸入部分：1. SZ => XY2. XY = > SZN ? 選擇計(jì)算方式，輸入 1 表示進(jìn)行由里程、邊距計(jì)算坐標(biāo) ；輸入2 表示由坐標(biāo)反算 里程和邊距。X0 ？線(xiàn)元起點(diǎn)的 X 坐標(biāo)Y0 ？線(xiàn)元起點(diǎn)的 Y 坐標(biāo)S0 ？線(xiàn)元起點(diǎn)里程F0 ？線(xiàn)元起點(diǎn)切線(xiàn)方位角LS ？線(xiàn)元長(zhǎng)度R0 ？線(xiàn)元起點(diǎn)曲率半徑RN ？線(xiàn)元止點(diǎn)曲率半徑Q ？ 線(xiàn) 元左右偏標(biāo)志（左偏Q=-1 ，右偏Q=1 ，直線(xiàn)段Q=0 ）S ？

8、 正算時(shí)所求點(diǎn)的里程Z ？正算時(shí)所求點(diǎn)距中線(xiàn)的邊距 （ 左側(cè)取負(fù)，值右側(cè)取正值，在中線(xiàn)上取 零）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=× ×× 正算時(shí)，所求點(diǎn)對(duì)應(yīng)的中線(xiàn)點(diǎn)的切線(xiàn)方位角S=× ×× 反算時(shí)，計(jì)算得出的所求點(diǎn)的里程Z=× ×× 反算時(shí)，計(jì)算得出的所求點(diǎn)的邊距四、算例 某匝道的由五段線(xiàn)元（直線(xiàn)+完

9、整緩和曲線(xiàn)+ 圓曲線(xiàn)+ 非完整緩和曲 線(xiàn)+直線(xiàn)）組成，各段線(xiàn)元的要素（起點(diǎn)里程S0 、起點(diǎn)坐標(biāo)X0 Y0 、起點(diǎn)切線(xiàn)方位角 F0 、線(xiàn)元長(zhǎng)度 LS 、起點(diǎn)曲率半徑 R0 、止點(diǎn)曲率半徑 RN、 線(xiàn) 元左右偏標(biāo)志 Q ）如下：S0 X0 Y0 F0LS R0 RN Q500.000 19942.837 28343.561 12516 .00 269.2561E45 1E45 0769.25619787.34028563.3781251631.0037.4921E45 221.75 -1806.74819766.56628594.5741202554.07112.779221.75 221.75

10、-1919.52719736.07228701.893911730.6380.285221.75 9579.228-1999.81219744.03828781.659804050.00100.0001E45 1E45 01 、正算注意：略去計(jì)算方式及線(xiàn)元要素輸入，請(qǐng)自行根據(jù)所求點(diǎn)所在的線(xiàn)元輸入線(xiàn)元要素）S=700 Z=-5得 XS=19831.41785YS=28509.72590S=700 Z=0得 XS=19827.33592YS=28506.83837S=700 Z=5 計(jì)得 XS=19823.25398YS=28503.95084S=780 Z=-5算算算算得 XS=19785.25

11、749 YS=28575.02270S=780 Z=0計(jì)算得 XS=19781.15561YS=28572.16358S=780 Z=5計(jì)算得 XS=19777.05373YS=28569.30446S=870 Z=-5計(jì)算得 XS=19747.53609YS=28654.13091S=870 Z=0計(jì)算得 XS=19742.68648YS=28652.91379S=870 Z=5計(jì)算得 XS=19737.83688YS=28651.69668S=940 Z=-5.123計(jì)算得XS=19741.59118 YS=28722.05802S=940 Z=0計(jì)算得 XS=19736.47687YS=

12、28722.35642S=940 Z=3.009計(jì)算得 XS=19733.47298YS=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=699.9999985 Z= 5.000003137X=19785.25749 Y=28575.02270計(jì)算得 S=780.0000035 Z= -5 .000001663X=19781.15561Y=28572.16358計(jì)算得 S=780.0000025 Z=- 0.000002979X=19777.05373Y=28569.30446計(jì)算得 S=780.0000016 Z= 4.99999578X=19747.536 Y=28654.131計(jì)算得 S=870.0001137 Z= -4.99941049X=19742.686 Y=28652.914計(jì)算得 S=870.0003175 Z=- 0.00041814X=19737.837Y=28651.69