EduKG
← Назад

Графическое представление измерений с неопределённостью: моделирование и вычислительное исследование

Цели занятия

  • Итог

    Постройте график табличных данных с усами неопределённости по обеим осям, указав числовые значения uxu_x и uyu_y.

  • Итог

    Определите наклон k±Δkk \pm \Delta k по центральной и крайним допустимым прямым, проходящим через усы.

  • Итог

    Обоснуйте принятие или отклонение линейной модели y=kx+by = kx + b, сравнив её предсказание с коридором допустимых прямых.

Понять идею

Графическое представление измерений с неопределённостью

Текст

Каждое измеренное значение несёт неопределённость uu, заданную разрешением прибора и разбросом повторных измерений. На графике значение изображают точкой, а неопределённость — отрезком («усом») полной длины 2u2u: истинное значение с высокой вероятностью лежит внутри этого интервала. Линейная модель y=kx+by = kx + b считается допустимой, если существует прямая, проходящая одновременно через все интервалы неопределённости.

Наклон такой прямой даёт kk, а крайние допустимые прямые ограничивают коридор, ширина которого равна 2Δk=kmax⁡−kmin⁡2\Delta k = k_{\max} - k_{\min}; половину этой ширины и принимают за неопределённость наклона Δk\Delta k. Точки на краях диапазона определяют kk сильнее всего, поэтому xx планируют так, чтобы охватить возможно более широкий интервал. Если через все усы нельзя провести ни одной прямой, линейная модель отвергается — и это тоже полноценный результат измерения.

Разобранные примеры

Напрямую

Задача

Измерения растяжения пружины дали таблицу:

xx, см 1.0 2.0 3.0
yy, Н 2.2 4.1 5.9

Неопределённости: ux=0.1u_x = 0.1 см, uy=0.2u_y = 0.2 Н. Постройте эти точки с усами и найдите наклон k±Δkk \pm \Delta k.

Решение

Шаг 1 — усы. Вокруг каждой точки рисуем прямоугольник 2ux×2uy2u_x \times 2u_y, то есть 0.20.2 см на 0.40.4 Н, или эквивалентные крестообразные усы.

Шаг 2 — центральная прямая. Через первую и последнюю точки:

kц=5.9−2.23.0−1.0=3.72.0=1.85 Н/см.k_{\text{ц}} = \frac{5.9 - 2.2}{3.0 - 1.0} = \frac{3.7}{2.0} = 1.85 \ \text{Н/см}.

Шаг 3 — крайние допустимые прямые. Самая крутая прямая проходит через низ первой точки и верх последней:

kmax⁡=6.1−2.03.0−1.0=2.05 Н/см.k_{\max} = \frac{6.1 - 2.0}{3.0 - 1.0} = 2.05 \ \text{Н/см}.

Самая пологая — через верх первой и низ последней:

kmin⁡=5.7−2.43.0−1.0=1.65 Н/см.k_{\min} = \frac{5.7 - 2.4}{3.0 - 1.0} = 1.65 \ \text{Н/см}.

Ответ: k=1.85±0.20k = 1.85 \pm 0.20 Н/см, где Δk=(kmax⁡−kmin⁡)/2\Delta k = (k_{\max} - k_{\min})/2.

Применение

Задача

Для тех же данных (xx от 1.0 до 3.0 см, uy=0.2u_y = 0.2 Н) найдите значение yy при x=0x = 0 (свободный член bb) вместе с его неопределённостью. Почему Δb\Delta b больше отдельного uyu_y?

Решение

Центральная прямая имеет kц=1.85k_{\text{ц}} = 1.85 Н/см, поэтому

bц=2.2−1.85×1.0=0.35 Н.b_{\text{ц}} = 2.2 - 1.85 \times 1.0 = 0.35 \ \text{Н}.

