代做SEMT20001-R: Principles of Computational Modelling Coursework Reassessment: Question 4代写Python编程

SEMT20001-R: Principles of Computational Modelling

COURSEWORK REASSESSMENT: QUESTION 4

Question 4:    Equation-based modelling

(a)    (i)  Write Python code that numerically solves an ODE system

(1)

given a right-hand-side function f, timespan [t0, tf], initial condition y0, and number of steps N / step-size h = (tf - t0)/N, using the time-stepping algorithm

(2)

where

(2 marks)

(ii) Using both the forward Euler method, and the method defined in part (a)(i), solve the ODE initial value problem

(3)

for t ∈ [1, 3], with step-size h = 0.5.  Plot a graph of the numerical solutions, together with the exact solution

(4)

and comment on the results.           (4  marks)

(iii)  Plot a graph of the absolute error of the numerical solution of equation (3), using the method given in part (a)(i), versus step-size h, for a suitably-chosen range of h.  You should use the Euclidean distance between the numerical solution

(5)

to define the error. Use the graph to find the order of the truncation error as h → 0.    (4  marks)

(b) The ODE system dt/dy = f(t, y), written in terms of components y = (S, E, I, R),

represents an equation-based model of the transmission of a non-fatal disease with four states:

in which the population is divided into four groups:

•  Susceptible: healthy and not yet infected with the disease,

•  Exposed: infected with the disease, but not yet infectious,

•  Infectious: both infected and infectious (i.e., able to pass on the disease to others),

•  Recovered: healthy, with long-lasting immunity.

The variables in the model, S(t), E(t), I(t), and R(t), are the (non-dimensional) proportions of the total population in the four states at any time t.  The parameters are b the birth rate, d the natural-causes death rate (assumed equal across the population groups), β the infection rate, γ the recovery rate, and α the infected to infectious rate (so 1/α is the average latency period). All the parameters are assumed to be constant and positive, and measured per simulated day.

(i) Use solve_ivp to solve the SEIR model (6), with parameter values b = d = 0.01, α = 0.25, β = 0.2 and γ = 0.1, for t ∈ [0, 3 × 365], with initial conditions S(0) = 0.99, E(0) = 0.01, I(0) = 0 and R(0) = 0, and plot a graph of the solution versus time.  Explain what the results mean in the context of the infection in the population.

How long does it take for the solution to reach a steady state, where the magnitude of the rate of change of the solution, |dt/dy|, is less than 10-5 ?   (4  marks)

(ii)  Plot a graph of the steady-state  solution of the SEIR model (6) versus β, for β ∈ [0, 0.5], by using solve_ivp to solve the SEIR model (6) until the dynamics reach a steady state for a range of values of β .  The other parameters should be fixed to the values given in part (b)(i).

Explain the graph you have obtained, in the context of the infection in the population.      (4  marks)

(c)  As described in the lecture notes, the Robertson system is an equation-based model of a reaction between three chemicals X , Y , Z, with concentrations x(t), y(t), z(t):

where the reaction rates are k1  = 0.04, k2  = 3 × 107 , k3  = 104 .

Write a Python function that measures the performance of a  solve_ivp solver method,  when solving the Robertson system (7) for t ∈ [0, 5] with initial conditions x(0) = 1, y(0) = z(0) = 0, by returning:

the number of right-hand-side function evaluations of solve_ivp,

the average run time of solve_ivp,

given a solve_ivp method  (the lectures used method= ' RK45 '  and method= ' Radau ' ), and desired relative and absolute tolerances, rtol and atol.

Use your code to plot graphs of

the number of right-hand-side function evaluations of solve_ivp,

the average run time of solve_ivp,

versus relative error tolerance rtol, for at least four solve_ivp solver methods, with values of rtol between 10-4  and 10-8  and atol=1e-3*rtol.

Interpret the results,  and explain how they might influence your choice of method when using solve_ivp.    (7 marks)


热门主题

课程名

mktg2509 csci 2600 38170 lng302 csse3010 phas3226 77938 arch1162 engn4536/engn6536 acx5903 comp151101 phl245 cse12 comp9312 stat3016/6016 phas0038 comp2140 6qqmb312 xjco3011 rest0005 ematm0051 5qqmn219 lubs5062m eee8155 cege0100 eap033 artd1109 mat246 etc3430 ecmm462 mis102 inft6800 ddes9903 comp6521 comp9517 comp3331/9331 comp4337 comp6008 comp9414 bu.231.790.81 man00150m csb352h math1041 eengm4100 isys1002 08 6057cem mktg3504 mthm036 mtrx1701 mth3241 eeee3086 cmp-7038b cmp-7000a ints4010 econ2151 infs5710 fins5516 fin3309 fins5510 gsoe9340 math2007 math2036 soee5010 mark3088 infs3605 elec9714 comp2271 ma214 comp2211 infs3604 600426 sit254 acct3091 bbt405 msin0116 com107/com113 mark5826 sit120 comp9021 eco2101 eeen40700 cs253 ece3114 ecmm447 chns3000 math377 itd102 comp9444 comp(2041|9044) econ0060 econ7230 mgt001371 ecs-323 cs6250 mgdi60012 mdia2012 comm221001 comm5000 ma1008 engl642 econ241 com333 math367 mis201 nbs-7041x meek16104 econ2003 comm1190 mbas902 comp-1027 dpst1091 comp7315 eppd1033 m06 ee3025 msci231 bb113/bbs1063 fc709 comp3425 comp9417 econ42915 cb9101 math1102e chme0017 fc307 mkt60104 5522usst litr1-uc6201.200 ee1102 cosc2803 math39512 omp9727 int2067/int5051 bsb151 mgt253 fc021 babs2202 mis2002s phya21 18-213 cege0012 mdia1002 math38032 mech5125 07 cisc102 mgx3110 cs240 11175 fin3020s eco3420 ictten622 comp9727 cpt111 de114102d mgm320h5s bafi1019 math21112 efim20036 mn-3503 fins5568 110.807 bcpm000028 info6030 bma0092 bcpm0054 math20212 ce335 cs365 cenv6141 ftec5580 math2010 ec3450 comm1170 ecmt1010 csci-ua.0480-003 econ12-200 ib3960 ectb60h3f cs247—assignment tk3163 ics3u ib3j80 comp20008 comp9334 eppd1063 acct2343 cct109 isys1055/3412 math350-real math2014 eec180 stat141b econ2101 msinm014/msing014/msing014b fit2004 comp643 bu1002 cm2030
联系我们
EMail: 99515681@qq.com
QQ: 99515681
留学生作业帮-留学生的知心伴侣!
工作时间:08:00-21:00
python代写
微信客服:codinghelp
站长地图