Главная Основа метрологии Проверка нормальности распределения результатов наблюдений


Проверка нормальности распределения результатов наблюдений

Метрология - Основа метрологии

В предыдущих разделах было показано, что результаты наблюдений можно оценить наиболее полно, если их распределение является нормальным. Поэтому исключительно важную роль при обработке результатов наблюдений играет проверка нормальности распределения.

Эта задача представляет собой частный случай более общей проблемы, заключающейся в подборе теоретической функции распределения, в некотором смысле наилучшим образом согласующейся с опытными данными.

При большом числе результатов наблюдений (n>40) данная задача решается в следующем порядке.

Весь диапазон полученных результатов наблюдений Xmax ... Xmin разделяют на r интервалов шириной и подсчитывают частоты mi, равные числу результатов, лежащих в каждом i-м интервале, т. е. меньших или равных его правой и больших левой границы.

Отношения

(50)


где n - общее число наблюдений, называются частостями и представляют собой статистические оценки вероятностей попадания результата наблюдений в i-й интервал. Распределение частот по интервалам образует статистическое распределение результатов наблюдений.

Если теперь разделить частость на длину интервала, то получим величины

(51)


являющиеся оценками средней плотности распределения в интервале .

Отложим вдоль оси результатов наблюдений (рис.11) интервалы в порядке возрастания индекса i и на каждом интервале построим прямоугольник с высотой, равной . Полученный график называется гистограммой статистического распределения.

Площадь суммы всех прямоугольников равна единице:

При увеличении числа наблюдений число интервалов можно увеличить. Сами интервалы уменьшаются, и гистограмма все больше приближается к плавной кривой, ограничивающей единичную площадь, - к графику плотности распределения результатов наблюдений.

При построении гистограмм рекомендуется пользоваться следующими правилами:

1. Число интервалов выбирается в зависимости от числа наблюдений согласно рекомендациям табл.6.

Таблица 6

n

r

40 – 100

7 – 9

100 – 500

8 – 12

500 – 1000

10 – 16

1000 – 10000

12 – 22

2. Длины интервалов удобнее выбирать одинаковыми. Однако если распределение крайне неравномерно, то в области максимальной концентрации результатов наблюдений следует выбирать более узкие интервалы.

3. Масштабы по осям гистограммы должны быть такими, чтобы отношение ее высоты к основанию составляло примерно 5:8.

Пример. Было выполнено 100 измерений среднего диаметра резьбового калибра. Результаты наблюдений лежат в диапазоне 8.911 - 8.927 мм, т. е. зона распределения результатов составляет 0.016 мм. Весь диапазон удобно разделить на восемь равных интервалов через 0.002 мм. В табл.7 приведены частоты mi, частости и плотности статистического распределения.

Таблица 7

, мм

, мм

, 1/мм

1

8.911

8.913

1

0.01

5

2

8.913

8.915

5

0.05

25

3

8.915

8.917

14

0.14

70

4

8.917

8.919

27

0.27

135

5

8.919

8.921

24

0.24

120

6

8.921

8.923

18

0.18

90

7

8.923

8.925

9

0.09

45

8

8.925

8.927

2

0.02

10

После построения гистограммы надо подобрать теоретическую плавную кривую распределения, которая, выражая все существеные черты статистического распределения, сглаживала бы все случайности, связанные с недостаточным объемом экспериментальных данных. Принципиальный вид теоретической кривой выбирают заранее, проанализировав метод измерения, или хотя бы по внешнему виду гистограммы. Тогда определение аналитического вида кривой распределения сводится к выбору таких значений его параметров, при которых достигается наибольшее соответствие между теоретическим и статистическим распределением. Одним из методов решения этой задачи является метод моментов. При его использовании параметрам теоретического распределения придают такие значения, при которых несколько важнейших моментов совпадают с их статистическими оценками. Так, если статистическое распределение, определяемое гистограммой, приведенной на рис.11, мы хотим описать кривой нормального распределения, то естественно потребовать, чтобы математическое ожидание и дисперсия последнего совпадали со средним арифметическим и оценкой дисперсий, вычисленным по опытным данным. В предыдущем примере мм, мм и уравнение кривой нормального распределения, лучше всего согласующегося со статистическим распределением, должно иметь вид:

Далее законно возникает вопрос, объясняются ли расхождения между гистограммой и подобранным теоретическим распределением только случайными обстоятельствами, связанными с ограниченным числом наблюдений, или они вызваны тем, что результаты наблюдений в действительности распределены иначе?

