5800P线元计算.doc

FX-5800P曲线元计算 FX-5800P线元计算F[SIPSON]（主程序）Lbl 0（全局变量初始化）“F0:KZ=XY”:“F1:XY=KZ”?→MCls:90→C:“ANG=”?CIf M=0:Then Goto 1:Else If M=1:Then Goto 2:Else Goto 0:IfEnd: IfEndLbl 1（正算模式）Norm 1:DegCls:“Kp=”?SCls:“Zp=”?ZProg”DATP[XX]”Prog”TURN[YS]”Prog”HHB1”Cls:Fix 3“X[K]=”:Locate 6,2,XLocate 1,3,X “Y[K]=”:Locate 6,4,Y?Goto 1Lbl 2（反算模式）Norm 1:Deg“Xp=”?N: Cls:“Yp=”?E线路起点桩号如10660→SProg”DATP[XX]”Prog”TURN[YS]”Prog”HHB2”Cls:Fix 3“Kp=”:Locate 6,2,SLocate 1,3,“Zp=”:Locate 6,4,Z?Goto 2?HHB1（四等分复化辛卜生坐标正算子程序）Lbl 0（循环计算各节点值）For 1→M To 8G+180Q(M(S-O)÷(8P)+M^2(P-R)(S-O)^2÷(128PR(H-O)))÷π→ListX[M] Next List X[8]→T360(Frac((T+360)÷360)→TLbl 1（中边桩坐标计算）U+(S-O)(Cos(G)+4∑(Cos(ListX[2M+1]),M,0,3)+2∑(Cos(ListX[2M+2]),M,0,3)-Cos(List X[8]))÷24→XV+(S-O)(Sin(G)+4∑(Sin(ListX[2M+1]),M,0,3)+2∑(Sin(ListX[2M+2]),M,0,3)-Sin(List X[8]))÷24→YX+ZCos(T+C)→XY+ZSin(T+C)→YReturn?HHB2（四等分复化辛卜生坐标反算子程序）Lbl 0（反算数据初始化）G-(180-C)→FO+Abs((E-V)Cos(F)-(N-U)Sin(F))→S0→Z(置零)Lbl 1（反算桩号偏距）Prog”HHB1”180(T-(180-C))÷π→L(E-Y)Cos(L)-(N-X)Sin(L))→ZIf Z＜10^(-5):Then Goto 2:ElseS+Z→S:Prog”DATP[XX]”:Prog”TURN[YS]”:Goto 1IfEndLbl 2（计算偏距）(N-Y)÷Sin(T+C)→ZReturn?DATP[XX]（平面数据库子程序）If S≤Kei:Then [ [ Ksi,XSi , YSi , TSi , RSi ,Rei,Kei,Qi] ]→Mat A :Return:IfEndIf S≤Ke i+1:Then [ [ Ksi+1,Xsi+1 , Ysi+1 , Tsi+1, Rsi+1 ,Rei+1,Kei+1,Qi+1] ]→Mat A Return:IfEnd.............?TURN[YS]（将存在矩阵中的线元数据赋给正算标准变量）Mat A[1,1]→OMat A[1,2]→UMat A[1,3]→VMat A[1,4]→GMat A[1,5]→PMat A[1,6]→RMat A[1,7]→HMat A[1,8]→QReturn?关于程序的一些说明：1.程序的功能：实现线路全线任意角度正反算，反算无需输入近似桩号。2.程序计算核心是复化辛卜生近似求积分公式，相对于高斯-勒让德公式来说，是比较好理解的，即用小段抛物线拟合曲线。但是精度在分段数相同的条件下要差些，当然程序中4等分能满足一般放样要求。在5800P中，运行的速度也慢些。在此主要是比较，大家可以选择适合自己计算器的数学模型来编制程序。 FX-5800P曲线元计算之高斯-勒戎德版 FX-5800P之曲线元计算（高斯-勒让德版）F[QXY]（曲线元主程序）Lbl 0（线元要素输入）Norm 1:Deg“KS=”?0:Cls:“XS=”？U: Cls:“YS=”？V:Cls“TS=”?G:Cls:“RS=”?P:Cls:“Re=”?R:Cls“Ke=”?H:Cls:“[-1,0,1]PX=”?QLbl 1（坐标及切线方位角计算）Cls:Norm1”Kp=”?S: IF S=-1:Then Cls:Goto 0:IfEndCls:”Zp=”?ZCls:“ANG=”?CAbs(S-O)→W(P-R)