Virtual Reconstruction-Based Robust Adaptive Beamforming for Distributed Digital Subarray Antennas

For the distributed digital subarray antennas (DDSA), the conventional beamforming may give rise to grating lobes, high sidelobes, and other problems. In this paper, the gaps between the subarrays are filled with virtual array elements, and then the DDSA can form a virtual contiguous array. More concretely, based on the direction-of-arrival (DOA) estimation of the signal sources, the interference components of the virtual elements and the interference-plus-noise covariance matrix (INCM) of the virtual contiguous array can be reconstructed. At low signal-to-noise ratio (SNR), the DOA estimation of the desired signal is implemented by subarray adaptive beamforming. Finally, with the steering vector of the desired signal and reconstructed INCM, the weight vector of the proposed beamformer can be obtained, which must be applied to the rearranged data matrix received by the actual and virtual elements. Simulation results are provided to demonstrate the effectiveness of the proposed algorithm.


Introduction
Compared with a single antenna, the distributed digital subarray antennas (DDSA) can provide better spatial resolution, wider beam coverage, and macro diversity gain [1].
At present, there are some research results about the DDSA, which can be generally divided into two categories, i.e., direction-of-arrival (DOA) estimation and digital beamforming. In terms of DOA estimation, the long baseline separation technique in [11] was proposed to achieve superresolution direction finding. In [12,13], the estimation of signal parameters via rotational invariance techniques-(ESPRIT-) based methods was applied to sparse arrays. In [14], by combining the multiple signal classification (MU-SIC) algorithm with sparse subarrays, the closed-form eigenstructure-based DOA estimation was proposed to eliminate ambiguity. In [15], by combining virtual filling of the gaps between the subarrays with the modified matrix pencil method, the grating lobes were suppressed, and then the DOA of the desired signal can be accurately determined. In [16], the virtual elements between inside and outside the subarrays were constructed with different approaches, and then the superresolution DOA estimation was realized with sparse subarrays. In the aspect of beamforming, to reduce the parameters storage and computational load, the distributed and parallel subarray beamforming was proposed in [7]. e main characteristic of this method is to reduce the hardware resources, but the anti-interference performance is not significantly improved. In [4], a Fresnel-based frequency domain adaptive beamforming was proposed to eliminate the envelope migration problem. e research object of this method is the large aperture distributed array radar. However, due to the negligible envelope migration in the small DDSA system, this beamforming algorithm is not suitable. In [17], the distributed linearly constrained minimum variance (LCMV) beamformer was proposed to suppress interference and optimize the signal-to-noise ratio (SNR), which greatly reduced the transmission requirements. But this method requires the accurate DOA of the desired signal in advance, especially at low SNR. In [6], a spatial multi-interference suppression method based on joint adaptive weight was proposed for the distributed array radar to improve the performance of the interference suppression in the case of massive sidelobe interferences. In addition to the DOA of the desired signal, this beamforming algorithm needs interference-plus-noise covariance matrix (INCM) reconstruction for each subarray, which greatly increases the computational complexity. In [5], a two-stage adaptive beamforming technique for distributed array radars was proposed to cancel multiple mainlobe and sidelobe jammings. However, this algorithm requires the full distributed array aperture to be sufficiently large. If not, the desired signal may be suppressed while cancelling the mainlobe jammings.
In this paper, the virtual filling technique is applied to the DDSA, and the gaps between the subarrays are filled with the virtual elements. By dividing the DDSA into two rotationinvariant subarray sets, the DOAs of the interfering signals can be accurately estimated, and the received data on the virtual elements can be reconstructed. en, a virtual contiguous array can be formed, which significantly increases the degrees of freedom (DOF) of the proposed beamformer, and then a novel beamforming algorithm is devised based on the INCM reconstruction. In addition, to obtain the DOA of the desired signal at low SNR, we propose a DOA estimation method by estimating the subarray steering vector. After that, simulation results demonstrate the effectiveness of the proposed algorithm. Our main contributions are summarized as follows: (1) We apply the ESPRIT-like method to DOA estimation of the incident signals by dividing the DDSA into two rotation-invariant subarray sets and estimate the DOA of the desired signal at low SNRs by solving the optimization problem. (2) Based on the DOA estimation, we estimate the received interfering signals on the virtual elements and apply the data to the INCM reconstruction of the virtual contiguous array. (3) We give the complete data treatment of the proposed beamforming involved in the weight vector and reconstructed sampling data set. e remainder of this paper is organized as follows. In Section 2, the DDSA signal model and background are described. A robust adaptive beamforming based on virtual array elements and INCM reconstruction is proposed in Section 3. Simulation examples are provided in Section 4. Finally, some conclusions are drawn in Section 5.

