Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
geokniga-геодезия-попов-вн-чекалин-ви-2007.pdf
Скачиваний:
3640
Добавлен:
16.03.2015
Размер:
39.37 Mб
Скачать

Если результат измерения заранее известен, то, казалось бы, зачем производить измерения? Однако такие действия часто приходится выполнять. Например, при проверке правильности работы или показаний прибора по эталону. Да и при самих непосредственных измерениях, например, углов в треугольнике, сумма углов треугольника (или многоугольника) является эталоном, известной величиной.

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

Рассмотрим ряд измерений одной и той же величины Х для случая, когда число измерений весьма большое (n → ∞). Составим ряд истинных погрешностей измерений, полагая, что измеряемая величина нам известна.

Сложим все разности в правых и левых частях формул (3.3) и разделим полученные результаты на n, получим

[

]

=

[x]

X

(3.6)

n

 

 

n

 

 

 

 

 

 

В соответствии со свойствами случайных погрешностей отношение [ ]/n стремится к нулю при n → ∞ . Отношение [х]/n = хо называется средним

арифметическим из результатов измерений.

 

С учетом сказанного можно записать, что

 

(хо Х) n ,

(3.7)

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

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

ями результатов измерений от среднего арифметического:

 

vi = хi - xо

(3.8)

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

§ 23. Средняя квадратическая погрешность

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

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

70

ками ординат при значениях + и − . Для нормального закона распределения вероятность появления погрешностей в установленных диапазонах равна следующим значениям :

для диапазона ± → Р = 68,3% (≈ 68%); для диапазона ±2Δ→ Р = 95,5% (≈ 95%);

для диапазона ±3Δ→ Р = 99,7% (практически 100%).

 

Таким образом, только в 3-х слу-

 

чаях из 1000 может появиться по-

 

грешность, превышающая значение

 

3 . Такие погрешности принято счи-

 

тать грубыми, и результаты измере-

 

ний, содержащие такие погрешности,

 

исключают из дальнейшей обработ-

 

ки. В некоторых случаях, для ужесто-

 

чения требований к точности измере-

Рис. 3.1. Нормальный закон распределения

ний, устанавливают предельную по-

случайных погрешностей

грешность до 2 (или до 2m).

 

Часто значение СКП указывают с

коэффициентом t (коэффициент Стьюдента), который и определяет доверительный вероятностный интервал (х ± tm) результата измерений при установленном уровне вероятности Р. Для этого удобно пользоваться табл. 3.1.

Например, необходимо определить доверительный интервал для величины Х с вероятностью 75%. По таблице интерполированием находим, что для

Р1 = 72,9 t1 = 1,1, для Р2 = 77,0 t2 = 1,2: tх 1,15.

Это значит, что результат измерений с вероятностью 75% находится в

пределах

(Х ± 1,15 m).

 

 

 

Таблица 3.1

 

 

 

 

 

t

P%

t

P%

t

P%

0,1

8,0

1,1

72,9

2,1

96,4

0,2

15,9

1,2

77,0

2,2

97,2

0,3

23,6

1,3

80,6

2,3

97,9

0,4

31,1

1,4

83,8

2,4

98,4

0,5

38,3

1,5

86,6

2,5

98,8

0,6

45,1

1,6

89,0

2,6

99,1

0,7

51,6

1,7

91,1

2,7

99,3

0,8

57,6

1,8

92,8

2,8

99,5

0,9

63,2

1,9

94,3

2,9

99,6

1,0

68,3

2,0

95,5

3,0

99,7

Если измеряемая величина Х известна, то значение СКП определяется по

формуле Гаусса:

 

 

 

 

m =

[ 2 ]

,

(3.9)

 

 

n

 

 

где

- истинные погрешности измерений.

 

71

Для случаев, когда измеряемая величина неизвестна, используется

формула Бесселя:

m =

 

2 ]

 

,

(3.10)

(n − 1)

 

 

 

 

