Spatial Spectral Enhancement of Broadband Signals in a Towed Array Using Deconvolved Subband Peak Energy Detection

Conventional energy detection is a robust method that is usually applied to underwater broadband acoustic signal processing for towed arrays. Due to its low resolution, the weak target detection performance of conventional energy detection is severely degraded in shallow sea environments with strong acoustical reverberation. Subband peak energy detection is an effective method to improve the display resolution of conventional energy detection. However, subband peak energy detection produces false alarms due to the presence of high sidelobe levels. In order to improve the underwater target detection performance, a deconvolved subband peak energy detection method for towed arrays is proposed in this paper. Compared with conventional beamforming, minimum-variance distortionless response with forward–backward averaging and diagonal loading algorithm and subband peak energy detection, the proposed method could robustly provide higher-resolution results and suppress the fake peaks induced by subband peak energy detection. The performance of the proposed method was evaluated with simulation results, and the sea experimental data processing results show that the proposed method is effective in engineering applications.


Introduction
Towed array sonars with a large number of acoustical sensors and large apertures are widely used for underwater weak signal detection [1], achieving high spatial resolution and signal amplification. Geometrically, a towed array sonar is a long flexible hose attached to the ship or the underwater vehicle with a connecting cable [2], as shown in Figure 1. During the motion of the tow vehicle, the sensors are thus subjected to hydrodynamic forces and under constant motion, and disorder of the array shape appears.
Conventional beamforming (CBF), also known as Bartlett [3], is the earliest spatial spectrum estimation technique. Because of its insensitivity to minor array shape uncertainty and ease of implementation, conventional beamforming has been frequently utilized in towed array sonar systems. Due to the inadequate resolution and high sidelobes level of conventional beamforming, the identification of low signal-to-noise ratio (SNR) signals in the presence of a noisy interferer using a towed array sonar is constrained [4]. Many adaptive beamforming methods have been developed to overcome the disadvantages of conventional beamforming such as the minimum-variance distortionless response (MVDR) and multiple signal classification (MUSIC) [3]. However, distortion of the array shape can cause significant degradation in the high-resolution beamforming process, resulting in Conventional beamforming (CBF), also known as Bartlett [3], is the earliest spatial spectrum estimation technique. Because of its insensitivity to minor array shape uncertainty and ease of implementation, conventional beamforming has been frequently utilized in towed array sonar systems. Due to the inadequate resolution and high sidelobes level of conventional beamforming, the identification of low signal-to-noise ratio (SNR) signals in the presence of a noisy interferer using a towed array sonar is constrained [4]. Many adaptive beamforming methods have been developed to overcome the disadvantages of conventional beamforming such as the minimum-variance distortionless response (MVDR) and multiple signal classification (MUSIC) [3]. However, distortion of the To lessen the impact of array disturbance, a tremendous number of array shape uncertainty estimation methods have been proposed, which can be classified into self-calibrating methods and calibrating source-dependent methods. For the self-calibrating method, the most common way is to place heading and depth measurement sensors at different points along the towed array, which provides localized horizontal and vertical information on the Remote Sens. 2022, 14, 3008 3 of 17 transverse displacements of the towed array sonar [17,18]. Alternatively, the positions can be determined by developing a water dynamical motion model of towed array [5,19,20]. Calibrating source-dependent method estimate array shapes can be conducted by processing known sound sources [21]. The disadvantages of these array shape estimation methods are that implementation is complicated, and the computational complexity is high, making real-time application difficult [22].
Another approach to enhancing the performance of spatial spectrum estimation is to exploit the broadband information of the signals. For a passive towed array sonar at low frequencies below 1000 Hz, broadband processing can be used for weak acoustic signal detection [1]. Methods such as conventional energy detection can be used to increase the array gain and narrow the beamwidth after the CBF process. Conventional energy detection can be defined as when the signal has a nonzero bandwidth, and the CBF output from the towed array can be incoherently averaged through the whole bandwidth to improve the performance of the beamforming results [23]. In shallow sea environments with strong acoustical reverberation, the acoustic energy of the target is weak, and the performance of conventional energy detection sharply decreases [24]. Subband peak energy detection (SPED) can be expected as a solution in such situations, and it can improve the display resolution of conventional energy detection [25]. SPED provides the possibility for detecting contacts of which the broadband signature will not be at minimum detectable levels but contain significant power in certain narrow spectral ranges. Although SPED can improve the direction resolution and provide better detection ability in real acoustic environments, SPED is still affected by its high sidelobe levels, and it generates false alarms [26]. To improve the performance of SPED, one can resort to upgrading the algorithms of energy detection or reducing the sidelobe levels of the beamformed output [27,28].
In this paper, a deconvolved subband peak energy detection (dCv-SPED) method is proposed to improve the performance of spatial-spectral estimation. A SPED algorithm is commonly used as a postprocessing algorithm for beamforming; however, the highresolution deconvolution beamforming and SPED were treated as a joint problem in our work. Compared with SPED, the dCv-SPED algorithm makes full use of the prior information of the beam pattern by using deconvolution; the peak points obtained at different subbands are thereby more accurate and spikier. The advantages of the proposed method are summarized as follows: (1) It achieves high resolution with a low sidelobe level and overcomes the suppression of weak signals near strong interferers in deconvolution beamforming; (2) Such an algorithm can suppress the fake peak induced by SPED while retaining a low SNR signal detection capability; (3) It is insensitive to array perturbation and the snapshot deficiency problem. This paper is organized as follows: the formulation of the problem is given in Section 2; Section 3 utilizes a simulation to study the performance of the proposed method, and the bearing time record results for 200-400 Hz broadband signals based on a towed array sonar is also introduced; Section 4 provides a conclusion of our research.

