Equilibria analysis for Biomath Problem

Click For Summary
SUMMARY

The discussion focuses on finding positive equilibria for the difference equations defined by x(t+1) = (ax(t)y(t))/(1+x(t)) and y(t+1) = (bx(t)y(t))/(1+y(t)), where a and b are positive constants. The equilibria identified are (0,0) and ((1+a)/(ab-1), (1+b)/(ab-1)) with the condition ab > 1. The Jacobian matrix is analyzed to determine local asymptotic stability, revealing that the trivial solution (0,0) is stable for any parameters a and b, while the second equilibrium requires further investigation into the conditions for stability, specifically the inequalities involving the trace and determinant of the Jacobian.

PREREQUISITES
  • Understanding of difference equations and equilibrium points
  • Familiarity with Jacobian matrices and their role in stability analysis
  • Knowledge of local asymptotic stability criteria
  • Proficiency in mathematical notation and LaTeX for expressing equations
NEXT STEPS
  • Study the properties of Jacobian matrices in the context of difference equations
  • Learn about local asymptotic stability and its implications for equilibria
  • Explore the concept of higher-order terms (H.O.T.) in stability analysis
  • Investigate the conditions for stability in nonlinear systems
USEFUL FOR

Students and researchers in mathematical biology, particularly those working on dynamical systems and stability analysis of difference equations.

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
 

Similar threads

  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 2 ·
Replies
2
Views
1K
Replies
3
Views
2K
Replies
12
Views
2K
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
2
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
Replies
8
Views
1K
  • · Replies 3 ·
Replies
3
Views
2K