Для ответа на этот вопрос используют методы проверки статистических гипотез. Идея их применения заключается в следующем. На основании гистограммы, полученной при обработке опытных данных, строится гипотеза, состоящая в том, что результаты наблюдений подчиняются распределению с плотностью .

Для того чтобы принять или опровергнуть эту гипотезу, выбирается некоторая величина U, представляющая собой меру расхождения теоретического и статистического распределений. В качестве меры расхождения можно принять сумму квадратов разностей частостей и теоретических вероятностей попадания результатов наблюдений в каждый интервал, взятых с некоторыми коэффициентами:

,

(52)


где – коэффициенты, называемые весами разрядов; – теоретические вероятности, определяемые как

,

(53)


Здесь – предполагаемая плотность распределения.

Мера расхождения U является случайной величиной и, независимо от исходного распределения подчиняется - распределению с k степенями свободы – см. формулу (44). Если значения всех частот , число измерений стремится к бесконечности, а веса выбираются равными . Число степеней свободы распределения k = r - s, где – число разрядов гистограммы статистического распределения, а s – число независимых связей, наложенных на частости .

Если проверяется гипотеза о нормальности распределения, то к числу этих связей относится равенство среднего арифметического математическому ожиданию, а точечной оценки дисперсии - дисперсии предполагаемого нормального распределения. Кроме того, всегда требуется, чтобы сумма частостей по всем интервалам была равна единице. Поэтому в данном случае s = 3.

По табл.П.6 можно при заданной доверительной вероятности найти тот доверительный интервал значений , в который мера расхождения может попасть по чисто случайным причинам.

Если вычисленная по опытным данным мера расхождения окажется в указанном интервале, то гипотеза принимается. Это, конечно, не значит, что гипотеза верна. Можно лишь утверждать, что она правдоподобна, т.е. не противоречит опытным данным. Если же она выходит за границы доверительного интервала, то гипотеза отвергается как противоречащая опытным данным.

Поскольку проверка гипотезы основывается на опытных данных, то при принятии решения всегда возможны ошибки. Отвергая в действительности верную гипотезу, мы совершаем ошибку первого рода. Вероятность ошибки первого рода называется уровнем значимости и составляет . Принимая в действительности неверную гипотезу, мы совершаем ошибку второго рода. Вычислить ее вероятность, вообще говоря, невозможно, поскольку для этого нужно рассмотреть все прочие возможные гипотезы, являющиеся альтернативой обсуждаемой гипотезы. Можно лишь утверждать, что при уменьшении ошибки первого рода ошибка второго рода увеличивается, поэтому не имеет смысла брать слишком высокие значения доверительных вероятностей.

Описанная процедура проверки гипотезы о том, что данное статистическое распределение является распределением с плотностью , называется критерием согласия . Проверка нормальности распределения согласно критерию сводится к следующему.

1. Данные наблюдений группируют по интервалам, как при построении гистограммы, и подсчитывают частоты . Если в некоторые интервалы попадает меньше пяти наблюдений, то такие интервалы объединяют с соседними. При этом число степеней свободы k, конечно, уменьшается.

2. Вычисляют среднее арифметическое и точечную оценку среднеквадратического отклонения результата наблюдений , которые принимают в качестве параметров теоретического нормального распределения с плотностью .

3. Для каждого интервала находят вероятности попадания в них результатов наблюдений либо по общей формуле (29), либо приближенно как произведение плотности теоретического распределения в середине интервала на его длину:

.

(54)

4. Для каждого интервала вычисляют величины и суммируют их по всем , в результате чего получают меру расхождения .

5. Определяют число степеней свободы и, задаваясь уровнем значимости , находят по табл.П.6 приложения значения и . Если , то распределение результатов наблюдений считают нормальным.

Критерий согласия , построенный на предельном переходе при , рекомендуется применять, если общее число наблюдений больше сорока.

При малом числе наблюдений нормальность распределения результатов наблюдений проверяется с помощью двух критериев.

Первый критерий основан на вычислении статистики

.

(55)

Гипотеза о нормальности распределения на основании первого критерия принимается, если при данном числе наблюдений и выбранном уровне значимости соблюдается условие

,


где и – квантили, выбираемые из табл.П.8.

На основании второго критерия гипотеза о нормальности распределения принимается, если не более разностей превосходят уровень , где – оценка среднеквадратического отклонения результатов наблюдения, – квантиль интегральной функции нормированного нормального распределения, определяемый по данным табл.П.2 приложения при значении

.

Величина находится при заданном уровне значимости второго критерия по данным табл.П.9.

Распределение результатов наблюдения считается отличным от нормального, если оно не соответствует хотя бы одному из этих двух критериев. Уровень значимости составного критерия

.

