Yu. A. Kuznetsov

Bio: Yu. A. Kuznetsov is an academic researcher from Utrecht University. The author has contributed to research in topics: Numerical continuation & Saddle-node bifurcation. The author has an hindex of 17, co-authored 27 publications receiving 1988 citations. Previous affiliations of Yu. A. Kuznetsov include University of Twente & Polytechnic University of Milan.

MATCONT: A MATLAB package for numerical bifurcation analysis of ODEs

Abstract: MATCONT is a graphical MATLAB software package for the interactive numerical study of dynamical systems. It allows one to compute curves of equilibria, limit points, Hopf points, limit cycles, period doubling bifurcation points of limit cycles, and fold bifurcation points of limit cycles. All curves are computed by the same function that implements a prediction-correction continuation algorithm based on the Moore-Penrose matrix pseudo-inverse. The continuation of bifurcation points of equilibria and limit cycles is based on bordering methods and minimally extended systems. Hence no additional unknowns such as singular vectors and eigenvectors are used and no artificial sparsity in the systems is created. The sparsity of the discretized systems for the computation of limit cycles and their bifurcation points is exploited by using the standard Matlab sparse matrix methods. The MATLAB environment makes the standard MATLAB Ordinary Differential Equations (ODE) Suite interactively available and provides computational and visualization tools; it also eliminates the compilation stage and so makes installation straightforward. Compared to other packages such as AUTO and CONTENT, adding a new type of curves is easy in the MATLAB environment. We illustrate this by a detailed description of the limit point curve type.

Cl_matcont: a continuation toolbox in Matlab

Abstract: CL_MATCONT is a Matlab continuation package for the numerical study of a range of parameterized nonlinear problems. In the case of ODEs it allows to compute curves of equilibria, limit point, Hopf points, limit cycles and period doubling bifurcation points of limit cycles. All curves are computed by the same function that implements a prediction-correction continuation algorithm based on the Moore - Penrose matrix pseudo-inverse. The continuation of bifurcation points of equilibria and limit cycles is based on bordering methods and minimally extended systems. Hence no additional unknowns such as singular vectors and eigenvectors are used and no artificial sparsity in the systems is created.The inherent sparsity of the discretized systems for the computation of limit cycles and their bifurcation points is exploited by using the standard Matlab sparse matrix methods.CL_MATCONT furthermore allows to compute solution branches to underdetermined systems of nonlinear equations and parameterized boundary value problems.

Two-parameter discontinuity-induced bifurcations of limit cycles: classification and open problems

Abstract: This paper proposes a strategy for the classification of codimension-two discontinuity-induced bifurcations of limit cycles in piecewise smooth systems of ordinary differential equations. Such nonsmooth transitions (also known as C-bifurcations) occur when the cycle interacts with a discontinuity boundary of phase space in a nongeneric way, such as grazing contact. Several such codimension-one events have recently been identified, causing for example, period-adding or sudden onset of chaos. Here, the focus is on codimension-two grazings that are local in the sense that the dynamics can be fully described by an appropriate Poincare map from a neighborhood of the grazing point (or points) of the critical cycle to itself. It is proposed that codimension-two grazing bifurcations can be divided into three distinct types: either the grazing point is degenerate, or the grazing cycle is itself degenerate (e.g. nonhyperbolic) or we have the simultaneous occurrence of two grazing events. A careful distinction is drawn between their occurrence in systems with discontinuous states, discontinuous vector fields, or that with discontinuity in some derivative of the vector field. Examples of each kind of bifurcation are presented, mostly derived from mechanical applications. For each example, where possible, principal bifurcation curves characteristic to the codimension-two scenario are presented and general features of the dynamics discussed. Many avenues for future research are opened.

Numerical Methods for Two-Parameter Local Bifurcation Analysis of Maps

Abstract: We discuss new and improved algorithms for the bifurcation analysis of fixed points and periodic orbits (cycles) of maps and their implementation in matcont, a MATLAB toolbox for continuation and bifurcation analysis of dynamical systems. This includes the numerical continuation of fixed points of iterates of the map with one control parameter, detecting and locating their bifurcation points (i.e., limit point, period-doubling, and Neimark-Sacker) and their continuation in two control parameters, as well as detection and location of all codimension 2 bifurcation points on the corresponding curves. For all bifurcations of codim 1 and 2, the critical normal form coefficients are computed, both numerically with finite directional differences and using symbolic derivatives of the original map. Using a parameter-dependent center manifold reduction, explicit asymptotics are derived for bifurcation curves of double and quadruple period cycles rooted at codim 2 points of cycles with arbitrary period. These asymptotics are implemented into the software and allow one to switch at codim 2 points to the continuation of the double and quadruple period bifurcations. We provide two examples illustrating the developed techniques: a generalized Henon map and a juvenile/adult competition model from mathematical biology.

Numerical Normalization Techniques for All Codim 2 Bifurcations of Equilibria in ODE's

