Программа статистической обработки результатов испытаний строительных материалов

1 REM ПРОГРАММА `МИКС` - МЕТОДЫ ИССЛЕДОВАНИЯ

10 REM СТРОИТЕЛЬНЫХ МАТЕРИАЛОВ

20 PRINT` ПРОГРАММА СТАТИСТИЧЕСКОЙ ОБРАБОТКИ

РЕЗУЛЬТАТОВ`

30 PRINT` ИСПЫТАНИЙ СТРОИТЕЛЬНЫХ МАТЕРИАЛОВ, ВКЛЮЧАЮЩАЯ`

40 PRINT` 1) УПОРЯДОЧЕНИЕ МАССИВА ДАННЫХ ПО ВОЗРАСТАНИЮ;`

50 PRINT` 2) ПРОВЕРКУ ПО КРИТЕРИЮ ДИКСОНА ПРИНАДЛЕЖНОСТИ ЭКСПЕРИМЕНТАЛЬНЫХ

60 PRINT` ДАННЫХ ОДНОЙ И ТОЙ ЖЕ ГЕНЕРАЛЬНОЙ СОВОКУПНОСТИ;`

70 PRINT` 3) ОЦЕНКУ ОСНОВНЫХ ПАРАМЕТРОВ РАСПРЕДЕЛЕНИЯ;`

80 PRINT` 4) ВЫЧИСЛЕНИЕ 95% ДОВЕРИТЕЛЬНЫХ ИНТЕРВАЛОВ ДЛЯ МАТЕМАТИЧЕСКОГО`

90 PRINT` ОЖИДАНИЯ.

100 REM S – СУММА РЕЗУЛЬТАТОВ НАБЛЮДЕНИЙ;

S₁ – ОЦЕНКА МАТЕМАТИЧЕСКОГО

110 REM ОЖИДАНИЯ; D₁ – ОЦЕНКА ДИСПЕРСИИ; S₂ – ОЦЕНКА СРЕДНЕГО

120 REM КВАДРАТИЧЕСКОГО ОТКЛОНЕНИЯ; S₃ – СУММА КВАДРАТОВ РЕЗУЛЬТАТОВ

130 REM НАБЛЮДЕНИЙ; V – ОЦЕНКА КОЭФФИЦИЕНТА ИЗМЕНЧИВОСТИ;

I₁, I₂ -

140 REM СООТВЕТСТВЕННО, НИЖНЯЯ И ВЕРХНЯЯ ГРАНИЦЫ ДОВЕРИТЕЛЬНОГО ИНТЕРВАЛА

150 INPUT ВВЕДИТЕ ЧИСЛО ТОЧЕК N=`N

160 DIM X(N), T(30), D(30)

170 POR I=1 TO N: PRINT!3.3! `X(`1;:INPUT`)=`X(1)

180 NEXT 1

190 GOSUB 1000

195 N1=N

200 GOSUB 2000

210 S=0; S₃=0

220 POR I=1 TO N1: PRINT !3.3!`X(`1;`)=`X(1)

230 S=S+X(1): S₃=S₃+X(1)X(1)

240 NEXT 1

250 S₁=S/N₁

260 D₁=(S₃-N₁S₁-2)/W₁-1)

270 S₂-SQR(D₁)

280 V=(S₂/S₁)100

290 PRINT `S<sub>1</sub>=`!3.2!S₁,`D<sub>1</sub>=`D₁, `S<sub>2</sub>=`S₂, V=`V

300 K=N1-1

310 GOSUB 2500

320 I₁=S1-T(K)S₂/SQR(N1)

330 I₂=S1+T(K)S₂/SQR(N1)

340 PRINT `(I<sub>1</sub>=`I₁;`,``I<sub>2</sub>=`I₂`)`

350 END

1000 REM УПОРЯДОЧЕНИЕ МАССИВА ПО ВОЗРАСТАНИЮ

Продолжение прил. 2

1010 POR I=1 TO N-1

1020 M=1; P=X(1)

1030 POR I=I+1 TO N

1040 IF X(I)P GOTO IOGO

1050 P=X(I) : N=1

1060 NEXT I

1070 X(M)=X(1); X(1)=P

1080 NEXT I

1090 POR I=1 TO N; PRINT !3.3! `X(`I;`)=`X(I)

1100 NEXT I

1110 RETURN

2000 REM РАСЧЕТ КРИТЕРИЯ ДИКСОНА НА ПРИНАДЛЕЖНОСТЬ ЭКСПЕРИМЕНТАЛЬНЫХ

2005 REM ДАННЫХ ГЕНЕРАЛЬНОЙ СОВОКУПНОСТИ

2010 REM ОБОЗНАЧЕНИЯ: R₁ – КРИТЕРИЙ ДЛЯ НАИМЕНЬШЕГО ЭКСПЕРИМЕНТАЛЬНОГО

