4800改5800全线坐标正反算程序.doc

线路中边桩任意里程坐标正、反算改编程序（CASIO fx-4800改CASIO fx-5800）1. 正算子程序(SUB1)0.1739274226Z3：0.3260725774Z4：0.0694318442Z5：0.3300094782Z6：1-Z6Z7：1-Z5Z8：A+Z10 (Z3cos(F +57.2958 O Z5 Z10 (1/C+ Z5 Z10 Z9)+Z4cos(F +57.2958 O Z6 Z10 (1/C+ Z6 Z10 Z9)+Z4cos(F +57.2958 O Z7 Z10 (1/C+ Z7 Z10 Z9)+ Z3cos(F +57.2958 O Z8 Z10 (1/C+ Z8 Z10 Z9) Z11：B+ Z10 (Z3sin(F +57.2958 O Z5 Z10 (1/C+ Z5 Z10 Z9)+Z4sin(F +57.2958 O Z6 Z10 (1/C+ Z6 Z10 Z9)+Z4sin(F +57.2958 O Z7 Z10 (1/C+ Z7 Z10 Z9)+Z3sin(F +57.2958 O Z8 Z10 (1/C+ Z8 Z10 Z9) Z12：F +57.2958 O Z10 (1/C+ Z10 Z9)+90Z7：Z11+Scos(Z7-90+M) Z11：Z12+Ssin(Z7-90+ M) Z12：2. 反算子程序(SUB2)Abs(Z12-B)cos(F-90)-( Z11-A)sin(F-90)Z10：0S：Lbl 0：Prog SUB1：(F-90)+57.2958 O Z10 (1/C+ Z10 Z9)Z6：(J- Z12)cos Z6-(I-Z11)sin Z6S：If AbsZ XY：2.XY = SZ：“K=? K：“N=?N ：If K 下一线元起点里程：then 本线元起点里程Z1：本线元起点XA：本线元起点YB：本线元起算方位角F：本线元长度Z2：起点曲率半径C：终点曲率半径R：0或1、-1O： Goto0：Ifend： (第一线元数据要素)If K 下一线元起点里程：then 本线元起点里程Z1：本线元起点XA：本线元起点YB：本线元起算方位角F：本线元长度Z2：起点曲率半径C：终点曲率半径R：0或1、-1O： Goto0：Ifend ：(第二线元数据要素)If K 下一线元起点里程：then 本线元起点里程Z1：本线元起点XA：本线元起点YB：本线元起算方位角F：本线元长度Z2：起点曲率半径C：终点曲率半径R：0或1、-1O： Goto0：Ifend ：(第三线元数据要素)If：。 Goto0：Ifend：(第N-1线元数据要素)If：。 Goto0：Ifend：(第N线元数据要素)LB1 0：(C-R)(2 Z2 C R) Z9：If N=1：Then Goto 1：Else Goto 2：IfendLbl 1：“S=?S：“M =? M ： Abs(K- Z1) Z10：Prog SUB1：XS=：Z11YS=：Z12FS：Z7-90(可以删除) “XC=” ?P：“YC=” ?Q0I：0J：Pol（Z11-P，Z12-Q）：“I=” ：IIf J XY2. XY = SZ N ?选择计算方式，输入1表示进行由里程、边距计算坐标 ； 输入2表示由坐标反算里程和边距。S ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)M？为斜交角度，规定M为涵轴右侧方向与“线路前进方向切线”之间的夹角，当涵轴与线路正交时，T=90XS=正算时，计算得出的所求点的X坐标YS=正算时，计算得出的所求点的Y坐标FS=正算时，所求点对应的中线点的切线方位角“XC、YC、”为测站坐标“”、“”为放样距离和方位角 A ？线元起点的X坐标 B？线元起点的Y坐标 Z1 ？线元起点里程 F ？线元起点切线方位角Z2 ？线元长度 C ？线元起点曲率半径 R ？线元止点曲率半径 O ？线 元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0）K ？ 正算时所求点的里程X ？反算时所求点的X坐标Y ？反算时所求点的Y坐标 显示部分：K=反算时，计算得出的所求点的里程Z=反算时，计