The invention provides a numerical simulation analysis method for multi-field coupling of a loudspeaker. The method includes: firstly, establishing a finite element model, to be more specific, (1) establishing a physical model, namely adding corresponding physical fields and adding a solving step according to different structures of the loudspeaker and different physical principles referred in different working stages, (2) establishing a geometric model, (3) defining material properties, modifying the physical fields and defining boundary conditions, load and the like, (4) performing meshing on the model; secondly, performing solving and post-processing, to be more specific, solving the established model via a finite element method, and subjecting results to post-processing to finally obtain loudspeaker characteristic parameters under mutual influences among the different physical fields, wherein the loudspeaker characteristic parameters include loudspeaker electromagnetic characteristics such as Bl(x) and Z(f), structural mechanical characteristics such as f0 and K(x) of a loudspeaker vibration system, loudspeaker acoustic characteristics of SPL, directivity curves and the like and temperature field characteristics and the like of the loudspeaker. In the process of establishing the finite element model, the material properties of variation with temperature and the special boundary conditions and load modes are defined to couple the different physical fields.