Next Article in Journal
Development of an On-Board Measurement System for Railway Vehicle Wheel Flange Wear
Next Article in Special Issue
A Novel Unitary ESPRIT Algorithm for Monostatic FDA-MIMO Radar
Previous Article in Journal
Channel Prediction Based on BP Neural Network for Backscatter Communication Networks
Previous Article in Special Issue
SPICE-ML Algorithm for Direction-of-Arrival Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Sparse Bayesian Learning-Based Off-Grid DOA Estimation Method for Vehicle Localization

Electronic Information School, Wuhan University, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(1), 302; https://doi.org/10.3390/s20010302
Submission received: 28 October 2019 / Revised: 26 December 2019 / Accepted: 1 January 2020 / Published: 5 January 2020

Abstract

:
With the rapid development of the Internet of Things (IoT), autonomous vehicles have been receiving more and more attention because they own many advantages compared with traditional vehicles. A robust and accurate vehicle localization system is critical to the safety and the efficiency of autonomous vehicles. The global positioning system (GPS) has been widely applied to the vehicle localization systems. However, the accuracy and the reliability of GPS have suffered in some scenarios. In this paper, we present a robust and accurate vehicle localization system consisting of a bistatic passive radar, in which the performance of localization is solely dependent on the accuracy of the proposed off-grid direction of arrival (DOA) estimation algorithm. Under the framework of sparse Bayesian learning (SBL), the source powers and the noise variance are estimated by a fast evidence maximization method, and the off-grid gap is effectively handled by an advanced grid refining strategy. Simulation results show that the proposed method exhibits better performance than the existing sparse signal representation-based algorithms, and performs well in the vehicle localization system.

1. Introduction

The emergence of the Internet of Things (IoT) and its applications have greatly affected the society and improved the quality of life [1,2]. For example, the IoT has been widely used in the autonomous vehicles [3,4], which are viewed as a potential solution to some critical problems, such as pollution, traffic congestion, and road accidents [5]. Almost all autonomous vehicles have five basic functions, including localization, perception, control, planning, and system management [6]. In order to implement correct driving decisions and actions, accurate knowledge of the vehicle’s location is required for the perception, planning, and control functional systems. Therefore, vehicle localization technology is critically important for autonomous vehicles. As the global positioning system (GPS) can offer a cheap and easily accessible solution for global localization, it is the first choice for vehicle localization. However, the accuracy and availability of the GPS signal cannot always meet the requirements of autonomous vehicle localization, due to satellite visibility interruption and signal multipath. To address this issue, many vehicle localization systems based on robust location algorithms with advanced sensors, such as RADAR [7], LiDAR [8,9], and camera [10], have been invented and become the focus of research.
In general, vehicle localization systems can work based on four principles, namely, time of arrival (TOA), time difference of arrival (TDOA), radio signal strength (RSS), and direction of arrival (DOA). In TOA approaches [11,12], a relative position to the reference point is obtained based on the time from a signal being sent to the reference point to the response signal back at the vehicle. The accuracy of the TOA approaches is susceptive to a synchronous clock. As a result, the perfect synchronization is hard to realize in practice. For TDOA approaches [13,14], the relative position of the vehicle can be calculated based on the time difference of signal arrivals. However, the accuracy of TDOA methods is limited due to their high sensitivity to measurements of time difference. In RSS methods [15,16], the attenuation of signal strength is measured to estimate the travel distance of the signal. However, the spatial fading characteristic of the signal is firstly required. In DOA approaches [17,18,19], the angle of the signal is measured via antenna arrays and the position can be calculated by the angle information. The accuracy of DOA approaches is only dependent on the accuracy of the angle estimation, which can be easily realized by the existing high accurate DOA estimation algorithms. DOA is thus the most efficient one among the four approaches.
Various DOA estimation methods have emerged during the past decades, among which subspace-based algorithms are best known with a wide range of applications [20]. However, they suffer from certain well-known constraints: prior knowledge of source number, uncorrelated circumstance, and a sufficient number of snapshots are required. The rapid development of sparse representation and its recent application in various fields [21,22,23,24,25,26,27,28,29,30,31,32] has offered a new perspective for DOA estimation. A great number of sparse representation-based DOA methods have emerged, which can be applied in several demanding scenarios that the subspace-based methods cannot. The recent DOA estimation based on sparse representation can be summarized into three categories. The first category concerns on-grid methods, including l 1 -SVD [33] , SPICE [34], SBL [35], and MSBL [36]. For these kinds of methods, true DOAs are assumed to lie on a set of fixed grid points, therefore the existing sparse representation techniques can be applied directly. However, the gap between true DOA and its nearest grid point always exists in practice. The other two categories, namely, gridless methods and off-grid methods, can eliminate or narrow this gap. Gridless methods [37,38,39,40,41,42] operate in the continuous domain directly so that they can avoid the grid mismatch problem. Gridless methods have strong theoretical guarantees and can only be applied to the uniform or sparse linear arrays.
For the off-grid methods, a sampling grid is still required but the true DOAs are not restricted to be on the grid [43]. Han et al. (2015) [44] introduce a sparse total least-squares (STLS) approach, which can yield a maximum a posteriori (MAP) optimal estimate if the matrix perturbation caused by the basis mismatch is Gaussian. Off-grid sparse Bayesian inference (OGSBI) [45] takes a Bayesian perspective on the off-grid methods. The model of OGSBI is based on a first-order Taylor series expansion, and a Laplace prior is assumed to exploit the spatial sparsity of signals. Wu et al. (2016) [46] present an off-grid model based on a perturbed sparse Bayesian learning, in which a linear interpolation between two adjacent grid points is adopted. The off-grid methods mentioned above are still limited in such a way: On the one hand, a denser grid set is required for higher accuracy, which, however, may cause massive computation. On the other hand, a coarser grid can greatly reduce the computational workload but may introduce larger model errors. In order to address this issue, Dai et al. (2017) [47] give a root SBL (RSBL) method, which decreases computational workload by using a root method. It also maintains high accuracy with a coarse grid. However, RSBL is still quite time consuming at other steps in the iterative updating. In GRDOA [48], the grid points get reconfigurated via a grid refining and fission strategy, which makes the method more efficient and accurate than RSBL.
Vehicle localization based on DOA estimation is suitable and efficient in practice [17,18,49]. Since many recent DOA estimation algorithms are very fast, the moving distance of a vehicle during algorithms execution can be ignored. In other words, the moving vehicle can be assumed to be motionless in DOA estimation-based vehicle localization. The location of the vehicle can be obtained by the cross localization which merely depends on the accuracy of DOA estimation. Therefore, DOA estimation provides a suitable and efficient solution to the localization of a moving vehicle. For example, [17] provides a vehicle localization model based on two cooperating Multi-input multi-output (MIMO) radars and a cross localization method with a real-valued covariance vector-based SBL DOA estimation algorithm. In [18], three collaborative MIMO radars are used for vehicle localization. A robust SBL based DOA estimation approach is proposed to address the non-uniform noise as well as the off-grid issue. The location of the vehicle can be obtained by averaging the cross-localization results from each pair of the three MIMO radars.
The radar-based vehicle localization systems mentioned above [17,18,49] all use the active radars. In recent years, the passive radar systems have been considered to be promising due to their advantages over active radars, such as low-cost surveillance, reduced pollution of the electromagnetic environment, as well as covert operation [50,51,52,53,54]. The DOA estimation for the bistatic passive radar is investigated in [50], where a joint sparse Bayesian model is established to combine the measurements from both stations under the framework of SBL. Colone et al. (2012) [51] analyze the practical feasibility of a WiFi transmissions-based passive bistatic radar and demonstrates the potential of passive bistatic radar for local area surveillance applications, where vehicles and people can be detected and tracked. More ways of detecting and tracking vehicular and human targets based on passive bistatic radar systems can be found in [52,53,54]. However, to our best knowledge, vehicle localization systems using passive bistatic radar with SBL-based DOA estimation have rarely been studied in the existing literature, despite their high usefulness in practice.
A vehicle localization system based on a passive bistatic radar is low cost, low power usage, and reduced pollution of the electromagnetic environment, since passive bistatic radars have no transmitting stations and only receive data from existing source emitters [51]. The system is also robust as it is not affected by the weather and light. In addition, many existing SBL based DOA estimation methods can achieve higher accuracy than conventional ones [45,46,47,48]; the SBL-based DOA method is thus an efficient scheme for vehicle localization, as the localization method is merely dependent on the accuracy of DOA estimation. Therefore, it is worthwhile to investigate the vehicle localization systems using passive bistatic radar with an advanced SBL based DOA estimation algorithm.
In general, the main novelties and contributions of this paper can be summarized as follows:
(1)
We provide a vehicle localization system consisting of a passive bistatic radar and the vehicle localization method via DOA estimation. The proposed system and method is practically feasible and valuable since it can achieve high efficiency, great robustness, low cost, as well as reduced electromagnetic pollution.
(2)
Under the framework of SBL, we propose a novel off-grid DOA estimation approach. We introduce a fast evidence maximization method to update the hyperparameters to the off-grid DOA estimation domain, which can achieve faster convergence than the expectation maximization method used in many other SBL-based off-grid DOA approaches.
(3)
Based on the conception of grid refining, we propose a new grid refining strategy consisting of a root process and a step process, which effectively reduces the off-grid errors. The grid refining strategy is less sensitive to the location of grid points and DOAs, which makes the proposed method more accurate than other DOA methods. The improvement is strengthened in good conditions with a large number of snapshots and high SNR. Furthermore, the grid refining strategy greatly speeds up the procedure of grid refining, making the proposed method more effective than other DOA methods.
(4)
We conduct adequate simulations to verify the superiority and effectiveness of the proposed method. Then, the influence of the parameter θ s t e p , the number of active grid points, and the different conditions on DOA estimation are examined. Then, we compare the proposed DOA method with GRDOA. Finally, a vehicle localization simulation is conducted to verify the feasibility of the proposed vehicle localization method.

