System of second-order autonomous ODE

In summary, the conversation discusses a system of second-order autonomous ODEs arising from a population genetics model. The question is whether there is a way to solve this system, and the response mentions the possibility of linearizing it and the need for additional details. The conversation then delves into the values of the parameters and the form of the solution being sought. There is also mention of a way to decouple the equations, but it may not yield a simple solution. The conversation then turns to discussing the values of x and y at equilibrium points and the possibility of reducing the system to first-order equations. The conversation ends with a discussion of a specific example and the search for steady-state solutions for a pair of PDEs.
  • #1
dr29
4
0
Hello,
I have a system of two second-order autonomous ODEs arising from a population genetics model:
(a-b y)x(1-x) + x''=0
(a-b x)y(1-y) + y''=0
where a and b are constants, and x, y are fonctions of t.
Is there any hope to solve this system?
Thank you for your help.
 
Last edited:
Physics news on Phys.org
  • #2
One can linearise this if you provide additional details about how fast x,y grow & whether xy is a small quantity.
 
  • #3
Thank you for your response!
Sorry I should have been more specific; the parameters satisfy 0<a<b<1; attached is the form of the solution I'm looking for (with x/y in blue/red). Here the variable is set at 0 when x=y.
Even if there's no simple solution, the value of x=y at 0, and the slopes at 0 would already be quite interesting...
 

Attachments

  • solution.jpg
    solution.jpg
    5.9 KB · Views: 440
  • #4
Well... there is a way to de-couple this set of equations, so one gets one ODE for x and one ODE for y, but i think it does more harm than good, as the equations one gets are of fourth-order, enormous, and probably impossible to solve analitically.