Methods
Consider a uniform linear array that has N equally spaced hydrophones, with d being their spacing. Sources at the far field of the array radiate narrowband acoustic signals with an angular direction θ i The array receives these signals as the incident plane waves. The received signal in the frequency domain can be expressed as: where F = F 1 F 2 · · · F D T is the frequency-domain snapshot vector of D source signals; N n denotes the isotropic white noise snapshot. vector can also be defined in u-space as v I = 1 e −jkdu i · · · e −jk(N−1)du i T , where u i = sin θ i .

Conventional Beamforming
In order to estimate the direction of arrival, the spatial-spectral obtained by the CBF in u-space is given by: where a H (u) = 1 N 1 e −jkdu · · · e −jk(N−1)du T is the steering vector; E(·) denotes the average over snapshots; and H denotes the conjugate transpose. We assumed that no two signals were coherent, and the noise component at each sensor had a spectral height σ 2 n .
Then, using , we can rewrite Equation (2) as: where the sinc function is given by sin(πx)/πx; f is the narrowband frequency; and c is the sound speed. If the plane waves are broadband signals, the M primarily frequency bins that we are interested in are f m = f 1 , · · · , f M , which can be rewritten as a spatial-spectral function of f m as: The average spatial-spectral across the primarily frequency bins is obtained by: (5)

MVDR-DL-FB Beamforming
In this subsection, we present an improved version of the MVDR algorithm, which is commonly used in towed array signal processing. Despite the fact that this algorithm has no relation to the algorithm proposed in this paper, we used it for performance comparison purposes. Based on the signal model above, the MVDR beamformer design is given by: where R X = E XX H is the covariance matrix and u is the desired signal direction in the u-space. Usually, the Lagrange multiplier method is used to solve the constrained minimization problem in Equation (6). The objective function to be minimized is: where λ is the Lagrange multiplier. Determining the complex gradient of L with respect to w H and letting it be zero, the optimal weight vector is: According to the constraint conditions of Equation (6), and solving it gives: Substituting Equation (9) into Equation (8), we can rewrite the optimal weight vector as: Thus, the spatial-spectral in the u-space of the MVDR is obtained by: By adjusting the weight factor, w, the MVDR algorithm minimizes the total outputting power, while ensuring the outputting signal from the desired direction is distortionless.
In many situations, the covariance matrix, R X , is unknown and has to be estimated from multiple sample snapshots. Assuming we obtain K copies of the array frequencydomain snapshots, say X[1], · · · , X[K], we may use the sample covariance matrix

