Pure 9.5 - Numerical methods in context

Pure 9.5 - Numerical methods in context

Turn a model into a numerical target, then interpret the result with its units and limits. You will find a model's maximum by iteration, estimate totals from rate readings and separate calculation accuracy from confidence in the model.

From context to numerical target

When does a vehicle reach its greatest speed between two stops? For the model v=(12−0.5t)ln⁡(t+1)v=(12-0.5t)\ln(t+1), with tt in seconds and 0≤t≤240\le t\le24, the answer needs more than a calculator root: we must choose the right equation and explain what its root means.

At t=0t=0 and t=24t=24, the model gives zero speed. At t=7t=7, it gives about 17.7 m s−117.7\,\mathrm{m\,s^{-1}}. The greatest speed occurs somewhere between the stops. We will find its time shortly.

You already know

A threshold asks for a root of “model minus target”. An interior maximum or minimum is a candidate where the derivative is zero. Integrating a rate gives an accumulated quantity, with units of rate multiplied by time.

Choose the mathematical target before choosing a numerical method.

Context asks forMathematical targetAppropriate numerical approach
Time a value is reachedSolve F(t)=0F(t)=0Root-finding or a suggested iteration
Greatest or least valueSolve the derivative equation and check candidatesIteration or Newton-Raphson if needed
Total from rate readingsEstimate an integralTrapezium rule
Proof of rounded rootTrap it inside rounding boundariesSign change and continuity

The same model can lead to different targets. Here is a pressure example where differentiation would answer the wrong question.

Worked example 1: a threshold target

Question. A cabin pressure is modelled by P=11+0.08tcos⁡(2t)P=11+0.08t\cos(2t) psi, where tt is hours and 0≤t≤90\le t\le9. Form a numerical target for the first time P=10.5P=10.5.

Plan. Equate the model to the required pressure, then move the target to the left.

11+0.08tcos⁡(2t)=10.5threshold equationF(t)=0.5+0.08tcos⁡(2t)=0root-finding form\begin{aligned} 11+0.08t\cos(2t)&=10.5 && \text{threshold equation} \\ F(t)=0.5+0.08t\cos(2t)&=0 && \text{root-finding form} \end{aligned}

Why this step? A time of reaching a level is a root problem; a stationary pressure would instead involve a derivative.

Answer. Find the smallest root of F(t)=0F(t)=0 in [0,9][0,9], using radians. Report the time in hours after take-off.

Check. F(0)=0.5F(0)=0.5, so take-off is not the requested threshold time. A root outside the stated interval would not answer the model's question.

Change only the threshold from 10.510.5 to 10.610.6 psi. The model and variable remain the same.

Worked example 2: a changed threshold

Question. Form the corresponding numerical target for the first time P=10.6P=10.6.

Plan. Subtract the new target from the model.

11+0.08tcos⁡(2t)=10.6G(t)=0.4+0.08tcos⁡(2t)=0new residual\begin{aligned} 11+0.08t\cos(2t)&=10.6 \\ G(t)=0.4+0.08t\cos(2t)&=0 && \text{new residual} \end{aligned}

Why this step? The residual must vanish at the target pressure, rather than at the old threshold.

Answer. Seek the smallest root of G(t)=0G(t)=0 in [0,9][0,9], again using radians and hours.

Check. At a root, the rearranged equation gives 11−0.4=10.611-0.4=10.6, as required.

A numerical routine may find a later crossing. Use a sketch or small-interval search to locate the first relevant crossing before refining it.

Watch out

Wrong: a converged root is automatically the requested time. Right: check the domain, the variable's meaning and whether the question asks for the first crossing.

Now recognise a different target from its units.

Iteration for contextual equations

Return to the vehicle. To locate its greatest speed, differentiate the model first. The resulting equation mixes a logarithm and an algebraic expression, so a suggested iteration is useful.

