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

Matobrabotka / Matobr20_

.pdf
Скачиваний:
96
Добавлен:
16.04.2015
Размер:
565.26 Кб
Скачать

5. Практические способы расчета случайных погрешностей

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

§ 1. Обработка прямых измерений (алгоритм прямых измерений)

Пусть искомая величина x измерена n раз, для нахождения x и x , рекомендуется записать данные в следующую таблицу и производить расчеты в указанном порядке.

Таблица 1

 

 

 

 

xi

 

xi x

 

(xi x )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

x1

 

x1 x

 

(x1 x )2

 

 

2

 

 

x2

 

x2 x

 

(x2 x )2

 

3

 

 

x3

 

x3 x

 

(x3 x )2

 

 

 

 

…..

 

……

 

……….

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

xn

 

xn x

 

(xn x )2

 

 

 

 

 

 

n

 

 

n

 

 

 

2

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

 

 

i

x

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

i

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

1.

Найти сумму всех xi

xi

 

 

 

 

 

 

 

 

 

 

i

1

 

 

 

 

 

2.

Найти x =

 

 

xi

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.Заполнить третий и четвертый столбцы таблицы.

21

n

x x 2

4. Сосчитать сумму в четвертом столбце

i

i

1

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

 

 

 

 

 

 

 

 

 

 

 

n

 

x)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sx

i

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n(n

1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.

Найти в «таблице коэффициентов Стьюдента» tn 1,P

для данного

 

числа измерений и выбранной вероятности.

 

6.

Определить x

tn 1,P Sx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.

Записать окончательный результат x

x

x x tn

1,P Sx .

 

Пример. Пять раз измерен диаметр проволоки с помощью мик-

рометра. Получены следующие результаты (столбец 2).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

di, мм

 

di

 

,мм

 

 

 

 

(di

 

)2,

 

 

 

d

 

 

 

d

 

 

1

3,90

 

 

0,01

 

 

 

 

 

1

104

 

 

 

2

3,85

 

 

0,06

 

 

 

 

36

104

 

 

 

3

3,88

 

 

0,03

 

 

 

 

9

104

 

 

 

4

3,97

 

 

0,06

 

 

 

 

 

36

104

 

 

 

5

3,95

 

 

0,04

 

 

 

 

 

16

104

 

 

 

 

 

19,55

 

 

 

 

 

 

 

 

 

 

98

104

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=19,55 / 5 =3,910 мм

 

 

 

 

 

 

 

d

 

 

 

 

 

S

 

 

98

10 4

 

 

10 2

 

98

 

0, 022 мм

 

 

d

5

4

 

 

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для доверительной вероятности

 

Р

0,95 и числа измерений n = 5,

коэффициент Стьюдента tn 1,P = 3,2, тогда

d =3,2 0,022= 0,070 мм.

Окончательный результат: d = (3,91 ± 0,07 P=0,95 ) мм.

22

Относительная погрешность d = (0,07 / 3,91) 100 % = 1,8 %.

Возможны другие способы расчета, смотри ПРИЛОЖЕНИЕ §§ 4, 5.

§2. Обработка косвенных измерений. Функция одной переменной. (Формулы переноса ошибок)

Пусть искомая физическая величина Y является функцией измеряемой величины x.

Y f (x)

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

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

Как определить погрешность функции, если известна погрешность аргумента?

Для этого пользуются известным соотношением между дифференциалом функции df (x) и бесконечно малым приращением аргу-

мента dx:

df (x)

df (x)

dx f (x)dx

 

 

dx

Полагая x dx , а Y dY , получаем выражение для погрешности функции:

Y

df (x)

x f (x) x

(17)

dx

 

 

 

где x tp,n 1Sx ,…… fX (x) производная функции f (x) при x

x .

Иногда оказывается удобнее (проще) вычислить сначала относительную погрешность, а уже зная ее, определить доверительный

интервал. Учитывая то, что:

d ln f (x)

1

 

df (x)

легко видеть, что

dx

 

f (x)

 

dx

 

 

 

 

 

23

 

 

 

 

 

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

Y

Y

100%

1

 

df (x)

x 100% ln f (x)

x 100%

(18)

Y

Y

 

dx

 

 

 

 

 

 

§3 Обработка косвенных измерений. Функция многих переменных. (Формулы переноса ошибок)

Вобщем случае искомая физическая величина может быть функцией не одной, а нескольких измеряемых величин, то есть:

Y f (X1,X2 , Xn ) Каждая из величин

X1,X2 ,

Xn

определяется с

соответствующей погрешностью X1,

X2 ,

Xn .

В этом случае

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

 

f

 

2

f

 

2

 

f

 

2

S

S

 

S

 

...

S

(19)

 

X1

 

X2

 

Y

X1

 

X2

 

 

Xn

Xn

 

 

 

 

 

 

 

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

 

f

 

2

f

 

2

f

 

2

Y

X1

 

X2

...

X n

(20)

X1

 

X2

Xn

 

 

 

 

 

 

 

Относительная погрешность функции может быть вычислена по формуле:

Здесь величины X1,X2,…Xn не разные значения одной и той же физической величины, а разные физические величины.

24

 

ln f (X1

..X n )

2

ln f (X1

..X n )

 

2

ln f (X1

..X n )

2

Y

X1

X 2

...

X n

X1

 

X2

 

Xn

 

 

 

 

 

 

 

 

 

(21)

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

Погрешность алгебраической суммы

Пусть функция имеет вид:

Y , Y X1 X2 X3 ...

X n , тогда среднеквадратичная погреш-

ность такой функции будет определяться:

 

 

SY

SX2

1

SX2

2 ...

SX2 n

(22)

а выборочная дисперсия:

 

 

 

 

 

 

 

 

SY2

SX2

1

 

SX2

2 ...

SX2 n

(23)

То есть выборочная дисперсия алгебраической суммы равна сумме выборочных дисперсий отдельных независимых переменных. Обратите внимание, на то, что в выражение для выборочной дисперсии функции все слагаемые входят со знаком «+», незави-

симо от того, с каким знаком соответствующая величина входила в алгебраическую сумму.

Погрешность произведения

Пусть функция имеет вид:

Y X1

X 2 .... X n

или

Y

X1

X 2

...

Xn

Z1

Z2

...

Zm

 

 

 

 

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

25

Y

Y

 

X

1

2

X

2

2

..

X

n

2

Z

2

..

Z

m

2

 

 

 

 

 

 

 

 

1

 

 

Y

 

X1

 

X 2

 

X n

 

Z1

 

Zm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(24)

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

тите внимание, на то, что все слагаемые под корнем берутся со знаком «+», независимо от того в числитель или знаменатель вы-

ражения функции они входили.

Производить расчет по этой формуле обычно гораздо проще, чем по формуле (19), а доверительный интервал искомой величины легко найти: Y Y Y .

Погрешности некоторых элементарных функций

1.

Y

С X , где С=const;

Y С X

2.

Y

X n ;

Y

n X

 

3.

Y

C X ;

Y

ln C

X

§ 4. Два способа оценки погрешности при косвенных измерениях

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

Первый способ

Для каждой измеряемой физической величины проводят серию измерений. Эти измерения (прямые) обрабатывают и находят средние значения Xi и соответствующие им доверительные интервалы

Xi tp,n 1 SXi с одной и той же доверительной вероятностью. Затем по выражению функциональной зависимости Y f (X1,X2 , Xn )

находят среднее значение искомой физической величины, используя найденные средние значения всех измеренных величин Y f (X1 , X2 ,....Xn ) . Затем по формулам переноса ошибок (формула

19 или ее частные случаи) рассчитывают доверительный интервал величины Y.

26

Пример

Задача — определить ускорение свободного падения с помощью математического маятника.

Формула для расчета — g

4 2l

, где l — длина маятника, T

T 2

 

 

период его колебаний.

Измерили 5 раз длину маятника и 5 раз его период. Получили следующие значения, которые обработали, как прямые измерения:

Таблица 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l(м)

l l (м)

l

l

 

2

 

l

0,9644

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,965

0,0006

 

3,6 ×10–7

 

 

 

 

 

 

 

 

 

 

 

2

 

 

10 7

 

2

0,966

0,0016

 

25,6×10–7

 

 

S l

 

 

l

l

52

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n(n

1)

5

4

 

3

0,964

–0,0004

1,6×10–7

 

 

 

 

 

 

 

 

 

 

 

=0,00051

 

 

 

 

 

 

4

0,963

–0,0014

19,6×10–7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3,2×0,00051=0,0016

 

 

5

0,964

–0,0004

1,6×10–7

 

 

 

 

 

l

 

 

Σ

4,822

 

 

 

 

 

52×10–7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т(сек)

Т

 

 

 

 

 

 

)2

 

 

 

 