2. Signal Model

2.1. System Model

Consider a vehicle localization system with a passive bistatic radar consisting of two stations that only receive data from source emitters, as shown in Figure 1. Point A and Point B are assumed to be the two receiving stations of the passive bistatic radar, which are placed on the same plane. The localization of a vehicle can be calculated via two DOAs, which can be estimated by these two stations, respectively. As shown in Figure 1, the location of one station is set as the reference point A ( 0 , 0 ) and another station is set as B ( b , 0 ) . V ( x , y ) is the coordinate of the target vehicle, and its azimuth angles related to A and B are set as θ 1 and θ 2 , respectively, which can be obtained via a DOA estimation method. It is easy to know that
tan θ 1 = x y ,
tan θ 2 = b x y ,
and, based on Equations (1) and (2), the location of the target vehicle can be calculated as
x = b tan θ 1 tan θ 1 + tan θ 2 ,
y = b tan θ 1 + tan θ 2 .
In the vehicle localization system, two stations are assumed to be configured and work in the same way, meaning that only one station needs to be described for DOA estimation because the other one works in the same way. Since the proposed DOA method is very fast, the moving distance of the vehicle during the method execution can be ignored. Namely, the moving vehicle can be assumed to be motionless.

2.2. Data Model

A uniform linear array (ULA) with M sensors and half-wavelength spacing is considered in each station of the vehicle localization system. Assume that there are K target vehicles in the same plane range, θ k represents the DOA of the kth target where k = 1 , 2 , , K , with K < M . In a passive radar system, the illuminating signal is emitted by a transmitter that is independent of the radar itself. The echo signal reflected by the target vehicles is considered as narrow-band and far-field. The signal received at the output of the array can be written as
y ( t ) = A s ( t ) + e ( t ) , t ( t 1 , t 2 , , t L ) ,
where y ( t ) = [ y 1 ( t ) , y 2 ( t ) , , y M ( t ) ] T , s ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s K ( t ) ] T , ( · ) T is the transpose, L is the number of snapshots, A = [ a ( θ 1 ) , a ( θ 2 ) , , a ( θ K ) ] is an M × K matrix of steering vectors with a ( θ k ) = [ 1 , v θ k , , v θ k M 1 ) ] T , v θ k = e j 2 π d / λ sin ( θ k ) , d is the distance between adjacent sensors, and λ is the wavelength of the source. e ( t ) = [ e 1 ( t ) , e 2 ( t ) , , e M ( t ) ] T is an unknown noise vector. The signals and the noise are assumed to be zero-mean, stationary complex Gaussian random processes and independent of each other. Furthermore, the noise is considered to be uncorrelated from sensor-to-sensor with common variance σ 2 . Equation (5) can be compactly written as
Y = A S + E ,
with the definitions of Y = [ y ( t 1 ) , y ( t 2 ) , , y ( t L ) ] , S = [ s ( t 1 ) , s ( t 2 ) , , s ( t L ) ] , and E = [ e ( t 1 ) , e ( t 2 ) , , e ( t L ) ] .
To cast the DOA estimation as a sparse representation problem, the sparse signal model is constructed. Uniform sampling over DOA range is used in the conventional off-grid methods. Let θ ˜ = [ θ ˜ 1 , θ ˜ 2 , , θ ˜ N ] be a fixed sampling grid in the range [ π 2 , π 2 ] , where N denotes the grid number ( N M ) and r = θ ˜ 2 θ ˜ 1 denotes the grid interval. If the grid is fine enough, the true DOAs will lie on (or, practically, close to) the grid. An M × N overcomplete matrix A = [ a ( θ ˜ 1 ) , a ( θ ˜ 2 ) , , a ( θ ˜ N ) ] is constructed to represent the steering vectors of the potential DOAs where a signal may or may not be present. An N × L matrix X is set as a zero-padded extension of S whose non-zero rows correspond to the true DOAs θ k , k = 1 , 2 , , K . In the presence of few stationary signal sources, X is row-sparse with K N . Then, the data model can be written as
Y = A ˜ X + E .
As the observations Y and the overcomplete matrix A ˜ are given, our purpose is to recover the row-sparse matrix X . It is a sparse signal recovery problem with multiple measurement vectors (MMVs) in noisy circumstances. However, the DOAs are not usually exactly on the grid, which leads to the off-grid gap problem, as shown in Figure 2 [55]. In general, there are two kinds of off-grid methods to handle this problem: One is to use a fixed grid with joint estimation of the sparse signal and the grid offset [45,46]; the other one is to use a dynamic grid [18,47,48]. Compared with fixed grid methods, dynamic grid methods have the advantages of lower computational complexity and higher accuracy.

3. The Proposed DOA Estimation Algorithm

