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

лабы по оптике

.pdf
Скачиваний:
184
Добавлен:
04.06.2015
Размер:
3 Мб
Скачать
☆

x

1

х

х

 

...х

.

(34)

 

 

10

1

 

2

10

 

 

Оно не совпадает со средним значением средних по сериям измерений в два приёма:

х 1

х

х

 

.

(35)

2

1

 

2

 

 

Величину x можно выразить через

х1

и х2 :

 

x 7x3 3x2

.

36)

 

7 3

 

 

 

Значение x из (36) можно рассматривать средним по десяти измерениям, в результате которых было получено семь раз одинаковое значение x1 и три

раза –

значение

x2 . Числа 7 и 3 называют статистическим весом величин

х1 и х2 .

 

случае, если имеется N величин z (z1, z2,…, zN)

со

В

общем

статическими весами ω1, ω2,…, ωN, наилучшее значение получается из

формулы

 

zi

 

 

 

 

~

 

 

 

 

z

 

 

.

(37)

 

 

 

i

 

 

 

 

 

 

На практике чаще встречается случай непрерывного распределения величины z и при N измерениях получается N значений zi, каждое со своей среднеквадратичной ошибкой:

z1 (z1 ),

z2 (z2 ), ...,

zN (zN ).

(38)

Чтобы получить наилучшее значение ~z из этих измерений, необходимо приписать каждому значению zi статистический вес ωi, выраженный через σ

(zi).

Предположим, что значение zi есть среднее из ni измеренных значений. В этом случае вес ωi пропорционален ni. Будем считать, что первоначальные значения являются равноточными, следовательно, они обладают одинаковыми весами и принадлежат одному и тому же распределению. Пусть среднеквадратичная ошибка распределения есть σ, тогда по (7) имеем:

(zi i )

 

 

 

.

 

(39)

 

ni

 

 

 

 

 

 

 

 

 

 

 

 

 

Для веса ωi с учётом (39) имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

.

(40)

 

(zi

)

 

 

i

i

 

 

 

 

По (37) и (40) найдём наилучшее значение

z :

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

490

 

 

 

2

 

z

 

 

 

 

 

 

 

z

 

 

 

 

 

 

~

 

 

 

i

 

2

(zi )

 

 

 

 

i

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(zi )

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

2

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2

(zi )

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(zi )

z найдем:

По (39) и (40) для среднеквадратичной ошибки величины

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

(z )

 

 

 

 

 

.

ni

1 2 (zi )

1 2 (zi ) 12

(41)

(42)

При одинаковых погрешностях измерений σ(z1) = σ(z2) = =…= σ(zN)

~

наилучшее значение z из (41), называемое средневзвешенным, переходит в среднее арифметическое, а среднеквадратичная ошибка величины σ~z переходит в среднеквадратичную ошибку среднего значения σz .

Приборные погрешности

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

приб

1

.

(43)

 

3

 

 

Стрелочные электроизмерительные приборы по величине допустимой погрешности делятся на классы точности, которые обозначаются на шкале прибора цифрами 0,1; 0,2; 0,5; 1,0; 1,5; 2,5; 4,0. Класс точности прибора показывает величину допустимой погрешности в процентах от значения измеряемой величины, при отклонении стрелки прибора до последнего деления шкалы.

Цифровые приборы (электроизмерительные и прочие) имеют, как правило, погрешность, составляющую 1-2 единицы последнего индицируемого разряда.

Если сведений о допустимой приборной погрешности не имеется, то в качестве нее можно принять половину наименьшего деления шкалы прибора. Например, при измерении длины линейкой с миллиметровыми делениями за допустимую погрешность принимается 0,5 мм.

Окончательный результат многократного измерения содержит в себе как случайную, так и приборную погрешность:

491

x

2 случ (x) 2 пр .

(44)

С увеличением числа измерений случайная погрешность все время уменьшается, а приборная погрешность остается неизменной. Отсюда следует правило: количество измерений должно быть таким, чтобы выполнялось неравенство

(x)случ2

.

(45)

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

Погрешности табличных значений

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

