Section 2.1: Fixed Points and Trajectories
Section 2.1: Fixed Points and Trajectories
Flows on the Line - Fixed Points and Trajectories
We’re actually going to focus on the simplest system of differential equations - a single differential equation. We’re going to see what interesting discoveries we can make in such simple settings. For this chapter we will be dealing with equations of the form:
x
Please note:
- We use as the shorthand for
x
dx
dt
- We leave out the functional dependence on t (i.e. write rather than (t) most of the time)
x
x
Recall, we can have equations where is itself a function of and explicitly ( =x+), but these non-autonomous systems can always be rewritten as autonomous systems of one more equation (we won't be dealing with these).
f
x
t
x
2
t
Here, is a function: (For our purposes it’s always going to be differentiable)
x
x:⟶
This is an equation on a line, or a flow on a line, because you can imagine that is the position (let’s say of a particle) on a line at time t, and the dynamical equation tells you how you move along the line. In particular, your velocity (which is just ) at any particular position (velocity is change in displacement, x, over time).
x
x
We’re going to start with an example that is possible to solve exactly, but we’ll see that we can get a lot out of it without doing so.
Our first equation is:
x
Think of our particle on a line, and this just tells us that the velocity along the line at any one point is given by the sin of that point’s position.
It is actually possible to solve this, but it is a pain... and we can be lazy. I’ll write out the solution, but you don’t need to know this!!
x=2cot+2πn()
-1
cot
-t
e
x
0
2
x
0
where is the position at time , and we have a function which we won’t go into, but it’s always an integer. As an exercise, you can show that this is a solution by plugging it into both sides of our differential equation above:
x
0
x=0
n()
x
0
d2cot+2πn()
-1
cot
-t
e
x
0
2
x
0
dt
-1
cot
-t
e
x
0
2
x
0
If you take the derivative, and do some serious manipulation you can show that this is indeed a solution.
Now frankly, that’s not easy to look at and know what is going on for any particular starting position. We are going to start developing some graphical methods for understanding such systems, instead of this complex route.
The first thing that we’re going to do is to plot against . Well, in this case is just , so we have:
x
x
x
sin(x)
In[]:=
equation1===Sin[x]pl1=Plot[equation1[[2]],{x,-10,10},AxesLabel(Style[#,20]&/@{"x",""})]
x
x
Out[]=
x
Out[]=
Well what does this mean?
- tells us how fast we are moving along the line, for a given (the velocity). If we want to ask how fast we are going at, say the position , we just read off the value of :
x
x
x=2
x
In[]:=
Show[pl1,Graphics[Line[{{2,0},{2,Sin[2]},{0,Sin[2]}}]],Graphics[Text[Style["x=2",16],{2,-0.1}]],Graphics[Text[Style["=sin(2)",16],{-2,Sin[2]-0.05}]]]
x
Out[]=
So the velocity of our particle, when we are at position is =sin(2)≈0.91 (i.e. moving right, with speed roughly 0.91). Now let's ask the same question at :
x=2
x
x=4
In[]:=
Show[pl1,Graphics[Line[{{4,0},{4,Sin[4]},{0,Sin[4]}}]],Graphics[Text[Style["x=4",16],{4.1,0.1}]],Graphics[Text[Style["=sin(4)",16],{2,-0.9}]]]
x
Out[]=
So the velocity of our particle, when we are at position is =sin(4)≈-0.76 (i.e. moving left, now at a slightly slower speed of 0.76)
x=4
x
In fact, we can see that whenever is positive, we are moving to the right, and whenever is negative we would be moving to the left. We put arrows on our line to show direction of motion in each region:
x
x
In[]:=
Show[pl1,Table[Graphics[Arrow[{{i,0},{i+#,0}}]]&/@{2.5,-2.5},{i,-2π,2π,2π}]]
Out[]=
We can see something else from the plots. The closer we get to an intersection point of the sinusoidal curve with the x-axis, the smaller the magnitude of our velocity will be.
In fact, we can see that when we are at an intersection point, our velocity will be zero (i.e. =0). These are called Fixed Points. If a particle (a flow) starts at one of these special points, it never moves from it, hence it is fixed to the point.
x
For now, let's put a black dot at every fixed point. In a bit we will see that not all fixed points are equal...
In[]:=
Show[pl1,ListPlot[{#,0}&/@Range[-3π,3π,π],PlotStyle{Black,PointSize[0.02]}]]
Out[]=
If I put you at a particular point. You could work out how fast and in which direction you would be moving, you can imagine going in that direction and then the next moment you will be somewhere else, and moving at a slightly different speed. Let’s zoom into one region (between and ) and see what happens if I dropped you at just a little after , let’s say at . Our velocity at this point is going to be =sin(0.1)≈0.1 to the right.
x=0
x=π
x=0
x=0.1
x
In[]:=
Show[Plot[equation1[[2]],{x,0,π},AxesLabel(Style[#,20]&/@{"x",""})],Graphics[{Arrowheads[Medium],Arrow[{{0.1,0},{0.1+Sin[0.1],0}}]}],ListPlot[{{0.1,0}},PlotStyle{Red,PointSize[0.02]}]]
x
One really important thing to keep in mind, is that we are only moving along the x-axis. The vertical axis gives us information about how we are moving, but nothing more.
This same thing happens when starting at any point we want. There are certain points that you will always move away from, and certain points that you will always move towards (but never arrive at). The points that we move away from are called unstable fixed points, and those that you move towards are called stable fixed points.
Stability: The idea that small shakes or perturbations don’t cause large changes. If we start a particle/flow exactly on a fixed point, it never moves.
If we perturb the flow and it moves away from a fixed point then it is an unstable fixed point and if it moves back towards it then it is a stable fixed point.
Stability: The idea that small shakes or perturbations don’t cause large changes. If we start a particle/flow exactly on a fixed point, it never moves.
If we perturb the flow and it moves away from a fixed point then it is an unstable fixed point and if it moves back towards it then it is a stable fixed point.
Think of balancing a ball at the top of a hill.
- In theory you can balance perfectly there, but any little wobble and you will fall down (unstable)
Conversely, imagine a ball at the bottom of a valley.
- Any little push away and it will roll back down to the bottom of the valley (stable)
Note that a system like a ball on a hill is not the same as our one-dimensional system. However, to get an idea of stable and unstable fixed points it’s a useful analogue.
We can label the stable and unstable fixed points separately:
- A solid black disk at a stable fixed point.
- An open circle at the unstable fixed point.
What are the trajectories going to look like starting with a particular initial condition?
We can plot different trajectories on the same graph:
Spot the stable and unstable fixed points here!
And we can follow along with one of the trajectories to see what the actual particle motion looks like:
Now let’s move on to a more general investigation of these methods.