DDSA Signal Model and Background
Consider a linear DDSA system of M identical subarrays, which belongs to the uniform linear array (ULA) of N omnidirectional sensors, with a half wavelength spacing d. As shown in Figure 1, the solid circles are actual elements, and the nonfilled circles are virtual sensors. e notch widths between the subarrays are D i � h i d (i � 1, 2, . . . , M − 1), where D i is an integer multiple of the interelement spacing d. If the first sensor of subarray 1 is considered as the reference one of the whole antenna, the location of the n-th element in subarray m can be represented by p(n, m)d, where p(n, m) is written as (1) In Figure 1, there are M − 1 virtual subarrays, and the location of the l-th element in virtual subarray m can be expressed as q(l, m)d, where q(l, m) is represented by (2) Accordingly, the element positions of subarray 1, 2, . . ., erefore, all of the sensors can form a virtual contiguous ULA with Z � MN + m−1 i�1 h i elements. Assume that the array receives S + 1 narrowband signals from far-field signal sources. e desired signal impinges upon the DDSA from direction θ 0 , and S interfering signals come from the directions θ 1 , θ 2 , . . . , θ S , respectively. e DDSA observation vector x(k) ∈ C MN×1 at time k can be modeled as where x d (k), x i (k), and n(k) are the desired signal, interference, and noise, respectively, which are statistically independent of each other. e element of the noise vector n(k) is assumed to be independent and identically distributed complex Gaussian variables with zero mean and variance δ 2 n . e beamforming output at time k can be expressed as where w � [w 0 , w 1 , . . . , w MN− 1 ] T denotes the weight vector of the beamformer and (·) H denotes the Hermitian transpose operation. e output SINR can be obtained from where δ 2 d denotes the power of the desired signal, λ denotes the carrier wavelength, and R i+n ∈ C MN×MN denotes the theoretical INCM, which can be reconstructed by where δ 2 i denotes the power of the i-th interfering signal and I MN stands for an MN × MN identity matrix. By maximizing the output SINR (5), the optimal weight vector can be expressed as which is a function of the interference energy δ 2 i , noise power δ 2 n , and the DOAs θ s (s � 0, 1, . . . , S) of the desired signal and interferers. From (7), the MN × 1 weight vector w opt can be obtained, which corresponds to MN physical elements of the DDSA. However, the sparse distribution of M subarrays may result in the grating lobes and high sidelobe levels. erefore, to suppress grating lobes and high sidelobes, it is a practical way to reconstruct the received signals of H � M−1 i�1 h i virtual elements and form a virtual contiguous ULA. By reconstructing the INCM of the virtual ULA, the adaptive beamforming is used for suppressing the interfering signals. To make the desired signal keep phase congruency among multiple channels, only the interfering signals received by the virtual elements are reconstructed. In addition, it is difficult to accurately estimate the DOA of the desired signal at low SNR, which may result in the drift of beam direction.

