Next Article in Journal
Autonomous Ride-Sharing Service Using Graph Embedding and Dial-a-Ride Problem: Application to the Last-Mile Transit in Lyon City
Previous Article in Journal
Investigating Behavior of Slider–Crank Mechanisms with Bearing Failures Using Vibration Analysis Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Sparse Bayesian Learning-Based Direction of Arrival Estimation Method with Array Position Errors

1
School of Information and Control Engineering, Qingdao University of Technology, Qingdao 266520, China
2
School of Marine Science and Technology, Northwest Polytechnical University, Xi’an 710129, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(4), 545; https://doi.org/10.3390/math12040545
Submission received: 5 January 2024 / Revised: 30 January 2024 / Accepted: 4 February 2024 / Published: 9 February 2024
(This article belongs to the Section Engineering Mathematics)

Abstract

:
In practical applications, the hydrophone array has element position errors, which seriously degrade the performance of the direction of arrival estimation. We propose a direction of arrival (DOA) estimation method based on sparse Bayesian learning using existing array position errors to solve this problem. The array position error and angle grid error parameters are introduced, and the prior distribution of these two errors is determined. The joint probability density distribution function is established by means of a sparse Bayesian learning model. At the same time, the unknown parameters are optimized and iterated using the expectation maximum algorithm and the corresponding parameters are solved to obtain the spatial spectrum. The results of the simulation and the lake experiments show that the proposed method effectively overcomes the problem of array element position errors and has strong robustness. It shows a good performance in terms of its estimation accuracy, meaning that the resolution ability can be greatly improved in the case of a low signal-to-noise ratio or small number of snapshots.

1. Introduction

Direction of arrival (DOA) estimation is still an important research direction in array signal processing. Its main principle is analyzing the characteristic information of signals (such as the angle of arrival, frequency, etc.) by using various algorithms in the data received by the sensor array. It has been widely used in the fields of radar, sonar, microphones, wireless communication, and other fields [1,2,3,4,5]. Traditional subspace class methods such as the multiple signal classification (MUSIC) method [6,7], a form of signal parameter estimation based on a rotation invariant technique [8], have attracted more attention from researchers. These methods are simple and can achieve real-time processing, but these methods also have some shortcomings, e.g., their estimation performance is severely degraded in cases where there is a low signal-to-noise ratio (SNR) or small snapshot number [9].
The array aperture is one of the key factors affecting DOA estimation. Compared with a uniform linear array with the same number of elements, the coprime array has a larger array aperture, which will provide a better estimation performance. Pal et al. proposed the spatial smoothing MUSIC (SS-MUSIC) method based on the coprime array, which could detect more signal sources than the elements of the coprime array and had a better estimation performance [10]. However, this method used the coprime array differential operation to form the augmented virtual array; the virtual elements of the discontinuous part were ignored in the smoothing process, and the information of the virtual array was not fully utilized. Zhou and Weng decomposed the coprime array into two uniform linear arrays for DOA estimation [11,12], and a unique estimation of the signal source was obtained from the results of phase ambiguity. This method reduced the number of identifiable targets by at least 1/2 when compared with the traditional uniform linear array [13]. Li proposed an extended aperture rooting MUSIC method in [14], but this is still a derivation and application of the traditional subspace class method, so its estimation accuracy is still reduced by factors such as the SNR or snapshot number.
In recent years, researchers have proposed many DOA estimation methods [15,16,17,18,19,20,21,22,23] with sparse spatial characteristics. These sparse signal representation methods not only discretize the angle range of interest into a spatial angle grid, but also assume that the real signal angle falls on a predefined grid, such as the L1-norm singular value decomposition (L1-SVD) method [15]. This method can fully exploit the advantage provided by the high degree of freedom (DOF) of the coprime array, so the number of detected targets can be significantly increased in comparison to traditional methods. In practice, the real incident angle may not be accurately located on the sampling grid no matter how dense the grid is, which results in a decrease in the accuracy of incident angle estimation. To solve this problem, an off-grid sparse Bayesian inference (OGSBI) method is proposed in [16], which introduces the angle grid error parameter and determines the prior distribution, and the joint probability density distribution function was established with the help of a sparse Bayesian learning model [17,18,19]. At the same time, the maximum expectation algorithm was used to optimize the iteration expression of each unknown parameter and solve the corresponding parameters; as a result, the DOA estimation performance was improved for the case of an off-grid incident signal. Sparse Bayesian learning has good flexibility in sparse signal modeling, so it is suitable for scenarios with a limited number of snapshots and an unknown number of sources [20,21]. An off-grid root sparse Bayesian learning (Root-SBL) method was proposed in [22], which greatly reduces the computational complexity of the OGSBL method. A new off-grid sparse Bayesian learning (SBL) method for nested arrays was also proposed in [23], which takes the noise variance as a part of the unknown signal of interest, and directly learns the unknown parameters through the sparse Bayesian learning model, which greatly improves the estimation performance.
However, the above methods do not consider various errors existing in practical applications, such as mutual coupling and array element position errors, which will seriously affect their effectiveness. To address the above problems, an SBL method was proposed to solve the problem of uncertain mutual coupling in [24,25,26] that deduces the theoretical expressions of all unknown parameters, e.g., noise variance, mutual coupling vector, and off-grid error vector, to achieve a good estimation performance. A sparse Bayesian reconstruction technique has also been considered to solve this estimation problem, with the joint impact of the simultaneous existence of amplitude phase and mutual coupling errors, and its results show that it has great engineering application value [27,28]. When considering the effect of impulsive noise, an off-grid sparse Bayesian learning method was also proposed in [29], which can effectively suppress impulsive noise while estimating off-grid DOA. In [30], Zhou proposes a reparameterized gamma process with random effects. Compared with the classical gamma process, the proposed model has a more intuitive physical interpretation. In addition, statistical inference for this model can be readily generated through the variational Bayesian algorithm.
This paper mainly studies how to reduce the impact of array element position errors on DOA estimation. One new DOA estimation method based on sparse Bayesian learning is proposed to overcome the presence of array position errors. The proposed method establishes the prior distribution of grid errors and array element position errors, as well as the joint probability density distribution function, and employs the expectation maximum (EM) method to iterate each unknown parameter. Finally, the effectiveness of this method is validated using the simulation results and the lake test.
The rest of this paper is organized as follows: In Section 2, the signal model is introduced. In Section 3, the theory derivation and algorithm flow are given. In Section 4, the simulation and lake test results are provided to validate the proposed method. Finally, the conclusion for this paper is provided in Section 5.

