任意三角网平差计算 [pc1500]

标签：工程测量规范,工程测量技术, 任意三角网平差计算 [pc1500],http://www.gong66.com
任意三角网平差计算 [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$="

[1] [2]  下一页

，任意三角网平差计算 [pc1500]