Proposed Algorithm
In this section, a robust adaptive beamforming algorithm using virtual array elements is proposed. Based on the DOA estimation of the desired signal and interferers, the received data on the virtual elements and the INCM of the virtual contiguous array can be reconstructed. Introducing the virtual array elements can greatly increase the DOF of the proposed beamformer and improve the anti-interference performance. With the INCM of the virtual ULA and the steering vector of the desired signal, the weight vector can be obtained, which must be used for processing the rearranged data set received by the actual and virtual elements.
3.1. DOA Estimation of the DDSA. As shown in Figure 2, the DDSA is divided into two subarray sets A z1 and A z2 . e subarray set A z1 consists of the first N − 1 sensors of each actual subarray, and the subarray set A z2 is composed of the last N − 1 elements of each actual subarray. Between A z1 and A z2 , there is a space-translation invariant d � (λ/2). erefore, the ESPRIT algorithm can be used for DOA estimation, and there is no ambiguity in the DOA estimation results.
If the DDSA observation vector at time k can be expressed as x(k) ∈ C MN×1 , the received data of the subarray sets A z1 and A z2 can be, respectively, written as , I denotes the identity matrix, and 0 (N−1)×1 denotes the (N − 1) × 1 vector of all zeros. Assume that the interfering signals are much stronger than noise, and we can neglect the effect of noise in (8). For the s-th interfering signal, the rotation-invariant relationship can be expressed as where θ s denotes the DOA of the s-th interferer. Consider that there are S incident signals received by the DDSA, and the DOAs are θ 1 , . . . , θ S , respectively. en, the rotationinvariant relationship can be further written as where denotes the array manifold of the DDSA. Because both the subspace spanned by the steering vectors and the signal subspace belong to the same subspace, there is an invertible matrix Φ, satisfying where U s denotes the signal subspace. In practice, U s can be estimated by performing eigendecomposition of the sample covariance matrix, and the number of the incident signals S can be correctly estimated by the methods in [18,19]. From (11), equation (10) can be updated to where E z1 � J z1 U s , E z2 � J z2 U s , and Ψ � Φ −1 ΓΦ. erefore, we can obtain the rotation-invariant relation matrix: International Journal of Antennas and Propagation By performing eigendecomposition of matrix Ψ, the DOA estimation can be obtained by where ρ s denotes the s-th eigenvalue of the rotation-invariant relation matrix Ψ. With the ESPRIT-like method [20,21], the DOAs of the interfering signals can be accurately estimated, and it is the same to the desired signal at high SNR.

Data Reconstruction of the Virtual Contiguous Array.
In general, the desired signal energy has a very large range of variation, and the input SNR can vary from −30 dB to +50 dB in the simulations. Due to the impact of DOA estimation errors, it is impossible to accurately estimate the phases and amplitudes of the desired signal, which can bring about bad impacts on the later data processing. However, the reconstruction errors of the interfering signals have different effects, which only change the positions and depths of the nulls. Moreover, the interference energy is usually much stronger than noise, which is helpful to improve the accuracy of DOA estimation. rough the virtual expansion, the DOF of the proposed beamformer is dramatically increased, and the nulls in the beampattern can be significantly deepened. erefore, the interference reconstruction on the virtual elements can enhance the anti-interference capability of the proposed beamformer.
Owing to the narrowband assumption, there are only phase differences among the received signals on different array elements. For the m-th virtual subarray, the received interfering signals at the elements can be expressed as where and v i exp(jϕ i ) denotes the i-th interfering signal. Casting the above equations in matrix form can get where To reconstruct the received interfering signals on the virtual subarrays, the DOAs, magnitudes, and phases are needed. Assume that all interference energy is much stronger than noise, and then the DOAs of the interfering signals can be accurately estimated by the ESPRIT-like method in the last subsection. e angular sector Θ in which the desired signal is located can be obtained in advance, and the DOA estimation result can be written as θ � [θ 0 , θ 1 , . . . , θ S ] T . On the actual subarray m, the energy of the interfering signals is much stronger than that of noise, and then the received data can be represented as where erefore, the signal vector r can be calculated by [16] r where r is the estimation of the signal vector r in (18) In the data processing, the mean values of the estimation r using multiple actual subarrays independently are adopted to cancel the effect of noise and improve the reconstruction accuracy. From (16), (18), and (20), the received data on the m-th virtual subarray can be constructed by By continuously processing multiple snapshots, the received data matrix X m and interference data matrix X i,m on the m-th virtual subarray can be obtained. erefore, two new data matrices of the virtual contiguous array can be, respectively, constructed as where X is used for the power estimation of signal sources and X p is applied to the beamforming output with the weight vector.