2. Signal Model

A model diagram of the coprime array is shown in Figure 1. Two uniform linear arrays with different spacings are nested into a coprime array: where subarray 1 has M2 sensors and its spacing is evenly distributed by M1; and subarray 2 has M1 sensors, and its spacing is evenly distributed by M1d. d represents the unit array element spacing, and this satisfies d = λ / 2 , where λ represents the signal wavelength. Additionally, the total number of elements is M = M 1 + M 2 1 . Based on the above assumption, the index for the collection of element positions is represented as follows:
Ω = [ M 1 m 2 | 0 m 2 M 2 ] [ M 2 m 1 | 0 m 1 M 1 ] .
Assuming that there are unknown errors at each position of the element Δ P = [ Δ P , 1 , Δ P , 2 , , Δ P , M ] T , the actual position of the element becomes Ω P = [ Ω 1 + Δ P , 1 , Ω 2 + Δ P , 2 , , Ω M + Δ P , M ] T . Given that K signals occur in different directions on the coprime array, and the incident direction is θ = [ θ 1 , θ 2 , , θ K ] , the receiving signal model of the coprime array is then as follows:
Y ( t ) = A ( θ ) s ( t ) + N ( t ) t = 1 , 2 , , T ,  
where Y ( t ) = [ y 1 ( t ) , y 2 ( t ) , , y M ( t ) ] T represents the data received by the coprime array. s ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s K ( t ) ] T represents K signal vectors. N ( t ) = [ n 1 ( t ) , n 2 ( t ) , , n M ( t ) ] T represents the noise vector received by the coprime array, which obeys the Gaussian distribution with a mean value of 0 and covariance σ n 2 I , where σ n 2 represents the power of noise. T represents the number of snapshots, and A ( θ , Ω P ) = [ a ( θ 1 , Ω P ) , , a ( θ K , Ω P ) ] is the manifold matrix of the array model; a ( θ k , Ω P ) is represented by:
a ( θ k , Ω P ) = [ e j 2 π ( Ω 1 + Δ P , 1 ) d sin θ k / λ , , e j 2 π ( Ω M + Δ P , M ) d sin θ k / λ ] T = ψ ( θ k , Δ P ) a ˜ ( θ k ) ,  
where a ˜ ( θ k ) = [ e j 2 π Ω 1 d sin θ k / λ , , e j 2 π Ω M d sin θ k / λ ] T denotes the array-steering vector of the k-th incident signal for the ideal condition; ψ ( θ k , Δ P ) = [ e j 2 π Δ P , 1 d sin θ k / λ , , e j 2 π Δ P , M d sin θ k / λ ] T represents the element position error vector factor of the coprime array; and denotes the Hadamard product operation. Therefore, the manifold matrix of the coprime array can be expressed as A ( θ , Ω P ) = Ψ ( θ , Δ P ) A ˜ ( θ ) , A ˜ ( θ ) = [ a ˜ ( θ 1 ) , , a ˜ ( θ K ) ] , Ψ ( θ , Δ P ) = [ ψ ( θ 1 , Δ P ) , , ψ ( θ K , Δ P ) ] .
The spatial interval [ π / 2 , π / 2 ] is evenly divided into N grid points, where N M > K , and the corresponding overcomplete angle set can be expressed as θ ¯ = [ θ ¯ 1 , , θ ¯ n , , θ ¯ N ] , so the expression under the sparse model is obtained as:
Y ( t ) = A ( θ ¯ , Ω P ) X ( t ) + N ( t ) t = 1 , 2 , , T ,  
where A ( θ ¯ , Ω P ) can be denoted as A ( θ ¯ , Ω P ) = [ a ( θ ¯ 1 , Ω P ) , a ( θ ¯ 2 , Ω P ) , , a ( θ ¯ K , Ω P ) ] , A ( θ ¯ , Ω P ) = Ψ ( θ ¯ , Δ P ) A ˜ ( θ ¯ ) , Ψ ( θ ¯ , Δ P ) = [ ψ ( θ ¯ 1 , Δ P ) , , ψ ( θ ¯ N , Δ P ) ] , A ˜ ( θ ¯ ) = [ a ˜ ( θ ¯ 1 ) , , a ˜ ( θ ¯ N ) ] . X ( t ) is the zero extension of the original signal s ( t ) , which has a value near the actual incident angle, and a value of 0 in other cases. However, the incident angle cannot be accurately incident on the preset grid point in the actual received data. Therefore, we introduced the off-grid error to the model, which was obtained with the help of the Taylor expansion as follows:
Φ ( Ω P , δ ) A ( θ ¯ , Ω P ) + B ( θ ¯ , Ω P ) diag ( δ ) ,  
where B ( θ ¯ , Ω P ) = [ b ( θ ¯ 1 , Ω P ) , , b ( θ ¯ n , Ω P ) , , b ( θ ¯ N , Ω P ) ] , 1 n N , b ( θ ¯ n , Ω P ) = a ( θ ¯ n , Ω P ) / θ ¯ n . δ = [ δ 1 , δ 2 , , δ N ] represents the angle grid error, and it satisfies δ [ ϑ / 2 , ϑ / 2 ] , where ϑ represents the grid spacing. diag ( ) represents the operation transforming the vector into a diagonal matrix. Therefore, the new discrete grid model can be expressed as:
Y ( t ) = Φ ( Ω P , δ ) X ( t ) + N ( t ) t = 1 , 2 , , T .  

3. Sparse Bayesian Learning Model

3.1. Model Parameter Distribution

