Обработка результатов равноточных многократных измерений.
1.1. Определить основные статистические показатели.
1.2. Построить гистограмму.
1.3. Построить теоретическую кривую нормального распределения на том же графике.
1.4. Определить доверительный интервал рассеивания случайных погрешностей вокруг среднего значения.
1.5. Определить суммарную погрешность обработки.
1.6. Ответить на вопросы: Для чего необходимо строить на одном графике теоретическую кривую нормального распределения и гистограмму? Какие выводы можно сделать по данному графику?
ВАРИАНТ 06. Обработка результатов равноточных многократных измерений с получением среднего арифметического X, среднеквадратичного отклонения Sx, и определением суммарной погрешности измерения в виде доверительного интервала ±∆Pд.
Исходные данные
Цена деления прибора С, мм 0,010
Результаты измерений, мм
1 80,090 12 80,150 23 80,130 34 80,090 45 80,120 56 80,120
2 80,030 13 80,130 24 80,170 35 80,050 46 80,090 57 80,190
3 80,080 14 80,110 25 80,090 36 80,030 47 80,010 58 80,070
4 80,170 15 80,190 26 80,030 37 80,160 48 80,070 59 80,090
5 80,070 16 80,100 27 80,110 38 80,080 49 80,050 60 80,110
6 80,130 17 79,970 28 80,050 39 80,230 50 80,150 61 80,090
7 80,090 18 79,990 29 80,110 40 80,130 51 80,110
8 80,110 19 80,070 30 80,130 41 80,110 52 80,150
9 80,060 20 80,230 31 80,150 42 80,110 53 80,050
10 80,130 21 80,130 32 80,210 43 80,090 54 80,100
11 80,070 22 80,070 33 80,130 44 80,110 55 80,170
Доверительная вероятность Рд = 0,95 – показывает вероятность нахождения истинного значения в рассчитанном интервале.
Уровень значимости q = 0,02 – показывающий, что принятый закон рассеивания размеров не будет соответствовать реальному закону.
Решение
1. Построение гистограммы.
Определяем величину размаха R (поле рассеяния):
R = Xmax – Xmin=80,230 -79,970=0,26 мм.
Xmax = 80,230 мм – наибольшее из измеренных значений;
Xmin = 79,970 мм – наименьшее из измеренных значений.
Определяем число интервалов разбиения n, в соответствии с рекомендациями:
n =N=61=7,81.
Количество интервалов принимается ближайшим большим нечетным.
Принимаем n = 7.
Определяем ширину интервала h:
h ==0,267=0,037 мм.
Определяем границы интервалов Xmini – Xmaxi
1 интервал: Xmin1 – Xmax1
Xmin1 = Xmin=79,970 мм
Xmax1 = Xmin1 + h = 79,970+0,037=80,007 мм
2 интервал: Xmin2 – Xmax2
Xmin2 = Xmax1 = 80,007 (мм)
Xmax2 = Xmin2 + h = 80,044 (мм)
3 интервал: Xmin3 – Xmax3
Xmin3 = Xmax2 = 80,044 (мм)
Xmax3 = Xmin3 + h = 80,081 (мм)
4 интервал: Xmin4 – Xmax4
Xmin4 = Xmax3 = 80,081 (мм)
Xmax4 = Xmin4 + h = 80,118 (мм)
5 интервал: Xmin5 – Xmax5
Xmin5 = Xmax4 = 80,118 (мм)
Xmax5 = Xmin5 + h = 80,155 (мм)
6 интервал: Xmin6 – Xmax6
Xmin6 = Xmax5 = 80,155 (мм)
Xmax6 = Xmin6 + h = 80,192 (мм)
7 интервал: Xmin7 – Xmax7
Xmin7 = Xmax6 = 80,192 (мм)
Xmax7 = Xmin7 + h = 80,229~80,230 мм.
Определяем середины интервалов Xoi
1 интервал: Xo1 = Xmin1 + =79,9885 (мм)
2 интервал:Xo2 = Xmin2 + = 80,0255 (мм)
3 интервал:Xo3 = Xmin3 + = 80,0625 (мм)
4 интервал:Xo4 = Xmin4 + = 80,0995 (мм)
5 интервал:Xo5 = Xmin5 + = 80,1365 (мм)
6 интервал:Xo6 = Xmin6 + = 80,1735 (мм)
7 интервал:Xo7 = Xmin7 + = 80,211 (мм)
Определение количества размеров попадающих в каждый интервал mi.
Используя заданную выборку, подсчитываем количество размеров попадающих в каждый интервал (если размер совпадает с границей интервала то его относят в интервал, находящийся слева по числовой оси)
Результаты выполненных выше расчетов занесем в таблицу:
Таблица 1.1
Номер интервала Границы интервала Середина интервала Xoi (мм) Число размеров в интервале, mi
Xmin (мм) Xmax (мм)
1 79,97 80,007 79,9885 2 0,033
2 80,007 80,044 80,0255 4 0,066
3 80,044 80,081 80,0625 13 0,213
4 80,081 80,118 80,0995 19 0,311
5 80,118 80,155 80,1365 14 0,230
6 80,155 80,192 80,1735 6 0,098
7 80,192 80,23 80,211 3 0,049
Используя табличные данные, строим гистограмму рассеивания единичных замеров и теоретическую кривую нормального распределения:
3387090521970001504950501650X=80,104 мм
00X=80,104 мм
158686577089000
Рисунок 1.1
2
. Проверка выборки на соответствие нормальному закону распределения
При числе измерений свыше 50 проверка распределения на соответствие нормальному закону может выполняться по критерию Пирсона. При использовании этого критерия определяется параметр хи-квадрат по следующей формуле:
,
где Noi – теоретическая частота попадания в интервал.
Теоретическая частота попадания в интервал определяется по формуле:
φ(z) – плотность вероятности появления размеров в каждом интервале;
σx – среднеквадратичное отклонение размеров (СКО) выборки.
Считая, что СКО практически совпадает с его оценкой (σx ≈ Sx) приведем формулу, по которой определяется оценка СКО:
В данную формулу входит величина , которая представляет среднеарифметическое значение измеряемой величины и определяется по формуле:
X=161×(79,9885×2+80,0255×4+80,0625×13+80,0995×19++80,1365×14+80,1735×6+80,211×3).
80,104 мм.
После подстановки 80,104 мм получим численные значения среднеарифметического и оценки СКО:
Sx=2*79,9885-80,1042+…+3*80,211-80,104261=0,1522261.
Sx=0,05 мм.
Кроме полученных величин, для определения теоретической частоты попадания в интервал Noi необходимо знать плотность вероятности попадания размеров в каждом интервале.
Эту величину можно определить по формуле: ϕ(z)=e(-z22)2π
Так как расчеты по данной формуле достаточно сложны, значения плотности вероятности выбирают из таблицы в зависимости от безразмерного параметра Z, который для каждого интервала определяется по формуле: Zoi=Xoi-XSx.
Для 1 интервала:
Zo1 =-2,31, что соответствует величине φ(z) = 0,028.
Для 2 интервала:
Zo2 = -1,57, что соответствует величине φ(z) = 0,116.
Для 3 интервала:
Zo3 = -0,83, что соответствует величине φ(z) = 0,283.
Для 4 интервала:
Zo4 = -0,09, что соответствует величине φ(z) = 0,397.
Для 5 интервала:
Zo5 = 0,65, что соответствует величине φ(z) = 0,322.
Для 6 интервала:
Zo6 = 1,39, что соответствует величине φ(z) = 0,152.
Для 7 интервала:
Zo7 = 2,184, что соответствует величине φ(z) = 0,040.
Определяем теоретические значения количества деталей для каждого интервала Noi.
Noi=N×hSxϕZoi.
Для 1 интервала:
No1 =1,250.
Для 2 интервала:
No2 = 5,251.
Для 3 интервала:
No3 = 12,761.
Для 4 интервала:
No4 =17,935.
Для 5 интервала:
No5 = 14,579.
Для 6 интервала:
No6 = 6,854.
Для 7 интервала:
No7 = 1,824.
На основании результатов измерений и расчета теоретических данных определяем фактическую и теоретическую частоту попадания размеров в интервал:
Таблица 1.2
№ интервала Фактическая частота Теоретическая частота
1 0,033 0,02
2 0,066 0,086
3 0,213 0,209
4 0,311 0,294
5 0,230 0,239
6 0,098 0,112
7 0,049 0,03
Полученные результаты позволяют получить расчетную величину параметра χ-квадрат:
Таблица 1.3
№
интервала mi
Noi
mi-Noi
(mi-Noi)2
(mi-Noi)2Noi
1 2 1,25 0,75 0,5625 0,45
2 4 5,251 -1,251 1,565001 0,298039
3 13 12,761 0,239 0,057121 0,004476
4 19 17,935 1,065 1,134225 0,063241
5 14 14,579 -0,579 0,335241 0,022995
6 6 6,854 -0,854 0,729316 0,106407
7 3 1,824 1,176 1,382976 0,758211
χ2=(mi-Noi)2Noi=1,703.
Для совпадения фактического закона распределения с теоретическим законом нормального распределения необходимо, чтобы выполнялось следующее условие: ,
где – теоретическое граничное значение параметра хи-квадрат, которое определяется по таблице (таблица 2 задания к контрольной работе).
Для получения табличного значения необходимо определиться с двумя параметрами:
- уровнем значимости q, который показывает вероятность того, что законы не совпадут.
В нашем случае, в соответствии с заданием, q = 0,02;
- числом степеней свободы , которое определяется в зависимости от числа интервалов n и числа определяемых по статистике параметров, необходимых для совмещения модели и гистограммы r.
Для нормального закона распределения r = 2, так как закон однозначно характеризуется двумя параметрами – СКО и МО (математическим ожиданием).
Число степеней свободы определяется по формуле:
Таким образом, табличное значение χq2=11,67.
В силу того, что выполняется условие
χ2=1,703<χq2=11,67,
делаем вывод о том, что исходные данные соответствуют нормальному закону распределения.
3