曲線(xiàn)任意里程中邊樁坐標(biāo)正反算(終版)

1、CASIO4800、4850坐標(biāo)正反計(jì)算程序坐標(biāo)正反計(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)

2、的程序中，用于對(duì)曲線(xiàn)任意里程中邊樁坐標(biāo)進(jìn)行正反算。本程序也可以在CASIOfx4500P計(jì)算器及CASIOfx4850P計(jì)算器上運(yùn)行。特別申明：(1).適用于弧長(zhǎng)小于2倍半徑的各種線(xiàn)元坐標(biāo)正反算，精度優(yōu)于1mm；(2).在引用該核心計(jì)算部分時(shí)，請(qǐng)注明來(lái)源。二、源程序1.主程序(TYQXJS)Lb14：“1.SZ=XY“：“2.XY=SZ“：｛NS｝：S≤下一線(xiàn)元起點(diǎn)里程＝＞O=本線(xiàn)元起點(diǎn)里程：U=本線(xiàn)元起點(diǎn)X：V=本線(xiàn)元起點(diǎn)Y：G=本線(xiàn)

3、元起算方位角：H=本線(xiàn)元長(zhǎng)度：P=起點(diǎn)曲率半徑：R=終點(diǎn)曲率半徑：Q=0或1、1：Goto0Δ←┘(第一線(xiàn)元數(shù)據(jù)要素第一線(xiàn)元數(shù)據(jù)要素)S≤下一線(xiàn)元起點(diǎn)里程＝＞O=本線(xiàn)元起點(diǎn)里程：U=本線(xiàn)元起點(diǎn)X：V=本線(xiàn)元起點(diǎn)Y：G=本線(xiàn)元起算方位角：H=本線(xiàn)元長(zhǎng)度：P=起點(diǎn)曲率半徑：R=終點(diǎn)曲率半徑：Q=0或1、1：Goto0Δ←┘(第二線(xiàn)元數(shù)據(jù)要第二線(xiàn)元數(shù)據(jù)要素)S≤下一線(xiàn)元起點(diǎn)里程＝＞O=本線(xiàn)元起點(diǎn)里程：U=本線(xiàn)元起點(diǎn)X：V=本線(xiàn)元起點(diǎn)Y：G=本

4、線(xiàn)元起算方位角：H=本線(xiàn)元長(zhǎng)度：P=起點(diǎn)曲率半徑：R=終點(diǎn)曲率半徑：Q=0或1、1：Goto0Δ←┘(第三線(xiàn)元數(shù)據(jù)要第三線(xiàn)元數(shù)據(jù)要素)。。。。。。。。。。Goto0Δ←┘(第N1線(xiàn)元數(shù)據(jù)要素線(xiàn)元數(shù)據(jù)要素)。。。。。。。。。。Goto0Δ←┘(第N線(xiàn)元數(shù)據(jù)要素線(xiàn)元數(shù)據(jù)要素)LB10：C=1PC=1P：D=(PR)(2HPR)D=(PR)(2HPR)：E=180πE=180π：N=1=GotoN=1=Goto1：≠Goto≠Goto2Δ←┘

5、2Δ←┘LblLbl1：ZZ：Z：W=Abs(SO)W=Abs(SO)：ProgProg“SUB1““SUB1“：“X=“:X◢“X=“:X◢“Y=“:Y◢“Y=“:Y◢GotoGoto4←┘4←┘LblLbl2：XYXY：XYXY：I=XI=X：J=YJ=Y：ProgProg“SUB2““SUB2“：“S=“:S=OW◢“S=“:S=OW◢“Z=“:Z◢“Z=“:Z◢GotoGoto42.正算子程序(SUB1)A=0.17392742