In this section, a robust off-grid DOA estimation algorithm is proposed to accurately estimate the azimuth angles of the target vehicles in vehicle localization. Based on the SBL formulation, all the observed and unknown variables are assumed to be stochastic and their joint prior probability distribution is specified. The joint distribution can be factored into individual prior or conditional distributions of parameters which can be estimated via a fast evidence maximization method [36]. Furthermore, an advanced grid refining strategy is designed to address the off-grid gap problem.

3.1. Sparse Bayesian Framework

Since the additive noise E is regarded to be complex Gaussian with noise variance σ 2 , p ( Y X ; σ 2 ) is also complex Gaussian. Thus, for each y ( t ) , x ( t ) pair, we have the likelihood of the array output as
p ( y ( t ) x ( t ) ; σ 2 ) = ( π σ 2 ) M exp ( 1 σ 2 y ( t ) A ˜ x ( t ) 2 2 ) ,
and hence
p ( Y X ; σ 2 ) = Π t = 1 L p ( y ( t ) x ( t ) ; σ 2 ) .
Following the commonly used SBL model [35], the source amplitudes X are assumed to be independent in terms of both snapshots and DOAs and follow a zero-mean complex Gaussian distribution. We have
p ( X ; α ) = Π t = 1 L CN ( x ( t ) 0 , Λ ) ,
where α = α 1 , α 2 , , α N T with α i being the inverse variance, and Λ = d i a g ( α ) . The hyperparameters α control the sparsity of the model. When α i = 0 , the corresponding source amplitude equals zero with probability 1. Thus, we estimate the hyperparameters α rather than the complex source amplitudes X , leading to significant reduction of parameters needed to be estimated.
Combining the likelihood of the array observations and the prior, the posterior distribution of X can be easily found by Bayes rule, and it is also a complex Gaussian with hyperparameters α and σ 2 :
p ( Y X , α , σ 2 ) = p ( Y X ; σ 2 ) p ( X ; α ) p ( Y ; α , σ 2 ) = Π t = 1 L p ( x ( t ) μ ( t ) , Σ ) ,
with posterior mean and covariance given by [35]
M = [ μ 1 , μ 2 , , μ L ] = Λ A ˜ H Σ y 1 Y ,
Σ = Λ Λ A ˜ H Σ y 1 A ˜ Λ ,
where Σ y 1 denotes the inverse of the array data covariance Σ y , and Σ y = E { y ( t ) y ( t ) H } = σ 2 I N + A Λ A H .

3.2. Hyperparameter Estimation

The hyperparameters α and σ 2 are estimated by an evidence maximization procedure [56], in which the unknown source amplitudes X are treated as nuisance parameters. The evidence is the product of the likelihood (9) and the prior (10) integrated over X , which is a function of α and σ 2 :
p ( Y ; α , σ 2 ) = p ( Y X ; σ 2 ) p ( X ; α ) d X = exp { tr ( Y H Σ y 1 Y ) } ( π N det Σ y ) L ,
and the marginal log-likelihood is
p ( Y ; α , σ 2 ) tr ( Y H Σ y 1 Y ) L ln ( det Σ y ) .
The hyperparameters α and σ 2 are obtained by maximizing Formula (15). The expectation maximization algorithm [35] is commonly used to maximize Formula (15), which is adopted by lots of methods [17,18,45,46,47]. Recently, a fast evidence maximization method has been proposed in an on-grid DOA estimation method named MSBL [36], in which the hyperparameters are updated alternatively. Since the fast evidence maximization method is proved to be faster than the expectation maximization algorithm and its local convergence is proved by [50], we adopted it for hyperparameters estimation in this paper. Readers can refer to [36] for more derivative details that are omitted here for brevity. We have
α m n e w = α m o l d L Y H Σ y 1 a m 2 / a m H S y 1 a m ,
( σ 2 ) n e w = 1 N K tr ( ( I N A M A M + ) S y ) ,
where S y = Y Y H / L , A M is composed of active columns from A that correspond to the η largest peaks in α , and A M + = ( A M H A M ) 1 A M H . If the number of sources K is available, η = K ; otherwise, η < M is recommended. As stated in [36], different choices of η only influence the convergence, and any choice 0 < η < M provides a better estimation than the method in [35].

3.3. Advanced Grid Refining Strategy

As mentioned in Section 2, the assumption that all DOAs are located on the predefined grid may not be met in reality. To handle this problem, there are two kinds of methods: the fixed grid methods and the dynamic grid methods. Since the dynamic grid methods have many advantages over the fixed grid methods, we provide a dynamic grid refining strategy consisting of a root process and a step process to eliminate the off-grid gap.

3.3.1. The Root Process

Considered as parameters, the grid points get refined by maximizing Formula (15). The maximization process is converted to a root seeking process [47]:
[ υ θ ˜ i , 1 , υ θ ˜ i 1 , , υ θ ˜ i ( M 2 ) ] M ( M 1 ) 2 ϕ ( i ) φ 2 ( i ) 2 φ 3 ( i ) ( M 1 ) φ M ( i ) = 0 ,
with
ϕ ( i ) t = 1 L ( | μ t i | 2 + γ i i ) ,
φ ( i ) L j i γ j i a j t = 1 L μ t i H y t i ,
where υ θ ˜ i e j 2 π d / λ sin ( θ ˜ i ) , a i , μ t i , and γ i j denote the ith column, the ith element, and the ( i , j ) th element of A ˜ , μ ( t ) , and Σ , respectively. In addition, y t i y ( t ) j i μ t j a j .
Note that Equation (18) has M 1 roots in the complex plane since its order is M 1 . In practice, the closest root to the unit circle is selected (which is denoted by z i * ) in the noisy environment. Then, the candidate point for grid refining is
θ ˜ i * n e w = arcsin ( λ 2 π d angle ( z i * ) ) ,
and θ ˜ i * n e w will be accepted if it falls into the set of [ θ ˜ i * 1 + θ ˜ i * 2 , θ ˜ i * + θ ˜ i * + 1 2 ] .

3.3.2. The Step Process

The root method [47] is susceptible to the location of grid points and DOAs, especially in a coarse grid condition. If a DOA is closed to the middle of two adjacent grid points, the total consumed time will increase and the mean square error (MSE) will become the largest [57]. To address this issue, a step process is designed in our method. θ ˜ i * n e w is compared with its value in the last iteration firstly, and then the grid point for grid refining is obtained as follows:
If θ ˜ i * n e w > θ ˜ i * ,
θ ˜ i * n e w = θ ˜ i * + θ s t e p ,
otherwise,
θ ˜ i * n e w = θ ˜ i * θ s t e p ,
where θ ˜ i * is the value of θ ˜ i * n e w in the last iteration, θ s t e p is a parameter with θ s t e p < 1 . The influence of different vaules of θ s t e p will be discussed in the next section. Then, we discuss the stop criterion for the step process. The extreme case is that a DOA is right in the middle of two adjacent grid points. The initial grid interval is assumed to be r. The number of iterations may be the round of r 2 θ s t e p for θ ˜ i * n e w to get very close to the DOA under the noiseless condition. Therefore, we set the round of r 2 θ s t e p as the maximum iteration time for the step process. When the step process reaches its stop criterion, the grid will continue to get updated via the root process until convergency.
As suggested in [47], we update active grid points rather than all grid points in each iteration. The first η maxima mean power of rows of X are regarded as the active grid points. The mean power estimation of the ith row of X can be calculated by P ( i ) = t = 1 L | μ t i | 2 , i = 1 , 2 , , N , where μ t i is estimated by Equation (12) in each iteration.
The stop criterion of the estimation is: α i + 1 α i 2 α i 2 < τ or the iteration reaches the predefined maximum time, where superscript i represents iteration time and τ is a settled tolerance.
The intuitive illustration of the grid refining procedures of RSBL and the proposed method is given in Figure 3. Figure 3a shows the root method of RSBL, while Figure 3b exhibits the advanced grid refining strategy of the proposed method. A true DOA is assumed to be in the middle of two adjacent initial grid points. For RSBL, an initial grid point is getting closer and closer to the true DOA via the root method after each iteration until the root method is convergent. Since the distance of the grid movement via the root method during one iteration is very small, it will lead to a large error in the condition of the coarse grid with a middle located DOA and also cost much time. The proposed method can effectively address these issues. The grid refining strategy contains two processes, i.e., the step process and the root process. The step process is first conducted to greatly reduce the distance between the grid point and the true DOA and to speed up the grid refining procedure because θ s t e p is much larger than the grid movement distance via the root method. Then, the root process is implemented to further reduce the estimation error. The influence of the location of grid points and DOAs of the proposed method is therefore much less than that of RSBL. Note that, although the grid point may move in the opposite direction of the true DOA in some iterations under the bad conditions, the general trend of movement is to get closer to the true DOA.

