книги / Методы испытаний и диагностики кабелей и проводов
..pdfностеи, попавших в этот ин |
|
|||||||
тервал; |
п — общее |
число |
из |
|
||||
мерений |
(рис. 12.1). |
измерений |
|
|||||
Когда |
число |
|
||||||
стремится |
к |
бесконечности, |
|
|||||
а 8— к |
нулю, |
величина р,(А) |
|
|||||
стремится |
к функции |
плотно |
|
|||||
сти |
распределения |
вероятно |
|
|||||
стей |
р (А), |
которая |
сможет |
|
||||
быть |
аппроксимирована |
ана |
|
|||||
литической функцией. Эта фу |
|
|||||||
нкция называется законом рас |
|
|||||||
пределения вероятностей. Если |
|
|||||||
на результаты |
измерений ока |
Рис. 12.1. Гистограмма погрешности |
||||||
зывает |
слабое |
влияние |
мно- |
измерений |
жество независимых факторов, то эта функция соответствует нормальному закону распределения.
Для построения гистограммы число измерений должно быть велико. Значение и,- для каждого интервала должно быть не меньше 5— 10. Найденная гистограмма соответствует кон кретному значению ха. Если р (А) зависит от х, то для нахождения этой зависимости необходимо производить мно гократные измерения при различных значениях х. Трудоемкость таких экспериментов очевидна.
Функция р (А) обладает следующими свойствами:
р(А)^0; |
J р(Д)</Д=1; |
|
- |
00 |
|
Д 2 |
|
|
J р (А)^А = Р [ - А 1<А ^А 2] = Рд. |
(12.2) |
Наиболее важными числовыми характеристиками закона распределения вероятностей является математическое ожидание Ас (см. рис. 12.1) и среднеквадратическое отклонение с (или дисперсия а 2):
Ас= J Ар (A)d А; |
(12.3) |
-00
а 2= J ( A - A c) 2p(A)dA. |
(12.4) |
-00
Величину Лс в метрологии принято называть систе матической погрешностью. Систематическую погрешность можно рассматривать как среднее значение погрешностей, которые получаются при многократных измерениях одной и той же величины. В общем случае величина Дс может изменяться как детерминированная функция некоторых
171
аргументов (измеряемой величины А, времени, температуры и других влияющих факторов).
При любом законе распределения погрешность измерения
можно представить в виде суммы |
|
Д= ДС+А'. |
(12.5) |
Величину А' в метрологии принято называть случайной погрешностью. Математическое ожидание величины А' равно нулю.
Если величина Ас известна (например, определена экс периментально), то ее целесообразно исключить из результата измерения, переходя от него к исправленному результату хисп:
*исп=*-Ас. (12.6)
Для достаточно точной оценки систематической погреш ности обычно достаточно произвести 10—20 измерений, тогда
|
Дс« - Z |
(12.7) |
|
п i=i |
|
Отметим следующее принципиальное обстоятельство. Так |
||
как любой результат |
измерений х х случаен, то |
значение Дс, |
полученное. по (12.7), |
также случайно. |
|
Среднеарифметическое значение из п измерений для ве личины х
* с = - i X , . |
( 1 2 . 8 ) |
П1=1 |
|
Тогда оценка среднеквадратического отклонения от среднего
значения |
|
*2= -Г Г £ (*.—Хер)2- |
(12.9) |
П1 1 = 1
Причем s стремится к а, если п стремится к бесконечности. Величина хср является случайной, и ее среднеквадратическое отклонение от математического ожидания Jc (среднего значения)
(тс=стl^/ntts/y/n. |
|
(12.10) |
|
Математическое ожидание |
величины |
х |
|
00 |
|
|
|
х — J |
xp{x)dx, |
|
(12.11) |
—оо |
|
|
|
где р(х)— плотность распределения величины |
х. |
||
Если число измерений п составляет |
5— 10 |
или более, то |
закон распределения величины хср стремится к нормальному независимо от того, какому закону распределения подчиняются единичные измерения величины х .
172
В дальнейшем, если это специально не оговаривается, мы будем считать, что систематическая погрешность в основном исключена (в результате измерений включены поправки) или достаточно мала, т. е. значение Лс много меньше а. Принципи ально важно, что систематические погрешности не могут быть исключены полностью по следующим причинам:
действительное значение ха никогда не равно точно ис тинному я:и;
величина Дс определяется приближенно, так как количество опытов ограничено;
измерительную аппаратуру практически можно проверить только в отдельных точках шкалы х.
Случайную погрешность А' в (12.5) можно уменьшить, повторяя измерения п раз и вычисляя лсср по (12.8). Тогда среднеквадратическое отклонение для хср по (12.10) уменьшается с увеличением п. При использовании автоматических измери тельных систем число измерений п можно получить очень большим. Но не следует делать вывод о том, что таким путем случайную погрешность можно снизить до любого малого значения. При малых промежутках времени между измерениями и по другим причинам может существовать статистическая зависимость случайных погрешностей в каждом измерении между собой.
В качестве количественной меры статистической зависимости между двумя случайными величинами используется коэффици
ент корреляции (см. § 12.3). В |
этом |
случае вместо |
(12.10) |
необходимо использовать более |
общее |
выражение |
|
0 . - 4 = /1 + - 1 г у> |
(12.12) |
||
xAV |
ni<j |
|
|
где г,7— коэффициент корреляции между погрешностями на блюдения i и j (суммирование ведется по всем /</).
При независимых наблюдениях все гц равны нулю и из (12.12) следует (12.10). Наиболее сильной зависимости между погрешностями соответствуют г1у равные единице. В этом случае можно показать, что значение стс будет максимальным и равным ст, т. е. увеличение числа наблюдений не увеличивает точности измерения. В большинстве практических случаев погрешности наблюдений в той или иной степени коррелированы, при этом справедливы неравенства
о/у/п<(5с<а- |
(12.13) |
Различают две составляющие погрешности измерения: ин струментальную, определяемую погрешностью средств измере ния, и методическую, связанную с несовершенством метода измерений. Например, при измерении температуры образца
173
термопарой методическая погрешность определяется рядом факторов (см. гл. 9), которые приводят к тому, что температура термопары отличается от действительной температуры образца. Инструментальные погрешности определяются точностью при бора, измеряющего ЭДС термоэлектродов, и точностью их градуировки.
12.2. ВЫЧИСЛЕНИЕ ПОГРЕШНОСТИ ИЗМЕРЕНИЙ
По способу нахождения числового значения измеряемой величины у измерения подразделяются на прямые, косвенные, совместные и совокупные:
при прямых измерениях искомая величина у равна опытному значению х, т. е. у = х ;
при косвенных измерениях величину у находят на основании имеющейся математической зависимости от наблюдаемых
опытных значений |
л'ь х 2) ..., х , |
т. е. y = f( x и |
х 2, ..., х); |
при совместных |
измерениях |
одновременно |
производятся |
измерения нескольких неодноименных величин для нахождения зависимости между ними;
при совокупных измерениях измеряется несколько одноимен
ных |
величин у и |
у 2, ..., у„, |
которые |
находят |
в результате |
|||
решения |
системы |
уравнений: |
|
|
|
|
|
|
|
|
Л {Уи Уг> •••* Ут Яь |
Ьи ки |
1и |
•••> Pi)=0; |
|||
|
|
Р2(У1,У2, •••5 Ут аъ ^2> к2, |
^2» |
•••, |
Pi) = 0; |
|||
|
|
Рт(У1,У2, |
•••» Ую |
ьт, К , |
|
•, |
Рт) = 0, |
|
где |
у и |
у ......... Уп--измеряемые величины; ah |
bt— величины, |
|||||
определяемые путем прямых измерений; kh |
pt— постоянные |
|||||||
величины. |
|
|
|
|
|
|
При совместных и совокупных измерениях число уравнений обычно больше, чем число искомых величин, т. е. т>п. Тогда систему решают методом наименьших квадратов на ЭВМ.
Примером совместных измерений может служить измерение сопротивления изоляции в зависимости от температуры.
Рассмотрим вычисление погрешности измерений при прямых измерениях.
На практике наиболее часто встречается нормальный закон распределения погрешностей, для которого плотность вероят
ностей /?(Д) выражается |
зависимостью |
|
р(Д)= _ 1 е -(Л-Л.)2/(2а)2, |
(12.14) |
|
|
у/2ка |
|
где Аа— математическое |
ожидание А. |
|
174
Всоответствии с (12.2) вероятность попадания погрешности
винтервал (— Д2)
(12.15)
где 0(z) = |
- t2l2du |
|
2у/п, |
|
о |
Функция Ф (z) называется функцией Лапласа или интегралом вероятностей, для нее составлены таблицы (см. табл. П4
приложения |
1). |
доверительного |
интервала |
|
Обычно |
принимают границы |
|||
Ai = A2, |
тогда |
|
|
|
|
|
Рд= /> [|Д |< Д 1] = Ф(Д1/а). |
(12.16) |
|
Для |
А1/<з= 3 значение Ф равно |
0,997. Таким образом, при |
нормальном законе распределения вероятность того, что по грешность однократного измерения превысит по абсолютному значению 3 ст будет не больше 0,3%. Этой вероятностью часто пренебрегают и называют величину За максимально возможной погрешностью.
Формулы, аналогичные (12.14)— (12.15), можно записать и для самой измеряемой величины х (вместо Д будет л\ вместо Аср—хср величины Д2 и Дх будут отклонениями х от математического ожидания х).
Обычно доверительную вероятность Рд принимают равной 0,95. В особенно ответственных случаях принимают Рд равной
0,99 или |
еще больше приближающейся к единице. Величину |
а = 1 —Рд |
называют уровнем значимости. |
Если задано значение Рд, то из (12.16) по таблице функции Ф можно найти интервал для погрешности Дь в который попадают все измерения с вероятностью Рд, но при этом должно быть известно среднеквадратическое отклонение для нормального закона распределения, которому подчиняется данная совокупность измерений.
Для уменьшения интервала случайных погрешностей по (12.8) вычисляют среднее значение х из п измерений. Сред неквадратическое отклонение для среднего вычисляется по (12.10). Тогда доверительный интервал для среднего вычисля
ется |
по формуле |
|
|
|
|
б = Д= /Ра с, |
или |
Б—/(Рд, /г)——, |
(12.17) |
|
|
|
V” |
|
где |
tp является функцией |
от Рд |
и определяется |
по таблицам |
функции Ф (/р— функция, обратная Ф); а с = а /Л/й .
175
На основании экспериментальных данных из п измерений по (12.9) определяется s, которая является лишь приближенной оценкой а. В этом случае коэффициент /(Рд, п) является функцией Рд и числа измерений п и выбирается по таблицам статистики Стьюдента (см. табл. П1 приложения 1). При увеличении п этот коэффициент стремится к значению, опре деляемому по интегралу вероятностей Ф.
При увеличении числа измерений вероятность распределения случайной величины хср стремится к нормальному закону распределения. Если л = 5-И 0 или более, то для л*ср пред полагается нормальный закон распределения. Методы проверки этой гипотезы показаны ниже.
В практике измерений встречаются и другие законы рас пределения, которые оговариваются в ГОСТ 8.011—72.
Для равномерного распределения функция плотности веро
ятности распределения измерения |
величины л: |
/7(JC)== 0 при х 2<х |
и Xi>x; |
|
(12.18) |
р(х)= \/(х2 — х 1) при x 1^ix ^ x 2. |
При этом плотность вероятности р(х) постоянна в интервале между х х и х2 и равна нулю за пределами этого интервала. Примерами случайной погрешности, имеющей равномерное распределение, является погрешность отсчета по шкале прибора и погрешность квантования измеряемой величины по уровню в цифровых измерительных приборах. Обычно равномерное распределение в пределах допускаемых границ принимают для погрешности измерительного прибора. Иногда равномерное распределение принимают в тех случаях, когда закон рас пределения неизвестен.
Дисперсия равномерного распределения |
|
а 2=Д 2/3, |
(12.19) |
где A„=x2—x = x —xi; х = ( х 1 + х 2) / 2 . |
равна Д/Дм. |
Вероятность попадания в интервал А = х —х |
Соотношение (12.19) использовано в ГОСТ 8.207—76 «Прямые измерения с многократными наблюдениями. Методы обработки результатов измерений» для учета неисключенных систематичес ких погрешностей.
Если закон распределения неизвестен, то |
при A2 = At по |
неравенству Чебышева можно получить |
|
Рд = Р [|Д К Д 1] ^ 1 - а 2/Д?. |
(12.20) |
Полагая, например, Д равным За, получим, что Рд превыша ет 0,9. Неравенство (12.20) дает грубую оценку величины Рд снизу. Например, как уже отмечалось, для нормального закона распределения при Д= 3а значение Рд составляет 0,997.
176
Остановимся на способе исключения из результатов измере ния грубых погрешностей. При этом принимается гипотеза о нормальном законе распределения величины х. Необходимо вычислить хср и s по (12.8) и (12.9). Затем для наибольших отклонений измерений от среднего значения вычислить от носительное отклонение t= \х,—xcp\/s. Далее необходимо вы брать уровень значимости (обычно а = 2,5-г 10%) и по табл. ПЗ приложения 1 найти значение /кр, отвечающее данным а и п. Если t> tKр, то данное измерение можно отбросить. С умень шением а растет tKp, и условие t>tKp выполняется труднее. После того как грубые ошибки будут отброшены, статистичес кая обработка оставшихся результатов измерений производится снова обычным порядком.
Последовательность операций при обработке данных можно определить следующим образом:
1) производится п наблюдений величины х = х{;
2)производится исключение систематической погрешности,
идальнейшие операции производятся с исправленными значе ниями xf;
3) находятся хср и s по (12.8) и (12.9);
4)производится исключение грубых ошибок;
5)снова вычисляются хср и s;
6)производится оценка закона распределения вероятностей;
7)по заданной доверительной вероятности Рл находится доверительный интервал е=А1. Если принят нормальный закон распределения, то в вычисляется по (12.17). При других законах распределения коэффициент tp вычисляется для соответст вующего закона распределения.
При достаточно большом числе измерений коэффициент Стьюдента стремится к коэффициенту tp, получаемому по интегралу вероятностей Ф (например, при Рд=0,95 разница между tp несущественна при я > 15). На рис. 12.2 показаны зависимости коэффициента tp от Рл равномерного закона распределения (кривая 1) и для нормального закона рас пределения (кривая 2). Следует отметить, что для большинства других законов распределения (треугольный, трапециевидный, арксинусный и др.) кривая tp=f(Pa) будет располагаться между
Рис. 12.2. Кривые для определения коэффициента tp при расчете довери
тельного интервала
12 Заказ 1841 |
177 |
кривыми 1 и 2. Поэтому для приближенной оценки коэф фициента tp можно использовать среднюю кривую 3.
При косвенных измерениях необходимое значение у вычис ляют по результатам непосредственных измерений х и х 2, ..., хп:
|
|
|
У=/(хь |
*2> |
•*,«)• |
|
(12.21) |
Оценка среднеквадратического |
отклонения для |
величины |
|||||
у и |
случайных |
погрешностей |
величины |
у |
|
||
|
|
т |
1 |
2 ” |
d f д/ |
(12.22) |
|
|
|
|
|||||
|
|
|
sf + 1 р и ^ г т - а д . |
||||
|
- |
ш |
|
к Ф 1 |
O X k O X i |
|
|
|
|
|
|
|
|
||
где |
Si— оценка |
среднеквадратического |
отклонения |
величины |
xh получаемой методом непосредственного наблюдения; pki— оценка коэффициента корреляции между случайными погреш
ностями аргументов |
хк и |
хг |
[кф1)\ |
Рkl~ |
1 " |
{^kj |
Хк,ср){%ij XjCp)j |
^ |
|||
|
ЗД j=i |
|
|
где xktCp и xr/>cp— средние значения из п проведенных измерений.
Коэффициенты |
— называют |
коэффициентами |
влияния, |
|
df |
|
погрешностями. |
|
|
а величины -^-Si— частными |
|
|||
O X i |
разделяют |
на |
сильнокоррелированные |
|
Погрешности |
||||
(| р |= 1 -т-0,7) и слабокоррелированные (р=0-н0,7). На |
практике |
для сильнокоррелированных зависимостей принимается р = ± 1 , а для слабокоррелированных р = 0.
Зависимые погрешности обычно обусловлены одной общей причиной. Например, если в измерительном устройстве имеется ряд измерительных усилителей, которые питаются от общего источника питания, то при уменьшении напряжения питания коэффициенты усиления всех усилителей будут падать. Воз никающие при этом погрешности отдельных усилителей будут сильнокоррелированными (р= +1).
Если Ptt=pkj=±l, то два слагаемых для к и / под корнем
в (12.22) объединяются в |
одно, |
т. е. |
|
|
s |
f ± 2 |
f f SkSl= |
( f Sk± |
^ Sl |
|
uXjt дх\ |
\dxji |
иХ[ |
Из этого следует, что сильнокоррелированные частные погрешности необходимо сначала сложить алгебраически, а за тем уже сумму включить в общую формулу для оценки s.
Если принимать для слабокоррелированных погрешностей р = 0, то из (12.22) следует
178
s= I |
2 |
(12.23) |
|
i=l |
|
В (12.23) сильнокоррелированные частные погрешности |
|
предварительно просуммированы |
алгебраически, и уже их |
|
Я f |
суммы вошли как отдельные величины — j.. При этом общее |
||
количество слагаемых |
уменьшилось, |
дх-, |
т. е. тх<т. |
||
Доверительный интервал погрешности г определяется по |
||
формуле, аналогичной |
(12.17): |
|
|
г -tpS. |
(12.24) |
Закон распределения оценки я, вычисленной по (12.23), является композицией законов распределения измерений от дельных величин X в формуле (12.21). Эта композиция зависит также от соотношения величин слагаемых частных погреш ностей в (12.23). Если слагаемых больше, чем четыре-пять, то закон распределения погрешностей у близок к нормальному и коэффициент определяется для нормального закона рас пределения. В общем случае приближенную оценку tp можно произвести по рис. 12.2. Для более точных вычислений необ ходим обширный статистический анализ.
Ранее, в § 12.1, было показано, что систематическая погреш
ность, вычисленная по (12.7), является случайной |
и имеет |
свое значение среднеквадратического отклонения. |
|
Если систематическая погрешность не учтена, то интервал |
|
суммарной погрешности составит Ас + ^ра(.г). При |
учете си |
стематической погрешности этот интервал составит tp1>/сг2(.т)+ + ст2(Дс). Следовательно, поправку целесообразно вводить, если
Ас+ tp<j(*) > tp i v/a 2(j:)+a2(Ac), |
(12.25) |
где а(х )— среднеквадратическое отклонение измеряемой вели чины JC; сх(Дс)— среднеквадратическое отклонение систематичес кой погрешности Дс. Коэффициенты tp и tpl выбираются при одинаковой доверительной вероятности в соответствии с за конами распределения величин х и (л*—Дс).
Следовательно, вводить поправку на систематическую по грешность не всегда целесообразно.
Для обнаружения, оценки и исключения систематических погрешностей требуется тщательное изучение конкретных при меняемых методов, средств и условий измерения.
Систематическая погрешность может быть определена путем применения более точного метода и средств измерений и будет равна Ас=л:Ср—хТ, где хт— результат точного измерения.
Оценку систематической погрешности можно произвести расчетным путем. Например, если измеряют вольтметром
171)
напряжение между двумя точками в электрической цепи, то влияние входного сопротивления вольтметра на напряжение при подключении вольтметра можно учесть расчетным путем. Разность напряжений без вольтметра и с вольтметром будет систематической погрешностью.
После анализа методов измерений, возможного влияния различных внешних факторов можно исключить систематичес кие погрешности либо техническим совершенствованием уста новки, либо применением специальных приемов измерения.
Систематическую погрешность можно определить, если в качестве измеряемого значения использовать точную меру, например, если с помощью моста измерить емкость эталонного конденсатора. Тогда разность измеряемого значения и значения меры будет равна систематической погрешности.
После исключения ряда систематических погрешностей и по лучения исправленного значения х остаются неисключенные систематические погрешности, которые по ряду причин рас сматриваются как случайные величины, подчиняющиеся своим законам распределения и имеющие свои границы интервалов. Если известны границы систематической погрешности и нет оснований для конкретной оценки закона распределения, то обычно принимается равномерное распределение, для которого используется формула (12.19).
Оценки среднеквадратического отклонения случайных по грешностей должны быть учтены в качестве слагаемых в фор муле (12.22). Если производятся прямые измерения (jp= JC), то в этой формуле одно слагаемое будет учитывать случайные погрешности, а другие— систематические.
Рекомендации для определения границ суммарной погреш ности даются в ГОСТ 8.207—76.
Если известны только границы отдельных неисключенных систематических погрешностей 0 £, а распределение предполага ется равномерным, то границы неисключенной систематической погрешности при т>4 рекомендуется вычислять по формуле
@ = к |
(12.26) |
i
При доверительных вероятностях 0,9; 0,95; 0,99 коэффициент к соответственно равен 0,95; 1,1; 1,4. При 4 используется спе циальный график, рекомендуемый соответствующим стандартом.
Далее все операции производятся с исправленными величинами измерений х, (введена поправка на систематическую погрешность). Грубые погрешности должны быть исключены (ГОСТ 11002— 73).
Величины х ср и s вычисляют по (12.8) и (12.9). Доверитель ный интервал определяется по (12.24) и коэффициенту tp ста тистики Стьюдента.
180