Feigenbaum constants

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

In mathematics, specifically bifurcation theory, the Feigenbaum constants are two mathematical constants which both express ratios in a bifurcation diagram for a non-linear map. They are named after the mathematician Mitchell Feigenbaum.

History

Feigenbaum originally related the first constant to the period-doubling bifurcations in the logistic map, but also showed it to hold for all one-dimensional maps with a single quadratic maximum. As a consequence of this generality, every chaotic system that corresponds to this description will bifurcate at the same rate. It was discovered in 1978.[1]

The first constant

The first Feigenbaum constant is the limiting ratio of each bifurcation interval to the next between every period doubling, of a one-parameter map

x_{i+1} = f(x_i)

where f(x) is a function parameterized by the bifurcation parameter a.

It is given by the limit:[2]

\delta = \lim_{n\rightarrow \infty} \dfrac{a_{n-1}-a_{n-2}}{a_n-a_{n-1}} = 4.669\,201\,609\,\cdots

where an are discrete values of a at the nth period doubling.

According to (sequence A006890 in OEIS), this number to 30 decimal places is: δ = 4.669 201 609 102 990 671 853 203 821 578(...).

Illustration

Non-linear maps

To see how this number arises, consider the real one-parameter map:

f(x)=a-x^2.

Here a is the bifurcation parameter, x is the variable. The values of a for which the period doubles (e.g. the largest value for a with no period 2 orbit, or the largest a with no period 4 orbit),are a1, a2 etc. These are tabulated below:[3]

n Period Bifurcation parameter (an) Ratio \dfrac{a_{n-1}-a_{n-2}}{a_n-a_{n-1}}
1 2 0.75 N/A
2 4 1.25 N/A
3 8 1.3680989 4.2337
4 16 1.3940462 4.5515
5 32 1.3996312 4.6458
6 64 1.4008286 4.6639
7 128 1.4010853 4.6682
8 256 1.4011402 4.6689

The ratio in the last column converges to the first Feigenbaum constant. The same number arises for the logistic map

 f(x) = a x (1-x)

with real parameter a and variable x. Tabulating the bifurcation values again:[4]

n Period Bifurcation parameter (an) Ratio \dfrac{a_{n-1}-a_{n-2}}{a_n-a_{n-1}}
1 2 3 N/A
2 4 3.4494897 N/A
3 8 3.5440903 4.7514
4 16 3.5644073 4.6562
5 32 3.5687594 4.6683
6 64 3.5696916 4.6686
7 128 3.5698913 4.6692
8 256 3.5699340 4.6694

Fractals

Self similarity in the Mandelbrot set shown by zooming in on a round feature while panning in the negative-x direction. The display center pans from (−1, 0) to (−1.31, 0) while the view magnifies from 0.5 × 0.5 to 0.12 × 0.12 to approximate the Feigenbaum ratio

In the case of the Mandelbrot set for complex quadratic polynomial

 f(z) = z^2 + c

the Feigenbaum constant is the ratio between the diameters of successive circles on the real axis in the complex plane (see animation on the right).

n Period = 2n Bifurcation parameter (cn) Ratio = \dfrac{c_{n-1}-c_{n-2}}{c_n-c_{n-1}}
1 2 -0.75 N/A
2 4 -1.25 N/A
3 8 -1.3680989 4.2337
4 16 -1.3940462 4.5515
5 32 -1.3996312 4.6458
6 64 -1.4008287 4.6639
7 128 -1.4010853 4.6682
8 256 -1.4011402 4.6689
9 512 -1.401151982029
10 1024 -1.401154502237
infinity -1.4011551890 ...

Bifurcation parameter is a root point of period = 2n component. This series converges to the Feigenbaum point c = −1.401155 The ratio in the last column converges to the first Feigenbaum constant.

Other maps also reproduce this ratio, in this sense the Feigenbaum constant in bifurcation theory is analogous to pi (π) in geometry and Euler's number e in calculus.

The second constant

The second Feigenbaum constant (sequence A006891 in OEIS),

\alpha = 2.502907875095892822283902873218...,

is the ratio between the width of a tine and the width of one of its two subtines (except the tine closest to the fold). A negative sign applied to \alpha when the ratio between the lower subtine and the width of the tine is measured.[5]

These numbers apply to a large class of dynamical systems (for example, dripping faucets to population growth).[5]

Properties

Both numbers are believed to be transcendental, although they have not been proven to be so.[6]

The first proof of the universality of the Feigenbaum constants carried out by Lanford[7] (with a small correction by Eckmann and Wittwer,[8]) was computer assisted. Over the years, non-numerical methods were discovered for different parts of the proof aiding Lyubich in producing the first complete non-numerical proof.[9]

Approximations

Though there is no known closed form equation or infinite series that can exactly calculate either constant, there are closed form approximations for several digits. One of the most accurate, up to six digits, is (sequence A094078 in OEIS)

 \pi + \tan^{-1}(e^{\pi})

which is accurate up to 4.669202. Two closely related expressions that accurately estimate both \delta and \alpha to three decimal places are given in [10]

 \frac{2\varphi}{\log 2} \approx 4.669
 \frac{2\varphi+1}{\log 2+1} \approx 2.502

where \varphi is the golden ratio and \log 2 is the natural logarithm of 2.

See also

Notes

  1. Chaos: An Introduction to Dynamical Systems, K.T. Alligood, T.D. Sauer, J.A. Yorke, Textbooks in mathematical sciences ,Springer, 1996, ISBN 978-0-38794-677-1
  2. Non-Linear Ordinary Differential Equations: Introduction for Scientists and Engineers (4th Edition), D.W. Jordan, P. Smith, Oxford University Press, 2007, ISBN 978-0-19-920825-8
  3. Alligood, p. 503.
  4. Alligood, p. 504.
  5. 5.0 5.1 Nonlinear Dynamics and Chaos, Steven H. Strogatz, Studies in Nonlinearity ,Perseus Books Publishing, 1994, ISBN 978-0-7382-0453-6
  6. Lua error in package.lua at line 80: module 'strict' not found.
  7. Lua error in package.lua at line 80: module 'strict' not found.
  8. Lua error in package.lua at line 80: module 'strict' not found.
  9. Lua error in package.lua at line 80: module 'strict' not found.
  10. Lua error in package.lua at line 80: module 'strict' not found.

References

  • Alligood, Kathleen T., Tim D. Sauer, James A. Yorke, Chaos: An Introduction to Dynamical Systems, Textbooks in mathematical sciences Springer, 1996, ISBN 978-0-38794-677-1
  • Lua error in package.lua at line 80: module 'strict' not found.
  • Lua error in package.lua at line 80: module 'strict' not found.
  • Lua error in package.lua at line 80: module 'strict' not found.

External links