где v - уклонения результатов измерений от среднего арифметического.

Как видно из формул (3.9) и (3.10), в случае, когда измеряемая величина известна, для оценки точности достаточно уже одного измерения (оно и является необходимым). Как уже указывалось выше, чаще всего формулу Гаусса используют при оценках точности эталонируемых приборов при измерении известных величин (эталонов). Для оценки точности по формуле Бесселя необходимыми являются как минимум два измерения. Формула Бесселя используется при оценках точности результатов массовых (многократных) измерений одной величины.

При возрастании числа измерений значения СКП, полученные по формулам Гаусса и Бесселя, становятся практически одинаковыми (примерно с n ³ 20). При этом значение СКП одного измерения стремится к пределу mпред, который определяется точностью прибора, точностью метода или программы измерений. Очевидно, что на практике невозможно, да и нецелесообразно по ряду причин, обеспечивать весьма большое число измерений одной величины. При этом практическое число измерений должно обеспечивать получение результата измерения с заданной точностью при установленном уровне доверительной вероятности.

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

mm =

 

m

 

.

(3.11)

 

 

 

 

 

 

2n

 

§ 24. Средние квадратические погрешности функций измеренных величин

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

Предположим, что имеется функция F аргументов х1 , х2 , ... , хn:

F = f (x1, х2, …, хn). (3.12) Величины хi известны из непосредственных измерений, а также известны

и их СКП: m1 , m2 , ... , mn . В этом случае СКП функции определяется по следующей формуле:

72

 

é

æ

f

ö 2

 

2

 

æ

f

ö

2

 

 

ê

ç

 

 

 

ç

 

 

2

 

mF =

x

m1

 

+

x

 

m2

+ ... +

ê

ç

 

 

ç

 

 

 

ë

è

1

ø

 

 

 

è

 

2 ø

 

 

æ

f

ö 2

 

2

ù

(3.13)

ç

 

 

mn

ú ,

ç

x

 

 

 

 

ú

 

è

 

n ø

 

 

û

 

где (f/хi) - частная производная функции по аргументу хi . Правила определения СКП функций следующие.

1.Выполнить последовательно дифференцирование функции отдельно по каждому из аргументов, считая остальные аргументы постоянными числами (коэффициентами).

2.Полученные выражения умножить на СКП аргументов, по которым производилось дифференцирование функции и возвести полученные выражения каждое отдельно в квадрат.

3.Записать полученные выражения в виде суммы под знаком квадратного корня.

Рассмотрим несколько примеров определения СКП функций.

Пример 3.1. Средняя квадратическая погрешность среднего арифметического. Очевидно, что значение среднего арифметического является функцией суммы

измеренных величин хi (3.1). Представим это выражение в виде

хо = (х1 + х2 + … + хn ) / n (3.14) Поскольку 1/n является постоянным коэффициентом, то при почленном

дифференцировании и после умножения на mi и возведения в квадрат пролучим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

m

x

= M =

 

[(m / n)

2 + (m

2

/ n)2

+ ...+ (m

n

/ n)2

]

(3.15)

 

o

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mxo = M =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.16)

 

 

(m1

2

+ m2

2

 

+ ... +

mn

2 )/ n

 

 

 

 

Полагая измерения равноточными, т.е. m1

 

= m2

= ... = mn

=

m, выражение (3.16)

преобразуем к виду

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M =

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т.е. СКП среднего арифметического в корень из числа измерений меньше СКП одного измерения.

С учетом (3.10)

M

=

 

2 ]

 

(3.18)

 

n(n - 1)

 

 

 

 

 

Очевидно, что, если при увеличении числа измерений значение СКП одного измерения стремится к предельному значению, отличному от нуля, то значение СКП среднего арифметического стремится при увеличении числа измерений к нулю, а само среднее арифметическое – к истинному значению.

Пример 3.2. Объем пирамиды,

основанием которой является прямоугольник,