X[k]X[k]
H as an estimate of R X . The estimation accuracy ofR X has a major influence on the performance of the MVDR algorithm. While facing the snapshot deficiency problem,R X does not represent statistical characteristics of signals properly and, furthermore, causes an ill-conditioned problem for the matrix inverse. Assuming that the uniform linear array steering vector that we are dealing with is conjugate symmetric, say: where * denotes a conjugate. J is an exchange or reflection matrix represented by: Then, we can introduce the forward-backward (FB) averaging technique to improve the snapshot deficiency problem by replacingR X with the forward-backward averaged sample covariance matrix,R X,fb , as: References [12][13][14] proved that forward-backward average brings a performance improvement and has computational advantages.
Manifold mismatch problems caused by array perturbations also deteriorate the behavior of the MVDR. In some mismatched cases, the MVDR beamformer nullifies the signal in the desired direction and leads to unacceptable performance loss. To deal with this problem, a technic called diagonal loading (DL) is adopted. Diagonal loading abates the sensitivity of the MVDR beamformer to array perturbations by adding a diagonal matrix toR X . Combining forward-backward and diagonal loading, a new covariance matrix is reconstructed as:R X,fb,dl =R X,fb + σ L I (15) where σ L is the loading versus, and it can be estimated by quadratic constraint: where the range of the constraint value ξ [7] is: Initially, by using the equality constraint in Equation (16) and combining Equation (6), the objective function to be minimized is obtained as: The solution in Equation (18) can be written as: Reference [3] proves that theŵ H optŵ opt decreases with the increase in σ L . Since the analytical solution of σ L cannot be directly obtained, an iterative process can be used to solve the constraint equation: with the initial value σ The iterative process of Equation (20) is repeated and w H optŵ opt is calculated at each step until the constraint in Equation (16) is satisfied. Thus, the spatial-spectral in the u-space of the MVDR-FB-DL is obtained by: For broadband signals, the spatial-spectral can be estimated by the average beam power across primarily frequency bins such as in Equation (5).

Formulations of SPED
Reconsider beamforming while the signal of interest is broadband. In Section 2.1, we simply averaged the beam power across frequency bins, the so-called conventional energy detection. Different from conventional energy detection, SPED only sums the beam peak across frequency bins. The first step of SPED is to estimate the floor of the beam power at each subband, sayσ(u, f m ). This is conducted since the SNR is used as the energy value, defined by: (22) can be written as: Then, the spatial-spectral of SPED in u-space is obtained as: where: Equations (24) and (25) ensure that only the energy at the peak accumulates across the subband. When there is a broadband signal coming from a certain azimuth in space, there should be some subband azimuth spectra that show local peak values at that azimuth. When there is only noise in a certain azimuth, although some subband azimuth spectra may show local peaks at that azimuth, these local peaks are random, and it is impossible for most subband azimuth spectra to show local peaks at the same azimuth. This property is called spatial coherence. By performing local peak accumulation, the spectral value of the signal azimuth is strengthened, and the noise spectral value of other azimuths is weakened; thus, the azimuth resolution can be improved, and the background level can be reduced.
SPED can be regarded as a frequency selection algorithm that selects the frequency of the beam response at the peak when the beam power is accumulated and discards the frequencies of those non-peak points. However, the extracted peaks cannot make full use of the prior information on the beam pattern; therefore, the extracted peaks have a large error. At the same time, when the broadband signal has a strong line spectrum, the peaks extracted from the sidelobes of the individual frequencies are difficult to suppress, resulting in false alarms.

