Proving stability of linear system equations

Click For Summary
SUMMARY

The discussion centers on proving the stability of linear system equations through various mathematical approaches. The characteristic equation derived is ##\lambda^2 + 2f \lambda + f^2 + 1 = 0##, yielding eigenvalues ##\lambda_1 = -f - i## and ##\lambda_2 = -f + i##. Participants suggest using polar coordinates for analysis, indicating that if ##\dot r < 0##, the origin is asymptotically stable, while ##\dot r > 0## indicates instability. The use of a Liapunov function, specifically ##V(x,y) = dx^2 + dy^2##, is proposed as a method for stability analysis.

PREREQUISITES
  • Understanding of characteristic equations in linear systems
  • Familiarity with eigenvalues and their significance in stability analysis
  • Knowledge of polar coordinates and their application in dynamical systems
  • Experience with Liapunov functions for stability proofs
NEXT STEPS
  • Study the derivation and implications of characteristic equations in linear systems
  • Learn about the application of polar coordinates in analyzing dynamical systems
  • Research the construction and use of Liapunov functions for stability analysis
  • Explore the Jacobian matrix and its role in linearization of non-linear systems
USEFUL FOR

Mathematicians, engineers, and students focused on control theory, dynamical systems, and stability analysis of differential equations.

member 731016
Homework Statement
Please see below
Relevant Equations
Please see below
For this problem,
1716877322765.png

My solution is to find the characteristic equation of the system by putting the system into a matrix. This gives ##\lambda^2 + 2f \lambda + f^2 + 1 = 0##
Then each eigenvalue is ##\lambda_1 = -f - i## and ##\lambda_2 = -f + i##

I then want to find the Jacobian, however, I would need to find the partial derivatives (with respect to x and y) of ##F(x,y) = y - xf(x,y)## and ##G(x,y) = - x - yf(x,y)##, however, I'm not sure how to do that with the ##f(x,y)## in there.

Does anybody please know what I should do?

Thanks!
 
Physics news on Phys.org
You don't need to look at the Jacobian; the fact that the question tells you nothing about the partial derivatives of f is a strong suggestion that this is the wrong way to proceed.

When I see scalar multiples of (x,y) and (-y,x) on the right hand side, I immediately think of polar coordinates (x,y) = (r \cos \theta, r \sin \theta). This is because \begin{split}<br /> r\frac{dr}{dt} &amp;= x\frac{dx}{dt} + y\frac{dy}{dt} \\<br /> r^2\frac{d\theta}{dt} &amp;= x\frac{dy}{dt} - y\frac{dx}{dt} \end{split} so that the coefficient of (x,y) tells you how r behaves and the coefficient of (-y,x) tells you how \theta behaves. If \dot r &lt; 0 the origin is asymptotically stable; if \dot r &gt; 0 the origin is unstable.
 
Last edited:
  • Love
Likes   Reactions: member 731016
pasmith said:
You don't need to look at the Jacobian; the fact that the question tells you nothing about the partial derivatives of f is a strong suggestion that this is the wrong way to proceed.

When I see scalar multiples of (x,y) and (-y,x) on the right hand side, I immediately think of polar coordinates (x,y) = (r \cos \theta, r \sin \theta). This is because \begin{split}<br /> r\frac{dr}{dt} &amp;= x\frac{dx}{dt} + y\frac{dy}{dt} \\<br /> r^2\frac{d\theta}{dt} &amp;= x\frac{dy}{dt} - y\frac{dx}{dt} \end{split} so that the coefficient of (x,y) tells you how r behaves and the coefficient of (-y,x) tells you how \theta behaves. If \dot r &lt; 0 the origin is asymptotically stable; if \dot r &gt; 0 the origin is unstable.
Thank you for your reply @pasmith!

That is a interesting idea that I have not seen before. They have taught me to me so far to methods to find the stability of a non-linear system of equations, to either use the Jacobian matrix for linearization to find a equivalent linear DE system to the non-linear DE system or use the direct method.

If I think about for the later method, we could also try to solve this problem by using a Liapunov function of the form ##V(x,y) = dx^2 + dy^2 = d(x^2 + y^2)##? If that does not work then I could generalize to more coefficients, ##V(x,y) = dx^2 + gy^2##?

Thanks!
 

Similar threads

Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
1K
Replies
2
Views
2K
Replies
2
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
8
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
3
Views
2K