Received 6 September 2015; accepted 26 February 2016; published 29 February 2016
Fractional calculus attracted the attention of researchers because of its application in physics as the nonlinear oscillation of earthquake can be modeled with fractional derivatives  , and the fluid-dynamic traffic model with fractional derivatives  can eliminate the deficiency arising from the assumption of continuum traffic flow. Based on experimental data fractional, partial differential equations for seepage flow in porous media are suggested in  , and differential equations with fractional order have recently proved to be valuable tools to the modeling of many physical phenomena  . A review of some applications of fractional derivatives in continuum and statistical mechanics is given by Mainardi  . The analytic results on the existence and uniqueness of solutions of the fractional differential equations have been investigated by many authors see as an example  . During the last decades, several methods have been used to solve fractional differential equations, fractional partial differential equations, fractional integro-differential equations and dynamic systems containing fractional derivatives, such as Adomian’s decomposition method  , variational iteration method  -  , homotopy perturbation method  , homotopy analysis method  , spectral methods  , and other methods  .
Recently, Wang  studied the synchronized motions in a star network of coupled fractional order systems in which the major element is coupled to each of the non-interacting individual elements and Kilicman and Al- Zhour  studied several operational matrices for fractional integration and differentiation and expanded the Kronecker convolution product to the Riemann-Liouville fractional integral of matrices. Al-Zhour  introduced the exact solution of coupled fractional order systems by using Kronecker structure.
In the present paper, the exact solutions of coupled and uncoupled systems of matrix fractional differential equations for diagonal unknown matrices are presented by using a new attractive method and some illustrated examples are also given to show our new approach.
2. Basic Results and Preliminaries
In this section, we recall some basic results and definitions associated to Hadamard product, Mittage-Leffler function and Caputo fractional derivative that will be used to get our results later.
Definition 2.1. Let and.Then the Hadamard product of A and B is defined by  -  .
Definition 2.2. Let be a diagonal matrix. Then the diagonal extraction operator of A is defined by   .
Theorem 2.3. Let be diagonal matrices. Then
Definition 2.4. The one parameter Mittage-Leffler functions and Mittage-Leffler matrix functions of matrix are defined, respectively, for by  .
Note that the Mittage-Leffler matrix function of can be represented by using spectral decomposition method by  .
where and are the eigenvectors corresponding to the eigenvalue of A and AT, respectively.
Theorem 2.5. Let is a diagonal matrix and. Then 
Definition 2.6. The Caputo fractional derivative of with order, , is defined by  - .
Theorem 2.7. The relationship between the Mittage-Leffler function and Caputo derivative are given by:
3. Main Results
In this section, we present the general exact solutions of the coupled and uncoupled system of fractional differential equations for diagonal unknown matrices by using the using vector extraction operators and Hadamard product.
Lemma 3.1. Let be a given scalar matrix, be a given scalar vector, be a given vector function and be an unknown vector function to be solved. Then the exact solution of the following non-homogenous linear fractional system of order is given by  - .
is given by:
Theorem 3.2. Let and be given diagonal scalar matrices, be a given diagonal matrix function and be an unknown diagonal matrix function. Then the general vector extraction solution of the following non-homogeneous matrix fractional differential equation
is given by:
Proof. By using (2-3), then (3.3) can be represented by:
Hence, the vector extraction solution of (3.3) is given by:
Theorem 3.3. Let be given diagonal scalar matrices, and be an unknown diagonal matrix functions. Then the general vector extraction solution of the following general system of linear matrix fractional differential equations of order:
is given by:
Proof. By using (2-3), then (3.5) can be represented by:
Now by letting
. Then (3.7) can be written as:
Hence by using Lemma 3.1 and simple computations, then we get the solution as in (3-6).
Below we will discuss some important special cases of the general system as in Theorem 3.3.
Theorem 3.4. Let be given scalar diagonal matrices, be diagonal matrix functions, and be unknown diagonal matrix functions. Then the general solutions of the following coupled matrix fractional differential equations of order:
are given by:
Proof. By multiplying the second equation in (3-9) by, we get:
Then (3-9) can be written as
Now, by using of (3.12), then we get the following equivalent system:
Now by using (3-6), then the solution of (3.13) is given by:
Now we deal with
Now from (3-13), (3-18) and (3-19), we get
Then, we get the vector extraction solution as in (3-11).
Corollary 3.5. Let be given scalar diagonal matrices and be an unknown diagonal matrix functions. Then the general vector extraction solutions of the following coupled matrix fractional differential equations of order:
are given by:
Proof. The proof is straightforward by applying Theorem 3.4 by letting and by using the following fact:
4. Illustrated Examples
In the section, we give some illustrated examples to show our new approach as discussed in above section.
Example 4.1. Consider the following matrix linear fractional differential equation:
where and is diagonal matrix. Then the exact solution of (4-1) by applying Theorem 3.2 is given by:
Example 4.2. Consider the following system of order:
where. Then the system (4-2) can be rewritten as:
Now the exact solution of (4-3) by applying Theorem 3.2 is given by:
Example 4.3. Consider the following matrix fractional differential equation:
where and is diagonal matrix. Then the exact solution of (4-4) by applying Theorem 3.2 is given by:
Example 4.4. Consider the following matrix fractional differential equations of order:
where are diagonal matrices. So
Then the exact solution of (4-5) by applying Corollary 3.5 is given by:
Example 4.5. Consider the following coupled matrix fractional differential equations:
Then the exact solution by applying Corollary 3.5 is given by:
The general exact solutions of coupled system of matrix fractional differential equations with diagonal matrices coefficients by using vector extraction operators and Hadamard product in Caputo sense are presented with some illustrated examples. How to find the complexity of this method requires further research.
The authors express their sincere thanks to referees for very careful reading and helpful suggestion of this paper.
 He, J. (1998) Approximate Analytical Solution for Seepage Flow with Fractional Derivatives in Porous Media. Computer Methods in Applied Mechanics and Engineering, 167, 57-68.
 Mainardi, F. (1997) Fractional Calculus: Some Basic Problems in Continuum and Statistical Mechanics, Fractals and Fractional Calculus in Continuum Mechanics. Springer-Verlag, New York, 291-348.
 Momani, S. and Shawagfeh, N. (2006) Decomposition Method for Solving Fractional Riccati Differential Equations, Applied Mathematics and Computation, 182, 1083-1092.
 Momani, S. and Odibat, Z. (2006) Analytical Approach to Linear Fractional Partial Differential Equations Arising in Fluid Mechanics. Physics Letters A, 355, 271-279.
 Momani, S. and Odibat, Z. (2007) Homotopy Perturbation Method for Nonlinear Partial Differential Equations of Fractional Order. Physics Letters A, 365, 345-350.
 Odibat, Z. and Momani, S. (2008) Modified Homotopy Perturbation Method: Application to Quadratic Riccati Differential Equation of Fractional Order. Chaos, Solitons and Fractals, 36, 167-174.
 Odibat, Z. and Momani S. (2006) Application of Variational Iteration Method to Nonlinear Differential Equations of Fractional Order. International Journal of Nonlinear Sciences and Numerical Simulation, 7, 271-279.
 Hashim, I., Abdulaziz, O. and Momani, S. (2009) Homotopy Analysis Method for Fractional IVPs. Communications in Nonlinear Science and Numerical Simulation, 14, 674-684.
 Rawashdeh, E.A. (2006) Numerical Solution of Fractional Integro-Differential Equations by Collocation Method. Applied Mathematics and Computation, 176, 1-6.
 Wang, Q. (2006) Numerical Solutions for Fractional KdV-Burgers Equation by Adomian Decomposition Method. Applied Mathematics and Computation, 182, 1048-1055.
 Kilicman, A. and Al-Zhour, Z. (2007) Kronecker Operational Matrices for Fractional Calculus and Some Applications. Applied Mathematics and Computation, 187, 250-265.
 Al-Zhour, Z. (2014) The General (Vetor) Solutions of Such Linear (Coupled) Matrix Fractional Differential Equations by Using Kronecker Structures. Applied Mathematics and Computation, 232, 498-510.
 Balachandran, K., Kokila, J. and Trujillo, J. (2012) Relative Controllability of Fractional Dynamicsystems with Multiple Delays in Control. Applied Mathematics and Computation, 64, 3037-3045.
 Kiliçman, A. and Al-Zhour, Z. (2009) On the Connection between Kronecker and Hadamard Convolution Products of Matrices and Some Applications. Journal of Inequalities and Applications, 2009, Article ID: 736243.
 Al-Zhour, Z. (2008) A Computationally-Efficient Solutions of Coupled Matrix Differential Equations for Diagonal Unknown Matrices. Journal of Mathematical Sciences: Advances and Applications, 1, 373-387.
 Kiliçman, A. and Al-Zhour, Z. (2007) Vector Least-Squares Solutions of Coupled Singular Matrix Equations. Journal of Computational and Applied Mathematics, 206, 1051-1069.
 Al-Zhour, Z. and Kiliçman, A. (2007) Some New Connections between Matrix Products for Partitioned and Non-Partitioned Matrices. Computers & Mathematics with Applications, 54, 763-784.