Крайние прямые. Крутая прямая: kmax⁡=2.05k_{\max} = 2.05, b=2.0−2.05×1.0=−0.05b = 2.0 - 2.05 \times 1.0 = -0.05 Н. Пологая прямая: kmin⁡=1.65k_{\min} = 1.65, b=2.4−1.65×1.0=0.75b = 2.4 - 1.65 \times 1.0 = 0.75 Н.

Ответ: b=0.35±0.40b = 0.35 \pm 0.40 Н.

Почему Δb\Delta b большое. Экстраполяция к x=0x = 0 удалена от измеренного диапазона на «плечо» x=1.0x = 1.0 см. Ошибка наклона умножается на это плечо: Δb≈Δk⋅xmin⁡=0.20×1.0=0.20\Delta b \approx \Delta k \cdot x_{\min} = 0.20 \times 1.0 = 0.20 Н у каждого конца плюс вклад uy=0.2u_y = 0.2 Н, что и даёт ±0.40\pm 0.40 Н. Чем дальше экстраполяция, тем больше неопределённость.

Перенос

Задача

Модель предсказывает наклон 2.22.2 Н/см. Набор A (данные выше) даёт k=1.85±0.20k = 1.85 \pm 0.20 Н/см. Набор B, полученный более грубой шкалой, даёт k=2.0±0.6k = 2.0 \pm 0.6 Н/см. Для каждого набора решите, согласуется ли модель с данными, и объясните, какой набор информативнее.

Решение

Интервал допустимых наклонов — это и есть усы неопределённости, пересчитанные в наклон.

Набор A: [1.85−0.20; 1.85+0.20]=[1.65; 2.05][1.85 - 0.20;\ 1.85 + 0.20] = [1.65;\ 2.05] Н/см. Предсказание 2.22.2 лежит вне интервала, значит при таком разрешении данные противоречат модели: либо модель неверна, либо есть неучтённая систематическая погрешность.

Набор B: [2.0−0.6; 2.0+0.6]=[1.4; 2.6][2.0 - 0.6;\ 2.0 + 0.6] = [1.4;\ 2.6] Н/см. Предсказание 2.22.2 лежит внутри, значит данные не противоречат модели — но и не подтверждают её: коридор слишком широк.

Вывод. Информативность определяется не числом точек, а шириной коридора. Набор A позволяет опровергнуть модель, набор B — нет; грубые измерения «безопасны», но почти ничего не проверяют.

Практика

Текст

Вычислительное исследование в электронной таблице или в Python — безопасно, никаких нагреваний и опасных веществ не требуется.

  1. На отрезке x=0.5…4.0x = 0.5\ldots 4.0 см сгенерируйте 8 значений по модели y=1.85x+0.35y = 1.85x + 0.35 и добавьте к каждому небольшой случайный шум.
  2. Задайте ux=0.1u_x = 0.1 см, uy=0.2u_y = 0.2 Н и постройте точечный график с крестообразными усами.
  3. Проведите центральную прямую и две крайние допустимые прямые; запишите k±Δkk \pm \Delta k и b±Δbb \pm \Delta b.
  4. Повторите всё, увеличив uyu_y вдвое. Отметьте, во сколько раз выросло Δk\Delta k и изменился ли вывод о совместимости с моделью y=2.2xy = 2.2x.

Сохраните оба графика рядом и сравните ширину коридоров.

Визуальная подсказка

Назначение

Показать, как неопределённость отдельной точки превращается в коридор допустимых прямых, и как из этого коридора читаются Δk\Delta k и интервал для bb.

Представление

Подписанный эскиз-график: оси xx и yy; три точки с прямоугольными усами 2ux×2uy2u_x \times 2u_y; сплошная центральная прямая; две штриховые крайние прямые, образующие заштрихованный коридор; горизонтальные скобки с подписями 2Δk2\Delta k (полная ширина коридора по наклону; Δk\Delta k — её половина) и 2Δb2\Delta b (полная ширина коридора при x=0x = 0).

Примеры из жизни

Спринт, снятый на телефон

Сценарий

