中边桩坐标计算程序:
ZBJS主程序
X:Y(这是立仪点坐标):D“X1”:E“Y1”(前一交点坐标):F“X2”:G“Y2”(曲线所在交点坐标):H“X3”:V“Y3”(后一交点坐标): L“LS”(缓和曲线长):W“J”(所在曲线交点桩号) ¿ I=D-F:J=E-G:Prog “C”(计算前一直线的极坐标):O(字母O)=J(前一直线方位角):I=H-F:J=V-G:Prog “C”:U=J(后一直线极坐标及方位角¿ M=O(字母O)-U(计算交点夹角):M<0(数字O)ÞM=M+360◣M<180ÞM=180-M:C=1: ¹ÞM=M-180:C=-1◣¿ P=L^2/(24R):Q=L/2-L^3/(240R^2):T=Q+(R+P)TAN(M/2) ◢ S=RMπ/180+L◢ M=MC:M“J0”◢ M=W-T:N=M:N“ZH”◢ N=M+L:N“HY”◢ N=M+S/2:N“QZ”◢ N=M+S-L:N“YH”◢ N=M+S:N“HZ”◢ (计算和显示曲线主点桩号,计算坐标时,可以不用显示,将黑三角形改为:就行了)Prog “B”
B(子程序、主要是进行桩号判断,进行分段处理)
Lbl 0:Fixm::I=K-W+T:J=Z:B(输入K:计算的桩号,Z:与中线的夹角、B:边桩的距离,左负右正,这是为计算边桩坐标用的)¿ I<S-LÞI≤0ÞM=I+BcosZ:N=BsinZ:Goto 2◣C=1ÞJ=180-Z◣I≤LÞ Prog “D”: ¹Þ Prog “F” ◣
Lbl 2:J=O(字母O): Prog “G”:Goto 0◣ I=S-I:I≤0ÞM=I-BcosZ:N=-BsinZ: ¹ÞC¹1ÞJ=180-Z◣Prog “D”:N=-N◣ J=U:Prog “G”:Goto 0
C(子程序、主要是进行坐标转换,转换为极坐标)
Pol(I,J):J<0ÞJ=J+360 ◣
D(子程序、主要计算缓和曲线上的X和Y)
A=RL:M=I-I^5/(40A^2)+I^9/(3456A^4):N=I^3/(6A)-I^7/(336A^3)+I^11/(42240A^5):
A=90I^2/(Aπ): Prog “E”
E(子程序,主要是计算边桩独立坐标)
A=A-J:M=M-CBcosA:N=CN-BsinA
F子程序(计算圆曲线上独立坐标)
A=90(2I-L)/(Rπ):M=RsinA+Q:N=P+R(1-cosA):Prog “E”
G(子程序,计算和显示坐标)
A=cosJ:J=sinJ:I=F+AT-MA+NJ:J=G+TJ-MJ-NA:M=I:N=J:I=I-X:J=J-Y¿
Prog “C”:M”X=” ◢ N”Y=” ◢ (计算和显示坐标) I◢ (放样时显示的距离) J◢(放样时显示的角度,在计算坐标时,可以不用显示,将黑三角形改为:就行了)
说明:1、由于显示的原因,所有“◣”均为空心的,不知怎么回事,全变成实心的了。
2、程序源于网络,只进行了部分修改和完善。
3、所有()中的均为说明,程序中不要出现。
把坐标一得X,Y分别输入到A1,B1;坐标二的X,Y分别输入到A2,B2距离,在A3输入公式 =SQRT((A2-A1)^2+(B2-B1)^2)
方位角,在A4输入公式 =ASIN((B2-B1)/A3)*180/PI()
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)