Основы метрологии и контроль качества в строительстве
.pdf
сможем только со значительной степенью неопределенности. Теория вероятностей и математическая статистика предлагают вероятностный подход к определению погрешностей измерений. В инженерной практике случайную составляющую погрешности и общую погрешность измерений выражают доверительным интервалом (см. рис. 3.3). Доверительный интервал – интервал значений ФВ, который с выбранной доверительной вероятностью Рд накрывает истинное значение измеряемой величины.
Рис. 3.2. Схема влияния различных типов погрешностей измерений на результат измерений и оценку его неопределенности
51
Доверительный интервал
д1 |
|
д2 |
|
|
|
с |
|
xист |
x |
xизм |
X |
xд1 |
x |
|
xд2 |
Левая доверительная |
|
Правая доверительная |
|
|
|
||
граница |
|
|
граница |
Рис. 3.3. Погрешности и доверительный интервал величины:
x – точечная оценка результата измерения; xи – измеренное значение величины X; xист – истинное (действительное) значение величины X; x – абсолютная погрешность измерения величины Х; с – систематическая погрешность измерения величины Х; д1, д2 – левая и правая составляющие доверительного интервала
Законы распределения вероятностей случайной величины. Нормальный закон распределения
Как уже указывалось, случайные погрешности проявляют себя при многократных измерениях. Для описания случайной величины, а таковой будет являться случайная погрешность, в теории вероятностей используется закон распределения плотностей вероятности случайной величины, в соответствии с которым случайная величина характеризуется несколькими детерминированными параметрами. Получить закон распределения плотностей вероятности случайной величины можно только на основе генеральной совокупности случайной величины20. Приблизительное представление о законе распределения можно получить, используя часть генеральной совокупности – выборку, в которой проявляются закономерности, присущие генеральной совокупности.
Во многих практических случаях закон распределения плотностей вероятности случайной величины принимается априори (без проверки) из перечня некоторых наиболее распространенных. К таким относятся: нормальный закон распределения (Гаусса), равномерный и треугольный законы распределения, законы распределения Пуас-
20 Генеральная совокупность – совокупность всех значений случайной величины. В большинстве практических случаев ее получить невозможно, так как она будет представлять собой бесконечное множество значений.
52
сона, Пирсона, Бернулли, экспоненциальный и некоторые другие законы распределения21.
Для наглядного отображения поведения исследуемой случайной величины, представленной выборкой своих значений, строят график – гистограмму или полигон частот, по которым можно ориентировочно определить эмпирический закон распределения и его основные параметры.
Рассмотрим пример построения гистограммы на примере обработки данных испытаний на разрыв партии образцов стальной арматуры. Испытано на разрыв 60 образцов арматурной стали; для
каждого образца определен условный предел текучести Ri = 0,2. Все значения полученного вариационного ряда попали в интервал от Rmin = 1040 МПа до Rmax = 1250 МПа.
1.Производится ранжирование вариационного ряда размером n по возрастанию – все результаты располагаются последовательно по мере их возрастания; первым будет значение 1040 МПа, последним – 1250 МПа.
2.Определяется количество интервалов гистограммы (карманов) k. Существуют разные рекомендации для определения k. Некоторые
из них носят эвристический характер, иные имеют теоретическое обоснование. Известно эмпирическоеправило (формула) Стерджеса
k 1 log2 n 1 3,322lg n, |
(3.6) |
правило Райс |
|
k 23 N . |
|
Часто используют следующую формулу для расчета k: |
|
k n, |
(3.7) |
где n – размер вариационного ряда, в примере n = 60. |
|
21 Существуют статистические способы проверки соответствия выбранного эмпирического закона распределения и имеющихся экспериментальных данных измерений, в частности, по критерию Пирсона [4] и критерию Колмогорова.
53
Используя формулу (3.6), |
получаем k = 7,75; по формуле (3.7) |
|||
k = 6,9. Примем компромиссный вариант – k = 7. |
|
|||
3. Рассчитывается размер интервала (кармана) h: |
|
|||
h |
Rmax Rmin |
1250 1040 30 МПа. |
(3.8) |
|
k |
||||
|
7 |
|
||
4. Рассчитывается заполнение интервалов (карманов) nj, подсчитывается количество значений вариационного ряда, попавших в соответствующий интервал. Значения, оказавшиеся на границе смежных интервалов, включают в левый интервал.
5. Рассчитывается относительная частота значений k и плотность распределения событий Pk:
|
j |
|
n j |
, |
P |
j |
. |
(3.9) |
|
|
|||||||
|
|
n |
|
j |
h |
|
|
|
|
|
|
|
|
|
|
В табл. 3.2 приведены данные примера расчетов параметров гистограммы, на рис. 3.4 по этим данным построена гистограмма. Высота прямоугольников гистограммы принимают равной плотности распределения событий (или относительной частоте), ширина основания – размеру интервала. Площадь Ф фигуры, образованной гистограммой, равна сумме площадей каждого прямоугольника:
|
|
k |
|
k |
|
|
|
Ф |
Pj h j 1. |
(3.10) |
|
|
|
j 1 |
|
j 1 |
|
|
|
|
|
|
Таблица 3.2 |
|
Пример расчетных данных для построения гистограммы |
||||
|
|
|
|
|
|
№ |
Интервал, |
Число значений |
Относительная |
Плотность распре- |
|
частота |
деления событий |
||||
|
МПа |
в интервале nj |
j n j n |
Pj j h 103 |
|
1 |
2 |
3 |
|
4 |
5 |
1 |
1040–1070 |
1 |
|
0,0167 |
0,557 |
2 |
1070–1100 |
4 |
|
0,067 |
2,233 |
54
|
|
|
|
Окончание табл. 3.2 |
|
|
|
|
|
1 |
2 |
3 |
4 |
5 |
3 |
1100–1130 |
9 |
0,15 |
5,000 |
4 |
1130–1160 |
15 |
0,25 |
8,333 |
6 |
1160–1190 |
16 |
0,267 |
8,900 |
6 |
1190–1220 |
10 |
0,167 |
5,567 |
7 |
1220–1250 |
5 |
0,083 |
2,766 |
|
|
nj n 60 |
j 1 |
|
Рис. 3.4. Гистограмма распределения значений условного предела текучести стали (штриховой линией показан график «эквивалентного» нормального
закона распределения)
Если увеличивать размер n вариационного ряда, то размер карманов h гистограммы будет стремиться к нулю, а количество карманов k – к бесконечности. «Ступени» гистограммы будут стано-
55
виться все более мелкими, и в пределе график гистограммы примет вид гладкой кривой, отображающей дифференциальный закон распределения f(x) плотности вероятности анализируемой величины X. Площадь фигуры Ф, ограничиваемой сверху указанной кривой, а снизу осью абсцисс, рассчитываемая по формуле
Ф |
f x dx 1, |
(3.11) |
|
|
|
будет равна вероятности появления любого значения величины, то есть – единице.
Интегрирование функции распределения в пределах аргумента от x1 до x2 дает значения вероятности попадания величины x в указанный интервал, значение доверительной вероятности для доверительного интервала x1 … x2:
x2 |
f (x)dx P |
|
x2 . |
(3.12) |
|
||||
x1 |
|
|
x1 |
|
|
|
|
|
Если какую-либо случайную величину x можно рассматривать как сумму большого числа других случайных величин, каждая из которых вносит в сумму небольшой вклад, то независимо от того, какому закону распределения подчиняются отдельные слагаемые xi , их сумма будет иметь распределение, близкое к нормальному, описываемому формулой Гаусса22:
|
|
1 |
|
|
1 |
x xm 2 |
|
||
f x |
|
|
|
|
|
|
|
||
|
|
2 |
|
|
|||||
|
|
e |
|
|
. |
(3.13) |
|||
|
2 |
|
|||||||
|
|
|
|
|
|
|
|
||
Как видно из формулы (3.13), нормальный закон распределения плотности вероятности случайной величины имеет два детерминированных параметра: xm – математическое ожидание величины X,
и – среднее квадратическое отклонение (СКО) величины X. График функции нормального распределения (рис. 3.5) симметричен относительно оси, проходящей через ее максимум, для аргумента
22 Приведенная формулировка получила название центральной предельной теоремы.
56
x = xm. Точкам перегиба скатов функции соответствуют значения аргумента <xm – > и <xm + >.
Рис. 3.5. Характерный вид графиков функций дифференциальной формы нормального закона распределения плотности вероятности величины x
В метрологической практике применяется функция F(x) интегральной формы нормального закона распределения, описываемая соотношением
|
|
|
1 |
x |
|
1 |
x xm 2 |
|
|||
F x x |
f x dx |
|
|
|
|
|
|
|
|||
|
2 |
|
|
|
|||||||
|
|
e |
|
|
dx. |
(3.14) |
|||||
|
|
||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Характерный вид графиков функции интегральной формы нормального закона распределения плотности вероятности приведен на рис. 3.6.
Функцию закона распределения плотности вероятности удобнее анализировать в безразмерном виде. Для этого все размерные величи-
ны, входящие в формулу (3.12), выразим в единицах , то есть
x xm / |
|
x |
t. |
(3.15) |
|
/ |
|
||||
|
|
|
57
Нормальный закон распределения приобретает вид
t |
|
1 |
|
|
t 2 |
|
1 |
|
|
t 2 |
|
||
|
e |
|
|
|
e |
|
|
(3.16) |
|||||
|
2 |
|
2 . |
||||||||||
|
2 / |
|
|
2 |
|
||||||||
Рис. 3.6. Графики функций интегральной формы нормального закона распределения плотности вероятности величины x
Рассмотрим два важных для практики измерений случая определения доверительной вероятности: для симметричного доверительного интервала (рис. 3.7, а) и для несимметричного (одностороннего) доверительного интервала (рис. 3.7, б).
В первом случае tд2 = –tд1. Доверительная вероятность Pд для данного доверительного интервала будет отображаться на рис. 3.7, а площадью, ограниченной кривой функции φ(t), осью аргумента t и границами доверительного интервала. Симметрия графика функции φ(t) и границ доверительного интервала делает справедливыми соотношения
Pд tд2 |
|
tд2 |
1 |
e |
t 2 |
dt 2Ф tд2 , |
|
||
(t)dt 2 |
|
2 |
|
(3.17) |
|||||
|
|||||||||
tд1 |
|
0 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
58
где Ф(tд) – интеграл вероятностей; tд – коэффициент доверия.
а |
б |
Рис. 3.7. Графическое представление симметричного (а) и несимметричного (б) доверительных интервалов
Левая и правая части д доверительного интервала в единицах величины x определяются по формуле
д = tд . |
(3.18) |
Интеграл вероятностей Ф(t) не может быть выражен через элементарные функции. Его значения, полученные численными методами, приводятся в математических справочниках. Некоторые значения Ф(t), часто используемые на практике, приведены в табл. 3.3.
Таблица 3.3
Значения интеграла вероятности и доверительной вероятности для нормального закона распределения
Доверительный |
|
|
|
|
|
|
интервал tд |
1,64 |
1,96 |
2,58 |
2,81 |
3,00 |
3,29 |
в значениях |
|
|
|
|
|
|
Ф(tд) |
0,4495 |
0,4750 |
0,4951 |
0,4975 |
0,4987 |
0,4995 |
Pд = 2Ф(tд) |
0,9 |
0,95 |
0,99 |
0,995 |
0,997 |
0,999 |
59
Из табл. 3.3 следует, что для доверительной вероятности Pд = 0,95, которая часто принимается в практике инженерных измерений, ширина доверительного интервала составляет д = tд = 1,96 .
Для несимметричного доверительного интервала tд2 ≠ –tд1. Характерным для практики оценки свойств строительных материалов является случай, когда правая граница доверительного интервала отсутствует, то есть tд2 = ∞ (рис. 3.7, б). В этом случае доверительную вероятность можно представить в виде суммы двух составляющих, одна – доверительная вероятность для интервала xд–0, вторая – для интервала 0–∞:
Pд P |
|
0 |
P |
|
P |
|
0 |
0,5. |
(3.19) |
|
|
|
|||||||
|
|
Xд |
|
|
0 |
|
Xд |
|
|
Доверительная вероятность для интервала xд–0 составляет ½ доверительной вероятности Pд для симметричного доверительного интервала xд–xд* (см. рис. 3.7, б).
Пример. Определить левую границу несимметричного доверительного интервала для доверительной вероятности Pд = 0,95. Из формулы (3.19) следует
P |
|
0 |
Pд 0,5 0,95 0,5 0,45. |
(3.20) |
|
||||
|
|
Xд |
|
|
|
|
|
|
|
Таким образом, левая граница доверительного интервала tд будет совпадать с левой границей симметричного интервала для довери-
тельной вероятности Pд = 2·0,45 = 0,90. С учетом данных табл. 3.3 и формулы (3.18) для абсолютных значений величины x получим
xд = xm – tд . |
(3.21) |
Границы доверительных интервалов иногда называют квантильными оценками величины. Квантиль – значение, которое величина x не превышает с заданной вероятностью. Например, 2,5 %-й квантиль совпадает с левой границей, а 97,5 %-й квантиль – с правой границей симметричного доверительного интервала для Pд = 0,95. С левой границей несимметричного (одностороннего) интервала
60
