Биномиальный закон действует, когда число объектов, из которого проводят статистические выборки, ограничено, и число возможных значений результата равно двум (например, вылет частицы в переднюю или заднюю полусферу). При этом нарушается одно из условий применимости закона Пуассона, так как число событий в данном интервале будет зависеть от числа событий в предыдущих интервалах.
Параметрами биномиального распределения являются вероятность одного события p и полное число возможных событий N0.
Вероятность осуществления N событий
|
P |
= C N |
P N (1 − p) N0 −N . |
(П.2.3) |
||
|
N |
N0 |
|
|
|
|
Здесь C N |
– число сочетаний из N0 по N: |
|
||||
N |
0 |
|
|
|
|
|
|
|
C N |
= |
N0 ! |
. |
|
|
|
|
|
|||
|
|
N0 |
|
N !(N0 − N )! |
|
|
|
|
|
|
|
||
Распределение имеет такой же вид, как и общий член разложения бинома Ньютона. При N0 → ∞ биномиальное распределение переходит в распределение Пуассона. Среднее значение N для биномиального распределения равно N0 p, дисперсия D = = N0 p (1 – p). Вероятность одного события может быть как постоянной величиной (например, при бросании игральной кости), так и переменной (для многих задач ядерной физики). Так, вероятность распада радиоактивного ядра за время t, как известно, имеет вид
p = 1 – exp(–λt),
где λ – константа радиоактивного распада. Вероятность распада N ядер из N0 ядер за время t определяется соотношением (П.2.3). Такой же вид, как вероятность распада радиоактивного ядра, имеет вероятность регистрации γ-кванта высокой энергии в пузырьковой камере
p = 1 – exp(–σ n r),
где r – расстояние от точки образования γ-кванта до точки конверсии; σ – сечение конверсии на атом; n – число атомов в единице объема. Пусть в одном акте образуется N0 γ-квантов (например,
115
исследуется распад K20 → 3π0 → 6γ ) при средней эффективности
регистрации γ-квантов p. Тогда вероятность наблюдения в одном акте N квантов можно рассчитать по формуле (П.2.3).
В случае радиоактивного распада среднее число не распавшихся ядер к моменту t равно
N = N0 exp(−λt) ,
где N0 – полное число ядер в момент времени t = 0. Дисперсия в случае биномиального закона меньше, чем в случае закона Пуассона, так как N не может быть больше N0 при биномиальном распределении, а при пуассоновском значения N не ограничены.
Следует отметить, что рассмотренные статистические законы дают распределение вероятности при известных параметрах. В экспериментальной практике чаще приходится решать обратную задачу – например, по наблюдаемым на опыте функциям распределения находить константу распада исследуемого вещества. Можно показать, что в случае биномиального распределения дисперсия параметра p равна p (1 – p) N0.
χ2-Распределение. Если хi – независимые нормально распределенные случайные величины со средним значением, равным нулю, и средним квадратичным отклонением, равным единице, то сумма квадратов этих величин
n
m2 = ∑xi2 i=1
подчиняется χ2-распределению с n степенями свободы:
|
1 |
|
|
χ2 |
|
|
P(χ2 ) = |
|
(χ2 )n / 2−1 exp |
− |
|
. |
(П.2.4) |
2n / 2 Γ(n / 2) |
|
|||||
|
|
|
2 |
|
|
Здесь Γ(n/2) – гамма-функция. Среднее значение χ2 = n , дисперсия
D(χ2) = 2 n.
χ2-Распределение используется при рассмотрении согласия между теоретическими расчетами и результатами эксперимента, при аппроксимации экспериментальных данных аналитическими функциями.
116
П.2.2. Расчет среднеквадратичных погрешностей
Критерием точности выполняемых измерений является среднеквадратичная погрешность, или стандартное отклонение σ, равное положительному значению квадратного корня из дисперсии. Отношение погрешности к истинному значению измеряемой величины (если последняя не равна нулю) называется относительной погрешностью. Относительная погрешность лучше характеризует достоверность результата, чем абсолютная.
Необходимо различать прямые и непрямые измерения, когда значение исследуемой величины вычисляется на основании результатов измерений других величин. Погрешность непрямых измерений вычисляется с помощью соотношений, связывающих искомую физическую величину с непосредственно измеряемыми величинами. Эти соотношения могут выражать известные физические законы, в некоторых случаях они должны быть установлены на основании полученных экспериментальных данных. Пусть искомая величин Z связана с несколькими непосредственно измеренными вели-
чинами х1, ..., хm функцией Z = f(х1, ..., хm).
Разлагая Z в ряд Тейлора в окрестности истинных значений Z0 и x0i, составляя соотношение Z – Z0, возводя его в квадрат и усредняя по распределениям величин х0i, можно получить при пренебрежении членами второго порядка, что
|
|
n |
|
∂ f |
|
2 |
|
m |
∂f |
|
∂f |
|
|
(Z − Z0 )2 |
|
|
|
|
|||||||||
= σ2 |
= ∑ |
|
σ2x |
+ |
∑ |
|
ρijσxi σx j , |
||||||
|
|
|
|||||||||||
|
|
i=1 |
|
∂xi |
|
i |
i, j=1 |
∂xi ∂x j |
|||||
где коэффициент корреляции
ρij = (xi − x0i ) (x j − x0 j ) [σ(xi )σ(x j )]−1.
Если погрешности измерений хi не коррелированы друг с другом, то ρij = 0, тогда стандартное отклонение непрямых измерений будет иметь хорошо известный вид:
m |
|
∂f |
|
2 |
|
|
|
σZ = ∑ |
|
σ2x |
. |
(П.2.5) |
|||
|
|||||||
i=1 |
|
∂xi |
|
i |
|
||
|
|
|
|
|
|
||
|
117 |
|
|
|
|
||
Практическое следствие этого соотношения: для создания оптимальных условий основные усилия должны быть направлены не на дальнейшее уточнение тех результатов измерений, которые и так являются наиболее точными, а на совершенствование наименее точных измерений.
В табл. П.2.1 приведены соотношения между среднеквадратичными погрешностями для некоторых функций.
|
|
|
Таблица П.2.1 |
|
|
||
Функция |
Соотношения между среднеквадратичными |
||
|
|
погрешностями |
|
Z = A ± B |
σ2Z = σ2A + σ2B |
|
|
Z = A×B, |
2 |
2 |
2 |
Z = A/B |
(σZ / Z ) |
= (σA / A) |
+ (σB / B) |
Z = An |
σZ / Z = nσ A / A |
|
|
Z = eA |
σZ/Z = σA |
|
|
Z = lnA |
σZ = σA/A |
|
|
П.2.3. Учет фона
Практически всегда при работе со счетными физическими приборами приходится исключать фон. Фон находится как разность двух измерений: суммарного N в течение времени t (полезный эффект плюс фон) и фонового N в течение времени tф (без источника излучения, при введенном в прибор экранирующем фильтре и пр.).
Исследуемый эффект представляет собой разность результатов двух измерений, отнесенных к одинаковым интервалам времени, например к единице времени:
|
A = |
N |
− |
N |
ф |
= n − nф, |
(П.2.6) |
|
t |
tф |
|||||
|
|
|
|
|
|||
σA = σ2n |
+ σn2 = |
|
N /t 2 + Nф /tф2 = |
n/t + nф /tф . (П.2.7) |
|||
|
ф |
|
|
|
|
|
|
Наиболее удобно выразить фон в единицах времени измерения суммарного эффекта, например если времена t и tф кратны, то эффект
118
A = N − tфt Nф ,
его погрешность
σA =
N + (t / tф)2 Nф .
Предполагается, что время измеряется с точностью много лучшей, чем число импульсов. Поэтому величина t не имеет дисперсии.
Следует заметить, что соотношение σN = σN =
N применимо
только к величинам, непосредственно измеренным на опыте. Погрешность скорости счета, например, зависит от того, каким способом получена эта величина, и может быть, в принципе, сделана как угодно малой при увеличении времени измерения.
П.2.4. Рациональный выбор времени измерения
К этой задаче можно подойти с двух сторон. Во-первых, можно найти то наименьшее время всех измерений, которое необходимо для получения заданной относительной погрешности δ = σА/A окончательной расчетной величины.
Во-вторых, при заданном общем времени, которое отводится для проведения всех измерений, можно найти распределение времени между измерениями, дающее наименьшую относительную погрешность исследуемой величины А.
В реальных экспериментах обычно ставят пробные опыты, в которых проверяется работа отдельных элементов установки, определяются интервалы значений каждой из величин и оцениваются их возможные погрешности. Последнее оказывает непосредственное влияние на проведение всего эксперимента. Большое внимание следует уделять измерению тех величин, погрешности которых вносят основной вклад в погрешность конечного результата. Поэтому при проведении эксперимента следует, априори, расчетным путем оптимизировать времена каждого из измерений с точки зрения их допустимых погрешностей, по возможности провести предварительные измерения, а затем составить план с указанием величин, которые необходимо измерить, и времени, отводимого на каждое измерение.
119