[DIAGRAM: asset_slug: p1_9_5_use_numerical_methods_to_solve_problems_in_context__diagram_01; asset_name: Maximum speed found by iteration; caption: The model's greatest speed occurs about 7.09 seconds after departure; description: Exact speed v=(12-0.5t)ln(t+1) on 0 to 24 seconds, with stationary endpoints and the unique maximum at t approximately 7.08922. A separate small table gives the first five iterates from t1=7 without covering the curve.]
The model's greatest speed occurs about 7.09 seconds after departure

Worked example 3: the time of maximum speed

Question. Find by iteration when v=(12−0.5t)ln⁡(t+1)v=(12-0.5t)\ln(t+1), 0≤t≤240\le t\le24, is greatest. Start with t1=7t_1=7.

Plan. Set the derivative to zero, rearrange as a fixed-point equation and iterate.

dvdt=−0.5ln⁡(t+1)+12−0.5tt+1product rule24−t=(t+1)ln⁡(t+1)set derivative to zero25=(t+1)[1+ln⁡(t+1)]add t+1tn+1=251+ln⁡(tn+1)−1isolate t\begin{aligned} \frac{\mathrm{d}v}{\mathrm{d}t}&=-0.5\ln(t+1)+\frac{12-0.5t}{t+1} && \text{product rule} \\ 24-t&=(t+1)\ln(t+1) && \text{set derivative to zero} \\ 25&=(t+1)[1+\ln(t+1)] && \text{add }t+1 \\ t_{n+1}&=\frac{25}{1+\ln(t_n+1)}-1 && \text{isolate }t \end{aligned}

Starting at t1=7t_1=7, the next values are 7.118355118…7.118355118\ldots, 7.079822017…7.079822017\ldots, 7.092265173…7.092265173\ldots, 7.088236347…7.088236347\ldots, 7.089539676…7.089539676\ldots, eventually approaching 7.089221031…7.089221031\ldots.

Why this step? Differentiation finds the stationary-time equation; iterating the speed formula itself would solve a different problem.

Answer. The model predicts greatest speed after 7.09 s7.09\,\mathrm s, to three significant figures.

Check. Throughout [0,24][0,24],

d2vdt2=−1t+1−12−0.5t(t+1)2<0\frac{\mathrm{d}^2v}{\mathrm{d}t^2}=-\frac1{t+1}-\frac{12-0.5t}{(t+1)^2}<0

so the derivative decreases through zero once. Together with zero endpoint speeds and positive interior speeds, this establishes the unique global maximum.

Change only the initial coefficient 1212 to 1010. The stop time becomes 2020 seconds to keep the same physical interpretation.

Worked example 4: changing the vehicle model

Question. For v=(10−0.5t)ln⁡(t+1)v=(10-0.5t)\ln(t+1), 0≤t≤200\le t\le20, derive a similar iteration and estimate the maximum-speed time. Start with t1=6t_1=6.

Plan. Apply the same derivative and rearrangement to the changed coefficient.

0=−0.5ln⁡(t+1)+10−0.5tt+120−t=(t+1)ln⁡(t+1)tn+1=211+ln⁡(tn+1)−1\begin{aligned} 0&=-0.5\ln(t+1)+\frac{10-0.5t}{t+1} \\ 20-t&=(t+1)\ln(t+1) \\ t_{n+1}&=\frac{21}{1+\ln(t_n+1)}-1 \end{aligned}

The iterates are 6.128526987…6.128526987\ldots, 6.084770174…6.084770174\ldots, 6.099517644…6.099517644\ldots, 6.094530254…6.094530254\ldots, eventually 6.095789440…6.095789440\ldots.

Why this step? The changed model alters the numerator of the recurrence; the earlier value 2525 must become 2121.

Answer. The greatest speed occurs after approximately 6.10 s6.10\,\mathrm s, to three significant figures.

Check. The analogous second derivative is negative on [0,20][0,20], and the speed is zero at both ends, so the stationary point is the unique maximum.

The numerical answer is a prediction under the stated speed model. Its digits do not establish that a real vehicle follows the model perfectly.

Exam tip

Give the requested iterate to the requested accuracy, then state a contextual limiting answer with units. Three decimal places and three significant figures are different instructions.

Use a supplied recurrence with a fresh coefficient.

Trapezium rule in context

A rate table asks for accumulation rather than a root. Multiply rate units by the time units before calculating, converting units first when needed. Speed in miles per hour and time in hours produce distance in miles.

[DIAGRAM: asset_slug: p1_9_5_use_numerical_methods_to_solve_problems_in_context__diagram_02; asset_name: Estimating distance from speed readings; caption: The three trapezia estimate the distance travelled during 0.3 hours; description: Recorded speeds of 48,54,50,46 miles per hour at times 0,0.1,0.2,0.3 hours are joined by straight segments. Shaded trapezia show distance as the area under the speed-time graph.]
The three trapezia estimate the distance travelled during 0.3 hours

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account