(1) Distribution of noise and its precision: We can assume that the noise is Gaussian white noise with a distribution function as follows:
p ( N | 0 , α n ) = t = 1 T C N ( N t | 0 , α n I M ) ,  
where N t represents the t-th column of signal noise. α n represents noise precision and is defined as α n 1 σ n 2 . Noise precision follows gamma distribution, and its probability density distribution function can be denoted as:
p ( α n ) = Γ ( α n | 1 , b ) ,  
where b is the hyperparameter of the gamma distribution, Γ ( v ) = 0 + t v 1 e t d t , Γ ( x | v , z ) = Γ 1 ( v ) z v x v 1 exp ( z x ) .
Therefore, sparse model Y ( t ) is subject to C N ( Y ( t ) | Φ ( Ω P , δ ) X ( t ) , α n I M ) , and its probability distribution function can be given as:
p ( Y | X , Δ P , δ , α n ) = t = 1 T C N ( Y t | Φ ( Ω P , δ ) X t , α n I M ) ,  
where Y t represents the t-th column of the received signal Y , and X t represents the t-th column of sparse signal X .
(2) Distribution of signal and its precision: Assuming that the sparse signal X ( t ) is independent in different snapshots, and has a normal distribution with a mean of 0 and variance γ = [ γ 1 , , γ n , , γ N ] , its probability distribution function is:
p ( X | γ ) = t = 1 T C N ( X t | 0 , γ ) ,  
where γ n is signal precision, γ = [ γ 1 , , γ n , , γ N ] represents the sparse signal precision vector, which is subject to joint gamma distribution, and the probability distribution function is as follows:
p ( γ ) = n = 1 N Γ γ n ; 1 , c ,  
where c represents the hyperparameter of the gamma distribution.
(3) Grid error distribution: The grid error satisfies δ [ ϑ / 2 , ϑ / 2 ] , which follows a uniform distribution, and its probability distribution function is expressed as:
p ( δ ) = U [ ϑ / 2 , ϑ / 2 ] N .
(4) Distribution of array error and its precision: We can suppose that the errors in the array position Δ P are independent of each other, and they obey a normal distribution with mean 0 and variance ρ = [ ρ 1 , , ρ m , ρ M ] ; therefore, their probability distribution function is:
p ( Δ P | ρ ) = m = 1 M C N ( Δ P , m | 0 , ρ m ) ,  
where ρ m represents the array element error precision, ρ = [ ρ 1 , ρ 2 , , ρ M ] follows the joint gamma distribution, and its probability distribution function is:
p ( ρ ) = m = 1 M Γ ρ m ; 1 , e ,  
where e represents the hyperparameter of the gamma distribution.
According to Equations (8)–(14), the joint probability density function can be obtained as:
p ( X , Y , δ , Δ P , α n , γ , ρ ) = p ( Y | X , Δ P , δ , α n ) p ( X | γ ) p ( γ ) p ( α n ) p ( Δ P | ρ ) p ( ρ ) p ( δ ) .

3.2. Proposed Method

