PC1500任意高程网平差计算

2&nbspY=VAL&nbspP$(I):Z=VAL&nbspP$(1,I):S=S(I):RETURN&nbsp
3&nbspM=INT&nbsp(X/10):N=X-10*M:IF&nbspS=0RETURN&nbsp
4&nbspN(M,N)=N(M,N)+J/S:RETURN&nbsp
5&nbsp"A"CLEAR&nbsp:WAIT&nbsp0:USING&nbsp:U$="#####.###":INPUT&nbsp"Network&nbspKind?(L,T)";N$
10&nbspB$="&nbspPoints":M$="cm":IF&nbspN$="L"LET&nbspB$="&nbspBM":U$=U$+"#":M$="mm"
20&nbspPRINT&nbsp"Known";B$;"=";:INPUT&nbspB:CLS&nbsp:PRINT&nbsp"Unknown";B$;"=";:INPUT&nbspU:CLS&nbsp:INPUT&nbsp"Tota1&nbspdh=";T
25&nbspDIM&nbspH(B+U),D(T),S(T),P$(1,T)*2:U=U-1:DIM&nbspL(U),D$(U)*4:IF&nbspN$="T"LET&nbspB$="&nbspP"
30&nbspFOR&nbspI=1TO&nbspB:PRINT&nbsp"H";I;"=";:INPUT&nbspH(I):CLS&nbsp:NEXT&nbspI
35&nbspFOR&nbspI=1TO&nbspT:PRINT&nbsp"dh";I;"=";:INPUT&nbspD(I):CLS&nbsp:PRINT&nbsp"S";I;"(km)=";:INPUT&nbspS(I):CLS&nbsp
40&nbspPRINT&nbsp"From";B$;:INPUT&nbspP$(I):CLS&nbsp:PRINT&nbsp"To";B$;:INPUT&nbspP$(1,I):CLS&nbsp:NEXT&nbspI:END&nbsp
45&nbsp"B"FOR&nbspI=1TO&nbspB:PRINT&nbsp"H";I;"=";H(I);"&nbsp";:INPUT&nbspH(I)
50&nbspCLS&nbsp:NEXT&nbspI:FOR&nbspI=1TO&nbspT:PRINT&nbsp"dh";I;"=";D(I);"&nbsp";:INPUT&nbspD(I)
55&nbspCLS&nbsp:PRINT&nbsp"S";I;"=";S(I);"&nbsp";:INPUT&nbspS(I)
60&nbspCLS&nbsp:PRINT&nbsp"From";B$;P$(I);"&nbsp";:INPUT&nbspP$(I)
65&nbspCLS&nbsp:PRINT&nbsp"To";B$;P$(1,I);"&nbsp";:INPUT&nbspP$(1,I)
70&nbspCLS&nbsp:NEXT&nbspI:END&nbsp
80&nbsp"C"GOSUB&nbsp1:LPRINT&nbsp"&nbspObserved":LPRINT&nbsp"No.&nbspdh&nbspS":GOSUB&nbsp1
85&nbspFOR&nbspI=1TO&nbspT:GOSUB&nbsp2:LPRINT&nbspSTR$ I;TAB 3;USING&nbspU$;D(I);TAB 12;USING&nbsp"###.##";S:USING&nbsp:IF&nbspN$="T"LET&nbspS(I)=S*S
90&nbspIF&nbspY>BAND&nbspVAL&nbspD$(Z-B-1)<Z-YLET&nbspD$(Z-B-1)=STR$&nbsp(Z-Y)
95&nbspH(Z)=H(Y)+D(I):NEXT&nbspI:GOSUB&nbsp1:IF&nbspUFOR&nbspI=1TO&nbspU:D$(I)=STR$&nbsp(VAL&nbspD$(I-1)+VAL&nbspD$(I)+1):NEXT&nbspI
100&nbspDIM&nbspN(INT&nbsp(VAL&nbspD$(U)/10),9):FOR&nbspI=1TO&nbspT:GOSUB&nbsp2:H=Z-B-1:X=VAL&nbspD$(H):J=1:GOSUB&nbsp3
105&nbspD=H(Z)-H(Y)-D(I):L(H)=L(H)+D/S
110&nbspIF&nbspY>BLET&nbspX=X+Y-Z:J=-1:GOSUB&nbsp3:H=Y-B-1:X=VAL&nbspD$(H):J=1:GOSUB&nbsp3:L(H)=L(H)-D/S
115&nbspL=L+D*D/S:NEXT&nbspI:S=0:USING&nbspU$:FOR&nbspH=0TO&nbspU:X=VAL&nbspD$(H):GOSUB&nbsp3:P=M:Q=N:A=X-VAL&nbspD$(H-SGN&nbspX)+(H=0)
120&nbspW=X:FOR&nbspI=HTO&nbspU:Y=VAL&nbspD$(I)-I:Z=VAL&nbspD$(I-SGN&nbspI)-Y-(I=0):IF&nbspH<=ZTHEN&nbsp155
125&nbspX=Y+H:GOSUB&nbsp3:F=M:G=N:IF&nbspA=1THEN&nbsp145
130&nbspFOR&nbspJ=H-A+1TO&nbspH-1:IF&nbspJ<=ZTHEN&nbsp140
135&nbspX=W-H+J:GOSUB&nbsp3:C=M:R=N:X=Y+J:GOSUB&nbsp3:N(F,G)=N(F,G)-N(C,R)*N(M,N)
140&nbspNEXT&nbspJ
145&nbspIF&nbspI=HLET&nbspO=SQR&nbspN(P,Q)
150&nbspN(F,G)=N(F,G)/O
155&nbspNEXT&nbspI:IF&nbspA>1FOR&nbspJ=H-A+1TO&nbspH-1:X=W-H+J:GOSUB&nbsp3:L(H)=L(H)-N(M,N)*L(J):NEXT&nbspJ
160&nbspL(H)=L(H)/O:L=L-L(H)^2:NEXT&nbspH:FOR&nbspH=UTO&nbsp0STEP&nbsp-1:W=-L(H)
165&nbspIF&nbspH<UFOR&nbspI=H+1TO&nbspU:X=VAL&nbspD$(I)-I+H:IF&nbspX>VAL&nbspD$(I-1)GOSUB&nbsp3:W=W-N(M,N)*L(I):NEXT&nbspI
170&nbspX=VAL&nbspD$(H):GOSUB&nbsp3:W=W/N(M,N):L(H)=W:H(B+1+H)=H(B+1+H)+W:NEXT&nbspH:C=.00005:O=.005
175&nbspGOSUB&nbsp1:LPRINT&nbspTAB&nbsp10;"Adjusted":LPRINT&nbspTAB 13;"dh":GOSUB&nbsp1:P=1000:IF&nbspN$="T"LET&nbspC=.0005:P=100
180&nbspFOR&nbspI=1TO&nbspT:GOSUB&nbsp2:D=H(Z)-H(Y):LPRINT&nbspD+C*SGN&nbspD:D(I)=D:NEXT&nbspI:GOSUB&nbsp1:GOSUB&nbsp1:LPRINT&nbspTAB 13;"1/Ph":S=0
185&nbspLPRINT&nbspB$;TAB&nbsp7;"H&nbspMh(";M$;")":GOSUB&nbsp1:V=SQR&nbspL/SQR&nbsp(T-U-1)*P:FOR&nbspH=-BTO&nbspU
190&nbspD=H+B+1:USING&nbspU$:LPRINT&nbspSTR$ D;:COLOR&nbsp(H<0)*3:LPRINT&nbspTAB 2;H(D)+C*SGN&nbspH(D);:IF&nbspH<0LPRINT&nbsp:LPRINT&nbsp:GOTO&nbsp220
195&nbspQ=0:FOR&nbspI=HTO&nbspU:L(I)=0:NEXT&nbspI:L(H)=1:FOR&nbspI=HTO&nbspU:Z=VAL&nbspD$(I):IF&nbspI=HTHEN&nbsp210
200&nbspA=I-Z+VAL&nbspD$(I-1)+1:E=H:IF&nbspA>HLET&nbspE=A
205&nbspFOR&nbspJ=ETO&nbspI-1:X=Z-I+J:GOSUB&nbsp3:L(I)=L(I)-N(M,N)*L(J):NEXT&nbspJ
210&nbspX=Z:GOSUB&nbsp3:L(I)=L(I)/N(M,N):Q=Q+L(I)^2:NEXT&nbspI:COLOR&nbsp2:USING&nbsp"###.##":LPRINT&nbspTAB 12;Q+O
215&nbspLPRINT&nbspTAB 12;V*SQR&nbspQ+O
220&nbspCOLOR&nbsp0:NEXT&nbspH:GOSUB&nbsp1:LPRINT&nbsp:COLOR&nbsp2:LPRINT&nbsp"Mo=";V+O;M$:COLOR&nbsp0:USING&nbsp:END&nbsp
10001&nbspGRAPH&nbsp:LINE&nbsp-(216,0):TEXT&nbsp:LPRINT&nbsp:RETURN&nbsp
OR&nbsp0:NEXT&nbspH:GOSUB&nbsp1:LPRINT&nbsp:COLOR&nbsp2:LPRINT&nbsp"Mo=";V+O;M$:COL

发表评论

您的邮箱不会被公开。* 为必填项