This problem considers a cardiovascular model due to Ottesen
[16] involving the arterial pressure, Pa(t) = y1(t), the
venous pressure, Pv(t) = y2(t), and the heart rate, H(t) = y3(t). Ottesen studies conditions under which the delay causes
qualitative differences in the solution and in particular,
oscillations in Pa(t). Delays t = 1.0, 1.4, 3.9, 5.0, 7.5,10 are considered in [16]. There are a number of
parameters in the model, so you might wish simply to declare them
as global values rather than pass them through dde23 to
the function for evaluating the equations, or even to hard code
them. Plot y1(t) = Pa(t) for several values of t. You
should find that the solutions obtained for different values of t differ dramatically. Solve on [0,350] the equations
| ||||||||||||||||||||||||||||||
| |||||||||||||||||||||
| ||||||||||||||||||||||||||||||
One of the figures of [16] shows the solution components
when the peripheral pressure R is reduced exponentially from its
value of 1.05 to 0.84 beginning at t = 600. For this computation
the delay was 4 and the interval [0,1000]. You can easily
modify the previous program to solve this problem. All you have to
do is inform the solver of the low-order discontinuity at a known
time by setting the value of the 'Jumps' option to 600,
modify the function for evaluating the DDEs to include
if t <= 600
R = 1.05;
else
R = 0.21 * exp(600-t) + 0.84;
end
and use the specified delay and interval. All the solution
components are of interest. The figure shows the sharp change in
the heart rate due to the change in R at t = 600.
