MHB Differential equations stability

Click For Summary
SUMMARY

The discussion focuses on the stability of equilibrium points in a one-dimensional dynamical system defined by the differential equation $x' = f(x)$, where $f(x) = x^4 \sin(1/x)$ for $x \neq 0$ and $f(0) = 0$. The equilibrium points identified are $x = 0$ and $x = \frac{1}{n\pi}$ for all integers $n$. The stability of these points is analyzed using Lyapunov stability criteria, where it is established that if $f'(q) < 0$, the equilibrium is stable, while $f'(q) > 0$ indicates instability. The discussion emphasizes the importance of correctly applying the definitions of stability and utilizing direction fields to visualize the behavior of the system.

PREREQUISITES
  • Understanding of differential equations, specifically first-order equations.
  • Familiarity with concepts of equilibrium points and stability in dynamical systems.
  • Knowledge of Lyapunov stability criteria and how to apply them.
  • Ability to analyze direction fields and phase portraits for dynamical systems.
NEXT STEPS
  • Study the application of Lyapunov's direct method for stability analysis in nonlinear systems.
  • Learn how to construct and interpret direction fields for various types of differential equations.
  • Explore the concept of asymptotic stability and its implications in dynamical systems.
  • Investigate the behavior of solutions near equilibrium points using perturbation methods.
USEFUL FOR

Mathematicians, physicists, and engineers involved in dynamical systems analysis, particularly those focusing on stability analysis of nonlinear differential equations.

Harambe1
Messages
5
Reaction score
0
A one-dimensional dynamical system is given by
$x′ = f(x), t \in [0,+\infty)$,
where $f : \mathbb{R} \to \mathbb{R}$ is the smooth function defined as follows:

$$f(x) = \begin{cases}
x^4 \sin \left(\frac{1}{x}\right) & x \neq 0\\ 0 & x = 0.
\end{cases}.$$

Find all the equilibrium points and determine the stability properties of each equilibrium point.

I've managed to find all the equilibrium points and their stability properties (although using the derivative test I'm not sure if less than zero implies 'stable' or 'asymtopically stable') with the exception of zero. I'm struggling to formally prove its stability using the definition:

The equilibrium point $p \in \mathbb{R}$ is stable if for all $\varepsilon > 0$ there exists a $\delta > 0$ such that $\|x(t)-p\| \leq \delta$ whenever $\|x(t)-p\| \leq \varepsilon$ for all $t \geq t_0 \geq 0$.
I'm not really sure where to start with choosing epsilon or delta or how to use these. Thanks for any tips.
 
Physics news on Phys.org
I think you have not quoted the definition of stability correctly for your application. It's looking like you're trying to apply Lyapunov stability:

Given the DE $\dot{x}=f(x)$, and solution $x=\varphi(t,p)$, a point $q$ is stable if given $\varepsilon > 0$, there is a $\delta > 0$ such that $\|\varphi(t,p)-q\|<\varepsilon$ for all $t>t_0$ and for all $p$ such that $\|p-q\|<\delta$. The idea in my mind is this: if you start the system close enough to a stable equilibrium point, then the system stays in an envelope about that equilibrium point for all time after the initial time.

Now your equilibria are:
$$x=0, \; \frac{1}{n\pi}, \; \forall \, n\in\mathbb{Z}.$$
I would do the $\dot{x}$ versus $x$ plot as done in http://www.math.psu.edu/tseng/class/Math251/Notes-1st%20order%20ODE%20pt2.pdf. You just draw arrows in the various regions, corresponding to whether you're above or below the $x$ axis. When you're above, you draw arrows to the right, and when you're below, draw them to the left. If an equilibrium has arrows on both sides going into it, you're stable. If going out of it, you're unstable. If the arrows are going through it, then you're neither unstable nor stable, but semistable (stable from one direction only). Does that make sense?

Another way of thinking about it is this: if $f'(q)<0$ at equilibrium point $q$, then it's stable. If $f'(q)>0$ at equilibrium point $q$, then it's unstable. If $f'(q)=0$, then $q$ is neither.
 
Ackbach said:
I think you have not quoted the definition of stability correctly for your application. It's looking like you're trying to apply Lyapunov stability:

Given the DE $\dot{x}=f(x)$, and solution $x=\varphi(t,p)$, a point $q$ is stable if given $\varepsilon > 0$, there is a $\delta > 0$ such that $\|\varphi(t,p)-q\|<\varepsilon$ for all $t>t_0$ and for all $p$ such that $\|p-q\|<\delta$. The idea in my mind is this: if you start the system close enough to a stable equilibrium point, then the system stays in an envelope about that equilibrium point for all time after the initial time.

Now your equilibria are:
$$x=0, \; \frac{1}{n\pi}, \; \forall \, n\in\mathbb{Z}.$$
I would do the $\dot{x}$ versus $x$ plot as done in http://www.math.psu.edu/tseng/class/Math251/Notes-1st%20order%20ODE%20pt2.pdf. You just draw arrows in the various regions, corresponding to whether you're above or below the $x$ axis. When you're above, you draw arrows to the right, and when you're below, draw them to the left. If an equilibrium has arrows on both sides going into it, you're stable. If going out of it, you're unstable. If the arrows are going through it, then you're neither unstable nor stable, but semistable (stable from one direction only). Does that make sense?

Another way of thinking about it is this: if $f'(q)<0$ at equilibrium point $q$, then it's stable. If $f'(q)>0$ at equilibrium point $q$, then it's unstable. If $f'(q)=0$, then $q$ is neither.

Thanks for the reply.

I'm not really sure how I would go about drawing a direction field for this particular function so have opted to use the "If $f'(q)>0$ at equilibrium point $q$, then it's unstable. If $f'(q)=0$, then $q$ is neither." Although having obtained $f'(q)=0$ using the equilibrium q=0, the test fails. As I also am unable to sketch the direction field is it possible to use the definition to prove stability at q=o (by choosing an arbitrary epsilon?)

Thanks again.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
1
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K