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

книги / Обработка и представление результатов эксперимента

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
2.37 Mб
Скачать

в)

отбрасываем число 58 и с оставшимися 9 числами проводим расчет

j =44,4;

а х = 2,9.

 

7. Задачи

Ниже приведен ряд задач с решениями. Погрешности рассчитываются

методом непосредственных вычислений.

Задача 1. Вычислить выражение

q = (5 ± 1) + (8 ± 2) - (7 ± 1).

Р е ш е н и е

 

 

Наилучшая оценка: q = 5 + 8 - 7 = 6.

 

Абсолютная погрешность:

Ддг=1+2 + 1=4.

Ответ: <7 = 6 ±4.

 

 

Задача 2. Вычислить выражение

q = (5 ± 1) • (8 ±.2).

Р е ш е н и е

 

 

Наилучшая оценка: q = 5 • 8 = 40.

 

Относительная погрешность: zq = 1/5 + 2/8 = 0,2 + 0,25 = 0,45.

Абсолютная погрешность: tsq = 40 • 0,45 = 18 « 20.

Ответ: q = 40 ± 20.

 

 

Задача 3. Вычислить выражение

q = (10 ± 1) / (20 ± 2).

Р е ш е н и е

 

 

НаилУЧшая оценка: q = 10 / 20 = 0,5.

 

Относительная погрешность: eq = 1/10 + 2/20 = 0,1 + 0,1 = 0,2.

Абсолютная погрешность:

Ад = 0,5 • 0,2 = 0,1.

Ответ: q = 0,5 ± 0,1.

 

 

Задача 4. Вычислить выражение q = 2я (10 ± 1). Р е ш е н и е

НаилУЧшая оценка: q - 2я-10 = 62,8 « 63. Абсолютная погрешность: Aq =-2ти*1 = 6,28 « 6. Ответ: q = 63 ± 6.

Задача 5. Посетитель средневекового замка решает определить глубину колодца, измеряя время падения брошенного в него камня. Он определил, что время падения равно / = (3,0 ± 0,5) с. Какова глубина колодца ?

Р е ш е н и е Для вычисления глубины колодца следует воспользоваться кинематиче-

ской формулой Л = g • / /

2.

 

 

 

 

 

Наилучшая оценка:

Л = 9,8 (3,0)2 / 2 -

44,1 м.

 

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

е/, = 2 ^г

= 2 • 0,5 / 3,0

= 0,3.

Абсолютная погрешность: АЛ = 44,1 • 0,3 = 14,6 «15 м.

Ответ: Л = (44 ± 15) м.

 

 

 

 

 

Задача 6. Угол 0 измерен как

(125 ± 2)°

Рассчитайте наилучшее зна­

чение sin 0

и его погрешность.

 

 

 

 

 

Ре ше н и е

 

 

 

 

 

 

 

Наилучшая оценка:

sin 0 = sin 125° = 0,819.

 

Выразим погрешность Д0 в радианах: 2° = 0,0349 рад.

Абсолютная

погрешность:

 

A(sin 0) =

d sin0

А0 = cos 0 • Д0

 

de 0

= |cos(125°)|- 0,0349 = 0,02.

 

 

 

 

 

Ответ:

sin 0

=0,82 ±0,02.

 

 

 

 

 

Задача 7. Вычислить значение

следующего выражения методом «шаг за

шагом»: q = Vl6±4 + (3,0± 0,1)3(2,0 ± 0,1) = (1) + (2)• (3).

 

Здесь пронумерованы (1),

члены выражения.

 

Р е ш е н и е

 

 

 

 

 

 

 

Вычисляем последовательно, записывая все шаги: 1. Вычисляем (1): q\ * Vl6±4 =(16±4)1/2

Наилучшая оценка: л/Гб = 4.

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

^ = 0,5-4/16 =0 125

 

 

 

 

41

2

х

 

 

Абсолютная погрешность: Aqx= 4 • 0,125 = 0,5.

 

Ответ: л/16 ± 4 = 4,0 ± 0,5.

 

 

 

 

 

 

2. Вычисляем (2): q2= (3,0 ± 0,1)3.

 

 

