精确坐标与切线方位角5800.doc

线路中桩坐标和切线方位角计算（CASIO-5800P程序） 适用于直线、第一缓和曲线（ZH-HY）、圆曲线、卵形曲线（YH-HY，可以从大半径到小半径，也可以从小半径到大半径）、第二缓和曲线（YH-HZ）。输入变量：A=起点桩号；B=终点桩号；D=起点切线方位角；F=起点X坐标；H=起点Y坐标；K=起点曲率；R=终点曲率（曲率为 1÷半径，直线的曲率为0）；I=判断因子（线路左转I=－1；线路右转I=1； 直线I=0）；C=桩号（起点和终点之间的任意桩号，A≦C≦B）。输出：X，Y=桩号为C点的坐标； P=桩号为C 点的切线方位角（单位：度）。UV为C点任意边距角度的XY坐标,O为边桩点到中线C点的距离(路左为负,路右为正)，Q为边桩点到路中线C点的偏角（顺时针）程序开始：？A：？B：？D：？F：？H：？K：？R：？I：Lbl 0:？C：If K-R=0：Then S=C-A:Goto 5:Else Goto 2 ：IfEnd Lbl 2: ABS((B-A)÷(R-K))→G K=0=L=0: C-A→S:Goto 1 R=0=B-A→L:L-C+A→S:Goto 1 KG→L:If KR:Then S=L+C-A: Else L-C+A→S: IfEnd: Goto 1 Lbl1: S-L-(S^5-L^5)/(40G^2)+(S^9-L^9)/(3456G^4)-(S^13-L^13)/(599040G^6)+(S^17-L^17)/175472640G^8)-(S^21-L^21)/(78033715200G^10)→N (S^3-L^3)/(6G)-(S^7-L^7)/(336G^3)+(S^11-L^11)/(42240G^5)-(S^15-L^15)/(9676800G^7)+(S^19-L^19)/(3530096640G^9)-(S^23-L^23)/(1880240947200G^11)→E If KR:Then N=-N:Z=-1: Else N=N：Z=1 : IfEnd:”X=”D-90IZL^2÷(πG)→TF+NcosT-IEsinT→ X◢“Y=“：H+NsinT+IEcosT→Y ◢”P=：“D+90I Abs(S^2-L^2)/(πG) →PP0=P=P+360 P360=P=P-360 P=P◢Goto 8Lbl 5:If K=0：Then Goto 6: Else Goto 7 ：IfEndLbl 6:“X=“： F+ScosD→X ◢? ? ? ”Y=“：H+SsinD→Y ◢? ? ”P=“?：D→P◢Goto 8Lbl 7:”X=“：F+2sin(90SR÷π)×cos(D+90ISR÷π) ÷R→X ◢ “Y=”：H+2sin(90SR÷π)×sin(D+90ISR÷π) ÷R→ Y◢? ? ?“ P=”：D+180ISR÷π→PP0=P=P+360⊿P360=P=P-360⊿P=P◢Goto 8 Lbl 8: ？O:？Q:“U=”:X+O*cos(P+Q)→U ◢ “V=”：Y+O*sin(P+Q)→V ◢ Goto 0 程序结束