Since p ( X , δ , Δ P , α n , γ | Y ) cannot be directly calculated, the EM algorithm will be used to perform Bayesian learning derivation. Using the above preset model, the posterior probability distribution of sparse signal X can be given as:
p ( X | Y , δ , Δ P , α n , γ ) = p ( Y | X , δ , Δ P , α n ) p ( X | γ ) p ( Y , δ , Δ P , α n , γ ) p ( Y | X , δ , Δ P , α n ) p ( X | γ ) .
We treat the sparse signal X as a hidden variable whose posterior probability distribution also follows the Gaussian distribution with the mean μ and covariance Σ x , as follows:
Σ x = α n Φ H ( Ω P , δ ) Φ ( Ω P , δ ) + Λ 1 1 ,  
μ = α n Σ x Φ H ( Ω P , δ ) Y ,  
where Λ = diag ( γ ) .
The calculation of μ and Σ x in (17) and (18) requires the estimation of hyperparameters δ , Δ P , α n , and γ . These hyperparameter estimation results are obtained by maximizing p ( X , δ , Δ P , α n , γ | Y ) using the MAP method; therefore, we have the following:
δ ^ , Δ ^ P , α ^ n , γ ^ , ρ ^ = max L ( δ , Δ P , α n , γ , ρ ) = max E ln p ( X , Y , δ , Δ P , α n , γ , ρ ) ,  
where E { } denotes the expectation operation.
(1) Noise precision α n : Ignoring the posterior probability distribution function unrelated to noise precision in Equations (15) and (19), the likelihood function expression of noise precision is:
L ( α n ) = E ln p ( Y | X , Δ P , δ , α n ) p ( α n ) = M T ln α n π α n t = 1 T Y t Φ ( Δ P , δ ) X t 2 2 + ln b b α n .
Taking the derivative of Equation (20) and setting the derivative result to 0, we can calculate the following:
ln L ( α n ) α n = M T α n t = 1 T Y t Φ ( Ω P , δ ) X t 2 2 b = M T α n t = 1 T Y t Φ ( Ω P , δ ) μ t 2 2 T × Tr Σ x Φ H ( Ω P , δ ) Φ ( Ω P , δ ) b .
Therefore, the iterative expression of noise precision can be obtained as follows:
α ^ n = M T t = 1 T Y t Φ ( Ω P , δ ) μ t 2 2 + T × Tr Σ x Φ H ( Ω P , δ ) Φ ( Ω P , δ ) + b ,
where μ t represents the t-th column of the mean μ , Tr { } represents the trace operation of the matrix, and 2 represents the two-norm operation of the vector.
(2) Signal precision γ : Ignoring any parameters unrelated to signal precision, the likelihood function expression regarding signal precision can be written as follows:
L ( γ ) = E ln p ( X | γ ) p ( γ ) = T ln Λ π N t = 1 T X t H Λ X t + N ln c c n = 1 N γ n .
Taking the derivative of Equation (23) and setting the derivative result to 0, we calculate:
ln L ( γ ) γ n = T γ n t = 1 T X t , n H X t , n c .
Therefore, the iterative expression of signal precision can be obtained as follows:
γ ^ n = T c + t = 1 T X t , n H X t , n = T c + T Σ x ( n , n ) + t = 1 T μ n , t 2 , 1 n N ,  
where Σ x ( n , n ) denotes the t-th row and t-th column of the covariance matrix Σ x , And μ n , t denotes the n-th row and t-th column of the mean matrix μ .
(3) Array position error Δ P : Ignoring the posterior probability distribution function unrelated to array position error in Equations (15) and (19), the likelihood function expression for array position error is:
L ( Δ P ) = E ln p ( Y | X , Δ P , δ , α n ) p ( Δ P | ρ ) = α n t = 1 T Y t Φ X t 2 2 Δ P H diag ( ρ ) Δ P = α n t = 1 T Y t Φ μ t 2 2 α n T × Tr Φ H Φ Σ x Δ P H diag ( ρ ) Δ P
Setting the derivative of (26) as 0, the iterative expression of the array position error is obtained as follows:
Δ ^ P , m = α n Re t = 1 T Υ t H diag ( e m ) Ξ μ t α n T × Tr Re Φ ˜ H diag ( e m ) Ξ Σ x ρ m + α n t = 1 T μ t H Ξ H diag ( e m ) Ξ μ t + α n T × Tr Ξ H diag ( e m ) Ξ Σ x ,
where m = 1 , 2 , , M , Re represents taking the real part of a matrix. Please see Appendix A for the proof process.
(4) Array position error precision ρ : Ignoring the posterior probability distribution function unrelated to array position error precision, the likelihood function expression of array position error precision is:
L ( ρ ) = E ln p ( Δ P | ρ ) p ( ρ ) = ln diag ( ρ ) π M Δ P H diag ( ρ ) Δ P + M ln e e m = 1 M ρ m .
The derivative of Equation (28) is set to 0, and the iterative expression of the array position error precision is obtained as:
ρ ^ m = 1 Δ P , m 2 + e , m = 1 , 2 , , M .
(5) Grid error δ : Ignoring the posterior probability distribution function unrelated to grid error, the likelihood function expression for grid error is:
L ( δ ) = E ln p ( Y | X , Δ P , δ , α n ) p ( δ ) = α n T t = 1 T Y t Φ ( Ω P , δ ) μ t 2 2 + α n Tr Σ x Φ H ( Ω P , δ ) Φ ( Ω P , δ ) ,
where:
Y t Φ ( Ω P , δ ) μ t 2 2 = Y t A μ t B diag ( δ ) μ t 2 2 = δ H B H B μ t μ t H δ 2 Re diag ( μ t ) B H Y A μ t H δ + const ,
Tr Σ x Φ H ( Ω P , δ ) Φ ( Ω P , δ ) = 2 Re Tr { B H A Σ x diag ( δ ) } + Tr { diag ( δ ) Σ x diag ( δ ) B H B } + const = 2 Re Tr { B H A Σ x } δ + δ H Σ x B H B δ + const .
If we let the derivative of Equation (30) be equal to 0, the iterative expression of the grid error is obtained as follows:
δ ^ = P 1 v ,  
where δ ^ [ ϑ / 2 , ϑ / 2 ] , P N × N and v N × 1 are, respectively, represented by:
P = Re 1 T t = 1 T ( B H B ) ( μ t μ t H + Σ x )
v = Re 1 T t = 1 T diag ( μ t ) B H Y t A μ t H Re diag { B H A Σ x }
The estimated grid points are as follows:
θ ¯ = θ ¯ + δ ^ .
Based on the above theoretical analyses, the process of the proposed DOA estimation method based on sparse Bayesian learning presented in this paper is shown in Algorithm 1.
Algorithm 1. Sparse Bayesian learning with array position errors
Input: received data Y, Φ ( Ω P , δ )
Initialization: b = c = e = 10−2, ε = 10 3 , Iter = 0
Repeat:
  Iter = Iter + 1.
  Update Σ x and μ using Equations (17) and (18).
  Update α ^ n using Equation (22).
  Update γ ^ n using Equation (25).
  Update Δ ^ P , m using Equation (27).
  Update ρ ^ m using Equation (29).
Until convergence
Update δ using Equation (33).
Output: estimated θ ¯ , Σ x , and μ .

4. Simulation Analysis

To verify the excellent performance of the proposed method in a DOA estimation, we compared it with existing estimation methods such as the MUSIC method, OGSBL method [16], and SunFG method [31]. The simulation parameters used in this simulation experiment are shown in Table 1.

4.1. Simulation Results

We assume that the coprime array is composed of seven arrays. According to Equation (1), the ideal positions corresponding to the array elements are [ 0 , 3 , 5 , 6 , 9 , 10 , 12 ] λ / 2 , respectively, and the errors between the actual placement position and the ideal position of each array element are [ 0.12 , 0.17 , 0.13 , 0.09 , 0.11 , 0.14 , 0.06 ] λ / 2 , respectively.
Figure 2 shows the spatial spectrum image of the simulation experiment. As can be seen from Figure 2, all four methods can estimate the orientation of the signal source better. In the signal-independent direction, the curve changes of the proposed method and the MUSIC method are relatively gentle, while the curve changes of the OGSBL method and the SunFG method change sharply, which will cause greater interference with the actual signal incident direction and affect the judgment. Compared with the MUSIC method, the proposed method improves the power spectrum of the signal-independent direction by 5~10 dB, which shows that it can achieve a better performance than the MUSIC method. When we review the amplification details of the signal incidence direction, we find that the angle estimated by the proposed method is closer to the actual incidence angle, which is superior to the OGSBL, SunFG, and MUSIC methods. Therefore, the performance of the proposed method is optimal.

4.2. RMSE Performance Analysis