Наилучшая оценка: (3,0)3 = 27,0.

 

 

 

 

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

га. =3 ~

=

3 - 0,1 / 3,0

= 0,1.

 

 

 

 

41

 

х

 

 

Абсолютная погрешность: bq2= 27,0 • 0,1 = 2,7.

 

Ответ: (3,0±0,1)3

=27,0 ±2,7.

 

 

 

 

3. Вычисляем произведение (2) • (3): q3= (27,0 ± 2,7) • (2,0 ± 0,1).

Наилучшая оценка:

27,0 • 2,0 = 54,0.

 

 

 

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

е^3 = 2,7 / 27,0 + 0,1 / 2,0 = 0,15.

Абсолютная погрешность: Aq3= 54,0 • 0,15 = 8,1.

 

Ответ: (27,0 ± 2,7) • (2,0 ± 0,1) = 54,0 ± 8,1.

 

 

4. Вычисляем сумму

qx+ q3: q =(4,0 ± 0,5) + (54,0 ± 8,1).

Наилучшая оценка:

4,0 + 54,0 = 58,0.

 

 

 

Абсолютная погрешность: Aq = 0,5 + 8,1 = 8,6.

 

Окончательный ответ задачи: q = 58,0 ± 8,6 « 58 ± 9.

Задача

8. Покажите,

что

среднее

арифметическое N отклонений

d\, d2, ..., dff

равно нулю.

 

 

 

 

 

 

Р е ш е н и е

 

 

 

 

 

 

 

N i=l

N ,

 

N

i

N

 

 

Задача 9. Чтобы вычислить стандартное отклонение ах для N измерений

х 1, JC2, , xN,

необходимо рассчитать сумму Х(-г/ - ^)

Докажите, что ее

можно представить в виде

£(*,- - х)2 = Х (х/)2 - N - x 2

 

Р е ш е н и е

Z(*, - -О2 = к *2 - 2х ■Xj + х 2) = Y .xf - 2xY.Xj + Nx2 =

= ^ 4 - 2 х Ы х + Ы х 2 = ^ х} - И - х 2

Полученный результат очень полезен практически, так как упрощает расчеты.

Задача 10. После нескольких измерений скорости звука стандартное от­ клонение среднего имело значение ст^ = 10 м/с. Сколько измерений необходимо

сделать, чтобы конечная погрешность результата составляла ± 3 м/с ?

Р е ш е н и е

 

 

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

=GU/ J N

Пусть количество измерений необходи­

мо увеличить в к раз,

тогда

= ои/^ к Отношение погрешностей

Ответ: количество измерений необходимо будет увеличить в 11 раз.

8. Аппроксимация данных методом наименьших квадратов

s.1. Метод наименьших квадратов

До сих пор мы имели дело с многократным измерением одной и той же физической величины и с последующей статистической обработкой получен­ ных данных. Чаще всего в экспериментах требуется совместно измерять значе­ ния двух различных физических переменных х и у для установления в после­

дующем математической связи этих

переменных: у = у(х) (например, для

получения вольт-амперной характеристики

I = I(U) какого-либо элемента).

Наиболее простая и важная связь - линейная

вида у = Л + Вх, которая является

частью более общего полинома типа

у = А + Вх + Сх2 + Dx3 +

Итак, мы измерили N различных значений х и JC2, ..., Хдг и соответствую­

