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

книги / Основы научных исследований

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
12.77 Mб
Скачать

Для проведения опытов с заданной точностью и досто­ верностью необходимо знать то количество измерений, при котором экспериментатор уверен в положитель­ ном исходе. В связи с этим одной из первоочередных задач при статических методах оценки является уста­ новление минимального, но достаточного числа изме­ рений для данных условий. Задача сводится к уста­ новлению минимального объема выборки (числа изме­ рений) Nmin при заданных значениях доверительного интервала и доверительной вероятности. При выпол­ нении измерений необходимо знать их точность:

А = о0/х,

(Ю-5)

где оь — среднеарифметическое значение

среднеквадра­

тичного отклонения сг, равное oo=o/V~n.

Значение оо часто называют средней ошибкой. Дове­ рительный интервал ошибки измерения А определяется аналогично для измерений р=/(То. С помощью t легко определить доверительную вероятность ошибки измере­ ний из табл. 10.1.

В исследованиях часто по заданной точности А и до­ верительной вероятности измерения определяют мини­ мальное количество измерений, гарантирующих требуе­ мые значения А и рд.

Аналогично уравнению (10.3) с учетом (10.5) можно

получить

 

_

 

|х = о arg <р (рд) =

o j Y п t.

(Ю.6)

При Nmin= n получаем

 

 

 

tfmin = <j2*2/ao =

^

2/A2,

(10.7)

здесь kB— коэффициент вариации

(изменчивости), %;

А — точность измерений, %.

Для определения JVmin может быть принята такая по­ следовательность вычислений: 1) проводится предвари­ тельный эксперимент с количеством измерений п, кото­ рое составляет в зависимости от трудоемкости опыта от 20 до 50; 2) вычисляется среднеквадратичное отклоне­ ние о по формуле (10:1); 3) в соответствии с поставлен­ ными задачами эксперимента устанавливается требуе­ мая точность измерений А, которая не должна превышать точности прибора; 4) устанавливается нор­ мированное отклонение t, значение которого обычно за­ дается (зависит также от точности метода); 5) по фор­ муле (10.7) определяют Nmin и тогда в дальнейшем

281

впроцессе эксперимента число измерений не должно быть меньше JV m in.

Пусть, например, при приемке сооружений комиссия

вкачестве одного из параметров замеряет их ширину. Согласно инструкции требуется выполнять 25 измере­ ний; допускаемое отклонение параметра ±0,1 м. Если

предварительно вычисленное значение о = 0,4м, то мож­ но определить, с какой достоверностью комиссия оцени­ вает данный параметр.

Согласно инструкции Д=0,1 м.

Из формулы (10.7)

можно записать t = У~п. — =

=1,25. В соот-

а

0,4

ветствии с табл. 10.1 доверительная вероятность для t=

 

= 1,25рд=0,79. Это низкая ве­

 

роятность.

Погрешность,

пре­

 

вышающая

доверительный

ин­

 

тервал

2р=0,2 м,

согласно

 

выражению

(10.4) будет встре­

 

чаться один раз из 0,79/(1—

 

—0,79) =3,37, т. е. из

четырех

 

измерений.

Это

 

недопустимо.

 

В связи с этим необходимо вы­

 

числить

минимальное

количе­

 

ство

измерений

с

доверитель­

Рис. 10.1. Кривые распреде­

ной

вероятностью

рд,

равной

0,9 и 0,95.

По форйуле

(10.7)

ления Стьюдента для раз­

личных значений:

имеем

Мат= 0,42 • 1,652/0,12=

1 —п-+» ; 2— 72= 10; 3 — л -2

= 43

измерения

при

рд= 0,90

 

и 64

измерения

при

рд=0,95,

что значительно превышает ус­ тановленные 25 измерений,

иценки измерении с помощью сг и сто по приведенным методам справедливы при п > 30. Для нахождения гра­ ницы доверительного интервала при малых значениях применяют метод, предложенный в 1908 г. английским математиком В. С. Госсетом (псевдоним Стьюдент). Кривые распределения Стьюдента в случае п-*~оо (прак­ тически при я>20) переходят в кривые нормального рас­ пределения (рис. 10.1).

Для малой выборки доверительный интервал

 

рСт = О:0а ст.