1,9698

 

 

 

 

 

 

 

T

T

T

 

 

 

 

 

 

1

1,970

0,0002

 

0,04×10–6

 

S T

 

6, 8 10

6

 

0,000583

 

2

1,969

–0,0008

0,64×10–6

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

1,971

0,0012

 

1,44×10–6

 

 

 

T

 

3,2× 0,000583

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=0,0018660,0019

 

 

 

 

4

1,968

–0,0018

3,24×10–6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

1,971

0,0012

 

1,44×10–6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Σ

9,849

 

 

 

 

 

6,8×10–6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получили:

l

0,9644

 

 

 

 

 

 

 

и T

1,9698.

 

 

 

 

 

4 2

 

4 3,14162 0,9644

 

Затем рассчитывают

g

l

9,812342

 

 

 

 

 

 

 

 

 

T 2

1,96982

 

 

 

 

 

и, в соответствии с формулой (7) рассчитывают относительную погрешность g:

g

 

l 2

2

T 2

0, 0016

2

2

0,0019 2

0,002540 0,0025.

l

 

 

T

 

0,9644

 

1,9698

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

27

 

 

 

 

 

Тогда g = g ×δg =9,812342 × 0,0025=0,0249 0,025, и оконча-

тельный результат, соответственно: g (9,812±0,025) м/с2 при

Р=0,95

Второй способ

