曲線任意里程中邊樁坐標(biāo)正反算

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

2、ASIO fx-4500P計算器及 CASIO fx-4850P計算器上運行。 特別申明：(1). 適用于弧長小于2倍半徑的各種線元坐標(biāo)正反算，精度優(yōu) 于1mm； (2). 在引用該核心計算部分時，請注明來源。 二、源程序 1.主程序(TYQXJS) 1.SZ = XY：2.XY = SZ：N：UX0：VY0：OS0：GF0：HLS：PR0：R RN：Q：C=1P：D=(P-R)(2HPR)：E=180：N=1=Goto 1：Goto 2 Lbl 1：SZ：SZ：W=Abs(S-O)：Prog SUB1：XXS=X YYS=Y FFS=F-90 Goto 1 Lbl 2：XY：XY：I=X：J

3、=Y：Prog SUB2：SS=O+W ZZ=Z Goto 2 2. 正算子程序(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)

4、+90：X=X+ZcosF：Y=Y+ZsinF 3. 反算子程序(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：AbsZGoto1：W=W+Z：Goto 0 Lbl 1：Z=0：Prog SUB1：Z=(J-Y)sinF 三、使用說明 1、規(guī)定 (1) 以道路中線的前進(jìn)方向（即里程增大的方向）區(qū)分左右；當(dāng)線元往左偏時， Q=-1；當(dāng)線元往右偏時，Q=1；當(dāng)線元為直線時，Q=0。 (2) 當(dāng)所求點位于中線時，Z=0；當(dāng)位于中線左鍘時，Z取負(fù)值；當(dāng)位于中

5、線中線右 側(cè)時，Z取正值。 (3) 當(dāng)線元為直線時，其起點、止點的曲率半徑為無窮大，以10的45次代替。 (4) 當(dāng)線元為圓曲線時，無論其起點、止點與什么線元相接，其曲率半徑均等于圓 弧的半徑。 (5) 當(dāng)線元為完整緩和曲線時，起點與直線相接時，曲率半徑為無窮大，以10的45 次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半 徑為無窮大，以10的45次代替；與圓曲線相接時，曲率半徑等于圓曲線的半徑。 (6) 當(dāng)線元為非完整緩和曲線時，起點與直線相接時，曲率半徑等于設(shè)計規(guī)定的 值；與圓曲線相接時，曲率半徑等于圓曲線的半徑。止點與直線相接時，曲率半徑等 于設(shè)計規(guī)定的值；與

6、圓曲線相接時，曲率半徑等于圓曲線的半徑。 2、輸入與顯示說明 輸入部分： 1. SZ = XY 2. XY = SZ N ? 選擇計算方式，輸入1表示進(jìn)行由里程、邊距計算坐標(biāo) ；輸入2表示由坐標(biāo)反算 里程和邊距。 X0 ？線元起點的X坐標(biāo) Y0 ？線元起點的Y坐標(biāo) S0 ？線元起點里程 F0 ？線元起點切線方位角 LS ？線元長度 R0 ？線元起點曲率半徑 RN ？線元止點曲率半徑 Q ？ 線 元左右偏標(biāo)志(左偏Q=-1，右偏Q=1，直線段Q=0） S ？ 正算時所求點的里程 Z ？正算時所求點距中線的邊距(左側(cè)取負(fù)，值右側(cè)取正值，在中線上取零) X ？反算時所求點的X坐標(biāo) Y ？反算時所求點的Y坐標(biāo) 顯示部分： XS= 正算時，計算得出的所求點的X坐標(biāo) YS= 正算時，計算得出的所求點的Y坐標(biāo) FS= 正算時，所求點對應(yīng)的中線點的切線方位角 S= 反算時，計算得出的所求點的里程 Z= 反算時，計算得出的所求點的邊距 U=起點X坐標(biāo)：V=起點Y坐標(biāo)：O=起點樁號：G=起點方位角：H=起點至終點之間的距離：P=半徑：R=半徑： 注：（當(dāng)為圓曲線段時P，R均為該段的半徑，當(dāng)為直線時，P，R均以10