Для физических величин, приводимых в справочниках, иногда указываются погрешности. В этом случае под погрешностями подразумеваются стандартные отклонения для доверительной вероятности 0,68. Если погрешность величины, приводимой в справочнике, не указана, то следует, что погрешность не превышает половины единицы последнего разряда числа.

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

Учет погрешности в записи окончательного результата измерения

Завершением обработки данных многократного прямого измерения при заданной доверительной вероятности являются два числа: среднее значение измеряемой величины и его доверительный интервал. Оба этих числа и есть результат проведенного измерения, стандартная форма записи которого имеет следующий вид:

x x x,

,

n k,

(46)

где – величина заданной доверительной вероятности, а k – число проведенных измерений. Устно выражение (46) звучит так: истинное

492

значение величины х с вероятностью после k измерений лежит в интервалеx около среднего значения x .

Возможности вычислительной техники позволяют при определении величин x и x получить числа с любым количеством цифр. Но заблуждением было бы считать, что количество цифр, полученных вычислением, отражает точность проведенных измерений. Математические расчеты показывают, что погрешность вычисления среднеквадратичной ошибки (x) при n измерениях составляет примерно (n 1) 0,5 . Для n = 10

относительная ошибка (x) достигает 30%. При такой ошибке бессмысленно приводить в погрешности лишние цифры, так как они окажутся

ненадежными. Исходя

из

ошибки (x) в (25÷50)% следует сначала

округлить число (x) ,

а

затем, соответственно, округлить и число x ,

руководствуясь следующими правилами:

1)погрешность округляется до одной значащей (ненулевой) цифры слева, если эта цифра больше 2, или до двух первых цифр в противном случае;

2)если цифра, расположенная за оставляемой меньше 5, то ее просто отбрасывают, иначе цифру увеличивают на единицу;

3)если же отбрасываемая цифра равна 5, то округляют до ближайшего четного числа (правило Гаусса);

4)при округлении среднего значения x справа оставляют цифры тех

разрядов, которые сохранились в погрешности x после ее округления;

5) порядковые множители (вида 10k, где k – целое число), одинаковые для среднего и погрешности среднего, выносятся за общую скобку.

Например, стандартная форма записи окончательных результатов

измерений длины волны излучения с вероятностью 95 % имеет вид

 

(576 4)нм,

0,95,

n 8 ,

(47)

или

мм, 0,95,

n 8.

 

(5,76 0,04) 10 4

 

Порядок действий при вычислении окончательных результатов измерений

Прямые многократные измерения

1.По результатам эксперимента найти среднее значение измерений величины по (2).

2.Вычислить выборочное среднеквадратичное отклонение S результатов отдельных измерений от среднего по (10).

3.Отбросить измерения, в которых отклонение от среднего больше 3δ.

4.Повторить вычисления пунктов 1 и 2 для оставшихся результатов.

5.Оценить среднеквадратичное отклонение окончательного результата по (12), если n≥5. При числе измерений n<5 сначала нужно найти

493

выборочное отклонение среднего, затем, используя коэффициенты Стьюдента из табл.1 для доверительной вероятности α = 0,68, вычислить стандартное отклонение среднего.

6.Определить полное стандартное отклонение с учётом приборной погрешности по (34) и (33).

7.Для заданной доверительной вероятности найти коэффициенты Стьюдента по табл.1 и вычислить границы доверительного интервала.

8.Записать окончательный результат в стандартной форме (36) и (37), предварительно проведя его округление.

Косвенные измерения

1.Рассчитать и записать в стандартной форме окончательные результаты прямых измерений величин, необходимых для искомой величины. Пользоваться доверительной вероятностью α = 0,68.

2.Найти среднее значение искомой величины по средним значениям прямых измерений.

3.Вывести формулу для расчёта погрешности результата косвенного измерения и вычислить значения погрешности по (31) или воспользоваться табл. 2.

4.Привести окончательный результат косвенного измерения к стандартной форме записи.

Пример расчёта погрешностей

В качестве примера рассчитаем погрешность при определении длины волны излучения, пропускаемого светофильтром в интерференционном опыте с бипризмой Френеля. Расчётная формула для длины волны λ имеет вид

 

xd

.

(1*)

(a b)

 

 

 

