Лекции по метрологии. Часть 4: Погрешности измерений

Вероятностный подход к описанию погрешностей

       

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

        Пусть величину А измеряли п раз и наблюдали при этом значе­ния а1, a2, а3,…,аi,...,аn. Случайная абсолютная погрешность еди­ничного измерения определяется разностью

                                                              Di = ai — A .                                             (4.1)

        Графически результаты от­дельных измерений представ­лены на рис. 4.2.

        При достаточно большом числе п одни и те же погрешности, если они имеют ряд дискретных значений, повторяются и поэтому можно установить относительную частоту (частость) их появления, т.е. отношение числа полу­ченных одинаковых данных mi к общему числу проведенных изме­рении п. При продолжении измерений величины А эта частота не изменится, поэтому ее можно считать вероятностью появления по­грешности при данных измерениях:  p(Ai) = mi / n.

clip_image014       

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

Интегральным законом, или функцией распределения вероятностей F(D) случайной погрешности Di в iопыте, называют функцию, значение которой для каждого D является вероятностью события Р(D), заключающегося в том, что случайная погрешность Di принимает значения, меньше некоторого значения D, т.е. функцию F(D) = Р[Di < D]. Эта функция при изменении D от -¥ до +¥ принимает значения от 0 до 1 и является неубывающей. Она существует для всех случайных  величин как дискретных, так и непрерывных (рис 4.3 а).

Если F(D) симметрична относительно точки А, соот­ветствующей вероятности 0,5 , то распределение результатов наблюдения будет симметрично относительно ис­тинного значения А. В этом случае целесообразно F(D) сдвинуть по оси абсцисс на значение DA, т.е. исключить систематическую составляющую погрешность (DА = Dс) и получить функцию распределения случайной составляющей погрешности D = clip_image016  (рис. 4.3 б). Функция рас­пределения вероятности погрешности D  отличается от функции распределения вероятности случайной состав­ляющей погрешности только сдвигом по оси абсцисс на значение систематической составляющей  погрешно­сти Dс .

Дифференциальным законом распределения вероятностей для случайной погрешности с непрерывной и дифференцируемой функцией распределения F(D) называют функцию clip_image018.

Эта зависимость есть плотность распределения вероятностей. График плотности распределения вероятностей может иметь различную форму в зависимости от закона распределения погрешностей. Для F(D), изображенной на рис. 4.3 б, кривая распределения f(D) имеет форму, близкую к форме колокола (рис. 4.3 в).

Вероятность появления случайных погрешностей определяется площадью, ограниченной кривой f(D) или её частью и осью абсцисс (рис. 4.3 в). В зависи мости от рассматриваемого интервала погрешности    clip_image020.


 

clip_image022
 

Значение f(D)dD есть элемент вероятности, равный площади прямоугольника с основанием dD и абсциссами D1 , D2, называемыми квантилями. Так как F(+¥)=1, то справедливо равенство            clip_image024,

т.е. площадь под кривой f(D) согласно правилу нормирования равна единице и отражает вероятность всех возможных событий.

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

        Математическое выражение нормального закона имеет вид

clip_image026,

где f(D) — плотность вероятности случайной погрешности D = аi — As — среднее квадратическое отклонение. Среднее квадратическое отклонение может быть выражено через случайные отклонения результатов наблюдений Di (см. формулу (4.1)):

clip_image028.

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

clip_image030 