(Ю.8)

где

а ст — коэффициент Стьюдента,

принимаемый

по

табл.

10.2 в зависимости от значения

доверительной

ве­

роятности рд.

282

Зная Цст. можно вычислить действительное значение изучаемой величины для малой выборки

хд = X ± цот.

(10.9)

Возможна и иная постановка задачи. По л известных измерений малой выборки необходимо определить довери­ тельную вероятность рл при условии, что погрешность среднего значения не выйдет за пределы ±p.CT. Задачу решают в такой последовательности: вначале вычисля­

ется среднее значение х, оо и а Ст=м.ст/ао. С помощью величины <хст, известного л и табл. 10.2 определяют до­ верительную вероятность.

В процессе обработки экспериментальных данных сле­ дует исключать грубые ошибки ряда. Появление этих ошибок вполне вероятно, а наличие их ощутимо влияет на результат измерений. Однако прежде чем исключить то или иное измерение, необходимо убедиться, что это действительно грубая ошибка, а не отклонение вследст­ вие статистического разброса. Известно несколько ме­ тодов определения грубых ошибок статистического ряда. Наиболее простым способом исключения из ряда резко выделяющегося измерения является правило трех сигм: разброс случайных величин от среднего значения не дол­ жен превышать

Xmax.mln = X i 3 d .

(Ю .10)

Более достоверными являются методы, базируемые на использовании доверительного интервала. Пусть име­ ется статистический ряд малой выборки, подчиняющий­ ся закону нормального распределения. При наличии грубых ошибок критерии их появления вычисляются по

формулам

_______

 

Pi *=* (Xmax— x )h V (n —Т)Щ

 

Рз =

(х — XmIn) / o У(П— 1 )/л,

(10.11)

где Xmax, xmin — наибольшее и наименьшее значения из

лизмерений.

Втабл. 10.3 приведены в зависимости от доверитель­ ной вероятности максимальные значения Ртах, возника­

ющие вследствие статистического разброса. Если р |> > р т а х , то значение хтах необходимо исключать из ста­ тистического ряда как грубую погрешность. При Рг< <Р шах исключается величина xmin. •После исключения грубых ошибок определяют новые значения х и а из (л—1) или (л—2) измерений.

283

 

 

 

 

 

 

Т а б л и ц а

10.3

 

 

Критерий появления грубых ошибок

 

п

 

Ршах ИР"

 

 

 

Рmax ПР"1

 

0,90

0,95

0,99

 

0,90

0,95 |

0,99

 

 

3

1,41

1,41

1,41

15

2,33

2,49

2,80

4

1,64

1,69

1,72

16

2,35

2,52

2,84

5

1,79

1,87

1,96

17

2,38

2,55

2,87

6

1,89

2 ,0 0

2,13

18

2,40

2,58

2,90

7

1,97

2,09

2,26

19

2,43

2,60

2,93

8

2,04

2,17

2,37

20

2,45

2,62

2,96

9

2 ,1 0

2,24

2,46

25

2,54

2,72

3,07

10

2,15

2,29

2,54

30

2,61

2,79

3,16

11

2,19

2,34

2,61

35

2,67

2,85

3,22

12

2,23

2,39

2 ,6 6

40

2,72

2,90

3,28

13

2,26

2,43

2,71

45

2,76

2,95

3,33

14

2,30

2,46

2,76

50

2,80

2,99

3,37

Второй метод установления грубых ошибок основан на использовании критерия В. И. Романовского и приме­ ним также для малой выборки. Методика выявления грубых ошибок сводится к следующему. Задаются дове­ рительной вероятностью рд и по табл. 10.4 в зависимости от п находится коэффициент q. Вычисляют предельно до­ пустимую абсолютную ошибку отдельного измерения

евр = стд.

(1 0 .1 2 )

Если X * т а х > е Пр, ТО ИЗМвреНИе * т а х

исключают из

ряда наблюдений. Этот метод более требователен к очи­ стке ряда.

При анализе измерений можно применять для при­ ближенной оценки и такую методику: вычислить по

(10.1) среднеквадратичное отклонение а; определить

с помощью (10.5) сто; принять доверительную

вероят­

ность рд и найти доверительные интервалы рСт из

(10.8);

