книги / Обработка и представление результатов эксперимента
..pdfв) |
отбрасываем число 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\
(*мУи)- Хотя результаты измеренийх н у содержат некоторые погрешности, для
упрощения будем предполагать, что Ах « |
Ау, поэтому величинами Ах будем |
||||||
пренебрегать, кроме того, считаем, что все Ау одинаковы по величине. |
|
||||||
Задача состоит |
в том, чтобы определить |
наилучшую |
прямую |
линию |
|||
у = А + Вх, аппроксимирующую набор |
экспериментальных |
точек |
(,гь у \\ |
||||
(*2,У2\ . |
(*ы,Ун\ иными словами, требуется рассчитать коэффициенты А и В. |
||||||
Зная А |
и Ву |
мы могли |
|
|
|
|
|
бы для любого значения х, вы |
|
|
|
|
|||
числить |
«истинное» |
значение |
|
|
|
|
|
= А + Вх( (рис.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. -