Broadband DOA Estimation Based on Nested Arrays

Direction of arrival (DOA) estimation is a crucial problem in electronic reconnaissance. A novel broadband DOA estimation method utilizing nested arrays is devised in this paper, which is capable of estimating the frequencies and DOAs of multiple narrowband signals in broadbands, even though they may have different carrier frequencies. The proposed method converts the DOA estimation of multiple signals with different frequencies into the spatial frequency estimation. Then, the DOAs and frequencies are pair matched by sparse recovery. It is possible to significantly increase the degrees of freedom (DOF) with the nested arrays and the number of sources can bemore than that of sensor array. In addition, the method can achieve high estimation precision without the two-dimensional search process in frequency and angle domain. The validity of the proposed method is verified by theoretic analysis and simulation results.


Introduction
Direction of arrival (DOA) has been an active research area, playing an important role in electronic reconnaissance [1], as it is a crucial parameter for sorting and recognition of sources, directing jamming and passive location.Unlike radar receiver, there is no priori information of signals for the reconnaissance receiver.Therefore, the reconnaissance receiver tends to have large instantaneous bandwidth in order to cover a wide spectrum [2][3][4].This leads to two problems in DOA estimation: one is that the array aperture must be increased to satisfy multiple signals of different frequencies within the reconnaissance band range and the other is that it must be possible to achieve the frequency and DOA estimation of multiple signals simultaneously in broadband.In [5], a nested array structure was proposed to increase the degrees of freedom by vectorizing the covariance matrix of the received signals among different sensors.Then, Pal and Vaidyanathan increase degrees of freedom of the coarray by extending the nesting strategy to multiple levels in [6], namely, the 2thorder nested array, whose 2th-order difference coarray is proved to contain a uniform linear array with ( 2 ) sensors.It can be viewed as a virtual array with a wider aperture, which is capable of being utilized to improve the DOA estimation performance of the multiple signals with same frequency.However, when the frequencies of the signals are different and unknown, this technique cannot be applied directly.By exploiting the four-level nested array, a novel estimation approach is proposed in [7], which can be extended to wideband scenarios.However, for each subband signal, the vectorization is performed twice to construct a fourth order difference coarray and the spectrum peak search process is conducted in the whole spatial domain, which suffers from the complex computation.Signal subspace techniques such as MUSIC [8] need two-dimensional search process in frequency and angle domain, which leads to large amount of computation.In [9], an angle-frequency joint dictionary is established to achieve the frequency and DOA estimation by sparse recovery.The results suggest that sparse recovery becomes more difficult and complicated while the length of dictionary increases.
A broadband DOA estimation method is investigated in this paper.A nested array is used to extend the array aperture and the spatial frequency, which contains the information of both frequency and angle, is defined.In this way, the DOA estimation of multiple signals with different frequencies is converted into the spatial frequency estimation.Then, the received signals of sensor array are transformed into 2 International Journal of Antennas and Propagation frequency domain.By sparse recovery approach, each of the signals with different frequency can find its corresponding spatial frequency in the dictionary composed of spatial frequencies that have been estimated.Once the frequency and spatial frequency are pair matched, DOA estimation can be implemented simply by the definition of spatial frequency.This paper is organized as follows.In Section 2, we introduce the nested array model.In Section 3, we propose the broadband DOA estimation method based on nested array and analyze the performance.Section 4 presents the simulation results, and conclusions are made in Section 5.
Assume that there are  sources imping on the array from directions  1 ,  2 , . . .,   with carrier frequency  1 ,  2 , . . .,   , respectively.The steering vector a(  ,   ) corresponding to the th source signal is expressed as where  is the velocity of wave propagation.The data vector received at the nested array is expressed as where where A denotes the array manifold matrix and s() denotes the source signal vector.n() is the white noise, which is uncorrelated from the sources.

The Nested Array Model
3.1.Spatial Frequency Estimation.From ( 1)-( 4), it is obvious that there is two-dimensional information in broadband DOA estimation, namely, frequency and angle.Therefore, we attempt to estimate spatial frequency instead of frequency and angle, and the spatial frequency Ω is defined as Combining ( 1) and ( 3), we obtain where Here, we assume that the sources are mutually uncorrelated and each source signal   () is wide-sense quasistationary [11] with frame length ; that is, Then, the autocorrelation matrix of the th frame received signal is where ) is the source covariance matrix at frame . 2  is the power of the th source and  2  is the power of noise.Vectorizing R  , we get the following  2 × 1 vector: where T with e T 1 being a column vector of all zeros except a 1 at the th position.Comparing it with (2), we can say that y  in (10) behaves like the data vector received at an array whose manifold is given by A * ⊙ A, where ⊙ denotes the Khatri-Rao (KR) product.
By stacking [y 1 , y 2 , . . ., y  ] ≜ Y, we have that where The noise subspaces U  can be obtained by applying singular value decomposition (SVD) on Y, which satisfies that Then, the conventional subspace-based DOA estimation approach as MUSIC can be exploited as follows: It is interesting to note that spectrum peak search process is applied in the spatial frequency domain instead of the angle domain.The source signals with different frequencies and angles correspond to respective spatial frequency; therefore they can be separated in the spatial frequency domain.Since the spatial frequency includes both frequency and angle information, we will show how to determine their matching relation in the following part.

