任意三角网平差计算 [pc1500]
由本人在网络上收集整理


2 LPRINT STR$ H;".P";STR$ PEEK (Q-3);"-P";STR$ PEEK (Q-2):RETURN 
3 R=196+2*(K+U):RETURN 
4 Z=0:GOSUB 1:LPRINT TAB 4;"Adjusted":LPRINT :RETURN 
5 N=PEEK (196+I):RETURN 
6 Q=2*(I-K)-1:RETURN 
7 X=X(O)-X(I):Y=Y(O)-Y(I):RETURN 
8 C=ATN (Y/X)+90*(2-SGN Y-SGN X*SGN Y):RETURN 
9 S=SQR (X*X+Y*Y):C(G)=V:E=SIN C/S*P:V=COS C/S*P
10 IF I>KLET C(Q-1)=E:C(Q)=-V:C(1,Q-1)=C(1,Q-1)+E:C(1,Q)=C(1,Q)-V
11 IF O>KLET C=2*(O-K)-1:C(C-1)=-E:C(C)=V:C(1,C-1)=C(1,C-1)-E:C(1,C)=C(1,C)+V
12 RETURN 
13 A=R+4*J:I=PEEK A:GOSUB 6:O=PEEK (A+1):GOSUB 7:GOSUB 8:RETURN 
16 V=1+G+5E-7:G$=STR$ (INT V-1):H$=STR$ (V-INT V)
17 IF LEN G$<3LET G$=X$+G$:GOTO 17
18 IF LEN H$<8LET H$=H$+"0":GOTO 18
19 LPRINT G$;X$;MID$ (H$,3,2);X$;MID$ (H$,5,2);".";MID$ (H$,7,2):RETURN 
23 R=PEEK (196+K+U+W):RETURN 
24 P=PEEK (196+2*(K+U)+R+H):RETURN 
25 FOR A=0TO G-1:C=C(M,A):IF C=0THEN 29
26 FOR B=ATO G-1:O=C(M,B)
27 IF OLET S=VAL I$(B)-B+A:E=INT (S/10):S=S-E*10:N(E,S)=N(E,S)+C*O*H
28 NEXT B:GOSUB 30
29 NEXT A:B=A
30 L(A)=L(A)+C(M,A)*C(M,B)*H:C(M,A)=0:RETURN 
31 FOR J=1TO F:Q=W+4*J:A=PEEK Q:B=PEEK (Q+1)
32 IF (I=AOR I=B)AND (O=AOR O=B)COLOR 1+2*(PEEK (Q+2)+PEEK (Q+3)=0):RETURN 
33 NEXT J:RETURN 
39 GOSUB 1:FOR I=1TO K+U:GOSUB 5:IF N=0LF 2:GOTO 42
40 LPRINT " 0 00 00.00":D(Z+1)=0:FOR J=2TO N:Y=Z+J:IF D(0)LET D(Y)=DMS D(Y)
41 G=D(Y):GOSUB 16:D(Y)=DEG G:NEXT J:Z=Z+N
42 GOSUB 1:NEXT I:Z=0:RETURN 
43 FOR I=ATO 2*U-1:E=VAL I$(I):F=INT (E/10):R=C(I)*(H>0)+(H=0):IF A=ITHEN 47
44 IF ILET B=I-E+VAL I$(I-1)+1
45 R=R*(H>0):IF I<=BTHEN 47
46 FOR J=A*(A>B)+B*(B>=A)TO I-1:Q=E-I+J:C=INT (Q/10):R=R-N(C,Q-10*C)*C(V,J):NEXT J
47 R=R/N(F,E-10*F):C(V,I)=R:IF V=0LET X=X+R*R:GOTO 49
48 Y=Y+R*R:Z=Z+R*C(I)
49 NEXT I:RETURN 
50 COLOR 2:LPRINT "M=";M;CHR$ 34:COLOR 0:M$=STR$ M:DIM I$(2*U-1)*4:IF D(0)=0RETURN 
55 X(0)=0:FOR I=0TO 2*U:C(I)=0:NEXT I:GOSUB 3:GOTO 330
60 "A"CLEAR :INPUT "Known Points=";K,"Unknown Points=";U:U$="########.###":X=K+U:X$=" "
62 F$=":From P":D$="Directions":B$="Equal a&b":S$="=":A$="Azimuths":Q$="Total "
68 V$="value":WAIT 0:INPUT "M=";M:PRINT Q$;D$;S$;:INPUT D:CLS :DIM X(X),Y(X),D(D),C(1,2*U)
70 O$="Observed ":Y$="###.##":FOR I=1TO X:GOSUB 3:PRINT "X";I;S$;:INPUT X(I):CLS 
72 PRINT "Y";I;S$;:INPUT Y(I):CLS :PRINT D$;S$;:INPUT N:CLS :POKE 196+I,N:IF N=0THEN 85
75 FOR J=1TO N:Y=Z+J:PRINT Y;":Aim P";:INPUT G:CLS :POKE 196+2*X+Y,G:IF J>1PRINT O$;V$;S$;:INPUT D(Y):CLS 
80 NEXT J:Z=Z+N
85 NEXT I:Z=0:Z$=":Points=":PRINT Q$;O$;A$;S$;:INPUT T:CLS :R=193+2*X+D:IF T=0THEN 100
90 DIM A(T):FOR I=1TO T:PRINT "a";I;F$;:INPUT A:CLS :INPUT "To P";B:PRINT O$;V$;S$;:INPUT A(I):CLS :
INPUT "Ma=";C
95 POKE R+4*I,A,B,INT C,100*(C-INT C):NEXT I
100 PRINT Q$;O$;"Sides=";:INPUT L:CLS :R=R+4*T:IF L=0END 
105 IF L>1PRINT B$;"(Y/N)";:INPUT W$:CLS :IF W$<>"Y"AND W$<>"N"THEN 105
110 DIM S(L):FOR I=1TO L:PRINT "S";I;F$;:INPUT C:CLS :INPUT "To P";E:PRINT O$;V$;S$;:INPUT S(I):CLS 
115 IF W$="N"OR I=1INPUT "Ms:a(mm)=";A,"b(ppm)=";B
120 POKE R+4*I,C,E,10*A,10*B:NEXT I:END 
160 "B"PRINT "M=";M;X$;:INPUT M
165 CLS :FOR I=1TO X:N=PEEK (196+I):PRINT "X";STR$ I;S$;X(I);X$;:INPUT X(I)
170 CLS :PRINT "Y";I;S$;Y(I);X$;:INPUT Y(I)
175 CLS :PRINT D$;S$;N;X$;:INPUT N:POKE 196+I,N
180 IF N=0THEN 200
185 CLS :FOR J=1TO N:Y=Z+J:GOSUB 3:R=R+Y:PRINT Y;":Aim P";PEEK R;X$;:INPUT E:POKE R,E
190 CLS :IF J>1PRINT "D";Y;S$;D(Y);X$;:INPUT D(Y)
195 CLS :NEXT J:Z=Z+N
200 NEXT I:Z=0:IF T=0THEN 230
205 FOR I=1TO T:P=R-3+4*I:PRINT "a";I;F$;PEEK P;X$;:INPUT E:POKE P,E
210 CLS :PRINT "To P";PEEK (P+1);X$;:INPUT E:POKE P+1,E
215 CLS :PRINT "a=";A(I);X$;:INPUT A(I)
220 CLS :PRINT "Ma=";PEEK (P+2)+PEEK (P+3)/100;X$;:INPUT C:POKE P+2,INT C,100*(C-INT C)
225 CLS :NEXT I
230 IF L=0END 
235 IF L>1PRINT B$;": ";W$;X$;:INPUT W$:IF W$<>"Y"AND W$<>"N"THEN 235
240 CLS :FOR I=1TO L:P=R-3+4*(T+I):PRINT "S";I;F$;PEEK P;X$;:INPUT E:POKE P,E
245 CLS :P=P+1:PRINT "To P";PEEK P;X$;:INPUT E:POKE P,E
250 CLS :PRINT "S=";S(I);X$;:INPUT S(I)
255 CLS :P=P+1:IF W$="Y"AND I>1POKE P,10*A,10*B:GOTO 270
260 PRINT "Ms:a(mm)=";PEEK P/10;X$;:INPUT A:POKE P,10*A
265 CLS :P=P+1:PRINT "b(ppm)=";PEEK P/10;X$;:INPUT B:POKE P,10*B
270 CLS :NEXT I:END 
310 "C"GOSUB 50:GOSUB 1:LPRINT "Point Direc":GOSUB 1:LPRINT :LPRINT TAB 7;"No. Aim":GOSUB 1
315 GOSUB 3:FOR I=1TO X:GOSUB 5:LPRINT "P";I;:IF N=0LF 2:GOTO 325
320 D(Z+1)=0:FOR J=1TO N:Y=Z+J:LPRINT TAB 6;Y;TAB 11;"P";STR$ PEEK (R+Y):NEXT J:POKE 196+X+I,Z:Z=Z+N
325 GOSUB 1:NEXT I:Z=0:GOSUB 1:LPRINT "tions":GOSUB 1:LPRINT "Reduced ";V$:GOSUB 39:LPRINT 
330 P=180/PI :R=R+D:IF T=0THEN 345
335 LPRINT O$;A$:FOR H=1TO T:Q=R+4*H:A=PEEK (Q-1)+PEEK Q/100
340 GOSUB 2:LPRINT " Ma=";A:COLOR 3*(A=0):G=A(H):LPRINT " a=";:GOSUB 16:COLOR 0:LPRINT :NEXT H:LPRINT 
345 R=R+4*T:LPRINT :IF L=0THEN 360
350 LPRINT O$;"Sides:":FOR H=1TO L:Q=R+4*H:A=PEEK (Q-1)/10:B=PEEK Q/10:GOSUB 2
355 LPRINT " a=";A;",b=";B:COLOR 3*(A+B=0):LPRINT " S=";S(H):COLOR 0:LPRINT :NEXT H:LPRINT 
360 GOSUB 3:FOR I=1TO K+U:GOSUB 5:IF N=0THEN 410
365 A=2*(I-K)-2:FOR J=1TO N:O=PEEK (R+Z+J):IF O<=KTHEN 405
370 B=2*(O-K)-2:IF I<=KTHEN 385
375 IF O<IAND VAL I$(A)<I-OLET I$(A)=STR$ (I-O)
380 IF O>IAND VAL I$(B)<O-ILET I$(B)=STR$ (O-I)
385 FOR H=JTO N:Q=PEEK (R+Z+H):IF Q<=KTHEN 400
390 IF Q<OAND VAL I$(B)<O-QLET I$(B)=STR$ (O-Q)
395 C=2*(Q-K)-2:IF Q>OAND VAL I$(C)<Q-OLET I$(C)=STR$ (Q-O)
400 NEXT H
405 NEXT J:X(0)=X(0)+1:Z=Z+N
410 NEXT I:Z=0:FOR I=1TO U:N=2*I-2:C=VAL I$(N)*2+1:IF I=1LET I$(0)="0":GOTO 420
415 I$(N)=STR$ (VAL I$(N-1)+C)
420 I$(N+1)=STR$ (VAL I$(N)+C+1):NEXT I:M=1:G=2*U:DIM N(INT (VAL I$(N+1)/10),9),L(G)
425 FOR I=1TO K+U:W=0:GOSUB 5:GOSUB 6:GOSUB 3:IF N=0THEN 445
430 FOR J=1TO N:O=PEEK (R+Z+J):GOSUB 7:GOSUB 8:F=C-D(Z+J):IF F<0LET F=F+360
435 IF M=1LET W=W+F/N:NEXT J:M=0:GOTO 43

[1] [2] 下一页