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

Полезные файлы / Теория по погрешностям

.pdf
Скачиваний:
25
Добавлен:
16.12.2017
Размер:
733.03 Кб
Скачать
yi , …

ращение, ax =

f

 

 

 

 

 

 

 

, ay =

f

 

 

 

 

 

 

 

, az =

f

 

 

 

 

 

 

 

 

…–

частные произ-

 

 

 

x

 

x′, y′, z

 

y

 

 

 

 

 

 

 

 

z

 

 

x′, y′, z

 

 

 

 

 

 

x′, y′, z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

водные функции, вычисленные в точке (x′, y′, z).

Рассмотрим вычисление случайной погрешности косвенно определяемой величины f. Для этого вычислим дисперсию ее среднего значения. С учетом (3.2) получим

 

 

S 2f =

 

1

 

fi2 =

1

(ax

xi + ay

yi + az

zi )2 ,

 

 

N (N

1)

N (N −1)

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

S 2f

=

1

 

(ax2

xi2 + a2y

yi2 + az2 zi2 + 2axay xi yi + ...).

 

 

N (N

−1)

 

 

 

 

 

 

 

 

 

 

Если аргументы функции случайны и независимы, то их отклонения от средних значений xi , также независимы. Учитывая, что среднее значение произведения независимых случайных величин равно произведению средних значений сомножителей, получаем, что суммы вида xi yi

равны нулю. Тогда

S

2

= a2S

2

+ a2S

2

+ a2S

2

,

(3.3)

f

x

y

z

 

x

y

z

 

 

где Sx2 , S y2 , Sz2 – дисперсии средних значений аргументов функции. Умножив обе части (3.3) на квадрат коэффициента Стьюдента tP2, N , где N – объем вы-

борок, по которым рассчитываются x′, y′, z′ и Sx2 , S y2 , Sz2 , получим для слу-

чайной погрешности функции

f

= fx2 +

f y2 + fz2 ,

(3.4)

где fx = ax x , f y = ay y ,

fz = az z

частные случайные погрешности

функции.

Смещенное среднее значение функции в (3.2), используя выражение

(3.1), можно выразить через ее истинное среднее значение

 

 

 

 

 

 

 

 

 

 

 

=

 

+ θ( f ) =

 

+ axθ( x) + ayθ( y) + azθ( z) ,

(3.5)

 

 

 

 

 

 

 

 

 

f

f

f

 

 

= f (

 

,

 

,

 

) – истинное среднее значение функции; θ( f )

 

где

f

x

y

z

приборная

погрешность функции. Из (3.5) следует, что истинное среднее значение кос-

32

венно

 

 

измеряемой

 

величины

будет

равно

 

 

=

 

− θ( f ) =

 

 

 

f

f

=

 

a

 

θ

 

a

θ

a

θ

 

, где ни величина, ни знак постоянных прибор-

f

x

( x)

(z)

 

 

 

 

 

y ( y)

z

 

 

 

 

 

 

 

 

 

ных погрешностей θ(x),

θ(y), θ(z) аргументов,

а значит и θ( f ) ,

неизвестны.

Приборные погрешности θ(x), θ(y), θ(z) представляют собой независимые случайные величины. Поэтому, как и в случае нахождения случайной погрешности, для приборной составляющей погрешности получим

q( f ) = q2f( x) + q2f( y) + q2f( z) ,

где q f( x) , q f( y) , q f( z ) частные приборные погрешности косвенно измеряе-

мой величины, откуда для верхней границы приборной погрешности величи-

ны

f получим q f = q2f x + q2f y + q2fz , где q fx =

 

ax

 

qx , q f y =

ay

qy ,

 

 

 

 

 

q fz

=

 

az

 

qz представляют собой верхние границы частных приборных по-

 

 

 

 

грешностей косвенно измеряемой величины, а qx ³ |q(x)|, qy ³ |q(y)|, qz ³ |q(z)| –

верхние границы аргументов функции. Коэффициенты

ax , ay , az имеют

смысл

весовых множителей,

показывающих, с каким

весом случайные

x, y,

z или приборные θx ,

θ y , θz погрешности аргументов функции вхо-

дят, соответственно, в случайную и приборную погрешности функции. Производя суммирование случайной и систематической приборной погрешностей согласно (2.16) получим:

f = f 2 + θ2f = ( fx2 + θ2fx ) + ( f y2 + θ2f y ) + ( fz2 + θ2fz ) ,

откуда полная погрешность величины f будет определяться как

 

 

 

 

 

 

 

 

 

= ax2

 

 

2 + a2y

 

 

2 + az2

 

 