DOA Estimation of the Desired Signal at Low SNR.
If the input SNR of the desired signal is much greater than 0 dB, the DOA can be simultaneously estimated while processing the interfering signals, but the performance may decrease at low SNR. To solve this problem, we propose a novel DOA estimation method based on subarray adaptive beamforming. Using the m-th actual subarray, the first element is considered as the reference one. With the Capon spatial spectrum estimator, the INCM of the subarray m can be reconstructed by [22] where T is the steering vector, and Θ is the complement sector of Θ, which is the angular sector in which the desired signal is located. In order to improve the beam pointing accuracy, the signal-plus-noise covariance matrix of the subarray m can be reconstructed by [23] where δ 2 n is the noise power estimate, which can be replaced by the minimum eigenvalue of the sample covariance matrix R m . To estimate the actual steering vector a, we can establish the optimization problem as follows [23]: where θ 0 ′ is the presumed DOA of the desired signal. e optimization problem is a quadratically constrained quadratic programming problem, which can be solved by the convex optimization software, e.g., CVX [24]. Taking out the first and second items, we can get Hence, the DOA of the desired signal can be estimated by where real(·) and imag(·), respectively, denote the real part and imaginary part of a complex number.

Adaptive Beamforming Using Virtual Array Elements.
From (23), if the received data of the virtual contiguous array have been reconstructed, the covariance matrix can be calculated by To estimate the signal energy, we can formulate the convex optimization problem as follows [25,26]: where ‖ · ‖ F denotes the Frobenius norm, p(θ) � [p(θ 0 ), p(θ 1 ), . . . , p(θ S )] T denotes the signal power vector with corresponding to the power estimate at the angle

Subarray set Az1
Subarray set Az2 exp( − jμ( Z − 1 ) d sin( θ s ) )] T ∈ C Z×1 . In practice, the covariance matrix R is a positive semidefinite matrix; the inequality constraint in (31) can be neglected. Without the inequality constraint, the solution can be written as [25,26] , and vec(·) denotes the vectorization operator. Consequently, the INCM of the virtual contiguous array can be reconstructed by With the steering vector a z (θ 0 )�[1, exp(−jμdsinθ 0 ),..., exp(−jμ(Z−1) dsin θ 0 )] T corresponding to the estimated DOA θ 0 of the desired signal, the weight vector of the proposed beamformer can be written as erefore, the output data vector of the proposed beamformer can be obtained by For clarity, the flowchart of the proposed algorithm is shown in Figure 3, and the arithmetic steps can be summarized as follows: (1) Estimate the DOAs θ � [θ 0 , θ 1 , . . . , θ S ] of the desired signal and interferers (2) Calculate the approximate matrices C m (p m ) and C m (q m ) via (18) with the DOA estimation θ (3) Substitute C m (p m ) into (20), and estimate the signal vector r (4) Substitute C m (q m ) and r into (22), and reconstruct the received data X m of the m-th virtual subarray (5) Reconstruct the virtual data matrices X via (23) and X p via (24) (6) Estimate the signal energy vector p(θ) using (32) e RAB-WCPO beamformer has complexity of O(MN) 3.5 due to solving the optimization problem. e computational complexity of SMID-JAW is O(M 3 N 2 L 1 ), which is mainly dominated by the M times INCM reconstruction. erefore, the proposed beamformer has less computational complexity when the number of virtual elements is less than that of actual ones.
3.5. Discussion. So far, the DDSA has played an increasingly important role in radar, sonar, and communication systems. e research object of DDSA in this paper is the subarray antennas distributed in different parts of the same carrier, such as aircraft, vehicle, and ship. e DDSA has some desirable attributes: (1) DDSA can be applied to mainlobe interference suppression; (2) by forming a long baseline array, DDSA can achieve very accurate angular locations of targets; (3) with the distributed arrangement of subarrays, the expanded antenna aperture can provide larger detection range; (4) for a given desired performance, the signal transmitting power can be reduced significantly [6,27,28]. However, the distributed subarrays may bring about some negative impacts on the beampattern. If the subarrays are uniformly sparse, this can result in grating lobes. If the subarrays are randomly arranged, this may bring about high sidelobes. From the view of application, the bad effects may lead to ambiguity in the actual positions of targets and to excess clutter background coming from grating lobes and high sidelobes [28]. To solve these problems, the gaps between subarrays are filled with virtual elements, and then a virtual contiguous array can be formed, which can effectively suppress grating lobes and high sidelobes. Moreover, the proposed method can also be used to deal with failure of some elements in uniform arrays. By reconstructing the received signals on the faulty elements, a virtual contiguous array can be rebuilt, which is helpful to improve the spatial resolution and interference rejection capability.

