The invention belongs to the field of
sonar signal processing, and discloses a windowing two-dimensional unrolling multi-beam power spectrum
estimation algorithm, which comprises the following steps of: pre-calculating and storing a guide vector universal set, namely a
point spread function, of each scanning direction (theta, F) according to an array structure and a
processing frequency band; performing windowing
Fourier transform by adopting the k-lattice element
time domain data to obtain the nth batch of
array element frequency domain data; adopting a pre-stored guide vector set to carry out windowing beam former
processing to obtain a multi-beam power spectrum
estimation result; carrying out two-dimensional unrolling processing on the multi-beam power spectrum windowed in the space domain-
time domain by using an R-L iterative
algorithm in a discrete two-dimensional integral form, and estimating the real
energy distribution estimation PBout of the target; and after the unrolling processing is completed, carrying out two-dimensional spectrum peak search on the PBout, carrying out background estimation by taking the minimum value in the extreme points as background estimation, and carrying out background
smoothing on all non-extreme points to obtain new multi-beam power spectrum estimation. According to the method, the performance of passive
sonar narrowband warning is effectively improved, and the method has high
engineering practical significance.