опытных данных. Некоторые из наиболее часто встречающихся законов распределения приведены в ГОСТ 8.011-84 «Показатели точности измерений и формы представления результатов измерений».

    Основными характеристи- ками законов распределения являются математическое ожидание и дисперсия.

    Математическое ожидание случайной величины — это такое ее значение, вокруг которого группируются результаты отдельных наблюдений. Мате-матическое ожидание дискрет-ной случайной величины М[X] определяется как сумма произ-ведений всех воз­можных значений случайной величины на вероятность этих зна­чений                               clip_image032.

        Для непрерывных случайных величин приходится прибегать к интегрированию, для чего необходимо знать зависимость плот­ности вероятности от х , т. е.  f(х), где х=D.   Тогда     clip_image034.

        Это выражение означает, что математическое ожидание равно сумме бесконечно большого числа произведений всех возможных значений случайной величины х на бесконечно малые площади f(х)dх, где f(х) —ординаты для каждого х, a  dх — элементарные отрезки оси абсцисс.

        Если наблюдается нормальное распределение случайных по­грешностей, то математическое ожидание случайной погрешности равно нулю (рис. 4.4). Если же рассматривать нормальное распределение результатов, то математическое ожи­дание будет соответствовать истинному значению измеряемой ве­личины, которое мы обозначаем через A.

        Систематическая погрешность при этом представляет собой отклонение математического ожидания результатов наблюдений от истинного значения А измеряемой величины:  Dс = М[X] — A, а случайная погрешность — разность между результатом единичного наблюдения и математическим ожиданием:  clip_image036.

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

D[X] =Dx=M[(ai — mx)2].

Чем меньше дисперсия, тем меньше разброс отдельных результатов, тем точнее выполнены измерения. Однако дисперсия выражается в единицах в квадрате измеряемой величины. Поэтому в качестве характеристики точности ряда наблюдений наиболее часто применяют среднее квадратическое отклонение (СКО), равное корню квадратному из дисперсии :        clip_image038.

        Рассмотренное нормальное распределение случайных величин, в том числе и случайных погрешностей, является теоретическим, поэтому описанное нормальное распределение следует рассматри­вать как «идеальное», т. е. как теоретическую основу для изучения случайных погрешностей и их влияния на результат измере­ний.

        Далее излагаются способы применения этого распределения на практике с той или иной степенью приближения. Рассматривает­ся также еще одно распределение (распределение Стьюдента), применяемое при небольших количествах наблюдений.

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

Di = ai — A,

где  Di  -  погрешность  i-го измерения; ai -  результат i-го измерения.

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

                                   clip_image040.                              (4.2)

где аi  — результаты отдельных измерений; п — число измерений.

Теперь аналогично выражению (4.1) можно определить отклонение результата каж­дого измерения от среднего значения  clip_image042

                                                   clip_image044                                                      (4.3)

где vi  -  отклонение результата единичного измерения от среднего значения. Следует помнить, что сумма отклонений результата изме­рений от среднего значения равна нулю, а сумма их квадратов минимальна, т. е.

clip_image046       и          clip_image048 min.


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

Затем вычисляют оценку значения средней квадратической погрешности  для данного ряда измерений

                                clip_image050.                 (4.4)

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

                       clip_image052.            (4.5)

Ввиду того что среднее арифметическое значение clip_image042[1] также является случайной величиной, имеет смысл понятие среднеквадратического от­клонения среднего арифметического значения. Эту величину обозна­чим символом sср. Можно показать, что для независимых погрешностей

                                                 clip_image055.                                    (4.6)


Значение sср  характеризует степень разброса clip_image042[2]. Как указывалось вы­ше, clip_image042[3] выступает оценкой истинного значения измеряемой величины, т.е. является конечным результатом выполняемых измерений. Поэтому sср называют также средней квадратической погрешностью результата измерений.

На практике значением s, вычисляемым по формуле (4.5), пользуются в том случае, если необходимо дать характеристику точности применяемо­го метода измерения: если метод точен, то разброс результатов отдельных измерений мал, т.е. мало значение s. Значение же sср , вычисляемое по (4.6), используется для характеристики точности результата изме­рений некоторой величины, т.е. результата, полученного посредством математической обработки итогов целого ряда отдельных прямых из­мерений.

