Approximate Stability Analysis and Computation of Solutions of Nonlinear Delay Differential Algebraic Equations with Time Periodic Coefficients

Approximate stability analysis of nonlinear delay differential algebraic equations (DDAEs) with periodic coefficients is proposed with a geometric interpretation of evolution of the linearized system. Firstly, a numerical algorithm based on direct integration by expansion in terms of Chebyshev polynomials is derived for linear analysis. The proposed algorithm is shown to have deeper connections with and be computationally less cumbersome than the solution of the underlying semi-explicit system via a similarity transformation. The stability of time periodic DDAE systems is characterized by the spectral radius of a “monodromy matrix”, which is a finite-dimensional approximation of a compact infinite-dimensional operator. The monodromy matrix is essentially a map of the Chebyshev coefficients (or collocation vector) of the state from the delay interval to the next adjacent interval of time. The computations are entirely performed with the original system to avoid cumbersome transformations associated with the semi-explicit form of the system. Next, two computational algorithms, the first based on perturbation series and the second based on Chebyshev spectral collocation, are detailed to obtain solutions of nonlinear DDAEs with periodic coefficients for consistent initial functions.