Pure 2.11 - Functions in modelling

Pure 2.11 - Functions in modelling

A function model links a real input to a predicted output. Learn to interpret its constants, choose a suitable family, fit it to data and decide where its predictions can be trusted. A useful refinement changes a particular weakness in the model.

What a function model says

Can a cooling drink end up colder than the room in a model that is supposed to describe ordinary cooling? A straight-line prediction can say yes. The algebra may be correct while the model needs changing; by the end, you will be able to explain and repair that mismatch.

First consider T(t)=18+62e−0.08tT(t)=18+62e^{-0.08t}, where TT is temperature in degrees Celsius and tt is minutes since pouring. At the start it predicts 18+62=8018+62=80, while much later it approaches 1818.

You already know

Substitution evaluates a function, e0=1e^0=1, and a negative exponential e−kte^{-kt}, k>0k>0, approaches zero as tt increases. A horizontal translation or vertical shift changes a graph's position.

The equation alone is incomplete: its variables, units, allowed inputs and assumptions tell you what it describes.

Function model

A function model is a simplified rule predicting one quantity from another. Variables change within the model; parameters are constants controlling its behaviour. The contextual domain contains the input values for which the situation and assumptions make sense.

The cooling rule is mathematically defined for all real tt, but negative time is outside “minutes since pouring”. Even t≥0t\ge0 need not guarantee reliability: the room temperature might change later.

Worked example 1: interpret initial and limiting values

Question. For T(t)=18+62e−0.08tT(t)=18+62e^{-0.08t}, 0≤t≤600\le t\le60, find the initial temperature and interpret 1818.

Plan. Substitute zero, then consider what happens as the exponential becomes small.

T(0)=18+62e0initial input=80use e0=1e−0.08t→0 as t→∞exponential decayT(t)→18limiting value of the rule\begin{aligned} T(0)&=18+62e^0 &&\text{initial input}\\ &=80 &&\text{use }e^0=1\\ e^{-0.08t}&\to0\text{ as }t\to\infty &&\text{exponential decay}\\ T(t)&\to18 &&\text{limiting value of the rule} \end{aligned}

Why distinguish the two constants? The added 1818 is the limiting temperature; the 6262 is the initial excess above it.

Answer. The initial temperature is 8080 degrees Celsius; 1818 represents room temperature. The limiting behaviour describes the mathematical rule, beyond its stated 6060-minute use interval.

Check. The predicted starting temperature exceeds the room temperature by 6262 degrees Celsius.

Change only the initial excess from 6262 to 5252. The surrounding temperature and decay parameter stay fixed.

Worked example 2: a different initial excess

Question. Interpret the initial and limiting temperatures in S(t)=18+52e−0.08tS(t)=18+52e^{-0.08t}, 0≤t≤600\le t\le60.

Plan. Apply the same zero-input and long-term checks.

S(0)=18+52(1)=70initial temperatureS(t)−18=52e−0.08t>0positive for finite tS(t)→18exponential term tends to zero\begin{aligned} S(0)&=18+52(1)=70 &&\text{initial temperature}\\ S(t)-18&=52e^{-0.08t}>0 &&\text{positive for finite }t\\ S(t)&\to18 &&\text{exponential term tends to zero} \end{aligned}

Why is the curve above room temperature? The excess stays positive even as it becomes smaller.

Answer. Initially 7070 degrees Celsius, tending towards 1818 degrees Celsius under unchanged conditions.

Check. At t=0t=0, the excess 70−18=5270-18=52 matches the coefficient.

A pleasant mathematical surprise is that these two drinks have the same proportional cooling of their excess: after 1010 minutes each excess is multiplied by e−0.8≈0.449e^{-0.8}\approx0.449. The total temperatures are not multiplied by that factor, because the room-temperature offset stays.

Watch out

The wrong statement “the initial temperature is 1818” ignores the exponential term at zero. Always calculate T(0)T(0) before interpreting a parameter.