Ученик снимает бег одноклассника на дистанции 6060 м телефоном в режиме 240240 кадров в секунду. Метки через каждые 1010 м отмерены рулеткой с точностью ±0,02\pm 0{,}02 м, а момент, когда бегун поравнялся с меткой, виден с точностью одного кадра, то есть ±0,004\pm 0{,}004 с. Получены пары: s=0s = 0 м при t=0t = 0 с; 1010 м при 1,551{,}55 с; 2020 м при 2,952{,}95 с; 3030 м при 4,384{,}38 с; 4040 м при 5,815{,}81 с.

Связь

Каждая пара чисел становится точкой с крестом погрешностей: горизонтальная полоса ±0,004\pm 0{,}004 с, вертикальная ±0,02\pm 0{,}02 м. Прямая наилучшего приближения по участку 2020–4040 м даёт наклон 7,07{,}0 м/с; коридор, который допускают кресты, — от 6,956{,}95 до 7,057{,}05 м/с. Но на первых 2020 м точки ложатся систематически выше прямой и расходятся шире крестов, поэтому реалистичная оценка скорости — (7,0±0,2)(7{,}0 \pm 0{,}2) м/с. График с неопределённостями показывает границу между приборным разбросом и несовершенством модели.

Тормозной путь велосипеда

Сценарий

Велосипедист на сухом асфальте тормозит с разных скоростей: велокомпьютер показывает скорость с точностью ±0,5\pm 0{,}5 км/ч, а тормозной путь отмерен рулеткой с точностью ±0,1\pm 0{,}1 м. Данные: 2020 км/ч — 2,32{,}3 м; 2525 км/ч — 3,53{,}5 м; 3030 км/ч — 5,15{,}1 м; 3535 км/ч — 6,96{,}9 м.

Связь

Зависимость dd от vv — парабола, поэтому данные переносят на график dd от v2v^2, где точки ложатся на прямую. Наклон этой прямой 0,0730{,}073 с2^2/м; кресты погрешностей (по vv — ±0,5\pm 0{,}5 км/ч, что при 3030 км/ч даёт по оси v2v^2 около ±2\pm 2 м2^2/с2^2) расширяют допустимый наклон до 0,0660{,}066–0,0800{,}080 с2^2/м. Из d=v2/(2μg)d = v^2/(2\mu g) наклон равен 1/(2μg)1/(2\mu g), откуда μ=0,70±0,07\mu = 0{,}70 \pm 0{,}07: именно длина крестов на графике задаёт, с какой точностью вообще можно узнать коэффициент трения.

Остывание чая на кухне

Сценарий

Кружку чая при $85{,}0\ ^\circ$C оставляют в комнате, где $22\ ^\circ$C. Термометр имеет цену деления $0{,}5\ ^\circ$C, секундомер — ±1\pm 1 с. Измерения температуры: $85{,}0\ ^\circ$C при t=0t = 0; $72{,}0\ ^\circ$C через 55 мин; $62{,}5\ ^\circ$C через 1010 мин; $55{,}0\ ^\circ$C через 1515 мин; $49{,}0\ ^\circ$C через 2020 мин.

Связь

Разность с комнатой $\Delta T = T - 22\ ^\circ$C убывает как ΔT=ΔT0e−t/τ\Delta T = \Delta T_0 e^{-t/\tau}, поэтому строят график ln⁡(ΔT)\ln(\Delta T) от tt и получают прямую с наклоном −0,042-0{,}042 мин−1^{-1}. Вертикальный крест каждой точки — это цена деления термометра $\pm 0{,}5\ ^\circ$C; при $\Delta T \approx 30\ ^\circ$C это всего ±0,017\pm 0{,}017 в логарифме, и прямые, ещё проходящие через все кресты, дают наклон от −0,038-0{,}038 до −0,046-0{,}046 мин−1^{-1}. Значит, постоянная времени равна 2222–2626 мин, то есть (24±2)(24 \pm 2) мин. Переводя экспоненту в прямую, график делает погрешность наглядной.