3.4. Operating Instruction

The operating instruction of the proposed method can be summarized as Table 1.
Both RSBL [47] and the proposed method are dynamic grid DOA estimation methods based on SBL. However, they are different from each other. First, they employed different hyperparameters estimation methods. An expectation maximization method is used in RSBL to estimate the hyperparameters, whereas a fast evidence maximization is adopted in the proposed method, which can achieve fast convergence. Second, they implemented different grid refining strategy. RSBL uses a polynomial root method to update the grid points, whereas the proposed method provides an advanced grid refining strategy consisting of a root process and a step process. The advanced grid refining strategy accelerates the iterative convergence and is less sensitive to the location of grid points and DOAs, which makes the proposed method more effective, robust, and accurate than RSBL. The improvement of the proposed method will be verified in the next section.

4. Simulation

In this section, the performance of the proposed vehicle localization system is mainly evaluated by the performance of the proposed DOA estimation method. We compare the proposed method with OGSBI [45], RSBL [47], l 1 -SVD [33], and MSBL [36] on accuracy and computational complexity. All the numerical simulations were carried out by MATLAB (R2016b, MathWorks, Natick, MA, USA) on a PC with an Intel i3-7350k CPU (Santa Clara, CA, USA) and 8 GB of RAM. As a representative of the on-grid model, the grid refinement strategy was not performed in l 1 -SVD. OGSBI was conducted without the SVD process.
In the following simulations, a ULA composed of M = 10 sensors with d = λ / 2 was used to receive K = 2 signals. The convergence parameters were set the same as the settings in [45]: ρ = 0.01 , c = d = 10 4 . σ 2 and all the elements of α were set to be 1. The tolerance τ = 10 3 and the maximum iteration time was 1000. It is assumed that two sources were uniformly chosen from the intervals of [ 20 , 10 ] and [ 20 , 30 ] , respectively. The estimation errors were investigated by the MSE defined as:
M S E = 1 K ξ k = 1 K i = 1 ξ ( θ ^ i , k θ k ) 2 ,
where ξ = 200 is the total number of Monte Carlo trials and θ ^ i , k is the estimated result of θ k in the ith Monte Carlo simulation.

4.1. Performance Analysis

We verify the improvement of the proposed method in terms of the MSE and the computational time. For this simulation, we set η = 2 and θ s t e p = 0.5 for the proposed method, and η = 2 for RSBL. Figure 4 shows the MSE of DOA estimation versus SNR, where the grid interval r = 4 and the number of snapshots L = 200 . SNR varies from 10 dB to 10 dB. As shown in Figure 4, while SNR increases, the MSEs of l 1 -SVD and MSBL almost keep constant and the MSEs of the other three methods decrease. It is clear that the proposed method has the lowest MSE among the five methods. The main reason is that l 1 -SVD and MSBL are on-grid methods, whose performance in a coarse grid condition is mainly limited by the off-grid error. Therefore, the increase of SNR cannot effectively reduce this error. The other three methods are off-grid methods, which can effectively handle the off-grid gap. RSBL and the proposed method perform better than OGSBI. This is because the linear approximation used by OGSBI brings large model error in the coarse grid condition, while grid refining methods adopted by RSBL and the proposed method effectively alleviate the off-grid gap. Compared with RSBL, the proposed method can achieve higher accuracy. The performance of RSBL is seriously affected by the location of grid points and DOAs. If a DOA is right in the middle of two adjacent grid points, the performance of RSBL becomes the worst. On the contrary, if the DOA is close to a grid point, the MSE decreases nearly to a noise level. However, for the proposed method, the advanced grid refining strategy makes the refining grid points much closer to the true DOAs, especially in the case of coarse grid and middle located DOAs. Therefore, the influence of the location of grid points and DOAs of the proposed method is much less than that of RSBL.
Figure 5 shows the MSE of DOA estimation versus grid interval, where S N R = 10 dB and the number of snapshots L = 200 . The grid interval varies from 2 to 6 . As shown in Figure 5, off-grid methods outperform on-grid methods in general, and the proposed method outperforms the state-of-art methods. When the grid interval becomes larger, the MSEs of l 1 -SVD, MSBL, and OGSBI increase, while RSBL and the proposed method still maintains high accuracy. This is because grid refining methods, in which the grid points are considered as adjustable parameters, adopted by RSBL and the proposed method, can effectively handle the modeling error caused by a coarse grid. The proposed method effectively eliminates the influence of the location of grid points and DOAs, which greatly affects RSBL.
Figure 6 shows the MSE of DOA estimation versus the number of snapshots, where the grid interval is set as 4 and S N R = 10 dB. The number of snapshots varies from 50 to 300. As shown in Figure 6, the proposed method has the lowest MSE among the five methods. With the increased number of the snapshots, the MSEs of l 1 -SVD, MSBL, and OGSBI almost remain unchanged, while that of RSBL and the proposed method decrease. It is because the performances of l 1 -SVD, MSBL, and OGSBI in coarse grid condition are mainly limited by the model error, and the increased number of snapshots cannot effectively reduce this error. RSBL and the proposed method can effectively handle this model error via grid refining methods. Furthermore, the proposed method efficiently eliminates the influence of the location of grid points and DOAs while RSBL does not.
Figure 7 shows the total CPU time of the five methods versus grid interval, where S N R = 10 dB and the number of snapshots L = 200 . Simulation was carried out by MATLAB on a PC with an Intel i3-7350k CPU and 8 GB of RAM. The grid interval varies from 2 to 6 . As shown in Figure 7, MSBL is the fastest of all, and the proposed method is the fastest among the off-grid methods when the grid is coarse. In general, the off-grid methods are able to address the off-grid gap at the cost of much more computational workload. When the grid interval is greater than 3 , RSBL and the proposed method are still more efficient than l 1 -SVD, and the proposed method is faster than RSBL. For the proposed method, the advanced grid refining strategy greatly speeds up the procedure of grid refining, which makes it more efficient than RSBL.

4.2. Influence of Parameter Settings and Conditions

