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

книги из ГПНТБ / Слабкий Л.И. Методы и приборы предельных измерений в экспериментальной физике

.pdf
Скачиваний:
27
Добавлен:
25.10.2023
Размер:
12.68 Mб
Скачать

Если принятая нами оценка для £ является правильной, то она должна лежать в пределах

-

а

, г

. —

,

а

%

хр

s

%

/—

 

у п

 

 

 

у п

или

 

 

 

 

(1.23)

1—х:p - ^ = < x d

+ xp^ =

 

уп

 

 

 

уп

с достоверностью 0,95 или с вероятной погрешностью 0,05.

Если, таким образом, данные неравенства выполняются, то мож­

но сказать, что случайная величина х принадлежит к такой гене­ ральной совокупности < Х > , для которой генеральное среднее есть £. Иначе говоря', при увеличении числа измерений (я-ѵоо) величина х сходится к генеральному среднему

В качестве примера рассмотрим измерение мощности W, выде­ ляемой на сопротивлении при токе / 0 и напряжении на сопротив­

лении

Ö 0.

 

/„ и U0 величина

 

При

данных значениях

 

 

 

 

£=/о£/о.

(1-24)

Полагая,

например, /„ =

10-15 А, U0 = ІО-5 В

получим

І = І О -2“

Вт.

 

 

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

Согласно формуле (В. 1), величина о может быть найдена по тео­

реме Найквиста:

 

o = WwyMa = kTAf.

(1.25)

Примем А/ = 1 гц и число измерений я =

9.

Соотношение (1. 23) показывает, что с достоверностью 0,95 най­ денное из опыта значение х не должно выходить за пределы указан­ ного интервала.

 

Результат эксперимента, таким образом, должен дать

10-2ö_l,96 4,15-10-21

< х =

№изм< 1 0

- 2“ +

1.96

4,15-10-2і

т.

1 /9

 

 

 

 

 

УЁГ

е.

 

 

 

 

 

 

 

7,29- 10-21<И7„з„<12,71 • ІО"81.

 

(1.26)

Допустим, что результаты

измерений

дали

 

 

 

 

10_21Ц7г = дг£= 7 ;

10;

8,5; 8,5;

12;

8; 7,5;

9; 10,5.

В

этом случае W = Елу/я — 9 • 10~21, т. е. условие

(1. 26) выпол­

нено.

 

 

 

 

 

 

Итак, прогнозирование результатов измерений, согласно кото­ рому в 95 случаев из 100 должно выполняться (1. 26), проведено на­ ми для случая всего девяти измерений.

21

Если найденная в результате измерений величина х не удовлет­ воряет (1. 26), т. е. если мы вышли за пределы указанного довери­

тельного интервала, то это означает, что число х = Ц7пзМ 'значимо, т. е. с достоверностью 0,95 средняя мощность не равна ІО-20 Вт. В рассмотренном примере число 0,05 есть ошибка 1-го рода.

Величина коэффициента в (1. 26) (в нашем случае 1,96) зависит от выбранной заранее степени достоверности (мы приняли ее рав­ ной 0,95) и может быть найдена из таблиц (интеграл Ф (х)). Так, например, для значений вероятностей 0,99 и 0,999 будем иметь

£ - 2 ,5 8 - ^ < x < M - 2 ,5 8 - f = ;

(1,27)

у п

у п

 

£ - 3 ,2 9 - 7 = < л <

1+ 3,29-^=.

(1.28)

У п

1/ п

 

£

Отметим, что величина и — '^~ = : распределена нормально со

а/ у п

средним значением ё = 0 и а = 1. «-Критерий значимости позволяет установить различие между случайными ошибками того или иного метода измерений и его систематическими ошибками.

В отличие от случайных ошибок, обусловленных большим чи­ слом не поддающихся учету факторов, ошибка от которых распре­ делена нормально со средним ёо — 0 и данным о, систематиче­ ские ошибки бывают обусловлены лишь незначительным числом фак­ торов и действуют, как правило, в какую-то одну сторону, «смещая» распределение на величину 6 = Ё — t„. Необходимо отметить так­ же, что применение «-критерия значимости возможно, если зара­ нее известно стандартное отклонение о.

Рассмотренный выше «-критерий значимости дает возможность определить, при заранее выбранном нами уровне значимости а (а — малое число, равное, например, 0,05), ту критическую область