6、26A=0.1739274226：B=0.3260725774B=0.3260725774：K=0.0694318442K=0.0694318442：L=0.3300094782L=0.3300094782：F=1LF=1L：M=1KM=1K：X=UW(Acos(GQEKW(CKWD))Bcos(GQELW(CLWD))Bcos(GQEFWX=UW(Acos(GQEKW(CKWD))Bcos(GQELW(CLWD))Bcos(GQEFW

7、(CFWD))Acos(GQEMW(CMWD)))(CFWD))Acos(GQEMW(CMWD)))：Y=VW(Asin(GQEKW(CKWD))Bsin(GY=VW(Asin(GQEKW(CKWD))Bsin(GQELW(CLWD))Bsin(GQEFW(CFWD))Asin(GQEMW(CMWD)))QELW(CLWD))Bsin(GQEFW(CFWD))Asin(GQEMW(CMWD)))：F=GQEW(CF=GQEW(CWD)9

8、0WD)90：X=XZcosFX=XZcosF：Y=YZsinFY=YZsinF3.反算子程序(SUB2)T=G90T=G90：W=Abs((YV)cosT(XU)sinT)W=Abs((YV)cosT(XU)sinT)：Z=0Z=0：LblLbl0：ProgProg“SUB1““SUB1“：L=TQEW(CL=TQEW(CWD)WD)：Z=(JY)cosL(IX)sinLZ=(JY)cosL(IX)sinL：AbsZGoto1AbsZ

9、Goto1：≠W=WZ≠W=WZ：GotoGoto0Δ←┘0Δ←┘CASIO4800、4850坐標(biāo)正反計(jì)算程序坐標(biāo)正反計(jì)算程序XS=正算時(shí)，計(jì)算得出的所求點(diǎn)的X坐標(biāo)YS=正算時(shí)，計(jì)算得出的所求點(diǎn)的Y坐標(biāo)S=反算時(shí)，計(jì)算得出的所求點(diǎn)的里程Z=反算時(shí)，計(jì)算得出的所求點(diǎn)的邊距四、算例某匝道的由五段線(xiàn)元（直線(xiàn)完整緩和曲線(xiàn)圓曲線(xiàn)非完整緩和曲線(xiàn)直線(xiàn)）組成，各段線(xiàn)元的要素（起點(diǎn)里程S0、起點(diǎn)坐標(biāo)X0Y0、起點(diǎn)切線(xiàn)方位角F0、線(xiàn)元長(zhǎng)度LS、起點(diǎn)曲率半徑

10、R0、止點(diǎn)曲率半徑RN、線(xiàn)元左右偏標(biāo)志Q）如下：S0X0Y0F0LSR0RNQ500.00019942.83728343.5611251631.00269.2561E451E450769.25619787.34028563.3781251631.0037.4921E45221.751806.74819766.56628594.5741202554.07112.779221.75221.751919.52719736.07228701.8

11、93911730.6380.285221.759579.2281999.81219744.03828781.659804050.00100.0001E451E4501、正算（注意：略去計(jì)算方式及線(xiàn)元要素輸入，請(qǐng)自行根據(jù)所求點(diǎn)所在的線(xiàn)元輸入線(xiàn)元要素）S=700Z=5計(jì)算得XS=19831.41785YS=28509.72590S=700Z=0計(jì)算得XS=19827.33592YS=28506.83837S=700Z=5計(jì)算得XS=1982

12、3.25398YS=28503.95084S=780Z=5計(jì)算得XS=19785.25749YS=28575.02270S=780Z=0計(jì)算得XS=19781.15561YS=28572.16358S=780Z=5計(jì)算得XS=19777.05373YS=28569.30446S=870Z=5計(jì)算得XS=19747.53609YS=28654.13091S=870Z=0計(jì)算得XS=19742.68648YS=28652.91379S=87

13、0Z=5計(jì)算得XS=19737.83688YS=28651.69668S=940Z=5.123計(jì)算得XS=19741.59118YS=28722.05802S=940Z=0計(jì)算得XS=19736.47687YS=28722.35642S=940Z=3.009計(jì)算得XS=19733.47298YS=28722.531682、反算X=19831.418Y=28509.726計(jì)算得S=699.9999974Z=5.00018164X=1982

14、7.336Y=28506.838計(jì)算得S=699.9996493Z=0.000145136X=19823.25398Y=28503.95084計(jì)算得S=699.9999985Z=5.000003137X=19785.25749Y=28575.02270計(jì)算得S=780.0000035Z=5.000001663X=19781.15561Y=28572.16358計(jì)算得S=780.0000025Z=0.000002979X=19777.05