окончательно установить действительное значение изме­ ряемой величины хд по формуле (10.9).

В случае более глубокого анализа экспериментальных данных рекомендуется такая последовательность: 1) по­ сле получения экспериментальных данных в виде стати­ стического ряда его анализируют и исключают система­ тические ошибки; 2) анализируют ряд в целях обнару­ жения грубых ошибок и промахов: устанавливают подо­ зрительные значения хтах или хт\п\ определяют средне­ квадратичное отклонение а; вычисляют по (10.11) кри­

284

терии Рь Рг и сопоставляют с рт ах, Ртт, исключают при необходимости из статистического ряда дгтах или Xmin и получают новый ряд _из новых членов; 3) вычисляют

среднеарифметическое х, погрешности отдельных изме­

рений (*—xt) и среднеквадратичное очищенного ряда о; 4) находят среднеквадратичное <то серии измерений, ко­ эффициент вариации kB; 5) при большой выборке зада­ ются доверительной вероятностью рд=ср(0 или уравне­ нием значимости (Г—рд) и по табл. 10.1 определяют /; 6) при малой выборке («<30) в зависимости от приня­

та б л и ц а 10.4

Коэффициент для вычисления предельно допустимой ошибки измерения

Значение q при рд

п

0.95

0,98

0,99 |

0,995

 

2

15,56

3 8,97

7 7 ,9 6

7 79,7

3

4 ,9 7

8 ,0 4

11,46

3 6 ,5

4

3 ,5 6

5 ,0 8

6 ,5 8

14,46

5

3 ,0 4

4 ,1 0

5 ,0 4

9 ,4 3

6

2 ,7 8

3 ,6 4

4 ,3 6

7,41

7

2 ,6 2

3 ,3 6

3 ,9 6

6 ,3 7

8

2,51

3 ,1 8

3,71

5 ,7 3

9

2 ,4 3

3 ,0 5

3 ,5 4

5,31

10

2 ,3 7

2 ,9 6

3,41

5,01

12

2 ,2 9

2 ,8 3

3 ,2 3

4 ,6 2

14

2 ,2 4

2 ,7 4

3 ,1 2

4 ,3 7

16

2 ,2 0

2 ,6 8

3 ,0 4

4 ,2 0

18

2 ,1 7

2 ,6 4

3 ,0 0

4 ,0 7

20

2 ,1 5

2 ,6 0

2 ,9 3

3,9 8

 

1 ,96

2 ,3 3

2 ,5 8

3 ,2 9

 

 

Та б л и ц а

10.5

Результаты

измерений и их

 

обработки

м

 

 

 

 

 

I

•н»

I

1

н*

•К*

 

X

X

 

67

— 8

—7 ,8 3

64

67

—8

— 7 ,8 3

64

68

—7

—6 ,8 3

49

68

— 7

—6 ,8 3

49

69

- 6

- 5 , 8 3

36

70

— 5

— 4 ,8 3

25

71

—4

— 3 ,8 3

16

73

- 2

— 1,83

4

74

— 1

— 0 ,8 3

1

75

0

+ 0 ,1 7

0

76

+ 1

+ 1,17

I

77

+ 2

+ 2 ,1 7

4

78

+ 3

+ 3 ,1 7

9

79

+ 4

+ 4 ,1 7

16

80

+ 5

+ 5 .1 7

25

81

+ 6

+ 6 ,1 7

36

82

+ 7

+ 7 ,1 7

49

92

+ 17

+ 1 7 ,2 7

289

х =

2 =

Проверка

2 =

= 7 4 ,8 3 = — 3 —4 6 ,5

= 7 3 7

 

 

+ 4 6 ,5

 

той доверительной вероятности рА и числа членов ряда п принимают коэффициент Стьюдента аст; с помощью формулы (10.2) для большой выборки или (10.8) для малой выборки определяют доверительный интервал; 7) устанавливают по (10.9) действительное значение ис­ следуемой величины; 8) оценивают относительную по­

285

грешность (%) результатов серии измерений при задан­ ной доверительной вероятности рд:

б = _6^ст100

(10.13)

X

 

Если погрешность серии измерений соизмерима с по­ грешностью прибора Впр, то границы доверительного ин­

тервала

______________

 

 

14,- ] / о 5 < & + [*“ !=>]’