щих им значений у\, у2, ...» ум, то есть получили N точек (xit у \\ (х2, у2\

(*мУи)- Хотя результаты измеренийх н у содержат некоторые погрешности, для

упрощения будем предполагать, что Ах «

Ау, поэтому величинами Ах будем

пренебрегать, кроме того, считаем, что все Ау одинаковы по величине.

 

Задача состоит

в том, чтобы определить

наилучшую

прямую

линию

у = А + Вх, аппроксимирующую набор

экспериментальных

точек

(,гь у \\

(*22\ .

(*ы,Ун\ иными словами, требуется рассчитать коэффициенты А и В.

Зная А

и Ву

мы могли

 

 

 

 

бы для любого значения х, вы­

 

 

 

 

числить

«истинное»

значение

 

 

 

 

= А + Вх( (рис.8.1). Предполага­

 

 

 

 

ем, что результат измерения каж­

 

 

 

 

дого у,

подчиняется

распределе­

 

 

 

 

нию Гаусса с центром на истин­

 

 

 

 

ном значении У, и с шириной

 

 

 

 

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

су,

одинаковой

 

 

 

 

для всех измерений.

 

 

 

 

 

Вероятность получения единичного у,

определяется выражением

 

Ру. ( Л) *

 

шJ _ е- ( * - л-& ,)2/ 2о2 = рА в{у>)

 

1

G y

G y

 

 

 

 

Вероятность получения всего набора N точек

 

 

 

РаМУ1’ - ’Ун ) = раМ у^ ' - 'PAM ? N ) *

е~Р ^2.

 

 

 

 

 

 

Оу

 

 

 

где р2 = 'Ziyi - A - B x i)2/ <т2 .

/=1

В соответствии с принципом максимального правдоподобия, наилучшие оценки А и В такие, при которых вероятность Раув (У\> ->Уи ) будет шах или р2 - min. Отсюда и название метода - метод наименьших квадратов, то есть

сумма квадратов (или р2) должна быть минимальной. Это означает

ЭВ2

2 N

- л - Л / ) - о,

® г

= - 4 Е ( л

дА

Оу /=1

 

откуда определяется величина А. Для определения В вычисляем

ЭВ2

2

^

дВ

<3Ly

/=1

Перепишем два последних уравнения в виде системы двух уравнений для опре­ деления неизвестных А и В:

AN + B ^ X i - ^ y i

A Y xi + B 'Z xf= 'L xiyi .

Решение этой системы:

 

. 5 > ? Е Л - 5 > /5 > /Л

 

А ~

Z

*

3 N 'Z xjyi-'Z xj'Z yi

 

 

Z

 

где

 

 

Z = N Y x f - ( Z x i)2

(8.1)

Можно получить более компактный алгоритм расчета коэффициентов А и

В.Перепишем первое уравнение системы (поделим обе части уравнения на N):

А+ Вх = у ,

где

X = £ * ,/# ; y = 'Zyi/N

Таким образом, наилучшая прямая проходит через точку (Jc,J), то есть через «центр тяжести» всех экспериментальных точек. Выражение для коэффициента

В можно преобразовать:

Q - I Xrfj - x'Zyj

Z(Xj ~ х)у,

I.XjXj-x'EXi

'

Порядок расчета следующий: х, у, В, Л = у - Я г.

8.2.Погрешность в измерениях у

Задача состоит в отыскании выражения для величины ау. В соответствии с принципом максимального правдоподобия, наилучшая оценка для су такая, при которой вероятность Ра,в ъ—*Ум ) получения значений у и у2, ..., yN

будет максимальной. Это означает дРдв jdGy = 0 . Запишем выражение для вероятности в виде

Ра,в (У\,-,Уи ) а ~ Y e ~ a2f 2° y , ау

где

a2 = c2y {i2 = j:(yi - A - B x i)2,

1=1

тогда

дРА>в/да у = [- Na<N+l) +а 2а ^ +3> ] - = 0.

Поскольку нулю может быть равно только выражение в квадратных скобках, то

N =a 2

Gy2,

G2 = а 2 IN ,

О \

N

 

7

— L

(y i- A -B x j)2

N i=1

 

Ha практике наилучшей оценкой является выражение для Gy, в котором

произведена замена N —>N - 2, то есть выражение

9

1

*

9

N ~ 2 i=\

Вспомним, что в выражении стх для стандартного отклонения N измере­ ний одной величины х была произведена замена N -> N - 1. Такие замены можно обосновать следующим образом. Прежде чем вычислить G x, полученные

данные JCI, JC2,

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

чения

x = Y,xi/N

Это оставляло только (ЛМ)

независимых

значений

JCI, JC2,

или

(AM)

степеней свободы.

 

 

В

методе

наименьших квадратов произвели

измерение

N точек

 

(х2ууг),

(JCM yri) и перед расчетом о> вычислили коэффициенты А и В,

то есть установили 2 связи. Сделав это, мы оставили только (N-2) степеней сво­ боды. Таким образом, в формулах для <тх и ау в знаменателях стоит не число произведенных замеров N, а число степеней свободы.

8.3. Погрешности в величинах А и В

Коэффициенты А и В являются функциями х\, Хг, ..., хн и у и У2> уы.

Так как по принятому допущению у величин хи *2, ..., х^ погрешность не учитывается, то необходимо учесть погрешности у А и В простым расчетом ошибок в косвенных измерениях:

2

тг,ОдА

2.

2

 

у

ayi)\2

GA~

сту/)

 

 

7 > 1

 

 

 

i

tyi

После подстановок и преобразований получаем

2 _

2

2

2

 

К

 

 

 

<*В=°у

 

 

где Z определяется выражением (8.1).

 

 

 

 

 

8.4. Аппроксимация полиномом

Часто зависимость у = у (х)

не является линейной, в этом случае для по­

лучения функциональной зависимости у

= у (х) можно использовать полином

у = А + Вх + Сх2 + Dx? +

 

 

 

 

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

квадратичного

полинома у = А + Вх + Сх2, чтобы

определить

значения

коэффициентов Ау

В, С, необходимо минимизировать

выражение

f

= U y i - A - B Xi-C x f)2/ o 2y ,

 

 

 