Simulation Results
In the simulations, consider that the DDSA contains two uniform linear subarrays with 8 elements spaced half a wavelength, and the gap between the subarrays is h 1 � 16 d. Assume that the desired signal impinges upon the DDSA from the direction θ 0 � −10°, and three interfering signals come from the directions θ 1 � −40°, θ 2 � 30°, and θ 3 � 50°, respectively. e interference-to-noise ratio (INR) is set to 30 dB. e additive noise is presumed to be a complex circularly symmetric Gaussian zero-mean spatially and temporally white process. e number of snapshots is set to N c � 60 while comparing the output SINR versus the input SNR, and the input SNR is fixed at 10 dB in the performance comparison versus the number of snapshots. Note that the random DOAs of the desired signal and the interferers change from trial to trial but remain fixed from snapshot to snapshot. For each scenario, the simulation results are provided as the average of 200 Monte Carlo trials. e proposed adaptive beamformer is compared to the worst-case-based beamformer (RAB-WCPO) [29], distributed LCMV beamforming (DA-LCMV) [17], INCM reconstruction-based beamformer (REC-CMRSVE) [25], and the joint adaptive beamforming (SMIS-JAW) [6]. In addition, the optimal output SINR of the subarray is shown in the simulations to clearly demonstrate the advantages of the proposed beamformer. According to the MVDR criterion, the optimal output SINR is calculated by (5). e parameter ε � 0.3 is used in the RAB-WCPO, and the angular sector of the desired signal is set to be Θ � [θ 0 − 5°, θ 0 + 5°] while suppressing the sidelobe interference. e sampling grid is uniform in Θ and Θ with 0.1°increment between adjacent grid points.

Simulation for Beampattern.
In the first example, the mainlobe width and sidelobe levels of the beampattern are examined. Assume that the DOAs of the desired signal and three interferers are accurately known. As the desired signal components in the sample covariance matrix may degrade the performance of some tested beamformers, the input SNR is set to be 0 dB. In the simulations, the DOA of the first interferer is changed from − 40 ∘ to − 6 ∘ , which can form the mainlobe jamming. Note that, in this scenario, the angular sector of the desired signal is set to be Θ � [− 13 ∘ , −7 ∘ ]. Figures 4(a) and 4(b), respectively, correspond to beampattern with sidelobe and mainlobe interferences. It can be seen from Figure 4(a) that, as a result of employing more elements spaced half a wavelength, the mainlobe width and sidelobe levels of the proposed algorithm are better than those of the other tested beamformers while suppressing the sidelobe interference. Moreover, the proposed beamforming beampattern has deeper nulls with the normalized amplitudes below -120 dB. Figure 4(b) shows that the subarray mainlobe interference is effectively suppressed by the virtual contiguous array, while the other tested beamformers cannot resist the beam-pointing drift while suppressing the mainlobe jamming, including the optimal beamformer. From the beampatterns, the proposed beamformer possesses the obvious advantages in many aspects, such as mainlobe width, sidelobe levels, null depth, and beam direction.

