Operators for the representation in the time-domain of systems with relaxation of Havriliak–Negami type are studied.
Systems of this type are encountered in physics, mechanical engineering, electric circuits, electromagnetism, etc.
It is provided an explicit representation in terms of fractional differences of Grunwald–Letnikov type.
The main properties of these operators are studied.
The discretization of the operators is also shown. The comparison with an alternative approach is numerically tested.