To investigate the influence of the SNR and the number of snapshots on the precision of the proposed method, root mean square error (RMSE) is used in this section as the only criterion for evaluating the accuracy, and we define the RSME expression as follows:
θ RSME = 1 Q K q = 1 Q k = 1 K θ ^ q k θ k
where Q represents the number of Monte Carlo simulation experiments; θ ^ q k denotes the k-th signal angle estimation result estimated by the q-th Monte Carlo simulation experiment.
To show the relation between the SNR and angle error, we set the number of snapshots to 500, and changed the SNR from −10 dB to 10 dB with a step of 2 dB, with the other conditions remaining unchanged. A total of 300 Monte Carlo simulation experiments were carried out. Figure 3 shows the simulation results for the relationship between the SNR and angle error.
It can be seen from Figure 3 that the angle RSMEs of the above four methods decrease with the increments of SNR, and the proposed method can achieve a better estimation performance in the whole SNR range. When the SNR is 10 dB, the RSMEs of the four methods are about 0.032° (the proposed method), 0.039° (OGSBL method), 0.095° (SunFG method), and 0.102° (MUSIC method). The angle error of the proposed method is lower than that of the other three methods. Therefore, it can be seen from the experimental results that the DOA estimation accuracy of the proposed method is better than that of the other three methods.
To show the relationship between the snapshots and the angle error, we set the SNR to 10 dB, and changed the snapshots from 50 to 500 with a step of 50, with other conditions remaining the same. A total of 300 Monte Carlo simulation experiments were also carried out. The simulation results regarding the relationship between the snapshots and the angle error are shown in Figure 4.
We can see that the angle RSMEs of the four methods decrease with the incremental changes in snapshots from Figure 4, and that the proposed method can achieve a better estimation performance in the whole snapshot range. When the number of snapshots is 500, the RSMEs of the four methods are about 0.086° (the proposed method), 0.105° (OGSBL method), 0.18° (SunFG method), and 0.29° (MUSIC method). The angle error of the proposed method is lower than that of the other three methods. Therefore, the experimental results prove that the DOA estimation accuracy of the proposed method is better than that of the other three methods.

4.3. Resolution Probability Performance Analysis

The probability of successful angle resolution was defined as the absolute deviation between the estimated DOA value, and an actual incidence angle of less than 0.2°, which was considered a successful resolution: namely, θ ^ q k θ k 0.2 ° .
While conducting the 300 Monte Carlo simulation experiments, the signal-to-noise ratio (SNR) varied from −10 dB to 10 dB in 2 dB increments while keeping the other conditions constant, with a shutter speed of 500 frames per second. Figure 5 illustrates the simulation results for the relationship between the SNR and angle resolution probability.
From this graph, we can see that the angle resolution capabilities of the four methods mentioned above increase with increases in the signal-to-noise ratio. Moreover, the method proposed in this paper exhibits a good estimation performance across the entire signal-to-noise ratio range. At a signal-to-noise ratio of 8 dB, the probability of successful angle resolution for the four methods is approximately 1 for the method proposed in this paper, compared to 0.97 for the OGSBI method, 0.925 for the SunFG method, and 0.88 for the MUSIC method. The method proposed in this paper achieves a higher accuracy in terms of angle resolution compared to the other three methods. Therefore, the experimental results indicate that the proposed method has a superior DOA estimation resolution compared to the other three methods.
Keeping other conditions unchanged, with a signal-to-noise ratio of 5 dB, the number of snapshots varied from 50 to 500 in increments of 50. A Monte Carlo simulation experiment was conducted for 300 iterations, and the simulation results of the varying number of snapshots and the probability of angle resolution are shown in Figure 6.
From this figure, it can be seen that the angle resolution capabilities of the four methods mentioned above increase with the increases in the number of snapshots. Moreover, the method proposed in this paper demonstrates a good estimation performance across the entire range of snapshot numbers. When the number of snapshots is 500, the probabilities of successful angle resolution for the four methods are approximately 0.95 (proposed method), 0.89 (OGSBI method), 0.8 (SunFG method), and 0.72 (MUSIC method). The proposed method exhibits higher accuracy in angle resolution compared to the other three methods. Therefore, the experimental results indicate that the proposed method has better DOA estimation resolution than the other three methods.

4.4. Running Time Analysis

Considering the influence of the snapshot number on the operation time, we kept all other conditions unchanged, and changed the number of snapshots. The snapshot number started from 50 and increased in steps of 50 until the simulation reached 300 iterations. We ran it 100 times to take the average value and record the change in the operation time between the different methods. The run time curves versus the number of snapshots across the different methods are shown in Figure 7.
From Figure 7, we can see that the MUSIC algorithm has the lowest computation time and the highest computational efficiency among the four methods. The main reason for this is that this method does not require iterative operations and can be achieved by performing eigen decomposition of the covariance matrix. The computation times of the proposed method, OGSBI method, and SunFG method increase gradually with the increases in snapshot number, and the computation time of the proposed method is always smaller than that of the other two methods; so, the computational efficiency of the proposed method is superior to that of the SunFG method and OGSBI method.

4.5. The Lake Experimental Data Analysis

To verify the effectiveness and reliability of the above-mentioned methods in practical applications, we conducted a lake experiment using a fiber-optic hydrophone array. The structure of the coprime array in the lake experiment is shown in Figure 8. The received signal of the fiber-optic hydrophone array is shown in Figure 9. We processed the lake experimental data received by the fiber-optic hydrophone array to verify the resolution performance of the proposed method. During the actual processing, we selected 60,000 valid data points from Figure 9 for analysis. Then, the experimental data were processed using three methods, and the spatial spectrum images of the three methods in the underwater environment are shown in Figure 10. The SunFG method could not successfully estimate the direction of the incident signal in this underwater experiment, so the processing results for this method are not shown in Figure 10.
From Figure 10, we can see that the MUSIC method is relatively stable in the signal-independent direction, but its power spectrum curve has a wide-angle range in the incident angle direction; thus, the MUSIC method has a bad performance in the estimation of adjacent directions, and its resolution is relatively poor compared to the other two methods. The power of the OGSBL method and the proposed method is concentrated in the direction of the incident angle, so these two methods have better resolution. However, the OGSBL method generates four pseudo-peaks and the maximum value reaches −20 dB, which will interfere with its application in practical engineering, while the proposed method generates only two pseudo-peaks and its maximum peak only reaches −30 dB, which is smaller than that of the OGSBL method. Therefore, our proposed method is far superior to the OGSBL method. Overall, the performance of the proposed method is better than that of the other methods.
We framed the receiving signal and used the above three methods to estimate the DOA of the incident signal in this underwater experiment. The bearing time recordings are shown in Figure 11, Figure 12 and Figure 13.
From the bearing time recordings of the above three methods shown in Figure 11 Figure 12 and Figure 13, we can see that all of the methods can obtain a linear course trajectory of the target. Considering the overall effect, the MUSIC method has the widest course trajectory and the worst resolution. The OGSBL method can also estimate a clear trajectory curve, but produces three severely interfering pseudo-peaks, which require pseudo-peak identification and greatly reduce the application of this method in practical engineering. The proposed method can estimate a clear linear trajectory curve, and only generates two pseudo-peaks, and the pseudo-peak value reaches −30 dB, which is better than the OGSBL method. Therefore, the proposed method can obtain the best estimation performance among the three methods.

