
книги из ГПНТБ / Слабкий Л.И. Методы и приборы предельных измерений в экспериментальной физике
.pdfЕсли принятая нами оценка для £ является правильной, то она должна лежать в пределах
- |
а |
, г |
. — |
, |
а |
% |
хр |
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