возможных значений х или с, для которой вероятность неприня­ тия первоначальной гипотезы (если в действительности она спра­ ведлива) равнаа. Таким образом, данный критерий приводит к оши­ бочному решению об отказе от заведомо верной гипотезы (т. е., что £ = Іо) с вероятностью 0,05, т. е. примерно в 5 случаях из 100. Ошибка, связанная с непринятием заведомо правильной гипотезы, есть ошибка 1-го рода, а вероятность такой ошибки есть а. Наряду с ошибкой 1-го рода обычно вводят понятие ошибки 2-го рода ß, которая равна вероятности ошибочного принятия заведомо невер­ ной гипотезы.

Допустим, что наша гипотеза М{х} =

£0 неверна, т. е. что

Ф ЕоТогда функция я

(£х), определяемая в виде

 

я ( У = Р

X —

Іо

при М {х} = ^

(1.29)

а/V«

 

 

 

 

является мощностью критерия, причем

величина ß =

1 — я (х)

есть ошибка 2-го рода.

 

 

 

 

 

22

Таким образом, функция мощности я (gj) есть вероятность принятия правильной гипотезы для всех £ Ф £п и вероятность при­ нятия неправильной гипотезы при I = | 0.

Полагая

х - г 1 _ х - и

(1.30)

 

 

где

 

 

К

Ф — £о

(1.31)

СГ/і/н"

 

 

можно написать:

 

 

я (£х) = Р {|м + А,х| > ир} = Ф { — Up A.J + Ф ( — Up-I-

(1.32)

x — h

Здесь и — °У"[/л" — функция, введенная выше; Ф (и, X) — интеграл ошибок.

График функции мощности я (X) приведен на рис. 4.

Чем больше величина \ отличается от £0>тем больше значение функции мощности. Малость величины я (X) означает, что при ма­ лом значении X возрастает вероятность принятия неверной гипоте­ зы. Однако в силу малости X это сопровождается лишь незначитель-. ной ошибкой.

Рис. 4. Функция мощности я(Х)

Для экспериментов, в которых ставится задача установить, равен или не равен нулю полезный сигнал при данной степени до­ стоверности, имеет смысл строить график функции ß = 1 — я (£х) для отклонений измеряемой величины от ее нулевого значения. Значения функции я (X) и ß для разных X приведены ниже:

я

1 X 1

0,0

0,5

1,0

1,5

2,0

Ш

0,05

0,079

0,17

0,323

0,516

ß

0,95

0,921

0,83

0,677

0,484

 

1 X 1

2,5

3,0

3,5

4

5

я

(к)

0,705

0,851

0,938

0,979

0,999

ß

0,295

0,149

0,062

0,021

0,001

В качестве примера рассмотрим эксперимент со взвешиванием тела массы т = 10 г при механической добротности измерительной системы Q — ІО2, со 0 = 10 с-1, Т — 300° К и А/ = 1 гц. Допус­ тим, что необходимо определить, имеется ли воздействие малой силы F на измерительную систему.

Примем число проведенных измерений равным 9 и среднее зна­

чение F = X =

4,65 ■10-13 Я.

 

Величина о

для данного

случая равна:

 

с — "j/"4/еГ

А/ ~ 4 -10~ 13Я.

Имея данные значения X и ст, можно поставить вопрос: существует ли предполагаемый эффект изменения веса или не существует? Иначе говоря, необходимо сделать выбор между двумя возможными гипотезами: 1) эффект рагвен нулю и 2) эффект не равен нулю.

Как мы уже знаем, ответ на этот вопрос можно получить по «-критерию значимости, построив соответствующие доверительные

границы для X. Однако насколько можно доверять данному крите­ рию? Можно ли, используя его, допустить ошибку и каков'а вероят­ ность такой ошибки?

Пусть заведомо известно, что эффекта нет, т. е. генеральное сред­ нее Іо = 0. Вероятность принятия ложной гипотезы ( | 0 =^0)

есть ß = 1 — л (І) = 1 — л (х) — 1 — л (А). В рассматриваемом

случае А = —^ = = 3,5 и ß = 0,062.

Таким образом, вероятность ошибиться (приняв вместо правиль­ ной неверную гипотезу) при применении «-критерия в данном слу­ чае очень мала — всего 6,2%. Чем больше А, т. е. чем сильнее раз­ личаются «правильная» и «неправильная» гипотезы, тем более донтоверным является «-критерий значимости, т. е. тем меньше стасовится вероятность совершения ошибки принятия заведомо лож­ ной гипотезы.

 

§ 4. Распределение у2

До сих пор

рассматривались случаи, когда величины а2 и £

для множства

(xj были заранее известны. Используя соотношение

где Up — квантиль, можно находить доверительные границы для X. Квантиль Up для нормированного нормального распределения

— 00

24

может быть найдена из

уравнения

 

Ф(ир) = Р

( 0 < Р < 1 )

(1.33)

— заданное значение

вероятности) с помощью таблиц для ин­

тегральной функции распределения Ф (и). Квантиль ир изменяется

от — оо до poo при изменении Р от нуля до единицы. Для Р =

0,5

Up = 0. В частности, «р-квантиль для х или £ определялась

как

l - U p -у^п< X < ^ + U p ~у п

 

и

 

х — ир

 

где п — число измерении.

Перейдем теперь к случаю, когда величина о2 неизвестна. Параметр ст2 может быть найден с помощью оценочной величины

s2 = ~ l S

(1-34)

і-1

 

или, точнее, с помощью

 

4 - 2

(1-35)

І=І

 

Данные оценки для а2 являются наилучшими. При /г-»- оо s2 s2,

т.е. при большом числе измерений обе эти оценки совпадают. Случайные величины (х,- — |) распределены нормально и имеют '

параметры

распределения (0;

а 2).

 

множества {хД из­

Будем предполагать в дальнейшем, что для

вестна величина | и требуется

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

а2 по найденному из опыта значению

s2.

 

 

Преобразуя (1. 35),

получим

 

9

 

 

s2 =

пт2 —

У1

 

 

(1.36)

 

 

 

і = 1

 

 

 

 

где их — (хг — |)/сг — универсальная

функция, введенная выше.

Напомним,

что случайная величина

и является

нормированной

и имеет плотность

распределения

 

 

 

 

 

р(и)=(2п)

ч'- ехр

2

 

 

 

 

 

 

 

Полагая в

(1.36) S

 

= %2 .будем иметь

 

 

 

 

 

 

 

 

 

(1.37)

25

Наряду с величиной х2можно ввести в рассмотрение величину

= 2 ( ^ ) 2.

(1-38)

i=1

 

которая может быть преобразована таким образом, что для нее функ­ ция распределения будет иметь тот нее вид, что и для %2»но при этом

параметром будет (п— 1) вместо числа п для %2.

Отметим, что функция распределения для генеральной сово­ купности (х2) зависит только от числа п и не зависит, в противопо-

p{xzj

Рис. 5. Плотность вероятности Р(%2) при различных значе­ ниях п

ложность множеству {х}, от величины с. Иначе говоря, если парамет­ рами распределения для {*} являются £ и а2, т. е. если

{*}-*(£; ° 2)>

.

то

{X2}->-(«)•

График р (%2) приведен на рис. 5.

Для функции X2 среднее значение £ равно числу /г*, т. е.

 

 

 

І= М {Х 2} =

S

М{и}} = п*,

(1.39)

 

 

 

= 1.

 

£=і

 

 

поскольку М { и і\2

 

 

 

 

 

Доверительные

границы для

х2 определяются

соотношением

 

 

Р {Хр, < X2 <

Хр,} = Р2-

Рх,

(1.40)

для которого квантили

Хр, и Хр2 вычисляются из функции Р (х2)

и зависят от числа степеней свободы п*.

 

Задаваясь

некоторыми значениями Р х и Р 2>можно, используя

таблицы

X2 =

X {%2Ь

«*)> определить

квантили

Хр, и Хр,- На­

пример,

для

Р х =

0,99

и Р 2 =

0,01 при /г* = 9

Xj>t — 2,088 и

Хр, = 21,666.

При меньшей степени достоверности,

например, для

Р х = 0,95 и Р г =

0,05 при том же числе п* степеней свободы Хр, =

26

=3,325 и xl, = 16,919, т. е. доверительный интервал стал уже. Используя соотношения (1.37) и (1.40), можно написать

 

 

=

0-41)

или

 

 

 

p ( J i - s 2< a 2<

^ s

2 ) = P 2- P 1.

(1.42)

^ УР>

УРі

J

 

Поскольку обычно величина £, с помощью которой можно найти s2, неизвестна, то вместо s2 можно использовать оценку

* = 7Г±12 (* -*)'■ 1=1

Отличие величины s2 от s2 нетрудно уяснить, преобразуя s2 следую­ щим образом:

s2= £

(*і - Ю 2= Е

l(Xi-x) + (x— g)]2= Е (Х і- х )2 + п ( х - ^ ) 2,

і=1

І=1

і=1

 

 

(1.43)

поскольку в силу независимости величин (х; — х)

£ (Хі- х ) ( х - 1 ) = 0.

(1.44)

t=i

 

Из (1.43) следует, в частности, что если х и t мало отличается

друг от друга, то можно приближенно считать, что

 

s2~ S 2.

(1.45)

Отметим еще одно важное обстоятельство. Распределение %2 не яв­ ляется нормальным, но сходится к нему весьма медленно при /і оо.

Как показал Фишер, для случая п > 30 случайная величина -j/2%2 распределена приближенно нормально и имеет параметры

о2 = 1

и I = -j/2п — 1.

по Фишеру, приближенно

равна:

Универсальная

функция и,

 

и ~ У ^х2 — -]/"2/г — 1 (ц^ЗО ).

(1-46)

Можно

написать

также,

что

 

 

 

 

ХІ

~ - ^ - ( ѵ 2 / г — 1 4-Их-р)2

(1.47)

или, более точно,

 

 

 

 

 

 

 

ш

+ иі - » Ѵ w )