Performance Comparison of Different Arrays.
In the second example, the DDSA with the proposed algorithm is compared to the traditional ULA and nonuniform arrays, including coprime array and nested array. Assume that, for the arrays taken into the comparison, the DOA of the desired signal and INCM can be accurately obtained in advance, and then the optimal beamformer based on the minimum variance distortionless response (MVDR) criterion can be applied to the traditional ULA, coprime array, and nested array, respectively. Note that all the arrays involved in the experiments have 16 elements. e ULA is composed of 16 omnidirectional elements with interspacing d. e coprime array contains two subarrays, whose first element is shared.
e first subarray has 8 elements with interelement spacing 9 d, while the second subarray has 9 elements with interelement spacing 8 d. e nested array is a concatenation of an inner ULA and an outer ULA, and both of them have 8 elements. e spacing between two adjacent elements of the inner ULA is half a wavelength d, and the intersensor  International Journal of Antennas and Propagation spacing of the outer ULA is 9 d. Assume that the desired signal comes from the direction θ 0 � 0 ∘ , and three interfering signals impinge upon the arrays from θ 1 � −40 ∘ , θ 2 � 30 ∘ , and θ 3 � 50 ∘ , respectively. In the case of sidelobe interferences, the beampattern is shown in Figure 5(a), and Figure 5(b) compares the output SINR versus input SNR with the number of snapshots K � 200. Figure 5(a) demonstrates that the mainlobe width of the proposed scheme is smaller than that of the ULA, and the sidelobe levels are much lower than those of the nonuniform arrays. As shown in Figure 5(b), the proposed scheme can achieve the best performance of the other tested arrays. When the DOA of the first interference changes from −40 ∘ to 4 ∘ and other DOAs remain unchanged, the interfering signal in the main lobe may have different impacts on the tested arrays. Note that, in this simulation, the angular sector of the desired signal is set to be Θ � [−3 ∘ , +3 ∘ ]. In the case of mainlobe and sidelobe interferences, Figure 5(c) gives the beampattern, and Figure 5(d) depicts the output SINR versus input SNR with the number of snapshots K � 200. It is shown in Figure 5(c) that the beam pointing of the ULA deviates from the preset direction, while other schemes maintain the correct direction. From Figure 5(d), it is indicated that while suppressing the mainlobe interference, the proposed scheme outperforms the other tested arrays with the optimal beamformer. Simulation results demonstrate that the proposed scheme can achieve or even exceed the optimal beamformers of the traditional ULA, coprime array, and nested array.  Figure 6(a) that if the DOA of the desired signal is known, DA-LCMV, REC-CMRSVE, and the proposed algorithm have similar performance and are superior to other tested beamformers. Figure 6(b) shows that both DA-LCMV and REC-CMRSVE need more snapshots to achieve good performance. e figures demonstrate that the proposed beamformer enjoys a better performance improvement in this scenario. Compared with the subarray optimal beamformer, the proposed beamformer through multiple subarray synthesis significantly outperforms that of single subarray.   Figure 7(a). When the input SNR is set as 10 dB, Figure 7(b) demonstrates the performance curves versus the number of snapshots. It can be seen from the figures that the proposed method dramatically outperforms other tested beamformers, even the optimal beamformer.