Implementation of the dCv-SPED
To improve the performance of SPED, a deconvolution algorithm was introduced to obtain more accurate peaks while suppressing the sidelobes. Following the work of Yang [16], Equation (23) can be rewritten in the integral form as: is the source spatial distribution in u-space, and the beam pattern corresponding to frequency f m is given by: Since B SNR (u, f m ) can be regarded as the convolution of the beam pattern with the source spatial distribution, a method called deconvolution is applied to recover the source spatial distribution from B SNR (u, f m ). In our study, the deconvolution problem in Equation (26) was solved with a fast gradient projection method called the fast iterative shrinkage-thresholding algorithm [29].
Given Bp f m as a point-spread function, the iterative procedure of the fast iterative shrinkage-thresholding algorithm is given by: where the superscript (p) is the iteration number; P + is the Euclidean projection of variables onto the nonnegative orthant R n + ; f is a convex function with a Lipschitz continuous gradient ∇ f and a Lipschitz constant By setting p = P as the stopping criterion, the estimated source spatial distribution iŝ S f m (u) (P) . Then, B SNR (u, f m ) is replaced by S f m (u) (P) in Equation (25) to obtain: and the spatial-spectral of the dCv-SPED can be expressed as: The summarized steps of the dCv-SPED algorithm are listed in Algorithm 1. We applied numerical simulations and experimental data to illustrate the performance of the dCv-SPED. Algorithm 1. Summary of the dCv-SPED beamforming algorithm.
Given the received acoustic sensors data 1. Achieve spatial-spectral B(u, f m ) by conventional beamforming; 2. Estimate the noise floorσ(u, f m ) to obtain B SNR (u, f m ) using Equation (22); 3. Estimate the source spatial distribution S f m (u) (P) by setting the stopping criterion P in Equation (28); 4. Obtain B dCV−peaks (u) by Equation (18) and generate the spatial-spectral of the dCv-SPED using Equation (30).

