An efficient numerical method, the fractional Chebyshev collocation method, is proposed for obtaining the solution of a system of linear fractional order delaydifferential equations (FDDEs). It is shown that the proposed method overcomes several limitations of current numerical methods for solving linear FDDEs. For instance, the proposed method can be used for linear incommensurate order fractional differential equations and FDDEs, has spectral convergence (unlike finite differences), and does not require a canonical form. To accomplish this, a fractional differentiation matrix is derived at the Chebyshev-Gauss-Lobatto collocation points by using the discrete orthogonal relationship of the Chebyshev polynomials. Then, using two proposed discretization operators for matrix functions results in an explicit form of solution for a system of linear FDDEs with discrete delays. The advantages of using the fractional Chebyshev collocation method are demonstrated in two numerical examples in which the proposed method is compared with the Adams- Bashforth-Moulton method.