Начальные сведения из теории ошибок

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

По своей природе ошибки бывают грубые, систематические и случайные.

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

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

Случайные ошибки измерений обусловлены точностью способа измерений (строгостью теории), точностью измерительного прибора, квалификацией исполнителя и влиянием внешних условий. Закономерности случайных ошибок проявляются в массе, то-есть, при большом количестве измерений; такие закономерности называют статистическими. Освободить результат единичного измерения от случайных ошибок невозможно; невозможно также предсказать случайную ошибку единичного измерения. Теория ошибок занимается в основном изучением случайных ошибок.

Случайная истинная ошибка измерения Δ - это разность между измеренным значением величины l и ее истинным значением X:

(1.25)

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

при данных условиях измерений абсолютные значения случайных ошибок не превосходят некоторого предела; если какая-либо ошибка выходит за этот предел, она считается грубой,

положительные и отрицательные случайные ошибки равновозможны,

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

(1.26)

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

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

Средняя квадратическая ошибка одного измерения. Для оценки точности измерений можно применять разные критерии; в геодезии таким критерием является средняя квадратическая ошибка. Это понятие было введено Гауссом; он же разработал основные положения теории ошибок. Средняя квадратическая ошибка одного измерения обозначается буквой m и вычисляется по формуле Гаусса:

(1.27)

где: ;
n - количество измерений одной величины.

Средняя квадратическая ошибка очень чувствительна к большим по абсолютной величине ошибкам, так как каждая ошибка возводится в квадрат. В то же время она является устойчивым критерием для оценки точности даже при небольшом количество измерений; начиная с некоторого n дальнейшее увеличение числа измерений почти не изменяет значения m; доказано, что уже при n = 8 значение m получается достаточно надежным.

Предельная ошибка ряда измерений обозначается Δпред; она обычно принимается равной 3*m при теоретических исследованиях и 2*m или 2.5*m при практических измерениях. Считается, что из тысячи измерений только три ошибки могут достигать или немного превосходить значение Δпред = 3*m.

Отношение mx/X называется средней квадратической относительной ошибкой; для некоторых видов измерений относительная ошибка более наглядна, чем m. Относительная ошибка выражается дробью с числителем, равным 1, например, mx/X = 1/10 000.

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

F = f( X, Y, Z ... ), (1.28)

здесь: X, Y, Z ... - истинные значения аргументов,
F - истинное значение функции.

В результате измерений получены измеренные значения аргументов lX, lY, lZ, при этом:

(1.29)

где ΔX, ΔY, ΔZ - случайные истинные ошибки измерения аргументов.

Функцию F можно выразить через измеренные значения аргуметов и их истинные ошибки:

Разложим функцию F в ряд Тейлора, ограничившись первой степенью малых приращений ΔX, ΔY, ΔZ:

(1.30)

Разность является случайной истинной ошибкой функции с противоположным знаком, поэтому:

(1.31)

Если выполнить n измерений аргументов X, Y, Z, то можно записать n уравнений вида (1.31). Возведем все эти уравнения в квадрат и сложим их; суммарное уравнение разделим на n и получим

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

(1.32)

Как частные случаи формулы (1.32) можно написать выражения для средней квадратической ошибки некоторых функций:

Если функция имеет вид произведения нескольких аргументов,

F = x * y * z,

то для нее можно записать выражение относительной ошибки функции:

(1.33)

которое в некоторых случаях оказывается более удобным, чем формула (1.32).

Принцип равных влияний. В геодезии часто приходится определять средние квадратические ошибки аргументов по заданной средней квадратической ошибке функции. Если аргумент всего один, то решение задачи не представляет трудности. Если число аргументов t больше одного, то возникает задача нахождения t неизвестных из одного уравнения, которую можно решить, применяя принцип равных влияний. Согласно этому принципу все слагаемые правой части формулы (1.32) или (1.33) считаются равными между собой.

Арифметическая середина. Пусть имеется n измерений одной величины X, то-есть,

(1.34)

Сложим эти равенства, суммарное уравнение разделим на n и получим:

(1.35)

Величина (1.36)

называется средним арифметическим или простой арифметической серединой. Запишем (1.35) в виде

по третьему свойству ошибок (1.26) можно написать:

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

Запишем формулу (1.36) в виде

и подсчитаем среднюю квадратическую ошибку арифметической середины, которая обозначается буквой M. Согласно формуле (1.32) напишем:

или

Но ml1 = ml2 = ... = mln= m по условию задачи, так как величина X измеряется при одних и тех же условиях. Тогда в квадратных скобках будет n * m2, одно n сократится и в итоге получим:

M2 = m2/n

или

(1.37)

то-есть, средняя квадратическая ошибка арифметической середины в корень из n раз меньше ошибки одного измерения.

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

Пусть имеется ряд равноточных измерений величины X:

l1, l2 , ..., ln .

Вычислим арифметическую середину X0 = [1]/n и образуем разности:

(1.38)

Сложим все разности и получим [l] - n * X0 = [V]. По определению арифметической середины n * X0 = [l], поэтому:

[V] = 0. (1.39)

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

(1.40)

Приведем вывод этой формулы. Образуем разности случайных истинных ошибок измерений Δ и вероятнейших ошибок V:

(1.41)

Разность (X0 - X) равна истинной ошибке арифметической середины; обозначим ее Δ0 и перепишем уравнения (1.41):

(1.42)

Возведем все уравнения (1.42) в квадрат, сложим их и получим:

.

Второе слагаемое в правой части этого выражения равно нулю по свойству (1.39), следовательно,

.

Разделим это уравнение на n и учтя, что [Δ2]/n =m2, получим:

(1.43)

Заменим истинную ошибку арифметической середины Δ0 ее средней квадратической ошибкой ; такая замена практически не изменит правой части формулы (1.43). Итак,

,
откуда ;

после перенесения (n-1) в правую часть и извлечения квадратного корня получается формула Бесселя (1.40).

Для вычисления средней квадратической ошибки арифметической середины на основании (1.37) получается формула:

(1.44)

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

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

p = C/m2 (1.45)

где C - в общем случае произвольное положительное число.

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

(1.46)
или X0 = [l*p] / [p] .

Ошибку измерения, вес которого равен 1, называют средней квадратической ошибкой единицы веса; она обозначается буквой m. Из формулы (1.45) получаем

откуда (1.47)

то-есть, за число C принимают квадрат ошибки единицы веса.

Подсчитаем вес P средневесовой арифметической середины. По определению веса имеем:

(1.48)

Согласно (1.46) и (1.32) напишем:

Подставим сюда вместо mli2 их выражения через вес m2 = C/p , тогда:

Подставим это выражение в формулу (1.48) и получим,

P = [p], (1.49)

то-есть, вес средневесовой арифметической середины равен сумме весов отдельных измерений.

В случае равноточных измерений, когда веса всех измерений одинаковы и равны единице, формула (1.49) принимает вид:

P = n. (1.50)

При обработке больших групп измерений (при уравнивании геодезических построений по МНК) вычисляются значение ошибки единицы веса, веса измерений и других элементов после уравнивания, а ошибка любого уравненного элемента подсчитывается по формуле:

(1.51)

где pi - вес i-того элемента.