Equilibria analysis for Biomath Problem

end3r7
Messages
168
Reaction score
0

Homework Statement


Find all the positive equilibria for the system of the difference equations
x(t+1) = \frac{ax(t)y(t)}{1+x(t)};<br /> <br /> y(t+1) = \frac{bx(t)y(t)}{1+y(t)}
a, b > 0.
Then determine conditions on parameters so the equilibriums are locally asymptotically stable.

Homework Equations


An equilibrium point will have corresponding eigenvalues < 1 iff |Trace(J)| &lt; 1 + det(J) &lt; 2,
where J is the Jacobian matrix of the system evaluated at the equilibrium point.
In this case J = \left(<br /> \begin{array}{cc}<br /> \frac{ay(t)}{1+x(t)} - \frac{ax(t)y(t)}{(1+x(t))^2} &amp; \frac{ax(t)}{1+x(t)}\\<br /> \frac{by(t)}{1+y(t)} &amp; \frac{bx(t)}{1+y(t)} - \frac{bx(t)y(t)}{(1+y(t))^2}<br /> \end{array}<br /> \right)<br />

The Attempt at a Solution


I found the following equilibria
(0,0) and (\frac{1+a}{ab-1},\frac{1+b}{ab-1}); ab &gt; 1
The Jacobian Matrix evaluated at (0,0) gives me a zero matrix. So that means the trivial solution is always stable for any choice of parameters a,b, right? Or do I have to look at H.O.T.?

The Jacobian matrix for the second point is:
J = \left(<br /> \begin{array}{cc}<br /> \frac{ab-1}{ab+a} &amp; \frac{1+a}{1+b}\\<br /> \frac{1+b}{1+a} &amp; \frac{ab-1}{ab+b}<br /> \end{array}<br /> \right)

or, if I let x and y denote the equilibrium points

J = \left(<br /> \begin{array}{cc}<br /> \frac{1}{ay} &amp; \frac{x}{y}\\<br /> \frac{y}{x} &amp; \frac{1}{bx}<br /> \end{array}<br /> \right)

And since the trace is always positive, we have the following requirements
ax + by &lt; 1 &lt; 2aybx
Or (ab)(xy) &gt; \frac{1}{2}
and ay + bx &lt; 1

Does that make any sense?
If we substitute x and y back in though, we get nonsense (something like, pos number < -1).

What am I doing wrong?
 
Last edited:
Physics news on Phys.org
Sorry about the confusion fellas, first time using latex in a long time.
Any questions about my work, I'll be gladly to answer. And if possible, don't tell me the answer, but rather around where my thinking goes wrong.
Thanks in advance,
 
Does the zero matrix mean that I have to consider the nonlinear terms of the Jacobian? Cuz that would suck =P
 
end3r7 said:
Does the zero matrix mean that I have to consider the nonlinear terms of the Jacobian? Cuz that would suck =P

Oops, I'm thinking Differential Equations, where the nonhyperbolic/stationary case is when the Jacobian Matrix is the zero Matrix. For difference, I would have to consider H.O.T. if the Jacobian was the identity, which it clearly isn't.
So I guess (0,0) is asymptotically stable, which makes sense if you look a the equations I guess. For x, y small it should converge.

Am I right?
 
What is the formal definition of asymptotically stable?
 
Well, locally asymptotically stability implies that it's both stable and attracting.

That is, given a certain distance from equilibrium point, it will not ever exceed that distance (stable), and will eventually converge to the equilibrium (attracting).

For the vector function X(t+1) = F(X(t)), where X(t) is a vector depending on x,y, whatever...
An equilibrium T is locally stable if for any E > 0 there exists D > 0 such that if |X(0) - T| < D, then |X(t) - T| < E
It's attracting if lim X(t) = T.
 
Can you show formally that either of your equilibria satisfy the definition?
 
EnumaElish said:
Can you show formally that either of your equilibria satisfy the definition?

Yes. But it really doens't matter, I know I"m supposed to use it, because that's what the book tells me. =P But yea, I can show it does.
 
Anyway, the identity that's giving me problems is the following

ay + bx &lt; 1

which is equivalent to

a + 2ab + b < ab - 1
or a + ab + b + 1 < 0, which is impossible since a,b>0
 
Back
Top