№
|
Число
измерений n
|
Количество
интервалов r
|
1
|
40-100
|
7-9
|
2
|
100-500
|
8-12
|
3
|
500-1000
|
10-16
|
4
|
1000-10000
|
12-22
|
Б) Строят ступенчатое распределение результатов
измерений, откладывая по оси ординат среднее значение плотности распределения в
интервале ∆хi
(рисунок
1).
Площадь под гистограммой, являющейся
экспериментальной функцией плотности распределения вероятности результатов
измерений, должна ровняться единице.
Рисунок 1 - Гистограмма теоретическое
(нормальное) распределение результатов измерений.
В) По виду гистограммы подбирается
теоретическая кривая распределения. Если гистограмма подобна нормальному
распределению, то, вычисляя среднее арифметическое и дисперсию результатов измерений , строят теоретическую кривую
(рисунок 1)
;
Г) Проверяют выдвинутую гипотезу о
нормальности распределения результатов измерения, используя метод Пирсона.
Расхождение между экспериментальной гистограммой и выбранной теоретической
кривой в этом методе оценивается с помощью величины:
;
где - коэффициент весов разрядов,
выбранных для усиления веса составляющих с малой вероятностью в виде ;
- значение вероятности в i-м интервале
экспериментальной гистограммы;
- значение вероятности i-м
интервале, определенное по теоретической кривой как площадь под i-м
интервалом:
;
Мера расхождения - величина случайная, и плотность
вероятности распределения этой величины подчиняется распределению Пирсона при и (реально при ).
Заметим, что число степеней свободы
здесь k=r-3,
поскольку результаты измерений использованы для вычисления среднего, дисперсии
и общей площади под гистограммой.
Задаваясь уровнем значимости , находят по таблице 2 значение , соответствующее значению q или
вероятности P. Если , то распределение результатов
измерений принимают нормальным и гипотеза о нормальности распределения
считается верной. При этом возможны ошибки. В соответствии с математической
статистикой, если отвергается правильная гипотеза, то имеют дело с ошибкой 1-го
рода, а если принимается неверная гипотеза (распределение нормально, но это
отвергается в результате исследований), то имеет место ошибка 2-го рода.
Таблица 2 - Интегральная функция
распределения Пирсона. Значения для различных k и P
.2 Проверка с помощью составного
критерия
Вычисляется статистика:
;
квантили (квантиль - абсцисса,
соответствующая определенной вероятности) распределение которых приведены в
таблице 3.
Таблица 3 - Квантили распределения
статистики d
Если при данном числе измерений n и выбранном
уровне значимости q1 соблюдается
условие d1-0,5q< d ≤d0,5q то гипотеза
о нормальности распределения на основании первого критерия принимается, если -
нет, то отвергается. Гипотеза по второму критерию принимается, если не более m разностей измерений превышают уровень z0,5(1+P)Sx, где Sx - оценка
СКО результатов измерений, z0,5(1+P) - квантиль
интегральной функции нормированного нормального распределения, определяемая по
таблице 4 или таблице 5 при Ф(z0,5(1+P)) =
0,5(1+Р).
Таблица 4 - Интегральная функция
нормированного нормального распределения
Величина Р находится при заданном
уровне значимости q2 по данным
таблицы 6.
Таблица 6 - Значения n и a,
соответствующие различным m и q
Распределение результатов измерений
считается отличным от нормального, если не выполняется хотя бы один из
рассмотренных критериев. Уровень значимости составного критерия q ≤ q1 +q2 или
вероятность того, что распределение результатов измерений нормально, P = 1 - q.
.3 Применение коэффициентов
асимметрии и эксцесса для проверки нормальности распределения результатов
измерений
Известно, что для нормального
распределения коэффициенты асимметрии и эксцесса определяются по формулам:
;
;
Известно также, что дисперсии этих
величин равны:
;
;
Зная оценки дисперсий и , можно определить, значимо ли
выборочные коэффициенты асимметрии и эксцесса отличаются от нуля. Если
;
;
то наблюдаемое распределение можно
считать нормальным.
.4 Метод Колмогорова
Для применения критерия Колмогорова
необходимо определить наибольшее абсолютное отклонение выборочной функции распределения
от генеральной , то есть экспериментальной функции
интегрального распределения результатов измерения от теоретической функции,
которая предлагается для сравнения: . Затем вычисляется . Квантили распределения Колмогорова,
приведенные в таблице 7, сравниваются со значением л, вычисленным по
экспериментальным данным.
Таблица 7
Если , то гипотеза о совпадении
теоретического закона распределения с выборочным не отвергается. Про гипотеза отклоняется (или считается
сомнительной). Уровень значимости q при
применение критерия Колмогорова выбирают обычно равным 0,2 - 0,3.
.5 Непараметрический метод оценки
При малом числе измерений для
быстрой, но грубой оценки нормальности распределения результатов измерений
иногда используют так называемый непараметрический метод.
Для проверки нормальности
распределения результатов измерений вначале результаты располагают в
вариационный ряд в порядке их возрастания, а затем строится статистическая
интегральная функция распределения:
;
где ; - общее число измерений.
Функция - ступенчатая функция, скачки
которой равны. Если в этом ряду имеются одинаковые числа, то скачок возрастает
на количество этих чисел.
Затем теоретическую функцию
распределения приравнивают к статистической и для каждого по таблице 4 находят
соответствующие значения . Поскольку переменные и связаны соотношением , то зависимость при нормальном распределении должна
иметь вид прямой линии или, точнее говоря, точки с координатами и должны расположиться вдоль одной
прямой линии. Эта прямая по оси абсцисс отсекает отрезок, примерно равный
математическому ожиданию, а котангенс угла наклона прямой к оси абсцисс равен
примерно СКО ряда измерений. Если расположение точек существенно отличается от
прямолинейного, то распределения результатов измерений отличается от
нормального[1].
2. Расчеты
Рисунок 2 - Сбор данных
Таблица 8 - Периоды калибровочной
решетки “TGZ2”
№
|
Т,
мкм
|
№
|
Т,
мкм
|
№
|
Т,
мкм
|
№
|
Т,
мкм
|
1
|
2.49
|
26
|
2.66
|
51
|
2.54
|
76
|
2.58
|
2
|
2.45
|
27
|
2.58
|
52
|
2.58
|
77
|
2.58
|
3
|
2.45
|
28
|
2.54
|
53
|
2.58
|
78
|
2.66
|
4
|
2.45
|
29
|
2.54
|
54
|
2.58
|
2.58
|
5
|
2.49
|
30
|
2.62
|
55
|
2.58
|
80
|
2.54
|
6
|
2.45
|
31
|
2.58
|
56
|
2.54
|
81
|
2.54
|
7
|
2.5
|
32
|
2.62
|
57
|
2.58
|
82
|
2.54
|
8
|
2.45
|
33
|
2.62
|
58
|
2.54
|
83
|
2.58
|
9
|
2.41
|
34
|
2.58
|
59
|
2.58
|
84
|
2.54
|
10
|
2.45
|
35
|
2.58
|
60
|
2.58
|
85
|
2.62
|
11
|
2.49
|
36
|
2.54
|
61
|
2.54
|
86
|
2.62
|
12
|
2.49
|
37
|
2.66
|
62
|
2.54
|
87
|
2.54
|
13
|
2.49
|
38
|
2.66
|
63
|
88
|
2.54
|
14
|
2.58
|
39
|
2.58
|
64
|
2.54
|
89
|
2.58
|
15
|
2.54
|
40
|
2.58
|
65
|
2.54
|
90
|
2.62
|
16
|
2.49
|
41
|
2.58
|
66
|
2.5
|
91
|
2.62
|
17
|
2.5
|
42
|
2.54
|
67
|
2.58
|
92
|
2.58
|
18
|
2.5
|
43
|
2.58
|
68
|
2.5
|
93
|
2.58
|
19
|
2.41
|
44
|
2.58
|
69
|
2.54
|
94
|
2.58
|
20
|
2.45
|
45
|
2.58
|
70
|
2.58
|
95
|
2.58
|
21
|
2.49
|
46
|
2.54
|
71
|
2.58
|
96
|
2.58
|
22
|
2.45
|
47
|
2.58
|
2.58
|
97
|
2.58
|
23
|
2.49
|
48
|
2.62
|
73
|
2.58
|
98
|
2.66
|
24
|
2.45
|
49
|
2.62
|
74
|
2.5
|
99
|
2.58
|
25
|
2.58
|
50
|
2.54
|
75
|
2.62
|
100
|
2.62
|
.1 Грубая оценка нормальности распределения
периода калибровочной решетки “TGZ2”
непараметрическим методом
Располагаем в вариационный ряд период
калибровочной решетки в порядке их возрастания, а затем строим статистическую
интегральную функцию распределения:
;
где ; - общее число измерений. Если имеются
одинаковые значения, то k возрастает на количество этих
чисел.
Затем теоретическую функцию
распределения приравнивают к статистической и для каждого по таблице 4 находят
соответствующие значения .
Таблица 9
, мкмЧисло
одинаковых результатов
|
|
|
|
|
2.41
|
2
|
2
|
0.02
|
-2.4046
|
2.45
|
9
|
11
|
0.11
|
-1.2326
|
2.49
|
8
|
19
|
0.19
|
-0.8806
|
2.5
|
6
|
25
|
0.25
|
-0.6745
|
2.54
|
22
|
47
|
0.47
|
-0.0754
|
2.58
|
37
|
84
|
0.83
|
0.9585
|
2.62
|
11
|
95
|
1.2816
|
2.66
|
5
|
100
|
0.99
|
2.3267
|
Расчет:
;
;
;
;
;
;
;
.
Для определения нормальности
распределения периода калибровочной решетки “TGZ2”
непараметрическим методом построим график функции с координатами и .
График 1 - Зависимость от .
.2 Оценка нормальности распределения
периода калибровочной решетки “TGZ2” методом Колмогорова
Для применения критерия Колмогорова
определяю наибольшее абсолютное отклонение выборочной функции распределения от генеральной , то есть экспериментальной функции
интегрального распределения результатов измерения от теоретической функции,
которая предлагается для сравнения: .
Затем вычисляю . Квантили распределения Колмогорова,
приведенные в таблице 7, сравниваются со значением л, вычисленным по
экспериментальным данным.
Таблица 9
i
|
Границы
интервалов
|
mi
|
pi
|
npi
|
nFn(x)
|
nF(x)
|
|
1
|
2.41
- 2.46
|
11
|
0.0466
|
4.66
|
11
|
4.66
|
6.34
|
2
|
2.46
- 2.51
|
14
|
0.1756
|
17.56
|
25
|
22.22
|
2.78
|
3
|
2.51
- 2.56
|
22
|
0.3240
|
32.4
|
47
|
54.62
|
-7.62
|
4
|
2.56
- 2.61
|
37
|
0.2889
|
28.89
|
84
|
83.51
|
0.49
|
5
|
2.61
- 2.66
|
16
|
0.1289
|
12.89
|
100
|
96.4
|
3.6
А) Определяю ∑Т:
мкм.
Б) Определяю :
мкм.
В) Определяю ST:
мкм.
Г) Определяю параметр Ф:
Теперь по таблице 4 определяю
значение интегральной функции нормированного нормального распределения:
Ф1 (-2.4976) = 0.0062
Ф2 (-1.6212) = 0.0528
Ф3 (-0.7449) = 0.2284
Ф4 (0.1320) = 0.5524
Ф5 (1.0078) = 0.8413
Ф6 (1.8840) = 0.9702
Д) Определяю вероятность Pi:
;
;
;
;
;
.
И перемножаем вероятность Pi с
количеством измерений n:
;
;
;
;
.
E) Определяю
значения выборочной функции распределения и генеральной функции распределения
.
;
;
;
;
;
.
;
;
;
;
;
.
Ж) Определяю :
;
;
;
;
.
Данные занес в таблицу 9
Для принятия решения по критерию
Колмогорова определяю л:
.
По таблице 7 для уровня значимости q = 0.3
определяю л0.7 = 0.97.
3. Результаты и обсуждение
А) В непараметрическом методе для
грубой оценки нормальности распределения периода калибровочной решетки “TGZ2”использовался
график функции с координатами и .
Из графика 1 видно, что точки этой
зависимости довольно близко расположены к прямой линии, поэтому распределение
периода калибровочной решетки “TGZ2” можно признать нормальным.
Б) В методе Колмогорова для
определения нормальности распределения периода калибровочной решетки “TGZ2”
использовался критерий Колмогорова.
По полученным данным можно сказать,
что критерий Колмогорова позволяет считать распределение периода калибровочной
решетки “TGZ2”
нормальным, так как л = 0.634 < л0.7 = 0.97.
Заключение
В данной курсовой работе была
осуществлена проверка нормальности распределения результатов измерений периода
калибровочной решетки “TGZ2” с использованием непараметрического метода и
метода Колмогорова.
По итогам проверки, можно сказать,
что оба метода подтвердили нормальность распределения результатов измерений
периода калибровочной решетки “TGZ2”.
Список использованных источников
1. Пронкин Н.С., Основы
метрологии: практикум по метрологии и измерениям: учеб. пособие для вузов. -
М.: Логос, 2007. - 392 с.
Похожие работы на - Определение нормальности распределения методом Колмогорова
|