В (1*) x – ширина интерференционной полосы, определяемая экспериментально с помощью измерительного микроскопа; d – расстояние между мнимыми источниками, определяемое также с помощью микроскопа; a и b –расстояние от источника света до бипризмы и от бипризмы до экрана с интерференционной картиной соответственно. По табл. 2 определим, исходя из функциональной связи λ с измеренными величинами, что удобнее всего расчёт погрешности проводить по формуле для относительной погрешности в частном случае имеющей вид:

δ2 2 δ2 ( x) δ2 (d) δ2 (a b) .

(2*)

Определение среднего значения ширины интерференционной полосы х проводим следующим образом: совмещаем крест нитей измерительного микроскопа с серединой тёмной полосы (интерференционный минимум) в левой части поля зрения микроскопа и снимаем отсчёт А1.Обработку

494

интерференционной картины проводим по тёмным полосам, они в два-три раза уже светлых полос, поэтому совмещение на глаз креста нитей с серединой полосы (наиболее тёмным местом) можно осуществить более точно. Далее перемещаем крест нитей окулярного микроскопа в правую часть поля зрения, считая при этом количество пересеваемых тёмных полос, и через n полос снимаем отсчёт А2 (в соответствии с инструкцией по эксплуатации окулярного микроскопа). Значение n следует выбирать по возможности большим, так как относительная погрешность δ(n) при этом

уменьшается, что следует из её определения.

Измерив расстояние А1-А2,

занимаемое n интерференционными полосами, можно найти х:

x

A2 A1

,

(3*)

 

 

n

 

где α – цена деления измерительного микроскопа. Она определяется предварительно путём наблюдения через микроскоп объект-микрометра, представляющего собой эталонную линейку длиной 1 мм нанесённую на стеклянную пластинку. Совместив крест нитей микроскопа с нулевым делением объект-микрометра, получим отчёт В1, совместив же крест нитей с последующим делением линейки, получим отсчёт В2. Цену деления определим по соотношению

 

α

 

1мм

 

.

 

(4*)

 

В

2

 

В

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

С учётом (4*) х из (3*) примет вид

 

 

 

 

x

A2

A1

 

1мм

.

(5*)

 

n

 

 

B

2

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

Структура формулы (5*) показывает, что расчёт погрешности х удобнее начать с определения относительной погрешности. Тогда по

формуле (2) из табл.2 запишем

δ2 ( х) = δ2(А2-А1)+ δ2(n)+ δ2(1мм)+ δ2(В2-В1). (6*)

По определению

 

2 (A

 

 

 

σ2 (A

2

A )

 

 

δ

 

A

)

 

 

 

 

 

1

.

 

(7*)

 

 

 

 

 

 

 

 

 

 

 

2

1

 

(A

2

A )2

 

 

По формуле 1 из табл. 2:

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ2(А2-А1)= σ2(А2)+ σ 2(А1).

(8*)

Так как микроскоп один, то ошибки отсчётов А1и А2 одинаковы:

 

Тогда (7*) имеет вид

σ (А1)=σ(А2)=σ(А);

 

σ2(А2-А1)=2 σ2(А).

(9*)

 

 

 

 

 

 

 

 

2 2 (А)

 

 

 

 

2 (А

А )

 

 

.

(10*)

 

(

А А )2

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

 

 

 

495

Точно так же получим

 

 

2 2 (А)

 

 

2 (В

 

В )

.

(11*)

 

(В

 

В )2

 

2

1

2

 

 

 

 

 

 

1

 

 

Повторные измерения отрезков (А2-А1) и (В2-В1) показали, что разности (А2-А1) и (В2-В1) остаются неизменными, то есть погрешности отрезков определяются погрешностью микроскопа.

Определим предельную погрешность микроскопа величиной, равной половине минимального деления микроскопа Θ = 0,005, тогда среднеквадратичная погрешность по микроскопу будет σ(А) =1/3Θ ≈ 0,002 деления.

Обратимся к величине n. Можно предположить, что погрешность в определении числа n равна нулю, однако если считать точку А1 за начало отсчёта интерференционных полос, то точка А2 соответствует концу n-й полосы с какой-то погрешностью. Положение А2 определено на глаз, по минимуму освещённости. Глазом же более или менее уверенно можно определить 1/10 небольшого интервала, поэтому предположим σ(n)=0,1.

