Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Основы метрологии и контроль качества в строительстве

.pdf
Скачиваний:
2
Добавлен:
29.11.2025
Размер:
4.3 Mб
Скачать

сможем только со значительной степенью неопределенности. Теория вероятностей и математическая статистика предлагают вероятностный подход к определению погрешностей измерений. В инженерной практике случайную составляющую погрешности и общую погрешность измерений выражают доверительным интервалом (см. рис. 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 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