Section 3.5: Bifurcations - Combinations of bifurcations

Let’s look at a somewhat more complicated system

x
=rx+
3
x
-
5
x
This is similar to the subcritical pitchfork bifurcation, but now with an extra
5
x
term. Let’s start by looking at the roots of the right hand side (which of course correspond to the fixed points). It’s messy, but it’s possible to show that this has zeros at:
x=0,x=
1
2
(1-
1+4r
)
,x=-
1
2
(1-
1+4r
)
,x=
1
2
(1+
1+4r
)
,x=-
1
2
(1+
1+4r
)
So we have five different zeroes. However, for a given value of r, not all of these give us real numbers. For instance, for
r<-
1
4
, we can see that all of the terms with square roots are undefined, so for this region of
r
there will actually only be one fixed point. The bifurcation diagram is found by plotting all of these functions of r. We will check whether these are stable or unstable fixed points in a bit...
​
In[]:=
x/.Solve[rx+
3
x
-
5
x
0,x]​​Plot[%,{r,-1,1},PlotStyleThick,AxesLabel(Style[#,20]&/@{"r","x"})]
Out[]=
0,-
1-
1+4r
2
,
1-
1+4r
2
,-
1+
1+4r
2
,
1+
1+4r
2

Out[]=
Each color corresponds to one of the above roots of the equation
rx+
3
x
-
5
x
=0
ok, so we can see that we have changes in behaviour at
r=-
1
4
and r = 0. Let’s explore the phase portraits for the different regions. We’ll skip the fixed points at first and just plot the

x
. We will plot it for
r=-
1
2
(below the first bifurcation point), at
-
1
4
(at the first bifurcation point), at
-
1
8
(in between the first and second bifurcation point), at 0 (at the second bifurcation point) and at
1
2
(above the second bifurcation point).
In[]:=
ploptions={AxesLabel{Style["x",14],Style["

x
",14]},AxesOrigin{0,0},AspectRatio0.5};GraphicsGridPlot
-1
2
,
-1
4
,
-1
8
,0,
1
4
[[#]]x+
3
x
-
5
x
,{x,-1.15,1.15},Evaluate[ploptions]&/@Range[5],Spacings-50,ImageSize1200​​
Out[]=
As an exercise, draw on the fixed points and whether they are stable or unstable. From this, you should be able to get the following bifurcation plot:
In[]:=
pl1=ShowPlot0,-
1-
1+4r
2
,
1-
1+4r
2
,-
1+
1+4r
2
,
1+
1+4r
2
,{r,-1,0},AxesLabel(Style[#,20]&/@{"r","x"}),PlotStyle{{Thick,Black},{Thick,Black,Dashed},{Thick,Black,Dashed},{Thick,Black},{Thick,Black}},Plot0,-
1+
1+4r
2
,
1+
1+4r
2
,{r,0,1},PlotStyle{{Thick,Black,Dashed},{Thick,Black},{Thick,Black}},AxesLabel(Style[#,20]&/@{"r","x"}),PlotRangeAll
Out[]=
Notice that we have two saddle-node bifurcations appearing at
r=-
1
4
, and then a subcritical pitchfork bifurcation happening at r=0. Let’s just box these to make sure that it’s clear. The saddle-node bifurcations are in blue, the pitchfork is in red.
In[]:=
Show[pl1,Graphics[{Blue,Thick,Line[{{-0.4,#0.5},{-0.15,#0.5},{-0.15,#0.9},{-0.4,#0.9},{-0.4,#0.5}}]}]&/@{1,-1},Graphics[{Red,Thick,Line[{{-0.1,-0.2},{0.1,-0.2},{0.1,0.2},{-0.1,0.2},{-0.1,-0.2}}]}]]
Out[]=
Now look at the phase portraits and try and match up the different types of fixed points.
​
OK, let’s step back for a moment and try and get a bigger picture...this has gotten pretty messy!
​
We started here with a differential equation which may describe the dynamics of some physical system. We want to know some of the qualitative features of the behaviour of this system. One of the most important types of qualitative thing we can talk about are the fixed points, and their nature. We have a parameter,
r
and we know that for different values of this parameter there may be different numbers of fixed points. We figured out where the fixed points were, and you figured out what kinds of fixed points they were. Of course one can run the linear stability analysis on the fixed points and see whether they are exponential, or less than exponential in their behaviour.
​
Hysteresis
​
Let’s ask now, not just about where the fixed points are, but how you might move around them. Let’s say you start your system off with an r value (which you can tune by a dial on your equipment) of
r=-0.5
, and an x value of
x
0
=0.1
. Close to
x
= 0 you can do a linear stability analysis and show that the equation looks like:

x
≈rx
Which means that we will have a behaviour with these parameters of:
x=0.1
-0.5t
e
.
So we are getting closer and closer to x=0 over time. Let’s say that we now quickly increase our value of r to r=-0.01. Because we are close to x=0, the linear stability analysis around here remains the same and so we will remain close to x=0. Notice that for values of
-1
4
<r<0
we do have another couple of stable fixed points at positive and negative values of x but because we have started close to x=0, these don’t affect us...we can plot our journey with a red arrow on our bifurcation diagram.
In[]:=
Show[pl1,Graphics[{Thick,Red,Dashed,Arrow[{{-0.5,0},{-0.01,0}}]}]]
Because of our exponential decay, we are still at a positive (though small) value of x. It actually depends how quickly we move our dial and how long we’ve been waiting as to how small, but we know that we couldn’t have crossed over to a negative value of x. Let’s say that after waiting at r=-0.5 for 2 seconds, we suddenly altered the value of r to -0.01 and stay there for 2 seconds. We can first work out that after the first 2 seconds our x value would be:
Then our new r value is -0.01 and we wait there for 2 seconds:
Let’s think careful about where that 4-2 comes from. The solution to the new equation is:
but we started now not at t=0, but really at t=2, so everything should be measured from t=2. The amount of time that we have been sat there at t=4 is 2 seconds in the new r=-0.01 system, and the initial condition is the value that we were at the end of the r=-0.5 system, which was 0.0367879.
In terms of the trajectory, this looks like:
Note that in seconds 2 to 4, it looks flat but we are still decaying, just with a decay rate of only -0.01.
​
Now let's say that we move our dial to r=0.5. All of a sudden the value of x that we are sitting at is an unstable fixed point. In fact we will move away from it now with an exponential of the form:
This takes us to 14 seconds and to an x value of 1.16877
On the bifurcation diagram it looks like we have made this journey:
Dashed arrows are instantaneous jumps in r, the full red arrows are trajectories which take time to get to a new x value.
ok, so what happens if we now move the value of r back -0.01? We could do this by either slowly tuning r down, in which case it’s going to stick to the positive stable fixed point:
Then we find ourself at the upper fixed point again:
Let’s just add to this the values of r across time as well. These are in red:
Which corresponds to:
​
You could have done the above, but changed the value of r slowly, in which case you would have:
​
Here is a really strange example of hysteresis: If you take water at room temperature and freeze it, it will take a certain amount of time to freeze. However, if you boil water, and then freeze it, it actually freezes faster, despite the fact that in order to freeze, it has to pass through the room temperature state. There is information in the system about what happened to it in the past (you can think of it as remembering that it was in a high temperature state - though note very very importantly, this “memory” is not like the memory which people claim that water has in homeopathy. That isn’t memory, it’s wishful thinking). This phenomenon with water is called the Mpemba Effect, named after the Tanzanian schoolboy who discovered it. It’s actually not completely agreed upon what is the cause of this effect.
​
​
That was a lot! We’re going to look at something a bit different on the next page before returning to the study of these bifurcations.
​
Next page
Homepage