For multi-particle coupled systems, the effects of environmental fluctuations on each particle are often different in actual situations. To this end, this paper studies the collective dynamic behaviors in globally coupled harmonic oscillators driven by different frequency fluctuations, including synchronization, stability and stochastic resonance (SR). The statistical synchronicity between particles' behaviors is derived by reasonably grouping variables and using random average method, and then the statistical equivalence between behaviors of mean field and behaviors of single particle is obtained. Therefore, the characteristics of mean field's behaviors (that is, collective behaviors) can be obtained by studying behaviors of any single particle. Moreover, the output amplitude gain and the necessary and sufficient condition for the system stability are obtained by using this synchronization. The former lays a theoretical foundation for analyzing the stochastic resonance behavior of the system, and the latter gives the scope of adaptation of the conclusions in this paper. In terms of numerical simulation, the research is mainly carried out through the stochastic Taylor expansion algorithm. Firstly, the influence of system size
Nand coupling strength
$\varepsilon$
on the stability area and synchronization time is analyzed. The results show that with the increase of the coupling strength
$\varepsilon$
or the increase of the system size
N, the coupling force between particles increases, and the orderliness of the system increases, so that the stable region gradually increases and the synchronization time gradually decreases. Secondly, the stochastic resonance behavior of the system is studied. Noises provide randomness for the system, and coupling forces provide orderliness for the system. The two compete with each other, so that the system outputs about the noise intensity
$\sigma$
, the coupling strength
$\varepsilon$
and the system size
Nexhibit stochastic resonance behavior. As the coupling strength increases or the system size increases, the orderliness of the system increases, and greater noise intensity is required to provide stronger randomness to achieve optimal matching with it, so as to the resonance of the noise intensity
$\sigma$
, the peak gradually shifts to the right. Conversely, as the noise intensity
$\sigma$
increases, the resonance peak of the coupling strength
$\varepsilon$
and the system size
Nwill also shift to the right.