
- •1. Введение
- •2. Измерение физических величин
- •Типы величин
- •Типы погрешностей измерений
- •3. Случайные величины и их характеристики
- •4. Нормальное распределение и его свойства
- •Нормальное распределение
- •Правило «3 стандартов»
- •Коэффициент Стьюдента
- •5. Суммарная погрешность измерений
- •6. Погрешности косвенных измерений
- •7. Учет погрешности в записи окончательного результата измерения
- •Порядок выполнения округления
- •8. Линеаризация данных
- •9. Метод наименьших квадратов
- •10. Статистическая проверка гипотез
- •Общие положения
- •Гипотеза совпадения экспериментального среднего и известного значения
- •Гипотеза совпадения двух независимых средних величин
- •Гипотеза о линейности данных
- •11. Основы планирования эксперимента
- •Определение необходимого числа измерений
- •Ведение лабораторного журнала
- •Требования к оформлению научного отчета
- •12. Контрольные задания и вопросы

Порядок выполнения округления
Выполнить предварительную запись окончательного результата измерения в виде x = x ± x и вынести за общую
скобку одинаковые порядки среднего и погрешности, т.е. множитель вида 10k, где k – целое число. Числа в скобках переписать в десятичном виде с использованием запятой, убрав тем самым оставшиеся порядковые множители.
Округлить в скобках число, соответствующее погрешности, до одной значащей (ненулевой) цифры слева, если эта цифра больше 2, или до двух первых цифр в противном случае. При округлении используют правило: если цифра, расположенная за оставляемой, меньше 5, то ее просто отбрасывают, иначе оставляемую цифру увеличивают на единицу. Если же отбрасываемая цифра равна 5, то наименьшая ошибка достигается при округлении по правилу Гаусса до ближайшего четного числа. К примеру, 4,5 округляют до 4, в то время как 3,5 также округляют до 4.
Округлить в скобках число, соответствующее среднему значению: последними справа оставляют цифры тех разрядов, которые сохранились в погрешности после ее округления.
Окончательно записать x = x ± x с учетом выполнен-
ных округлений. Общий порядок и единицы измерения величины приводят за скобками – получена стандартная форма записи.
8. Линеаризация данных
Физические величины, определяющие результаты эксперимента, выступают в роли переменных и параметров некоторой функциональной зависимости, теоретически получаемой в рамках модели. После экспериментальной регистрации зависимости ее сравнивают с теоретической. Методом сравнения можно не только численно определить, т.е. измерить,
21

значения физических величин, не измеряемых другим способом, но и вывести заключение об адекватности применения модели к эксперименту.
Проще всего проверить линейную зависимость:
y = ax +b, |
(8.1) |
где x, y – измеряемые величины, a, b – параметры зависимости. Если зависимость нелинейная, в некоторых случаях ее можно преобразовать в линейную (табл. 4).
Таблица №4
Способы линеаризации зависимостей
Вид нели- |
Получаемая ли- |
|
|
|
|
|||
|
|
|
|
|||||
нейной |
нейная зависи- |
y |
x |
a |
b |
|||
зависимо- |
мость |
|
|
|
|
|||
|
сти |
|
|
|
|
|||
|
|
|
|
|
|
|||
v = k u z |
ln v = z ln u + ln k |
ln v |
ln u |
z |
ln k |
|||
v = k ezu |
ln v = zu + ln k |
ln v |
u |
z |
ln k |
|||
v = k ez u |
ln v = zu −1 +ln k |
ln v |
u–1 |
z |
ln k |
|||
v = |
|
u |
|
v−1 = ku −1 + z |
v–1 |
u–1 |
k |
z |
|
k + zu |
|||||||
|
|
|
|
|
|
|
Обращаем внимание на то, что после вычисления погрешностей величин a и b переход к погрешностям реальных физических величин k и z осуществляется по формулам для погрешностей косвенных измерений.
9. Метод наименьших квадратов
Этот метод является одним из наиболее распространенных приемов статистической обработки экспериментальных данных, относящихся к различным функциональным зависимостям физических величин друг от друга. В том числе он применим к линейной зависимости и позволяет получить
22

достоверные оценки ее параметров a и b, а также оценить их погрешности.
Рассмотрим статистическую модель эксперимента, в котором исследуют линейную зависимость. Пусть проведено n > 2 парных измерений величин x и y: xi, yi, где i = 1, ..., n. По экспериментальным данным необходимо найти оценки
параметров a и b, а также оценки их дисперсий σa2 и σb2 . О
природе экспериментальных погрешностей сделаем следующие предположения:
1.Значения xi известны точно, т.е. без погрешностей.
2.Распределения величин yi взаимно независимы, имеют
одну и ту же дисперсию σ2 и отвечают нормальному закону. Распределения yi имеют средние значения yi , которые сов-
падают с точным значением функции axi + b.
3. Систематические погрешности отсутствуют. В частности, все промахи, т.е. точки, выходящие за интервал 3σ, отброшены.
Запишем функцию правдоподобия как вероятность реализации набора полученных экспериментальных данных:
L = |
1 |
|
|
|
|
|
(y1 −(ax1 |
+b))2 |
|
|
|
|
|
1 |
|
(yn −(axn |
+b))2 |
|||||||
|
|
|
exp − |
|
|
|
|
2σ2 |
|
|
... |
|
|
|
|
|
exp − |
2σ2 |
= |
|||||
σ 2π |
|
|
|
|
|
σ |
|
|
2π |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
n |
|
|
|
1 |
|
n |
|
|
|
|
|
|
|
2 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
= |
|
|
|
exp |
− |
|
|
|
∑(yi −(axi +b)) |
|
. |
|
|
(9.1) |
||||||||||
σ 2π |
|
2σ |
2 |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
Логарифмируем обе части: |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
n |
|
|
|
n |
|
|
|
1 |
|
|
n |
|
|
|
||
|
|
ln L = − |
ln 2π − |
lnσ2 |
− |
|
|
∑(yi −(axi +b))2 . |
|
|||||||||||||||
|
|
2 |
|
2σ |
2 |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
i=1 |
|
|
|
Оценками a, b, σ2 будет правильным считать значения, при которых L и lnL максимальны, т.е. реализуется наибольшая вероятность получения набора экспериментальных данных. Экстремум функции lnL находят дифференцированием:
23