5. Conclusions

In this paper, we analyzed and studied a DOA estimation method in which we employed sparse Bayesian learning to overcome array position errors. The proposed method substitutes the problem of array position errors into the signal receiving model to obtain the joint probability density distribution function. We used the EM method to optimize and derive parameters such as the off-grid error and array position error and determined the DOA through the spatial spectral image. The simulation experiment results and lake experiment results show that the DOA estimation performance of the proposed method is better than that of three other classical methods, i.e., OGSBL, MUSIC, and SunFG. In cases with fewer snapshots and a low SNR, the proposed method can still obtain robust DOA estimation performance, and its estimation accuracy is better than that of the three other methods. The proposed method effectively solved the problem of element position errors and improved the robustness of the mutual prime array greatly, so it has great value in practical application.

Author Contributions

Conceptualization, Y.T. and X.W. (Xuhu Wang); methodology, Y.T.; validation, Q.Z. and X.W. (Xuhu Wang); formal analysis, Y.T.; investigation, Y.T.; resources, X.W. (Xuhu Wang); data curation, L.D.; writing—original draft preparation, Y.T.; writing—review and editing, X.W. (Xinjie Wang) and Q.F.; visualization, X.W. (Xuhu Wang); supervision, X.W. (Xuhu Wang); project administration, Q.Z.; funding acquisition, X.W. (Xuhu Wang). All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded in part by the National Natural Science Foundation of China, grant number 62171247, and in part by the Shandong Provincial Natural Science Foundation of China, grant number ZR2022MF273.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors wish to thank Qingdao University of Technology and their intelligent perception and wireless communication research team.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Proof of Equation (25).
The array position error is very small, so we performed a Taylor series expansion to the array position error vector factor, from which we can obtain
ψ ( θ k , Δ P ) = [ e j 2 π Δ P , 1 d sin θ k / λ , , e j 2 π Δ P , M d sin θ k / λ ] T 1 + j 2 π Δ P , 1 d sin θ k λ , , 1 + j 2 π Δ P , M d sin θ k λ T = 1 M × 1 + Δ P j 2 π d sin θ k / λ
The over-complete basis matrix Φ ( θ ¯ , Ω P ) can be transformed into
Φ ( Ω P , δ ) = A ( θ ¯ , Ω P ) + B ( θ ¯ , Ω P ) diag ( δ ) = Φ ˜ ( δ ) + diag ( Δ P ) Φ ˜ ( δ ) diag ( j 2 π d sin θ ¯ / λ ) + diag ( Δ P ) A ˜ ( θ ¯ ) diag ( δ ) diag ( j 2 π d cos θ ¯ / λ ) = Φ ˜ ( δ ) + diag ( Δ P ) Ξ
where Ξ = Φ ˜ ( δ ) diag ( j 2 π d sin θ ¯ / λ ) + A ˜ ( θ ¯ ) diag ( δ ) diag ( j 2 π d cos θ ¯ / λ ) , Φ ˜ ( δ ) A ˜ ( θ ¯ ) + B ˜ ( θ ¯ ) diag ( δ ) , B ˜ ( θ ¯ ) = [ b ˜ ( θ ¯ 1 ) , , b ˜ ( θ ¯ n ) , , b ˜ ( θ ¯ N ) ] , and b ˜ ( θ ¯ n ) = a ˜ ( θ ¯ n ) / θ ¯ n .
By means of the definition in Equation (35), the first term of Equation (24) can be converted into
t = 1 T Y t Φ ( Ω P , δ ) μ t 2 2 = t = 1 T Y t Φ ˜ ( δ ) + diag ( Δ P ) Ξ μ t 2 2 = t = 1 T Υ t diag ( Δ P ) Ξ μ t 2 2
where Υ t = Y t Φ ˜ ( δ ) μ t . Taking the derivative of the element position error of equation in Equation (35), the result can be obtained as
t = 1 T Υ t diag ( Δ P ) Ξ μ t 2 2 Δ P , m = Δ P , m t = 1 T 2 Re Υ t H diag ( Δ P ) Ξ μ t + μ t H Ξ H diag ( Δ P Δ P ) Ξ μ t = 2 Re t = 1 T Υ t H diag ( e m ) Ξ μ t + 2 Δ P , m t = 1 T μ t H Ξ H diag ( e m ) Ξ μ t
The second term of Equation (24) can also be converted into
Tr Φ H ( Ω P , δ ) Φ ( Ω P , δ ) Σ x = Tr Φ ˜ ( δ ) + diag ( Δ P ) Ξ H Φ ˜ ( δ ) + diag ( Δ P ) Ξ Σ x = Tr Φ ˜ H ( δ ) Φ ˜ ( δ ) Σ x + Tr Φ ˜ H ( δ ) diag ( Δ P ) Ξ Σ x + Tr Ξ H diag ( Δ P ) Φ ˜ ( δ ) Σ x + Tr Ξ H diag ( Δ P Δ P ) Ξ Σ x
If we take the derivative of the element position error of Equation (A1), the result can be derived as
Δ P , m Tr Φ H ( Ω P , δ ) Φ ( Ω P , δ ) Σ x = Tr Φ ˜ H ( δ ) diag ( Δ P ) Ξ Σ x Δ P , m + Tr Ξ H diag ( Δ P ) Φ ˜ ( δ ) Σ x Δ P , m + Tr Ξ H diag ( Δ P Δ P ) Ξ Σ x Δ P , m = 2 Tr Re Φ ˜ H ( δ ) diag ( Δ P ) Ξ Σ x Δ P , m + Tr Ξ H diag ( Δ P Δ P ) Ξ Σ x Δ P , m = 2 Tr Re Φ ˜ H ( δ ) diag ( e m ) Ξ Σ x + 2 Δ P , m × Tr Ξ H diag ( e m ) Ξ Σ x
Setting the result of Equation (24) as 0, the iterative expression of the array element position error can be obtained as Equation (25). The proof is thus completed. □

