himoto / bifurcation

Numerical Bifurcation & Stability Analysis of ODE Models

Home Page:https://github.com/biomass-dev/BioMASS.jl

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Bifurcation Analysis

A numerical study of the changes in the dynamics and stability of a system upon variations in its parameters


Procedure for stability analysis at fixed points

Consider the following system of ordinary differential equations: dx/dt = F(x)

  1. Determine the fixed point vector, x, solving F(x) = 0
  2. Construct the Jacobian matrix, J(x) = F(x)/x
  3. Compute eigenvalues of J(x): det |J(x) − λE| = 0
  4. Conclude on stability or instability of x based on the real parts of eigenvalues
    • All eigenvalues have real parts less than zero → x is stable
    • At least one of the eigenvalues has a real part greater than zero → x is unstable

Requirements

Plotting bifurcation diagram

# Change directory to examples/XXX
include("bifurcation.jl")
using .Bifurcation

Bifurcation.analysis()
Bifurcation.diagram()

Installation

$ git clone https://github.com/himoto/bifurcation.git

Acknowledgements

I would particularly like to thank Dr. Gouhei Tanaka (Graduate School of Engineering, The University of Tokyo) for valuable discussions.

References

License

MIT

About

Numerical Bifurcation & Stability Analysis of ODE Models

https://github.com/biomass-dev/BioMASS.jl

License:MIT License


Languages

Language:Julia 96.6%Language:Shell 3.4%