In this article, a systematic numerical study is described of the effect of the polydispersity of suspensions of spherical particles on their dielectric behavior, in both the frequency and time domains, starting from the model proposed by DeLacey and White (J. Chem. Sec., Faraday Trans. 277, 2007 (1981)) for monodisperse suspensions. The distribution function of relaxation times, characterizing the dielectric response of the systems, is also calculated. It is found that in both the frequency and time domains the predicted behavior does not differ in any essential way from the one obtained for a monodisperse suspension with particle radius close to the volume-averaged mean radius of the polydisperse system. Hence, no arguments related to polydispersity seem to be useful for explaining the discrepancies frequently found between measured and calculated dielectric increments in suspensions, namely, those concerning the magnitude of the dielectric constant of the suspension (its low-frequency value), the value of the characteristic or relaxation frequency, or the overall shape of the relaxation pattern.