ФЭА / АИТ / КУРСОВАЯ РАБОТА на тему «Обработка результатов наблюдений и оценка погрешностей»
(автор - student, добавлено - 25-04-2014, 14:07)
СКАЧАТЬ:
КУРСОВАЯ РАБОТА на тему
«Обработка результатов наблюдений и оценка погрешностей»
ВВЕДЕНИЕ. 3 Классификация измерений. 3 Оценка результата измерения. 3 Варианты оценки случайных погрешностей. 4 ОСНОВНАЯ ЧАСТЬ. 7 Измерения с однократными наблюдениями. 7 Прямые однократные измерения с точным оцениванием погрешности. 7 Однократные измерения с приближённым оцениванием погрешности. 10 Обработка прямых многократных измерений. 12 Равноточные измерения. 12 Неравноточные измерения. 17 Обработка результатов косвенных измерений. 19 Критерий ничтожных погрешностей. 22 Совокупные и совместные измерения. 24 ПРАКТИЧЕСКАЯ ЧАСТЬ. 27 ПРИЛОЖЕНИЕ. 30 ВЫВОДЫ.. 31 ЛИТЕРАТУРА.. 32
|
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).
Для суммирования систематической и случайной составляющих погрешностей рекомендуется следующий способ.
Если<0.8, то неисключенным остатком систематической погрешности θ(Р) пренебрегают и окончательно принимают ε(Р) за погрешность результата измерения ΔP при доверительной вероятности Р.
Если >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, и t2 no формулам
Значения t, и t2 сравнивают с табличным 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
Таблица З
Номер интервала |
Середина интервала |
|
Плотность нормального потока распределения Р(ti) |
Теоретическая частота |
Частота интервала |
= |
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 Определяют теоретические значения критерия dq1/2, d1-q1/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-q1/2≤d1≤ dq1/2. В противном случае отвергается.
6.3.1 Критерий 2. Критерий 2 введен дополнительно для проверки "концов" распределений. Считается, что результаты наблюдений соответствуют нормальному распределению, если не более m разностей превзойдет значение tp/2Sx, где 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.3.2 Гипотеза о _нормальности принимается, если число разностей больших tp/2Sx, не превышает т.
6.4 Гипотеза о нормальности принимается, если для проверяемой группы измерений выполняются оба критерия.
Уровень значимости составного критерия q=q1+q2,где уровень значимости для критерия 1 (d-критерия); q1 - то же для критерия 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 со средними квадратичными отклонениями s1 s2,...,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 -оценка истинного значения косвенно измеряемой величины и его случайная погрешность.
Таким образом, из уравнения (25) следует, что Z будет равна сумме оценок X иY, а случайные погрешности Δ Х и Δ Y в сумме дадут случайную погрешность Δ 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(Δx2+ Δy2+2 Δx Δy)=M(Δx2)+M(Δy2)+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) является приближенной, так как учитывает только линейную часть приращения функции, однако в большинстве практических случаев она обеспечивает удовлетворительную точность оценки погрешностей результатов косвенных измерений.
Систематические погрешности ΔmXp, если они определены или известны, используются для определения систематической погрешности Δ mZ с учетом их знаков подстановкой в (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=