References

  1. Ai, W.S.; Zheng, E.M.; Bao, K.K. A method of line spectrum extraction based on target radiated spectrum feature and its post-processing. J. Syst. Eng. Electron. 2021, 32, 1381–1393. [Google Scholar]
  2. Xie, N.; Guo, W.; Zhang, L.; Wang, H. High-Resolution DOA Estimation Method in MIMO Radar. Wirel. Pers. Commun. 2016, 91, 219–236. [Google Scholar] [CrossRef]
  3. Chen, T.; Huang, Q.; Zhang, L.; Fang, Y. Direction of Arrival Estimation Using Distributed Circular Microphone Arrays. In Proceedings of the 2018 14th IEEE International Conference on Signal Processing (ICSP), Beijing, China, 12–16 August 2018; pp. 182–185. [Google Scholar]
  4. Zhang, H.; Zhang, F.; Chen, M.; Fu, H. Non-uniform linear sonar array-based DOA estimation. In Proceedings of the 33rd Chinese Control Conference, Nanjing, China, 28–30 July 2014; pp. 7240–7243. [Google Scholar]
  5. Chen, Y.H.; Yan, L.F.; Han, C. Millidegree-level direction-of-arrival (DOA) estimation and tracking for terahertz wireless communications. In Proceedings of the 2020 17th Annual IEEE International Conference on Sensing, Communication, and Networking (SECON), Como, Italy, 22–25 June 2020; pp. 1–9. [Google Scholar]
  6. Li, X.Y.; Yang, G.Q.; Gu, Y.K. Simulation analysis of MUSIC algorithm of array signal processing DOA. In Proceedings of the International Conference on Automatic Control and Artificial Intelligence (ACAI 2012), Xiamen, China, 3–5 March 2012; pp. 1838–1841. [Google Scholar]
  7. He, Z.; Li, Y.; Xiang, J. A modified root-music algorithm for signal DOA estimation. J. Syst. Eng. Electron. 1999, 10, 42–47. [Google Scholar]
  8. Zheng, G.; Tang, J.; Yang, X. ESPRIT and Unitary ESPRIT Algorithms for Coexistence of Circular and Noncircular Signals in Bistatic MIMO Radar. IEEE Access 2016, 4, 7232–7240. [Google Scholar] [CrossRef]
  9. Lei, R.; Wang, J.; Li, H.; Wang, Z. DOA estimation algorithm of MIMO system. In Proceedings of the IET International Conference on Communication Technology and Application (ICCTA 2011), Beijing, China, 14–16 October 2011; pp. 150–154. [Google Scholar]
  10. Pal, P.; Vaidyanathan, P. Coprime sampling and the music algorithm. In Proceedings of the 2011 Digital Signal Processing and Signal Processing Education Meeting (DSP/SPE), Sedona, AZ, USA, 4–7 January 2011; pp. 289–294. [Google Scholar]
  11. Zhou, C.W.; Shi, Z.G.; Gu, Y.J.; Shen, X. DECOM: DOA estimation with combined MUSIC for coprime array. In Proceedings of the 2013 International Conference on Wireless Communications and Signal Processing, Hangzhou, China, 24–26 October 2013; pp. 1–5. [Google Scholar]
  12. Weng, Z.Y.; Peter, M.D. A search-free DOA estimation algorithm for coprime arrays. Digit. Signal Process. 2014, 24, 27–33. [Google Scholar] [CrossRef]
  13. Zhou, C.W.; Zheng, H.; Gu, Y.J.; Wang, Y.; Shi, Z. Research progress on coprime array signal processing: Direction-of arrival estimation and adaptive beamforming. J. Radars 2019, 8, 558–577. [Google Scholar]
  14. Li, J.F.; Li, D.; Jiang, D.F.; Zhang, X. Extended-aperture unitary root MUSIC-based DOA estimation for coprime array. IEEE Commun. Lett. 2018, 22, 752–755. [Google Scholar] [CrossRef]
  15. Malioutov, D.; Cetin, M.; Willsky, A.S. A sparse signal reconstruction perspective for source localization with sensor arrays. IEEE Trans. Signal Process. 2005, 53, 3010–3022. [Google Scholar] [CrossRef]
  16. 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]
  17. Wu, X.H.; 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]
  18. Das, A. Real-Valued Sparse Bayesian Learning for Off-Grid Direction-of-Arrival (DOA) Estimation in Ocean Acoustics. IEEE J. Ocean. Eng. 2021, 46, 172–182. [Google Scholar] [CrossRef]
  19. Liu, Y.; Dong, N.; Zhang, X.; Zhao, X.; Zhang, Y.; Qiu, T. DOA Estimation for Massive MIMO Systems with Unknown Mutual Coupling Based on Block Sparse Bayesian Learning. Sensors 2022, 22, 8634. [Google Scholar] [CrossRef] [PubMed]
  20. Ge, S.D.; Fan, C.Y.; Wang, J.; Huang, X. Robust adaptive beamforming based on sparse Bayesian learning and covariance matrix reconstruction. IEEE Commun. Lett. 2022, 26, 1893–1897. [Google Scholar] [CrossRef]
  21. Dai, J.S.; So, C.H. Real-valued sparse Bayesian learning for DOA estimation with arbitrary linear arrays. IEEE Trans. Signal Process. 2021, 69, 4977–4990. [Google Scholar] [CrossRef]
  22. Dai, J.S.; 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]
  23. Chen, F.F.; Dai, J.S.; Hu, N.; Ye, Z. Sparse Bayesian learning for off-grid DOA estimation with nested arrays. Digit. Signal Process. 2018, 82, 187–193. [Google Scholar] [CrossRef]
  24. 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] [CrossRef]
  25. Wang, X.H.; Bai, H.D.; Zhang, Q.F.; Yuan, N. Off-grid DOA estimation based on sparse Bayesian learning under interaction among array elements. J. Vib. Shock. 2022, 41, 303–312. [Google Scholar]
  26. Dai, J.; Hu, N.; Xu, W.; Chang, C. Sparse Bayesian Learning for DOA Estimation with Mutual Coupling. Sensors 2015, 15, 26267–26280. [Google Scholar] [CrossRef]
  27. Wang, D.; Gao, W.G.; Wu, Z.D. Array amplitude-phase and mutual coupling error joint correction method based on sparse Bayesian. J. Commun. 2022, 43, 112–120. [Google Scholar]
  28. Dai, F.; Fu, H.; Hong, L.; Li, L.; Liu, H. Off-Grid Error and Amplitude–Phase Drift Calibration for Computational Microwave Imaging with Meta surface Aperture Based on Sparse Bayesian Learning. In IEEE Transactions on Geoscience and Remote Sensing; IEEE: New York City, NY, USA, 2022; Volume 60, pp. 1–14. [Google Scholar]
  29. Ma, J.; Zhang, J.; Yang, Z.; Qiu, T. Off-Grid DOA Estimation Using Sparse Bayesian Learning for MIMO Radar under Impulsive Noise. Sensors 2022, 22, 6268. [Google Scholar] [CrossRef] [PubMed]
  30. Zhou, S.; Xu, A.; Tang, Y.; Shen, L. Fast Bayesian Inference of Reparameterized Gamma Process with Random Effects. In IEEE Transactions on Reliability; IEEE: New York City, NY, USA, 2023; pp. 1–14. [Google Scholar]
  31. Sun, F.G.; Wu, Q.H.; Sun, Y.M.; Lan, P. An iterative approach for sparse direction-of-arrival estimation in co-prime arrays with off-grid targets. Digit. Signal Process. 2017, 61, 35–42. [Google Scholar] [CrossRef]