определен по формуле

hab

 

V =

(3.19)

 

3

 

где h – высота пирамиды, а и b – стороны основания.

Требуется определить СКП объема пирамиды, вычисленного по формуле (3.19), если известно, что h = 12,34 м, а = 23,46 м, b = 39,63 м и их СКП равны соответственно: mh = 0,07 м, ma = 0,02 м, mb = =0,04 м.

Решение.

Выполняем последовательное дифференцирование по аргументам h, a и b:

73

-

по аргументу h:

V

=

ab mh ;

h

 

 

 

3

-

по аргументу а:

V

=

hb ma ;

a

 

 

 

3

-по аргументу b: Vb = ah3 mb .

Возводим в квадрат полученные части и записываем в виде суммы в подкоренном выражении:

 

 

 

 

 

 

 

 

 

mV =

1

[(ab)2 mh

2

+ (hb)2 ma

2 + (ha)2 mb

2 ]

(3.20)

 

9

 

 

 

 

 

 

 

Формулу (3.20) можно преобразовать к следующему виду. Разделим правую и левую части соответственно на ( h a b) и V, получим

 

m

=

1

æ

m

a

ö 2

æ

m

ö 2

æ

m

ö 2

 

 

(3.21)

 

V

3

ç

 

 

 

+ ç

 

 

b

+ ç

 

h

 

V

 

è

a

ø

 

è

 

b ø

è

h ø

 

или

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ V =

 

δ a

2

+ δ b

2

+ δ h

2

 

 

 

 

 

 

 

(3.22)

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где δ - относительные СКП аргументов и функции.

 

 

 

 

 

 

Выполним вычисления по формуле (3.20).

 

 

 

 

 

 

 

 

 

 

 

mV = 1

 

 

 

 

м3.

(23,46 ×39,63 ×0,70)2

+ (12,34 ×23,46 ×0,04)2

+ (12,34 ×39,63 ×0,02)2 = 7,42

3

V = (23,46 · 39,63 х12,34) : 3 = 1274,75 м3.

 

 

 

 

 

Значение

 

 

 

 

 

1

Относительная СКП определения объема равна δV = 7,42 /1274,75 = 0,00582 = 1718 .

Выполним проверку значения δV по формуле (3.22). Относительные СКП аргументов равны:

δh = 0,07 : 12,34 = 0,00567, δa = 0,02 : 23,46 = 0,00085, δb = 0,04 : 39,63 = 0,00101.

После подстановки в формулу (3.22) получим δV = 0,00582, что совпадает с предыдущим результатом.

Пример 3.3. Сторона а треугольника определена по теореме синусов по значению стороны b и двум углам треугольника А и В:

Известно: b = 140,12 м (δb = 1 : 2000; mb = 140,12 : 2000 = 0,07 м), А = 73о18,8' (mА = 0,4'), В = 63о05,6' (mВ = 0,3').

Необходимо определить СКП стороны а.

Решение.

Вычисление стороны а производится по формуле

 

a = b sin A

(3.23)

 

sin B

 

Запишем члены подкоренного выражения для СКП параметра а:

 

-

для аргумента b: (mb sin A / sin B )2;

 

-

для аргумента А: (b mA cos A ' sin B )2;

 

где ρ' = 3438′ (число минут в радиане; для выражения угловой меры СКП угла в меру радианную);

- для аргумента В: (b mB sin A cos B / ρ' sin2 B)2. Следовательно,

m

 

=

æ

mb sin A ö

2

æ

mAb cos A ö 2

æ

mBbsin Acos B ö

2

 

a

ç

 

 

+ ç

 

 

+ ç

 

2

 

 

(3.24)

 

 

è

sin B

ç

ρ ¢ sin B

 

ç

ρ ¢ sin

B

 

 

 

 

 

 

ø

 

è

ø

è

 

ø

 

 

После подстановки значений аргументов получим: mа= 0,096 м.

74