Change the initial excess once more for this short interpretation.

Choosing a function family

Look for the pattern before selecting a rule. Repetition suggests a trigonometric model; constant percentage change suggests an exponential; a constant product suggests inverse proportion.

PatternFamily and formMeaning to check
tide or annual daylight cycley=M+Acos⁡(2π(t−h)P)y=M+A\cos(\frac{2\pi(t-h)}P)midline MM, amplitude $
growth or decay by a fixed factor in equal intervalsy=abty=ab^t, a>0,b>0a>0,b>0, or aektae^{kt}initial value and growth/decay rate
inverse proportiony=k/xy=k/x, x>0x>0 in many applicationsconstant product xy=kxy=k

For the cosine form shown, take A>0A>0 when hh is the time of a maximum. All angles here are in radians. The midline, amplitude and period formulas used for these models are not in the formula booklet: learn them.

[DIAGRAM: asset_name: Three function model families; asset_slug: p1_2_11_use_functions_in_modelling_and_appreciate_their_limitations_and_refinements__diagram_01; caption: Repetition, constant percentage decay and inverse proportion; description: Three graphs show the constructed tide D=7+3cos(pi t/6), exponential C=20exp(-0.3t) and gas model P=240/V, with their important features labelled.]
Repetition, constant percentage decay and inverse proportion

These are constructed examples, rather than measurements from a named harbour or experiment. For a fixed amount of gas at constant temperature, the simple inverse model P=240VP=\frac{240}V predicts 80 kPa80\,\mathrm{kPa} at V=3V=3 litres and 100 kPa100\,\mathrm{kPa} at V=2.4V=2.4 litres. Both products PVPV equal 240240; the physical assumptions are essential.

Worked example 3: fit a repeating cycle

Question. A harbour depth is modelled with maximum 10 m10\,\mathrm m at midnight, minimum 4 m4\,\mathrm m six hours later and period 1212 hours. Find a cosine model for depth DD metres, tt hours after midnight, 0≤t≤120\le t\le12.

Plan. Average the extremes for the midline, halve their difference for the amplitude and use cosine at its maximum.

M=(10+4)/2=7midlineA=(10−4)/2=3amplitudeω=2π/12=π/6one cycle per 12 hours\begin{aligned} M&=(10+4)/2=7 &&\text{midline}\\ A&=(10-4)/2=3 &&\text{amplitude}\\ \omega&=2\pi/12=\pi/6 &&\text{one cycle per }12\text{ hours} \end{aligned}

Why use cosine without a shift? It equals 11 at the starting maximum.

Answer. D(t)=7+3cos⁡(πt6)D(t)=7+3\cos(\frac{\pi t}6), 0≤t≤120\le t\le12.

Check. D(0)=10D(0)=10, D(6)=4D(6)=4, D(12)=10D(12)=10.

Change only the time of high tide, delaying the same cycle by two hours. The heights and period remain fixed.

Worked example 4: move the cycle in time

Question. With the same extremes and period, high tide occurs at t=2t=2, low tide at t=8t=8. Find a model on 2≤t≤142\le t\le14.

Plan. Keep the midline and amplitude, replacing tt by t−2t-2.

M=7,A=3,ω=π/6unchanged shapeD(t)=7+3cos⁡(π(t−2)6)delay by two hours\begin{aligned} M&=7,\quad A=3,\quad\omega=\pi/6 &&\text{unchanged shape}\\ D(t)&=7+3\cos\left(\frac{\pi(t-2)}6\right) &&\text{delay by two hours} \end{aligned}

Why subtract inside? At t=2t=2, the cosine must receive zero to put the peak there.

Answer. D(t)=7+3cos⁡(π(t−2)6)D(t)=7+3\cos(\frac{\pi(t-2)}6), 2≤t≤142\le t\le14.

Check. D(2)=10D(2)=10, D(8)=4D(8)=4, D(14)=10D(14)=10.