При оценке результатов измерений иногда пользуются понятием максимальной или пре­дельной допустимой погрешности, значение которой определяют в долях s или S . В настоящее время существуют разные критерии установления максимальной погрешности, т. е. границы поля допу­ска ±D, в которые случайные погрешности должны уложиться. Об­щепринятым пока является определение максимальной погрешно­сти  D = 3s (или 3S). В последнее время на основании ин­формационной теории измерений профессор П. В. Новицкий рекоменду­ет пользоваться значением D = 2s.

Введем теперь важные понятия доверительной вероятности и доверитель­ного интервала. Как указывалось выше, среднее арифметическое зна­чение clip_image042[4], полученное в результате некоторого ряда измерений, является оценкой истинного значения  А и, как правило, не совпадает с ним, а отличается на значение погрешности. Пусть Рд есть вероятность того, что clip_image042[5] отличается от А  не более чем на D, т.е.    Р(clip_image042[6]-D< А<clip_image042[7]+D)=Рд . Вероятность Рд называется доверительной вероятностью, а интервал значений измеряемой величины от clip_image042[8]-D  до clip_image042[9]+D — доверительным интервалом.

Приведенные выше неравенства означают, что с вероятностью Рд доверительный интервал от clip_image042[10]-D  до clip_image042[11]+D заключает в себе истинное значение А. Таким образом, чтобы характеризовать случайную погреш­ность достаточно полно, надо располагать двумя числами — доверитель­ной вероятностью и соответствующим ей доверительным интервалом. Если закон распределения вероятностей погрешностей известен, то по заданной доверительной вероятности можно определить доверительный интервал. В частности, при достаточно большом числе измерений часто бывает оправданным использование нормального закона, в то время как при небольшом числе измерений (п < 20), результаты которых принадлежат  нормальному распределению, следует пользоваться распределением Стьюдента. Это распределение имеет плотность вероят­ностей, практически совпадающую с нормальной при больших п, но зна­чительно отличающуюся от нормальной при малых п.

В табл. 4.1 приведены так называемые квантили распределения Стью­дента  ½t(n)½Рд  для числа измерений п = 2 — 20 и  доверительных вероятностей Р = 0,5 — 0,999.

Укажем, однако, что обычно таблицы распределения Стьюдента приводятся не для значений п и Рд , а для значений m = n-1 и  a =1 — Рд , что следует учитывать при пользовании ими. Чтобы определить довери­тельный интервал, надо для данных п и Рд найти квантиль ½t(n)½Рд и вычислить величины   Ан = clip_image042[12] - sср× ½t(n)½Рд  и Ав = clip_image042[13] + sср× ½t(n)½Рд, которые будут являться нижней и верхней  границами доверительного интервала.

После нахождения доверительных интервалов для заданной до­верительной вероятности согласно выше приведенной методике делают запись результата измерения в виде                        clip_image042[14]; D = Dн  ¸ Dв ; Рд ,

где clip_image042[15] - оценка истинного значения результата измерения в единицах измеряемой   величины;   D — погрешность   измерения;   Dв = +sср× ½t(n)½Рд    и Dн = -sср× ½t(n)½Рд — верхняя и нижняя границы погрешности измерения; Рд — доверительная вероятность [14].

 

Таблица 4.1  

 

n 

Значения квантилей распределения Стьюдента t(n) при доверительной

 вероятности Рд

0,5

0,6

0,7

0,8

0,9

0,95

0,98

0.99

0,995

0,999

2

1.000

1.376

1,963

3,080

6,310

12,71

31,80

63,70

127,3 

637,2

3

0.816

1,061

1,336

1,886

2,920

4,30

6,96

9,92

14,10

31,60

4

0.765

0.978

1,250

1,638

2,350

3,18

4,54

5,84

7,50

12,94

5

0.741

0.941

1,190

1,533

2,130

2,77

3,75

4,60

5,60

8,61

6

0,727

0.920

1,156

1,476

2,020

2,57

3,36

4,03

4,77

6,86

7

0,718

0.906

1,134

1,440

1,943

2,45

3,14

3,71

