Совокупность операций по применению технического средства, хранящего единицу физической величины, заключающихся в сравнении измеряемой величины, с её единицей, с целью полученпия значения этой величины или информации о нём в формате, наиболее удобном для использования называется измерением.
Измерения могут быть классифицированы:
равноточные и неравноточные;
однократные
статические
технические
абсолютные
− по общим приёмам получения результатов измерений на прямые , косвенные, совместные и совокупные.
Оценка результата измерения
Задача оценки результата измерения состоит в том, чтобы по полученным экспериментальным путём результатам наблюдений, содержащим случайные погрешности, найти оценку истинного значения измеряемой величины – результат измерения. Будем полагать, что систематические погрешности в результатах наблюдений отсутствуют или исключены.
состоятельной
несмещённой
Способы нахождения оценок результата зависят от вида функции распределения и от имеющихся соглашений по этому вопросу, регламентируемых в законодательной метрологии. Общие соображения по выбору оценок заключаются в следующем.
Распределения погрешностей результатов наблюдений, как правило, являются симметричными относительно центра распределения, поэтому истинное значение измеряемой величины может быть определено как координата центра рассеивания хц, т.е. центра симметрии распределения случайной погрешности (при условии, что систематическая погрешность исключена).
Отсюда следует принятые в метрологии правило оценивания случайной погрешности в виде интервала, симметричного относительно результата измерения ( хц±∆х).
Координата хц может быть найдена несколькими способами. Наиболее общим является определение центра в симметрии из принципа симметрии вероятности, т.е. нахождение такой точки на оси х, слева и справа от которой вероятности появления различных значений случайных погрешностей равны между собой и составляют Р1=Р2=0,5.Такое значение хц называется медианой.
Методы и средства измерений и виды погрешностей
... погрешности измерений. Оценка достоверности результата измерений, т.е. определение погрешности измерений - одна из основных задач метрологии [1]. Погрешность -- это отклонение результата измерения от истинного значения измеряемой величины. Погрешности условно можно разделить на погрешности средств измерения и погрешности результата. измерений. Погрешности средств измерения, Погрешность результата ...
Координата хц может быть определена и как центр тяжести распределения, т.е. как математическое ожидание случайной величины.
При ассиметричной кривой плотности распределения вероятностей оценкой центра распределения может служить абсцисса моды распределения, т.е. координата максимума плотности. Однако есть распределения, у которых не существует моды (например, равномерное), и распределения, у которых не существует математического ожидания.
В практике измерений встречаются различные формы кривой закона распределения, однако чаще всего имеют дело с нормальным и равномерным распределением плотности вероятностей.
Учитывая многовариантность подходов к выбору оценок в целях обеспечения единства измерений, правила обработки результатов наблюдений обычно регламентируются нормативно-техническими документами (стандартами, методическими указаниями, инструкциями).
Варианты оценки случайных погрешностей
Для количественной оценки случайных погрешностей и установления границ случайной погрешности результата измерения могут использоваться: предельная погрешность, интервальная оценка, числовые характеристики закона распределения. Выбор конкретной оценки определяется необходимой полнотой сведений о погрешности, назначением измерений и характером использования их результатов. Комплексы оценок показателей точности установлены стандартами.
Предельная погрешность ∆m — погрешность, больше которой в данном измерительном эксперименте не может появиться. Теоретически, такая оценка погрешности правомерна только для распределений, границы которых чётко выражены и существует такое значение ±∆m, которое ограничивает возможные значения случайных погрешностей с обеих сторон от центра распределения ( например, равномерное).
На практике такая оценка есть указание наибольшей погрешности, которая может встретиться при многократных измерениях одной и той же величины.
Недостатком такой оценки является то, что она не содержит информации о характере закона распределения случайных погрешностей. При арифметическом суммировании предельных погрешностей получаемая сумма может значительно превышать действительные погрешности.
Более универсальным и информативным являются квантильные оценки. Площадь, заключённая под всей кривой плотности распределения погрешностей, отражает вероятность всех возможных значений погрешности и по условиям нормирования равна единице. Эту
площадь можно разделить вертикальными линиями на части. Абсциссы таких линий называются квантилями. Так, на рис. 1 ∆Х1 есть 25%-ная квантиль, так как площадь под кривой f(∆x) слева составляет 25% всей площади. Абсцисса ∆Х2 соответствует 75%-ной квантили. Между ∆Х1 и ∆Х2 заключено 50% всех возможных значений погрешности, а остальные лежат вне этого интервала.
Квантильная оценка погрешности представляется интервалом от -∆Х(Р) до +∆Х(Р), на котором с заданной вероятностью Р встречаются Р∙100% всех возможных значений случайной погрешности. Интервал с границами ±∆Х(Р) называется доверительным интервалом случайной погрешности, а соответствующая ему вероятность – доверительной вероятностью. Принято границы доверительного интервала указывать симметричными относительно результата измерения.
Погрешности и их классификации
... числового выражения делятся на три группы а) Абсолютные погрешности - ДX является оценкой абсолютной ошибки измерения. Величина этой погрешности зависит от способа её вычисления, который, всвою очередь, определяется распределением случайной величины Xmeas. При ...
Так как квантили, ограничивающие доверительный интервал погрешности могут быть выбраны различными, то при оценивании случайной погрешности доверительными границами необходимо одновременно указывать значение принятой доверительной вероятности ( например, ±0,3 при Р=0,95).
Доверительные границы случайной погрешности ∆Х(Р), соответствующие доверительной вероятности Р, находят по формуле ∆Х(Р)=tσ, где t – коэффициент, зависящий от Р и от формы закона распределения.
В целях единообразия в оценивании случайных погрешностей интервальными оценками при технических измерениях доверительная вероятность принимается равной 0,95. Лишь для особо точных и ответственных измерений допускается применение более высокой доверительной вероятности.
Недостатком оценивания случайной погрешности доверительным интервалом при произвольно выбираемых доверительных вероятностях является невозможность суммирования нескольких погрешностей, доверительный интервал суммы не равен сумме доверительных интервалов. В то же время необходимость в суммировании случайных погрешностей существует, когда нужно оценить погрешность суммированием её составляющих, подчиняющихся к тому же разным законам распределения.
В теории вероятностей показано, что суммирование статистически независимых случайных величин осуществляется путём суммирования их дисперсий
= или (1)
Таким образом, для того чтобы отдельные составляющие случайной погрешности можно было суммировать расчётным путём, они должны быть определены своим СКО, а не предельными или доверительными границами.
Формула (1) правомерна только для некорелированных случайных величин. В том случае, когда суммируемые составляющие погрешности коррелированы, расчётные соотношения усложняются, так как требуется учёт корреляционных связей. Методы выявления корреляционных связей и их учёт являются предметом изучения в теории вероятностей.
Рассотренные свойства распределений следует понимать как «идеальные», полученные на основе бесконечно большого числа опытов. В реальных условиях результат измерения получают либо путём обработки ограниченной группы наблюдений, либо на основе однократного измерения. Правила обработки данных для получения оценок результата и погрешности статистических измерений определены стандартами Государственной системы обеспечения единства измерений.
Измерения с однократными наблюдениями.
Прямые однократные измерения с точным оцениванием погрешности
Большинство измерений являются однократными. В обычных условиях их точность вполне приемлема. Результат однократного измерения Qi, записывается следующим образом
где Xi значение i-ro показания; θi — поправка.
При однократных измерениях для получения результата измерения используется одно значение отсчёта показаний прибора. Будучи по сути дела случайным, однократный отсчёт Х включает в себя инструментальную, методическую и личную составляющие погрешности измерения, в каждой из которых может быть выделены систематическая и случайная составляющая.
Сравнительно легко, путём проверки или по паспортным данным может быть получена оценка систематической погрешности прибора, а анализом метода измерения – оценка систематической погрешности методического происхождения. При наличии в документации на прибор сведений о дополнительных систематических погрешностях, обусловленных влияющими величинами, эти погрешности также оцениваются и учитываются. После исключения из отсчета всех известных систематических погрешностей можно полагать, что погрешность исправленного результата хиспр состоит из неисключённых остатков систематических погрешностей и случайных составляющих погрешностей.
Необходимым условием проведения однократного измерения служит наличие априорной информации. К ней относится, например, информация о виде закона распределения вероятности показания и мере его рассеивания, полученная из опыта предшествующих измерений. Такой информацией может служить класс точности средства измерения. Без априорной информации выполнение однократного измерения бессмыслено.
За результат измерения в этом случае принимают результат однократного наблюдения х (с введением поправки, если она имеется), используя предварительно полученные (например, при разработке методики выполнения измерений) данные об источниках, составляющих погрешность.
Доверительные границы неисключенного остатка систематической погрешности результата измерения 0(Р) вычисляют по формуле:
, (2)
т
Таблица 1
m |
≥5 |
4 |
3 |
2 |
k(P) |
1.45 |
1.40 |
1.30 |
1.20 |
Если составляющие неисключепного остатка систематической погрешности распределены равномерно и заданы доверительными границами θi ( Pi ), то доверительную границу результата измерения вычисляют по формуле:
(2)
где к и к i — те же, что и в предыдущем случае, коэффициенты, соответствующие доверительной вероятности Р и Р i соответственно; т — число составляющих неисключенного остатка систематической погрешности.
Среднее квадратическое отклонение результата измерения с однократным наблюдением вычисляют одним из следующих способов:
- Если в технической документации на средство измерения или в методике выполнения измерения, указаны нормально распределенные составляющие случайной погрешности результата наблюдения (инструментальная, методическая, из-за влияющих факторов, оператора и т.д.), то среднее квадратическое отклонение определяют по формуле:
(3)
где п — число составляющих случайной погрешности Si -значения среднего квадратического отклонения этих составляющих.
Доверительную границу случайной погрешности результата измерения ε(Р) в этом случае вычисляют по формуле:
ε(Р)= Z p /2 · S ( x ) (4)
где zp /2 — значение нормированной функции Лапласа в точкеР/2 при доверительной вероятности Р (табл. 2):
Таблица 2
Р |
0,90 |
0,95 |
0,96 |
0,97 |
0,98 |
0,99 |
Z p /2 |
1,65 |
1,96 |
2,06 |
2,17 |
2,33 |
2,58 |
2. Если в тех же документах случайные составляющие погрешности результата наблюдения представлены доверительными границами ε(Р) при одной и той же доверительной вероятности Р, то доверительную границу случайной погрешности результата измерения с однократным наблюдением при доверительной вероятности Р вычисляют по формуле:
(5)
3. Если случайные составляющие погрешности результата наблюдения определяют предварительно в реальных рабочих условиях экспериментальными методами при числе наблюдений «; <30, то:
(6)
где t-коэффициент Стыодента, соответствующий наименьшему числу наблюдений п из всех nj ; Si (x)— оценки средних квадратических отклонений случайных составляющих погрешности результата наблюдения, определяемых по формуле
(7)
Если в эксперименте невозможно или нецелесообразно определять среднее квадратическое отклонение составляющих случайной погрешности и определено сразу суммарное среднее квадратическое отклонение, то в формуле (6) n=1.
4. Если случайные составляющие погрешности результата наблюдений представлены доверительными границами ε(Pi ), соответствующими разным вероятностям Pi , то сначала определяют среднее квадратическое отклонение результата измерения с однократным наблюдением по формуле
(8)
где z pi /2 — значения функции Лапласа, Затем вычисляют ε(Р) по формуле (4).
Для суммирования систематической и случайной составляющих погрешностей рекомендуется следующий способ.
,
Если >8,то пренебрегают случайной погрешностью и принимают ΔP=θ(P)
Если 0.8≤≤8, то доверительную границу погрешности результата измерений вычисляют по формуле
ΔP=Kp (9)
где
K (10)
Однократные измерения с приближённым оцениванием погрешности
Для таких измерений в качестве результата измерения принимают значение отсчета х, а оценивание погрешностей проводится на основе нормативных данных о свойствах используемых средств измерений( пределов допускаемой основной погрешности, дополнительных погрешностей и др.).
Поскольку эти данные относятся к множеству средств данного типа, то у конкретного экземпляра прибора, используемого в измерении, действительные свойства могут значительно отличаться от нормированных. Тем не менее, не имея другой достоверной информации о реальных метрологических характеристиках, мы вынуждены проводить оценку погрешности измерения на основе предельных норм. Такие оценки, хотя и грубо, но всё же дают возможность оценить погрешность сверху; но для корректировки результата измерения, для введения поправок они недостаточно надёжны.
Общую схему оценивания погрешностей можно представить следующим образом. Выбрав, исходя из условий измерительной задачи, необходимое средство измерения, уточняют условия измерения и оценивают возможные дополнительные погрешности прибора, возникающие от воздействия влияющих величин.
В результате для оценивания погрешности измерения имеем сведения о погрешностях средства измерения:
предел допускаемой основной погрешности прибора Δпр ;
дополнительные погрешности ψ1, ψ2, ψ3 ….ψm .
Таким образом, задача сводится к суммированию составляющих погрешности Δпр, ψ1, ψ2, ψ3 ….ψm .
Верхняя оценка погрешности результата измерения Δ∑ может быть найдена суммированием составляющих по абсолютным значениям:
Δ∑ = (11)
Более реальная оценка погрешности может быть получена статистическим сложением составляющих погрешности. Поскольку основная и дополнительная погрешности средства измерения заданы границами, то, считая их случайными величинами с равномерным распределением, границы их суммы вычислим по формуле
(12)
Обработка прямых многократных измерений
Равноточные измерения
Рассмотрим прежде всего статистические измерения, при которых многократные измерения проводятся для уменьшения влияния случайных погрешностей. Результат каждого измерения при этом дает оценку измеряемой величины.
Результат наблюдения отличается от истинного значения измеряемой величины из-за случайной Δ и систематической Δс , составляющих погрешности
х`=Х+ Δ+ Δс (13)
Повторяя наблюдение, можно получить информацию о случайной погрешности. О систематической погрешности из этих наблюдений информацию извлечь нельзя. Для оценки систематической погрешности необходимо знать свойства используемых средств измерений, метод измерений и условия измерений.
Основная задача — найти оценку x=f(x)При нормальном распределении наблюдений, как уже отмечалось, оптимальной оценкой распределения Х является среднее арифметическое результатов измерений:
(14)
Сумма квадратов отклонений от среднего арифметического минимальна, поэтому среднее арифметическое значение является эффективной оценкой измеряемой величины.
В общем, алгоритм обработки результатов измерений может быть представлен в следующем виде.
1, Исключают из результатов наблюдений известные систематические погрешности. Если известно, что все результаты наблюдений отягощены одинаковой постоянной систематической погрешностью, ее исключают из результата измерений.
2. Если есть подозрение о наличии грубых погрешностей, то определяют среднее арифметическое и среднее квадратичное значения случайной величины s.
Определяют значения t , и t 2 no формулам
Значения t , и t 2 сравнивают с табличным tn , q
Если они больше tn , q , то результаты xmin и хтах , можно считать анормальными и исключить их из дальнейшей обработки.
3. Вычисляют среднее арифметическое X исправленных результатов наблюдений.
4. Вычисляют оценку среднего квадратичного отклонения результата измерений по формуле
5. Рассчитывают оценку среднеквадратического отклонения среднего арифметического значения по формуле
6. Определяют принадлежность результатов измерений нормальному распределению. При числе результатов измерений п>50 для проверки этой принадлежности используют критерии или
Если n<50, то используют составной критерий или критерий Стьюдента.
6.1 Проверка гипотезы с помощью критерия.
6.1.1 Определяют наименьшее xmin и наибольшее хтах значения результатов измерений.
6.1.2 Определяют размах варьирования R
R= хтах — xmin (15)
Количество интервалов, на которое следует разбивать совокупность результатов измерений, определяется из выражения
r=int (16)
6.1.4 Определяется цена деления интервала с
c=R/r (17)
Цена деления с должна быть больше цены деления прибора, с помощью которого производились измерения.
6.1.5 Данные измерений группируют по интервалам и подсчитывают частоты m i .Если в некоторые интервалы попадает меньше пяти наблюдений, то такие интервалы объединяют с соседними.
6.1.6 Для каждого интервала определяется вспомогательная величина ti , по формуле
6.1.7Определяется плотность нормированного распределения по формуле
или по таблицам нормированного нормального распределения.
6.1.8 Определяют теоретическую частоту mmi в середине каждого интервала по формуле mmi = n*c*/S(x) (18)
6.1.9 Для каждого интервала определяют значение по формуле
= (19)
6.1.10 Определяется значение критерия суммированием значений
= (20)
При обработке результатов измерений без применения ЭВМ, расчеты удобно вести в соответствии со столбцами табл. 3
Таблица З
Номер интервала |
Середина интервала |
Плотность нормального потока распределения Р(t i ) |
Теоретическая частота |
Частота интервала |
= |
|
1 2
r |
||||||
S |
S |
S |
6.1.11 Определяется число степеней свободы к= r -3. Если интервалы объединялись, то число степеней свободы уменьшается. Под r в таком случае понимается количество интервалов с учетом объединения.
6.1.12 Задаются уровнем значимости α=0.05; 0.10; 0.20 и т.д., определяют табличные значения и . Если выполняется условие
<< (21)
то распределение результатов измерений считают нормальным.
6.2 Проверка гипотезы с помощью составного критерия.
6.2.1 Определяется отношение d:
d= (22)
6.2.2 Выбирают уровень значимости критерия (обычно 0.02<q1 <0.1 или в % 2<q1 <10).
6.2.3 Определяют теоретические значения критерия dq 1/2 , d1- q 1/2 по табл. 4
Таблица 4
При q/2, % |
При 1-q/2,% |
|||||
Число наблюдений |
1 |
5 |
10 |
1 |
5 |
10 |
11 |
0,9359 |
0,9073 |
0,8899 |
0,7409 |
0,7153 |
0,6675 |
16 |
0,9137 |
0,8884 |
0,8733 |
0,7452 |
0,7236 |
0,6829 |
21 |
0,9001 |
0,8768 |
0,8631 |
0,7495 |
0,304 |
0,6950 |
26 |
0,8901 |
0,8686 |
0,8570 |
0,7530 |
0,7360 |
0,7040 |
31 |
0,8827 |
0,8625 |
0,8511 |
0,7559 |
0,7404 |
0,7110 |
36 |
0,8769 |
0,8578 |
0,8468 |
0,7583 |
0,7440 |
0,7167 |
41 |
0,8722 |
0,8540 |
0,8436 |
0,7640 |
0,7470 |
0,7216 |
46 |
0,8682 |
0,8508 |
0,8409 |
0,7621 |
0,7496 |
0,7256 |
51 |
0,8648 |
0,8481 |
0,8385 |
0,7636 |
0,7518 |
0,7291 |
6.2.4 Гипотеза о нормальности по критерию d принимается, если d1- q 1/2 ≤d1 ≤ dq 1/2 . В противном случае отвергается.
6.3.1 Критерий 2. Критерий 2 введен дополнительно для проверки «концов» распределений. Считается, что результаты наблюдений соответствуют нормальному распределению, если не более m разностей превзойдет значение tp /2 Sx , где tp /2 квантиль распределения нормированной функции Лапласа, отвечающий вероятности P/2. Вероятность Р определяют по п и q как корень уравнения
(23)
Для нахождения р по заданным п и q составлена табл.5. При 10<n<20 следует принимать т=1, а при 50>n>20 следует принимать т=2.
Таблица 5
Уровень значимости q, % |
||||
n |
m |
1 |
2 |
3 |
10 |
1 |
0.98 |
0.98 |
0.96 |
11-14 |
1 |
0.99 |
0.98 |
0.97 |
15-20 |
1 |
0.99 |
0.99 |
0.98 |
21-22 |
2 |
0.98 |
0.97 |
0.96 |
23 |
2 |
0.98 |
0.98 |
0.96 |
24-27 |
2 |
0.98 |
0.98 |
0.97 |
28-32 |
2 |
0.99 |
0.98 |
0.97 |
33-35 |
2 |
0.99 |
0.98 |
0.98 |
36-39 |
2 |
0.99 |
0.99 |
0.98 |
т.
6.4 Гипотеза о нормальности принимается, если для проверяемой группы измерений выполняются оба критерия.
Уровень значимости составного критерия q=q1 +q2 ,где уровень значимости для критерия 1 ( d — критерия); q 1 — то же для критерия 2.
7. Находят доверительную погрешность результата измерений и доверительный интервал для среднего квадратичного отклонений.
7.1 Нахождение доверительных интервалов при известной точности измерений. Если заранее известна средняя квадратичная погрешность σ*,то доверительный интервал имеет вид
(24)
Значение t = tp определяется по заданной доверительной вероятности Р из условия 2Ф(t)=Р.
7.2 Нахождение доверительного интервала при неизвестной точности измерений. В этом случае используют распределение Стьюдента. Доверительный интервал принимает вид
где к=п-1, а множитель tP ( k ) зависит от доверительной вероятности Р и числа измерений п
7.3 Нахождение доверительных интервалов для средней квадратичной погрешности.
Для нахождения доверительных интервалов для средней квадратичной погрешности используют распределение (табл.6)
7.3.1 Определяют Рв и Рн по формулам
Рв =(1+Р)/2, Рн =(1-P)/2
7.3.2 Определяют число степеней свободы по формуле к=n-1
7.3.3 Для полученных значений Рв и Рн по табл. 6 находят соответственно значение и (см. Приложение 1)
7.3.4 Определяют доверительный интервал для среднеквадратичного отклонения
8. Определяют границы θ неисключенной систематической погрешности. Если известно, что погрешность результата измерений определяется рядом составляющих неисключенных систематических погрешностей, каждая из которых имеет свои доверительные границы, то при неизвестных законах распределения их границы суммарной погрешности находят по формуле
где m— число неиеключенных систематических составляющих погрешностей результата измерения; к — коэффициент, принимаемый равным 1.1 при доверительной вероятности P=q=0.95 и зависящий от числа составляющих неисключенных систематических погрешностей.
9. Определяют соотношение .
Если<0.8 то неисключенными погрешностями пренебрегают и в качестве границы погрешности результата измерения принимают Δ=±tp (k)Sx */
Если , то пренебрегают случайной погрешностью и считают, чтоΔ=θ
Если 0.8<<8 то при определении погрешности D необходимо учитывать и случайную и систематическую составляющую.
- Определяют границу погрешности результата измерений по формуле
, где
,а ,
11. Представляют результат измерения и погрешности для случая симметричных доверительных границ в форме X±Δ.
Неравноточные измерения.
Неравноточные результаты измерений возникают, если заданная величина измерялась средствами измерений: различной точности; одинаковой точности, но при разном числе измерений; одинаковой точности при одинаковом числе измерений, но в различных условиях.
Задача обработки неравноточных измерений состоит в определении достоверного значения измеряемой величины и оценке воспроизводимости измерений.
Пусть некоторая величина, истинное значение которой X измерено многократно и неравноточно, и получены результаты х1 х2 …,х n со средними квадратичными отклонениями s 1 s 2 ,…, sn , то наиболее вероятное значение ХР может быть найдено по формуле
Для удобства вычислений по этой формуле, выбирают некоторый коэффициент с таким расчетом, чтобы отношение было по возможности близко к единице. Величину называют «весом» соответствующих результатов и обозначают через р,с учетом этого выражения формулу для определения ХР можно записать в виде
Величина ХР носит название весовое среднее или общей арифметической средней измерений.
Среднее квадратичное отклонение результатов измерений вычисляется по формуле
а для оценки среднего квадратичного отклонения Sxp весового среднего Хр используется формула
Если значения Si отсутствуют, а известно среднее значение измеряемой величины в i-й серии Х i и количество наблюдений nij „ то весовое среднее можно вычислять по формуле
Обработка результатов косвенных измерений
При косвенных измерениях значение искомой величины получают на основании известной зависимости, связывающей ее с другими величинами, подвергаемыми прямым измерениям.
Рассмотрим простейший случай, когда косвенно измеряемая величина является суммой или разностью величин, определяемых прямыми измерениями, т.е.
Z=X+Y
Так как в результате прямых измерении величин X и Y (после исключения систематических погрешностей) включают в себя некоторые случайные погрешности, то формулу косвенного измерения суммы можно переписать в виде:
Z+ΔZ=X+ΔX+Y+ΔY (25)
где X , У-средние арифметические (или средние взвешенные), полученные при обработке результатов прямых измерений величин X и У;Δ X ,Δ Y – случайные погрешности средних значений величин X и Y ; Z и Δ Z — оценка истинного значения косвенно измеряемой величины и его случайная погрешность.
Z :
Z=X+Y, ΔZ= ΔX+ ΔY
Математическое ожидание оценки Z равно, очевидно, значению искомой величины:
MZ=M(X+Y)=MX+MY=MZ=Z
и ее дисперсия соответственно равна:
Δ(Z)= Δ(Δz )= Δ(Δx + Δy )=M[(Δx + Δy )2 ]=M(Δx 2 + Δy 2 +2 Δx Δy )=M(Δx 2 )+M(Δy 2 )+2M(Δx Δy ) = ΔX+ ΔY+2M(Δ x Δy ) (26)
Из (26) следует, что дисперсия суммы двух слагаемых величин кроме суммы дисперсий этих величин включает еще удвоенное математическое ожидание произведения погрешностей, которое называют корреляционным моментом. Корреляционный момент определяет степень тесноты »линейной» зависимости между погрешностями. Через корреляционный момент выражается безразмерная величина, получившая название коэффициент корреляции rxy
(27)
С учетом формулы (27) уравнение (26) примет вид
Если погрешности измерения величин X и У не коррелированны, то
Очень часто теоретические дисперсии распределения прямых результатов измерений случайных величин X и У неизвестны. В этом случае оценка дисперсии результата косвенных измерений определяется через оценки дисперсий
Оценки коэффициента корреляции вычисляют на основании результатов наблюдений исходных величин:
Чем ближе значение коэффициента корреляции к единице, тем теснее связь между величинами X и Y .
Если rxy >0, то имеет место положительная корреляция, т.е величины X и Y изменяются согласованно в одном направлении — увеличение одной величины влечет за собой увеличение другой.
Если rx у <0, то имеет место отрицательная корреляция — увеличение одной величины сопровождается уменьшением другой.
Если rXY = 0, то величины 1и У некоррелированы.
Если требуется оценить истинное значение величины Z , которая связана со многими величинами Xj (j=l,2,…,m), измеряемыми прямым способом, некоторым уравнением
(в общем случае — нелинейным), то для этого поступают следующим образом.
Рассматривая Z как функцию m переменных Хj запишем ее полный дифференциал
(28)
Каждая из величин Xj измерена с некоторой погрешностью Δхj , Полагая, что погрешности Δх: малы, мы можем заменить dX на Δхj :
В выражении (28) каждое слагаемое представляет собой частную погрешность результата косвенного измерения, вызванную погрешностью Δхj , определения величины Xj . Частные производные носят названия коэффициентов влияния соответствующих погрешностей.
Формула (28) является приближенной, так как учитывает только линейную часть приращения функции, однако в большинстве практических случаев она обеспечивает удовлетворительную точность оценки погрешностей результатов косвенных измерений.
Систематические погрешности Δm Xp , если они определены или известны, используются для определения систематической погрешности Δ m Z с учетом их знаков подстановкой в (28).
Эта же формула используется и для определения предельной погрешности косвенно измеряемой величины по предельным погрешностям аргументов.
Рассмотрим оценки случайных погрешностей результатов косвенных измерений. Предположим, что величины Xj измерены со случайными погрешностями Δj, имеющими нулевые математические ожидания и дисперсии.Найдем выражения для математического ожидания и дисперсии погрешности Δ Z, принимая но внимание (28).
(29)
где rij — коэффициенты корреляции погрешностей всех испытаний; и i, кроме i = j .
Если погрешности ΔXj некоррелированы, что чаще всего имеет место при независимых измерениях, то
В качестве оценки косвенно измеряемой величины принимается величина, рассчитанная по следующей формуле:
Z=F(X1 ;X2 ;….Xm )
Дисперсия этой оценки определяется по формуле (29).
Коэффициенты влияния, приведенные в формулах, в случае нелинейной функции F зависят от значений величин Xj . Коэффициенты влияния определяются подстановкой в выражение частных производных оценок соответствующих параметров, что является дополнительным источником погрешности. При экспериментальном определении коэффициентов влияния также возникает погрешность их определения.
Критерий ничтожных погрешностей.
В высокоточных измерениях погрешности обычно округляют до двух значащих цифр. Например, если в результате многократного измерения получены оценки измеряемой величины X=212.3764мм,то невероятно, чтобы погрешность могла быть определена с пятью значащими цифрами. Очевидно, что запись результата, полученного с калькулятора, не повысит точности. Для экспериментатора, как правило, не будет иметь значения, если будут отброшены три последних разряда числа, а погрешность взята равной 0.018 мм. В этом случае запись результата Х=212.3764±0.018 также некорректна, так как оценка размера имеет, по крайней мере, один избыточный разряд, который подлежит округлению. Правильная запись будет Х=212,376±0,018.
Значит, наибольшее различие в двух значащих цифрах, которое может быть при округлении, составляет 5%. Например, число 104.5 будет округлено до 100, а число 105.5-до 110. При определении суммарной погрешности случайных погрешностей результат выражается формулой
Если в этом равенстве, например, к-я погрешность такая, что
(30)
то этой погрешностью можно пренебречь, так как полученное различие при округлении теряется (число 1.04999 принимается равным 1)
Возведем обе части неравенства (30) в квадрат
получим
получим
;
Эта формула в метрологии называется критерием ничтожных погрешностей, а сами погрешности, отвечающие условию (7.56), называются ничтожными или ничтожно малыми.
Формула (7.56) справедлива и для нескольких погрешностей, если
Использование критерия ничтожных погрешностей при анализе частных погрешностей косвенных измерений позволяет выделить те величины, которыми существенно влияют на погрешность результата. Повышение точности измерения этих величин позволит уменьшить суммарную погрешность.
Совокупные и совместные измерения
При совокупных и совместных измерениях искомые значения физических величин Х1 , Х2 ,…Хm и полученные в i- м опыте в результате прямых или косвенных измерений значения физических величин Xi ,Yj связаны между собой уравнениями вида
Fi (X1 ,X2 ,…Xi ,…Xm , Y1 ,Y2 ,….Ym )=0 (31)
i=1, 2,…. n.
После подстановки в каждое уравнение определённых экспериментальных значений Х1 , Х2 , Х3 ….Хm получаем уравнение
Fi (Y1 ,Y2 ,….Ym )=0. (32)
В уравнении (32) знак равенства имеет условный характер, так как полученные в результате эксперимента коэффициенты содержат погрешности. Поэтому уравнение вида (32) называют условным.
Если уравнения (32) составлены из одноимённых величин, то измерения называют совокупными, если физические величины, входящие в уравнения, имеют разные размерности, измерения называют совместными.
Для того чтобы рассчитать значения искомых величин, достаточно иметь m уравнений, т.е. столько, сколько имеется неизвестных. В этом случае результаты измерений и доверительные границы их погрешностей можно найти методами обработки результатов косвенных измерений. Практически для уменьшения погрешности результата делается значительно больше измерений, чем необходимо для определения m неизвестных.
Из-за ограниченной точности определения коэффициентов, условные уравнения (32) одновременно не обращаются в тождества ни при каких значениях искомых величин. Так как истинные значения искомых величин определить нельзя, то задача сводится к нахождению их оценок Yi , представляющих собой наилучшие приближения к истинным значениям.
Введем в каждое из уравнений (32) такие дополнительные слагаемые Δi , которые превращали бы их в тождества.
Значения ∆i называют невязками, или остаточными погрешностями уравнений. Тогда уравнения (32) можно записать следующим образом :
Fi(Y1,Y2,….Ym)+∆i=0 (33)
Одним из наиболее общих способов отыскания оценок истинных значений измеряемых величин является метод наименьших квадратов.
Согласно этому методу наилучшие оценки величин Yi , будут получены тогда, когда сумма квадратов невязок условных уравнений будет минимальна:
SS=