EduKG
← Back

Graphing measurements with uncertainty

Learning objectives

  • Outcome

    Plot tabular data with uncertainty bars on both axes, labelling the numerical values of uxu_x and uyu_y.

  • Outcome

    Determine the slope k±Δkk \pm \Delta k from the central and the extreme admissible straight lines drawn through the bars.

  • Outcome

    Justify accepting or rejecting a linear model y=kx+by = kx + b by comparing its prediction with the corridor of admissible straight lines.

Understand the idea

Graphing Measurements with Uncertainty

Text

Every measured value carries an uncertainty uu set by the instrument resolution and the scatter of repeated readings. On a graph the value is drawn as a point and the uncertainty as a bar ('whisker') of total length 2u2u: the true value most likely lies inside that interval. A linear model y=kx+by = kx + b is admissible if a single straight line can pass through all the uncertainty intervals at once.

The slope of that line gives kk, and the extreme admissible lines bound a corridor whose width is 2Δk=kmax⁡−kmin⁡2\Delta k = k_{\max} - k_{\min}; half of that width is reported as the slope uncertainty Δk\Delta k. Points at the ends of the range constrain kk most strongly, so the xx values should be spread as widely as possible. If no straight line can be drawn through all the bars, the linear model is rejected — and that too is a genuine measurement result.

Worked examples

Direct

Problem

Measurements of a spring extension gave the table:

xx, cm 1.0 2.0 3.0
yy, N 2.2 4.1 5.9

Uncertainties: ux=0.1u_x = 0.1 cm, uy=0.2u_y = 0.2 N. Plot these points with bars and find the slope k±Δkk \pm \Delta k.

Solution

Step 1 — bars. Around each point draw a 2ux×2uy2u_x \times 2u_y rectangle, that is 0.20.2 cm by 0.40.4 N, or the equivalent cross-shaped bars.

Step 2 — central line. Through the first and last points:

kc=5.9−2.23.0−1.0=3.72.0=1.85 N/cm.k_{\text{c}} = \frac{5.9 - 2.2}{3.0 - 1.0} = \frac{3.7}{2.0} = 1.85 \ \text{N/cm}.

Step 3 — extreme admissible lines. The steepest line passes through the bottom of the first point and the top of the last:

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

The shallowest passes through the top of the first point and the bottom of the last:

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

Answer: k=1.85±0.20k = 1.85 \pm 0.20 N/cm, where Δk=(kmax⁡−kmin⁡)/2\Delta k = (k_{\max} - k_{\min})/2.

Application

Problem

For the same data (xx from 1.0 to 3.0 cm, uy=0.2u_y = 0.2 N) find the value of yy at x=0x = 0 (the intercept bb) together with its uncertainty. Why is Δb\Delta b larger than a single uyu_y?

Solution

Central line: kc=1.85k_{\text{c}} = 1.85 N/cm, hence

bc=2.2−1.85×1.0=0.35 N.b_{\text{c}} = 2.2 - 1.85 \times 1.0 = 0.35 \ \text{N}.

Extreme lines. Steep line: kmax⁡=2.05k_{\max} = 2.05, b=2.0−2.05×1.0=−0.05b = 2.0 - 2.05 \times 1.0 = -0.05 N. Shallow line: kmin⁡=1.65k_{\min} = 1.65, b=2.4−1.65×1.0=0.75b = 2.4 - 1.65 \times 1.0 = 0.75 N.

Answer: b=0.35±0.40b = 0.35 \pm 0.40 N.

Why Δb\Delta b is large. The extrapolation to x=0x = 0 lies a lever arm of x=1.0x = 1.0 cm away from the measured range. The slope uncertainty is multiplied by that arm: Δb≈Δk⋅xmin⁡=0.20×1.0=0.20\Delta b \approx \Delta k \cdot x_{\min} = 0.20 \times 1.0 = 0.20 N at each end, plus the contribution uy=0.2u_y = 0.2 N, which gives ±0.40\pm 0.40 N. The farther the extrapolation, the larger the uncertainty.

Transfer

Problem

A model predicts a slope of 2.22.2 N/cm. Dataset A (the data above) gives k=1.85±0.20k = 1.85 \pm 0.20 N/cm. Dataset B, taken with a coarser scale, gives k=2.0±0.6k = 2.0 \pm 0.6 N/cm. For each dataset decide whether the model is compatible with the data, and explain which dataset is more informative.

Solution

The admissible interval for the slope is just the uncertainty bars converted into slope.

Dataset 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] N/cm. The prediction 2.22.2 lies outside the interval, so at this resolution the data contradict the model: either the model is wrong or an unaccounted systematic error is present.

Dataset 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] N/cm. The prediction 2.22.2 lies inside, so the data do not contradict the model — but they do not confirm it either: the corridor is too wide.

Conclusion. Informativeness is set not by the number of points but by the width of the corridor. Dataset A can falsify the model, dataset B cannot; coarse measurements are 'safe' but test almost nothing.

Practice

Text

