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/τ. Three measurements are taken, each with the same uncertainty σ=0.3V:
t (s)
0
1
2
V (V)
9.5
5.6
3.4
Starting from the guess V0=10.0V, τ=1.0s, carry out one Gauss-Newton step.
(a) Write the residual ri=Vi−V0e−ti/τ and the two partial derivatives ∂f/∂V0 and ∂f/∂τ, and evaluate all of these at the starting guess for the three data points.
(b) Assemble the normal equations JTJΔp=JTr, where J is the 3×2 matrix of the derivatives found in (a) (the common weight 1/σ2 cancels because σ is the same for every point, so it can be dropped from both sides).
(c) Solve the 2×2 system for ΔV0 and Δτ, and give the updated parameters V0+ΔV0 and τ+Δτ to three significant figures.