In this paper, we consider numerical solutions for a general form of fractional delay differential equations (FDDEs) with fractional derivatives defined in the Caputo sense. A fractional integration operational matrix, created using a fractional Taylor basis, is applied to solve these FDDEs. The main characteristic of this approach is, by utilizing the operational matrix of fractional integration, to reduce the given differential equation to a set of algebraic equations with unknown coefficients. This equation system can be solved efficiently using a computer algorithm. A bound on the error for the best approximation and fractional integration are also given. Several examples are given to illustrate the validity and applicability of the technique. The efficiency of the presented method is revealed by comparing results with some existing solutions, the findings of some other approaches from the literature and by plotting absolute error figures.