Stability of piecewise polynomial collocation for computing periodic solutions of delay differential equations

Summary. We prove numerical stability of a class of piecewise polynomial collocation methods on nonuniform meshes for computing asymptotically stable and unstable periodic solutions of the linear delay differential equation $\dot y(t) = a(t)y(t)+b(t)y(t-\tau) + f(t)$ by a (periodic) boundary value approach. This equation arises, e.g., in the study of the numerical stability of collocation methods for computing periodic solutions of nonlinear delay equations. We obtain convergence results for the standard collocation algorithm and for two variants. In particular, estimates of the difference between the collocation solution and the true solution are derived. For the standard collocation scheme the convergence results are “unconditional”, that is, they do not require mesh-ratio restrictions. Numerical results that support the theoretical findings are also given.

[1]  Kolmanovskii,et al.  Introduction to the Theory and Applications of Functional Differential Equations , 1999 .

[2]  Robert D. Russell,et al.  Numerical solution of boundary value problems for ordinary differential equations , 1995, Classics in applied mathematics.

[3]  Christopher A. H. Paul,et al.  A User-Guide to Archi - An Explicit Runge-Kutta Code for Solving Delay and Neutral Differential Equa , 1995 .

[4]  H. Keller,et al.  Approximation methods for nonlinear problems with application to two-point boundary value problems , 1975 .

[5]  E. Doedel Some stability theorems for finite difference collocation methods on nonuniform meshes , 1980 .

[6]  Robert D. Russell,et al.  Collocation Software for Boundary-Value ODEs , 1981, TOMS.

[7]  G. Bader,et al.  Solving Boundary Value Problems for Functional Differential Equations by Collation , 1985 .

[8]  Uri M. Ascher,et al.  Collocation Software for Boundary Value Differential-Algebraic Equations , 1994, SIAM J. Sci. Comput..

[9]  Alfredo Bellen,et al.  One-step collocation for delay differential equations , 1984 .

[10]  Dirk Roose,et al.  Collocation Methods for the Computation of Periodic Solutions of Delay Differential Equations , 2000, SIAM J. Sci. Comput..

[11]  D. Sarafyan,et al.  DKLAG6: a code based on continuously imbedded sixth-order Runge-Kutta methods for the solution of state-dependent functional differential equations , 1997 .

[12]  Thomas F. Fairgrieve,et al.  AUTO 2000 : CONTINUATION AND BIFURCATION SOFTWARE FOR ORDINARY DIFFERENTIAL EQUATIONS (with HomCont) , 1997 .

[13]  H. Keller,et al.  Analysis of Numerical Methods , 1969 .

[14]  K. J. in 't Hout,et al.  A new interpolation procedure for adapting Runge-Kutta methods to delay differential equations , 1992 .

[15]  Yunkang Liu,et al.  Runge–Kutta–collocation methods for systems of functional–differential and functional equations , 1999, Adv. Comput. Math..

[16]  K. J. In 't Hout Convergence of Runge-Kutta Methods for Delay Differential Equations , 2001 .

[17]  Fotios Giannakopoulos,et al.  Local and Global Hopf Bifurcation in a Scalar Delay Differential Equation , 1999 .

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

[19]  Dirk Roose,et al.  COMPUTATION, CONTINUATION AND BIFURCATION ANALYSIS OF PERIODIC SOLUTIONS OF DELAY DIFFERENTIAL EQUATIONS , 1997 .

[20]  C. Travis,et al.  Approximation methods for boundary value problems of differential equations with functional arguments , 1974 .

[21]  H. Stech A numerical analysis of the structure of periodic orbits in autonomous functional differential equations , 1987 .

[22]  C. D. Boor,et al.  Collocation at Gaussian Points , 1973 .