First, we investigate the influence of the parameter θ s t e p on the DOA estimation performance. SNR is set as 10 dB, the number of snapshots L = 200 and η = 2 . Figure 8 shows the MSE of DOA estimation versus θ s t e p with different grid intervals, and Figure 9 shows the CPU time of DOA estimation versus θ s t e p with different grid intervals, where the grid intervals are set as r = 3 , 4 , 5 , and 6 , respectively. Simulations were carried out by MATLAB on a PC with an Intel i3-7350k CPU and 8 GB of RAM. In short, different values of θ s t e p do not affect the performance much. A relatively large value of θ s t e p brings light efficiency in coarse grid conditions. In general, θ s t e p = 0.3 or 0.5 is recommended.
Then, we investigate the influence of the number of active grid points on the DOA estimation performance. SNR is set as 10 dB, the number of snapshots L = 200 , and θ s t e p = 0.5 . Figure 10 shows the MSE of DOA estimation versus the number of active grid points with different grid intervals, and Figure 11 shows the CPU time of DOA estimation versus the number of active grid points with different grid intervals, where the grid intervals are set as r = 3 , 4 , 5 , and 6 , respectively. Simulations were carried out by MATLAB on a PC with an Intel i3-7350k CPU and 8 GB of RAM. It is shown that the number of active grid points barely affects the performance in coarse grid conditions. However, if a relatively large value is chosen, there is a slight increase in computational cost.
Finally, we illustrate the influence of different conditions on the proposed method. Figure 12 shows the MSE of DOA estimation versus grid interval with S N R = 0 dB and L = 30 . Compared with Figure 5, which shows the MSE versus grid interval with S N R = 10 dB and L = 200 , the improvement of the proposed method in a bad condition becomes less significant compared with a good condition with higher SNR and a larger number of snapshots. Since the advanced grid refining strategy of the proposed method is an adaptive process, it thus will be affected by the SNR as well as the number of snapshots. Despite that, the proposed method still performs better than other methods even in bad conditions and the improvement of the proposed method will be strengthened in good conditions with a large number of snapshots and high SNR.

4.3. Compared with GRDOA

Both GRDOA [48] and the proposed method are off-grid DOA estimation methods based on SBL. The discrete grid points are regarded as the dynamic parameters in both methods to reduce the model error, and the performance degradation in RSBL [47] caused by middle located DOAs is effectively alleviated in both methods. However, the essential differences between GRDOA and the proposed method are:
(1)
Different grid structures: the overall grid structure of GRDOA varies during the algorithm execution. Some new grid points are generated by the fission process, whereas some are discarded according to a certain criterion in the initial estimation of GRDOA. However, the overall grid structure of the proposed method stays unchanged.
(2)
Different hyperparameters estimation methods: an expectation maximization method is used in GRDOA to estimate the hyperparameters, whereas a fast evidence maximization method is used in the proposed method.
(3)
Different grid refining strategies: the polynomial root strategy [47] is adopted in the grid update process of GRDOA, whereas an advanced grid refining strategy consisting of a root process and a step process is provided in the proposed method.
Simulations are conducted to compare the proposed method with GRDOA and RSBL. Assume that two sources were uniformly chosen from the intervals of [ 20 , 10 ] and [ 20 , 30 ] , respectively. ρ = 0.01 , c = d = 10 4 . σ 2 and all the elements of α were set to be 1. The tolerance τ = 10 3 and the maximum iteration time was 1000. The number of active grid points η = 2 . θ s t e p = 0.5 for the proposed method. Simulations are based on 200 Monte Carlo trials.
Figure 13 shows the MSE of DOA estimation versus the number of snapshots, where the grid interval is set as 4 . As shown in Figure 13, both GRDOA and the proposed method achieve higher accuracy than RSBL, since they can effectively reduce the performance degradation caused by middle located DOAs. The proposed method performs better than GRDOA. There are three main reasons. First, the fission process of GRDOA may introduce additional errors. When a new grid point is generated by the fission process, a new element corresponding to that grid point should be added to some parameters. Because the value of this new element is calculated based on some assumptions and principles, it will lead to a gap between the calculated value and the real value. Second, the performance degradation problem caused by middle located DOAs still exists in GRDOA because the root method [39] is used in the grid update process of GRDOA. In other words, even though the fission process can address this issue to some extent, the fundamental problem has not been fully solved. For the proposed method, an advanced grid refining strategy is used for grid refining which fully solves the middle located DOA problem. Finally, since the final grid number of GRDOA is quite small, additional model error may be introduced to some extent. This is because, for off-grid methods, a small number of grid points brings high efficiency, but, at the same time, it causes a large model error in general. However, the grid number of the proposed method is much larger than that of GRDOA with the same initial grid interval, so the model error of the proposed method is much less than that of GRDOA.
Figure 14 shows the total CPU time versus grid interval, where S N R = 10 dB and the number of snapshots L = 200 . Simulation was carried out by MATLAB on a PC with an Intel i3-7350k CPU and 8 GB of RAM. GRDOA is faster than the proposed method when r < 8 . As the grid number of GRDOA is much smaller than that of the proposed method under a relatively fine grid condition, GRDOA achieves more efficiency than the proposed method. However, this superiority of GRDOA decreases as the grid interval grows. The proposed method is faster than GRDOA when r > 8 . Furthermore, the proposed method shows more efficiency than RSBL when the grid interval is greater than 3 , as the advanced grid refining strategy greatly speeds up the procedure of grid refining.
In general, the proposed method is more suitable than GRDOA for the vehicle localization system, as the accuracy of the vehicle localization system is highly dependent on the accuracy of the DOA estimation and the proposed DOA method shows a great improvement in the accuracy.

4.4. Vehicle Localization

The vehicle localization performance of different DOA estimation methods is tested based on the proposed localization model, where the grid interval is set as 4 . For this simulation, we set η = 2 and θ s t e p = 0.5 for the proposed method, and η = 2 for RSBL. As shown in Figure 15, there are two vehicles with position coordinates V 1 (120 m, 300 m) and V 2 (320 m, 335 m), respectively. Two stations of a bistatic passive radar are located at A ( 0 , 0 ) and B(500 m, 0), respectively. Theoretically, the DOAs of these two vehicles with respect to A and B are θ 1 A = 21.8 , θ 1 B = 51.71 and θ 2 A = 43.69 , θ 2 B = 28.25 . Assume that θ 1 A * , θ 1 B * , θ 2 A * , and θ 2 B * are the estimations of θ 1 A , θ 1 B , θ 2 A , and θ 2 B , respectively, obtained by the DOA estimation algorithm. The location of V 1 can be calculated by putting θ 1 A * and θ 1 B * into (3) and (4). Similarly, the location of V 2 can be calculated by putting θ 2 A * and θ 2 B * into (3) and (4). The vehicle localization results and errors versus different DOA estimation methods with different SNR are given in Table 2, Table 3 and Table 4. By comparing these results, it can be obviously found that the localization error of the proposed method is the smallest among the four methods, which is consistent with the performance of DOA estimation.

5. Conclusions

In this paper, we present a passive bistatic radar-based vehicle localization system as well as a robust SBL-based off-grid DOA estimation method for vehicle localization. With this proposed DOA estimation method, a fast evidence maximization method is used to iteratively update the source powers and the noise variance. In order to address the off-grid gap issue, an advanced grid refining strategy is designed to iteratively update the discrete grid points regarded as the dynamic parameters. The position of target vehicles can be easily obtained via the DOA estimation results. Simulation results show that the proposed method exhibits superiority and effectiveness in DOA estimation as well as vehicle localization. However, the experiment was only performed by simulations since we don’t have the real experimental data yet. In the future, the verification of the proposed method with the real experimental data set is needed. In addition, hybrid vehicle localization methods, combining the proposed method with TDOA based solutions, may be a potential solution to further achieve excellent robustness and accuracy.