Anyway, it seems possible to get some useful informations about this system from the graph you sent. Looking at it, we see that the slope of y (which is y') starts (at t=-100) at almost zero, then increases near t=0, then decreases again to almost zero as t increases past 10. If y' starts at almost zero, increases then returns to almost zero, it reaches a maximum at somewhere between. At the maximum of y', its derivative (y'') is zero.

If we put y''=0 at the second equation, assuming that 0<y<1, we find that this maximum is reached when x=a/b.

The same applies to the curve of x, the only difference is that x' reaches a minimum instead a maximum. Anyway, x'' still is zero, and so we find that at this point y=a/b.

It is very tempting to assume that, as we found equal values for x and y, it means that this value is the one where both curves intercept (at t=0), but there is no reason why it should be. This values only represent the points where x'' and y'' are zero, but nothing ensures that both x'' and y'' are zero at the same instant of time.

Anyway, if you know the values of a and b, you may verify if x (and y) are actually a/b at t=0. If it were, it may be just a coincidence, or it may be a "right" result, even if by pure luck, or guessing.
 
  • #5
You can always reduce a system of higher order d.e. s to first order.

Let u= x' and v= y' then x''= u'= (a- by)x(1- x) and y''= v'= (a- bx)y(1- y). Those, together with the equationx x'= u and y'= v give a system of four non-linear equations. It is easy to see that the equilibrium solutions are given by (a- by)x(1- x)= 0, (a- bx)y(1- y)= 0, u= 0, v= 0. The first equation has solutions x= 0, x= 1 and y= a/b. If x=0, the second equation is ay(1- y) which hs solutions y= 0 and y= 1. If x= 1, the second equation is (a- b)y(1- y)= 0 which has the same solutions. if y= a/b then the equation is (a- bx)(a/b)(1- a/b)= 0 which has x= a/b as solution. That is, the equilbrium solutions are
(x, y, u, v)= (0, 0, 0, 0)
(x, y, u, v)= (0, 1, 0, 0)
(x, y, u, v)= (1, 0, 0, 0)
(x, y, u, v)= (1, 1, 0, 0)
(x, y, u, v)= (a/b, a/b, 0, 0).

The Jacobian for the system is
[tex]\begin{bmatrix}(a- by)(1- 2x) & -bx(1- x) & 0 & 0 \\ -by(1- y) & (a- bx)(1- 2y) \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1\end{bmatrix}[/tex]

Finding the eigenvalues and eigenvectors for that system at each equilibrium point should tell you about the behaviour of the system around that point.
 
  • #6
OK thanks for your help!
Unfortunately x and y don't reach a/b at the same time, as can be seen in this new figure with a=0.05, b=0.2 (so a/b=0.25).
So with this solution, I have the equilibria (x,y,u,v)=(1,0,0,0), (x,y,u,v)=(0,1,0,0), but not (a/b, a/b, 0,0)...
What would be great would be to have the values of x and x' when y=a/b (and hence x''=0), but I'm not sure if it's feasible...
 

Attachments

  • sol2.jpg
    sol2.jpg
    6.8 KB · Views: 471
  • #7
That's not true. x= a/b, y= a/b, is, as I said before, an equi8librium point- if x= a/b, y= a/b, initially, then they have those values for all t. At the point (a/b, a/b), -a(1- a/b) is an eigenvalue. Whether solutions that start near that will move toward or away from it depends upon whether -a(1- a/b) is positive or negative. In your example, you chose initial values well away from a/b.
 
  • #8
Yes sorry, it is an equilibrium, but it is never reached in the case I'm interested in. Let me be more explicit: I'm looking for a steady-state solution for the pair of PDEs:
dx/dt=(a-b y)x(1-x)+d^2 x/dz^2
dy/dt=(a-b x)y(1-y)+d^2 y/dz^2
The initial condition (at time t=0) is x=y=0 everywhere (for all values of z, which is a spatial variable), except x=epsilon, y=0 at z=-z_init < 0, x=0, y=epsilon at z=+z_init > 0. Two traveling waves develop for x and y starting from -z_init and z_init, and finally collide at z=0, to take the form shown on the previous figures. I'm trying to describe x(z) and y(z) at this equilibrium.
This new figure shows the system at different times (the x-axis is the spatial z-variable... that was called t in previous posts -- sorry for the confusion)
 

Attachments

  • fig3.jpg
    fig3.jpg
    11.7 KB · Views: 519

1. What is a system of second-order autonomous ODE?

A system of second-order autonomous ODE (ordinary differential equations) is a set of two or more equations that describe the change in a set of variables over time. These equations are autonomous, meaning they do not explicitly depend on time, and are of second order, meaning they involve the second derivative of the variables.

2. How is a system of second-order autonomous ODE different from a single second-order ODE?

A single second-order ODE only involves one dependent variable, while a system of second-order autonomous ODE involves multiple dependent variables. Additionally, a system of ODEs is a set of equations, while a single ODE is a single equation. Finally, the solutions to a system of ODEs are a set of functions, while the solution to a single ODE is a single function.

3. What are some real-world applications of a system of second-order autonomous ODE?

A system of second-order autonomous ODE can be used to model a variety of physical and biological systems, such as the motion of a pendulum, the growth of a population, or the spread of a disease. It can also be used in engineering to model control systems, electrical circuits, and chemical reactions.

4. How is a system of second-order autonomous ODE solved?

There is no general method for solving a system of second-order autonomous ODE. However, there are various techniques such as substitution, elimination, and linearization that can be used to simplify the system and find a solution. Numerical methods, such as Euler's method or Runge-Kutta methods, can also be used to approximate solutions.

5. What are some challenges in solving a system of second-order autonomous ODE?

One of the main challenges in solving a system of second-order autonomous ODE is finding an analytical solution, as it may not always be possible. Another challenge is ensuring the stability and accuracy of numerical solutions, especially for complex systems. Additionally, interpreting the solutions and understanding the behavior of the system can be challenging, as the equations may involve multiple variables and parameters.

Similar threads

  • Differential Equations
Replies
3
Views
1K
  • Differential Equations
Replies
2
Views
1K
  • Differential Equations
Replies
7
Views
2K
Replies
3
Views
740
  • Differential Equations
Replies
3
Views
1K
  • Differential Equations
Replies
1
Views
2K
  • Differential Equations
Replies
16
Views
814
Replies
1
Views
900
  • Differential Equations
Replies
1
Views
1K
Replies
28
Views
2K
Back
Top