After in-depth investigation of literature, three possible approximations seem to be feasible:
This approximation includes the inflection point, i.e. it is very usefull to describe the virgin characteristic of a material.
However, calculating the parameters used in the approximation formula is a challenging nonlinear optimization problem.
This approximation neglects the inflection point, i.e. it is suitable to describe the center line of the (narrow) hysteresis loops of a soft magnetic material.
Calculating the parameters used in the approximation formula is a nonlinear optimization problem but much easier to solve.
Since the approximation formula describes the characteristic of magnetic flux density B dependent on magnetic field strength H, and in many cases the derivative of B is needed to calculate induced voltage, it is recommended to describe the derivative dB/dH additionally.
This approximation is mainly intended to describe roughly magnetization characteristic for educational purposes. It is sufficient to specify several points Bk(Hk) along the original characteristic. Especially for designing exam questions, it is easier to calculate the intersection between a linear interpolation segment and the linear charcteristic of a magnetic source than to determine the intersection between an nonlinar characteristic and the linear charcteristic of a magnetic source graphically.
Since the approximation formula describes the characteristic of magnetic flux density B dependent on magnetic field strength H, and in many cases the derivative of B is needed to calculate induced voltage, it is recommended to describe the derivative dB/dH additionally.