Черный А.А. Математическое моделирование
Подождите немного. Документ загружается.
121
Таблица 51
План 3
3
(Х = 27)
Факторы
Номер опы-
та
F(J) H(J) L(J)
Показатель
Y(J)
1 A1 A2 A3 Y(1)
2 B1 A2 A3 Y(2)
3 A1 B2 A3 Y(3)
4 B1 B2 A3 Y(4)
5 A1 A2 B3 Y(5)
6 B1 A2 B3 Y(6)
7 A1 B2 B3 Y(7)
8 B1 B2 B3 Y(8)
9 A1 E2 E3 Y(9)
10 B1 E2 E3 Y(10)
11 E1 A2 E3 Y(11)
12 E1 B2 E3 Y(12)
13 E1 E2 A3 Y(13)
14 E1 E2 B3 Y(14)
15 A1 A2 E3 Y(15)
16 B1 A2 E3 Y(16)
17 A1 B2 E3 Y(17)
18 B1 B2 E3 Y(18)
19 A1 E2 A3 Y(19)
20 B1 E2 A3 Y(20)
21 A1 E2 B3 Y(21)
22 B1 E2 B3 Y(22)
23 E1 A2 A3 Y(23)
24 E1 B2 A3 Y(24)
25 E1 A2 B3 Y(25)
26 E1 B2 B3 Y(26)
27 E1 E2 E3 Y(27)
122
Таблица 52
План 4 ⋅ 5 (Х = 20)
Факторы
Номер опыта
F(J) H(J)
Показатель
Y(J)
1 A1 A2 Y(1)
2 B1 A2 Y(2)
3 A1 B2 Y(3)
4 B1 B2 Y(4)
5 A1 E2 Y(5)
6 B1 E2 Y(6)
7 A1 C2 Y(7)
8 B1 D2 Y(8)
9 A1 D2 Y(9)
10 B1 C2 Y(10)
11 С1 A2 Y(11)
12 С1 C2 Y(12)
13 С1 E2 Y(13)
14 С1 D2 Y(14)
15 C1 B2 Y(15)
16 D1 A2 Y(16)
17 D1 C2 Y(17)
18 D1 E2 Y(18)
19 D1 D2 Y(19)
20 D1 B2 Y(20)
123
Таблица 53
План 3 ⋅ 4 (Х = 12)
Факторы
Номер опыта
F(J) H(J)
Показатель
Y(J)
1 A1 A2 Y(1)
2 B1 A2 Y(2)
3 A1 B2 Y(3)
4 B1 B2 Y(4)
5 E1 A2 Y(5)
6 E1 B2 Y(6)
7 A1 C2 Y(7)
8 B1 D2 Y(8)
9 A1 D2 Y(9)
10 B1 C2 Y(10)
11 E1 C2 Y(11)
12 E1 D2 Y(12)
Обозначение: Х – количество опытов по плану.
124
ПРОГРАММА МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ
VL0 ДЛЯ СЛУЧАЕВ ПЛАНИРОВАНИЯ 2
1
(Х = 2), 2
2
(Х = 4),
2
3
(Х = 8), 2
4
( Х = 16), 2
5
(Х = 32)
125
5 PRINT "ПРОГРАММА VL0,РАЗРАБОТКА А.А.ЧЕРНОГО"
6 CLS
7 PRINT "РЕЗУЛЬТАТЫ ВЫПОЛНЕНИЯ ПРОГРАММЫ VL0"
8 PRINT "ЗАНОСЯТСЯ В ФАЙЛ,ИМЯ КОТОРОГО НАДО ВВЕСТИ,"
9 PRINT "НАПРИМЕР, ВВЕСТИ ИМЯ ФАЙЛА VL01"
10 INPUT "ВВОД ИМЕНИ ФАЙЛА ", FA$
11 OPEN "O", #1, FA$
17 PRINT #1, "РЕЗУЛЬТАТЫ ВЫПОЛНЕНИЯ ПРОГРАММЫ VL0 ЗАНО-
СЯТСЯ В ФАЙЛ "; FA$
40 PRINT "X=2,X=4,X=8,X=16,X=32"
41 PRINT #1, "РЕЗУЛЬТАТЫ ВЫПОЛНЕНИЯ ПРОГРАММЫ
VL0,РАЗРАБОТАННОЙ А.А. ЧЕРНЫМ"
42 DIM F(50), H(50), L(50), K(50), M(50), Y(32), I(50), P(50)
44 DIM U(50), Q(50), V(50), O(32), B(32), Z(50), G(20), T(32)
46 DIM K6(50), K7(50), K8(50), J7(50), J8(50), J9(50)
55 PRINT "КОЛИЧЕСТВО ОПЫТОВ ПО ПЛАНУ "
56 PRINT #1, "ВВОД X-КОЛИЧЕСТВО ОПЫТОВ
ПО ПЛАНУ "
60 INPUT X: PRINT #1, "X="; X
61 PRINT "ВВОД ВЕЛИЧИН ФАКТОРОВ И ПОКАЗАТЕЛЕЙ СТЕПЕНИ"
62 PRINT #1, "ВЕЛИЧИНЫ ФАКТОРОВ И ПОКАЗАТЕЛЕЙ СТЕПЕНИ"
70 PRINT "ВВОД A1,B1,J1": INPUT A1, B1, J1
75 PRINT #1, "A1="; A1; "B1="; B1; "J1="; J1
80 A = A1: B = B1: N = J1: GOSUB 2480
90 V1 = V0: PRINT #1, "V1="; V1: IF X = 2 GOTO 220
100 PRINT "ВВОД A2,B2,J2": INPUT A2, B2, J2
105 PRINT #1, "A2="; A2; "B2="; B2; "J2="; J2
110 A = A2: B = B2: N = J2: GOSUB 2480
120 V2 = V0: PRINT #1, "V2="; V2: IF X = 4 GOTO 220
130 PRINT "ВВОД A3,B3,J3": INPUT A3, B3, J3
135 PRINT #1, "A3="; A3; "B3="; B3; "J3="; J3
140 A = A3: B = B3: N = J3: GOSUB 2480
150 V3 = V0: PRINT #1, "V3="; V3: IF X = 8 GOTO 220
160 PRINT "ВВОД A4,B4,J4": INPUT A4, B4, J4
165 PRINT #1, "A4="; A4; "B4="; B4; "J4="; J4
170 A = A4: B = B4: N = J4: GOSUB 2480
180 V4 = V0: PRINT #1, "V4="; V4: IF X = 16 GOTO 220
190 PRINT "ВВОД A5,B5,J5": INPUT A5, B5, J5
195 PRINT #1, "A5="; A5; "B5="; B5; "J5="; J5
200 A = A5: B = B5: N = J5: GOSUB 2480
210 V5 = V0: PRINT #1, "V5="; V5: GOTO 220
215 REM ПЛАНЫ
220 F(1) = A1: F(2) = B1: IF X = 2 GOTO 580
230 H(1) = A2: H(2) = A2: F(3) = A1: H(3) = B2
240 F(4) = B1: H(4) = B2: IF X = 4 GOTO 580
250 L(1) = A3: L(2) = A3: L(3) = A3: L(4) = A3
260 F(5) = A1: H(5) = A2: L(5) = B3: F(6) = B1: H(6) = A2
270 L(6) = B3: F(7) = A1: H(7) = B2: L(7) = B3: F(8) = B1
280 H(8) = B2: L(8) = B3: IF X = 8 GOTO 580
290 K(1) = A4: K(2) = A4: K(3) = A4: K(4) = A4: K(5) = A4
126
300 K(6) = A4: K(7) = A4: K(8) = A4: F(9) = A1: H(9) = A2
310 L(9) = A3: K(9) = B4: F(10) = B1: H(10) = A2: L(10) = A3
320 K(10) = B4: F(11) = A1: H(11) = B2: L(11) = A3: K(11) = B4
330 F(12) = B1: H(12) = B2: L(12) = A3: K(12) = B4: F(13) = A1
340 H(13) = A2: L(13) = B3: K(13) = B4: F(14) = B1: H(14) = A2
350 L(14) = B3: K(14) = B4: F(15) = A1: H(15) = B2: L(15) = B3
360 K(15) = B4: F(16) = B1: H(16) = B2: L(16) = B3: K(16) = B4
370 IF X = 16 GOTO 580
380 M(1) = A5: M(2) = A5: M(3) = A5: M(4) = A5: M(5) = A5
390 M(6) = A5: M(7) = A5: M(8) = A5: M(9) = A5: M(10) = A5
400 M(11) = A5: M(12) = A5: M(13) = A5: M(14) = A5: M(15) = A5
410 M(16) = A5: F(17) = A1: H(17) = A2: L(17) = A3: K(17) = A4
420 M(17) = B5: F(18) = B1: H(18) = A2: L(18) = A3: K(18) = A4
430 M(18) = B5: F(19) = A1: H(19) = B2: L(19) = A3: K(19) = A4
440 M(19) = B5: F(20) = B1: H(20) = B2: L(20) = A3: K(20) = A4
450 M(20) = B5: F(21) = A1: H(21) = A2: L(21) = B3: K(21) = A4
460 M(21) = B5: F(22) = B1: H(22) = A2: L(22) = B3: K(22) = A4
470 M(22) = B5: F(23) = A1: H(23) = B2: L(23) = B3: K(23) = A4
480 M(23) = B5: F(24) = B1: H(24) = B2: L(24) = B3: K(24) = A4
490 M(24) = B5: F(25) = A1: H(25) = A2: L(25) = A3: K(25) = B4
500 M(25) = B5: F(26) = B1: H(26) = A2: L(26) = A3: K(26) = B4
510 M(26) = B5: F(27) = A1: H(27) = B2: L(27) = A3: K(27) = B4
520 M(27) = B5: F(28) = B1: H(28) = B2: L(28) = A3: K(28) = B4
530 M(28) = B5: F(29) = A1: H(29) = A2: L(29) = B3: K(29) = B4
540 M(29) = B5: F(30) = B1: H(30) = A2: L(30) = B3: K(30) = B4
550 M(30) = B5: F(31) = A1: H(31) = B2: L(31) = B3: K(31) = B4
560 M(31) = B5: F(32) = B1: H(32) = B2: L(32) = B3: K(32) = B4
570 M(32) = B5
580 PRINT "IF I0=6 GOTO 40-НАЧАЛО"
585 PRINT "IF I0=7 GOTO 610-ПРОДОЛЖЕНИЕ"
590 INPUT I0: IF I0 = 6 GOTO 40
600 IF I0 = 7 GOTO 610
610 PRINT "ВВОД ВЕЛИЧИН ПОКАЗАТЕЛЕЙ ПО ПЛАНУ Y(J)"
611 PRINT #1, "ВВОД ВЕЛИЧИН ПОКАЗАТЕЛЕЙ В СООТВЕТСТВИИ С
ПЛАНОМ Y(J)"
615 FOR J = 1 TO X: PRINT "ВВОД Y("; J; ")": INPUT Y(J)
620 PRINT #1, "Y("; J; ")="; Y(J): NEXT J
630 PRINT "IF I0=1 GOTO 610-ПОВТОРЕНИЕ ВВОДА ПОКАЗАТЕЛЕЙ"
635 PRINT "IF I0=2 GOTO 660-ПРОДОЛЖЕНИЕ"
640 PRINT "ВВОД I0": INPUT I0: IF I0 = 1 GOTO 610
650 IF I0 = 2 GOTO 660
660 IF X = 2 GOTO 710
670 IF X = 4 GOTO 720
680 IF X = 8 GOTO 730
690 IF X = 16 GOTO 740
700 IF X = 32 GOTO 760
710 GOSUB 2490: GOTO 780
720 GOSUB 2490: GOSUB 2500: GOTO 780
730 GOSUB 2490: GOSUB 2500: GOSUB 2510: GOTO 780
740 GOSUB 2490: GOSUB 2500: GOSUB 2510: GOSUB 2520
127
750 GOTO 780
760 GOSUB 2490: GOSUB 2500: GOSUB 2510: GOSUB 2520
770 GOSUB 2530: GOTO 780
780 S = 0: O(1) = 0: FOR J = 1 TO X: S = S + Y(J): O(1) = O(1) + 1: NEXT J
790 B(1) = S / O(1): S = 0: O(2) = 0: FOR J = 1 TO X: S = S + I(J) * Y(J)
800 O(2) = O(2) + I(J) ^ 2: NEXT J: B(2) = S / O(2): IF X = 2 GOTO 1440
810 S = 0: O(3) = 0: FOR J = 1 TO X: S = S + P(J) * Y(J): O(3) = O(3) + P(J) ^ 2
820 NEXT J: B(3) = S / O(3): S = 0: O(4) = 0: FOR J = 1 TO X
830 S = S + I(J) * P(J) * Y(J): O(4) = O(4) + (I(J) * P(J)) ^ 2: NEXT J
840 B(4) = S / O(4): IF X = 4 GOTO 1440
850 S = 0: O(5) = 0
860 FOR J = 1 TO X: S = S + U(J) * Y(J): O(5) = O(5) + U(J) ^ 2
870 NEXT J: B(5) = S / O(5): S = 0: O(6) = 0: FOR J = 1 TO X
880 S = S + I(J) * U(J) * Y(J): O(6) = O(6) + (I(J) * U(J)) ^ 2: NEXT J
890 B(6) = S / O(6): S = 0: O(7) = 0: FOR J = 1 TO X: S = S + P(J) * U(J) * Y(J)
900 O(7) = O(7) + (P(J) * U(J)) ^ 2: NEXT J: B(7) = S / O(7): S = 0: O(8) = 0
910 FOR J = 1 TO X: S = S + I(J) * P(J) * U(J) * Y(J)
920 O(8) = O(8) + (I(J) * P(J) * U(J)) ^ 2: NEXT J: B(8) = S / O(8)
930 IF X = 8 GOTO 1440
940 S = 0: O(9) = 0: FOR J = 1 TO X
950 S = S + Q(J) * Y(J): O(9) = O(9) + Q(J) ^ 2: NEXT J: B(9) = S / O(9)
960 S = 0: O(10) = 0: FOR J = 1 TO X: S = S + I(J) * Q(J) * Y(J)
970 O(10) = O(10) + (I(J) * Q(J)) ^ 2: NEXT J: B(10) = S / O(10): S = 0
980 O(11) = 0: FOR J = 1 TO X: S = S + P(J) * Q(J) * Y(J)
990 O(11) = O(11) + (P(J) * Q(J)) ^ 2: NEXT J: B(11) = S / O(11): S = 0
1000 O(12) = 0: FOR J = 1 TO X: S = S + I(J) * P(J) * Q(J) * Y(J)
1010 O(12) = O(12) + (I(J) * P(J) * Q(J)) ^ 2: NEXT J: B(12) = S / O(12)
1020 S = 0: O(13) = 0: FOR J = 1 TO X: S = S + U(J) * Q(J) * Y(J)
1030 O(13) = O(13) + (U(J) * Q(J)) ^ 2: NEXT J: B(13) = S / O(13): S = 0
1040 O(14) = 0: FOR J = 1 TO X: S = S + I(J) * U(J) * Q(J) * Y(J)
1050 O(14) = O(14) + (I(J) * U(J) * Q(J)) ^ 2: NEXT J: B(14) = S / O(14): S = 0
1060 O(15) = 0: FOR J = 1 TO X: S = S + P(J) * U(J) * Q(J) * Y(J)
1070 O(15) = O(15) + (P(J) * U(J) * Q(J)) ^ 2: NEXT J: B(15) = S / O(15): S = 0
1080 O(16) = 0: FOR J = 1 TO X: S = S + I(J) * P(J) * U(J) * Q(J) * Y(J)
1090 O(16) = O(16) + (I(J) * P(J) * U(J) * Q(J)) ^ 2: NEXT J: B(16) = S / O(16)
1100 IF X = 16 GOTO 1440
1110 S = 0: O(17) = 0: FOR J = 1 TO X
1120 S = S + V(J) * Y(J): O(17) = O(17) + V(J) ^ 2: NEXT J: B(17) = S / O(17)
1130 S = 0: O(18) = 0: FOR J = 1 TO X: S = S + I(J) * V(J) * Y(J)
1140 O(18) = O(18) + (I(J) * V(J)) ^ 2: NEXT J: B(18) = S / O(18): S = 0
1150 O(19) = 0: FOR J = 1 TO X: S = S + P(J) * V(J) * Y(J)
1160 O(19) = O(19) + (P(J) * V(J)) ^ 2: NEXT J: B(19) = S / O(19): S = 0
1170 O(20) = 0: FOR J = 1 TO X: S = S + I(J) * P(J) * V(J) * Y(J)
1180 O(20) = O(20) + (I(J) * P(J) * V(J)) ^ 2: NEXT J: B(20) = S / O(20)
1190 S = 0: O(21) = 0: FOR J = 1 TO X: S = S + U(J) * V(J) * Y(J)
1200 O(21) = O(21) + (U(J) * V(J)) ^ 2: NEXT J: B(21) = S / O(21)
1210 S = 0: O(22) = 0: FOR J = 1 TO X: S = S + I(J) * U(J) * V(J) * Y(J)
1220 O(22) = O(22) + (I(J) * U(J) * V(J)) ^ 2: NEXT J: B(22) = S / O(22)
1230 S = 0: O(23) = 0: FOR J = 1 TO X: S = S + P(J) * U(J) * V(J) * Y(J)
1240 O(23) = O(23) + (P(J) * U(J) * V(J)) ^ 2: NEXT J: B(23) = S / O(23)
128
1250 S = 0: O(24) = 0: FOR J = 1 TO X: S = S + I(J) * P(J) * U(J) * V(J) * Y(J)
1260 O(24) = O(24) + (I(J) * P(J) * U(J) * V(J)) ^ 2: NEXT J: B(24) = S / O(24)
1270 S = 0: O(25) = 0: FOR J = 1 TO X: S = S + Q(J) * V(J) * Y(J)
1280 O(25) = O(25) + (Q(J) * V(J)) ^ 2: NEXT J: B(25) = S / O(25)
1290 S = 0: O(26) = 0: FOR J = 1 TO X: S = S + I(J) * Q(J) * V(J) * Y(J)
1300 O(26) = O(26) + (I(J) * Q(J) * V(J)) ^ 2: NEXT J: B(26) = S / O(26)
1310 S = 0: O(27) = 0: FOR J = 1 TO X: S = S + P(J) * Q(J) * V(J) * Y(J)
1320 O(27) = O(27) + (P(J) * Q(J) * V(J)) ^ 2: NEXT J: B(27) = S / O(27)
1330 S = 0: O(28) = 0: FOR J = 1 TO X: S = S + I(J) * P(J) * Q(J) * V(J) * Y(J)
1340 O(28) = O(28) + (I(J) * P(J) * Q(J) * V(J)) ^ 2: NEXT J: B(28) = S / O(28)
1350 S = 0: O(29) = 0: FOR J = 1 TO X: S = S + U(J) * Q(J) * V(J) * Y(J)
1360 O(29) = O(29) + (U(J) * Q(J) * V(J)) ^ 2: NEXT J: B(29) = S / O(29)
1370 S = 0: O(30) = 0: FOR J = 1 TO X: S = S + I(J) * U(J) * Q(J) * V(J) * Y(J)
1380 O(30) = O(30) + (I(J) * U(J) * Q(J) * V(J)) ^ 2: NEXT J: B(30) = S / O(30)
1390 S = 0: O(31) = 0: FOR J = 1 TO X: S = S + P(J) * U(J) * Q(J) * V(J) * Y(J)
1400 O(31) = O(31) + (P(J) * U(J) * Q(J) * V(J)) ^ 2: NEXT J: B(31) = S / O(31)
1410 S = 0: O(32) = 0: FOR J = 1 TO X: S = S + I(J) * P(J) * U(J) * Q(J) * V(J) *
Y(J)
1420 O(32) = O(32) + (I(J) * P(J) * U(J) * Q(J) * V(J)) ^ 2: NEXT J
1430 B(32) = S / O(32): GOTO 1440
1440 PRINT #1, "B(J) ДО АНАЛИЗА": FOR J = 1 TO X
1445 PRINT #1, "B("; J; ")="; B(J)
1450 NEXT J: GOTO 1454
1454 PRINT "ВВОД N0-КОЛИЧЕСТВО ОПЫТОВ НА СРЕДНЕМ УРОВНЕ
ФАКТОРОВ"
1456 PRINT #1, "N0-КОЛИЧЕСТВО ОПЫТОВ НА СРЕДНЕМ УРОВНЕ ФАК-
ТОРОВ"
1455 INPUT N0
1460 PRINT #1, "N0="; N0
1470 IF X = 2 GOTO 1520
1480 IF X = 4 GOTO 1530
1490 IF X = 8 GOTO 1540
1500 IF X = 16 GOTO 1550
1510 IF X = 32 GOTO 1560
1520 GOSUB 2540: GOTO 1570
1530 GOSUB 2550: GOTO 1570
1540 GOSUB 2570: GOTO 1570
1550 GOSUB 2600: GOTO 1570
1560 GOSUB 2670: GOTO 1570
1570 PRINT #1, "РАСЧЕТНЫЕ.ВЕЛИЧИНЫ ПОКАЗАТЕЛЯ Z(J) ДО АНАЛИ-
ЗА B(J)"
1580 FOR J = 1 TO X: PRINT #1, "Z("; J; ")="; Z(J): NEXT J
1590 PRINT #1, "ПРОВЕРКА ПО РАЗНОСТИ Y(J)-Z(J)"
1600 PRINT #1, "В
ПРОЦЕНТАХ (Y(J)-Z(J) * (100 / Y(J))"
1601 FOR J = 1 TO X: PRINT #1, "Y("; J; ")-Z("; J; ")="; Y(J) - Z(J)
1602 PRINT #1, "(Y("; J; ")-Z("; J; ")) * (100 / Y("; J; ")) = "; (Y(J) - Z(J)) * (100 /
Y(J))
1603 NEXT J
1604 PRINT "ВВОД F8=N0-1": INPUT F8: PRINT #1, "F8=N0-1="; N0; "-1="; F8
1620 PRINT "IF I0=3 GOTO 1710-ВВОД РЕЗУЛЬТАТОВ ОПЫТОВ"
129
1622 PRINT " НА СРЕДНЕМ УРОВНЕ ФАКТОРОВ"
1625 PRINT "IF I0=4 GOTO 1760-ВВОД ДИСПЕРСИИ ОПЫТОВ"
1630 PRINT "IF I0=5 GOTO 2150-ПРОВЕРКА ТОЧНОСТИ И"
1633 PRINT " РАСЧЕТЫ ПО МОДЕЛИ"
1635 PRINT "IF I0=6 GOTO 40-НАЧАЛО"
1640 PRINT "IF I0=20 GOTO 6830-КОНЕЦ"
1641 PRINT "IF I0=25 GOTO 2820-МАТЕМАТИЧЕСКАЯ МОДЕЛЬ"
1642 PRINT "IF I0=27 GOTO 7000-ВЫЧИСЛЕНИЯ ПОКАЗАТЕЛЕЙ Z(K5)"
1646 PRINT " С ИСПОЛЬЗОВАНИЕМ ЦИКЛОВ И"
1647 PRINT " ПОСТРОЕНИЕ ГРАФИКОВ"
1650 PRINT "ВВОД I0": INPUT I0: IF I0 = 3 GOTO 1710
1660 IF I0 = 4 GOTO 1760
1670 IF I0 = 5 GOTO 2150
1680 IF I0 = 6 GOTO 40
1690 IF I0 = 20 GOTO 6830
1700 IF I0 = 25 GOTO 2820
1705 IF I0 = 27 GOTO 7000
1710 PRINT "ВВОД РЕЗУЛЬТАТОВ ОПЫТОВ "
1711 PRINT "НА СРЕДНЕМ УРОВНЕ "
1712 PRINT #1, "РЕЗУЛЬТАТЫ ОПЫТОВ ";
1713 PRINT #1, "НА СРЕДНЕМ УРОВНЕ ФАКТОРОВ "
1715 FOR J = 1 TO N0: PRINT "ВВОД G("; J; ")": INPUT G(J)
1720 PRINT #1, "G("; J; ")="; G(J): NEXT J: S = 0: FOR J = 1 TO N0: S = S + G(J)
1730 NEXT J: S0 = S / N0: PRINT "S0="; S0: S = 0: FOR J = 1 TO N0
1740 S = S + (G(J) - S0) ^ 2: NEXT J: U9 = S / F8
1745 PRINT #1, "ДИСПЕРСИЯ ОПЫТОВ U9="; U9
1750 GOTO 1770
1760 PRINT "ВВОД U9": INPUT U9: PRINT #1, "ДИСПЕРСИЯ ОПЫТОВ U9=";
U9
1770 PRINT #1, "РАСЧЕТНЫЕ ВЕЛИЧИНЫ T(J)": FOR J = 1 TO X
1780 T(J) = ABS(B(J) / SQR(U9 / O(J))): PRINT #1, "T("; J; ")="; T(J): NEXT J
1781 PRINT " ДЛЯ УРОВНЯ ЗНАЧИМОСТИ 5% "
1782 PRINT " ПРИ F8 2 ; 3 ; 4 ; 5 ; 6"
1783 PRINT "СООТВЕТСТВЕННО T0 4.303 ;3.182 ;2.776 ;2.571 ;2.447"
1784 PRINT "F8="; F8
1785 PRINT "ВВОД T0-ТАБЛИЧНЫЙ Т-КРИТЕРИЙ"
1790 PRINT "ВВОД T0": INPUT T0: PRINT #1, "ТАБЛИЧНЫЙ Т-КРИТЕРИЙ
T0="; T0
1800 PRINT #1, "B(J) ПОСЛЕ АНАЛИЗА": FOR J = 1 TO X
1810 IF T(J) < T0 GOTO 1830
1820 IF T(J) >= T0 GOTO 1840
1830 B(J) = 0
1840 PRINT #1, "B("; J; ")="; B(J): NEXT J
1850 K9 = 0: FOR J = 1 TO X: IF B(J) = 0 GOTO 1852
1851 K9 = K9 + 1: NEXT J
1852 PRINT #1, "КОЛИЧЕСТВО СТАТИСТИЧЕСКИ ЗНАЧИМЫХ"
1853 PRINT #1, " КОЭФФИЦИЕНТОВ РЕГРЕССИИ K9="; K9
1862 PRINT "F9=X-1": F9 = X - 1: PRINT #1, "F9=X-1="; F9: CLS
1863 PRINT "! ! ЗНАЧЕНИЯ F-КРИТЕРИЯ F7 ДЛЯ 5%-ГО УРОВНЯ ЗНА-
ЧИМОСТИ"
130
1864 PRINT "! !----------------------------------------------"
1865 PRINT "! F8! F9 "
1866 PRINT "! !----------------------------------------------"
1867 PRINT "! ! 1 ! 3 ! 7 ! 15...16 ! 30...31 !"
1868 PRINT "---------------------------------------------------"
1869 PRINT "! 2 ! 18.51 ! 19.16 ! 19.36 ! 19.43 ! 19.46 !"
1870 PRINT "! 3 ! 10.13 ! 9.28 ! 8.88 ! 8.69 ! 8.62 !"
1871 PRINT "! 4 ! 7.71 ! 6.59 ! 6.09 ! 5.84 ! 5.74 !"
1872 PRINT "! 5 ! 6.61 ! 5.41 ! 4.88 ! 4.6 ! 4.5 !"
1873 PRINT "! 6 ! 5.99 ! 4.76 ! 4.21 ! 3.92 ! 3.81 !"
1874 PRINT
"!=================================================!"
1887 PRINT "F8="; F8; "F9="; F9
1890 PRINT "ВВОД ТАБЛИЧНОГО F7-ТАБЛИЧНЫЙ F-КРИТЕРИЙ"
1891 INPUT F7: PRINT #1, "ТАБЛИЧНЫЙ F-КРИТЕРИЙ F7="; F7
1900 IF X = 2 GOTO 1950
1910 IF X = 4 GOTO 1960
1920 IF X = 8 GOTO 1970
1930 IF X = 16 GOTO 1980
1940 IF X = 32 GOTO 1990
1950 GOSUB 2540: GOTO 2000
1960 GOSUB 2550: GOTO 2000
1970 GOSUB 2570: GOTO 2000
1980 GOSUB 2600: GOTO 2000
1990 GOSUB 2670: GOTO 2000
2000 PRINT #1, "РАСЧЕТНЫЕ ВЕЛИЧИНЫ ПОКАЗАТЕЛЯ ";
2005 PRINT #1, "Z(J) ПОСЛЕ АНАЛИЗА B(J)"
2010 FOR J = 1 TO X: PRINT #1, "Z("; J; ")="; Z(J): NEXT J
2011 PRINT #1, "ПРОВЕРКА ПО РАЗНОСТИ Y(J)-Z(J)"
2012 PRINT #1, "В ПРОЦЕНТАХ (Y(J)-Z(J) * (100 / Y(J))"
2013 FOR J = 1 TO X: PRINT #1, "Y("; J; ")-Z("; J; ")="; Y(J) - Z(J)
2014 PRINT #1, "(Y("; J; ")-Z("; J; ")) * (100 / Y("; J; ")) = "; (Y(J) - Z(J)) * (100 /
Y(J))
2015 NEXT J
2020 S = 0: FOR J = 1 TO X: S = S + (Z(J) - Y(J)) ^ 2: NEXT J
2025 F6 = S / (F9 * U9)
2030 PRINT #1, "РАСЧЕТНАЯ ВЕЛИЧИНА F-КРИТЕРИЯ F6="; F6
2040 IF F6 <= F7 GOTO 2060
2050 IF F6 > F7 GOTO 2070
2060 PRINT #1, "АДЕКВАТНО,ТАК КАК F6<=F7": GOTO 2080
2070 PRINT #1, "НЕАДЕКВАТНО,ТАК КАК F6>F7"; ""
2080 PRINT "I0=7 GOTO 2150-ПРОВЕРКА ТОЧНОСТИ И"
2081 PRINT " РАСЧЕТЫ ПО МОДЕЛИ"
2085 PRINT "I0=8 GOTO 40-НАЧАЛО"
2090 PRINT "I0=17 GOTO 2820-МАТЕМАТИЧЕСКАЯ
МОДЕЛЬ"
2091 PRINT "I0=9 GOTO 6830-КОНЕЦ"
2095 PRINT "IF I0=22 GOTO 7000-ВЫЧИСЛЕНИЯ ПОКАЗАТЕЛЕЙ Z(K5)"
2096 PRINT " C ИСПОЛЬЗОВАНИЕМ ГРАФИКОВ И"
2097 PRINT " ПОСТРОЕНИЕ ГРАФИКОВ"
2100 PRINT "ВВОД I0": INPUT I0