2015 REM ЗНАЧЕНИЯ; R₂ – КРИТЕРИЙ ДЛЯ НАИБОЛЬШЕГО ЭКСПЕРИМЕНТАЛЬНОГО ЗНАЧЕНИЯ

2020 DIM R₁(N), R₂(N)

2030 IFN3 GOTO 2235

2040 IFN1 GOTO 2060

2050 R₁(N)=(X(2))-(X(1))/X(N)-X(I): R₂(N)=(X(N)-X(N-I))/(X(N)-X(I))

2055 GOTO 2110

2060 IFN10 GOTO 2060

2070 R₁(N)=(X(2)-X(1))/(X(N-1)-X(1)): R₂(N)=(X(N)-X(N-1))/(X(N)-X(2))

2075 GOTO 2110

2080 IFN13 GOTO 2100

2090 R₁(N)=(X(3)-X(1))/(X(N-1)-X(1)): R₂(N)=(X(N)-X(N-2))/(X(N)-

2X(2))

2095 GOTO 2110

2100 IFN30 GOTO 2235

2105 R₁(N)=(X(3)-X(1))/(X(N-2)-X(1)): R₂(N)=(X(N)-X(N-2))/(X(N)-Х(2))

2110 REM КОЭФФИЦИЕНТЫ ДИКСОНА, УЧИТЫВАЮЩИЕ ЭКСПЕРИМЕНТАЛЬНЫЕ

2115 REM ЗНАЧЕНИЯ ПРИ УРОВНЕ ЗНАЧИМОСТИ, РАВНОМ 0,05

2120 ДАТА,941,.765,.642,.560,.507,.554,.512,.477,.576,.546,.521

2125 ДАТА,546,.525,.507,.490,.475,.462,.450,.440,.430,.421,.413,.406,.399

2130 ДАТА,393,.387,.381,.376

2135 POR I=3 TO 30

2140 READ D(I); NTXT I

2145 IF R₁(N)D(N) THEN 2175

2150 PRINT !2.3! R₁(N)=`R<sub>1</sub>(N), `D(N), TO`

2160 PRINT `НАИМЕНЬШЕЕ ЭКСПЕРИМЕНТАЛЬНОЕ ЗНАЧЕНИЕ НЕ ПРИНАДЛЕЖИТ ГЕНЕРАЛЬНОЙ СОВОКУПНОСТИ

2165 POR I=2 TON:A=X(1):X(1-1)=A: NEXT I

2170 N₁=N1-1; GOTO 2190

2175 PRINT !2.3! `R<sub>1</sub>(N)=`R₁(N), `D(N)=`D(N)

2180 PRINT` ТАК КАК R₁(N)D(N), ТО НАИМЕНЬШЕЕ ЭКСПЕРИМЕНТАЛЬНОЕ ЗНАЧЕНИЕ

2185 PRINT` ПРИНАДЛЕЖИТ ГЕНЕРАЛЬНОЙ СОВОКУПНОСТИ

Окончание прил. 2

2190 IF R2(N)D(N) THEN 2220

2195 PRINT !2,3! R₂(N)=`R<sub>2</sub>(N), `D(N)=`D(N)

2200 PRINT` ТАК КАК R<sub>2</sub>(N)D(N), ТО НАИБОЛЬШЕЕ ЭКСПЕРИМЕНТАЛЬНОЕ ЗНАЧЕНИЕ`

2205 PRINT` НЕ ПРИНАДЛЕЖИТ ГЕНЕРАЛЬНОЙ СОВОКУПНОСТИ`

2210 N₁-N1-1

2215 GOTO 2235

2220 PRINT !2,3! `R<sub>2</sub>(N)=`R₂(N), `D(N)=`D(N)

2225 PRINT ТАК КАК R₂(N)D(N), ТО НАИБОЛЬШЕЕ ЭКСПЕРИМЕНТАЛЬНОЕ ЗНАЧЕНИЕ`

2230 PRINT ПРИНАДЛЕЖИТ ГЕНЕРАЛЬНОЙ СОВОКУПНОСТИ`

2235 RETURN

2500 REM КРИТИЧЕСКИЕ ТОЧКИ РАСПРЕДЕЛЕНИЯ СТЬЮДЕНТА ПРИ Р=0,95

2510 REM И ЧИСЛЕ СТЕПЕНЕЙ СВОБОДЫ K=N-1

2530 ДАТА 12.7,4.30,3.18,2.78,2.57,2.45,2.36,2.31,2.26,2.23,2.20,2.18

2540 ДАТА 2.16,2.14,2.13,2.12,2.11,2.10,2.09,2.08,2.07,2.06,2.06

2550 ДАТА 2.06,2.05,2.05,2.05,2.04

2560 POR M=1 TO 30

2570 READ T(M)

2580 NEXT M

2590 RETURN