Numerical Simulations
Now, the role of the dCv-SPED in the broadband spatial-spectral estimation of a uniform linear array is demonstrated, and its performance in the presence of array perturbations was studied using numerical examples. A uniform linear array with N = 64 elements receiving signals with a frequency band of [200 Hz, 400 Hz] was considered. The number of snapshots processed was 32. The half-wavelength inter-element spacing of the uniform linear array was d = λ h /2 = 1.875 m, where λ h was the wavelength corresponding to the highest frequency, f h = 400 Hz.
In the first simulation, the presence of signals with a strong line spectrum was considered first to examine the false alarm suppression and weak target detection ability of the dCv-SPED algorithm. Two ship noise signals arriving at the directions of 0 • [u = sin(0 • ) = 0] and 11.54 • [u = sin(11.54 • ) ≈ 0.2] were generated. The SNRs of the two targets were 10 dB (at 0 • ) and 0 dB (at 11.54 • ). Figure 2 shows the spatial spectra estimated by CBF, SPED, MVDR-FB-FL, and dCv-SPED. In order to generate the same dynamic range, all spectra plots were normalized with the max value set at 0 dB. As shown in the figure, for the strong signal at u = 0, the beamwidth of the MVDR-FB-DL algorithm had the highest resolution, and its beamwidth was ∼ 0.28 • (measured at −3 dB point), followed by the dCv-SPED with a beamwidth ∼ 1.10 • . For the weak signal at u = 0.2, the dCv-SPED algorithm had the highest resolution, and its beamwidth was ∼ 1.24 • , followed by SPED with a beamwidth at ∼ 1.58 • . The dCv-SPED generated the highest dynamic range (target beam power in dB minus noise floor power in dB) of the weak target spatial spectra at u = 0.2, and its dynamic range was ∼ 18 dB, followed by SPED with a dynamic range of ∼ 15.5 dB. Compared with the SPED algorithm, dCv-SPED effectively suppressed the false peak caused by the SPED algorithm. 9 of 17 . u 02, the dCv-SPED algorithm had the highest resolution, and its beamwidth was .°1 24 , followed by SPED with a beamwidth at .°1 58 . The dCv-SPED generated the highest dynamic range (target beam power in dB minus noise floor power in dB) of the weak target spatial spectra at . u 02, and its dynamic range was 18 dB , followed by SPED with a dynamic range of . 15 5dB . Compared with the SPED algorithm, dCv-SPED effectively suppressed the false peak caused by the SPED algorithm.  For each loop, the two most energetic peaks in the spatial spectrum were chosen as the estimated direction in the u-space,û 1 andû 2 (û 1 <û 2 ), and the RMSE was: Eventually, the RMSE was averaged across all 200 Monte Carlo trials. Figure 3 presents the RMSEs versus the input noise powers for the four algorithms. The dCv-SPED reached the lowest RMSE over the other algorithms for a noise power lower than 25 dB. The RMSE of all four algorithms decreased as the noise power decreased. The proposed algorithm achieved the minimum RMSE for a noise power below 10 dB, while the SPED algorithm achieved the minimum RMSE for a noise power below 5 dB. The CBF and MVDR-FB-DL had relatively high RMSEs. Figure 4 displays the probability of resolution versus noise power. We can say that a strong signal and a weak signal were successfully resolved if: The probability resolution was computed by dividing the number of successful resolutions by the Monte Carlo trials. The resolution probability of all four algorithms increased as the noise power decreased. For noise power above 15 dB, the proposed algorithm obtained the highest resolution probability, and its resolution probability reached 1 at a noise power of 20 dB, while the SPED's resolution probability reached 1 at a noise power of 15 dB. Figure 5 provides the spatial spectra estimated by the four algorithms when the noise was 10 dB. For the proposed algorithm, the two targets were clearly resolved, and the peak position at the weak target was accurate. The SPED algorithm could also resolve two targets, but the peaks' position near weak targets was somewhat inaccurate, which waŝ u 2 = 0.095.     For the third simulation, the actual and hypothetical array element positions are displayed in Figure 6. The position errors were determined by random perturbations in the x-and y-directions. We assumed the perturbation was statistically independent zeromean Gaussian random variables with standard deviation p σ , where p h σ = λ 0.1 It was assumed that the actual position of the array was unknown and beamformed according to a uniform linear array. The parameter setting of the receiving signals and noise was the same as in the first simulation.  Figure 7 shows the spatial spectra estimated by the four algorithms under array perturbations. As the figure shows, the proposed algorithm provided the highest resolution for the spatial spectrum estimation of both strong and weak targets, and its beam width was .°1 53 and .°1 81 , respectively. When the array was perturbed, the beam widths Next, the presence of array perturbations was considered. In towed array beamforming, the array shape is usually assumed as a uniform linear array, and the position of the array elements is known. However, in most situations, the array manifold is subject to random perturbations due to the flexible structure of the towed array. The ability to suppress and discriminate the interference of the beamforming algorithms decreased sharply. Therefore, it is important to study the performance of the algorithm under array disturbance.
For the third simulation, the actual and hypothetical array element positions are displayed in Figure 6. The position errors were determined by random perturbations in the xand y-directions. We assumed the perturbation was statistically independent zero-mean Gaussian random variables with standard deviation σ p , where σ p = 0.1λ h It was assumed that the actual position of the array was unknown and beamformed according to a uniform linear array. The parameter setting of the receiving signals and noise was the same as in the first simulation. For the third simulation, the actual and hypothetical array element positions are displayed in Figure 6. The position errors were determined by random perturbations in the x-and y-directions. We assumed the perturbation was statistically independent zeromean Gaussian random variables with standard deviation p σ , where ph σ = λ 0.1 It was assumed that the actual position of the array was unknown and beamformed according to a uniform linear array. The parameter setting of the receiving signals and noise was the same as in the first simulation.  Figure 7 shows the spatial spectra estimated by the four algorithms under array perturbations. As the figure shows, the proposed algorithm provided the highest resolution for the spatial spectrum estimation of both strong and weak targets, and its beam width was .°1 53 and .°1 81 , respectively. When the array was perturbed, the beam widths  Figure 7 shows the spatial spectra estimated by the four algorithms under array perturbations. As the figure shows, the proposed algorithm provided the highest resolution for the spatial spectrum estimation of both strong and weak targets, and its beam width was ∼ 1.53 • and ∼ 1.81 • , respectively. When the array was perturbed, the beam widths of SPED for the strong and weak targets were ∼ 1.56 • and ∼ 1.87 • , close to those of the proposed algorithm, but SPED produced strong false alarms at approximately u = 0.2 and u = 0.6. Due the array perturbations, the dynamic range and spectra resolution of the MVDR-FB-DL declined sharply. The explanation for this is that the MVDR was affected by the array perturbations, although the improved algorithms of the forward-backward averaging and diagonal loading were adopted to obtain spatial spectra, the cost was a decrease in the dynamic range, which led to the reduction in the weak target detection effect. of SPED for the strong and weak targets were .°1 56 and .°1 87 , close to those of the proposed algorithm, but SPED produced strong false alarms at approximately . u 02 and .
u 06. Due the array perturbations, the dynamic range and spectra resolution of the MVDR-FB-DL declined sharply. The explanation for this is that the MVDR was affected by the array perturbations, although the improved algorithms of the forward-backward averaging and diagonal loading were adopted to obtain spatial spectra, the cost was a decrease in the dynamic range, which led to the reduction in the weak target detection effect. In the final simulation, we evaluated the performance of the four algorithms in terms of RMSE and resolution probability when the array was perturbed. The signal and noise parameters were identical to those used in the second simulation, and the array perturbation parameters were identical to those used in the third simulation. Figure 8 presents the root mean square errors of the estimated direction in the u-space vs. noise power for the four algorithms under array perturbations. As shown in the figure, the effect of the array perturbation increased the RMSEs of the four algorithms. Below a 25 dB noise power, the dCv-SPED outperformed the other algorithms and reached the minimum RMSE for a noise power below 15 dB. Figure 9 provides the probability of resolution versus noise power under array perturbations. For noise power above 10 dB, the proposed algorithm reached the highest resolution probability, and its resolution probability reached 1 at a noise power of 20 dB, while SPED's resolution probability reached 1 at a noise power of 10 dB. As Figure 10 shows, SPED's performance was insufficient in that it may treat false peaks around strong targets as real targets; meanwhile, the poor performance of the CBF and MVDR-FB-DL was attributed to insufficient weak target detection performance under strong interference. In the final simulation, we evaluated the performance of the four algorithms in terms of RMSE and resolution probability when the array was perturbed. The signal and noise parameters were identical to those used in the second simulation, and the array perturbation parameters were identical to those used in the third simulation. Figure 8 presents the root mean square errors of the estimated direction in the u-space vs. noise power for the four algorithms under array perturbations. As shown in the figure, the effect of the array perturbation increased the RMSEs of the four algorithms. Below a 25 dB noise power, the dCv-SPED outperformed the other algorithms and reached the minimum RMSE for a noise power below 15 dB. Figure 9 provides the probability of resolution versus noise power under array perturbations. For noise power above 10 dB, the proposed algorithm reached the highest resolution probability, and its resolution probability reached 1 at a noise power of 20 dB, while SPED's resolution probability reached 1 at a noise power of 10 dB. As Figure 10 shows, SPED's performance was insufficient in that it may treat false peaks around strong targets as real targets; meanwhile, the poor performance of the CBF and MVDR-FB-DL was attributed to insufficient weak target detection performance under strong interference.

