Magnetization of permanent magnet in the strong magnetic field is changeable. So the change of magnetization must be considered for the accurate analysis of the system with permanent magnets. But the hysteresis characteristice of permanent magnet cannot be represented in simple mathematical form. In this paper, Preisach model combined with finite element method which can describe the hysteresis phenomena is applied to analyze the permanent magnet system. To validate the method, it is applied to the force calculations between two magnets with different coercivities and the numerical results are compared with measured data.