1= 1

 

 

 

 

 

то есть рассчитать производные

 

 

^

 

= 0;

^

= 0;

3^ = 0

 

а 4

 

дВ

 

дС

 

 

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

коэффициентов А, В, С полинома у = А + Вх + Сх2

 

 

 

 

 

8.5. Экспоненциальные функции

 

В

физике

ряд

зависимостей

между двумя переменными

имеет вид

п

 

где А и В - постоянные, подлежащие определению. В этом случае

у = А е ,

поступают

следующим

образом.

Логарифмируют исходное

выражение

In у = In А + Вх, тем самым сводя его к линейной зависимости z = w + Вх, где

z = In у, w = In А. Находят наилучшие оценки для w и В методом наименьших квадратов, а затем переходят к исходным переменным х, у, А, В.

8.6. Метод наименьших квадратов с весами

Пусть мы измерили N пар значений (х\9у \), (*2, Уг), (*м Улг), которые,

как полагаем, связаны линейной зависимостью у =А+Вх. Для упрощения бу­

дем предполагать, что Ах, « Ау„ поэтому величинами

А*/

будем пренебре­

гать. В отличие от рассмотренного случая

считаем, что

Ау,

разные. Иными

словами, у\ имеет погрешность сгь уг - с ъ

; (ty. Задача состоит в том, что­

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

Вероятность получения всего набора точек

РА,в(Уь-’Уы) = рА,в(УО - - рА,в(Ум),!‘ е р /2, где теперь

р2 = 1(л - A - fix,)2А,2 = Z®,Ov - Л - fix,)2,

1= 1

/=1

соу = \/o f - веса.

В соответствии с принципом максимального правдоподобия минимизи­

руем выражение для р2, то есть рассчитываем выражения

ав2

N

~ г

=”2 Y,®i(yi - A - B x j) = 0 ,

дА

/=1

ав2

N

~ г

= “2 ХсО;(д'/ - A-Bxj)-Xj = 0.

дВ

/=1

Запишем последние выражения в виде системы двух уравнений для опре­

деления коэффициентов полинома А и В:

 

ЛХ>/ + 5 2 > ,х , =

;

Л2> ,х , + fi2>,-x,2 = 1> ,х ,л

Решение этой системы:

 

 

 

1 _ Е СО/Х,21

- £ <о,х,XЮ/Х/У, . д

Z Ю, X^,Х,\У/ - Z СО/Х,X С0,У, ,

 

 

Z

Z

z = Е “ , 1 “ ,х2 - ( 1 со,х,.)2

Метод наименьших квадратов с весами можно использовать, когда по­ грешности ст,, а значит, и веса со, известны.

9. Сглаживание данных

Сглаживание данных эксперимента является специальной операцией ус­ реднения с помощью интерполяционных многочленов. В результате получается уточненная величина р/, при расчете которой использовались близлежащие значения: ...,Л -2,Л - 1,Л ,Л +ь й +2. -

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