摘要:In this work, we derive the operational matrix using poly-Bernoulli polynomials. These polynomials generalize the Bernoulli polynomials using a generating function involving a polylogarithm function. We first show some new properties for these poly-Bernoulli polynomials; then we derive new operational matrix based on poly-Bernoulli polynomials for the Atangana–Baleanu derivative. A delay operational matrix based on poly-Bernoulli polynomials is derived. The error bound of this new method is shown. We applied this poly-Bernoulli operational matrix for solving fractional delay differential equations with variable coefficients. The numerical examples show that this method is easy to use and yet able to give accurate results.