5800.docVIP

5800,坐标反算,缓和段反算用时小于5秒自编单交点坐标反算 技术支持 作者：涂振宇QQ：285507730 E-Mail: tuxhyu-yu@163.com 电话转载及引用核心技术请保留作者信息说明：单交点坐标反算桩号，缓和段反算速度小于5秒。 完成于2007-11-23晚含DJDZBFS????(主核心程序)SET????????(设置程序)TL???????? (点到直线，求垂足，及垂足至‘原点’距离)X　待求点N坐标，当输入-1值时可以进入要素输入程序Y　待求点E坐标Z[1] 起点NZ[2] 起点EZ[3] 交点NZ[4] 交点EZ[5] 终点NZ[6] 终点EZ[7]，L 　第一缓和段长Z[8]，R 　半径Z[9]，M 　第二缓和段长Z[10]，V 　交点桩号T第一切线长（SET程序内部计算，不需输入）U第二切线长（SET程序内部计算，不需输入）H第一切线方位角（SET程序内部计算，不需输入）C第二切线方位角（SET程序内部计算，不需输入）Z[11] 用于判断程序要是否被其它程序扰乱（SET程序内部计算，不需输入）先编SET程序再TL　最后编djdZBFSDJDZBFS????(主核心程序)Z[1]+Z[2] +Z[3] +Z[4] +Z[5] +Z[6] +L +R +M +V +T+U+H+C=Z[11]= Porg”SET”??;判断要素是否被其它程序扰乱，是则调用SET程序重设要素Lbi SRad:Norm 1　　设置计算机角度单位为弧度及小数点显示模式?XIF　X=-1：Then Porg”SET”:Goto S :IfEnd　　 ;如要想改变要素可以在X坐标输入”-1”?YZ[3]→NZ[4]→EH+π→WProg TLIf GT:Then V-G→G ;反算点在ZH（ZY）点前直线段上-I→IGoto E Else??　　;反算点在 在ZH（ZY）点后If L≠0:Then Rec(T,H)　　;有缓和段Z[3]-J→N??;计算ZH点X坐标Z[4]-I→E??;计算ZH点Y坐标Pol(L-L^(3)÷40÷R2,L2÷6÷R-L^(4)÷336÷R^(3))K0=-J→JH+J→JRec(I,J)N+J→N???;计算HY点X坐标E+I→E???;计算HY点Y坐标K0=-B→BH+B→WProg TLIf G0 :Then Rec(T,H)　　;反算点在第一缓和段上Z[3]-J→N??;计算ZH点X坐标Z[4]-I→E??;计算ZH点Y坐标Pol(L-L^(3)÷40÷R2,L2÷6÷R-L^(4)÷336÷R^(3))K0=-J→JH+J→WProg TLN→AE→SG→QLbi 1Pol(Q-Q^(5)÷40÷R2÷L2,Q^(3)÷6÷R÷L)K0=-J→JH+J→JRec(I,J)A+J→N??;计算假设垂足点X坐标S+I→E??;计算假设垂足点Y坐标Q2÷2÷R÷L→W　;计算假设垂足点切线方向K0=-W→W　　　H+W→WProg TLIf Abs(G)0.001 :Then??Q+G→Q:Goto 1?????? 如果反算桩号精度要求不高的话可以把 Abs(G)0.001??改为Abs(G)0.01 CM级别，这样反算桩号速度就快，但如设得过小的值 这可能导致无法反算成功。ElseV-T+Q+G→GGoto E:IfEnd??;IfEndIfEnd点不在HY点前继续执行以下语句Rec(U,C)Z[3]+J→N??;计算CZ(YZ)点X坐标Z[4]+I→E??;计算CZ(YZ)点Y坐标M=0=Goto 2 ;如果没有缓和段则跳过Pol(M-M^(3)÷40÷R2,M2÷6÷R-M^(4)÷336÷R^(3))K0=-J→JC+π+J→JRec(I,J)N+J→N??;YC点坐标E+I→ELbi 2K0=-O→OC+O→WProg TLIf G0 :Then???? ;反算点在圆曲线上π÷2→WK0=-W→WC+π+O+W→WIF Abs(W) π:Then IF W0:Then W-2π→W:Else W+2π→W:IfEnd:IfEndRec(R,W)N+JN??;圆心坐标E+I→EPol(E-Y,N-X)I-R→IK0=-I→IK0=W-J→JK0=J-W→JIF Abs(J) π:Then If J0:Then J-2π→J:Else J+2π→J:IfEnd:IfEndV-T+L+P