In existing imaging algorithms for wide-beam multiple-receiver synthetic aperture sonar (SAS) systems, the double-square-root (DSR) range history of each receiver is generally converted into the sum of a single-square-root (SSR) range history and an error term using displaced phase center aperture (DPCA) approximation. Therefore, before imaging, each receiver’s error term needs to be individually compensated in the azimuth frequency domain, which is computationally expensive. As a result, a novel wide-beam multiple-receiver SAS system algorithm with low complexity and high precision is suggested. First, the translation relationship between the range histories of the reference receiver and other receivers is used to derive an SSR approximation range history that takes into account the azimuth variance of the non-stop-hop-stop time while ignoring differential range curvature (DRC) between the range histories from different receivers. Then, using the principle of stationary phase (POSP), the two-dimensional (2-D) spectrum of the point target is obtained. Finally, the multiple-receiver data are transformed into monostatic SAS-equivalent data for imaging after phase correction, time delay correction, and azimuth reconstruction. The range-Doppler (RD) algorithm is used as an example to explain the specific steps of the proposed approach. Simulation data and ChinSAS data experiments verify that the proposed algorithm achieves an imaging performance that is comparable to that of the existing wide-beam algorithm, but with much higher computational efficiency, making it suitable for real-time imaging.

When a synthetic aperture sonar (SAS) system operates under low-frequency broadband conditions, the azimuth range coupling of the point target reference spectrum (PTRS) is severe, and the high-resolution imaging range is limited. To solve the above issue, we first convert multi-receivers’ signal into the equivalent monostatic signal and then divide the equivalent monostatic signal into range subblocks and the range frequency subbands within each range subblock in order. The azimuth range coupling terms are converted into linear terms based on piece-wise linear approximation (PLA), and the phase error of the PTRS within each subband is less than π/4. Then, we use the chirp-z transform (CZT) to correct range cell migration (RCM) to obtain low-resolution results for different subbands. After RCM correction, the subbands’ signals are coherently summed in the range frequency domain to obtain a high-resolution image. Finally, different subblocks are concatenated in the range time domain to obtain the final result of the whole swath. The processing of different subblocks and different subbands can be implemented in parallel. Computer simulation experiments and field data have verified the superiority of the proposed method over existing methods.