2 ,

(3.6)

 

 

 

 

 

 

 

f

x

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

= x2 + θ2x ,

 

= y2 + θ2y ,

 

=

 

 

z2 + θ2z

 

 

 

x

y

z

 

 

 

 

полные погрешности

аргументов.

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

f = f ± f с вероятностью P = P0 , δf = f f ,

33

где P0 – вероятность определения случайной составляющей погрешности измерения; δf – относительная погрешность косвенно измеряемой величины f.

Числовые значения f и f округляются по тем же правилам, которые сформулированы для прямо измеряемых величин.

Замечание 1. Полученное выражение для полной погрешности величины f (3.6) остается справедливым также в случае различного объема выборок величин x, y, z, ..., являющихся аргументами функции f. При этом пол-

ные погрешности аргументов x , y , z , входящие в f , должны быть определены для одной и той же доверительной вероятности P0.

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

их необходимо складывать по модулю линейно θ f = θ fx + θ f y + θ fz . В

этом случае случайная (3.3) и приборная погрешности функции складываются (объединяются) в полную погрешность функции линейно f = f + θ f .

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

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

ln (axn/(ymtg x)) = ln a + n ln x – m ln у – ln tg x.

В таком случае, используя тождество (ln f )= f f и вводя новые ве-

совые множители b

=

ax

 

, b

 

=

ay

, b

=

az

, получим

 

 

 

 

 

 

 

 

x

 

 

f

y

 

 

f

 

z

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

34

f = f (bx x)2 + (by y )2 + (bz z )2 ,

где b =

d ln f

 

 

 

 

 

 

 

, b

 

=

d ln f

 

 

 

 

 

 

 

, b =

d ln f

 

 

 

 

 

 

 

 

 

,

 

.

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в точке

x

,

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

dx

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

z

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x , y ,z

 

 

 

 

 

 

 

 

 

 

 

 

 

x , y ,z

 

 

 

 

 

 

 

 

 

 

x , y ,z

 

 

 

 

3.2. Выборочный метод

Выборочный метод расчета погрешностей применяется в тех случаях, когда значения каждой из совместно измеренных величин x, y, z, ... не обра-

зуют выборок, но значения функции fi= f (xi, yi, zi, ...) образуют выборку,

т. е. величина f является некоторой физической константой. Штрих у аргументов означает, что они содержат неизвестные постоянные приборные по-

грешности: x

= x

+ θ

( x)i

, y

= y

+ θ

( y)i

, z

= z + θ

( z)i

. Здесь учтено, что

i

i

 

i

i

 

i

i

 

приборные погрешности измеряемых величин могут быть разными в разных опытах, поскольку зависят от отсчетов не образующих выборок величин x, y, z, ... по шкалам приборов.

Статистическая обработка полученных значений fiпроизводится так же, как и в анализе данных прямых измерений, которые позволяют найти ее смещенное среднее значение и СКО среднего значения (либо размах выборки

R = f

f

):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

min

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

N

 

 

 

2

N ( N −1) ,

 

 

 

 

f

=

 

f

,

S

 

=

