On Stability of LMS Methods and Characteristic Roots of Delay Differential Equations

We investigate the use of linear multistep (LMS) methods for computing characteristic roots of systems of (linear) delay differential equations (DDEs) with multiple fixed discrete delays. These roots are important in the context of stability and bifurcation analysis. We prove convergence orders for the characteristic root approximations and analyze under what condition for the steplength the discrete integration scheme retains certain delay-independent stability properties of the original equations. Unlike existing results, we concentrate on the recovery of both stability and instability. We illustrate our findings with a number of numerical test results.

[1]  Ben P. Sommeijer,et al.  Stability in linear multistep methods for pure delay equations , 1983 .

[2]  G. Strang,et al.  THE OPTIMAL ACCURACY OF DIFFERENCE SCHEMES , 1983 .

[3]  Toshiyuki Koto,et al.  A stability property ofA-stable natural Runge-Kutta methods for systems of delay differential equations , 1994 .

[4]  J. H. van Lint,et al.  Functions of one complex variable II , 1997 .

[5]  Dirk Roose,et al.  Numerical stability analysis and computation of Hopf bifurcation points for delay differential equations , 1996 .

[6]  Guang-Da Hu,et al.  Estimation of numerically stable step-size for neutral delay-differential equations via spectral radius , 1997 .

[7]  Neville J. Ford,et al.  Numerical Hopf bifurcation for a class of delay differential equations , 2000 .

[8]  Stability , 1973 .

[9]  K. J. in 't Hout,et al.  On the stability of adaptations of Runge-Kutta methods to systems of delay differential equations , 1996 .

[10]  K. Meerbergen,et al.  Matrix transformations for computing rightmost eigenvalues of large sparse non-symmetric eigenvalue problems , 1996 .

[11]  D. S. Watanabe,et al.  THE STABILITY OF DIFFERENCE FORMULAS FOR DELAY DIFFERENTIAL EQUATIONS , 1985 .

[12]  J. Conway,et al.  Functions of a Complex Variable , 1964 .

[13]  Christian Lubich,et al.  Periodic orbits of delay differential equations under discretization , 1998 .

[14]  D. Sorensen Numerical methods for large eigenvalue problems , 2002, Acta Numerica.

[15]  J. Hale Theory of Functional Differential Equations , 1977 .

[16]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

[17]  Tian Hong-jiong,et al.  The numerical stability of linear multistep methods for delay differential equations with many delays , 1996 .

[18]  Sue Ann Campbell,et al.  Stability, Bifurcation, and Multistability in a System of Two Coupled Neurons with Multiple Time Delays , 2000, SIAM J. Appl. Math..

[19]  Christopher T. H. Baker,et al.  Issues in the numerical solution of evolutionary delay differential equations , 1995, Adv. Comput. Math..

[20]  Neville J. Ford,et al.  The Use of Boundary Locus Plots in the Identification of Bifurcation Points in Numerical Approximati , 1998 .

[21]  G. Hu,et al.  Stability of Runge-Kutta methods for delay differential systems with multiple delays , 1999 .

[22]  V. Kolmanovskii,et al.  Applied Theory of Functional Differential Equations , 1992 .

[23]  G. Strang Trigonometric Polynomials and Difference Methods of Maximum Accuracy , 1962 .

[24]  Neville J. Ford,et al.  How do numerical methods perform for delay differential equations undergoing a Hopf bifurcation , 2000 .

[25]  Dirk Roose,et al.  Numerical computation of stability and detection of Hopf bifurcations of steady state solutions of delay differential equations , 1999, Adv. Comput. Math..

[26]  M. N. Spijker Stiffness in numerical initial-value problems , 1996 .

[27]  William L. Kath,et al.  Bifurcation with memory , 1986 .