В соответствии с (6*) необходимо определить погрешность эталонной линейки в 1 мм. Минимальное деление эталонной линейки 0,01 мм. Примем за предельную погрешность 0,5 деления, тогда среднеквадратичная погрешность 1 мм линейки будет σ(1мм) = ⅓·0,5·0,01 = 0,002 мм, а относительная погрешность составит

 

δ(1 мм)=0,002 / 1= 2·10 -3.

(12*)

Предположим, что измерения дали следующие значения:

 

А1=0,28 дел; А2=7,95 дел; n=12; В1=2,38 дел; В2=7,60 дел.

 

Проведём расчёты.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Среднее значение цены деления:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

1мм

 

 

 

0,192 мм

дел

.

(13*)

7,60 2,38

 

 

 

 

 

 

 

 

Относительная погрешность:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ2 (α) δ2 (1мм) δ

2 (В

2

В )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

(14*)

(2 10 3)2

 

 

 

2(2 10 3) 2

4 10

 

 

 

6.

 

(7,60 2,38)2

 

 

 

 

 

 

 

 

Относительная погрешность числа полос:

 

 

 

 

 

 

 

 

δ

 

 

0,1

2

7 10 8.

 

(15*)

 

2 (n)

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Относительная погрешность отсчёта:

2(2 10 3 )2

 

 

 

 

 

2 (А А )

 

 

 

7 10 8.

(16*)

 

 

 

 

 

 

 

 

 

 

2

1

 

 

(7,95 0,28)2

 

 

 

 

 

 

 

 

 

 

 

 

496

Относительная погрешность ширины интерференционной полосы:

 

 

δ2 (

х) = 7·10-8+4·10-6+7·10-5≈7·10-5.

(17*)

Среднее значение ширины интерференционной полосы:

 

 

х

А2 А1

 

 

(7,95

0,28)

0,192 0,122

мм.

(18*)

n

 

2

 

 

 

 

 

 

Из вычислений видно, что основной вклад в ошибку определения ширины интерференционной полосы вносит неопределённость в измерение числа интерференционных полос n.

Рассчитаем d – среднее расстояние между мнимыми источниками и ошибку этой величины. По заданию нужно было несколько раз менять положение фокусирующей линзы и определять d для каждого положения линзы. Было получено 3 значения d: d1=3,47 мм; d2=3,60 мм; d3=3,55 мм.

Определим среднее значение d:

 

 

d=1/3(d1+d2+d3)=3,54 мм.

(19*)

Ошибка величины d содержит в себе как приборную погрешность

σпр(d), так и случайную погрешность многократного измерения σсл

(d).

За

величину приборной погрешности можно взять значение погрешности σ(

х),

полученное при расчёте ошибки ширины интерференционной полосы. Такая оценка σпр(d)= σ( х), следует из того, что измерения d и х проводились одним и тем же прибором. За полную погрешность d примем среднеквадратичную ошибку σ(d), определяемую по правилу

 

(d)

сл

2

(d) пр

2 (d) .

(20*)

Рассчитаем σсл(d). Найдём выборочное среднеквадратичное

отклонение d:

 

 

 

 

 

 

S

(d1 d)2 (d2

d)2

(d3 d)2 3,79 10 2 мм

.(21*)

 

3 2

 

 

 

 

По табл. 1 найдём коэффициент Стьюдента для трёх измерений и доверительной вероятности α=0,68. Этот коэффициент равен t =1,4. Найдём

для заданной

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

по определению равен среднеквадратичной ошибке среднего значения d:

Найдём

σсл (d) = tn··S = 1,4 · 3,79·10-2 =5,3 ·10-2мм.

(22*)

относительную

погрешность

δ(d).

Определим

полную

погрешность d по (17*),(18*),(20*):

 

 

 

 

( d )

( 5,31 10 2 )2 7 10 5 ( 0,122 ) 2 5,31 10 2 мм; (23*)

 

 

2 ( d )

2 ( d )

 

( 5,31 10 2 )2

2,25 10 4.