∂ln L |
= 0 , |
∂ln L |
= 0 , |
∂ln L |
= 0. |
(9.2) |
∂a |
|
∂b |
|
∂σ2 |
|
|
После дифференцирования система уравнений относительно искомых параметров примет вид
n |
|
|
|
|
|
|
|
∑[xi (yi −axi −b)]= 0 , |
|
|
|||||
i=1 |
|
|
|
|
|
|
|
n |
|
−axi −b)= 0 , |
|
|
|||
∑(yi |
|
(9.3) |
|||||
i=1 |
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
nσ 2 = ∑(yi −axi −b)2 . |
|
|
|||||
|
i=1 |
|
|
|
|
|
|
|
|
z = |
1 |
n |
|
|
|
Вводя средние величины |
∑zi , получаем |
|
|||||
|
|
|
n i=1 |
|
|
||
a = |
xy |
− x |
|
y , |
|
|
|
|
x2 |
− x 2 |
|
|
|||
b = y −a x |
= |
x2 |
y |
− x xy |
, |
(9.4) |
|
x2 |
− x 2 |
||||||
|
|
|
|
σ 2 = y2 − y 2 −a2 ( x2 − x 2 ).
Для получения несмещенной относительно точного значения оценки дисперсии (вследствие конечного числа измерений полученные значения параметров a и b отличаются от предельных) последнее выражение надо умножить на n(n −2):
σ 2 = |
n |
y2 |
− |
y 2 −a2 ( x2 − |
x |
2 ) . |
(9.5) |
|
n −2 |
|
|
|
|
|
|
Оценим теперь дисперсии параметров. Преобразуем выражение для a:
24

n |
|
xj − x |
|
|
|
a = ∑k j y j , где k j |
= |
|
|
. |
(9.6) |
n |
) |
||||
|
|
∑(xj − x |
|
|
|
j=1 |
|
|
2 |
|
|
|
|
|
|
|
j=1
После преобразования видно, что a получается как линейная комбинация взаимно независимых величин yj, так как коэффициенты kj заданы точно – согласно пункту 1 предположений о статистике изучаемых величин. Следовательно, параметр a распределен нормально, а его дисперсия σa2 представляет собой линейную комбинацию дисперсий величин yj с коэффициентами kj2 – это свойство сложения нормальных распределений уже встречалось при рассмотрении погрешностей косвенных измерений:
2 |
n |
2 |
2 |
|
|
nσ 2 |
|
|
|
|
|
σ 2 |
|
|
|
|
σa |
= ∑k j σ |
|
= |
|
|
|
|
2 |
= |
|
|
|
|
|
(9.7) |
|
|
j=1 |
|
|
n |
2 |
|
n |
|
|
|
n |
( |
x2 − |
x |
) |
|
|
|
|
|
n∑xi |
− |
∑xi |
|
|
|
|
|
|
||||
|
|
|
|
i=1 |
|
i=1 |
|
|
|
|
|
|
|
|
|
Подставляя выражение для дисперсии из формулы (9.5), получаем оценку для стандартного отклонения параметра a:
σ |
a |
= |
1 |
y2 |
− |
y 2 |
−a2 . |
(9.8) |
|
|
n −2 |
x2 |
− |
x 2 |
|
|
Преобразуем выражение для b:
|
1 |
n |
1 |
n |
|
|
b = |
∑yi −a |
∑xi . |
(9.9) |
|||
|
|
|||||
|
n i=1 |
n i=1 |
|
Параметр b также нормально распределен. Его стандартное отклонение:
σb = |
σ 2 |
+ |
x |
2 |
σa2 =σa |
x2 . |
(9.10) |
n |
|
||||||
|
что для перехода от величин σa и |
||||||
Обращаем внимание, |
|||||||
σb к погрешностям a и |
b их следует умножить на коэф- |
||||||
фициент Стьюдента: |
|
|
|
|
|
|
|
25

a = t(α, n −1)σa , |
(9.11) |
b = t(α, n −1)σb , |
|
где α – уровень значимости, n – количество парных измерений (см. табл. 2). Значение n – 1 берется в связи с тем, что в методе наименьших квадратов из экспериментальных данных находят не одну величину, а две – a и b. Связь между ними уменьшает количество независимых случайных переменных, складывающихся в распределение Стьюдента.
Если линейная зависимость была получена из модельной путем линеаризации, переход к погрешностям реальных физических величин осуществляется по формулам для погрешностей косвенных измерений, после чего идет окончательная запись результата.
Описанный выше способ определения наилучшей линейной аппроксимации данных носит название метода наименьших квадратов, поскольку при данных значениях параметров достигается минимум величины «отклонения» прямой от экспериментальных данных:
n |
|
S = ∑(yi −(axi +b))2 . |
(9.12) |
i=1
Частным является случай, когда предполагается, что теоретическая линейная зависимость проходит через начало координат, т.е. имеет функциональный вид y = kx . Для оп-
ределения величины параметра k и его погрешности используют формулы:
|
k = |
xy |
, |
|
|
|
x2 |
|
|
σ 2 = n −n |
2 ( y2 |
−k 2 x2 ), |
(9.13) |
26