При малом числе наблюдений для оценки нормальности можно воспользоваться понятием статистической функции распределения результатов наблюдений. Для ее построения полученные в процессе эксперимента результаты группируют в так называемый вариационный ряд , члены которого располагаются в порядке их возрастания, так что всегда . Статистическую функцию распределения определяют по формуле

(56)

представляет собой ступенчатую линию, скачки которой соответствуют значениям членов вариационного ряда. Каждый скачок равен , если все членов ряда различны. Если же для некоторого , то в точке возрастает на , где i – число равных между собой членов ряда.

Если число наблюдений безгранично увеличивать, то статистическая функция распределения сходится по вероятности к истинной функции .

Для проверки нормальности распределения результатов наблюдений по табл.3 приложения находят значения , соответствующие полученным значениям статистической функции распределения . Но переменная определяется через результаты наблюдений как


и если в координатах нанести точки , то при нормальном распределении они должны расположиться вдоль одной прямой линии. Если же в результате такого построения получится некоторая кривая линия, то гипотезу о нормальности распределения придется отвергнуть как противоречащую опытным данным.

Пример. Даны результаты девятнадцати измерений длины детали (см. табл.3). Проверить нормальность распределения результатов наблюдений.

Вычисления по изложенной методике сведены в табл.8.

Таблица 8

, мм

18.303

0.05

-1.6449

18.304

0.10

-1.2816

18.305

0.20

-0.8416

18.306

0.30

-0.52.44

18.307

0.40

-0.2533

18.308

0.60

0.2533

18.309

0.75

0.6745

18.310

0.85

1.0364

18.311

0.90

1.2816

18.312

0.95

1.6449

На рис.12 представлена зависимость . Отдельные точки располагаются очень близко к прямой, поэтому распределение результатов наблюдений можно считать нормальным.

 


Читайте:


Добавить комментарий


Защитный код
Обновить

Учёные первооткрыватели:

Геккель, Эрнст Генрих

News image

Эрнст Ге нрих Фили пп А вгуст Ге ккель (нем. Ernst Heinrich Philipp August Haeckel; 16 февраля 1834, Потсдам — 9 ...

ДЖОН ФЛЕМСТИД

News image

Джон Флемстид родился 9 (19) августа 1646 г. в г. Дерби графства Дербишир (Derbyshire) в северной Англии. Заинтересовавшись астрономией, он на...

Авторизация



Единицы измерений:

Гигабайт

News image

Гигабайт  (Гбайт, Г, ГБ) — кратная единица измерения количества информации, равная 109 стандартных (8-битным) байтов или 1000 мегабайтам. Неправильность названия Читая нижеизложенный те...

Единицы измерения количества информации

News image

Единицы измерения информации служат для измерения объёма информации — величины, исчисляемой логарифмически. Это означает, что когда несколько объектов рассматриваются как од...

Ом

News image

Ом (обозначение: Ом, Ω) — единица измерения электрического сопротивления в СИ. Ом равен электрическому сопротивлению проводника, между концами которого возникает на...

Атмосфера (единица измерения)

News image

Атмосфера — внесистемная единица измерения давления, приблизительно равная атмосферному давлению на поверхности Земли на уровне Мирового океана. Существуют две примерно равные др...

Открыватели:

Вавилов, Сергей Иванович

News image

Серге й Ива нович Вави лов (12 (24) марта 1891, Москва — 25 января 1951, Москва) — советский физик, академик (1932), основатель научной школы физической оп...

Универсальный конвертер
Conversion Type:
Quantity:

converts to:

Construction Unit converter provided by: EcoLog Homes

Интересные факты:

Таблица Менделеева

News image

В конце августа 1875 г. в кабинет акад. Вюрца входит его ученик, молодой французский химик Лекок-де-Буабодран. н долго не решается об...

О звуке

News image

Звук с давних пор считался одним из самых загадочных явлений природы. В самом деле, что порождает звук? Что заставляет его не...

Эйнштейн и квантовая теория света

News image

Эйнштейн является одним из основателей новой, квантовой теории света и основателем теории относительности. Согласно квантовой теории свет представляет поток своеобразных ча...

Как происходит кристаллизация жидкости

News image

В настоящее время можно считать твердо установленным, что жидкость может затвердевать после ее охлаждения до температуры плавления только при наличии в ...

Атом и время

News image

Трудно себе представить более простое и вместе с тем более сложное понятие, чем время. Старая пословица говорит: «нет ничего в ми...

Ньютон и Марат о притяжении лучей света

News image

Что такое свет?— На этот вопрос Ньютон, очень много поработавший над изуче­нием световых явлений, отвечал так: свет — это поток бы...