( f

 

(3.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

f

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N i=1

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

а затем вычислить ее случайную погрешность

f

= tP, N S

 

, или

f = βP, N R .

f

Для определения приборной погрешности θf представим i-е смещенное

значение величины

f

в окрестности точки x , y , z ,

координаты которой не

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

i

 

i

 

i

 

 

 

 

 

зависят от приборных погрешностей,

fi= f (xi′, yi′, zi) = f (xi + θ(x)i , yi + θ( y)i , zi + θ(z)i )

в виде суммы истинного значения этой величины и малого конечного приращения, определяемого выражением (3.1):

f = f

i

+ θ

= f

i

+ a

xi

θ

(x)i

+ a

yi

θ

( y)i

+ a

zi

θ

(z)i

,

(3.8)

i

( f )i

 

 

 

 

 

 

 

 

 

35

где f

i

= f ( x , y , z

) , a

xi

= f

 

 

, a

yi

= f

 

 

, a

zi

= f

 

 

.

 

 

 

i

i i

 

x

 

xi′, yi′, zi

 

y

 

xi′, yi′, zi

 

z

 

xi′, yi′, zi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ввиду

малости

приборных

погрешностей θ( x)i , θ( y)i , θ(z)i

значения

производных в точке xi , yi , zi можно считать совпадающими с их значения-

ми в экспериментальной точке x

¢, y ¢, z ¢. Смещенное среднее значение кос-

 

 

 

 

 

 

 

 

 

i

i

i

 

 

венно определяемой величины с учетом (3.8) будет иметь вид

 

 

 

 

 

 

1

 

N

1

N

 

1

N

 

 

 

 

 

 

¢ =

 

fi¢=

fi

+

q( f )i =

 

 

+ q( f ) ,

(3.9)

 

 

 

f

f

 

 

 

 

 

 

 

 

 

 

 

 

 

N i=1

N i=1

 

N i=1

 

 

где q( f ) =

1

q( f )i =

1

 

(axiq( x)i + ayiq( y)i + aziq( z)i )

– приборная погреш-

 

N

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ность функции.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно

(3.9)

 

несмещенное

значение величины будет

равно

fi = fi¢ - q( f )i , где ввиду неизвестности величин и знаков приборных по-

грешностей θ(x)i , θ( y)i , θ(z)i приборная погрешность функции θ( f ) также неизвестна. Поэтому заменим приборную погрешность функции θ( f ) ее верхней границей q f ³ q( f ) . Тогда

 

 

 

N

 

N

 

q f

=

1

q fi =

1

(

 

axi

 

qxi +

 

ayi

 

qyi +

 

azi

 

qzi ),

(3.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N i=1

N i=1

 

 

 

 

 

 

 

 

 

где θxi, θyi, θzi – верхние границы приборных погрешностей аргументов. Выражение для верхней границы приборной погрешности функции можно также записать в виде, удобном в ряде приложений: q f = axqx + a yqy + azqz , где

 

x =

1

 

axi

 

,

 

 

y =

1

 

ayi

 

,

 

 

z =

1

 

azi

 

; θx = max θxi , θ y = max θ yi ,

 

 

 

 

 

 

 

a

a

a

 

 

 

 

N

N

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θz = max θzi

наибольшие значения верхних границ приборных погрешно-

стей аргументов в серии опытов.

Несмещенное среднее значение функции можно найти как f = f ¢ ± q f .

Тогда результат косвенного измерения с учетом его случайной погрешности можно записать в виде f = f ± Df = f ¢ ± (q f + Df ) = f ¢ ± Df , где Df = q f + Df представляет собой полную погрешность функции.

36

При практических расчетах штрихи у аргументов функции и самой функции опускают.

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

Dfi = fi¢ - f ¢ = ( fi - f ) + (axi - a x )q( x) + (ayi - a y )q( y) + (azi - a z )q( z) ,

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

axi ¹ ax , ayi ¹ ay , azi ¹ az и Dfi ¹ fi - f .

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

(ln f )¢ =

f ¢

f

и вводя новые весовые множители b

 

= a

x

f

i

, b

y

= a

y

f

i

,

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

i

 

 

 

 

 

i

 

i

 

 

 

bz = az

fi , получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где q fi

 

 

 

q fi = fi (

 

bx (xi , yi , zi )

 

qx +

 

by (xi , yi , zi )

 

qy +

 

bz (xi , yi , zi )

 

qz ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

fi

приборная погрешность и значение косвенно определяемой ве-

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

b (x, y, z) =

d ln f (x, y, z)

, b

y

(x, y, z) =

d ln f (x, y, z)

, b (x, y, z) =

d ln f (x, y, z)

.

 

 

 

x

dx

 

 

dy

z

dz

 

 

 

 

 

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

37

3.3. Алгоритм обработки данных косвенных измерений методом переноса погрешностей

Данный метод используется в случае, когда каждая из величин x, y, z , представляющих собой аргументы функции, измеряется независимо от остальных в своей серии опытов, и эти величины образуют выборки (близки друг к другу). Число опытов в сериях, вообще говоря, не обязано быть одинаковым, требуется только неизменность условий для прямо измеряемой величины в своей серии, неизменность условий для f во всех сериях и взаимная независимость всех опытов.

1. По формулам прямых измерений определить величины x , x ; y ,

y ; z , z (с учётом приборных погрешностей).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

 

Рассчитать значение функции

f = f ( x , y ,

 

z ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

 

Вычислить

 

частные

производные

от

 

функции

ax = f

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

ay =

f

 

 

 

 

 

 

 

az = f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x, y, z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или, для легко логарифмируемой функции f, от ее

,

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x, y, z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x, y, z

=

d ln f

 

 

 

=

d ln f

 

 

 

 

 

 

 

 

 

 

=

d ln f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, b

 

 

 

 

 

 

 

 

, b

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

логарифма b

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в точке

x

,

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

z

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x, y, z

 

 

 

 

 

x, y, z

 

 

 

 

 

 

 

 

 

x, y, z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4. По формуле переноса погрешностей вычислить полную погрешность функции f = (ax x)2 + (ay y )2 + (az z )2 или по эквивалентной формуле для легко логарифмируемой функции: f = f (bx x)2 + (by y )2 + (bz z )2 .

5.Записать результат измерения и округлить его.

6.Свести результаты обработки эксперимента в табл. 3.1.

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xi

 

 

 

 

 

 

 

 

θx=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yi

 

 

 

 

 

 

 

 

θy=

x↑i

 

 

 

 

 

 

 

 

 

 

= , Rx=x↑ N x↑1=

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

xi+1 xi

 

 

 

 

 

 

 

 

UP, N Rx =

xi = xi

x

 

 

 

 

 

 

 

 

 

ΣΔxi = 0

 

 

 

 

 

 

 

 

 

 

xi)2

 

 

 

 

 

 

 

 

Σ(Δxi)2=

 

 

 

 

38

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончание табл. 3.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

=

 

 

( xi )2 N (N −1)

=

,

 

 

 

 

 

x= tP,N S

 

 

=

,

 

 

 

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

x2 2x

=

,

 

 

 

x =

 

 

 

±

 

 

=

 

 

 

 

 

, P = %, N =

 

 

 

 

x

 

 

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y↑i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

= , Ry=y↑N y↑1=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yi+1 yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

UP, N Ry =

 

yi = yi

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ΣΔyi =0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yi)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Σ(Δyi)2 =

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

y= tP, N S

 

 

=

 

 

 

 

 

 

 

 

 

S

 

 

 

( yi )2 N (N −1)

,

 

 

 

 

 

 

,

 

 

 

y

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

y2 + θ2y

=

,

 

 

 

 

y =

 

±

 

=

 

 

, P = %, N =

 

 

 

 

y

 

 

 

 

y

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

= f (

 

 

 

,

 

 

 

) =

 

, f =

 

(ax

 

 

)2 + (ay

 

)2 + (az

 

)2

 

f

x

y

 

x

y

z

 

f =

 

±

 

 

 

=

, P = %, N =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В качестве примера обработки данных косвенных измерений методом переноса погрешностей рассмотрим эксперимент по определению ускорения свободного падения g по 5 измерениям периода колебания математического маятника T = 2πl g и его длины l. Выражая g через период колебаний и

длину, получим: g = 4π2lT 2 . Результаты расчетов будут иметь вид табл. 3.2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

li , м

 

 

 

 

0.782

0.810

 

0.795

 

 

0.801

 

 

0.787

θl = 5.10– 4 м

Тi, с

 

 

 

 

1.776

1.798

 

1.789

 

 

1.794

 

 

1.780

θT = 10– 4 c

l↑i

 

 

 

 

0.782

0.787

 

0.795

 

 

0.801

 

 

0.810

Rl = l↑ N – l ↑1 = 0.028,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l = 0.795

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

li+1–l i

 

 

 

 

0.005

 

 

 

0.008

 

 

 

 

 

0.006

 

 

0.009

 

 

UP, N Rl =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0.64. 0.028 = 0.018

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

–0.013

 

0.015

 

0

 

 

 

 

 

0.006

 

 

 

–0.008

 

 

ΣΔli = 0

li = li

– l

 

 

 

 

 

 

li)2

 

 

 

 

169·10– 6

 

 

225·10– 6

 

0

 

 

 

 

36·10– 6

 

 

64·10– 6

 

 

Σ(Δli)2 = 494·10– 6

Sl =

 

 

 

 

 

 

 

 

 

 

= 0.00497 ,

 

 

l= tP, N Sl = 0.013915 ,

 

 

 

 

 

 

(

li )2 N (N −1)

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

l =

 

±

 

 

= 0.795 ± 0.014 м, P = 95 %, N = 5

 

 

 

l 2 + θ2

= 0.013925,

 

 

 

l

l

l

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T↑i

 

 

 

 

1.776

 

1.780

 

1.789

 

 

1.794

 

 

1.798

 

 

 

RT=T↑N –T ↑1=0.022,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T =1.7874

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ti+1–T i

 

 

 

 

0.004

 

0.009

 

 

 

 

 

0.005

 

 

0.004

 

 

 

UP, N RT =0.0141

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

39

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончание табл. 3.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ti= Ti – T

 

–0.0114

0.0106

0.0016

 

0.0066

–0.0074

ΣΔTi=0

Ti)2

 

 

 

 

1.300·10– 4

1.124·10– 4

2.56·10– 6

 

 

4.356·10– 5

5.476·10– 5

Σ(ΔTi)2=3.432·10– 4

 

 

 

=

 

 

 

 

 

 

 

 

= 0.004142, T = tP, N S

 

= 0.011516,

 

S

 

 

(

 

Ti )2

N (N −1)

 

 

T

T

 

 

 

 

=

 

 

 

= 0.011517,

T =

 

±

 

 

= 1.787 ± 0.012 c , P = 95 %, N = 5

 

 

 

 

 

T 2 + θ2

 

 

 

 

 

T

T

T

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

g = g(l,T ) = 4π2 lT 2 = 9.82388.

Для определения погрешности используем метод полного дифференциала.

 

 

 

 

 

 

 

 

 

a =

dg

 

=

2

,

 

 

 

 

a =

dg

= −2

2

l

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dl

 

T

2

 

 

 

 

 

 

 

 

 

dT

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

π

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g = (al l )

+ (aT

T )

 

 

 

l

 

 

 

4

 

l

 

 

 

 

 

l

2 T

 

 

=

 

 

 

 

 

l

+

 

 

 

T

=

 

 

 

 

 

 

 

 

 

+

 

= 0.2136,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

T

3

 

 

 

 

 

 

2

 

 

 

 

 

l

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

g= 9.8 ± 0.2 м/ с2 , P = 95 %, N = 5 .

3.4.Алгоритм обработки данных косвенных измерений

выборочным методом

Выборочный метод применяется в том случае, если совместно измеренные значения аргументов функции xi, yi и zi не образуют выборок, но можно создать выборку значений функции {fi}.

1.По каждому набору совместно измеренных значений аргументов рассчитать значения функции fi = f(xi, yi, zi).

2.Обработать полученную выборку {fi} согласно алгоритму обработки

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

3. Вывести выражения для частных производных от функции

 

ax (x, y, z) =

df (x, y, z)

, ay (x, y, z) =

df (x, y, z)

, az (x, y, z) =

df (x, y, z)

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

dy

 

 

 

 

dz

или для легко логарифмируемой функции f

от ее логарифма

 

 

 

b (x, y, z) =

d ln f (x, y, z)

, b

y

(x, y, z) =

d ln f (x, y, z)

, b (x, y, z) =

d ln f (x, y, z)

.

 

 

 

x

 

 

dx

 

 

 

 

dy

 

 

z

dz

 

 

 

 

 

 

 

 

 

 

 

 

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

40

θ fi = ax (xi , yi , zi ) θxi + ay (xi , yi , zi ) θ yi + az (xi , yi , zi ) θzi ,

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

θ fi = fi ( bx (xi , yi , zi ) θxi + by (xi , yi , zi ) θ yi + bz (xi , yi , zi ) θzi ),

где fi – соответствующее данному набору аргументов значение функции (не путать со строкой таблицы упорядоченных по возрастанию значений f↑i).

5. Если приборные погрешности аргументов одинаковы во всех опытах или при нахождении максимальных по всей серии опытов значений приборных погрешностей θx = max θxi , θ y = max θ yi , θz = max θzi , для определения

приборной погрешности величины f

можно использовать выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ f =

 

 

xθx +

 

yθ y +

 

zθz ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

a

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1

 

 

 

 

 

 

 

 

(x , y , z )

 

,

 

y =

1

 

 

 

 

 

(x , y , z )

 

,

 

 

 

=

1

 

 

 

(x , y , z )

 

.

где ax

 

 

 

 

a

 

 

a

 

 

a

 

 

 

az

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

N

 

 

 

 

 

6.

 

 

 

 

 

x

i

i i

 

 

 

 

 

 

 

y

i i i

 

 

 

 

 

 

 

N

 

z

i

i i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислить

среднюю

приборную

погрешность

функции

 

 

 

1

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ f

=

 

 

θ fi

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

f + θ f .

 

 

 

 

7.

 

Вычислить полную погрешность функции

 

 

f

 

 

 

 

8.Записать результат измерения и округлить его.

9.Свести результаты обработки эксперимента в табл. 3.3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 3.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θxi

 

 

 

 

 

 

 

 

θx = max θxi =

yi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θyi

 

 

 

 

 

 

 

 

θy = max θyi =

fi

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f↑i

 

 

 

 

 

 

 

 

Rf = f↑N f↑1 =

Ufi = fi+1 – f i

 

 

 

 

 

 

 

 

Ufi < UP, N Rf =

fi = fi – f

 

 

 

 

 

 

 

 

ΣΔfi = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fi)2

 

 

 

 

 

 

 

 

Σ(Δfi)2 =

θfi

 

 

 

 

 

 

 

 

θ f = (

 

θ fi

 

) N =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

41