Элементы теории ошибок (погрешностей)

Цельюлюбогоизмерения некоторой физической величины является получение её истинного значения. Однако это весьма непростая задача из-за различных ошибок (погрешностей), неизбежно возникающих при измерениях.

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

Ошибки (погрешности) измерений принято делить на систематические и случайные.

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

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

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

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

Допустим, измеряется величина Х и мы хотим найти её истинное значение – х_ист. Результатом n измерений, проведенных соответствующим прибором, является ряд её значений: х₁, х₂, х₃ ,…, х_n.

Разность между полученным х_i и истинным х_ист значениями представляет собой случайную абсолютную погрешность отдельного измерения D х_{i =} х_i- х_ист. Причём из теории ошибок следует, что при большом числе измерений (большом n) ошибки одной и той же величины, но разного знака встречаются одинаково часто. Посмотрим, к чему это приводит. Представим полученные нами значения х_i через х_ист и D х_i и сложим получившиеся соотношения:

х₁ = х_ист. + D х₁;

х₂ = х_ист. + D х₂;

……………….

х_n = х_ист. + D х_n;

_____________

= nx_ист. + .

Отсюда найдем истинное значение измеряемой величины:

x_ист = - .

Поскольку при большом числе измерений n ошибки равные по величине, но разные по знаку встречаются одинаково часто, то сумма абсолютных ошибок не растет с увеличением n, а лишь колеблется вблизи нуля, поэтому с увеличением n слагаемое уменьшается и стремится к нулю при n ® ¥. Следовательно, при очень большом количестве измеренийистинное значение измеряемой величины практически совпадает со средним арифметическим всех полученных значений:

x_ист = = .

Однако при любом ограниченном количестве проведенных измерений n истинное значение х_ист будет отличаться от найденного среднего арифметического значения – х ¹ х_ист. –, необходимо оценить величину этого различия.

К решению данного вопроса можно подойти следующим образом. В связи с влиянием случайных ошибок на результаты измерений некоторой физической величины Х ряд полученных в эксперименте её значений х₁, х₂, х₃ …, х_n можно рассматривать как выборку из генеральной совокупности, которой соответствует
n ® ¥ и математическое ожидание которой – М_г(Х) = _г = х_ист. – надо найти (предполагается, и теория ошибок это подтверждает, что результаты измерений в генеральной совокупности распределены по нормальному закону).

Полученной выборке, естественно, соответствует свое среднее арифметическое значение:

= .

Тогда с определенной доверительной вероятностьюg можно утверждать, что х_ист. лежит в доверительном интервале, построенном около , а полуширина этого интервала при n < 30 рассчитывается по известной формуле:

Dх = t_{g, n} . (36)

Следовательно х_ист. = ± Dх, или - Dх < х_ист. < + Dх. (37)

В теории ошибок величину

S = (38)

называют средней квадратичной ошибкой прямо измеряемой величины х, величину Dх (см. (36)) – её абсолютной ошибкой, а величину e = × 100 % – относительной ошибкой, оценивающей точность измерений.

При косвенных измеренияхискомую величинуZ вычисляют по некоторой формуле

Z = f(x, y),

где x и y – прямо измеряемые величины.

Число значений x и y, полученных при измерении каждого из них, равно n:

x₁, х₂, х₃, …., х_n ;

у₁, у₂, у₃, … , у_n.

Теперь можно найти их средние арифметические значения:

= , = (39)

и средние квадратичные ошибки:

S_x = ; S_у = , (40)

Среднее арифметическое значение косвенно измеряемой величины вычисляют по формуле

= f( ). (41)

Истинное значение Z – Z_ист. лежит в доверительном интервале:

– DZ < Z_ист. < + DZ или Z_ист.= ± DZ. (3.7.5)

Полуширина данного интервала для нормально распределенной величины Z рассчитывается по формуле:

DZ = t_{g, n} . (43)

В (43) средняя квадратичная ошибка S_z косвенно измеряемой величины, равна:

= , (44)

где =Z_x´ и =Z_y´ – частные производные величины Z=f(x, y), соответственно, по x и по у, вычисляемые при их средних значениях, S_x и S_у – средние квадртичные ошибки величин х и у, значения которых получаются по формулам (40).

Окончательный результат обычно записывается в виде: Z_ист. = ± DZ, с указанием выбранного значения g. Приводится так же относительная ошибка косвенно измеряемой величины:

e = × 100 %.

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

h = h₀ ,

где h, r, t – вязкость, плотность и время истечения исследуемой жидкости из капилляра вискозиметра; h₀, r₀, t₀ – соответственно вязкость, плотность и время истечения эталонной жидкости (воды).

Величины h₀, r₀ и r считаем точно известными, t и t₀ измеряем секундомером, вязкость исследуемой жидкости – косвенно измеряемая величина.

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

для исследуемой жидкости t= 79, 2с;80,4с;78,0с; 83,6с; 80,2 с;

для воды t₀ = 51,0с; 48,4с; 50,6с; 47,4с; 44,2с.

2. Найдем по (39) средние арифметические значения t и t₀:

= = 80,28 с,

= = 48,32 с.

Определим по (41) среднее арифметическое значение вязкости исследуемой жидкости при: r = 790 , r₀ = 998,2 , h₀ = 1,0 × 10^-3 Па × с:

= h₀ ; = 1,0 × 10^-3 × = 1,31 × 10^-3 Па× с = 1,31 мПа× с.

3.Рассчитаем среднюю квадратичную ошибку вязкости по (44):

S_h = .

Для этого по (40) определим средние квадратичные ошибки времени истечения исследуемой жидкости S_t и воды :

S_t = =2,09 с

= = 2,75 с.

Найдем частные производные при t = и t₀ = ₀:

= = = 16,38 × 10^-6 Па ,

= - = – = -27,21 × 10^-6 Па.

Тогда S_h = = 82,2 × 10^-6 Па × с.

4. Определим полуширину доверительного интервала или абсолютную ошибку вязкости Dh по (43). Для этого, приняв доверительную вероятность g = 0,95, и, зная число измерений непосредственно определяемых величин (n = 5), найдем коэффициент Стьюдента, [cм. табл., напр. в (4, 9)], t_{g, n} = 2,78, тогда:

Dh = 2,78 × = 0,1× 10^-3 Па × с = 0,1 мПа × с.

Следовательно, с доверительной вероятностью g = 0,95 = 95% истинное значение вязкости исследуемой жидкости лежит в интервале

η = ± Dh = (1,31 ± 0,1) × 10^-3 Па × с = (1,31 0,1) мПа × с.

Относительная ошибка равна

e = × 100 % = » 7,6 %