“…Due to the spatial nonlinearity and the temporal nonhomogeneity, there is normally no analytic or numerical technique which can calculate the response of general stochastic resonant systems apart from bistable systems. For a single standard or biased overdamped bistable system, several theoretical methods, including two-state theory [14]- [16], linear response approximation [17]- [20], matrix continuation fraction [21]- [23], eigenfunction perturbation [24,25], and the method of moments [26]- [28], have been developed to explore SR within or beyond the linear response range. Among these methods, the linear response approximation and the method of moments have been extended to uncoupled bistable arrays [29] and mean-field coupled bistable arrays [30], respectively; however, the effect of correlation of the internal noise of the consisting elements has not been taken into account in these extensions.…”