(10.14)

Формулой (10.14) следует пользоваться при а Ст<7о^ЗВПр. Если же а СтСГо>ЗВпр, то доверительный интервал вычис­ ляют с помощью (10.1) или (10.9).

Пусть, например, имеется 18 измерений (табл. 10.5). Если анализ средств и результатов измерений показал, что систематических ошибок в эксперименте не обнару­ жено, то можно выяснить, не содержат ли измерения грубых ошибок. Если воспользоваться первым методом (критерий ртах), то надо вычислить среднеарифметиче­

ское* и отклонение о^При этом удобно пользоваться фор­ мулой x=x'-\-{Xi х')/п, где х' — среднее произвольное число. Для вычисления х, например, принять произволь­ но *'= 75 . Тогда х=7Ъ—3/18=74,83. В формуле (10.1) значение- (*—х{)2 можно найти упрощенным методом:

Х[)2= 2 (xt х ) П ^ .

В данном случае (*— xt)2 = 737—32/18 = 736,5. По (10.1) а=736,5/(18—1) =6,58, коэффициент вариации

кв — 100 = 8,8%. Следовательно,

74,83

р ,------ _ 2,68.

Как видно из табл. 10.3, при доверительной вероятно­ сти рд=0,99 и /г=18 ршах—2,90. Поскольку 2,68<рППХ) измерение 92 не является грубым промахом. Если рд= =0,95, Ртах= 2,58, то значение 92 следует исключить.

Если применить правило Зо, то ■^max, min —74,83±ЗХ Х6,58=94,6...55,09, т. е. измерение 92 следует оставить.

В случае, когда измерение 92 исключается, *=73,8, а=5,15. Среднеквадратичное отклонение для всей серии

286

измерений при п= 18 Оо=6,58/18= 1,55; при очищенном ряде ао=5,15/17=1,25.

Поскольку л<30, ряд следует отнести к малой вы­ борке и доверительный интервал вычисляется с приме­ нением коэффициента Стьюдента а ст. По табл. 10.2 при­

нимается доверительная

вероятность 0,95 и тогда а Ст =

=2,11 в случае п= 18;

а ст=2,12,

если

я=17. Довери­

тельный

интервал

 

при

л= 18

цст=

± 1,55-2,11=3,2;

при

л= 17

цСт = ±

1,25-2,12=2,7.

Действительное

значение

изучаемой

 

величины:

при

л= 18

хя=

= 74,8± 3,2;

при

я= 1 7

хд=73,8±2,7.

Относитель­

ная

погрешность

результатов серии

измерений:

при

п=18

б = (3,2-100)/74,8=4,3 %;

при

л = 17 б=(2,7Х

ХЮ0)/73,8=3,7 %. Таким образом, если принять * j= 92 за грубый промах, то погрешность измерения уменьша­ ется с 4,3 до 3,7 %, т. е. на 14 %.

Если необходимо определить минимальное количест­ во измерений при их заданной точности, проводят серию опытов, вычисляют а, затем с помощью формулы (10.7) определяют Nmщ.

В рассмотренном случае ст=6,58; 6В=8,91 %. Если задана точность Д = 5 и 3% при доверительной вероят­

ности рд= 95% , а Ст = 2 ,1 1 .

Следовательно, при Д = 5 %

Nmin= (8,912-2,112)/52 = 14,

а при Д = 3 % Nmi„ =

= (8,912-2,112)/З2=40.

 

Таким образом, требование повышения точности из­ мерения (но не выше точности прибора) приводит к зна­ чительному увеличению повторяемости опытов.

Во многих случаях в процессе экспериментальных ис­ следований приходится иметь дело с косвенными изме­ рениями. При этом неизбежно в расчетах применяют те или иные функциональные зависимости типа

y = f(xltx2,...,*„)•

(10.15)

Так как в данную функцию подставляют не истинные, а приближенные значения, то и окончательный резуль­ тат также будет приближенным. В связи с этим од­ ной из основных задач теории случайных ошибок явля­ ется определение ошибки функции, если известны ошиб­ ки их аргументов.

При исследовании функции одного переменного пре­

дельные абсолютные епр и относительные бПр

ошибки

(погрешности) вычисляют так:

 

епР = ± е * № ,