Author Contributions

Y.L. and H.G. conceived and designed the experiments; Y.L. and S.Z. performed the experiments; Y.L., L.Y., and F.R. analyzed the data; S.Z., L.Y., and F.R. contributed reagents/materials/analysis tools; Y.L. and H.G. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This paper is supported by the National Natural Science Foundation of China (No.61671333), the Natural Science Foundation of Hubei Province (2014CFA093), and the Fundamental Research Funds for the Central Universities (2042019K50264, 2042019gf0013).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, R. A new communication-efficient privacy-preserving pange query scheme in fog-enhanced IoT. IEEE Internet Things J. 2018, 6, 2497–2505. [Google Scholar] [CrossRef]
  2. Lu, R.; Heung, K.; Lashkari, A.; Ghorbani, A. A lightweight privacy-preserving data aggregation scheme for fog computing-enhanced IoT. IEEE Access 2017, 5, 3302–3312. [Google Scholar] [CrossRef]
  3. Rathee, G.; Sharma, A.; Iqbal, R.; Aloqaily, M.; Jaglan, N.; Kumar, R. A Blockchain Framework for Securing Connected and Autonomous Vehicles. Sensors 2019, 19, 3165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Huang, C.; Lu, R.; Lin, X.; Shen, X. Secure automated valet parking: A privacy-preserving reservation scheme for autonomous vehicles. IEEE Trans. Veh. Technol. 2018, 67, 11169–11180. [Google Scholar] [CrossRef]
  5. Kuutti, S.; Fallah, S.; Katsaros, K.; Dianati, M.; Mccullough, F.; Mouzakitis, A. A survey of the state-of-the-art localization techniques and their potentials for autonomous vehicle applications. IEEE Internet Things J. 2018, 5, 829–846. [Google Scholar] [CrossRef]
  6. Jo, K.; Kim, J.; Kim, D.; Jang, C.; Sunwoo, M. Development of autonomous car—Part I: Distributed system architecture and development process. IEEE Trans. Ind. Electron. 2014, 61, 7131–7140. [Google Scholar] [CrossRef]
  7. Vivet, D.; Gérossier, F.; Checchin, P.; Trassoudaine, L.; Chapuis, R. Mobile ground-based radar sensor for localization and mapping: Anevaluation of two approaches. Int. J. Adv. Robot. Syst. 2013, 10, 307–318. [Google Scholar] [CrossRef]
  8. Hsu, C.M.; Shiu, C.W. 3D LiDAR-Based Precision Vehicle Localization with Movable Region Constraints. Sensors 2019, 19, 942. [Google Scholar] [CrossRef] [Green Version]
  9. Ren, R.; Fu, H.; Wu, M. Large-Scale Outdoor SLAM Based on 2D Lidar. Electronics 2019, 8, 613. [Google Scholar] [CrossRef] [Green Version]
  10. Cai, H.; Hu, Z.; Huang, G.; Zhu, D.; Su, X. Integration of GPS, Monocular Vision, and High Definition (HD) Map for Accurate Vehicle Localization. Sensors 2018, 18, 3270. [Google Scholar] [CrossRef] [Green Version]
  11. Ferreira, B.; Matos, A.; Cruz, N. Optimal positioning of autonomous marine vehicles for underwater acoustic source localization using toa measurements. In Proceedings of the 2013 IEEE International Underwater Technology Symposium (UT), Tokyo, Japan, 5–8 March 2013; pp. 1–7. [Google Scholar]
  12. Wahab, A.A.; Khattab, A.; Fahmy, Y.A. Two-way TOA with limited dead reckoning for GPS-free vehicle localization using single RSU. In Proceedings of the 2013 13th International Conference on ITS Telecommunications (ITST), Tampere, Finland, 5–7 November 2013; pp. 244–249. [Google Scholar]
  13. Jin, B.; Xu, X.; Zhu, Y.; Zhang, T.; Fei, Q. Single-Source Aided Semi-Autonomous Passive Location for Correcting the Position of an Underwater Vehicle. IEEE Sens. J. 2019, 19, 3267–3275. [Google Scholar] [CrossRef]
  14. Huang, B.; Xie, L.; Yang, Z. TDOA-based source localization with distance-dependent noises. IEEE Trans. Wirel. Commun. 2015, 14, 468–480. [Google Scholar] [CrossRef]
  15. Sallouha, H.; Azari, M.M.; Chiumento, A.; Pollin, S. Aerial anchors positioning for reliable rss-based outdoor localization in urban environments. IEEE Wirel. Commun. Lett. 2017, 7, 376–379. [Google Scholar] [CrossRef]
  16. Sun, W.; Xue, M.; Yu, H.; Tang, H.; Lin, A. Augmentation of Fingerprints for Indoor WiFi Localization Based on Gaussian Process Regression. IEEE Trans. Veh. Technol. 2018, 67, 10896–10905. [Google Scholar] [CrossRef]
  17. Wang, X.; Huang, M.; Shen, C.; Meng, D. Robust vehicle localization exploiting two based stations cooperation: A MIMO radar perspective. IEEE Access 2018, 6, 48747–48755. [Google Scholar] [CrossRef]
  18. Wang, H.; Wan, L.; Dong, M.; Ota, K.; Wang, X. Assistant vehicle localization based on three collaborative base stations via SBL-based robust DOA estimation. IEEE Internet Things J. 2019, 6, 5766–5777. [Google Scholar] [CrossRef]
  19. Zhou, B.; Yao, X.; Yang, L.; Yang, S.; Wu, S.; Kim, Y.; Ai, L. Accurate Rigid Body Localization Using DoA Measurements from a Single Base Station. Electronics 2019, 8, 622. [Google Scholar] [CrossRef] [Green Version]
  20. Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, AP-34, 276–280. [Google Scholar] [CrossRef] [Green Version]
  21. Zhou, C.; Shi, Z.; Gu, Y.; Shen, X. DECOM: DOA Estimation with Combined MUSIC for Coprime Array. In Proceedings of the 2013 International Conference on Wireless Communications and Signal Processing (WCSP), Hangzhou, China, 24–26 October 2013; pp. 1–5. [Google Scholar]
  22. Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A.; Qu, F. Source estimation using coprime array: A sparse reconstruction perspective. IEEE Sens. J. 2016, 17, 755–765. [Google Scholar] [CrossRef]
  23. Gu, Y.; Goodman, N.A. Information-theoretic compressive sensing kernel optimization and Bayesian Cramér–Rao bound for time delay estimation. IEEE Trans. Signal Process. 2017, 65, 4525–4537. [Google Scholar]
  24. Zhou, C.; Zhou, J. Direction-of-arrival estimation with coarray ESPRIT for coprime array. Sensors 2017, 17, 1779. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zhou, C.; Gu, Y.; Zhang, Y.D.; Shi, Z.; Jin, T.; Wu, X. Compressive Sensing based Coprime Array Direction-of-Arrival Estimation. IET Commun. 2017, 11, 1719–1724. [Google Scholar] [CrossRef]
  26. Gong, P.C.; Wang, W.Q.; Li, F.; Thus, H.C. Sparsity-Aware Transmit Beamspace Design for FDA-MIMO Radar. Signal Process. 2018, 144, 99–103. [Google Scholar] [CrossRef]
  27. Zhou, C.; Gu, Y.; Fan, X.; Shi, Z.; Mao, G.; Zhang, Y.D. Direction-of-arrival estimation for coprime array via virtual array interpolation. IEEE Trans. Signal Process. 2018, 66, 5956–5971. [Google Scholar]
  28. Zhou, C.; Gu, Y.; He, S.; Shi, Z. A Robust and Efficient Algorithm for Coprime Array Adaptive Beamforming. IEEE Trans. Veh. Technol. 2018, 67, 1099–1112. [Google Scholar] [CrossRef]
  29. Gu, Y.; Zhang, Y.D. Information-theoretic pilot design for downlink channel estimation in FDD massive MIMO systems. IEEE Trans. Signal Process. 2019, 67, 2334–2346. [Google Scholar]
  30. Chen, P.; Cao, Z.; Chen, Z.; Wang, X. Off-Grid DOA Estimation Using Sparse Bayesian Learning in MIMO Radar With Unknown Mutual Coupling. IEEE Trans. Signal Process. 2019, 67, 208–220. [Google Scholar]
  31. Wang, X.; Huang, M.; Wu, X.; Bi, G. Direction of arrival estimation for MIMO radar via unitary nuclear norm minimization. Sensors 2017, 17, 939. [Google Scholar] [CrossRef] [Green Version]
  32. Zhang, X.; Zhou, M.; Li, J. A PARALIND decomposition-based coherent two-dimensional direction of arrival estimation algorithm for acoustic vector-sensor arrays. Sensors 2013, 13, 5302–5316. [Google Scholar] [CrossRef] [Green Version]
  33. Malioutov, D.; Cetin, M.; Willsky, A. A sparse signal reconstruction perspective for source localization with sensor arrays. IEEE Trans. Signal Process. 2005, 53, 3010–3022. [Google Scholar] [CrossRef] [Green Version]
  34. Stoica, P.; Babu, P.; Li, J. SPICE: A sparse covariance-based estimation method for array processing. IEEE Trans. Signal Process. 2011, 59, 629–638. [Google Scholar] [CrossRef]
  35. Wipf, D.P.; Rao, B.D. An empirical Bayesian strategy for solving the simultaneous sparse approximation problem. IEEE Trans. Signal Process. 2007, 55, 3704–3716. [Google Scholar] [CrossRef]
  36. Gerstoft, P.; Mecklenbräuker, C.F.; Xenaki, A.; Nannuru, S. Multisnapshot sparse Bayesian learning for DOA. IEEE Signal Process. Lett. 2016, 23, 1469–1473. [Google Scholar] [CrossRef] [Green Version]
  37. Yang, Z.; Tang, J.; Eldar, Y.C.; Xie, L. On the sample complexity of multichannel frequency estimation via convex optimization. IEEE Trans. Inf. Theory 2018, 65, 2302–2315. [Google Scholar] [CrossRef]
  38. Wu, X.; Zhu, W.P.; Yan, J. A fast gridless covariance matrix reconstruction method for one-and two-dimensional direction-of-arrival estimation. IEEE Sens. J. 2017, 17, 4916–4927. [Google Scholar] [CrossRef]
  39. Wu, X.; Zhu, W.P.; Yan, J. A Toeplitz covariance matrix reconstruction approach for direction-of-arrival estimation. IEEE Trans. Veh. Technol. 2017, 66, 8223–8237. [Google Scholar] [CrossRef]
  40. Wu, X.; Zhu, W.P.; Yan, J.; Zhang, Z. A Spatial Filtering Based Gridless DOA Estimation Method for Coherent Sources. IEEE Access 2018, 6, 56402–56410. [Google Scholar] [CrossRef]
  41. Cui, Y.; Wang, J.; Qi, J.; Zhang, Z.; Zhu, J. Underdetermined DOA Estimation of Wideband LFM Signals Based on Gridless Sparse Reconstruction in the FRF Domain. Sensors 2019, 19, 2383. [Google Scholar] [CrossRef] [Green Version]
  42. Han, M.; Dou, W. Atomic Norm-Based DOA Estimation with Dual-Polarized Radar. Electronics 2019, 8, 1056. [Google Scholar] [CrossRef] [Green Version]
  43. Zhou, C.; Gu, Y.; Shi, Z.; Zhang, Y.D. Off-Grid Direction-of-Arrival Estimation Using Coprime Array Interpolation. IEEE Signal Process. Lett. 2018, 25, 1710–1714. [Google Scholar] [CrossRef]
  44. Han, K.; Yang, P.; Nehorai, A. Calibrating nested sensor arrays with model errors. IEEE Trans. Antennas Propag. 2015, 63, 4739–4748. [Google Scholar] [CrossRef]
  45. Yang, Z.; Xie, L.; Zhang, C. Off-grid direction of arrival estimation using sparse Bayesian inference. IEEE Trans. Signal Process. 2013, 61, 38–43. [Google Scholar] [CrossRef] [Green Version]
  46. Wu, X.; Zhu, W.P.; Yan, J. Direction of arrival estimation for off-grid signals based on sparse Bayesian learning. IEEE Sens. J. 2016, 16, 2004–2016. [Google Scholar] [CrossRef]
  47. Dai, J.; Bao, X.; Xu, W.; Chang, C. Root sparse Bayesian learning for off-grid DOA estimation. IEEE Signal Process. Lett. 2017, 24, 46–50. [Google Scholar] [CrossRef] [Green Version]
  48. Ling, Y.; Gao, H.; Ru, G.; Chen, H.; Li, B.; Cao, T. Grid Reconfiguration Method for Off-Grid DOA Estimation. Electronics 2019, 8, 1209. [Google Scholar] [CrossRef] [Green Version]
  49. Han, G.; Wan, L.; Shu, L.; Feng, N. Two novel DOA estimation approaches for real-time assistant calibration systems in future vehicle industrial. IEEE Syst. J. 2017, 11, 1361–1372. [Google Scholar] [CrossRef]
  50. Zhang, X.; Huo, K.; Liu, Y.; Li, X. Direction of Arrival Estimation via Joint Sparse Bayesian Learning for Bi-static Passive Radar. IEEE Access 2019, 7, 72979–72993. [Google Scholar] [CrossRef]
  51. Colone, F.; Falcone, P.; Bongioanni, C.; Lombardo, P. WiFi-based passive bistatic RADAR: Data processing schemes and experimental results. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 1061–1079. [Google Scholar] [CrossRef]
  52. Raja, R.A.; Noor, A.A.; Nur, A.R.; Asem, A.S.; Fazirulhisyam, H. Analysis on target detection and classification in lte based passive forward scattering radar. Sensors 2016, 16, 1607. [Google Scholar] [CrossRef] [Green Version]
  53. Falcone, P.; Colone, F.; Macera, A.; Lombardo, P. Two-dimensional location of moving targets within local areas using WiFi-based multistatic passive radar. IET Radar Sonar Navig. 2014, 8, 123–131. [Google Scholar] [CrossRef]
  54. Chetty, K.; Smith, G.E.; Woodbridge, K. Through-the-wall sensing of personnel using passive bistatic wifi radar at standoff distances. IEEE Trans. Geosci. Remote Sens. 2011, 50, 1218–1226. [Google Scholar] [CrossRef]
  55. Wang, H.; Wang, X.; Wan, L.; Huang, M. Robust Sparse Bayesian Learning for Off-Grid DOA Estimation With Non-Uniform Noise. IEEE Access 2018, 6, 64688–64697. [Google Scholar] [CrossRef]
  56. MacKay, D. Bayesian interpolation. Neural Comput. 1992, 4, 415–447. [Google Scholar] [CrossRef]
  57. Wang, Q.; Zhao, Z.; Chen, Z.; Nie, Z. Grid evolution method for DOA estimation. IEEE Trans. Signal Process. 2018, 66, 2374–2383. [Google Scholar] [CrossRef]
