The present article proposes a dynamical model to obtain ferroelectric hysteresis dynamics based on fractional derivatives. The consideration of a fractional derivative term widely increases the frequency bandwidth of the accuracy of the traditional hysteresis models. As a consequence, the model is suited for successfully taking into account the well-known scaling relations of the ferroelectric hysteresis area, < A > , versus the frequency, f , and field amplitude, E₀ . Under low frequency excitation, simulation tests provided good results regarding the comparison of the fractional model, experimental results and the well-known nonentire power law < A >∞f^1/3E₀^2/3 (where < A > represents the hysteresis loop area). These results were followed by comparing the hysteresis area obtained from the fractional model with that from the well known scaling relations as f→∞ , and the results were proposed as validation of the high frequency behavior. Next, the model was tested on large frequency bandwidths ( >6 decades) and validated with success using the comparison between simulation tests and the only experimental results available in literature obtained in such conditions by Liu etal [J. Phys.: Condens. Matter 16, 1189 (2004)] for BNT thin film samples.