In this work, we develop a homotopy continuation method to find the characteristic roots of delay differential equations with multiple delays. We introduce a homotopy parameter μ into the characteristic equation in such a way that for μ = 0 this equation contains only one exponential term (corresponding to the largest delay) and for μ = 1 the original characteristic equation is recovered. For μ = 0, all the characteristic roots can be expressed in terms of the Lambert W function. Pseudo-arclength continuation is then used to trace the roots as a function of μ. We demonstrate the method on several test cases. Cases where it may fail are also discussed.
BortzDM (2015) Eigenvalues for a two-lag linear delay differential equation. IFAC-PapersOnLine48(12): 13–16.
4.
Cepeda-GomezRMichielsW (2015) Some special cases in the stability analysis of multi-dimensional time-delay systems using the matrix Lambert W function. Automatica53: 339–345.
5.
ChowSNMallet-ParetJYorkeJA (1978) Finding zeroes of maps: homotopy methods that are constructive with probability one. Mathematics of Computation32(143): 887–899.
6.
CorlessRMGonnetGHHareDEJeffreyDJKnuthDE (1996) On the Lambert W function. Advances in Computational mathematics5(1): 329–359.
7.
ErgencAFOlgacNFazeliniaH (2007) Extended Kronecker summation for cluster treatment of LTI systems with multiple delays. SIAM Journal on Control and Optimization46(1): 143–155.
8.
GargNKMannBPKimNHKurdiMH (2007) Stability of a time-delayed system with parametric excitation. Journal of Dynamic Systems, Measurement, and Control129(2): 125–135.
9.
GuKChenJKharitonovVL (2003) Stability of time-delay systems, Springer Science & Business Media.
10.
InspergerTStépánG (2002) Semi-discretization method for delayed systems. International Journal for Numerical Methods in Engineering55(5): 503–518.
11.
InspergerTStépánGBaylyPMannB (2003) Multiple chatter frequencies in milling processes. Journal of Sound and Vibration262(2): 333–345.
12.
Kalmár-NagyTStépánGMoonFC (2001) Subcritical Hopf bifurcation in the delay equation model for machine tool vibrations. Nonlinear Dynamics26(2): 121–142.
13.
KaneDMShoreKA (2005) Unlocking dynamical diversity: optical feedback effects on semiconductor lasers, New York: John Wiley & Sons.
14.
KarimiHR (2011) Robust delay-dependent control of uncertain time-delay systems with mixed neutral, discrete, and distributed time-delays and Markovian switching parameters. IEEE Transactions on Circuits and Systems I: Regular Papers58(8): 1910–1923.
15.
LiTY (1997) Numerical solution of multivariate polynomial systems by homotopy continuation methods. Acta numerica6: 399–436.
16.
MeijerHKalmár-NagyT (2012) The Hopf–van der Pol System: failure of a homotopy method. Differential Equations and Dynamical Systems20(3): 323–328.
17.
NayfehAHBalachandranB (2008) Applied nonlinear dynamics: analytical, computational and experimental methods, New York: John Wiley & Sons.
18.
NiculescuSI (2001) Delay effects on stability: a robust control approachVolume 269.Springer Science & Business Media.
19.
OlgacNSipahiR (2002) An exact method for the stability analysis of time-delayed linear time-invariant (LTI) systems. IEEE Transactions on Automatic Control47(5): 793–797.
20.
OlgacNSipahiR (2004a) The cluster treatment of characteristic roots and the neutral type time-delayed systems. In: ASME 2004 International Mechanical Engineering Congress and Exposition, American Society of Mechanical Engineers, pp. 1359–1368.
21.
OlgacNSipahiR (2004b) A practical method for analyzing the stability of neutral type LTI-time delayed systems. Automatica40(5): 847–853.
22.
OroszGStépánG (2004) Hopf bifurcation calculations in delayed systems with translational symmetry. Journal of Nonlinear Science14(6): 505–528.
23.
OzoegwuCG (2016) Estimating stability boundaries of distributed-delay oscillators via two-stage numerical integration. Journal of Vibration and Control. DOI: 10.1177/1077546316668466.
24.
Pitcher AB and Corless RM (2005) Quasipolynomial root-finding: A numerical homotopy method. In: Proceedings of the Canadian Undergraduate Mathematics Conference.
25.
RodriguesSBartonDSzalaiRBenjaminORichardsonMPTerryJR (2009) Transitions to spike-wave oscillations and epileptic dynamics in a human cortico-thalamic mean-field model. Journal of Computational Neuroscience27(3): 507–526.
26.
SalleseLGrossiNScippaACampatelliG (2016) Numerical investigation of chatter suppression in milling using active fixtures in open-loop control. Journal of Vibration and Control. DOI: 10.1177/1077546316668686.
27.
SchreiberHMeerKSchmittB (2002) Dimensional synthesis of planar Stephenson mechanisms for motion generation using circlepoint search and homotopy methods. Mechanism and Machine Theory37(7): 717–737.
28.
SieberJKrauskopfB (2008) Control based bifurcation analysis for experiments. Nonlinear Dynamics51(3): 365–377.
29.
SommeseAJVerscheldeJWamplerCW (2004) Advances in polynomial continuation for solving problems in kinematics. Journal of Mechanical Design126(2): 262–268.
30.
SunJSongB (2009) Control studies of time-delayed dynamical systems with the method of continuous time approximation. Communications in Nonlinear Science and Numerical Simulation14(11): 3933–3944.
31.
TariHSuJJLiTY (2010) A constrained homotopy technique for excluding unwanted solutions from polynomial equations arising in kinematics problems. Mechanism and Machine Theory45(6): 898–910.
32.
VyasarayaniCP (2012) Galerkin approximations for higher order delay differential equations. Journal of Computational and Nonlinear Dynamics7(3): 031004.
33.
VyasarayaniCPSubhashSKalmár-NagyT (2014) Spectral approximations for characteristic roots of delay differential equations. International Journal of Dynamics and Control2(2): 126–132.
34.
VyasarayaniCPUchidaTMcPheeJ (2012) Nonlinear parameter identification in multibody systems using homotopy continuation. Journal of Computational and Nonlinear Dynamics7(1): 011012.
35.
VyhlidalTZitekP (2009) Mapping based algorithm for large-scale computation of quasi-polynomial zeros. IEEE Transactions on Automatic Control54(1): 171–177.
36.
WatsonLTBillupsSCMorganAP (1987) Algorithm 652: HOMPACK: A suite of codes for globally convergent homotopy algorithms. ACM Transactions on Mathematical Software (TOMS), 13(3, 281–310.
37.
Yi S, Yu S and Kim JH (2011) Analysis of neural networks with time-delays using the Lambert W function. In: American Control Conference (ACC), 2011. IEEE, pp. 3221–3226.