Figure 1. The DOA estimation-based vehicle localization model with passive bistatic radar.
Figure 1. The DOA estimation-based vehicle localization model with passive bistatic radar.
Sensors 20 00302 g001
Figure 2. The off-grid DOA estimation model.
Figure 2. The off-grid DOA estimation model.
Sensors 20 00302 g002
Figure 3. The grid refining procedures of the root method and the advanced grid refining strategy.
Figure 3. The grid refining procedures of the root method and the advanced grid refining strategy.
Sensors 20 00302 g003
Figure 4. Mean square error (MSE) of DOA estimation versus signal-to-noise ratio (SNR) with r = 4 and L = 200 .
Figure 4. Mean square error (MSE) of DOA estimation versus signal-to-noise ratio (SNR) with r = 4 and L = 200 .
Sensors 20 00302 g004
Figure 5. MSE of DOA estimation versus grid interval with S N R = 10 dB and L = 200 .
Figure 5. MSE of DOA estimation versus grid interval with S N R = 10 dB and L = 200 .
Sensors 20 00302 g005
Figure 6. MSE of DOA estimation versus the number of snapshots with S N R = 10 dB.
Figure 6. MSE of DOA estimation versus the number of snapshots with S N R = 10 dB.
Sensors 20 00302 g006
Figure 7. Computational time versus grid interval with S N R = 10 dB and L = 200 .
Figure 7. Computational time versus grid interval with S N R = 10 dB and L = 200 .
Sensors 20 00302 g007
Figure 8. MSE of DOA estimation versus θ s t e p with S N R = 10 dB and L = 200 .
Figure 8. MSE of DOA estimation versus θ s t e p with S N R = 10 dB and L = 200 .
Sensors 20 00302 g008
Figure 9. Computational time of DOA estimation versus θ s t e p with S N R = 10 dB and L = 200 .
Figure 9. Computational time of DOA estimation versus θ s t e p with S N R = 10 dB and L = 200 .
Sensors 20 00302 g009
Figure 10. MSE of DOA estimation versus the number of active grid points with S N R = 10 dB and L = 200 .
Figure 10. MSE of DOA estimation versus the number of active grid points with S N R = 10 dB and L = 200 .
Sensors 20 00302 g010
Figure 11. Computational time of DOA estimation versus the number of active grid points with S N R = 10 dB and L = 200 .
Figure 11. Computational time of DOA estimation versus the number of active grid points with S N R = 10 dB and L = 200 .
Sensors 20 00302 g011
Figure 12. MSE of DOA estimation versus grid interval with S N R = 0 dB and L = 30 .
Figure 12. MSE of DOA estimation versus grid interval with S N R = 0 dB and L = 30 .
Sensors 20 00302 g012
Figure 13. MSE of DOA estimation versus the number of snapshots with the grid interval r = 4 .
Figure 13. MSE of DOA estimation versus the number of snapshots with the grid interval r = 4 .
Sensors 20 00302 g013
Figure 14. Computational time versus grid interval with S N R = 10 dB and L = 200 .
Figure 14. Computational time versus grid interval with S N R = 10 dB and L = 200 .
Sensors 20 00302 g014
Figure 15. Localization diagram for two vehicles.
Figure 15. Localization diagram for two vehicles.
Sensors 20 00302 g015
Table 1. Operating instruction of the proposed method.
Table 1. Operating instruction of the proposed method.
The Proposed Method.
(1) Input: Y ;
(2) Initialize: α , σ 2 , and θ s t e p ;
(3) while not converge do
  (a) Calculate the posterior moments μ and Σ using Equations (12) and (13);
  (b) Update the hyperparameters α and σ 2 using Equations (16) and (17);
  (c) Implement the advanced grid refining strategy to update θ ˜ :
    Conduct the root process to calculate the candidate point using Equations (19)–(21);
    if not reach the maximum iteration time of the step process
    Conduct the step process using Equation (22) or Equation (23) for grid refining;
    else
    Directly use the candidate point for grid refining;
    end if
