Fundamental solution matrix and Cauchy properties of quaternion combined impulsive matrix dynamic equation on time scales

In this paper, we establish some basic results for quaternion combined impulsive matrix dynamic equation on time scales for the first time. Quaternion matrix combined-exponential function is introduced and some basic properties are obtained. Based on this, the fundamental solution matrix and corresponding Cauchy matrix for a class of quaternion matrix dynamic equation with combined derivatives and bi-directional impulses are derived.

On representation formulas for solutions of linear differential equations with Caputo fractional derivatives

In the paper, a linear differential equation with variable coefficients and a Caputo fractional derivative is considered. For this equation, a Cauchy problem is studied, when an initial condition is given at an intermediate point that does not necessarily coincide with the initial point of the fractional differential operator. A detailed analysis of basic properties of the fundamental solution matrix is carried out. In particular, the Hölder continuity of this matrix with respect to both variables is proved, and its dual definition is given. Based on this, two representation formulas for the solution of the Cauchy problem are proposed and justified.

Stability Assessment of Stochastic Differential-Algebraic Systems via Lyapunov Exponents with an Application to Power Systems

In this paper, we discuss stochastic differential-algebraic equations (SDAEs) and the asymptotic stability assessment for such systems via Lyapunov exponents (LEs). We focus on index-1 SDAEs and their reformulation as ordinary stochastic differential equations (SDEs). Via ergodic theory, it is then feasible to analyze the LEs via the random dynamical system generated by the underlying SDEs. Once the existence of well-defined LEs is guaranteed, we proceed to the use of numerical simulation techniques to determine the LEs numerically. Discrete and continuous QR decomposition-based numerical methods are implemented to compute the fundamental solution matrix and use it in the computation of the LEs. Important computational features of both methods are illustrated via numerical tests. Finally, the methods are applied to two applications from power systems engineering, including the single-machine infinite-bus (SMIB) power system model.

LYAPUNOV'S EXPONENTS FOR NONSMOOTH DYNAMICS WITH IMPACTS: STABILITY ANALYSIS OF THE ROCKING BLOCK

The plane dynamics of a rigid block simply supported on a harmonically moving rigid ground is a problem which still needs investigating, although the matter has been the subject of much research since the last century. Unilateral contacts, Coulomb friction and impacts make the system hybrid as it reveals a mixed continuous and discontinuous nature. Thus, stability analysis requires the extension and adaptation of concepts with regard to variational-perturbative procedures. In particular, discontinuous systems exhibit discontinuities or "saltations" in the fundamental solution matrix which must be analyzed carefully. In this paper, the adaptation of numerical methods that permit us to obtain characteristic multipliers and Lyapunov's exponents for the rocking mode of the block will be tackled. Analytical methods are used for the linearized equations of motion; the results are compared with those in the scientific literature.

Chebyshev Expansion of Linear and Piecewise Linear Dynamic Systems With Time Delay and Periodic Coefficients Under Control Excitations

In this paper, a new efficient method is proposed to obtain the transient response of linear or piecewise linear dynamic systems with time delay and periodic coefficients under arbitrary control excitations via Chebyshev polynomial expansion. Since the time domain can be divided into intervals with length equal to the delay period, at each such interval the fundamental solution matrix for the corresponding periodic ordinary differential equation (without delay) is constructed in terms of shifted Chebyshev polynomials by using a previous technique that reduces the problem to a set of linear algebraic equations. By employing a convolution integral formula, the solution for each interval can be directly obtained in terms of the fundamental solution matrix. In addition, by combining the properties of the periodic system and Floquet theory, the computational processes are simplified and become very efficient. An alternate version, which does not employ Floquet theory, is also presented. Several examples of time-periodic delay systems, when the excitation period is equal to or larger than the delay period and for linear and piecewise linear systems, are studied. The numerical results obtained via this method are compared with those obtained from Matlab DDE23 software (Shampine, L. F., and Thompson, S., 2001, “Solving DDEs in MATLAB,” Appl. Numer. Math., 37(4), pp. 441–458.) An error bound analysis is also included. It is found that this method efficiently provides accurate results that find general application in areas such as machine tool vibrations and parametric control of robotic systems.

Fold Bifurcations in Discontinuous Systems

This paper treats discontinuous fold bifurcations of periodic solutions of discontinuous systems. It is shown how jumps in the fundamental solution matrix lead to jumps of the Floquet multipliers of periodic solutions. A Floquet multiplier of a discontinuous system can jump through the unit circle causing a discontinuous bifurcation. Numerical examples are treated which show discontinuous fold bifurcations. The discontinuous fold bifurcation can connect stable branches to branches with infinitely unstable solutions.