Pair Matching.
According to the definition of spatial frequency in ( 5), once we know the frequency of each signal and match it to the corresponding spatial frequency, DOA can be obtained by simple computation.Hence, the pair matching of frequency and angle is equivalent to that of frequency and spatial frequency.
Consider the frequency domain model of ( 2), which is given by where X(  ) is the received data at the array with frequency   , which is achieved by the Fourier transform of x().Our pair matching process is targeted at finding one or more spatial frequencies corresponding to X(  ) among those that have been already estimated.We can treat this problem as a sparse recovery problem [12]; that is, where A dict is an overcomplete dictionary composed of the steering vector generated by the estimated values of spatial frequency.v is the sparse vector whose position of nonzero elements represents the spatial frequency matched to X(  ).However, the  0 -minimization is NP-hard in general.Thus, it is often converted into  1 -minimization [13] as follows: The support of k can be obtained by convention sparse recovery algorithms, such as Basis Pursuit (BP) [14] and Orthogonal Matching Pursuit (OMP) [15].Considering the scenario that some of the source signals are with the same frequency, the sparse recovery problem in the case of high sparsity ratio may exist.Kwon et al. proposed a multipath matching pursuit (MMP) algorithm [16], which has been proved effective in high sparsity ratio case.Therefore, the MMP sparse recovery algorithm is selected.
As it is described above, it is possible to find the matched spatial frequency in the dictionary A dict for each X(  ) with the MMP algorithm.
where f()  denotes the frequency corresponding to Ω .As θ is obtained by the matched f and Ω , θ and f are pair matched automatically.

Performance Analysis.
The implementation steps of the proposed method are summarized as follows.
Step 2. Perform SVD on Y in order to obtain the noise subspace matrix U  .Then, compute the spatial frequency spectrum by (13) and obtain Ω , for  = 1, 2, . . ., .
Step 4. Compute the support of X(  ) in A Dict by sparse recovery and determine the corresponding Ω to f .
Step 5. Compute   by (18) with the matched Ω and f .
Step 6. Repeat Steps 4 and 5 for each X(  ) and obtain the set { θ }  =1 as the DOA estimation.
It is obvious that the performance of DOA estimation is related to the spatial frequency and frequency.The precision of frequency estimation is due to discrete Fourier transform (DFT) points   , whose maximum value is the number of snapshot.From (18), we have This means that the influence of frequency estimation error on the DOA estimation decreases with the growth of frequency.Therefore, the DOA estimation precision is mainly International Journal of Antennas and Propagation due to the spatial frequency estimation precision when the carrier frequency is high.As the nested array, which is capable of extending the array aperture effectively, is utilized in spatial frequency estimation, we can achieve high estimation precision compared to the ULA.If the spatial frequency is estimated accurately, the DOA estimation with high precision follows.
In fact, the broadband DOA estimation problem includes two-dimensional parameter: frequency and angle.It is normally implemented with a two-dimensional search process in frequency and angle domain.Compared with traditional method, the proposed method only needs to conduct spectral searching process once in the spatial frequency domain and a simple pair matching process.However, the traditional method, which uses band-pass filters for each element and obtains narrowband DOA estimation by MUSIC, needs multiple search processes for the signals with different frequencies.Therefore, the proposed method reduces the computational complexity to some extent.

Simulation Results
In this section, we provide several sets of simulation results to demonstrate the performance of the proposed method.In all the simulation examples below, we suppose that  ∈ (0.1, 1.5) GHz,  ∈ (−90 ∘ , 90 ∘ ), for all the source signals.The root mean square error (RMSE) of angle from  Monte Carlo trials is used as our performance measure and it is defined as where   and θ denote the true and estimation DOAs, respectively.

