資源描述：
《坐標(biāo)計(jì)算4800源程序》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、公路通用復(fù)化辛普森公式匝道點(diǎn)位坐標(biāo)計(jì)算4800源程序??????公路通用復(fù)化辛普森公式匝道點(diǎn)位坐標(biāo)計(jì)算4800源程序------------------杭浦高速臨平互通--------------------- 本文利用的是計(jì)算公路匝道點(diǎn)位坐標(biāo)的復(fù)化辛普森通用公式數(shù)學(xué)模型,集直線、圓曲線、回旋線通用，占字符內(nèi)存較小，計(jì)算精度不限的程序一、運(yùn)行變量名稱說明：????V=1、2分別進(jìn)入坐標(biāo)計(jì)算、樁號(hào)反算 K1、K2-------曲線起點(diǎn)、終點(diǎn)里程 F0-----------曲線起點(diǎn)方位角 R1、R2------曲
2、線起點(diǎn)、終點(diǎn)半徑(ρ左-右+，0為直線) X0、Y0-------曲線起點(diǎn)、終點(diǎn)坐標(biāo) M------------求和累積次數(shù)n的2倍(偶數(shù))，精度迭代次數(shù) K------------曲線待求點(diǎn)里程????BP-----------求點(diǎn)左右偏距(左-右+)????ANG---------求點(diǎn)的右斜交角 X、Y---------曲線求得坐標(biāo)????FW-----------待求點(diǎn)的即時(shí)切線方位角????XF、YF-------為需求樁號(hào)的點(diǎn)坐標(biāo)????DL、K+O、LP分別為樁號(hào)誤差、求得樁號(hào)、左右偏距(左-右
3、+) 當(dāng)曲線的設(shè)計(jì)半徑較小時(shí)，為保證點(diǎn)位計(jì)算精度，M(即程序中n的2倍)的取值可適當(dāng)?shù)拇笮?。M為偶數(shù)，直線時(shí)M=2即可，經(jīng)計(jì)算M=16即可滿足半徑為60的小半徑曲線精度。二、曲線計(jì)算程序名:?Prog"CURVE"Defm4V"V=12"Lbl0:{KLW}Lbl4Q"OPT:M0AB1C2D3E4FH5G6I7J8CR9"=0=>Prog"M"△Q=1=>Prog"AB"△Q=2=>Prog"C"△Q=3=>Prog"D"△Q=4=>Prog"E"△Q=5=>Prog"FH"△Q=6=>Prog"G"△Q=7=
4、>Prog"I"△Q=8=>Prog"J"△Q=9=>Prog"CR"△A"K1"B"K2"C"F0"D"R1"E"R2"F"X0"G"Y0"D≠0=>I=1/D:≠=>I=D△E≠0=>J=1/E:≠=>J=E△AbsD+AbsE=0=>M=2:≠=>M=16△V=2=>L=0:W=90△KL"BP"W"ANG"N=0：Z[1]=0：Z[2]=0：Z[3]=0：Z[4]=0??????進(jìn)入坐標(biāo)迭代計(jì)算Lbl2N=N+1：H=2(K-A)/M：R=NH/2+A：R=C+180/π*(I+(J-I)/2(B-A)*(
5、R-A))*(R-A)Int(N/2)=N/2=>Z[1]=Z[1]+cosR：Z[2]=Z[2]+sinR：≠=>Z[3]=Z[3]+cosR：Z[4]=Z[4]+sinR△N=M=>Goto3：≠=>Goto2Lbl3X=F+H/6*(cosC+4Z[3]+2Z[1]-cosR)+Lcos(R+W)Y=G+H/6*(sinC+4Z[4]+2Z[2]-sinR)+Lsin(R+W)V=2=>Goto6△X"X="◢Y"Y="◢R"FW"=R-360Intg(R/360◢Goto0Lbl6????????進(jìn)入樁號(hào)求
6、算Pol(T"XF"-X,U"YF"-YO=Icos(J-RAbsO≤1e-4=>O"DL"◢K=K+O◢????O"LP"=Isin(J-R◢{TU}Goto6:≠=>K=K+O:L=0:Goto4三、數(shù)據(jù)文件：???線元要素?cái)?shù)據(jù)文件每行為一個(gè)線元段，逐句執(zhí)行賦值，直至不滿足、運(yùn)行完成。主線"M":K>=5865=>A=5865:B=586.526:C=101`49`58`:D=0:E=0:F=3361488.152:G=449202.642K>B=>A=B:B=6235.526:C=101`49`58`:D=0
7、:E=-3250:F=3361488.044:G=449203.157K>B=>A=B:B=7674.84:C=98`34`16.9`:D=-3250:E=-3250:F=3361419.068:G=449566.617K>B=>A=B:B=8044.84:C=73`11`49.3`:D=-3250:E=0:F=3361521.529:G=450990.516K>B=>A=B:B=8273.243:C=69`56`7.9`:D=0:E=0:F=3361641.834:G=451340.355"8+273=8+100
8、S=12":{S}S=2=>A=8100:B=9610:C=C:D=0:E=0:F=3361720.194:G=451554.895注：主線中S=1、2為斷鏈1、2段(如K8+273.243=K8+100)匝"AB"K≥0:A=0:B=716.396:C=85`46`15`:D=800:E=800:F=3361721.181:G=450482.643K>B=
