线元法线路坐标正反算程序.docVIP

??经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高，功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。功能简介及特点：　1、选用高斯-勒让德公式作计算内核，保证精度，模块化设计，便于扩充功能。　2、线元数据可自动从数据库调用，也可手工输入。　3、可管理多条线路，如里程不在线路或线元范围，将警告里程偏大、偏小。　4、边桩计算设计为导线式递推方式，可用于由一个中桩推出结构物所有角点坐标。　5、反算实现了智能化操作，只需输入线路号（或手工输线元资料）、坐标，不需近似里程，即可自动从起点向后开始试算出里程、位置，如对算出里程、位置表示怀疑，还可以让计算器从终点起再向前试算下一个可能的位置（匝道、回头曲线同一坐标可能会有一个以上结果）。第三次及以后试算才要求输入近似里程。　6、程序代码规范简洁，便于阅读、理解。 ?? 完整程序清单:ZFS %正反算主程序B=.1739274226:C=.5-B:Lbl 1:U0 ZS 1 FS=0=Prog ZS:≠U=1=ProgFS:≠Goto 1ZS %正算子程序{K}:ProgZZ:I=0:{I}:IL≠0=ProgWY:≠ProgZBFS %反算子程序 {KVW}:VXCWYC:Lbl 2:Prog ZZ:I=V-S:J=W-T:Pol(I,J:J=J-F:K=K+Rec(I,J:AbsI1m=ProgWZ:≠Goto 2ΔM=0:{M}:M0 NEXT=0=U=U+1:Goto 2:≠U=1ZZ %高斯法中桩子程序 (4节点)ProgXL:M=K-L:O=(P-R)÷2PQR:D=.0694318442:E=.3300094782:F=1:G=1-E:H=1-D:I=5:Lbl 1:C[I]=A+MrC[I](1÷P+OMC[I]:Dsz I:Goto 1:S=X+M(BcosD+CcosE+CcosG+BcosH:T=Y+M(BsinD+CsinE+CsinG+BsinHWY %外移点计算子程序Lbl 1:J=90:{J}:J=F+J:F=J:S=S+Rec(I,J:T=T+J:ProgZB:I=0:{I}:IL≠0=Goto 1WZ %位置显示子程序KJ:K:Pause 1:JZB %坐标显示子程序XY:S:Pause 1:TYC %异常处理子程序U=1=K=L:U=2ΔU=3=K=M:U=4ΔU=5={K}:U=4ΔKL=!:Z=1ΔKM=!:Z=1 DL %断链处理子程序DL:K=L:I0=K=L+Q-------------------------以上为程序运算部分,以下为数据库部分-------------------------------XL %线路数据库选择子程序Lbl 1:Z=0:N0 SD=0=Prog0N=1=Prog1△N=2=Prog2△...有几条线路仿上行格式输几行Z=1={NLXYOPQRK}:Goto 10 %手工输入子程序LK0XYAQLSPR0RRN:M=L+Q:ProgYC1 %线路一数据库子程序Lbl B:L=线路起点里程:M=线路终点里程:ProgYC:Z=1=Goto EΔQ=线元长:P=起点半径:R=终点半径:K≤L+Q=X=起点X坐标:Y=起点Y坐标:A=起点方位角:Goto EΔL=L+Q:......④Q=短链长:KL+Q=Prog DL:Goto BΔL=L+Q:Q=线元长:P=起点半径:R=终点半径:K≤L+Q=X=起点X坐标:Y=起点Y坐标:A=起点方位角:Goto EΔL=L+Q:......⑦Q=线元长:P=起点半径:R=终点半径:X=起点X坐标:Y=起点Y坐标:A=起点方位角:Lbl E2 %线路二数据库子程序输入要求和线路一相同。数据输入要求：是线路总体资料格式是标准线元资料格式，按里程顺序依次排断链资料格式最后一个线元资料格式　注：1、加下划线的行是程序名及功能说明　　2、ZZ程序第三行的r是计算器的弧度转度标志,相当于*180/π。????3、请注意区别数字0及字母O　　每个子程序第一行为程序名和功能简介用 %隔开,本程序语句在行未统一用冒号分隔，不须再加换行符，除≠前外所有冒号也可用换行代替。　　程序未尾的Δ（一个程序中Δ个数本应和=的个数相等）、：前的)]、小数点前的0可省略,不影响程序运行结果，省略乘号的连乘运算优先级高于除法。=========================================================================================程序使用说明：??　数据库准备　　将各条线路的“线元要素”依次输入