end while
(4) Output: μ and the final grid;
(5) Achieve off-grid DOA estimation through 1D spectrum search on the final grid.
Table 2. Vehicle localization results and errors with S N R = 10 dB.
Table 2. Vehicle localization results and errors with S N R = 10 dB.
MethodLocation V 1 Error V 1 Location V 2 Error V 2
MSBL(126.59, 313.33)14.87(304.65, 338.35)15.71
OGSBI(126.57, 313.02)14.58(304.57, 338.25)15.77
RSBL(126.38, 309.46)11.41(306.34, 338.06)14.00
Proposed(121.81, 298.77)2.19(317.72, 337.13)3.12
Multisnapshot sparse Bayesian learning (MSBL) [36], Off-grid sparse Bayesian inference (OGSBI) [45], Root sparse Bayesian learning (RSBL) [47].
Table 3. Vehicle localization results and errors with S N R = 0 dB .
Table 3. Vehicle localization results and errors with S N R = 0 dB .
MethodLocation  V 1 Error  V 1 Location  V 2 Error  V 2
MSBL(126.59, 313.33)14.87(304.65, 338.35)15.71
OGSBI(125.77, 311.60)12.96(306.85, 338.51)13.61
RSBL(122.89, 305.56)6.27(314.15, 335.66)5.89
Proposed(119.87, 299.16)0.85(319.13, 333.77)1.51
Table 4. Vehicle localization results and errors with S N R = 10 dB .
Table 4. Vehicle localization results and errors with S N R = 10 dB .
MethodLocation  V 1 Error  V 1 Location  V 2 Error  V 2
MSBL(126.59, 313.33)14.87(304.65, 338.35)15.71
OGSBI(122.02, 304.16)4.62(315.04, 337.36)5.49
RSBL(120.81, 301.29)1.52(318.44, 336.06)1.89
Proposed(120.03, 299.84)0.16(320.12, 334.84)0.20

Share and Cite

MDPI and ACS Style

Ling, Y.; Gao, H.; Zhou, S.; Yang, L.; Ren, F. Robust Sparse Bayesian Learning-Based Off-Grid DOA Estimation Method for Vehicle Localization. Sensors 2020, 20, 302. https://doi.org/10.3390/s20010302

AMA Style

Ling Y, Gao H, Zhou S, Yang L, Ren F. Robust Sparse Bayesian Learning-Based Off-Grid DOA Estimation Method for Vehicle Localization. Sensors. 2020; 20(1):302. https://doi.org/10.3390/s20010302

Chicago/Turabian Style

Ling, Yun, Huotao Gao, Sang Zhou, Lijuan Yang, and Fangyu Ren. 2020. "Robust Sparse Bayesian Learning-Based Off-Grid DOA Estimation Method for Vehicle Localization" Sensors 20, no. 1: 302. https://doi.org/10.3390/s20010302

APA Style

Ling, Y., Gao, H., Zhou, S., Yang, L., & Ren, F. (2020). Robust Sparse Bayesian Learning-Based Off-Grid DOA Estimation Method for Vehicle Localization. Sensors, 20(1), 302. https://doi.org/10.3390/s20010302

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop