坐标计算4800源程序.docVIP

公路通用复化辛普森公式匝道点位坐标计算4800源程序 ?????? 公路通用复化辛普森公式匝道点位坐标计算4800源程序------------------杭浦高速临平互通---------------------　 本文利用的是计算公路匝道点位坐标的复化辛普森通用公式数学模型,集直线、圆曲线、回旋线通用，占字符内存较小，计算精度不限的程序一、运行变量名称说明：???? V=1、2分别进入坐标计算、桩号反算　　 K1、K2-------曲线起点、终点里程 　　 F0-----------曲线起点方位角　　 R1、 R2------曲线起点、终点半径 (ρ左-右+，0为直线)　　 X0、Y0-------曲线起点、终点坐标　　 M------------求和累积次数n的2倍(偶数)，精度迭代次数　　 K------------曲线待求点里程???? BP-----------求点左右偏距(左-右+) ???? ANG-------- -求点的右斜交角　　 X、Y---------曲线求得坐标???? FW-----------待求点的即时切线方位角???? XF、YF-------为需求桩号的点坐标???? DL、K+O、LP分别为桩号误差、求得桩号、左右偏距(左-右+)　　当曲线的设计半径较小时，为保证点位计算精度，M(即程序中n的2倍)的取值可适当的大些。M为偶数，直线时M=2即可，经计算M=16即可满足半径为60的小半径曲线精度。 二、曲线计算程序名:? Prog CURVEDefm 4VV=1 2Lbl 0:{KLW}Lbl 4QOPT:M0AB1C2D3E4FH5G6I7J8CR9=0=Prog MQ=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=Prog I△Q=8=Prog J△Q=9=Prog CR△AK1BK2CF0DR1ER2FX0GY0D≠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△KLBPWANGN=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)*(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]+ sinRN=M=Goto3：≠=Goto 2Lbl3X=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=Goto 6XX=◢YY=◢ RFW=R-360Intg(R/360◢Goto 0 Lbl 6???????? 进入桩号求算Pol(TXF-X,UYF-YO=Icos(J-RAbsO≤1e-4=ODLK=K+O◢???? OLP=Isin(J-R◢{TU}Goto 6:≠=K=K+O:L=0:Goto 4 三、数据文件：??? 线元要素数据文件每行为一个线元段，逐句执行赋值，直至不满足、运行完成。主线M:K=5865=A=5865:B=586.526:C=101`49`58`:D=0:E=0:F=3361488.152:G=449202.642KB=A=B:B=6235.526:C=101`49`58`:D=0:E=-3250:F=3361488.044:G=449203.157KB=A=B:B=7674.84:C=98`34`16.9`:D=-3250:E=-3250:F=3361419.068:G=449566.617KB=A=B:B=8044.84:C=73`11`49.3`:D=-3250:E=0:F=3361521.529:G=450990.516KB=A=B:B=8273.243:C=69`56`7.9`:D=0:E=0:F=3361641.834:G=451340.3558+273=8+100 S=1 2:{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) 匝ABK≥0:A= 0:B=716.396:C=85`46`15`:D=800:E=800:F=3361721.