Based on the fast Fourier wave superposition spectrum method, a new equivalent source method (ESM) with a sparse sampling technique is proposed. First, the equivalent source intensities are expanded on a rectangular virtual surface using a bidirectional Fourier series, resulting in a semi-analytic and half-numerical acoustic pressure expression. The Fourier coefficients result in good sparsity for continuous acoustic pressures from structural vibration sources, and the proposed sparse sampling method can further reduce correlation in the measurement matrix. Better results can be obtained by solving the l1 norm optimization problem. Finally, the method was verified using several examples. The proposed method offers two main advantages compared with the traditional compressive equivalent source method: (1) the unknown source intensity vector is expanded into a bidirectional Fourier series, thereby transforming an unknown source intensity vector into a sparse Fourier coefficient vector, which has better sparsity; (2) the proposed method constructs a random sampling matrix, which is expanded into a sparse sampling matrix by random distribution, thereby improving the reconstruction accuracy of planar near-field acoustic field compared with the traditional random position sampling method reducing correlation in the transfer matrix.