(1-48)

27

В частном случае, когда Р 2 и соответственно рзрны 0,025 и 0,975, то для разности двух квантилей можно написать:

= 5(0,025 5(о.975 ^ 2 - 1 ,9 6 /2 д — 1 ^ 5 ,5 / п .

(1.49)

Таким образом, с вероятностью 0,95 выполняется приближенное соотношение

Хо.025 — 5(0,975 ^ 5 , 5 / п .

(1.50)

§ 5. Критерий равенства теоретических дисперсий (критерий Бартлетта)

Часто возникает задача установить, обладают ли различные се­ рии измерений некоторой[величиных одной и той же дисперсией о2. Ответ на этот вопрос позволяет судить о том, были ли допущены систематические ошибки в каких-либо і из /е серий измерений, так как иногда приходится сравнивать серии измерений, проведенных в разное время, в разных условиях, а иногда и различными авто­ рами.

Поскольку теоретическое значение дисперсии для подобного рода измерений должно быть постоянным (например, дисперсия а2 = kTAf для случая измерения шумовой мощности), то естествен­

но, что найденные в каждой серии измерений оценки для af, сг|, ..., а* должны иметь своим генеральным средним истинное значение а 2, равное теоретическому.

Критерий, позволяющий судить о том, имеют ли все оценки сво­ им генеральным средним о2, т. е. позволяющий проверить гипотезу

Gl = 02 = .. . = ак = а2,

был установлен Бартлеттом. Пусть в результате эксперимента име­ ются k серий различных измерений какой-либо величины, каждое из которых состоит из п = /; + 1 замеров, причем для каждой се­

рии известна величина sj — [2 (лу — х)2]/2/;. Введем обозначения:

*і

/ =

*2 = / і ------

(1-51)

*=>

2 и

 

 

1=1

 

Если о? = a I = . . . —а, то случайная' величина s2 имеет парамет­ рами распределения величины / и о2, т. е.

а2 = ^ { / , а2}.

28

Как показал Бартлетт, оценка для распределения %2 для данного случая может быть записана в виде

(1.52)

i=i

где

Это означает, что так же, как

и для х2, имеющей

параметры рас­

пределения /и а2, данная оценка для х2 также имеет

своими пара­

метрами величины f и о2 при

 

10.

 

Поскольку, в соответствии с (1.52), каждому значению k (в дан­

ном случае k = п) соответствуют

вполне определенные значения

Хр, и x L то при условии, что

X2 <

УІ, гипотеза а? = о2 = . . . =

= ok = а2 может быть принята.

Другими словами, для выяснения того, имеют ли все серии про­ веденных измерений одну и ту же дисперсию а2 (или в предположе­ нии, что ст2 должна быть одной и той же для разных измерений, т. е. что в процессе измерений не допущено каких-то систематических ошибок для некоторых из k серий), необходимо:

1.Вычислить X2 по формуле (1.52).

2.Выбрать заданный (например, 5%-ный) уровень значимости для границы между случайным и существенным расхождением срав­ ниваемых величин s?.

3. Если для найденного значения х2 = %1 величина Р (хд) при

а = 0,05 будет меньше числа 5,99

(число 5,99

есть табличное зна­

чение X2 при / = k —- 1 = 2 для а — 0,05), то гипотеза of =

а2 =

= оI принимается с вероятностью

0,95. Если

же Р (хд) >

5,99,

то эта гипотеза отвергается, т. е. различие между найденными из опыта оценками s2, s2 и s2 не является случайным, т. е. данные ве­

личины, по крайней мере, с вероятностью 0,95 не могут быть объеди­ нены.

Другими словами, если Хд — «наблюденное» значение (т. е. зна­ чение, вычисленное по формуле (1.52)), то при 5%-ном уровне зна­ чимости между случайным и существенным расхождением наблю­ денное значение Хд считается не случайным, если Р (Хд) < 0,05

при данном числе степеней свободы /* =

k — 1. В этом случае рас­

хождение между ох, . . ., ok является

значимым, т. е. данные at

не могут быть объединены.

 

Если же Р (Хд) > 0,05, то гипотезу о равенстве можно, с вероят­ ностью 0,95, считать согласующейся с наблюдениями.

Рассмотрим три серии замеров веса образца (k — 3, /* = k — 1 =

29

/

— 2) при числах а замеров в каждой серии соответственно равных 16, 13 и 21:

*<!>

 

5

7

 

3

8

и

9

7

6

9

4

7

5

2

8

 

( 9 )

 

4

8

 

12

7

5

3

4

9

4

8 ' 7

6

14 —

x

i

 

 

*<?>

 

9

11

 

7

10

8

9

11

5

3

7

_ 6

5

10

8

x

( [ )

 

7

6

— —

 

 

 

 

 

 

 

*<?)

 

7

8

 

5

6

10

9

8

 

 

 

 

 

 

 

А

c

 

 

 

 

 

 

 

 

 

“ (1)

 

 

—(2) .

 

 

—(3).

 

„2

1

 

V

,

 

 

9

 

= 6,5;

 

=

7;

 

 

 

 

S L —

n

1 2 J

I

- I <

‘ > ) 2 = 5 , 4 7 ;

 

1 0 , 6 6 ;

2

5,

 

 

 

 

 

 

 

 

 

 

 

=

s 3

=

4 .

 

 

 

 

 

 

 

 

 

 

SJ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим

s2,

f

и X2 =

Хо

по формулам

 

 

 

 

 

 

 

к

f.o2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

Дч + /2s| + fs4 .

 

 

 

 

 

 

 

s 2 =

t’*

1

 

 

 

 

 

 

 

 

 

 

 

к

 

 

h ^ h + h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i = l

f = L + h + /з;

^ 2’3"T“

+ /а + /з) g

c= 1

1

 

3 ( 6 - 1 )

I f:

 

(здесь fi

= n-t

1.)

Подставив значения

= 1,03; xg- = 2,12.

Д і + Д і + Д5

/і + /з + /з

 

— I/l lg S1+ /2 s2 + /3 lg S32] 1 '>

 

1

h

3(6- 1) [L + h + h

/г и sf, полѵчнм s2 = 6,75; / = 47; с =

В соответствии с таблицей для ^-распределения при числе степе­ ней свободы /* = V = k — 1 = 2 и Р (х2) = 0,05 находим %2 = = 5,99. Поскольку XQ = 2,12 < 5,99, то с вероятностью 0,95 гипо­

тезу об объединении этих серий измерений можно принять. Предпо­ ложим теперь, что в результате эксперимента получены сколь угодно близкие значения s-, s|, . . ., s2 и заранее известно, что

°1 = СТ2 = °1-

Для этого случая, как нетрудно видеть, формула для %о Дает

при /1 =

ft =

. . . =

fi

поскольку

S2

«rf sf.Хо

А [kfi lg s2- fik lg s] = 0,

30

Соседние файлы в папке книги из ГПНТБ