坐标计算4800源程序

《坐标计算4800源程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、公路通用复化辛普森公式匝道点位坐标计算4800源程序      公路通用复化辛普森公式匝道点位坐标计算4800源程序------------------杭浦高速临平互通---------------------　本文利用的是计算公路匝道点位坐标的复化辛普森通用公式数学模型,集直线、圆曲线、回旋线通用，占字符内存较小，计算精度不限的程序一、运行变量名称说明：    V=1、2分别进入坐标计算、桩号反算　　K1、K2-------曲线起点、终点里程　　F0-----------曲线起点方位角　　R1、R2------曲

2、线起点、终点半径(ρ左-右+，0为直线)　　X0、Y0-------曲线起点、终点坐标　　M------------求和累积次数n的2倍(偶数)，精度迭代次数　　K------------曲线待求点里程    BP-----------求点左右偏距(左-右+)    ANG---------求点的右斜交角　　X、Y---------曲线求得坐标    FW-----------待求点的即时切线方位角    XF、YF-------为需求桩号的点坐标    DL、K+O、LP分别为桩号误差、求得桩号、左右偏距(左-右

3、+)　　当曲线的设计半径较小时，为保证点位计算精度，M(即程序中n的2倍)的取值可适当的大些。M为偶数，直线时M=2即可，经计算M=16即可满足半径为60的小半径曲线精度。二、曲线计算程序名: 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      进入坐标迭代计算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        进入桩号求

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三、数据文件：   线元要素数据文件每行为一个线元段，逐句执行赋值，直至不满足、运行完成。主线"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=