A capacitor discharges through a resistor

PHYSICS · Subjective · Fermi University Physics · Scientific Practice & Experimental Methods · Mathematical Methods for Physics · Data analysis and model fitting in the laboratory

Computational Physics, Ch Finding the best parameters: minimisation and fitting

Problem

A capacitor discharges through a resistor, and its voltage is modelled as V(t)=V0e−t/τV(t) = V_0 e^{-t/\tau}. Three measurements are taken, each with the same uncertainty σ=0.3 V\sigma = 0.3\,\mathrm{V}:
tt (s)012
VV (V)9.55.63.4
Starting from the guess V0=10.0 VV_0 = 10.0\,\mathrm{V}, τ=1.0 s\tau = 1.0\,\mathrm{s}, carry out one Gauss-Newton step. (a) Write the residual ri=Vi−V0e−ti/τr_i = V_i - V_0 e^{-t_i/\tau} and the two partial derivatives ∂f/∂V0\partial f/\partial V_0 and ∂f/∂τ\partial f/\partial\tau, and evaluate all of these at the starting guess for the three data points. (b) Assemble the normal equations JTJ Δp=JTrJ^{\mathsf T}J\,\Delta\mathbf{p} = J^{\mathsf T}\mathbf{r}, where JJ is the 3×23\times2 matrix of the derivatives found in (a) (the common weight 1/σ21/\sigma^2 cancels because σ\sigma is the same for every point, so it can be dropped from both sides). (c) Solve the 2×22\times2 system for ΔV0\Delta V_0 and Δτ\Delta\tau, and give the updated parameters V0+ΔV0V_0 + \Delta V_0 and τ+Δτ\tau + \Delta\tau to three significant figures.

Browse the problem library