Dynamical system maps
Alexei Gilchrist
1 Introduction
There a neat way of visualising the dynamics of small-dimension deterministic systems, especially nonlinear dissipative systems where we may have no hope of finding closed form solutions. Recall that in all the solutions to one-dimensional systems we have found there are two free parameters to be determined by the initial conditions. In effect the freedom we have is in specifying the initial position and velocity of the system. This is a generic feature of second order differential equations. Since these two quantities are independent conditions we can plot them against each other—a point on this graph will represent a particular starting value and it will determine the motion for all future times. That is, from that point on we can then track the evolution of the systems as a path or trajectory through the space of states. This graph is a plot of the phase-space of the system (well technically it should be a co-ordinate vs the canonical momentum for that co-ordinate but velocity will do for our purposes).
With this new tool at our disposal we could sample the phase space and draw a small arrow at each point indicating the direction the system’s state would take if it where to have that initial value, and in this way we can draw a `map’ of the dynamics and see how the system would behave with any initial condition at a glance!
In practice though, we can just map out some key features in this space, such as the location of fixed points (starting points where the system does not change in time), the general flow of trajectories in different regions, and the curves that separate different flows which are called a separatrix. This map would then give a very good idea of the dynamics without ever having solved the equations. It should also be easy to see how to extend this description to the multidimensional phase spaces describing more complex systems.
2 Damped harmonic oscillator
As a first example we will determine the phase space of a damped harmonic oscillator (see damped-harmonic-oscillator) for which we already know an explicit general solution. The equation of motion for this system is
The general solution to a damped harmonic oscillator is
This plot should be familiar from the phasor diagram for sinusoid motion, the only difference here is that the shape is slightly squashed in one direction (see caption for explanation), and the phasor is slowly getting shorter. The complicated looking solutions above are the \(x\) and \(v\) projections of this curve.
What we would really like to do is be able to determine this motion without solving the equation of motion. Luckily there is a technique for doing just this. In order to draw an arrow pointing in the direction of the flow at some point \((x,v)\) in phase space, we need to know the direction of the change in time at that point. In other words, we need to know the time derivatives \((\dot{x},\dot{v})\).
It’s always possible to reduce an \(n^\mathrm{th}\) order differential equation to \(n\) coupled first order differential equations by introducing extra variables. In the case of \eqref{eq:damped-eqn-motion} an obvious new variable to introduce is the velocity: \(\dot{x}=v\). Then rewriting the remainder to use \(v\) gives two coupled first order equations,
The first thing we should do is find any fixed points for the system. By definition these points don’t change with time so \(\dot{x}=\dot{v}=0\). The two equations above imply that \(v=0\) and consequently \(x=0\). As this is the only solution, the point \((0,0)\) is the only fixed point.
To figure out the direction of flow at a particular point we just use the coupled differential equations. Say we want to know in which direction a trajectory at the position \((x_0,0)\) will move. Plugging these values in gives
You can explore this map below.
3 Damped Pendulum
Consider now the equation of motion for a damped rigid pendulum (pendulum):
Following the prescription above, we can introduce a new variable \(v=\dot{\theta}\), then split the equation into coupled differential equations
Setting \(\dot{\theta}=\dot{v}=0\) we find two fixed points at \((\theta=0,v=0)\), and \((\theta=\pi,v=0)\). The first corresponds to the pendulum hanging straight down and the second the the pendulum balanced pointing straight up. With any amount of damping (\(\gamma\ge 0\)) all the trajectories flow into the first fixed point.