Experiments
In this section, the dCv-SPED algorithm was evaluated with real data collected during a sea trial on 22 October 2017, in the east sea of China. The array had 32 elements uniformly spaced at 2 m (design frequency: 375 Hz). The acoustic source was located at 35 m below the sea's surface, and the ship with the towed array was 2~3 km away from the acoustic source. The sampling frequency was 8 kHz. The time window length was 600 s, during which the acoustic source signal starting from t 0 at . u 06 and several noise signals radiated by many ships were recorded. Twenty snapshots of the data were used for beamforming. The processing frequency band of the beamforming was 200-400 Hz. Figure 11 shows the bearing time records obtained by the CBF, SPED, MVDR-FB-FL, and dCv-SPED. It can be seen from the figure that the dCv-SPED improved the performance of the increased spatial resolution substantially; therefore, weak contact trajectories next to many strong targets could be detected. Compared to SPED, the dCv-SPED generated a clearer bearing time record with fewer false artifacts. In order to more clearly reflect the performance of the algorithms, we give the output of the beam power of the first second in Figure 12. The trajectories of the bearing time records indicate that there were three targets intersecting at the location marked by the ellipse in Figure 12. As can be seen from the figure, the SPED and dCv-SPED could effectively distinguish three close targets, while the other two algorithms could not. However, it can be seen from the labels in Figure 12c that SPED produced spurious peaks with respect to the sound source signal, while the proposed algorithm did not.

