MT-BCS-Based DoA and Bandwidth Estimation of Unknown Signals through Multiple Snapshots Data

,


Introduction
Direction-of-Arrival (DoA) estimation has important applications in many traditional fields such as smart antenna array, radar signal processing, geophysical or seismic sensing, acoustics, multiple-input multiple-output (MIMO), and other applications related to finding the direction of the incoming signals or sources [1][2][3][4].In fact, the DoAs knowledge of receiving signal can be properly used to locate the corresponding source, and the adaptive beamforming of the receiving antenna pattern can improve the sensitivity of the system to the desired signal direction or suppress unnecessary interference.
In some scientific literature, several effective methodologies have been proposed about DoA estimation, such as the multiple signal classification (MUSIC) [5] and the signal estimation parameter via rotational invariance technique (ESPRIT) [6].For the limitations that the need of a priori knowing, the number of incoming signals and the calculation of the correlation matrix are usually computationally expensive, and their standard implementations are rarely available especially nowadays with the huge proliferation of wireless devices/services and the presence of noncollaborative users [7].More recently, DoA estimation strategies based on the Compressive Sensing (CS) [8,9] are proposed and have shown promising features and results.CS is an enabling paradigm for many applications where there is the need of overcoming the Shannon's limit in data acquisition and to recover sparse signals from far fewer measurements [10].
anks to these features that the computational efficiency, the accuracy, and the robustness to the noise, CS-based strategies have already been applied to a variety of applications in electromagnetic engineering [11,12].But for guaranteeing reliable estimations, the sampling matrix must satisfy the restricted isometry property (RIP) [13] when applying CS.Unfortunately, such a condition cannot be easily verified because of it results are computationally demanding [14].But innovative approaches using the Bayesian compressive sensing (BCS) [15,16] have been proposed, and they have improved this problem.However, most of the works about DoA estimation are mainly focused on narrowband signals.In [17], the narrowband algorithm is simply extended to broadband.As for the DoAs estimation of wideband signals, because of the guiding vector of broadband signal frequency, the signal space of different sub-bands is inconsistent, which makes it difficult to separate the signal space from noise space.General, subspace-based estimation techniques [18] based on the root-MUSIC [19] and maximum likelihood [20] are used to estimate the DoAs of the wideband signal, but they also have limitations.
is problem has been shown in [21,22], where the MT-BCS [23] has been customized to deal with wideband signals while exploiting the correlation among different frequency samples taken from single snapshot data and multiple snapshots, respectively.
In this paper, the DoAs and BW estimation problem of wideband signals are formulated within the MT-BCS framework based on Laplace priors.Starting from the key observation that the wideband signal affecting the antenna array is essentially sparse in spatial domain.Because the signal is a wideband signal, we sample the signal with a certain range of frequencies in order to estimate the bandwidth of the signal.Set the sampling frequency range to [f min , f max ], the DoAs and BW are estimated at the same time by which the range of frequencies has measured data.e rest of the paper is organized as follows.e Bayesian modeling is mathematically formulated in Section 2, where the wideband signal model of MT-BCS using Laplace priors is described.A set of representative numerical results is then reported and discussed in Section 3, where reference DoAs and bandwidth estimation methods are performed.Finally, some conclusions are drawn in Section 4.

Bayesian Modeling
CS theory can cover certain signals from far fewer samples or measurements than traditional methods, and it is a very important step to restore the signal we want.In some work, narrowband signals are recovered by MT-BCS [24,25] e input data of the problem are the voltages v (p) n , n � 1, . . ., N, measured by each sensor expressed as where s k , k � 1, . . ., K, is the amplitude of the impinging signals, h is the antenna effective length, λ (p) being the free-space length at pth frequency, and η (p) n , n � 1, . . ., N, the contribution of the noise at each sensor.
en, in order to employ the proposed MTBCS-based methodology, the problem is reformulated by sampling the angular domain of interest over a very fine grid of I >> K angular location, being θ i � − (π/2) + π(i − 1)/(I − 1), i � 1, 2, . . ., I.Under this assumption and considering the multiple-snapshots case, (1) can be rewritten in a matrix form as where L is the number of snapshots, s M×N , is the measurement matrix, we will use the Gaussian matrix in this paper, and M (M << I) is the number of samples [16].A (p) is the matrix of array manifold: where x n � (n − 1)d, n � 1, . . ., N, is the distance between the ith array element and the reference array element and d is the array element spacing.as a Laplace priors about parameter λ, and it is defined as

Model for MT-BCS
And the mean and covariance of signals are given by where Λ � diag(1/c i ), i � 1, . . ., I, c i , is the signal hyperparameter and estimated as with 2 International Journal of Antennas and Propagation where ϕ i , i � 1, . . ., I, expresses the ith column of Φ, β is the noise hyperparameter and Gamma prior obeying parameters a and b.Parameter λ realizes the following Gamma hyperprior: e parameters estimation of λ, β, and υ is expressed as following: i > λ and c i � 0 then add c i to the model 7: Else if q 2 i > λ and c i > 0 then re-estimate c i using (7) 8: Else if q 2 i < λ then prune ith column of Φ from the model and set c i � 0 End if 9: Update μ, L, Σ using ( 5) and (6) 10: Update s i , q i using ( 8 From the above formula, we can see that it is known that c, λ, β, and ] can find μ and Σ using ( 5) and (6).It is known that μ and Σ can also find c, λ, β, and ] using (7)∼( 15).e procedure is summarized in Algorithm 1.

