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
, where
represents the signal wavelength. Additionally, the total number of elements is
. Based on the above assumption, the index for the collection of element positions is represented as follows:
Assuming that there are unknown errors at each position of the element
, the actual position of the element becomes
. Given that
K signals occur in different directions on the coprime array, and the incident direction is
, the receiving signal model of the coprime array is then as follows:
where
represents the data received by the coprime array.
represents
K signal vectors.
represents the noise vector received by the coprime array, which obeys the Gaussian distribution with a mean value of 0 and covariance
, where
represents the power of noise. T represents the number of snapshots, and
is the manifold matrix of the array model;
is represented by:
where
denotes the array-steering vector of the
k-th incident signal for the ideal condition;
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
,
,
.
The spatial interval
is evenly divided into
N grid points, where
, and the corresponding overcomplete angle set can be expressed as
, so the expression under the sparse model is obtained as:
where
can be denoted as
,
,
,
.
is the zero extension of the original signal
, 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:
where
,
,
.
represents the angle grid error, and it satisfies
, where
represents the grid spacing.
represents the operation transforming the vector into a diagonal matrix. Therefore, the new discrete grid model can be expressed as:
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 , respectively, and the errors between the actual placement position and the ideal position of each array element are , 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:
where
Q represents the number of Monte Carlo simulation experiments;
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, .
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.