Simulation and Setting
Example 1.This simulation example considers an underdetermined case where 13 uncorrelated source signals ( = 13) impinge on the nested array with 6 sensors ( 1 =  2 = 3).The number of snapshots  ( = ) and SNR are set to  ×  ( = 32,  = 1024) and 15 dB, respectively.The true DOAs are given by  2 and 3 show the spatial frequency and frequency estimated results, respectively, and Figure 4 shows the pair matching results of frequency and angle by MMP algorithm.
Example 2. In this simulation example, we consider the spatial frequency estimation performance with respect to SNR.Five source signal angles and frequencies are set to {1.0, 0.9, 1.1, 1.2, 1.3} GHz and {−30 ∘ , −15 ∘ , 5 ∘ , 15 ∘ , 30 ∘ }, respectively.The SNR varies from −10 dB to 15 dB in 2.5 dB intervals.Figure 5 shows the RMSE of the spatial frequency estimation as a function of SNR, averaged over 1000 Monte Carlo trials.the SNR and snapshots.The SNR varies from −5 dB to 20 dB in 2.5 dB intervals.The RMSE of angle as a function of SNR with different snapshots is shown in Figure 7.
Example 5.In this example, we compare with traditional method, which uses band-pass filters for each element and obtains narrowband DOA estimation by MUSIC.The SNR varies from −5 dB to 15 dB in 2.5 dB intervals.We choose five source signals with different frequencies and DOAs.Every fixed SNR conducts 1000 Monte Carlo trials.Figure 8 compares the RMSE of the proposed method and MUSIC method based on band-pass filtering.The CPU runtime of MUSIC and the proposed method as a function of source number are shown in Figure 9.  2 and 3, we see that the proposed method is effective for spatial frequency and frequency estimation.In Figure 4, e and I denote the true values and estimated values, respectively, which prove that the frequency and angle can be pair matched correctly by MMP algorithm.

Discussion and Analysis. From Figures
The RootCRB() in Figure 5 denotes the square root of Cramer-Rao bound with  sensors.It is interesting to note that the RootCRB is for spatial frequency instead of angle.Since the nested array increases the DOF from  to ( 2 − 2)/2 +  ( is even) [4], we compare it with the RootCRB (6) and RootCRB(23).From Figure 5, we observe that the nested array with 6 sensors yields better RMSE of spatial frequency than the ULA with 6 sensors, but still not as good as the ULA with 23 sensors.
The results of Figure 6 show that the smaller   is, the better RMSE we can obtain.However, with the same   , the RMSE of angle will decrease when the frequency varies from 0.1 GHz to 1.5 GHz.Therefore, we can conclude that the frequency estimation error has less influence at higher frequencies, and vice versa for lower frequencies, which is equal to the analysis in (19).
From Figure 7, we see that better RMSE of angle can be obtained by increasing  with fixed  or increasing  with fixed , as it increases the number of snapshots  =  in both scenarios.However, when the number of snapshots is fixed at  = 16384, the case of  = 512 and  = 32 provides better RMSE than that of  = 1024 and  = 16.Hence, this set of empirical results suggests that, for a fixed , trying to obtain more frames  by decreasing the frame length  ( should not be overly small) tends to be a better option for performance improvement than the opposite.
From Figure 8, we see that the RMSE of the proposed method is better than MUSIC.This is because it extends the array aperture by nested array, which improves the DOA estimation precision.The results of Figure 9 show that CPU time of the proposed method is less than MUSIC and varies little with the growth of source number.The reason is that the proposed method can estimate the spatial frequencies of all the signals in a single spectral search process, which is less affected by source number.However, the times of spectral search in MUSIC are related to the source number, which equals the number of signals with different frequencies.Therefore, the CPU runtime becomes longer significantly when the source number is large.

Conclusion
This paper has addressed the broadband DOA estimation problem with nested array.The proposed method can obtain the frequency and angle estimation of multiple source signals by spatial frequency estimation and pair matching with sparse recovery.The pair matching process is conducted in the set of spatial frequency that have been estimated for reducing the computation.Moreover, the proposed method can work in broadband, underdetermined, and low SNR cases, which is very suitable for electronic reconnaissance system.

3. 3 .
Frequency and DOA Estimation.The set of { f }  =1 achieved by Fourier transform on the rows of x() is regarded as the estimation of frequency.Since the pair matching for the set of { f }  =1 and { Ω }  =1 has been already implemented, the DOA estimation is given by

Example 4 .
In this example, we examine the DOA estimation performance with respect to SNR and snapshots.The simulation conditions are similar to those of Example 1 except
Figure 1: A 2-level nested array with  1 sensors in the inner ULA and  2 sensors in the outer ULA, with spacings  1 and  2 , respectively.