(10.16)

6DP = ±dln(x),

(10.17)

887

где f'(x)

— производная функции f(x); din \х) — диффе­

ренциал натурального логарифма функции.

Если

исследуется функция многих

переменных, то

 

П

 

 

(10.18)

 

 

 

 

 

^пр — ±

d | In (х^у х2у ч

я71)|.

(10.19)

В (10.18) и (10.19)

выражения

под

знаком суммы

и дифференциала принимают абсолютные значения. Ме­ тодика определения ошибок с помощью этих уравнений следующая: вначале определяют абсолютные и относи­ тельные ошибки аргументов (независимых переменных). Обычно величина хд± е каждого переменного измерена, следовательно, абсолютные ошибки для аргументов из­ вестны, т. е. e*i, еХ2,...,гХп. Затем вычисляют относитель­

ные ошибки независимых переменных:

Находят частные дифференциалы функции и по фор­ муле (10.18) вычисляют еПр в размерностях функции

f(y) И С ПОМОЩЬЮ (10.19) ВЫЧИСЛЯЮТ бпр, % .

Одной из задач теории измерений является установ­ ление оптимальных, т. е. наиболее выгодных, условий из­ мерений. Оптимальные условия измерений в данном экс­ перименте имеют место при 6 np = 6 npmin. Методика ре­ шения этой задачи сводится к следующему. Если исследуется функция с одним неизвестным переменным, то вначале следует взять первую производную по х, при­ равнять ее нулю и определить х х. Если вторая производ­

ная по Х\ положительна, то функция

(10.15) в случае

х = х { имеет минимум. При наличии

нескольких пере­

менных поступают аналогичным образом, но берут про­ изводные по всем переменным х и ...,хп. В результате минимизации функций устанавливают оптимальную об­ ласть измерений (интервал температур, напряжений, силы тока, угла поворота стрелки на приборе и т. д.) каждой функции f(xu ...,xn), при которой относительная ошибка измерений минимальна, т.е. 6*,=m in.

В исследованиях часто возникает вопрос о достовер­ ности данных, полученных в опытах. Решение такой за­ дачи можно проиллюстрировать примером.

Пусть установлена прочность контрольных образцов бетона до виброперемешцвания R I = R 1±.OQ= 20±

288

±0,5 МПа и прочность_бетонных образцов после вибро­

перемешивания #2=^?2± оо= 23±0,6 МПа. Прирост прочности составляет 15%. Это упрочнение относитель­ но небольшое, его можно отнести за счет разброса опыт­ ных данных. В этом случае следует провести проверку на достоверность экспериментальных данных по условию

- L > 3 .

(10.21)

Оi

 

В данном случае проверяется разница x — R t — /?2= = 3,0 МПа. Ошибка измерения равна о0 = ]/"а| + а\, поэтому

(Я, — R2) / V o’? 4- о* = 3,0/(0,25 + 0,36) = 3,84 > 3.

( 10. 22)

Следовательно, полученный прирост прочности явля­ ется достоверным.

Выше были рассмотрены общие методы проверки экспериментальных измерений на точность и достовер­ ность. Ответственные эксперименты должны быть прове­ рены также и на воспроизводимость результатов, т. е. на их повторяемость в определенных пределах измерений с заданной доверительной достоверностью. Суть такой проверки сводится к следующему. Имеется несколько па­ раллельных опытов (серий). Для каждой серии вычис­

ляют среднеарифметическое значение xt (п — число из­ мерений в одной серии, принимаемое обычно равным 3...4). Далее вычисляют дисперсию А. Чтобы оценить воспроизводимость, рассчитывают критерий Кохрена (расчетный):

т

 

kKP= max D i/^D t,

(10.23)

1

 

где шах А — наибольшее значение дисперсий из числа

т

рассматриваемых параллельных серий т ; 2 А — сумма

дисперсий т серий. Рекомендуется принимать 2 ^ т ^ 4 . Опыты считают воспроизводимыми при

Акр < Лк,,

(Ю-24)

где йкт — табличное значение критерия Кохрена

(табл.

1 0 .6 ), принимаемое в зависимости от доверительной ве­ роятности рд и числа степеней свободы q— n 1 . Здесь т — число серий опытов; п — число измерений в серии.

289