文摘
The determination of the stability of systems with time delays is of high importance in many industrial and research applications. In this study, we improve the numerical integration method (NIM) by using the Lagrange form interpolating polynomial to approximate the delayed terms and construct a periodic discrete dynamical map for the damped Mathieu equation with time delays. Hence, we can obtain the Floquet transition matrices without matrix multiplication, which can reduce calculation time when the matrix inversed has a small bandwidth according to the computation based on MATLAB. To compare the calculation efficiency and computational accuracy between the updated and original methods, we compare the stability charts calculated by using three NIM methods for a damped Mathieu equation with multiple delays. To further confirm the efficiency of the presented method, we calculate the stability of the Mathieu equation with distributed delays and time-periodic delays, and then we compare the computational accuracy and calculation time with those of the semi-discretization method.