Modeling - Math 056 - Homework 2
1. In class, we determined that the model for annual plants could be condensed into a single equation:
pn+1 = ασγpn + βσ2(1 - α)γpn-1.
(a) Show that the model can be written as a single equation in s1n. (Write s1n+1 as a function of s1n and s1n-1).
(b) Explain the term ασγpn.
2. Suppose that we define the beginning of a generation at the time when seeds are produced.
(a) Define the variables and write new equations linking them.
(b) Derive the system as two equations for S0n+1 and S¯1n+1. Write the solution in matrix form.
i. Solve the system found in (b).
ii. Show that pn+2 - ασγpn+1 - βσ2γ(1 - α)pn = 0. Is it important when we define the start of a generation?
iii. Find the solution of (ii).
iv. Find the condition for survival of the adult plant population. Biologically interpret this condition.
(c) Consider the problem in which seeds S0n can now survive and germinate in generations n, n+ 1, n+ 2, and n+ 3, with germination fractions α, β, δ, and ∈. Derive the corresponding system of equations in matrix form,
Interpret the terms in the matrix.
3. Solve the first order equation xn+1 = axn + b (a, b constants).
4. Solve and describe the solutions to the following difference equations. (You can check solutions with matlab simulations.)
(a) xn+1 = 3xn + 2yn, yn+1 = xn + 4yn
(b) xn+1 = -xn + 3yn, yn+1 = yn/3
(c) xn+2 - xn+1 + xn = 0
5. Consider the difference equation xn+2 - 3xn+1 + 2xn = 0.
(a) Show that the general solution to this equation is
xn = A1 + 2nA2.
(b) Now suppose that x0 = 10 and x1 = 20. Then A1 and A2 must satisfy the system of equations
A1 + 20A2 = x0 = 10,
A1 + 21A2 = x1 = 20.
(c) Solve for A1 and A2 and find the solution to the above initial value problem.
6. In the circulatory system, the red blood cells (RBCs) are constantly being destroyed and replaced. Since these cells carry oxygen throughout the body, their number must be maintained at some fixed level. Assume that the spleen filters out and destroys a certain fraction of the cells daily and that the bone marrow produces a number proportional to the number lost of the previous day. What would be the cell count on the nth day?
(a) Derive a model for the number of RBCs in circulation and the number of RBCs produced by the marrow. Don't forget to define your variables and parameters.
(b) For homeostasis in the red cell count, the total number of RBCs should remain roughly constant. Explain why one way of achieving this is by letting λ1 = 1.
(c) Using the results of part (c), find λ2. What then is the behavior of the solution?
7. In class we discussed solving homogeneous equations (axn+2 + bxn+1 + cxn = 0). This problem will show you how to solve nonhomogeneous equations.
(a) Consider the difference equation
aXn+2 + bXn+1 + cXn = d.
If d ≠ 0, the equation is called nonhomogeneous. Show that Xn = K, where K is a particular constant, will solve this second order nonhomogenous equation (provided λ = 1 is not an eigenvalue) and find the value of K. This is called a particular solution.
(b) Suppose the solution to the corresponding homogeneous equation aXn+2 + bXn+1 + cXn = 0 is
Xn = c1λn1 + c2λn2,
where c1 and c2 are arbitrary constants. This is called the complementary or homogenous solution. Show that
Xn = K + c1λn1 + c2λn2
will be a solution to the nonhomogeneous problem (provided λ1 ≠ λ2 and λi ≠ 1). This solution is called the general solution.
8. In the blood there is a steady production of CO2 that results from the basal metabolic rate. CO2 is lost by way of the lungs at a ventilation rate governed by CO2-sensitive chemoreceptors located in the brainstem. In reality, both the rate of breathing and the depth of breathing (volume of a breath) are controlled physiologically. We will simplify the problem by assuming that breathing takes place at constant intervals t, t + τ , t + 2τ , ... and that the volume, Vn is controlled by the CO2 concentration in the blood in the previous time interval, Cn-1. Keeping track of the two variables Cn and Vn might lead to the following equations:
Cn+1 = Cn - ?(Vn, Cn) + m,
Vn+1 = ρ(Cn)
where ? is the amount of CO2 lost as a function of volume and concentration, ρ is the sensitivity to blood CO2 as a function of concentration, and m is the constant production rate of CO2 in the blood.
(a) As a first model, assume that the amount of CO2 lost, ?(Vn, Cn) is simply proportional to the ventilation volume Vn with constant factor β (and does not depend on Cn). Further assume that the ventilation at time n + 1 is directly proportional to Cn (with factor α), i.e., that ρ(Cn) = αCn. Write down the system of equations and show that it corresponds to a single equation
Cn+1 - Cn + αβCn-1 = m.
(b) For m ≠ 0 the equation in part (a) is a nonhomogeneous problem. Use the steps outlined in problem 7 to solve it.
i. Show that Cn = m/αβ is a particular solution.
ii. Find the general solution.
(c) Now consider the nature of this solution in two stages.
i. First assume that 4αβ < 1. Interpret this inequality in terms of the biological process. Give evidence for the assertion that under this condition, a steady blood CO2 level C equal to m/αβ will eventually be established, regardless of the initial conditions. What will the steady ventilation rate then be?
ii. Now suppose 4aαβ > 1. Show that the CO2 level will undergo oscillations. If αβ is large enough, show that the oscillations may increase in magnitude. Comment on the biological relevance of this solution.
(d) Suggest how the model might be made more realistic.