В тех случаях, когда наблюдаемые значения случайной величины задаются многозначными числами и объем выборки достаточно велик (n > 25), вначале целесообразно найти среднюю арифметическую по формуле а за тем перейти к вычислению центральных моментов порядка k (k = 2, 3, 4):
Интервалы
наблюдаемых значений СВ Х, МПа
Середины интервалов xi
Частоты mi
[18;19)
18,5
1
-4,44
19,71
-87,53
388,63
[19;20)
19,5
9
-30,96
106,50
-366,37
1260,31
[20;21)
20,5
20
-48,80
119,07
-290,54
708,91
[21;22)
21,5
41
-59,04
85,02
-122,43
176,29
[22;23)
22,5
56
-24,64
10,84
-4,77
2,10
[23;24)
23,5
60
33,60
18,82
10,54
5,90
[24;25)
24,5
38
59,28
92,48
144,26
225,05
[25;26)
25,5
16
40,96
104,86
268,44
687,19
[26;27)
26,5
7
24,92
88,72
315,83
1124,34
[27;28]
27,5
2
9,12
41,59
189,64
864,75
Итого
250
0
687,61
57,07
5443,47
Следовательно,
Для предварительного выбора закона распределения вычислим вначале средние квадратические ошибки определения асимметрии
и эксцесса
Критерием «нормальности» распределения прочности бетона на сжатие является равенство нулю асимметрии и эксцесса. Из приведенных расчетов видно, что выборочные коэффициенты асимметрии и эксцесса Э отличаются от нуля не более чем на удвоенные средние квадратические ошибки их определения, что соответствует нормальному распределению. Вид полигона и гистограммы частостей также напоминает нормальную кривую (кривую Гаусса).
Можно предположить, прочность бетона на сжатие (СВ Х) изменяется под влиянием большого числа факторов, примерно равнозначных по силе. Поэтому, исходя из «технологии» образования СВ Х, т. е. механизма образования отклонений прочности от некоторого номинального значения, можно предположить, что распределение прочности бетона на сжатие является нормальным.
Плотность вероятности нормального распределения имеет вид
Найдём точечные оценки параметров a и σ нормального распределения методом моментов:
Следовательно, плотность вероятности предполагаемого нормального распределения имеет вид
Функция распределения предполагаемого нормального распределения имеет вид
Используя нормированную функцию Лапласа , функцию нормального распределения можно записать в виде
Проведем проверку гипотезы о нормальном распределении СВ Х (прочности бетона на сжатие) с помощью критерия согласия для этого интервалы наблюдаемых значений нормируют, т.е. выражают их в единицах среднего квадратического отклонения s: , причем наименьшее значение полагают равным , наибольшее . Далее вычисляют вероятности попадания СВ Х, имеющей нормальное распределение, с параметрами а = 22,94, σ = 1,65 в частичные интервалы (хi-1; хi) по формуле
,
где
.
Например, вероятность того, что СВ Х (прочность бетона на сжатие) попадает в первый частичный интервал (;19) , равна
Аналогично
и т. д. После этого вычисляют теоретические (модельные) частоты нормального распределения и наблюдаемое значение критерия
Затем по таблицам квантилей распределения по уровню значимости q = 0,05 и числу степеней свободы ‚ (k — число интервалов; r — число параметров предполагаемого распределения СВ Х) находят критическое значение .
Если , то считают, что нет оснований для отклонения гипотезы о нормальном распределении прочности бетона на сжатие.
В противном случае, т. е. если , считается, что гипотеза нормального распределения прочности бетона на сжатие не согласуется с экспериментальными данными.
Вычисления, необходимые для определения наблюдаемого значения выборочной статистики приведем в таблице:
Интервалы наблюдаемых значений СВ Х, МПа
Нормированные интервалы [ui, ui-1]
pi
npi
(-∞;-2,39)
0,008
2,00
0,05
[-2,39;-1,78)
0,029
7,25
3,06
0,42
[-1,78;-1,18)
0,081
20,25
0,06
0,00
[-1,18;-0,57)
0,168
42,00
0,02
[-0,57;0,04)
0,231
57,75
[0,04;0,64)
0,223
55,75
18,06
0,32
[0,64;1,25)
0,154
38,50
0,25
0,01
[1,25;1,85)
0,074
18,50
6,25
0,34
[1,85;2,46)
0,025
0,56
0,09
[2,46;+∞)
0,007
1,75
0,03
1.000
250,0
Замечание. Наименьшее значение стандартизованной переменной
заменено , наибольшее значение заменено . Эта замена произведена для того, чтобы сумма теоретических (модельных) частот npi была равна объему выборки.
В результате вычислений получили . Найдем по таблице квантилей распределения по уровню значимости α = 0,05 и числу степеней критическое значение . Так как , то нет оснований для отклонения гипотезы о нормальном распределении прочности бетона на сжатие.
Построим нормальную кривую. Для этого из середин частичных интервалов восстании перпендикуляры высотой pi/h (pi — вероятность попадания СВ Х в частичный интервал; h — длина интервала). На рисунке 3 концы этих перпендикуляров отмечены кружками. Полученные точки соединены плавной кривой. Сравнение гистограммы и нормальной кривой наглядно показывает, что нормальная кривая хорошо сглаживает гистограмму относительных частот.
Найдем интервальные оценки параметров нормального распределения. Для вычисления доверительного интервала накрывающего математическое ожидание прочности бетона на сжатие (СВ Х), найдем по таблицам квантилей распределения Стьюдента по заданной доверительной вероятности и числу степеней свободы‚ квантиль.
Вычислим предельную погрешность интервального оценивания
Искомый доверительный интервал для математического ожидания
Смысл полученного результата: если будет произведено достаточно большое число выборок по 250 исследований прочности образцов бетона на сжатие, то в 95% из них доверительный интервал накроет математическое ожидание прочности бетона и только в 5% случаев математическое ожидание может выйти за границы доверительного интервала.
Для нахождения доверительного интервала, накрывающего неизвестное среднее квадратическое отклонение σ с заданной вероятностью , найдем по доверительной вероятности и числу степеней свободы ‚ два числа; . Искомый доверительный интервал
Полученный результат означает, что если будет произведено достаточно большое число выборок по 250 исследований прочности образцов бетона на сжатие, то в 95% из них доверительный интервал накроет среднее квадратическое отклонение σ и только в 5% среднее квадратическое отклонение σ можёт выйти за границы доверительного интервала.
Страницы: 1, 2