В случаях, когда по условиям опыта измерения делаются в не воспроизводимых условиях, значения функции, Y f (X1,X2 , Xn )

вычисляют для каждой отдельной серии измерений {Xi } , а затем

полученный ряд значений Yi , обрабатывают по алгоритму прямых измерений.

Пример

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

 

 

 

 

 

 

 

Таблица 4

 

 

 

 

 

 

 

 

 

 

L(м)

Т(сек)

g (м/с2)

g

 

 

g

 

2

g

 

g

маятника

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,965

1,970

9,816453

0,00355

0,126032×10–4

2

1,222

2,222

9,771086

–0,041817

17,48669×10–4

3

1,568

2,510

9,825584

0,01268

1,60795×10–4

4

0,559

1,498

9,834401

0,0215

 

4,62167×10–4

5

0,752

1,739

9,816991

0,00409

0,167152×10–4

Σ

 

 

49,06452

 

 

 

24,0095×10–4

 

 

 

 

 

 

g

9,812903

Sg

24,01

10 4

10

2

24, 01

 

0,0109566 , g = 3,2×0,01096 =0,035

5

4

20

 

 

 

 

 

и, соответственно, окончательный результат:

g (9,813±0,035) м/с2 при Р=0,95

28

§ 5 Метод наименьших квадратов (МНК)

Помимо двух вышеописанных способов оценки погрешности результата при косвенных измерениях, иногда применяют еще так называемый «метод наименьших квадратов» или сокращенно МНК. Этот метод можно использовать, если известен вид функциональной зависимости между измеряемыми физическими величинами, а требуется определить коэффициенты, входящие в эту функцию. В наших лабораторных работах предлагается применять этот метод для определения параметров линейной зависимости.

Пусть в эксперименте можно измерить ряд значений некоторой величины x и, соответствующие им значения, величины y. И пусть при этом известно, что между ними справедлива зависимость вида: y = ax + b. Как известно, такая зависимость графически представляется прямой линией (рис.4). Однако измеренные значения xi и yi включают в себя погрешность и, в результате, не лягут идеально на прямую линию.

Как по данным экспериментальных наблюдений наилучшим образом найти коэффициенты a и b? Графически эта задача сводится к построению прямой, ближе всего лежащей ко всем экспериментальным токам, так как прямая однозначно задается этими коэффициентами (рис.4.).

Для аналитического выражения коэффициентов применяется метод наименьших квадратов.

Утверждается, что наилучшей будет та прямая, сумма квадратов расстояний до которой, от всех экспериментальных точек будет минимальной. Расстояние (вдоль оси y) от точки с координатами xi, yi до искомой прямой

определяется

выражением:

(axi b

yi ) ,тогда

сумма

квадратов расстояний

будет

n

 

 

 

равна:

(axi

b yi )2 .

 

i

1

 

 

Решение задачи на нахож-

дение минимума этого выражения (см. приложение, § 7)

Рис. 4

29

приводит к следующим выражениям для коэффициентов a

и b.

 

 

 

n

 

n

 

n

 

 

 

 

 

 

a

n

xi yi

 

 

xi

2

yi

(25)

 

 

 

i 1

 

i

1

i

1

 

 

 

 

 

 

 

 

n

 

n

 

 

 

 

 

 

 

 

n

x2

 

 

x

 

 

 

 

 

 

 

 

i

i

i

i

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

n

x2

n

n

 

 

n

 

 

 

 

y

 

 

x y

 

x

 

b

 

 

i

i

 

 

i i

2

 

i

(26)

 

 

i

1

i 1

i

1

 

i

1

 

 

 

 

 

 

 

n

 

n

 

 

 

 

 

 

 

 

 

n

x2

 

 

x

 

 

 

 

 

 

 

 

i

i

i

i

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

Дисперсию отклонения экспериментальных точек от прямой — S02 и дисперсию коэффициентов a и b — Sa2 и Sb2 можно вычислить по формулам:

 

 

n

 

n

 

 

 

 

 

n

 

 

 

 

y2

a

 

x y

i

b

y

S02

=

 

i

 

 

 

i

 

 

 

i

i 1

 

i 1

 

 

 

 

 

i

1

 

 

 

n

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sa2

 

 

nS

2

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

n

 

 

 

n

 

 

2

 

 

 

 

 

n

x

2

 

 

 

 

 

x

 

 

 

 

i 1

i

 

i

 

1

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

S

2

 

x

2

 

 

 

 

 

 

Sb2=

 

a

 

i

 

 

 

 

 

 

 

 

i

1

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доверительные интервалы для коэффициентов a ются как обычно:

ap tn 1,PSa

b t

S

b

 

p n 1,P

 

(27)

(28)

(29)

и b определя-

Если график исследуемой зависимости проходит через начало координат, то есть b =0, формулы 25, 27 и 28 существенно упрощаются. В этом случае коэффициент a и его дисперсию можно рассчитать по следующим формулам:

 

n

xi

yi

a

 

i

n

(30)

 

1

 

 

 

 

x2

 

i

1

i

 

 

 

30

Соседние файлы в папке Matobrabotka