basic正反算.txt

用BASIC语言实现54坐标系参数下的的正反算。对于84坐标系

5：REM "GSTY--ZBHDJS:ZHS-FSHD"
7:"A"CLEAR
10:R=6367558.49686:RH=180/π
20:INPUT "K=";K,"DS=";DS
30:ON KGOTO 40,110,150
40:INPUT "L0=";L0,"D6=";D6,"A=";A
60:FOR I=1 TO DS
70:INPUT "B=";B:INPUT "L-L0=";L:PRINT "B=";B:PRINT "L=";DMS(DEG L+DEG L0)
80:GOSUB 370
90:GOSUB 410
100:NEXT I:END
110:INPUT "L0=";L0,"D6=";D6,"A=";A
120:GOSUB 430
140:END
150:INPUT "L0=";L0,"LT=";LT,"D6=";D6,"A1=";A1,"A2=";A2
161:FOR I=1 TO DS:A=A1
162:INPUT "X=";X:X1=X
163:INPUT "Y=";Y:Y1=Y 
164:PRINT "X1=";X:PRINT "Y1=";Y:GOSUB 410:L=DMS(DEG L+DEG L0-DEG LT)
170:LL=L0:L0=LT:LT=LL
185:GOSUB 370
195:A=A2
200:GOSUB 410:L=DMS(DEG L+DEG L0-DEG LT):LL=L0:L0=LT:LT=LL:GOSUB 370
210:NEXT I:END
230:T=TAN BF:C=COS BF:E=0.006738525415*C*C:H=T*T:Q=1+E:N=6399698.9018/sqr(Q):RETURN
240:K1=32005.78006:K2=133.92133:K3=0.70310:W=W*C:M=W*W:U=T*C:V=U*U
250:K4=(K1+V*(K2+V*K3))
260:X1=R*B/RH-U*C*K4+((((H-58)*H+61)*M/30+(4*E+5)*Q-H)*M/12+1)*N*T*M/2
270:Y1=((((H-18)*H-(58*H-14)*E+5)*M/20+Q-H)*M/6+1)*N*W
280:G=DMS ((T*W*(1+M*((Q+E)*Q/3+M*(2-H)/15)))*RH:RETURN
290:IF D6=6LET DH=((L0+3)/6)*1E6+5E5
293:IF D6=3LET DH=(L0/3)*1E6+5E5+5E5
296:IF D6=1LET DH=0
298:RETURN
300:BF=X1/R:BF=BF*RH:U=SIN BF:V=U*U
310:BF=BF+U*COS BF*(5.051773759E-3-V*(2.983718E-5-V*2.38209E-7))*rh
320:GOSUB 230:N=Y1/N:M=N*N
330:B=BF-((((45*H+90)*H+61)*M/30-(3-9*E)*H-5-E)*M/12+1)*M*T*Q/2*RH
340:B=DMS B
350:L=((((24*H+28)*H+(8*H+6)*E+5)*M/20-2*H-Q)*M/6+1)*N/C
360:L=L*RH:L=DMS L
365:RETURN
370:B=DEG B:L=DEG L:L0=DEG  L0
380:W=L/RH:BF=B:GOSUB 230:GOSUB 240
390:L0=DMS L0:GOSUB 290:Y1=Y1+DH
400:PRINT "X=";X1:PRINT "Y=";Y1:PRINT "G=";G:RETURN
410:Y1=Y1-A*1E5:GOSUB 300:PRINT "B=";B:PRINT "L=";DMS (DEG L+DEG L0):RETURN
430:FOR I=1 TO DS
440:INPUT "X=";X:INPUT "Y=";Y
445:X1=X:Y1=Y
446:PRINT "X=";X:PRINT "Y=":Y 
450:GOSUB 410
460:GOSUB 370
470:NEXT I:RETURN