The reconstruction of broadband sound sources is an important issue in industrial acoustics. In this paper, a model comprising multiple incoherent Gaussian random sources is considered. The aim is to estimate locations and powers of the sound sources using the pressures measured by an array of microphones. Each measured pressure is interpreted as a mixture of latent signals emitted by different sound sources. Then, an Iterative Beamforming (IB) method is developed to estimate the source parameters. This approach is based on the Expectation-Maximization (EM) algorithm, a well-known iterative procedure for solving maximum likelihood parameter estimation. More specifically, IB iteratively estimates the source contributions and performs beamforming on these estimates. In this work, experiments on real data illustrate the advantage of IB with respect to classical beamforming and Near-field Acoustical Holography (NAH). In particular, the proposed method is shown to work over a wider range of frequencies, to better estimate the source locations, and is able to quantify the powers of the sources. Furthermore, experiments illustrate that IB can not only localize the sources on a given surface, but also accurately estimate their 3D locations.