(24*)

 

 

 

 

 

 

 

d 2

3,542

 

 

Определим среднеквадратичную ошибку среднего значения (a+b). Измерения проводились металлической линейкой, следовательно, можно предположить, что среднеквадратичная ошибка отсчёта равна 1/3·0,3 = 0,1

мм. Тогда для a=b=400 мм:

 

σ2(a+b) = σ2(b)+σ2(a) = 2(0,1)2 .

(25*)

497

Найдём относительную погрешность величины (a+b):

 

2

(a b)

 

2 (a b)

 

( a2

b2 )

 

 

 

(a b)2

 

(a

b)2

 

 

 

 

 

 

 

 

 

 

(1 1)

 

 

3 10 10.

 

 

 

 

(400 400)

2

 

 

 

 

 

 

 

 

 

 

 

Найдём среднее значение длины волны λ по (1*):

0,122 3,54 540 10 6 мм 540 нм. 800

Найдём относительную погрешность δ(λ) по (2*):

δ2(λ) = δ2( х) + δ 2(d) + δ 2(a+b) =

(26*)

(27*)

= 7·10-5+2,25·10-4 +6·10-8 ≈ 30·10-5; δ(λ) = 1,7·10-2.

(28*)

Величины, по которым определили среднее значение λ, имеют погрешность с доверительной вероятностью 0,68, поэтому и абсолютная погрешность λ имеет ту же доверительную вероятность и является среднеквадратичной погрешностью:

( ) ( ) 1,7 10 2 540 9,18 10 6 мм 10 нм.

(29*)

Итак, с вероятностью 68 % истинное значение длины волны излучения

δлежит в интервале (540 ± 10) нм.

Свероятностью 99 % доверительный интервал составит ± 3σ и длина волны будет лежать в пределах (540 ± 30) нм.

Линейная зависимость

Часто встречаются ситуации, когда известна функциональная связь двух физических величин у и х, и требуется найти значение третьей величины, которая присутствует в функции в виде параметра. Например, при рассмотрении интерференционного опыта «Кольца Ньютона» диаметр интерференционного кольца Dк связан с порядком интерференции k соотношением

Dk A k , (48)

где А – известный коэффициент пропорциональности, а – длина волны излучения, используемого в опыте. Из (38) видно, что, возведя обе части в квадрат, мы получим линейную функцию относительно k типа у =Вk. Искомая величина будет входить в коэффициент при переменной k. Параметры линейной зависимости и их погрешности могут быть надежно оценены стандартным методом математической статистики, называемым

методом наименьших квадратов (МНК). Определив параметр В и его погрешность, можно найти и его ошибку.

Суть МНК в следующем. Допустим, имеется n пар измеренных значений (х1, у1), (х2, у2), …, (хn, уn). Предположим, что ошибки содержат

498

лишь величины у. Если величины у и х связаны с параметрами m и с линейной зависимостью у = mx + c, то отклонение в i-м измерении составляет

di yi y yi (mxi c) .

(49)

Наилучшие значения m и c выбираются так, чтобы сумма квадратов

отклонений

 

S ( yi mxi c)

(50)

была минимальной, отсюда и название: метод наименьших квадратов.

 

Таблица 3

Формулы для нахождения параметров прямой линии методом наименьших квадратов

 

 

 

 

Прямая

 

Прямая y mx

y mx c

 

 

 

 

 

 

 

 

 

 

 

m

1

(xi x)

y

 

 

m

xi yi

 

 

D

 

 

 

xi2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c y mx

 

 

 

di yi mxi

 

 

 

 

x

1

xi

 

 

2

(m)

di2

 

 

 

n

 

 

xi2 (n

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

xi yi

 

 

 

 

 

 

 

 

 

xi2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

(m)

 

1 di2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D n 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 (c)

1

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yi

1

yi

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

di

 

yi

mxi

c

 

 

 

 

 

 

Принцип минимизации суммы квадратов отклонений впервые предложен Лежандром в 1806 г.

При решение задачи на

экстремум

выражения (40) получаем,

что

оптимальная прямая проходит

через точку с координатами

 

x 1 xi ,

y

1 yi ,

(51)

 

n

 

n

 

499