A computational study in a spreadsheet or in Python — safe, no heating and no hazardous substances are required.

  1. On the interval x=0.5…4.0x = 0.5\ldots 4.0 cm generate 8 values from the model y=1.85x+0.35y = 1.85x + 0.35 and add a small random noise to each.
  2. Set ux=0.1u_x = 0.1 cm, uy=0.2u_y = 0.2 N and make a scatter plot with cross-shaped bars.
  3. Draw the central line and the two extreme admissible lines; record k±Δkk \pm \Delta k and b±Δbb \pm \Delta b.
  4. Repeat everything with uyu_y doubled. Note by what factor Δk\Delta k grew and whether the verdict on compatibility with the model y=2.2xy = 2.2x changed.

Save both graphs side by side and compare the widths of the corridors.

Visual reference

Purpose

Show how the uncertainty of a single point turns into a corridor of admissible lines, and how Δk\Delta k and the interval for bb are read off from that corridor.

Representation

A labelled sketch-graph: xx and yy axes; three points with rectangular 2ux×2uy2u_x \times 2u_y bars; a solid central line; two dashed extreme lines forming a shaded corridor; horizontal brackets labelled 2Δk2\Delta k (full corridor width in slope; Δk\Delta k is half of it) and 2Δb2\Delta b (full corridor width at x=0x = 0).

Real-world examples

A sprint filmed with a phone

Scenario

A student films a classmate running 6060 m, using a phone in slow-motion mode at 240240 frames per second. Marks every 1010 m are measured with a tape to ±0,02\pm 0{,}02 m, and the moment the runner passes a mark can be read to one frame, that is to ±0,004\pm 0{,}004 s. The data are: s=0s = 0 m at t=0t = 0 s; 1010 m at 1,551{,}55 s; 2020 m at 2,952{,}95 s; 3030 m at 4,384{,}38 s; 4040 m at 5,815{,}81 s.

Connection

Every pair of numbers becomes a point with an error cross: a horizontal bar of ±0,004\pm 0{,}004 s and a vertical bar of ±0,02\pm 0{,}02 m. A best-fit straight line over the 2020–4040 m section has a slope of 7,07{,}0 m/s; the corridor allowed by the crosses is 6,956{,}95 to 7,057{,}05 m/s. Over the first 2020 m, however, the points lie systematically above the line and spread wider than the crosses, so a realistic speed estimate is (7,0±0,2)(7{,}0 \pm 0{,}2) m/s. A graph with uncertainties marks the boundary between instrument scatter and an imperfect model.

Bicycle braking distance

Scenario

A cyclist brakes on dry asphalt from several different speeds: the bike computer reads speed to ±0,5\pm 0{,}5 km/h, and the braking distance is measured with a tape to ±0,1\pm 0{,}1 m. The data are: 2020 km/h — 2,32{,}3 m; 2525 km/h — 3,53{,}5 m; 3030 km/h — 5,15{,}1 m; 3535 km/h — 6,96{,}9 m.

Connection

The dependence of dd on vv is a parabola, so the data are moved to a graph of dd against v2v^2, where the points fall on a straight line. Its slope is 0,0730{,}073 s$^2$/m; the error crosses (a ±0,5\pm 0{,}5 km/h uncertainty in vv, which at 3030 km/h is about ±2\pm 2 m$^2$/s$^2$ along the v2v^2 axis) widen the admissible slope to 0,0660{,}066–0,0800{,}080 s$^2$/m. Since d=v2/(2μg)d = v^2/(2\mu g) makes the slope equal to 1/(2μg)1/(2\mu g), this gives μ=0,70±0,07\mu = 0{,}70 \pm 0{,}07: the size of the crosses on the graph fixes how accurately the friction coefficient can be known at all.

Tea cooling in the kitchen

Scenario

A mug of tea at $85{,}0\ ^\circ$C is left in a room at $22\ ^\circ$C. The thermometer has a scale division of $0{,}5\ ^\circ$C and the stopwatch reads to ±1\pm 1 s. The recorded temperatures are: $85{,}0\ ^\circ$C at t=0t = 0; $72{,}0\ ^\circ$C after 55 min; $62{,}5\ ^\circ$C after 1010 min; $55{,}0\ ^\circ$C after 1515 min; $49{,}0\ ^\circ$C after 2020 min.

Connection

The excess over room temperature, $\Delta T = T - 22\ ^\circ$C, falls as ΔT=ΔT0e−t/τ\Delta T = \Delta T_0 e^{-t/\tau}, so one plots ln⁡(ΔT)\ln(\Delta T) against tt and gets a straight line with slope −0,042-0{,}042 min$^{-1}$. The vertical cross of each point is the thermometer division $\pm 0{,}5\ ^\circ$C; at $\Delta T \approx 30\ ^\circ$C that is only ±0,017\pm 0{,}017 in the logarithm, and the lines still passing through all crosses give slopes from −0,038-0{,}038 to −0,046-0{,}046 min$^{-1}$. Hence the time constant is 2222–2626 min, that is (24±2)(24 \pm 2) min. Turning the exponential into a straight line makes the uncertainty visible.