Adaptive Interference Rejection
at is because the beam-pointing drift while suppressing the mainlobe jamming, which has been verified in Figure 4(b), seriously degrades the performance of the optimal beamformer. Figure 7(b) shows that, with the increase of the number of snapshots, the performance of the proposed beamformer approaches or even exceeds the optimal beamformer when the number of snapshots is greater than 90. erefore, the proposed algorithm significantly outperforms the other tested methods including the optimal beamformers in the case of mainlobe jamming.  Figure 8(a) compares the output SINR versus the input SNR, and the output SINR against the number of snapshots is displayed by Figure 8(b). Figure 8(a) demonstrates that when the input SNR is greater than 0 dB, the proposed algorithm is very close to the optimal beamformer. Despite the DOA estimation errors at low SNRs make the performance of the proposed beamformer degrade slightly, whose output SINR is close to that of the subarray optimal beamformer, but the output SINR gradually approaches the optimal beamformer when the input SNR is higher than −20 dB. It is seen from Figure 8(b) that, with the increase of the number of snapshots, the output SINR of the proposed beamformer is very close to the optimal beamformer, and the proposed algorithm outperforms the other tested beamformers.

Incoherent Local Scattering.
In the last example, we consider the influence of incoherent local scattering on performance of the tested beamformers. Assume that the desired signal has a time-varying spatial signature as follows [23,31]: where s p (k) ∼ N(0, 1), p � 0, 1, 2, 3, 4 are independently and identically distributed zero-mean complex Gaussian random variables changing from snapshot to snapshot; θ l , l � 1, 2, 3, 4 are independently normally distributed in N(θ 0 , 2°). Note that θ l changes both from trial to trial and from snapshot to snapshot. In this scenario, the signal covariance matrix R s is no longer a rank-one matrix, and the optimal SINR should be obtained by [23,28] which is maximized by where P · { } denotes the principal eigenvector of a matrix. Figures 9(a) and 9(b) display the beamforming performance curves versus the input SNR and the number of snapshots, respectively. It can be seen from Figure 9(a) that when the input SNR is lower than 20 dB, both the proposed algorithm and REC-CMRSVE outperform the other tested methods, but the performance of REC-CMRSVE degrades with the further increase of input SNR. e output SINR of the proposed method has always been very close to that of the optimal beamformer. Figure 9(b) shows that, with the increase of the number of snapshots, the performance of the proposed algorithm approaches that of the optimal beamformer. ere is only a little SINR loss for the proposed beamformer because the DOF has been greatly increased and the interfering signals can be effectively suppressed. erefore, the proposed beamforming algorithm outperforms other tested beamformers in the case of incoherent local scattering.

Conclusion
In this paper, a novel robust adaptive beamforming algorithm for DDSA is proposed by virtual reconstruction. e gaps between the distributed subarrays are filled with the virtual elements, and then a new contiguous array can be formed. By dividing the DDSA into two rotation-invariant subarray sets, the DOAs can be accurately estimated by the ESPRIT-like method, and the received interfering signals on the virtual elements can be successfully reconstructed. At low SNRs, the steering vector of the desired signal can be obtained by solving the optimization problem. Based on the power estimation of the interfering signals, the INCM of the virtual contiguous array can be reconstructed, and then the weight vector of the proposed beamformer can be obtained. Without increasing the number of actual elements, the proposed algorithm greatly improves the DOF to enhance the interference rejection capability. Note that the weight vector must be applied to the rearranged data matrix received by the actual and virtual elements. With the virtual elements forming a contiguous array, the proposed beamforming algorithm can effectively suppress the grating lobes, reduce the mainlobe width, and degrade sidelobe levels. Simulation results highlight the validity and effectiveness of the proposed algorithm.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare no conflicts of interest.

Authors' Contributions
Jian Lu and Jian Yang contributed to conceptualization and methodology; Jian Lu and Bo Hou were responsible for software, validation, and investigation; Jian Lu contributed to formal analysis; Jian Yang and Zhiyong Yu contributed to data curation; Jian Lu and Zhiyong Yu prepared the original draft ; Jian Lu, Jian Yang, and Guangbin Liu reviewed and edited the manuscript; Zhiyong Yu contributed to visualization; Guangbin Liu supervised the study; and Jian Yang contributed to project administration and funding acquisition.