Bifurcation diagram

related topics
{math, energy, light}
{math, number, function}
{rate, high, increase}

In mathematics, particularly in dynamical systems, a bifurcation diagram shows the possible long-term values (equilibria/fixed points or periodic orbits) of a system as a function of a bifurcation parameter in the system. It is usual to represent stable solutions with a solid line and unstable solutions with a dotted line.


Bifurcations in the 1D discrete dynamical systems

Logistic map

An example is the bifurcation diagram of the logistic map:

The bifurcation parameter r is shown on the horizontal axis of the plot and the vertical axis shows the possible long-term population values of the logistic function. Only the stable solutions are shown here, there are many other unstable solutions which are not shown in this diagram.

The bifurcation diagram nicely shows the forking of the possible periods of stable orbits from 1 to 2 to 4 to 8 etc. Each of these bifurcation points is a period-doubling bifurcation. The ratio of the lengths of successive intervals between values of r for which bifurcation occurs converges to the first Feigenbaum constant. The code in MATLAB can be written as:

r_values = (2:0.0002:4)';
iterations_per_value = 10;
y = zeros(length(r_values), iterations_per_value);
y0 = 0.5;
y(:,1) = r_values.*y0*(1-y0);
for i = 1:iterations_per_value-1
    y(:,i+1) = r_values.*y(:,i).*(1-y(:,i));
plot(r_values, y, '.', 'MarkerSize', 1);
grid on;

[edit] Real quadratic map

For x_{n+1}=x_n^2-c, the code in MATLAB can be written as:

c = (0:0.001:2)';
iterations_per_value = 100;
y = zeros(length(c), iterations_per_value);
y0 = 0;
y(:,1) = y0.^2 - c;
for i = 1:iterations_per_value-1
    y(:,i+1) = y(:,i).^2 - c;
plot(c, y, '.', 'MarkerSize', 1, 'MarkerEdgeColor', 'black');

[edit] Symmetry breaking in bifurcation sets

Symmetry breaking in pitchfork bifurcation as the parameter epsilon is varied. epsilon = 0 is the case of symmetric pitchfork bifurcation.

In a dynamical system such as

 \ddot {x} + f(x;\mu) + \epsilon g(x) = 0,

which is structurally stable when  \mu \neq 0 , if a bifurcation diagram is plotted, treating μ as the bifurcation parameter, but for different values of ε, the case ε = 0 is the symmetric pitchfork bifurcation. When  \epsilon \neq 0 , we say we have a pitchfork with broken symmetry. This is illustrated in the animation on the right.


Full article ▸

related documents
Lorenz attractor
Triangle wave
Grashof number
List of brightest stars
Log-periodic antenna
Knife-edge effect
Circular definition
Optical density
SN 1604
Vela (constellation)
Isaac Barrow
Rankine scale
Dactyl (moon)
Pulse duration
Ejnar Hertzsprung
Puck (moon)
Lupus (constellation)
Electromagnetic environment
Desdemona (moon)
Igor Tamm
Barn (unit)
Faraday constant
Primary time standard