【正文】 0(B+C) 247。2688247。5800公路線路坐標(biāo)正反算程序時間:20100801 18:59:01 R ^(3) →Z[6] ↙C2 247。2→Q↙(R+(Z[6]+Z[7])247。40247。 選擇線路，左幅=1，右幅=1，整體式=0 正算子程序 GSZS ((PR)247。π)+Z[4]cos(G+QZ[1]J(1247。π))) →Y： (G+QJ(1247。 3) 當(dāng)線元為直線時，其起點、止點的曲率半徑為無窮大，以10的45次代替。 80 40 π)) →Z[6]: (Z[4]cos(G+QZ[1]J(1247。π)) →Z[16]:( Z[12]+ Z[13]+ Z[14]+ Z[15]+ Z[16])→Z[17] : (JZ[17]) →Z[18]:(S+ Z[18]) →Y: (G+M) →Z[19]:(( Z[19]+QJ(1247。 輸入與顯示說明 （1）輸入部分： I ？線元起點的X坐標(biāo) S？線元起點的Y坐標(biāo) O？線元起點里程G ？線元起點切線方位角 H ？線元終點里程 P ？線元起點曲率半徑 R？線元止點曲率半徑 Q ？ 線 元左右偏標(biāo)志(左偏Q=1，右偏Q=1，直線段Q=0)K ？ 正算時所求點的里程L ？ 正算時所求點距中線的邊距(左側(cè)取負值，右側(cè)取正值，在中線上取零)ANG？正算邊樁時左右邊樁連線與線路中線的右交角 X ？ 反算時所求點的X坐標(biāo)Y ？ 反算時所求點的Y坐標(biāo) M ? 斜交右角 （2）顯示部分 X正算時，計算得出的所求點的X坐標(biāo)Y正算時，計算得出的所求點的Y坐標(biāo)K反算時，計算得出的所求點的里程 L反算時，計算得出的所求點的邊距 三、算例 某匝道的由五段線元（直線+完整緩和曲線+圓曲線+非完整緩和曲線+直線）組成，各段線元的要素（起點里程S0、起點坐標(biāo)X0 Y0、起點切線方位角F0、線元長度LS、起點曲率半徑R0、止點曲率半徑RN、線元左右偏標(biāo)志Q）如下：S0 X0 Y0 F0 LS R0 RN Q 125 16 1E45 1E45 0 125 16 1E45 1 120 25 1 91 17 1 80 40 1E45 1E45 0 （注：該算例中線元要素Ls為程序修改前須輸入的線元長度，程序修改后改為輸入線元終點里程KN）。π)) →F： (X+LcosF)→X：(Y+LsinF) →YReturn4. 線元要素數(shù)據(jù)庫：（DATMDZD ）If K≥(線元起點里程)And K＜(線元終點里程):Then (線元起點橫坐標(biāo)x) →I：(線元起點縱坐標(biāo)y)→S: (線元起點里程)→O: (線元起點方位角)→G：(線元終點里程)→H：(線元起點半徑)→P：(線元終點半徑)→R：Q(Z=1 ZX=0 Y=1)(偏轉(zhuǎn)系數(shù))→Q: IfEnd ……….(注：如有多個曲線元要素繼續(xù)添加入數(shù)據(jù)庫中) 說明： 一、程序功能及原理 ： 本程序由兩個主程序——正算主程序(GSZS)、反算主程序(GSFS)和兩個子程——正算子程序(SUB1)、線元數(shù)據(jù)庫(DATM)構(gòu)成，可以根據(jù)曲線段——直線、圓曲線、緩和曲線（完整或非完整型）的線元要素（起點坐標(biāo)、起點里程、起點切線方位角、終點里程、起點曲率半徑、止點曲率半徑）及里程邊距或坐標(biāo)，對該曲線段范圍內(nèi)任意里程中邊樁坐標(biāo)進行正反算。π)) →Z[7]:( Bcos(G+Q(1E)J(1247。 (4) 當(dāng)線元為圓曲線時，無論其起點、止點與什么線元相接，其曲率半徑均等于圓弧的半徑。π+M) →F： (X+LcosF)→X：(Y+LsinF) →Y4. 曲線元要素數(shù)據(jù)庫：ZX/YX/QX If K(起點里程):Then Goto 2：IfEnd If K（ ZDLC）: Then QDXO →I：QDY0→S:QDLC→O:QDFWJ →G：ZDLC→H：QDR→P：ZDR→R：Q(Z=1 ZX=