线路中桩坐标和切线方位角计算.doc

线路中桩坐标和切线方位角计算一友人整理了线路中桩坐标和切线方位角计算的两种方法，拿出来与大家探讨，目的是对您的测量工作能有所帮助。这两种方法尚未经过实践验证，如果您发现它们有错误或缺陷，别忘了告诉我。方法一：线路中桩坐标和切线方位角计算（CASIO-4800P程序）适用于直线、第一缓和曲线（ZH-HY）、园曲线、卵形曲线（YH-HY，可以从大半径到小半径，也可以从小半径到大半径）、第二缓和曲线（YH-HZ）。输入变量：A=起点桩号；B=终点桩号；D=起点切线方位角；F=起点X坐标；H=起点Y坐标；K=起点曲率；R=终点曲率（曲率为1半径，直线的曲率为0）；M=判断因子（线路左转M=1；线路右转M=1；直线M=0）；C=桩号（起点和终点之间的任意桩号，ACB）。输出：X，Y=桩号为C点的坐标；P=桩号为C点的切线方位角（单位：度）。程序开始：ABDFHKRMLbl0:CK-R=0=S=C-A:Goto5:=Goto2Lbl2:G=ABS(B-A)(R-K)K=0=L=0:S=C-A:Goto1R=0=L=B-A:S=L-C+A:Goto1L=KGKS=L+C-A:=S=L-C+ALbl1:N=S-L-(S5-L5)/(40G2)+(S9-L9)/(3456G4)-(S13-L13)/(599040G6)+(S17-L17)/175472640G8)-(S21-L21)/(78033715200G10)E=(S3-L3)/(6G)-(S7-L7)/(336G3)+(S11-L11)/(42240G5)-(S15-L15)/(9676800G7)+(S19-L19)/(3530096640G9)-(S23-L23)/(1880240947200G11)KR=N=-N:Z=-1:=N=N：Z=1T=D-90MZL2(G)X=F+NcosT-MEsinTY=H+NsinT+MEcosTP=D+90MAbs(S2-L2)/(G)PP=P+360P360=P=P-360P=PGoto0Lbl5:K=0=Goto6:=Goto7Lbl6:X=F+ScosDY=H+SsinDP=DGoto0Lbl7:X=F+2sin(90SR)cos(D+90MSR)RY=H+2sin(90SR)sin(D+90MSR)RP=D+180SRPP=P+360P360=P=P-360P=PGoto0程序结束 线路中桩坐标和切线方位角计算公式（方法一）A=起点桩号，B=终点桩号，C=AB上任意点桩号，D=起点切线方位角，X0=起点X坐标，Y0=起点Y坐标，M=左转为-1；右转为1；直线为0，K=起点曲率，R=终点曲率。直线：R，S=C-AX=X0+ScosDY=Y0+SsinDP=D园曲线：，S=C-AX=X0+2sin（90SR/）.cos（D+90MSR/）/RY=Y0+2sin（90SR/）.sin（D+90MSR/）/RP=D+180MSR/缓和曲线：曲线参数G=Abs（B-A）/（R-K）第一缓和曲线：，L=0，S=C-A第二缓和曲线：，L=B-A，S=L-（C-A）=B-C是卵形曲线：，0，L=KGR，S=L+(C-A)，Z=1R，S=L-(C-A)，Z=-1（S-L）（S5-L5）40G（S9-L9）3456G-（S13-L13）599040G6（S17-L17）175472640G8-（S21-L21）78033715200G10（SL）G（SL）336G（SL）42240G-（S15-L15）9676800G7（S19-L19）3530096640G9-（S23-L23）1880240947200G11KR，N=-NT=D-90MZL2/GX=X0+NcosT-MEsinTY=Y0+NsinT+MEcosTP=D+90MAbs（S2-L2）/G方法二：缓和曲线中桩坐标和切线方位角计算程序（PC-E500）适用于第一缓和曲线、第二缓和曲线、卵型曲线（R1R2或R2R1）。该方法无须知道起点切线方位角，但要已知终点坐标。如果想用它求起点切线方位角，只要输入的桩号C比A大1毫米就可以了。1.输入信息A=起点桩号；B=终点桩号；XA，YA=起点坐标；XB，YB=终点坐标；R1=起点半径；R2=终点半径；K=左偏为-1，右偏为1；C=AB之间任意点桩号，ACB。直缓（缓直）点半径输1E50；Y$输“Y”，结束，回车继续下一个桩号的计算。2.输出信息XP，YP=桩号为C点的坐标；AT=桩号为C点的切线方位角（度.分秒）。3.源程序：10:INPUTA=;A:INPUTB=;B:INPUTXA=;XA:INPUTYA=;YA:INPUTXB=;XB:INPUTYB=;YB15:INPUTR1=;R1:INPUTR2=;R2:INPUTK=;K20:G=ABS(R1*R2*(B-A)/(R1-R2):S=SQR(YB-YA)2+(XB-XA)2)21:P0=ASN(YB-YA)/S)23:IF(XB-XA)0ANDP00THENP0=P0+36025:IFR1R2THENZ=-1:M=G/R2+B-A:GOTO3026:Z=1:M=G/R130:L=B-A:GOSUB104031:P1=PP35:INPUTC=;C:L=C-A40:GOSUB104050:P2=P1-PP:P3=P0-K*P260:XP=XA+DP*COSP3:YP=YA+DP*SINP370:P4=180*(M*L+Z*L2/2)/G/PI:AT=P3+K*(P4-PP)80:PRINTXP=;XP:PRINTYP=;YP:PRINTAT=;DMSAT100:INPUTEND?;Y$110:IFY$=YTHENEND120:GOTO351040:X=L-L3*(M*M/3+Z*M*L/4+L*L/20)/(2*G2)+L5*(M4/5+Z*M3*L/3+3*M2*L2/14+Z*M*L3/16+L4/144)/24/G41050:Y1=L2*(M+Z*L/3)/2/G-L4*(M3/2+3*Z*M2*L/5+M*L2/4+Z*L3/28)/12/G31