180 likes | 408 Views
N. t. K/2. Logistic Differential Equation. K. K. K/2. Logistic Differential Equation. d/dt N = - d/dN [r N 2 (N/3k -1/2)] = - d/dN V(N). V(N). K. N. 0. K/2. K. 3K/2. What is a bifurcation?. Max & Moritz, W. Busch. A qualitative change in the solution !.
E N D
N t K/2 Logistic Differential Equation K K
K/2 Logistic Differential Equation d/dt N = - d/dN [r N2 (N/3k -1/2)] = - d/dN V(N) V(N) K N 0 K/2 K 3K/2
What is a bifurcation? Max & Moritz, W. Busch A qualitative change in the solution !
Saddle Node Bifurcation (1-dim) Prototypical example:
Transcritical Bifurcatoin Prototypical example:
Pitchfork Bifurcation Prototypical example:
Hopf-Bifurcation Prototypical example:
Bifurcations ? Catastrophes ? • Use simple model from classical mechanics • Particle (m=1) in 1D-potential V(x,λ) • Potential changes very slowly with time • (slow time scale represented by λ) • Equations of motion: Stability of fix points:
Loss-of-Equilibrium (Saddle-Node Bifurcation) Potential, Lyapunov function Loss-of-equilibrium or catastrophe at λ = λc
Transcritical Bifurcation Potential, Lyapunov function Equilibrium branches exchange stability at λ=λc Other possibility: supercitical pitchfork bifurcation
Construct the dynamical System ! d/dt x = - d/dx V(x) V(x) = x2 (b - x2)
Dynamics of Two Dimensional Systems • Find the fixed points in the phase space! • Linearize the system about the fixed points! • Determine the eigenvalues of the Jacobian.
Liapunov Exponents • to determine the phase space extent • These exponents describe the mean exponential rate of divergence of two trajectories initially close to each other in phase space and are defined as: • where and are the initial (t=0) position and deviation, and is the deviation at time t. Such limit is proven to exist and is finite. • the largest exponent is the most important because it results from almost all deviations and therefore is the one found in practice when the initial deviation is chosen randomly.
Liapunov Exponents • numerical computation of the Liapunov exponents • the finite initial deviation may rapidly grow as large as the size of the orbits themselves -> scaling • Noting the deviation before the first rescaling, this is done by normalising to : every time ti