This paper deals with the superharmonic resonance of a ferromagnetic thin rectangular plate in the air-gap magnetic field excited by armature magnetic potential. Electromagnetic forces applied on the plate by the air-gap magnetic field cause the plate to transversal vibrate, which affects the air-gap magnetic field in turn, resulting the vibration and magnetic field coupling. According to the basic electromagnetic field theory and considering the magnetoelastic coupling effect, the air-gap magnetic field intensity is obtained by solving the Laplace's equation satisfied the air-gap magnetic boundary conditions. The electromagnetic force model of soft ferromagnetic plates is determined based on theories of electromagnetic and elasticity. According to the large deflection theory of plates, basic energy relationships and variational equations of the elastic plate are given. Eventually, the nonlinear magnetoelastic vibration equation of ferromagnetic thin plates is derived using Hamilton's principle and Galerkin method. The multi-scale method is used to solve the superharmonic resonance to obtain the amplitude-frequency response equation and the stability discriminant of solutions. The topological analysis of amplitude-frequency equation is carried out using singularity theory, and the bifurcation characteristics of systems on physical parameter planes in different regions are obtained according to the transition set. The correctness of analytical solutions is verified by comparison with numerical solutions. Through numerical calculations, curves of the static deflection and equivalent magnetic force of plate with parameters are given, and the amplitude curves, dynamic phase plane trajectories and time history diagrams of system response with changes in electromagnetic and structural parameters are plotted. Results show that both the decrease of armature magnetic potential amplitude and the increases of plate thickness and initial air-gap thickness reduce the static deflection. The increase of armature magnetic potential amplitude increases the equivalent magnetic force. As the decrease of initial air-gap thickness and the increases of armature magnetic potential amplitude and excitation force amplitude, the amplitudes of the upper branch and lower branch curves representing stable solutions decrease and increase, respectively, and the single-value solution region increases.