三角网平差计算

三角网平差计算 任意三角网平差计算 [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 EXT I:Z=0:IF T=0THEN 230 200 N 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 OIAND VAL I$(B)OAND VAL I$(C)KLET C(Q-1)=-COS C:C(Q)=-SIN C 470 S=SQR (X*X+Y*Y):C(G)=S-S(J):IF O>KLET F=2*(O-K)-1:C(F-1)=COS C:C(F)=SIN C 475 H=10:Z=PEEK (A+3)/1E3:IF ZLET H=(VAL M$/(PEEK (A+2)+S(J)*Z))^2/.1296 480 GOSUB 25:NEXT J 485 FOR H=0TO W:Z=VAL I$(H):Q=INT (Z/10):A=Z-VAL I$(H-SGN Z)+(H=0):FOR I=HTO W:Y=VAL I$(I)-I 490 X=VAL I$(I-SGN I)-Y-(I=0):IF H<=XTHEN 525 495 E=Y+H:F=INT (E/10):M=E-10*F:IF A=1THEN 515 500 FOR J=H-A+1TO H-1:IF J<=XTHEN 510 505 B=Z-H+J:C=INT (B/10):V=Y+J:S=INT (V/10):N(F,M)=N(F,M)-N(C,B-10*C)*N(S,V-10*S) 510 NEXT J 515 IF I=HLET O=SQR N(Q,Z-10*Q) 520 N(F,M)=N(F,M)/O 525 NEXT I:IF A>1FOR J=H-A+1TO H-1:B=Z-H+J:F=INT (B/10):L(H)=L(H)-N(F,B-10*F)*L(J):NEXT J 530 L(H)=L(H)/O:L(G)=L(G)-L(H)^2:NEXT H:J=K+U 535 FOR H=WTO 0STEP -1:X=-L(H):FOR I=H+1TO W:B=VAL I$(H) 540 IF IVAL I$(I-1)LET E=INT (F/10):X=X-N(E,F-10*E)*L(I) 545 NEXT I:Q=INT (B/10):X=X/N(Q,B-10*Q):L(H)=X:IF INT (H/2)0AND Y<0)):IF I<=KAND O<=KCOLOR 3:GOTO 590 585 IF TGOSUB 31 590 GOSUB 16:COLOR 0:IF H=1LET E=G:NEXT H 595 C=DEG G-DEG E:IF C<0LET C=C+360 600 D(Z+H)=C:NEXT H:Z=Z+N 605 GOSUB 1:NEXT I:GOSUB 4:LPRINT TAB 8;"S":LPRINT :GOSUB 1:F=L 610 W=W+4*T:FOR I=1TO K+U:GOSUB 5:IF N=0LF 2:GOTO 630 615 FOR H=1TO N:R=Z:GOSUB 24:O=P:GOSUB 7:IF I<=KAND O<=KCOLOR 3:GOTO 625 620 IF LGOSUB 31 625 LPRINT TAB 1;SQR (X*X+Y*Y)+5E-4:COLOR 0:NEXT H:Z=Z+N 630 GOSUB 1:NEXT I:GOSUB 4:LPRINT TAB 2;D$:LPRINT :GOSUB 1:USING Y$:FOR I=1TO K+U:GOSUB 5:IF N=0LF 2:GOTO 645 640 LPRINT " 0 00 00.00":FOR J=2TO N:G=DMS D(Z+J):GOSUB 16:NEXT J:Z=Z+N 645 GOSUB 1:NEXT I:Z=0:GOSUB 1:LF 2:LPRINT " Ms/S":LPRINT :GOSUB 1:V=0:USING :GOSUB 3:G=R-3+D+4*T 650 FOR O=1TO K+U:IF PEEK (196+O)=0LF 2:GOTO 715 655 FOR H=1TO PEEK (196+O):R=Z:GOSUB 24:IF P<=KAND O<=KLPRINT :GOTO 710 660 IF OKLET A=2*(P-K)-2:C(A)=X/S:C(A+1)=Y/S 700 IF O>KLET C=2*(O-K)-2:C(C)=-X/S:C(C+1)=-Y/S:IF O0):TEXT :LPRINT :RETURN