Help ODE System Stability - Origin Analysis

Click For Summary
SUMMARY

The discussion focuses on analyzing the stability of the origin for a system of ordinary differential equations (ODEs) defined by the equations dx/dt = -2y, dy/dt = x + 2y, and dz/dt = -2z. The eigenvalues of the corresponding 3x3 matrix are -2, 1 + i, and 1 - i, indicating the presence of both stable and unstable behavior. The solution, confirmed using Mathematica, suggests that the origin is unstable due to the presence of terms with e^t in x(t) and y(t), while z(t) decays exponentially. A potential typo in z(t) was noted, as one eigenvalue is -2, which typically indicates a decaying exponential should be present.

PREREQUISITES
  • Understanding of ordinary differential equations (ODEs)
  • Familiarity with eigenvalues and eigenvectors
  • Knowledge of stability analysis in dynamical systems
  • Experience with Mathematica for confirming solutions
NEXT STEPS
  • Study the stability criteria for linear systems of ODEs
  • Learn about phase plane analysis and its application to stability
  • Explore the role of complex eigenvalues in stability analysis
  • Investigate the use of Mathematica for solving and analyzing ODEs
USEFUL FOR

Students and researchers in mathematics, engineering, or physics who are studying dynamical systems and stability analysis of ODEs.

Carl140
Messages
49
Reaction score
0

Homework Statement



Hello.


I want to study the stability of the origin of the following problem:


dx/dt = -2y


dy/dt = x + 2y


dz/dt = -2z


So the eigenvalues of this 3 x 3 matrix are -2, 1 + i, 1-i.


The eigenvectors are (0,0,1) , (2,-1-i,0), (-2,-1+i,0).


The solution (confirmed with Mathematica) is given by:


x(t) = 2*exp(t)cos(t) * C_1 + 2*exp(t)*sin(t) * C_2 - 2*exp(t)*sin(t) *C_3


y(t) = C_1 * (-exp(t)*¨sin(t) - exp(t)*cos(t)) + C_2 * (exp(t)cos(t) -exp(t)sin(t)) +
C_3 * (-exp(t)cos(t)+exp(t)sin(t) )


z(t) = 2*exp(2*t) *C_3


Where C_1,C_2,C_3 are constants.


How can I find (analytically, not by plotting) if the origin (0,0,0) is stable, asymptotically stable? unstable, a node, a center?


I'm having trouble figuring this out. Can you please help?
 
Physics news on Phys.org
Carl140 said:

Homework Statement



Hello.


I want to study the stability of the origin of the following problem:


dx/dt = -2y


dy/dt = x + 2y


dz/dt = -2z


So the eigenvalues of this 3 x 3 matrix are -2, 1 + i, 1-i.


The eigenvectors are (0,0,1) , (2,-1-i,0), (-2,-1+i,0).


The solution (confirmed with Mathematica) is given by:


x(t) = 2*exp(t)cos(t) * C_1 + 2*exp(t)*sin(t) * C_2 - 2*exp(t)*sin(t) *C_3


y(t) = C_1 * (-exp(t)*¨sin(t) - exp(t)*cos(t)) + C_2 * (exp(t)cos(t) -exp(t)sin(t)) +
C_3 * (-exp(t)cos(t)+exp(t)sin(t) )


z(t) = 2*exp(2*t) *C_3


Where C_1,C_2,C_3 are constants.


How can I find (analytically, not by plotting) if the origin (0,0,0) is stable, asymptotically stable? unstable, a node, a center?


I'm having trouble figuring this out. Can you please help?

Is there a typo in z(t)? One of your eigenvalues is -2, so I would expect to see e^(-2t) in one of your solution functions.

Because x(t) and y(t) both have terms with e^t, I would expect orbits that move away from the origin over time, which would make the origin unstable or a node (I don't recall exactly what these terms mean in the context of phase diagrams. And because z(t) is a decaying exponential function, whatever the orbits are doing, they are going to be heading down to the x-y plane over time. I'm just going off the top of my head here, so take what I'm saying with a grain of salt.
 

Similar threads

Replies
8
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
9
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
1
Views
2K
Replies
1
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K