Daylight can be modelled similarly: extremes 1616 and 88 hours, with period 1212 months and a maximum at t=0t=0, give 12+4cos⁡(πt6)12+4\cos(\frac{\pi t}6). A fixed sinusoid assumes identical cycles and ignores irregular weather effects on tides or gradual changes in the underlying pattern.

Watch out

A constant percentage increase is not a constant amount increase. For instance, multiplying by 1.11.1 each interval is exponential; adding 0.10.1 each interval is linear.

Change only the high-tide timing again.

Fitting constants and using a model

A proposed form becomes a usable model when observations determine its constants. Keep exact values during fitting, then round the prediction. For V=pektV=pe^{kt}, the parameter pp is the initial value; k>0k>0 gives growth and k<0k<0 gives decay.

You already know

Natural logarithms reverse exponentials: eu=ve^u=v gives u=ln⁡vu=\ln v, for v>0v>0. Divide by a multiplier before taking logarithms; 15000=24000e4k15000=24000e^{4k} first becomes e4k=15000/24000e^{4k}=15000/24000.

This inverse operation turns an unknown exponent into an ordinary linear equation.

Fit an exponential from two observations

If V(0)=V0>0V(0)=V_0>0 and V(T)=VT>0V(T)=V_T>0, T>0T>0, in the model V=pektV=pe^{kt}, then

p=V0,k=1Tln⁡(VTV0).p=V_0,\qquad k=\frac1T\ln\left(\frac{V_T}{V_0}\right).

This follows by dividing VT=V0ekTV_T=V_0e^{kT} by V0V_0 and taking logarithms.

Not in the formula booklet: learn it.

You can equivalently write V=V0(VTV0)t/TV=V_0(\frac{V_T}{V_0})^{t/T}; it makes the observed time interval particularly clear.

Worked example 5: fit depreciation

Question. A machine's value VV pounds, tt years after purchase, follows V=pektV=pe^{kt}. It costs £24,000\pounds24{,}000 and is worth £15,000\pounds15{,}000 after four years. Find the model and predict its value after six years.

Plan. Use the initial value for pp, then solve the four-year equation for kk.

p=24000set t=015000=24000e4kuse the second observationk=14ln⁡(5/8)divide and take logsV(6)=24000(5/8)6/4predict with exact constants=11858.541…evaluate\begin{aligned} p&=24000 &&\text{set }t=0\\ 15000&=24000e^{4k} &&\text{use the second observation}\\ k&=\tfrac14\ln(5/8) &&\text{divide and take logs}\\ V(6)&=24000(5/8)^{6/4} &&\text{predict with exact constants}\\ &=11858.541\ldots &&\text{evaluate} \end{aligned}

Why avoid rounding kk first? The rounded rate would be exponentiated, reducing the accuracy of the final value.

Answer. V=24000(58)t/4V=24000(\frac58)^{t/4}, predicting about £11,900\pounds11{,}900 after six years, to 33 significant figures.

Check. V(4)=24000⋅58=15000V(4)=24000\cdot\frac58=15000; the six-year value is smaller, consistent with decay.

Change only the observed four-year value, from £15,000\pounds15{,}000 to £12,000\pounds12{,}000. The model now loses half its value in each four-year interval.

Worked example 6: a changed decay rate

Question. Fit the same model with initial value £24,000\pounds24{,}000 and four-year value £12,000\pounds12{,}000, then predict the six-year value.

Plan. Replace the observed ratio, retaining the same fitting method.

p=24000initial observation12000=24000e4knew four-year observationk=14ln⁡(1/2)solve for the rateV(6)=24000(1/2)3/2=8485.281…predict at six years\begin{aligned} p&=24000 &&\text{initial observation}\\ 12000&=24000e^{4k} &&\text{new four-year observation}\\ k&=\tfrac14\ln(1/2) &&\text{solve for the rate}\\ V(6)&=24000(1/2)^{3/2}=8485.281\ldots &&\text{predict at six years} \end{aligned}

