A new identification algorithm for the Preisach model is presented. The algorithm treats the identification procedure of the Preisach model as an inverse problem where the independent variables are parameters of the distribution function and the objective function is constructed using only the initial magnetization curve or only tile major loop of the hysteresis curve as well as the whole reversal curves. To parameterize the distribution function, the Bezier spline and Gaussian function are used for the coercive and interaction fields axes, respectively. The presented algorithm is applied to the ferrite permanent magnets, and the distribution functions are correctly found from the major loop of the hysteresis curve or the initial magnetization curve.