Figure 1. Model diagram of the coprime array.
Figure 1. Model diagram of the coprime array.
Mathematics 12 00545 g001
Figure 2. Spatial spectra of various methods.
Figure 2. Spatial spectra of various methods.
Mathematics 12 00545 g002
Figure 3. Relationship between SNR and angle error.
Figure 3. Relationship between SNR and angle error.
Mathematics 12 00545 g003
Figure 4. Relationship between snapshots and angle error.
Figure 4. Relationship between snapshots and angle error.
Mathematics 12 00545 g004
Figure 5. Relationship between SNR and resolution probability.
Figure 5. Relationship between SNR and resolution probability.
Mathematics 12 00545 g005
Figure 6. Relationship between the number of snapshots and resolution probability.
Figure 6. Relationship between the number of snapshots and resolution probability.
Mathematics 12 00545 g006
Figure 7. Relationship between the number of snapshots and computation time.
Figure 7. Relationship between the number of snapshots and computation time.
Mathematics 12 00545 g007
Figure 8. Hydrophone array configuration structure in the lake experiment.
Figure 8. Hydrophone array configuration structure in the lake experiment.
Mathematics 12 00545 g008
Figure 9. Graph of data collected via four fiber-optic hydrophones in the lake experiment.
Figure 9. Graph of data collected via four fiber-optic hydrophones in the lake experiment.
Mathematics 12 00545 g009
Figure 10. Spatial spectrum of the data received via fiber-optic hydrophone in the lake experiment.
Figure 10. Spatial spectrum of the data received via fiber-optic hydrophone in the lake experiment.
Mathematics 12 00545 g010
Figure 11. The bearing time recording of the MUSIC method.
Figure 11. The bearing time recording of the MUSIC method.
Mathematics 12 00545 g011
Figure 12. The bearing time recording of the OGSBL method.
Figure 12. The bearing time recording of the OGSBL method.
Mathematics 12 00545 g012
Figure 13. The bearing time recording of the proposed method.
Figure 13. The bearing time recording of the proposed method.
Mathematics 12 00545 g013
Table 1. Experimental simulation conditions.
Table 1. Experimental simulation conditions.
ParameterValue
The number of arrays M7
The number of subarray1 M23
The number of subarray2 M15
The number of snapshots T500
Signal-to-noise ratio 10 dB
The number of signals K2
The incidence angle θ −27.32° and 17.75°
The space between arrays d0.5 wavelength
The grid interval ϑ
The search interval of the MUSIC algorithm0.5°
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tian, Y.; Wang, X.; Ding, L.; Wang, X.; Feng, Q.; Zhang, Q. A New Sparse Bayesian Learning-Based Direction of Arrival Estimation Method with Array Position Errors. Mathematics 2024, 12, 545. https://doi.org/10.3390/math12040545

AMA Style

Tian Y, Wang X, Ding L, Wang X, Feng Q, Zhang Q. A New Sparse Bayesian Learning-Based Direction of Arrival Estimation Method with Array Position Errors. Mathematics. 2024; 12(4):545. https://doi.org/10.3390/math12040545

Chicago/Turabian Style

Tian, Yu, Xuhu Wang, Lei Ding, Xinjie Wang, Qiuxia Feng, and Qunfei Zhang. 2024. "A New Sparse Bayesian Learning-Based Direction of Arrival Estimation Method with Array Position Errors" Mathematics 12, no. 4: 545. https://doi.org/10.3390/math12040545

APA Style

Tian, Y., Wang, X., Ding, L., Wang, X., Feng, Q., & Zhang, Q. (2024). A New Sparse Bayesian Learning-Based Direction of Arrival Estimation Method with Array Position Errors. Mathematics, 12(4), 545. https://doi.org/10.3390/math12040545

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