4,32

5,96

8

0.711

0.896

1,119

1,415

1,895

2,36

3,00

3,50

4,03

5,40

9

0.706

0.889

1,108

1,397

1,860

2,31

2,90

3,36

3,83

5,04

10

0,703

0.883

1,110

1,383

1,833

2,26

2,82

3,25

3,69

4,78

11

0.700

0.879

1,093

1,372

1,812

2,23

2,76

3,17

3,58

4,59

12

0.697

0.876

1,088

1.363

1,796

2,20

2,72

3,11

3,50

4,49

13

0.695

0.873

1,083

1,356

1,782

2,18

2,68

3,06

3,43

4,32

14

0.694

0.870

1,079

1,350

1,771

2,16

2,65

3,01

3,37

4,22

15

0.692

0,868

1,076

1,345

1,761

2,14

2,62

2,98

3,33

4,14

16

0.691

0,866

1,074

1,341

1,753

2,13

2,60

2,95

3,29

4,07

17

0.690

0.865

1,071

1,337

1,746

2,12

2,58

2,02

3,25

4,02

18

0,689

0.863

1,069

1,333

1,740

2,11

2,57

2,90

3,22

3,96

19

0.688

0,862

1,067

1,330

1,734

2,10

2,55

2,88

3,20

3,92

20

0.688

0.861

1,066

1,328

1,729

2,09

2,54

2,86

3,17

3,88

¥

0,674

0,842

1,036

1,282

1,645

1,96

2,33

2,58

2,81

3,29

       

        Оценка погрешностей результатов косвенных измерений. При косвенных измерениях искомая величина А функционально связана с одной или несколькими непосредственно измеряемыми ве­личинами: х, y,...,t. Рассмотрим простейший случай определения по­грешности при одной переменной, когда A=F(x). Обозначив абсо­лютную погрешность измерения величины х через  ±Dx , получим  A+DA = F(x±Dx).

   Разложив правую часть этого равенства в ряд Тейлора и прене­брегая членами разложения, содержащими Dх в степени выше пер­вой, получим

A+DA » F(x) ± clip_image057Dx             или             DA »  ± clip_image057[1]Dx.

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

clip_image059.


Если измеряемая величина А является функцией нескольких пе­ременных: A=F(x, y,...,t), то абсолютная погрешность результата косвенных измерений

                                  clip_image061  .


Частные относительные погрешности косвенного измерения опреде­ляются по формулам  clip_image063;   clip_image065 и т. д. Относительная погрешность результата измерений

clip_image067.


        Остановимся также на особенностях оценки результата косвенного измерения при наличии случайной погрешности.

Для оценки случайной погрешности результатов косвенных измерений величины А бу­дем полагать, что систематические погрешности измерений вели­чин x, y,…, t  исключены, а случайные погрешности измерения этих же величин не зависят друг от друга.

При косвенных измерениях значение измеряемой величины на­ходят по формуле                                          clip_image069,

где clip_image071 - средние или средние взвешенные значения величин  x, y,…, t .

Для вычисления среднего квадратического отклонения значе­ния clip_image042[16] измеряемой величины А целесообразно использовать сред­ние квадратические отклонения, полученные при измерениях x, y,…, t .

В общем виде для определения среднего квадратического от­клонения s косвенного измерения служит следующая формула:

                                                     clip_image073,                                   (4.7)


где Dx ; Dy ;…; Dt так называемые частные погрешности косвенного измерения  clip_image075; clip_image077; …; clip_image079;   clip_image081; clip_image083; … ; clip_image085 частные производные А по x, y,…, t ;   sx; sy ,…, st , …— средние квадратические отклонения результатов измере­ний величин x, y,…, t .

Рассмотрим некоторые частные случаи применения уравнения (4.7), когда функциональная зависимость между косвенно и непосредственно измеряемыми величинами выражается формулой A = k×xa×yb×zg , clip_image087где k — числовой коэффициент (безразмерный).