Abstract: Explicit computational formulas for the coefficients of the normal forms for all codim 2 equilibrium bifurcations of equilibria in autonomous ODEs are derived. They include second-order coefficients for the Bogdanov--Takens bifurcation, third-order coefficients for the cusp and fold-Hopf bifurcations, and coefficients of the fifth-order terms for the generalized Hopf (Bautin) and double Hopf bifurcations. The formulas are independent on the dimension of the phase space and involve only critical eigenvectors of the Jacobian matrix of the right-hand sides and its transpose, as well as multilinear functions from the Taylor expansion of the right-hand sides at the critical equilibrium. The normal form coefficients for the fold-Hopf bifurcation in the "new" Lorenz model are computed using the derived formulas, proving the existence of a nontrivial invariant set in the system.

MATCONT: A MATLAB package for numerical bifurcation analysis of ODEs

Abstract: MATCONT is a graphical MATLAB software package for the interactive numerical study of dynamical systems. It allows one to compute curves of equilibria, limit points, Hopf points, limit cycles, period doubling bifurcation points of limit cycles, and fold bifurcation points of limit cycles. All curves are computed by the same function that implements a prediction-correction continuation algorithm based on the Moore-Penrose matrix pseudo-inverse. The continuation of bifurcation points of equilibria and limit cycles is based on bordering methods and minimally extended systems. Hence no additional unknowns such as singular vectors and eigenvectors are used and no artificial sparsity in the systems is created. The sparsity of the discretized systems for the computation of limit cycles and their bifurcation points is exploited by using the standard Matlab sparse matrix methods. The MATLAB environment makes the standard MATLAB Ordinary Differential Equations (ODE) Suite interactively available and provides computational and visualization tools; it also eliminates the compilation stage and so makes installation straightforward. Compared to other packages such as AUTO and CONTENT, adding a new type of curves is easy in the MATLAB environment. We illustrate this by a detailed description of the limit point curve type.

Mixed-Mode Oscillations with Multiple Time Scales

Abstract: Mixed-mode oscillations (MMOs) are trajectories of a dynamical system in which there is an alternation between oscillations of distinct large and small amplitudes. MMOs have been observed and studied for over thirty years in chemical, physical, and biological systems. Few attempts have been made thus far to classify different patterns of MMOs, in contrast to the classification of the related phenomena of bursting oscillations. This paper gives a survey of different types of MMOs, concentrating its analysis on MMOs whose small-amplitude oscillations are produced by a local, multiple-time-scale “mechanism.” Recent work gives substantially improved insight into the mathematical properties of these mechanisms. In this survey, we unify diverse observations about MMOs and establish a systematic framework for studying their properties. Numerical methods for computing different types of invariant manifolds and their intersections are an important aspect of the analysis described in this paper.

New features of the software MatCont for bifurcation analysis of dynamical systems

Abstract: Bifurcation software is an essential tool in the study of dynamical systems. From the beginning (the first packages were written in the 1970's) it was also used in the modelling process, in particular to determine the values of critical parameters. More recently, it is used in a systematic way in the design of dynamical models and to determine which parameters are relevant. MatCont and Cl_MatCont are freely available matlab numerical continuation packages for the interactive study of dynamical systems and bifurcations. MatCont is the GUI-version, Cl_MatCont is the command-line version. The work started in 2000 and the first publications appeared in 2003. Since that time many new functionalities were added. Some of these are fairly simple but were never before implemented in continuation codes, e.g. Poincare maps. Others were only available as toolboxes that can be used by experts, e.g. continuation of homoclinic orbits. Several others were never implemented at all, such as periodic normal forms for codimension 1 bifurcations of limit cycles, normal forms for codimension 2 bifurcations of equilibria, detection of codimension 2 bifurcations of limit cycles, automatic computation of phase response curves and their derivatives, continuation of branch points of equilibria and limit cycles. New numerical algorithms for these computations have been published or will appear elsewhere; here we restrict to their software implementation. We further discuss software issues that are in practice important for many users, e.g. how to define a new system starting from an existing one, how to import and export data, system descriptions, and computed results.

Modelling the fear effect in predator–prey interactions

Abstract: A recent field manipulation on a terrestrial vertebrate showed that the fear of predators alone altered anti-predator defences to such an extent that it greatly reduced the reproduction of prey. Because fear can evidently affect the populations of terrestrial vertebrates, we proposed a predator-prey model incorporating the cost of fear into prey reproduction. Our mathematical analyses show that high levels of fear (or equivalently strong anti-predator responses) can stabilize the predator-prey system by excluding the existence of periodic solutions. However, relatively low levels of fear can induce multiple limit cycles via subcritical Hopf bifurcations, leading to a bi-stability phenomenon. Compared to classic predator-prey models which ignore the cost of fear where Hopf bifurcations are typically supercritical, Hopf bifurcations in our model can be both supercritical and subcritical by choosing different sets of parameters. We conducted numerical simulations to explore the relationships between fear effects and other biologically related parameters (e.g. birth/death rate of adult prey), which further demonstrate the impact that fear can have in predator-prey interactions. For example, we found that under the conditions of a Hopf bifurcation, an increase in the level of fear may alter the direction of Hopf bifurcation from supercritical to subcritical when the birth rate of prey increases accordingly. Our simulations also show that the prey is less sensitive in perceiving predation risk with increasing birth rate of prey or increasing death rate of predators, but demonstrate that animals will mount stronger anti-predator defences as the attack rate of predators increases.