Exact solutions of fractional oscillation systems with pure delay

In this paper we study the exact solutions of a class of fractional delay differential equations. We consider the fractional derivative of the order between 1 and 2 in the sense of Caputo. In the first part, we introduce two novel matrix functions, namely, the generalized cosine-type and sine-type delay Mittag-Leffler matrix functions. Then we obtain the explicit solutions for the linear homogeneous equations subjecting to corresponding initial conditions, by means of undetermined coefficients. In the second part, we first obtain a particular solution by means of the Laplace transform for the inhomogeneous equations with null initial conditions. Then we give an analytical representation of the general solution of the inhomogeneous equations through the sum of its particular solution and the general solution of the corresponding homogeneous equation.