Законы распределения случайных погрешностей

(на примере математического маятника).

Цель работы: изучить основные законы распределения случайных погрешностей.

Введение.

Измерить физическую величину абсолютно точно принципиально невозможно по следующим причинам:

а) несовершенство приборов;

б) влияние внешних условий на характеристики приборов;

в) влияние внешних факторов на параметры изучаемого объекта;

г) субъективные человеческие факторы.

Учитывая все выше изложенное, на практике возможно указать лишь диапазон значений, в который попадает истинное значение искомой величины с определенной вероятностью:

х ± Dх,

где Dх – абсолютная ошибка измерений.

Этот диапазон называют доверительным интервалом с соответствующей надежностью a. Абсолютную ошибку измерений определяют как разницу между истинным и измеренным значением определяемой величины. Но, так как, истинное значение не известно, то существуют разные способы оценки абсолютной погрешности. Для бесконечно большого числа измерений применяют следующие способы (на практике они применяются при условии n³ 30):

1) в качестве абсолютно погрешности Dх принимают среднеарифметическую погрешность r:

с надежностью a=0, 57, (1)

где: – средняя величина n измеренных значений искомой величины,

x_i – величина каждого измерения;

2) в качестве абсолютной погрешности Dх принимают среднеквадратичную погрешность s:

, с надежностью a=0, 68; (2)

3) в качестве абсолютной погрешности Dх принимают вероятную погрешность h:

Dх=h, с надежностью a=0, 5.

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

Для r, s и h существуют теоретические соотношения:

; ; (3)

Все абсолютные погрешности разделяют на две большие группы:

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

2. Случайные погрешности. Случайными называются погрешности, величина и знак которых изменяются при повторении измерений.

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

Таблица 1

	-0, 047	-0, 024	-0, 009		0, 013	0, 027
	-0, 041	-0, 021	-0, 007	0, 002	0, 014	0, 032
	-0, 039	-0, 018	-0, 006	0, 004	0, 017	0, 034
	-0, 034	-0, 016	-0, 004	0, 006	0, 019	0, 039
	-0, 032	-0, 014	-0, 003	0, 009	0, 021	0, 041
	-0, 026	-0, 013	-0, 001	0, 011	0, 023	0, 048
s=0, 025

Выбирается система координат DT0DN, причем ось DT строится как в положительном, так и в отрицательном направлении. Ось DT разбивается на отрезки равной длины, т.е. с одинаковым шагом. В качестве размера шага, как правило, выбирается половина разряда первой значащей цифры границы доверительного интервала рассчитанного по формуле (2). Таким образом, для данного примера шаг оси DT равен 0, 005с (в данном примере первая значащая цифра s равна 2, ее разряд равен 0, 01). Ось DN имеет единичный отрезок равный 1. По этой оси откладывается количество значений погрешностей, которые попадают в каждый выбранный интервал на оси DT.

Интервалы значений по оси DT выбираются равными величине разряда первой значащей цифры границы доверительного интервала рассчитанного по формуле (2). Для рассматриваемого примера величина интервала будет равна 0, 01. Первый интервал выбирается симметрично относительно оси DN. Таким образом, для данного случая первый интервал будет от –0, 005 до 0, 005. Из таблицы 1 следует, что в этот интервал попадает 6 значений погрешностей. Поэтому на этом шаге строится прямоугольник шириной равной выбранному интервалу и высотой в шесть единиц. В следующий интервал от 0, 005 до 0, 015 попадает 5 значения погрешностей, т.е. на этом интервале строится прямоугольник высотой в 5 единиц и т.д. Аналогично, строятся прямоугольники для отрицательных значений абсолютных погрешностей. Пример построенной гистограммы приведен на рис.1.

Рис.1. Гистограмма

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

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

, (4)

где f(Dx_i) – плотность вероятности появления ошибки Dх_i, рассчитываемая по формуле:

(5)

Если по оси Y откладывать не значения DN, а выражение , то получим кривую, которая называется графиком распределения Гаусса. На рисунке 2 показаны графики распределения Гаусса с разной дисперсией.

Рис.2

Анализируя функцию распределения Гаусса и гистограмму на рисунках 1, 2 можно установить основные закономерности распределения случайных погрешностей в зависимости от их величины:

1. При большом числе измерений случайные ошибки принимают непрерывный ряд значений. Это подтверждает сплошная линия на рис.1.

2. Ошибки одинаковые по величине, но разного знака встречаются одинаково часто. Из гистограммы рис.1 видно, количество ошибок в интервале от 0, 005 до 0, 015 на положительной ветви оси DТ равно количеству ошибок в таком же интервале на отрицательной ветви оси DТ.

3. Частота появления ошибок уменьшается с увеличением их величины. Из гистограммы рис.1 видно, высота прямоугольников убывает вдоль оси DТ как в положительную, так и в отрицательную сторону.

Функция распределения случайных погрешностей Гаусса (6) позволяет установить физический смысл понятия «доверительная вероятность» или «надежность».

Если в соответствии с данными рис.1 доверительный интервал выберем от –0, 02 до 0, 02, то доверительная вероятность вычисляется путем интегрирования функции (4) в пределах выбранного интервала погрешностей:

. (6)

Тогда ответ можно будет представить в виде:

Т_ист= ±0, 02 с надежностью a = Р.

Если интервал взять «шире» от –0, 05 до 0, 05, то доверительная вероятность, посчитанная по формуле (6), возрастет. Это и понятно: чем шире доверительный интервал, тем вероятнее истинное значение периода колебаний попадет в этот интервал. Если интеграл (6) посчитать в пределах от –µ до +µ, то получим цифру 1. Это означает, истинное значение периода колебаний абсолютно достоверно попадет в интервал от –µ до +µ. Равенство интеграла (6) единице называется условием нормировки вероятностей.

Функция распределения Гаусса в виде, представленном формулой (4) применима на бесконечно малом интервале dT. В этой работе интервалы DT нельзя считать бесконечно малыми, поэтому следует вести речь о дискретных значениях функции распределения Гаусса определенных для границ рассматриваемых интервалов.

Приравняем выражения (4) и (5):

Выразим DN и перейдем к дискретным значениям:

, (7)

где Dх – величина интервала, DT – граница интервала.

Систематические погрешности смещают гистограмму вправо или влево по оси DТ, не искажая характер распределения случайных погрешностей.

Описание установки

Установка включает в себя модель математического маятника 1 и электронный блок 2 для измерений периода колебаний (рис.3). Электронный блок измеряет число периодов (индикатор 3) и время колебаний (индикатор 4). На блоке имеются кнопки управления: «СБРОС»(5) – очищает индикаторы от предыдущих показаний; «ПУСК»(6) – запускает счетчик времени; «СТОП»(7) – останавливает счет времени на заданном числе периодов.

Рис.3 Внешний вид установки

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

, (8)

где L – длина математического маятника.

123	Поделиться: