обработки образцов, аппаратуру контроля и регулирования температуры. При работе с микроскопом записывают его марку, объектив, окуляр, общее увеличение, способ освещения.
Систематическая ошибка обнаруживается чаще всего при испытании эталонного образца. В редких случаях результат измерения исправляют введением поправки, чаще, что более правильно, производят юстировку или отладку прибора для исключения систематической ошибки. Другая трудно устранимая причина систематической ошибки – неконтролируемые изменения условий эксперимента во времени. Эту ошибку можно из систематической превратить в случайную, если последовательность проведения опытов задать с использованием таблицы случайных чисел.
Случайная ошибка присутствует практически во всех измерениях и обнаруживается при многократном измерении одной и той же величины xв идентичных условиях. Эта ошибка отражает рассеивание результатов из-за нерегистрируемых изменений условий эксперимента. Случайные ошибки изучают и оценивают с помощью законов распределения. Если по горизонтальной оси отложить значения измеряемой величины в интервале от x до x x, по вертикали – частоту появления этих значений, то получим распределение измеряемой величины (рис. 2.1). При
уменьшении |
интервала |
x и увеличении числа |
интервалов |
|||||||||||
|
|
|
|
ступенчатая |
гистограмма |
|
переходит |
|||||||
частота |
|
в плавную кривую распределения. |
||||||||||||
|
Рассеяние |
|
|
результатов |
||||||||||
|
измерений |
чаще всего описывается |
||||||||||||
|
нормальным законом распределения. |
|||||||||||||
|
а |
|
х |
Аналитическое |
|
выражение |
этого |
|||||||
|
|
|
||||||||||||
|
|
закона: |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
1 |
|
|
x a 2 |
|
|||||
Рис. 2.1. Кривая нормаль- |
|
|
|
|
||||||||||
|
|
|
|
|
|
|||||||||
ного распределения |
y |
|
|
|
|
exp |
|
|
, |
(1) |
||||
|
|
|
|
2 2 |
|
|||||||||
2 2 |
||||||||||||||
величины x |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
где |
x |
– |
|
значение |
|
случайной |
||||
|
|
|
|
9 |
|
|
|
|
|
|
|
|
|
|
величины; a – математическое ожидание случайной
величины; 2 – дисперсия случайной величины. Приблизительно величину a можно оценить по
ограниченному числу измерений:
a xi /n x,
где x– среднеарифметическое значение.
Дисперсия 2 характеризует размытие результатов измерения, ширину кривой распределения. Поэтому дисперсия может характеризовать точность методики, однородность результатов измерений. При ограниченном числе измерений
(ограниченной выборке) эквивалентом 2 служит выборочная дисперсия
S2 x xi 2 / n 1
Построение доверительных интервалов. Первая задача, которая возникает при оценке результатов измерений, – определение погрешности измерений математического ожидания a по ограниченной выборке. Значение x лишь приблизительно оценивает a , т. е. (x ) a x . Заметим, что также случайная величина и в различных сериях измерений она может быть реализована по-разному. Поэтому, оценивая , задаемся надежностью – доверительной вероятностью, с которой гарантируется появление погрешности, не выходящей за пределы . Если обозначить доверительную вероятность P , то степень риска f P , где– уровень значимости.
Значения всех статистических критериев выражаются через уровень значимости или доверительную вероятность. Обычно в технике принимают = 0,05, что соответствует 95 % надежности. Если число измерений велико и известно
значение дисперсии 2 , то нетрудно определить – доверительный интервал при различных значениях
10
доверительной вероятности [по формуле (1)]. Из формулы (1) следует правило «трех сигм», заключающееся в том, что в
интервале 3 находится 99,7 % |
всех результатов, в |
|||
интервале 2 |
– 95 % и в интервале |
– 68 %. |
||
Если число измерений ограничено, то доверительный |
||||
интервал определяется по формуле |
|
|
||
t mS / |
|
, |
|
(2) |
n |
|
|||
t m – коэффициент Стьюдента; |
S |
– корень квадратный |
||
из выборочной дисперсии; n – число измерений.
Отбраковка резко выделяющихся результатов. На практике часто возникает вопрос об отношении к измерению, результат которого сильно отличается от значений, представленных во всей совокупности. В этом случае руководствуются общим правилом, по которому подобный результат следует учитывать при анализе, если вероятность появления такого результата превышает уровень значимости; если же не превышает, то результат отбраковывают. В большинстве случаев задачи решают с помощью статистических критериев. На основе результатов всей совокупности измерений рассчитывают значение статистического критерия и сравнивают с табличным. Если рассчитанное значение меньше табличного, то гипотезу о принадлежности сомнительного результата к данной генеральной совокупности принимают и результат учитывают, а если нет, результат отбраковывают. Одним из таких критериев является критерий Ирвина:
расч xk 1 |
xk / S , |
(3) |
где xk 1 – |
сомнительный результат; xk |
– результат, |
ближайший по значению к сомнительному; – корень квадратный из выборочной дисперсии. Ниже приведены значения этого критерия для разных значений :
n ……….......... 2 |
3 10 |
20 50 100 400 1000 |
11
|
= 0,05 ......... |
2,8 |
2,2 |
1,5 |
1,3 |
1,1 |
1,0 |
0,9 |
0,8 |
|
= 0,01 ......... |
3,7 |
2,9 |
2 |
1,8 |
1,6 |
1,5 |
1,3 |
1,2 |
|
Если |
расч табл , |
то |
результат |
учитывают, если |
||||
расч табл , то результат отбраковывают.
Исследование связи между двумя рядами измерений
(парная корреляция). Одной из задач любого исследования является определение причинно-следственных связей явлений, т.е. установление закономерностей влияния одних измерений на другие. В условиях статистического разброса результатов измерения прямая функциональная связь через экспериментальные точки теряет смысл, так как и частные результаты, и средние значения «плывут» в некоторой области значений. Для исключения субъективных оценок при решении этих вопросов используют методы корреляционного анализа. В частном случае изучения связей между двумя измеряемыми величинами рассчитывают коэффициент парной корреляции:
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
( |
x |
i |
|
x |
) (yi |
|
y |
) |
|
|
|
|
|
|||
rxy |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
(4) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
n |
|
|
|
n |
|
||||||||||||||
|
|
|
( |
x |
i |
x |
)2 (yi |
y |
)2 |
|
|
|
|||||||
|
|
|
1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|||
где xi , yi – частные результаты измерения величин x и
y .
Заметим, что при условии +1 ≥ rxy ≥ -1 и когда rxy > 0
зависимость y f x возрастающая, а при rxy < 0 – убывающая.
Если rxy = 1, связь однозначная и не вызывает сомнений, если rxy = 0 – корреляция отсутствует. При промежуточных значениях rxy необходимо проверять, существенно ли
отличается от нуля найденный коэффициент корреляции, т. е. не являются ли признаки x и y независимыми. Для этого
12
проверяют статистическую гипотезу rxy = 0 путем построения функции
|
rxy |
n 2 |
|
|||
Tрасч |
|
|
|
|
. |
(5) |
|
|
|
|
|||
|
|
|||||
|
1 rxy2 |
|
||||
Эта функция подчиняется распределению Стьюдента с m n 2степенями свободы. Если Т расч ≥ t m , то гипотеза на данном уровне значимости отвергается, т. е. rxy существенно
отличается от нуля. При Т расч < t m основания отвергать нулевую гипотезу, а значит, и говорить о наличии связи между x и y .
Построение математической модели методом наименьших квадратов. После того как установлено, что связь между x и y существует, определяют эту связь в форме уравнения. Вид уравнения либо выбирают из физических соотношений, либо получают путем аппроксимации (чаще полиномом). В условиях множества экспериментальных точек и их разброса зависимость y f (x)должна проходить наиболее вероятным образом между экспериментальными точками. Анализ нормального закона распределения показывает, что наиболее вероятное положение аппроксимирующей функции такое, при котором сумма квадратов отклонений экспериментальных точек от значений y искомой функции при соответствующих значениях х будет минимальной. Это условие можно записать в виде
yi f (xi ) 2 min .
Рассмотрим простейший, но встречающийся случай аппроксимации зависимости линейной функцией
|
(6) |
достаточно |
часто |
экспериментальной y a bx. Двумя
неизвестными параметрами a и b задают все возможные множества прямых. Единственная прямая с определенными
13
параметрами a |
и b описывается зависимостью y f (x). Для |
|||
нее условие метода наименьших квадратов имеет вид |
|
|||
yi |
(a b(xi ) 2 |
min , |
(7) |
|
где |
xi и |
yi |
– числа, соответствующие |
точкам |
эксперимента. Необходимым условием наличия минимума является равенство нулю производных по a и b от указанной суммы:
|
|
|
|
yi |
|
|
(a bxi ) 2 |
|
0 ; |
|
2 yi |
(a bxi ) 0; |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
a |
||||||
|
|
|
|
yi |
|
|
(a bxi) 2 |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0; |
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
2 xi yi |
(a bxi ) 0. |
|
|
|||||||||||||||||||
|
|
|
|||||||||||||||||||||
|
b |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
Получаем систему уравнений: |
|
|||||||||||||||||||
|
na b |
|
x |
i |
|
|
y |
i |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
(8) |
||||||||||||
|
a |
|
x |
|
b |
|
x2 |
|
|
|
x |
|
|
|
|
||||||||
|
i |
|
y |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
i |
|
|
i |
|
|
i |
|
|
||||||||
из которой находим a и b.
Примеры использования методов математической обработки результатов эксперимента
1.Построение доверительного интервала. Ниже приведены результаты измерения твердости образца из конструкционной стали. Требуется определить доверительный интервал, в который попадает ее среднеарифметическое значение с 95 % – ой вероятностью.
Номера |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
измерений |
|||||||||||
|
|
|
|
|
|
|
|
|
|
||
Твердость |
380 |
470 |
470 |
410 |
450 |
700 |
390 |
420 |
510 |
400 |
|
HV |
|||||||||||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
14 |
|
|
|
|
|
Определяем выборочные значения среднеарифметического и дисперсии:
x xi /n 463; S2 xi x 2 /n 1 959.
Рассчитываем значения доверительного интервала
t mS / 
n .
Значение коэффициента определяем по таблице для уровня значимости = 0,05 и числа степеней свободы m n 1 9 .
t m = 2,26; S 
S 
959 30,97; 
n 
10 3,16
2,26 30,97 /3,16 22,15 22; |
HV |
x |
463 22. |
2. Отбраковка резко выделяющихся результатов. В качестве примера воспользуемся рядом измерений,
полученных в предыдущем примере, где x 463, S2 =959, S 30,97. Из всей совокупности результатов измерений результат №6, который составляет 700 HV , вызывает сомнения. Проверим, не является ли это значение результатом промаха.
Определяем x 463;
S 30,97 .
Определяем расчетное значение критерия Ирвина:
расч |
x6 |
x9 |
|
700 510 |
|
190 |
6,15. |
|
S |
|
|
||||
|
|
30,97 |
30,97 |
|
|||
Находим |
табличные |
значения критерия Ирвина |
|||||
(справочник):
15
табл 1,5.
Поскольку расч табл (6,15 > 1,5), результат № 6 следует отбросить как промах.
3. Исследование связи между двумя рядами измерений (парная корреляция). На образцах из конструкционной стали, обработанных по различным режимам, измерили твердость и предел прочности. Результаты измерений приведены ниже:
Номера |
1 |
2 |
3 |
4 |
5 |
|
6 |
|
|
|
7 |
8 |
|
|
|
9 |
10 |
|||||||
измерений |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
xi |
HV |
263 |
701 |
427 |
531 |
322 |
|
511 |
460 |
602 |
|
|
560 |
380 |
||||||||||
yi |
B , |
530 |
1407 |
867 |
920 |
644 |
|
1040 |
936 |
1270 |
|
1210 |
753 |
|||||||||||
МПа |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Рассчитаем средние значения: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
xi |
| xi HV 4757 |
|
|
|
|
HV |
4757 |
=476; |
|
|||||||||||
|
|
x |
|
|
x |
|||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
n |
10 |
|
|
|
|
|
|||
|
|
|
|
|
yi |
| yi bi |
|
|
|
|
|
|
bi |
958. |
|
|
|
|
|
|||||
|
|
y |
9577 |
|
y |
|
|
|
|
|
||||||||||||||
|
|
n |
|
n |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Определяем для каждого из значений отклонения от |
|||||||||||||||||||||||
среднего |
арифметического, |
|
суммы |
произведений |
||||||||||||||||||||
отклонений и |
суммы квадратов отклонений: xi |
|
; yi |
|
|
; |
||||||||||||||||||
x |
y |
|||||||||||||||||||||||
xi x yi y ; xi x 2 yi y 2 .
Результаты расчетов приведены ниже:
Номера измерений …… 1 |
2 |
3 |
4 |
5 |
|||||
xi |
|
|
……………… |
–213 |
225 |
-49 |
55 |
-154 |
|
|
x |
||||||||
xi |
|
2 |
…………… |
45369 |
50625 |
2401 |
3025 |
23716 |
|
x |
|||||||||
|
|
|
|
|
|
16 |
|
|
|
yi |
|
|
|
|
……………. |
-428 |
|
449 |
-91 |
-38 |
-314 |
|||||||
y |
|
|||||||||||||||||
уi |
|
|
|
2 |
…………… |
183184 |
201601 |
8281 |
1444 |
98596 |
||||||||
у |
||||||||||||||||||
xi |
|
|
y |
|
………. |
91164 |
101025 |
4459 |
-2090 |
48356 |
||||||||
|
x |
y |
||||||||||||||||
Номера измерений |
6 |
7 |
8 |
9 |
10 |
Суммы |
||||||||||||
xi |
|
|
……………. |
35 |
-16 |
126 |
84 |
-96 |
|
|||||||||
|
x |
|
||||||||||||||||
xi |
|
2 |
…………… |
1225 256 |
15876 |
7056 |
9516 |
159065 |
||||||||||
x |
||||||||||||||||||
yi |
|
……………. |
82 |
-22 |
312 |
252 |
-205 |
|
||||||||||
y |
|
|||||||||||||||||
уi |
|
2 |
…………… |
6724 |
484 |
97344 |
63504 |
42025 |
703187 |
|||||||||
у |
||||||||||||||||||
xi |
|
y |
|
……. |
2870 |
352 |
39312 |
21160 |
19680 |
326288 |
||||||||
x |
y |
|||||||||||||||||
В соответствии с методическими указаниями рассчитываем значение коэффициента корреляции:
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
( |
x |
i |
|
x |
) (yi |
y |
) |
|
|
|
|
326288 |
|
|
||||||
rxy |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
n |
|
|
|
|
|
|
n |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
159065 703187 |
|||||||
|
|
|
|
|
( |
x |
i |
|
x |
)2 (yi |
y |
)2 |
|
|
|
|
|
|
|||||||
|
|
|
1 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
326288 |
|
|
|
|
|
326288 |
0,977 . |
|
|
|
|
|||||||||||
398,83 838,56 |
|
|
|
|
|
||||||||||||||||||||
|
334442 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Проверяем статистическую гипотезу об отличии коэффициента корреляции от нуля, для чего рассчитываем значение коэффициента Стьюдента
|
rxy |
n 2 |
0,977 2,84 |
||||
Tрасч |
|
|
|
|
|
|
13,87. |
|
|
|
|
0,045 |
|||
|
|
||||||
|
1 rxy |
|
|||||
Табличное значение коэффициента Стьюдента, определенное для уровня значимости 0,05 и числа степеней
17
свободы m n 2 8, составляет |
t m = 2,31, т.е. при |
Tрасч >> t m на статистически значимом уровне следует
признать наличие корреляции между твердостью и пределом прочности.
4. Построение математической модели методом наименьших квадратов. Поскольку коэффициент парной корреляции весьма близок к единице, целесообразно выбрать линейную модель, описывающую изменение предела прочности от числа твердости:
B a0 bHV .
Коэффициент bможно рассчитать по формуле
|
xi |
|
yi |
|
|
|
|
|
||
b |
x |
y |
|
326288 |
2,0513; |
|||||
xi |
|
2 |
|
|||||||
|
159065 |
|||||||||
x |
|
|||||||||
a y bx 958 2,0513 476 958 976,42 18,42.
Таким образом, B 18,42 2,0513HV .
Лабораторная работа № 3
МЕТАЛЛОГРАФИЧЕСКИЕ ИССЛЕДОВАНИЯ СВАРНЫХ СОЕДИНЕНИЙ
Цель работы
1.Изучить методику подготовки металлографических шлифов сварных соединений.
2.Приобрести практические навыки в изготовлении макро- и микрошлифов.
3.Изготовить макро- и микрошлифы и выполнить их исследование
18