В этом случае формула (4.7) примет следующий вид:

clip_image089.


Если a = b = g = 1 и  A = k×x×y×z, то формула относительной по­грешности упрощается до вида                clip_image091.


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

 

4.5. Правила суммирования случайных и систематических погрешностей

      

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

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

        Арифметическое суммирование  clip_image093или clip_image095дает максимальную погрешность прибора, которая имеет ничтожно малую вероятность и поэтому редко используется для оценки точности ра­боты прибора в целом. Согласно теории ошибок результирующая погрешность sрез и Мрез определяется сложением по квадратическому закону   clip_image097 или  clip_image099.

Аналогично определяется и результирующая относительная по­грешность измерения:                                     clip_image101.                                              (4.8)

        Уравнение (4.8) можно использовать для определения допу­стимых погрешностей отдельных блоков разрабатываемых прибо­ров с заданной общей погрешностью измерения. При конструирова­нии прибора обычно задаются равными погрешностями для отдель­ных входящих в него блоков. Если существует несколько источни­ков погрешностей, которые на конечный результат измерения влияют неодинаково (или прибор состоит из нескольких блоков с разными погрешностями), в формулу (4.8) следует ввести весовые коэффициенты ki :          

                                          clip_image103,                  (4.9)             

где d1, d2, … , dm — относительные погрешности отдельных узлов (бло­ков) измерительного прибора; k1, k2, … , km  - коэффициенты, учиты­вающие степень влияния случайной погрешности данного блока на результат измерения.

        При наличии у измерительного прибора (или его блоков) также и систематических погрешностей общая погрешность определяется их суммой:

clip_image105,  где  dс(im) — система­тическая погрешность от воздейст- вия на i-й блок m-го фактора; di - случайные погрешности для i-го блока.

        Суммирование погрешностей, имеющих взаимную корреляцион­ную связь, основано на следующем положении теории вероятностей: дисперсия суммы двух коррелированных случайных величин, харак­теризующихся дисперсиями clip_image107 и  clip_image109 и коэффициентом корреля­ции r12, определяется выражением

      clip_image111.

Из  этого следует, что средняя квадратическая результирующая по­грешность вычисляется по формуле             

                                                clip_image113 .                                 (4.10)

        На практике обычно пользуются двумя крайними случаями формулы (4.10); при r12 » ±1, когда составляющие погрешности суммируют­ся алгебраически: clip_image113[1]= s1+ s2  и при   r12 » 0 , когда погрешности суммируются геометрически: clip_image116. Такой же подход справедлив и для большего числа составляю­щих.

При оценке влияния частных погрешностей следует учитывать, что точность измерений в основном зависит от погрешностей, боль­ших по абсолютной величине, а некоторые наименьшие погрешно­сти можно вообще не учитывать. Частная погрешность оценивается на основании так называемого критерия ничтожной погрешности, который заключается в следующем. Допустим, что суммарная по­грешность dрез определена по формуле (4.8) с учетом всех m частных погрешностей, среди которых некоторая погрешность di имеет ма­лое значение. Если суммарная погрешность d¢рез, вычисленная без учета   погрешности   di,   отличается   от   dрез   не   более   чем   на   5 %, т.е. dрез-d¢рез< 0,05×dрез или 0,95×dрез<d¢рез , то di можно считать ничтожной погрешностью. Принимая во внимание, что (d¢рез)2 = d2рез — d2i , легко установить критерий ничтожной погрешности: di £ 0,3×dрез . Это означает, что если частная погрешность меньше 30 % общей погрешности, то этой частной погрешностью можно пренебречь. В случае нескольких погрешностей критерий ничтожности их сово­купности имеет вид    clip_image118 .

В практике технических расчетов часто пользуются менее стро­гим критерием — в эти формулы вводят коэффициент 0,4.

 

Вы здесь: Главная Метрология Лекции по метрологии. Часть 4: Погрешности измерений