A Mathematical Model
A Mathematical Model
One division between models is based on the type of outcome they predict.
Deterministic models ignore random variation, and so always predict the
same outcome from a given starting point. On the other hand, the model may
be more statistical in nature and so may predict the distribution of possible
outcomes. Such models are said to be stochastic.1
A model may help to explain a system and to study the effects of different
components, and to make predictions about behavior
Empirical Mechanistic
Deterministic Predicting cattle growth Planetary motion, from a
regression relationship
based on Newtonian
mechanics with feed
intake(differential equations)
Stochastic Analysis of variance of Genetics of small
variety yields over sites populations based on
and years Mendelian inheritance,
(probabilistic equations)
One further type of model, the system model, is worthy of mention. This is
built from a series of sub-models, each of which describes the essence of
some interacting components. The above method of classification then refers
more properly to the sub-models: different types of sub-models may be used
in any one system model. Much of the modelling literature refers to ’simulation
models’. Why are they not included in the classification? The reason for this
apparent omission is that ’simulation’ refers to the way the model calculations
are done - i.e. by computer simulation. The actual model of the system is not
changed by the way in which the necessary mathematics is performed,
although our interpretation of the model may depend on the numerical
accuracy of any approximations . In the physical sciences, a traditional
mathematical model contains most of the following elements:
1. Governing equations
2. Supplementary sub-models
1. Defining equations
2. Constitutive equations
3. Assumptions and constraints
1. Initial and boundary conditions
2. Classical constraints and kinematic equations
Dp/dt=ap
where p(t)is population size at time -t, and a is a constant. Solution of this
equation by integration gives
p(t) =p(0)ea
where p(0) is population size at time zero. According to this solution,
populations grow in size at an exponential rate. Clearly, not all populations
grow exponentially fast. Since the differential equation arose from an
interpretation of the assumption, we must look to the assumption for an
explanation for this discrepancy. In this case, the explanation is the qualifier”
in the absence of limiting factors”. Most natural populations are subject to
constraints such as food supply or habitat which restrict the range of
sustainable population sizes. It is important that all assumptions are stated
clearly and concisely. This allows us to return to them later to assess their
appropriateness.
Where no information exists about the form of a relationship, the only way forwards
is to acquire a body of data and fit equations to it. This has the advantage that we are in
control of the analysis.
2.4.1 Analytically
There is much to be gained from obtaining an analytical solution to a model. This will
allow us to perform all of the manipulations implied by the model with the minimum of
fuss. Note that a full analytical solution for a stochastic model involves finding the
distribution of outcomes, but we may feel satisfied if we can solve the equations for the
mean and standard deviation.
2.4.2 Numerically
When analytical methods are unproductive we can use numerical methods to obtain
approximate solutions. Although they can never have the same generality as analytical
solutions, they can be just as good in any particular instance.
Numerical solution of model equations generally mimics the processes described in
the model. For difference equations, numerical solution is exact since we can use the
rules laid down in the equations to follow the evolution of the system. With a stochastic
model, we can repeatedly simulate out comes using a random number generator as
described earlier, and combine a large number of simulations to approximate the
distribution of outcomes. Differential equations provide a rather more difficult problem.
The basic method is to divide continuous time into discrete intervals, and to estimate the
state of the system at the start of each interval. Thus the approximate solution changes
through a series of steps. The crudest method for calculating the steps is to multiply the
step length by the derivative at the start of the interval. This is called Euler’s method.
More sophisticated techniques are used in performing the Runge-Kutte types of
integration. Fourth order Runge-Kutte is both commonly used and sufficiently accurate
for most applications. It is always worth treating numerical solutions to differential
equations with caution
Philosophical considerations