矩阵单交点曲线正反算程序.docVIP

矩阵单交点曲线正反算程序 主程序 Lbl A:1.NEW:2.OLD?M (选1调用矩阵数据进行初始化计算) M=2=Goto 0 (2执行上次数据) B-NUM.?N (选择MatB中的数据N为行数) 10→dimz (增加10个额外变量) Mat B[N,1]+(Mat B[N,2]i)→O (将起点坐标复数化并赋值于O) Mat B[N,3]→W (起点方位角赋值于W) cos(W)+sin(W)i→W (转化为向量) Mat B[N,4]→R (线路半径赋值于R) for 1→I to 5 (线元数据处理) Mat B[N,4+I]→Mat A[I,1] Mat B[N,4+I]→Mat A[I,2] Mat B[N,5+I]-Mat B[N,4+I]→V V→Mat A[I,3] O→Z[2I-1] W→Z[2I] (计算相关主点数据) I=1=progZZS (I为1时调用直线正算) If I=2:Then abs(R)abs(RV)^.5÷R→F:F→Mat A[I,4]: progHZS: IfEnd (为2时计算第一缓和参数并调用回旋正算) I=3=progYZS (为3调用圆正算) If I=4:Then -W→W:-F→F:F→Mat A[I,4]:progHDQ:IfEnd Next:Mat A[5,2]→Mat A[4,2] Z[9]→Z[7] Z[10]→Z[8] -Z[10]→Z[10] (为4时计算二缓起点要素，并调用相应子程序) Lbl 0:1.ZS:2.FS?N (正反算选择) N=0=Goto A (为0时返回到Lbl A) If N=1: Then Lbl 1:S?V:V=0=Goto 0 If VMat A[1,1]:Then cls:MIN◢ Goto 1:IfEnd If VMat A[5,1]+Mat A[5,3]:Then cls:MAX◢ Goto 1:IfEnd If V=Mat A[2,1]:Then 1→I:Goto C:IfEnd 6→I:do:I-1→I:lpwhile VMat A[I,1] Lbl C:Mat A[I,4]→F abs(V-Mat A[I,2])→V Z[2I-1]→O Z[2I]→W I=1=progZZS I=2=progHZS I=3=progYZS I=4=progHZS (根据输入里程调用相关子程序进行坐标正算) I=5=progZZS Deg:?X: ?Z:cos(X)+sin(X)i→X:cls MID:O◢ LEFT:O-WXZ◢ RIGHT:O+WXZ◢ Goto 1: IfEnd (输入夹角和边距计算边桩坐标) If N=2:Then do NUM.?I (确定反算所在线元，由1到5判断) I=0=Goto 0 (为零重新选择正反算) ?X:?Y:X+iY→V Mat A[I,2]→H Mat A[I,4]→F Z[2I-1]→O Z[2I]→W I=1=progZFS If I=2:Then 1→G:progHFS: IfEnd I=3=progYFS If I=4:Then -1→G:progHFS:IfEnd I=5=ProgZFS cls S:S◢ Z:Z◢ lpwhile 1:IfEnd (显示反算结果，) 直线正算子程序ZZS VW+O→O:return 直线反算ZFS (V-O)÷W+H→V:rep(V)→S:imp(V)→Z:return 圆正算YZS Rad:V÷R→V 2Rsin(.5V)(cos(.5V)+sin(.5V)i)W+O→O:cos(V)+sin(V)i→B:BW→W:return 圆反算YFS Rad:(V-O)÷W→V:V-Ri→A Ri+abs(R)A÷abs(A)→A arg(V)-arg(A)→B 2Rsin-1(.5abs(A)÷R)+H→S abs(V-A)abs(B)÷B→Z:return 回旋曲线正算HZS V+U→V:progHHHH:AW+O→O:BW→W:return 回旋曲线反算HFS (V-O)÷W→Q:rep(Q)→V do progHHHH rep((Q-A)÷B)+V→V lpwhile rep((Q-A)÷B)10^(-6) arg(Q)-arg(A)→B G(V-U)+H→S abs(Q-A)abs(B)÷B→Z:retu