Why is the new prediction lower? The observed retention factor 12\frac12 is smaller than 58\frac58.

Answer. V=24000(12)t/4V=24000(\frac12)^{t/4}, predicting £8,490\pounds8{,}490 to 33 significant figures.

Check. The model gives £12,000\pounds12{,}000 at four years and £6,000\pounds6{,}000 at eight years, bracketing the six-year value.

The same fitting works for mass: an initial 80 mg80\,\mathrm{mg} falling to 50 mg50\,\mathrm{mg} in three hours gives M=80(58)t/3M=80(\frac58)^{t/3}, predicting 36.6 mg36.6\,\mathrm{mg} after five hours. Context changes the interpretation, not the algebra.

Exam tip

Use the stated model and give the prediction in context. A modelled resale value is an estimate, even when the calculation is exact.

Change the four-year observation once more for a matched fitting task.

Limitations and validity

A model can fit its observations and still fail later. Check the contextual domain, physically possible outputs and assumptions before trusting a prediction. For a battery model B=100e−0.12tB=100e^{-0.12t}, the initial output is 100%100\%, but extending to negative time predicts more than a full charge and lies outside “hours after charging”.

Interpolation and extrapolation

Interpolation predicts between observed input values. Extrapolation predicts outside their range. For example, a six-year prediction fitted from ages zero and four is extrapolation, even if it is mathematically defined.

Extrapolation is often less reliable because the assumptions may stop holding. A possible output is only a basic check, not proof that the model is accurate.

Worked example 7: find a physical contradiction

Question. Algae area AA, in square metres, follows A=3e0.4tA=3e^{0.4t}, where tt is days. The lake surface is 200200 square metres. Predict the area at t=12t=12 and assess the model there.

Plan. Calculate the prediction and compare it with the available surface.

A(12)=3e4.8substitute the time=364.531…evaluate364.531…>200compare with the capacity\begin{aligned} A(12)&=3e^{4.8} &&\text{substitute the time}\\ &=364.531\ldots &&\text{evaluate}\\ 364.531\ldots&>200 &&\text{compare with the capacity} \end{aligned}

Why use the lake area? It supplies a real upper bound that the exponential model lacks.

Answer. The model predicts about 365365 square metres, which is impossible on a 200200-square-metre surface, so it is unsuitable at this time.

Check. The prediction is more than the entire lake, not merely a small rounding excess.

Change only the prediction time to ten days. A smaller prediction removes that particular contradiction, without establishing reliability.

Worked example 8: distinguish possible from validated

Question. Use the same algae model at t=10t=10. Does the lake-area check establish that the model is accurate?

Plan. Evaluate and compare, then state what this check can actually show.

A(10)=3e4substitute=163.794…evaluate163.794…<200physically possible area\begin{aligned} A(10)&=3e^4 &&\text{substitute}\\ &=163.794\ldots &&\text{evaluate}\\ 163.794\ldots&<200 &&\text{physically possible area} \end{aligned}

Why avoid concluding “accurate”? Passing one bound says nothing about whether the growth rate matches observations.

Answer. About 164164 square metres is possible, but accuracy still requires suitable observations and reasonable assumptions.

Check. 0<A(10)<2000<A(10)<200; this is the check passed, rather than a validation of the entire model.

Likewise, P=240VP=\frac{240}V needs V>0V>0, constant temperature and a fixed quantity of gas. Heating the gas changes the pressure relationship; very small volumes may expose the limits of the idealised inverse-proportion rule.

Key idea

Name the assumption or bound that fails, then connect it to the prediction. “The model is unrealistic” on its own does not explain anything.

Keep the same bound and change the prediction time once more.

Refining a model

A refinement changes a specific weakness. You might restrict the time interval, fit new parameters from more observations, add a limiting value, use different rules over different intervals or include an ignored variable. Choose the change that addresses the evidence.

