Lbl A:{SDO}:S"K=":D"DP=":O"ANG=":Z"NO.=1,2,SD=3"=1=>Goto X:≠>Z=2=>Goto Y:≠>Goto W△△
Lbl X:S≤593.652=>J=1:Prog "DAT5":Goto 0△
S≤852.072=>J=2:Prog "DAT5":Goto 0△
S≤1033.94=>J=3:Prog "DAT5":Goto 0△
S≤1170.162=>J=4:Prog "DAT5":Goto 0△
………………………………
S≤6275.818=>J=18:Prog "DAT5":Goto 0△
Lbl Y:S≤997.192=>J=1:Prog "DAT6":Goto 0△
S≤2208.389=>J=2:Prog "DAT6":Goto 0△
S≤5565.209=>J=3:Prog "DAT6":Goto 0△
S≤7113.552=>J=4:Prog "DAT6":Goto 0△
………………………………
S≤17189.541=>J=8:Prog "DAT6":Goto 0△
Lbl W:A"XJ"B"YJ"L"LS1"M"LS2"F"F0"E"PJ"K"JD"R
Lbl 0:P=L2÷24÷R-L^4÷2688÷R^3:Q=L÷2-L^3÷240÷R2+L^5÷34560÷R^4:G=πRAbs E÷180+L:I=M2÷24÷R-M^4÷2688÷R^3:C=(R+P)÷cos tan-1(((R+I)÷(R+P)-cos Abs E)÷sin Abs E)-R:C:C=M÷2-M^3÷240÷R2+M^5÷34560÷R^4:T=(R+I-(R+P) cos AbsE)÷sin Abs E+Q:T:N=(R+P-(R+I) cos Abs E)÷sin Abs E+C:N:H=K-T:H:I=H+L:I:I=H+G-M:I:I=H+G:I
Lbl 1:E<0=>C=-1:≠>C=1△S≤H=>Goto 3△
S≤H+L=>Goto 4△
S≤H+G-M=>Goto 5△
S≤H+G=>Goto 6:≠>Goto 7△
Lbl 2:Fix3:"X=":X=X:Pause 0:"Y=":Y=Y▲
"FW=":Prog "FW4":J→DMS▲ Goto A
Lbl 3:J=F:X=A-T cos J+(S-H) cos J+D cos (J+O):Y=B-T sin J+(S-H)sin J+Dsin (J+O):Goto 2
Lbl 4:V=S-H:Prog "H1.S":J=F:I=T:Prog "H2.S":J=F+90 CV2÷(πRL):Prog "FW4":X=X+Dcos (J+O):Y=Y+D sin (J+O):Goto 2
Lbl 5:V=2(S-H)-L:W=90 V÷π÷R:U=R sin W+Q:W=R(1-cos W)+p:J=F:I=T:Prog "H2.S":J=F+90 CV÷π÷R:Prog "FW4":X=X+D cos (J+O):Y=Y+D sin (J+O):Goto 2
Lbl 6:V=H+G-S:L=M:Prog "H1.S":U= -U:J=F+E:I= -N :Prog "H2.S":J=F+E-90 CV2÷(πRM):Prog "FW4":X=X+D cos (J+O):Y=Y+D sin (J+O):Goto 2
Lbl 7:V=S-H-G:J=F+E:X=A+(N+V) cos J+D cos (J+O):Y=B+(N+V) sin J+D sin (J+O):Goto 2
二、子程序:
1、文件名:H1.S
Q=6:G=1
U=∑((-1)^(G-1)×V ^ (4G-3)÷((2G-2)!×(4G-3)×2 ^ (2G-2)×(RL)^(2G-2)),G,1,Q)
W=∑((-1)^(G-1)×V ^ (4G-1)÷((2G-1)!×(4G-1)×2 ^ (2G-1)×(RL)^(2G-1)),G,1,Q)
2、文件名:H2.S
X=A-I cos J+U cos J - CW sin J:Y=B-I sin J+U sin J + CW cos J
3、文件名:FW4
J<0=>J=J+360:≠>J>360=>J=J-360:≠>J=J
4、数据库1:DAT5
J=1=>A=6136.073:B=3300.4972:L=35:M=35:F=171°36′6.9″:E=74°58′29″:K=529.067:R=87.72△
………………………………………
J=N=>A=****:B=****:L=****:M=***:F=****:E=****:K=****:R=****△
5、数据库2:DAT6
J=1=>A=63559.945:B=531452.29:L=0:M=0:F=286°16′56.9″:E=10°46′40.7″:K=528.306:R=5000△
………………………………………
J=N=>A=****:B=****:L=****:M=***:F=****:E=****:K=****:R=****△
三、说明:
1、输入部分:
"NO.=1,2,SD=3"—输入1进行第一条线路正算,输入2进行第二条线路正算,输入3进行人工正算
XJ—交点X坐标(利用数据库全线贯通计算时不出现提问,以下类同)
YJ—交点X坐标
LS1—前缓长(如仅为圆曲线时,输入零)
LS2—后缓长(如仅为圆曲线时,输入零)
F0—起始方位角
PJ—曲线偏角(左负右正)
JD—交点里程
R—圆曲线半径
2、显示部分:
X= ——所求点X坐标
Y= ——所求点Y坐标
FW=——所求点切线方位角
3、功能:本程序适用于对称或非对称的完整曲线的正算