# 仿真建模代写 | COM3001 MODELLING AND SIMULATION OF NATURAL SYSTEMS

本次英国代写主要为仿真建模的限时作业

COM3001 MODELLING AND SIMULATION OF NATURAL SYSTEMS

SECTION A

1. a) Red Foxes and European rabbits can be found in ecosystems all over the UK,

often forming a binary predator-prey relationship within the ecosystems. Based on the

assumptions below, write down a pair of differential equations (do not include code) that

describe the changes over time of fox population F, and rabbit population R. Clearly

define the state variables and parameters, and specify their units.

● Rabbits have access to unlimited food

● The only predator of rabbits are foxes

● Foxes fully depend on rabbits as their food supply, and die otherwise

[5%]

b) Discuss how you can find the equilibrium points of the system? Discuss how you

could check the stability properties of the equilibrium points, i.e. whether they are

asymptotically stable, Lyapunov stable or unstable, via simulation (only discuss, don’t

simulate).

[5%]

c) The system in part (a) was simulated using the Euler method and its trajectory

was plotted on the phase plane in Fig 1(i). This was repeated with a different time step Δt

in Fig 1(ii). Which of the two plots used a larger Δt? Justify your answer.

Fig 1(i) Fig 1(ii)

[5%]

d) Write 3 circumstances for which it would be applicable to use an agent-based

model, rather than an equation-based approach in order to simulate this system?

[5%]

2. a) Write the differential equation (do not include code) and the spiking procedure for

a Leaky Integrate and Fire (LIF) neuron with constant input. Define the state variables

and parameters and provide their units.

[5%]

b) Assume that the constant input remains greater than the threshold for the

following two questions.

i) If threshold increases, will the mean time between two consecutive spikes

increase or decrease and explain why?

ii) If τm (the time constant of leaky integration) increases, will the mean time

between two consecutive spikes increase or decrease and explain why?

[5%]

SECTION B

3. The equations of a predator-prey system are simulated using a 2nd order Runge Kutta

method with time step Δt=0.001 and simulated from time t=0 to t=100 (arbitrary time

units). At t=100 units, the simulated predator and prey values are each off from their ideal

ones, in absolute value by approximately 1e-5. With Δt=0.01, calculate how far off the

predator and prey values are expected to be from their ideal ones, in absolute value, at

t=100 units. Explain your working / calculation.

[15%]

4. Suppose that the volume of a plant bud is proportional to the amount of water it contains.

Take the bud to be spherical. Assume that the rate at which it absorbs water from its

parent stem is proportional to the radius of the bud. Further, assume that the bud loses

water at a rate proportional to its surface area.

Write a differential equation (do not include code) that describes the growth of the volume

of this plant bud over time, as a function of its radius. Explain how you constructed this

equation. What is the differential equation for the rate of change of the radius of the bud?

(Note: This is not necessarily a realistic model of bud growth.)

[15%]

5. Complete the code snippet below, to simulate the time evolution of the logistic growth

system: dN/dt = aN(1-N/Nc) using ode45, starting from t=0 to t=100 (arbitrary time units),

with an absolute tolerance of 1e-5, and a relative tolerance of 1e-6. Use a=0.2, Nc=5000,

N(t=0)=1. You only need to type/write the code on paper, not simulate in Matlab (though

you can, to verify it). Insert your own code at the <your code> locations.

Explain briefly, in your own words, how the ode45 function solves the initial value

problem to maintain the error within tolerance.