Consider a drink initially at 9090 degrees Celsius in a room at 2020 degrees Celsius. A line can describe the first few minutes, but its constant decrease eventually carries it through room temperature. A shifted exponential can build in the limiting temperature.

Worked example 9: add a sensible limiting value

Question. Compare L(t)=90−2.6tL(t)=90-2.6t and E(t)=20+70e−0.06tE(t)=20+70e^{-0.06t} for ordinary cooling, with tt in minutes. Calculate both at t=30t=30 and explain the better long-term model.

Plan. Compare predictions with room temperature and inspect their limiting behaviour.

L(30)=90−2.6(30)=12linear predictionE(30)=20+70e−1.8=31.570…exponential predictionE(t)→20excess tends to zero\begin{aligned} L(30)&=90-2.6(30)=12 &&\text{linear prediction}\\ E(30)&=20+70e^{-1.8}=31.570\ldots &&\text{exponential prediction}\\ E(t)&\to20 &&\text{excess tends to zero} \end{aligned}

Why add 2020? A model 90e−0.06t90e^{-0.06t} would tend to zero rather than to the room temperature.

Answer. LL predicts 1212 degrees Celsius, below the room temperature; EE predicts 31.631.6 degrees Celsius and tends towards 2020. Under fixed room conditions, EE has the more suitable long-term behaviour.

Check. E(0)=20+70=90E(0)=20+70=90; the initial temperature is preserved.

Change only the room temperature to 2525 degrees Celsius, keeping the initial temperature and assumed decay parameter. This changes the limiting value and therefore the initial excess.

Worked example 10: refine the surrounding conditions

Question. Construct the corresponding cooling model for a drink initially at 9090 degrees Celsius in a room at 2525 degrees Celsius, keeping decay parameter 0.060.06. Predict its temperature after 3030 minutes.

Plan. Use 2525 as the limit and choose the excess to preserve the initial value.

90−25=65initial excessR(t)=25+65e−0.06trefined modelR(30)=25+65e−1.8=35.744…predict\begin{aligned} 90-25&=65 &&\text{initial excess}\\ R(t)&=25+65e^{-0.06t} &&\text{refined model}\\ R(30)&=25+65e^{-1.8}=35.744\ldots &&\text{predict} \end{aligned}

Why change the coefficient too? Keeping 7070 would give an incorrect initial temperature of 9595.

Answer. R(t)=25+65e−0.06tR(t)=25+65e^{-0.06t}, predicting 35.735.7 degrees Celsius after 3030 minutes.

Check. R(0)=90R(0)=90, and R(t)>25R(t)>25 for every finite t≥0t\ge0.

Holding the decay parameter fixed was an explicit modelling assumption, not a claim that changing room conditions never affects it. New temperature readings could estimate it again. For a daily cycle between 33 and 99, the midline is 66, amplitude 33; if the daily maximum changes with season, fitting separate seasonal models or a varying amplitude would address that weakness.

The opening question is now resolved. The drink can fall below room temperature in the linear equation because the equation assumes a constant cooling amount forever. A limiting-value refinement removes that contradiction under the stated physical assumptions.

Problem-solving

Connect each refinement to its purpose: a new limit fixes long-term behaviour; new observations improve fitted constants; a restricted domain avoids an unsupported prediction.

Change the surrounding temperature once more while preserving the initial condition.

Quick Recap

  • Define variables, units, parameters and the contextual domain before using a model.
  • Match repeated cycles, percentage change and constant products to suitable function families.
  • Fit constants from observations and keep exact values until the final prediction.
  • Check domain, physical bounds and assumptions; a possible output is not proof of accuracy.
  • Refine the specific feature that causes a failure.

The opening drink falls below room temperature only because the linear rule keeps subtracting a fixed amount. A shifted exponential preserves the initial temperature and approaches the room temperature instead.