Finding the spiral sinks and spiral sources of a linear system

the7joker7
Messages
111
Reaction score
0

Homework Statement



Basically, the problem involves a linear system dx/dt = ax + by and dy/dt = -x - y, with a and b being parameters that can take on any real value. Basically, you go through this system for several values of a and b (I did -12 to 12) to find the state at various points. That is, whether or not the point is a saddle, source, sink, perodic, or whatever else. I've found all of that. The only thing I haven't found yet is which sources and sinks are spirals. I know that it's a spiral if a complex number is involved. I've been told two different ways to find this.


The Attempt at a Solution



A: Graph the trace of the system (a + d)*x and the determinant (ad - bc) and the area inside the parabola has the spirals. I'm having a hard time doing this on the graphing tools I've found on the internet.

B: Using the quadratic x^2 - (trace*x) + determinant = 0, find what values of the trace and determinant make is such that x = sqrt(negative number). This equation simplifies to x = sqrt(trace*x - determinant). An x on both sides, so that complicates things...

What would be the easiest way to go about this?
 
Physics news on Phys.org
The behavior of a linear system like you describe is pretty much determined by the eigenvalues of the matrix. Here's a link to an article on Equilibrium in dynamical systems, as these are often called. The section titled Two-Dimensional Space might answer some of your questions.
 
For the most part, I understand the relationship between behavior and eigenvalues. When you have two positive eigenvalues it's a source, two negative it's a sink, so on and so forth...

My issue is, what's the easiest way to find for which values of A and B are there complex values involved.
 
This is pretty hazy in my memory since it was a long time ago that I studied this stuff, but your quadratic (B) brings up some memories. It could be that this equation is the characteristic polynomial for the system, which gives you roots r1 and r2, and from which you get solutions e^(r1*t) and e^(r2 * t). If the solutions to the quadratic are complex, they will be conjugates, so you'll get r = a +/-bi. The imaginary parts will lead to solutions involving sin() and cos(), and that's why you get spiral behavior, either toward a source (if a < 0) or away to a sink (if a > 0).

Again, this is pretty fuzzy in my mind, but maybe I've given you some directions to go in.
Mark
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top