Experiments
In this section, the dCv-SPED algorithm was evaluated with real data collected during a sea trial on 22 October 2017, in the east sea of China. The array had 32 elements uniformly spaced at 2 m (design frequency: 375 Hz). The acoustic source was located at 35 m below the sea's surface, and the ship with the towed array was 2~3 km away from the acoustic source. The sampling frequency was 8 kHz. The time window length was 600 s, during which the acoustic source signal starting from t = 0 at u ≈ 0.6 and several noise signals radiated by many ships were recorded. Twenty snapshots of the data were used for beamforming. The processing frequency band of the beamforming was 200-400 Hz. Figure 11 shows the bearing time records obtained by the CBF, SPED, MVDR-FB-FL, and dCv-SPED. It can be seen from the figure that the dCv-SPED improved the performance of the increased spatial resolution substantially; therefore, weak contact trajectories next to many strong targets could be detected. Compared to SPED, the dCv-SPED generated a clearer bearing time record with fewer false artifacts. In order to more clearly reflect the performance of the algorithms, we give the output of the beam power of the first second in Figure 12. The trajectories of the bearing time records indicate that there were three targets intersecting at the location marked by the ellipse in Figure 12. As can be seen from the figure, the SPED and dCv-SPED could effectively distinguish three close targets, while the other two algorithms could not. However, it can be seen from the labels in Figure 12c that SPED produced spurious peaks with respect to the sound source signal, while the proposed algorithm did not.

Conclusions
In this paper, a robust algorithm, dCv-SPED, was proposed to obtain high-resolution and interference-free spatial spectra. Through the joint processing of deconvolution and frequency selection problems, the algorithm made full use of the prior information of the beam. The simulation results showed that the algorithm could accurately and stably estimate weak signals near strong interference. When the array was disturbed, the performance of the algorithm was slightly reduced, but it still had a higher resolution and antiinterference advantages than the traditional algorithm. Meanwhile, the proposed algorithm could effectively improve the dynamic range of the beam power for the weak target, which brought a better display effect for the bearing time record, which is beneficial to the discovery of the target of interest. Finally, experimental data were used to test the proposed algorithm. The results revealed strong performance characteristics that matched those obtained in the simulated data.

Conclusions
In this paper, a robust algorithm, dCv-SPED, was proposed to obtain high-resolution and interference-free spatial spectra. Through the joint processing of deconvolution and frequency selection problems, the algorithm made full use of the prior information of the beam. The simulation results showed that the algorithm could accurately and stably estimate weak signals near strong interference. When the array was disturbed, the performance of the algorithm was slightly reduced, but it still had a higher resolution and anti-interference advantages than the traditional algorithm. Meanwhile, the proposed algorithm could effectively improve the dynamic range of the beam power for the weak target, which brought a better display effect for the bearing time record, which is beneficial to the discovery of the target of interest. Finally, experimental data were used to test the proposed algorithm. The results revealed strong performance characteristics that matched those obtained in the simulated data.