Numerical Results
In this section, we present experimental results that demonstrate the performance of MT-BCS based on Laplace priors to recovery wideband signals.Different from narrow narrowband signals, the wideband signal exists only when sampling frequency f p is within the range of bandwidth.
We use the following default setup in the experimental results reported in this section.Signal matrix of length I is generated, where K coefficients are located at signal randomly, and the rest (I-K) of the coefficients are set equal to zero.As the measurement matrix Ψ, we chose a Gaussian matrix, where the columns ψ i are Gaussian distributed on the sphere R I .Moreover, we present results with noisy where  s p l are the estimated signal energy in lth snapshot at pth frequency sampling.Moreover, the angular grid has been discretized with I � 181 samples in order to obtain a resolution of Δθ � 1 deg.
e actual DoAs and those estimated values by the proposed spectral correlation MT-BCS based strategy are shown in Figure 1 and Table 1.As it can be observed in Table 1, there are some values that are small enough to be ignored.is is mainly caused by noise.We ignore these values below 0.  1.
e actual and estimated bandwidths of the K � 3 impinging signals are shown in Figure 2. It is possible to observe that both the directions of the incoming signals are correctly retrieved (Figure 1) as well as their band widths.
en, we will see reconstruction error by RMSE.e equation of RMSE is expressed as follows: where T is number of experiments and  θ f p ,k,t,l is the estimated value of kth signal at pth frequency.
In     is will have an impact on the accuracy of DoA estimation results.
In order to obtain better reconstruction performances, we will add signal energy  s p avg in (16) at different frequencies of the same signal and average it, and then get  s avg .( 16) can be rewritten as follows: In Figure 4, the reconstruction DoA results using (18) at different snapshots are shown.We have also given estimated DoAs and values at pth frequency, as shown in Table 2.  s avg is the average signal energy of all estimated angles from low to high.e signal energy that shows the unrelated angle is much lower than what we get in the first set experiments. is shows that using (18) to estimate DoA can get better reconstruction results than using (16).We also can note that as the number of snapshots increases, unrelated angles no longer exist.Frequency away from f 0 still has poor reconstruction results, but this is an improvement after the signal energy of all frequencies is averaged, that is, the signal energy of the unrelated angle becomes smaller.
We ignore values less than 0.1 as in the first set of experiments.It can be concluded that the DoAs and bandwidths of signals form the result shown in Figure 4 and Tables 2 and 3

Conclusion
In this paper, we formulated the DoAs and BW-estimated problem of the wideband signal using MT-BCS based on Laplace priors and presented a framework for multiple snapshots data.Using this framework, we first get the actual and estimated DoAs for multiple snapshots at different frequencies.At the same time, the signal BW is obtained by observing at which frequency sampling point there is a signal energy after ignoring values that are small enough.en, we use the RMSE to evaluate the reconstruction results in different snapshots and SNRs, finding with the increasing of snapshots or SNR, the reconstruction results are better.International Journal of Antennas and Propagation )∼(11) 11: Update λ, ], β using (12)∼(15) 12: end while ALGORITHM 1: Algorithm of MT-BCS.

Figure 1 :
Figure 1: Actual and estimated DoAs for multiple snapshots.

Figure 2 :
Figure 2: Actual and estimated signals bandwidth.

1 .
It can be observed that DoA θ 1 � − 9 deg and bandwidth b 1 � [f 1 : f 3 ] for s 1 , θ 2 � − 5 deg and b 2 � [f 1 : f 4 ] for s 2 , and θ 3 � 60 deg and b 3 � [f 1 : f 5 ] for s 3 , which is shown Figure 1 and Table the second set of experiments, let us consider K � 2 signals, designed to work in the frequency range [f min , f max ] � [0.6: 1.0] GHz.DoAs equal to θ 1 � − 55 deg and θ 2 � − 32 deg, and in order to facilitate the use of RMSE calculation error, we set the two signals which have the same bandwidths [0.6 : 1.0] GHz. e measured data are sampled at P � 5 frequencies (f 1 , f 2 , f 3 , f 4 , f 5 ) � (0.6, 0.7, 0.8, 0.9, 1.0) GHz, being Δf � 0.1 GHz.Average reconstruction errors of T � 500 runs are shown for the case in Figure 3 for all types of signals.It is noted that RMSE results in lower values than setting snapshots minor from Figure 3(b) and the frequency close to f 0 show better performance than others from Figure 3(a).

Figure 5 :
Figure 5: Actual and estimated signal bandwidth.

Table 1 :
Estimated DoAs and signal energy for different frequencies.

Table 3
International Journal of Antennas and Propagation e results of the same experimental setup with different SNR are shown in Figures 3(c) and 3(d).Similar performance decrease by the increase of SNR can be observed.ese suggest that the estimating DoA results can get lower reconstruction error in multiple snapshots and higher SNR.For wideband signals, different sampling frequencies can obtain different reconstruction results.