Основы научных исследований. Ч.2. Килов А.С. - 12 стр.

UptoLike

Составители: 

10 REM ПРОГРАММА ДЛЯ ОБРАБОТКИ РЕЗУЛЬТАТОВ
20 REM ПОЛНОГО ТРЕХФАКТОРНОГО ЭКСПЕРИМЕНТА
30 DIM X(8,3), Y(8), B(6), S(6), T1 (6)
40 DATA -1, 1, -1, 1, -1, 1, -1, 1
50 DATA -1, -1, 1, 1, -1, -1, 1, 1
60 DATA -1, -1, -1, -1, 1, 1, 1, 1
70 J=1 TO 3: FOR I=1 TO 8
80 READ X (I,J) : NEXT I : NEXT J
90 FOR I=1 TO 8
100 INPUT “экспериментальные значения Y”; Y(I)
110 NEXT I
120 INPUT “оценка дисперсии воспроизводимости”; S1
130 INPUT “критерий Стьюдента”; T
140 FOR I=0 TO 6: LET S(I)=0 : NEXT I
150 FOR I=1 TO 8: LET S(0)=S(0)+Y(I)
160 LET S(1)=S(1)+X(I,1)*Y(I)
170 LET S(2)=S(2)+X(I,2)*Y(I)
180 LET S(3)=S(3)+X(I,3)*Y(I)
190 LET S(4)=S(4)+X(I,1)*X(I,2)*Y(I)
200 LET S(5)=S(5)+X(I,1)* X(I,3)*Y(I)
210 LET S(6)=S(6)+X(I,2)* X(I,3)*Y(I)
220 NEXT I
230 LET S2 = SQR(S1/8)
240 PRINT “коэффициенты регрессии”;
250 FOR I=0 TO 6
260 LET B(I)=S(I)/8:LET T1(I)=ABS (B(I)/S2)
270 PRINT B(I):NEXT I:LET Z=0
280 PRINT “значимые коэффициенты
290 FOR I=0 TO 6
300 IF T1(I)=T THEN GOTO 320
310 LET B(I)=0 : GOTO 330
320 LET Z=Z+1
330 NEXT I
340 FOR I=0 TO 6: PRINT B(I)
350 FOR J=1 TO 8: LET Y2(J)=B(0)
360 FOR I=1 TO 3
370 LET Y2(J)=Y2(J)+B(I)*X(J,I)
380 NEXT I
390 LET Y2(J)=Y2(J)+X(J,1)*X(J,2)*B(4)
400 LET Y2(J)=Y2(J)+X(J,1)*X(J,3)*B(5)
410 LET Y2(J)=Y2(J)+X(J,2)*X(J,3)*B(6)
420 NEXT J : LET S3=0
430 PRINT “расчетные значения, Y”;
440 FOR I=1 TO 8 : PRINT Y2(I)
10 REM ПРОГРАММА ДЛЯ ОБРАБОТКИ РЕЗУЛЬТАТОВ
20 REM ПОЛНОГО ТРЕХФАКТОРНОГО ЭКСПЕРИМЕНТА
30 DIM X(8,3), Y(8), B(6), S(6), T1 (6)
40 DATA -1, 1, -1, 1, -1, 1, -1, 1
50 DATA -1, -1, 1, 1, -1, -1, 1, 1
60 DATA -1, -1, -1, -1, 1, 1, 1, 1
70 J=1 TO 3: FOR I=1 TO 8
80 READ X (I,J) : NEXT I : NEXT J
90 FOR I=1 TO 8
100 INPUT “экспериментальные значения Y”; Y(I)
110 NEXT I
120 INPUT “оценка дисперсии воспроизводимости”; S1
130 INPUT “критерий Стьюдента”; T
140 FOR I=0 TO 6: LET S(I)=0 : NEXT I
150 FOR I=1 TO 8: LET S(0)=S(0)+Y(I)
160 LET S(1)=S(1)+X(I,1)*Y(I)
170 LET S(2)=S(2)+X(I,2)*Y(I)
180 LET S(3)=S(3)+X(I,3)*Y(I)
190 LET S(4)=S(4)+X(I,1)*X(I,2)*Y(I)
200 LET S(5)=S(5)+X(I,1)* X(I,3)*Y(I)
210 LET S(6)=S(6)+X(I,2)* X(I,3)*Y(I)
220 NEXT I
230 LET S2 = SQR(S1/8)
240 PRINT “коэффициенты регрессии”;
250 FOR I=0 TO 6
260 LET B(I)=S(I)/8:LET T1(I)=ABS (B(I)/S2)
270 PRINT B(I):NEXT I:LET Z=0
280 PRINT “значимые коэффициенты”
290 FOR I=0 TO 6
300 IF T1(I)=T THEN GOTO 320
310 LET B(I)=0 : GOTO 330
320 LET Z=Z+1
330 NEXT I
340 FOR I=0 TO 6: PRINT B(I)
350 FOR J=1 TO 8: LET Y2(J)=B(0)
360 FOR I=1 TO 3
370 LET Y2(J)=Y2(J)+B(I)*X(J,I)
380 NEXT I
390 LET Y2(J)=Y2(J)+X(J,1)*X(J,2)*B(4)
400 LET Y2(J)=Y2(J)+X(J,1)*X(J,3)*B(5)
410 LET Y2(J)=Y2(J)+X(J,2)*X(J,3)*B(6)
420 NEXT J : LET S3=0
430 PRINT “расчетные значения, Y”;
440 FOR I=1 TO 8 : PRINT Y2(I)