(ProQuest: ... denotes non-US-ASCII text omitted.)
Academic Editor:Shen Yin
Department of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China
Received 17 January 2014; Accepted 15 February 2014; 27 March 2014
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The methods of the direction of arrival (DOA) estimation have received considerable attention because of their importance in a great variety of applications such as radar, sonar, and communication [1-5]. The subspace-based methods include multiple signal classification (MUSIC) algorithm [6] and estimation of signal parameters via rotational invariance techniques (ESPRIT) algorithm [7]. However, these high resolution methods rely on an exact characterization of the array manifold which is inevitably influenced by the array perturbations such as mutual coupling errors [8], gain and phase errors [9], and sensor position perturbation errors [10]. The performance of the direction finding methods will be greatly degraded with the array perturbations, and their applications are limited in practice. Therefore, it is essential to develop DOA estimation methods in the presence of array errors.
In recent decades, many researchers attempt to estimate the signal DOAs in the presence of array errors. To calibrate errors of mutual coupling, a number of calibration algorithms [11-18] have been proposed. A maximum-likelihood approach is developed in [11] by setting calibration sources with known locations. An alternative minimization procedure is presented in [12], but the multidimensional search of the associated nonlinear optimization is computationally complicated and convergence can not be guaranteed. DOAs and mutual coupling matrix are estimated in [13, 14] by utilizing rank-reduction (RARE) estimator based on the banded symmetric Toeplitz structural characteristics of the mutual coupling matrix. The methods developed in [15, 16] are based on auxiliary sensors for self-calibration of mutual coupling. Unfortunately, spectral peak search technique is required in these methods, which are limited in practical direction finding system because of multiple calculations. The methods in [17, 18] make use of auxiliary sensors similarly without requiring spectral peak search.
Besides researching the problem of mutual coupling errors, sensor gain and phase errors calibration methods [11, 19-25] have also been thoroughly investigated in the open literature. The active calibration method in [11] can also calibrate sensor gain and phase errors. The method in [19] can estimate the DOAs of signals and array parameters simultaneously based on the assumption that array perturbations are slight, but it suffers from suboptimal convergence caused by the joint iteration between DOA estimation and array parameter estimation. In [20], Liu et al. have proposed an eigenstructure DOA estimation method based on the eigendecomposition of a covariance matrix, which is constructed by dot product of array output and its conjugate. Shenghong et al. try to improve this eigenstructure DOA estimation method in [21]. However, at least two signals are spatially separated in the above two methods. In addition, the problem of DOA estimation using partially calibrated sensors or auxiliary sensors has been studied. The rank-reduction method [22-24] is typical, but it requires an additional spectral search. The root-RARE algorithm in [17] is computationally more efficient but also has its limitations. A simple and efficient method based on the conventional ESPRIT algorithm is proposed in [25], which can only calibrate the gain and phase errors.
In practical situations, various sources of array errors may be present simultaneously, but the aforesaid calibration methods are capable of analyzing statistical performance of circumstances where only one array error model presents. Assuming that sensor locations are known, array calibration can be considered as estimating the DOAs in the presence of sensor gain and phase errors and mutual coupling errors. Generally, calibration methods can be classified as active calibration methods [26-29] and self-calibration methods [30-33], which also have the above described problems.
In this paper, we propose a direction finding method in the presence of sensor gain and phase errors as well as mutual coupling errors. By applying a group of auxiliary sensors, we can eliminate mutual coupling errors by the inherent mechanism of the proposed method, assuming that DOA estimation is independent of mutual coupling coefficients. Then DOAs as well as gain and phase coefficients can be jointly estimated in closed form based on extended least-squares ESPRIT (LS-ESPRIT) algorithm. Finally, mutual coupling coefficients can also be estimated by utilizing a new decoupling method. The proposed method requires no spectral search or iteration, so the complexity is low, which makes it simpler but more efficient, and it can realize the estimation of DOAs, gain/phase error coefficients, and mutual coupling coefficients simultaneously.
The rest of the paper is organized as follows. Section 2 illustrates problem formulation. In Section 3, the proposed method for estimation of DOA, sensor gain and phase error coefficients, and mutual coupling coefficients is described, and some discussions are given. Simulation results are presented in Section 4 to evaluate advantages of the proposed method over some representative calibration methods. At last, conclusions are drawn in Section 5.
2. Problem Formulation
Consider a uniform linear array (ULA) consisting of M identical sensors with the space d between neighboring sensors and N uncorrelated narrowband far-field signals s1 (t),s2 (t),...,sN (t) that arrive at the array from directions θ1 ,θ2 ,...,θN . Assuming that signals and noises are stationary, zero mean, uncorrelated random processes and noises are uncorrelated with variance σ2 ; the output of the array is [figure omitted; refer to PDF] where x(t) is the M×1 array output vector, x(t)=[x1 (t),x2 (t),...,xM (t)]T , s(t)=[s1 (t),s2 (t),...,sN (t)]T is the N×1 vector of signal waveforms, and n(t)=[n1 (t),n2 (t),...,nM (t)]T is the M×1 sensor noise vector that is commonly assumed to be the additive white Gaussian noise (AWGN). Moreover, A=[a(θ1 ),a(θ2 ),...,a(θN )] is steering matrix in which the i th signals direction manifold vector is denoted by a(θi )=[1,βi1 ,βi2 ,...,βiN-1]T with βi =exp...{-j2dπsin(θi )/λ} , and i=1,2,...,N . λ denotes carrier wavelength.
According to (1), the array covariance matrix of the array output is [figure omitted; refer to PDF] where RS =E{s(t)sH (t)} is signal covariance matrix, E{·} represents the statistic expectation, and σ2 and I denote noise variance and M×M identity matrix, respectively.
The eigendecomposition of R is expressed as [figure omitted; refer to PDF] where US is signal subspace matrix containing eigenvectors with the largest eigenvalues and UV is noise subspace matrix containing eigenvectors with the smallest eigenvalues. ΣS and ΣV are diagonal matrices associated with signal power and noise power, respectively. In the case of finite snapshots, the array covariance matrix can be estimated as [figure omitted; refer to PDF] where L is snapshot number. After the eigendecomposition of array covariance matrix is accomplished, signal subspace matrix U^S and noise subspace matrix U^V can be estimated.
Considering mutual coupling errors and gain and phase uncertainties of the array, the corresponding array output vector is rewritten by [figure omitted; refer to PDF] where C and Γ denote mutual coupling matrix (MCM) and gain and phase matrix (GPM), respectively, which are generally considered independent of direction.
As we know, mutual coupling coefficients of sensors are inversely related to their distance. On the basis of ULA structure, MCM can be expressed as a banded symmetric Toeplitz matrix. Suppose there are p nonzero complex coefficients in MCM, with p<M . The 1×M coefficient vector can be written as [figure omitted; refer to PDF] where cm , m=1,...,p-1 is the mutual coupling coefficient. The corresponding M×M matrix C is given by [figure omitted; refer to PDF]
Actually, mutual coupling errors and gain and phase errors coexist. Assuming that the gain and phase of M sensor are uncertain, gain and phase coefficients can be presented by ρ1 ,...,ρM and [straight phi]1 ,...[straight phi]M , respectively. M×1 gain and phase vector is defined as [figure omitted; refer to PDF] and M×M GPM can be written as [figure omitted; refer to PDF]
Here, we are interested in the problem of how to estimate the DOAs {θi}i=1N in the presence of sensors gain and phase errors and mutual coupling errors and how to obtain their coefficients, respectively.
3. Proposed Method
3.1. Array Model
Auxiliary sensors with known gain and phase are utilized to estimate DOAs, gain and phase coefficients, and mutual coupling coefficients on the basis of original M sensors. First, we add q sensors in front of M sensors ULA and then add p-1 auxiliary sensors on each side of the M+q sensors array. Thus a new array with a total of M+q+2(p-1) sensors is formed, and the space between all neighboring sensors is d . Here p is the mutual coupling degree of freedom without generality loss as mentioned above, which is assumed to be fixed. A detailed analysis will be performed in Section 3.5 concerning how the number of accurate calibration sensors q affects algorithm performance. Besides, the array with a total of Mw =M+q+2(p-1) sensors is called whole array and that with Mm =M+q middle sensors is called middle array.
Assuming that q of sensors in the middle array are auxiliary sensors with known gain and phase, we can divide the middle array into two overlapping subarrays: one comprising first Mm -1 sensors where q of them are well calibrated and the other comprising last Mm -1 sensors where q-1 of them are well calibrated. The array geometry is shown in Figure 1. Thus DOAs and gain and phase coefficients can be estimated by extended LS-ESPRIT algorithm which will be introduced below. Meanwhile, mutual coupling coefficients can be estimated by using the receiving data of the whole array.
Figure 1: Array geometry. Solid points represent uncertain sensors and circles represent auxiliary sensors.
[figure omitted; refer to PDF]
Let C- and C_ represent MCM of whole array and middle array, respectively. Γ- and Γ_ denote their GPM, respectively. More precisely, gain and phase vectors are expressed as η- and η_ , respectively.
We define matrix P as [figure omitted; refer to PDF] where 0Mm ×(p-1) is an Mm ×(p-1) zero matrix and IMm ×Mm is an Mm ×Mm unit matrix. Thus the MCM C- of the whole array can be expressed as [figure omitted; refer to PDF]
The MCM C_ of the middle array is written by [figure omitted; refer to PDF]
Taking the unknown gain and phase uncertainties into account, their gain and phase vectors can be denoted by [figure omitted; refer to PDF] respectively, where 1 is a vector with all elements equal to one. Here gain and phase matrix can be written as [figure omitted; refer to PDF]
3.2. Joint Estimation of DOA and Gain and Phase Coefficient
Considering auxiliary sensors, the output of the middle array can be expressed as [figure omitted; refer to PDF] where A-=[a-(θ1 ),a-(θ2 ),...,a-(θN )] , a-(θ)=[1,β1 ,β2 ,...,βMw -1]T , and Mw =Mm +2(p-1) . Consequently, the steering matrix is given by [figure omitted; refer to PDF]
Substituting β=exp...{-j2dπsin(θi )/λ} into (17), we have [figure omitted; refer to PDF] where F=βp-1 (2∑m=1p-1cm cos...(2nπdsinθ/λ)+1) is a scalar function related to mutual coupling coefficient cm and the direction θ and a_(θ)=[1,β,...,βMm -1 ] is the ideal steering vector of the middle array.
According to the above analysis, MCM in the steering vector of the middle array is converted to the scalar function related to mutual coupling coefficients, thus realizing decoupling. Therefore, the steering vectors of two subarrays of the middle array are given by [figure omitted; refer to PDF] where A1 and A2 denote the ideal steering vectors of these two subarrays, respectively, and Γ1 and Γ2 are GPM of these two subarrays, respectively, which can be written as [figure omitted; refer to PDF]
The relationship between A1 and A2 is [figure omitted; refer to PDF] where Φ=diag...{ej2πdsinθ1 /λ ,ej2πdsinθ2 /λ ,...,ej2πdsinθN /λ } .
In the case that gain and phase errors and mutual coupling errors exist simultaneously, the received data of two subarrays are stacked into a column [figure omitted; refer to PDF] where A~=[Γ1 FA1Γ2 FA2 ] . Then the covariance matrix R~=E[x~(t)x~H (t)] can be obtained. By taking the eigendecomposition of R~ , we arrive at the following equation: [figure omitted; refer to PDF] where U~S is the signal subspace spanned by eigenvectors corresponding to the big eigenvalues and U~V is the noise subspace spanned by eigenvectors corresponding to the small eigenvalues; Σ~S and Σ~V are the diagonal matrixes associated with signal power and noise power, respectively. Since the signal subspace spanned by eigenvectors corresponding to the big eigenvalues is equal to the signal subspace spanned by array manifold, we have [figure omitted; refer to PDF]
As it is known, N×N unique nonsingular matrix T satisfies the following formula: [figure omitted; refer to PDF]
Obviously, (25) can be applied to the two subarrays equally, so we can get [figure omitted; refer to PDF] where US1 and US2 are signal subspaces of the two subarrays, respectively. Based on (26), the relationship between their signal subspaces can be expressed as [figure omitted; refer to PDF] where Ψ=T-1 ΦT . For clarity, we rewrite the formula above as [figure omitted; refer to PDF] where Γ~=Γ1Γ2-1 . Combining the expression of Γ1 and Γ2 , we can get [figure omitted; refer to PDF]
Thus η~ can be expressed as [figure omitted; refer to PDF]
As it is known that Ψ is the similarity matrix of N×N matrix Φ , it is easy to see that diagonal elements of Φ are eigenvalues of Ψ and every column of T is the eigenvector of the matrix Ψ . In this case, we only need to estimate the matrix Ψ , and then N eigenvalues e1 ,e1 ...,eN can be obtained with eigendecomposition. Finally, DOAs of incident signals can be estimated by [figure omitted; refer to PDF]
Considering the finite sampling data, the two signal subspace U^S1 and U^S2 are easily estimated. According to LS-ESPRIT algorithm and based on (27), we only need to minimize the following function to estimate the DOAs: [figure omitted; refer to PDF] where ||·||F is the Frobenius norm. Now, Γ~ and Ψ are still unknown. By letting Ψ be a variable and minimizing the function above, we can obtain the least-squares solution [figure omitted; refer to PDF]
Substituting formula (33) into (32), (32) can be expressed as [figure omitted; refer to PDF] where the (Mm -1)×(Mm -1) matrix K is [figure omitted; refer to PDF] and KH K=K . Substituting (29) into (34) and applying the conclusion proved in [34], we get the following equation: [figure omitted; refer to PDF] where [ecedil]9; is the Hadamard product. Here we assume that [figure omitted; refer to PDF]
In that way, ||Γ~U^S2 -U^S1 Ψ||F2 =η~HQd η~ . According to (30), it can be known that front the q-1 elements of η~ are defined as 1, so without loss of generality, the (M-1)×(q-1) matrix [straight epsilon]d can be used to constrain η~ : [figure omitted; refer to PDF] Thus solving the problem by minimizing the objective function (32) is transformed into [figure omitted; refer to PDF]
It should be noted that Qd is a singular matrix for infinite sample data. Actually, sample data is finite, which implies that Qd is a nonsingular matrix; thus the inverse matrix of Qd exists, and the inverse matrix of QdH also exists. This is the necessary condition for the Lagrange method, which means the Lagrange method can be applied to solve the problem in (39). In this case, we can derive the Lagrange function [figure omitted; refer to PDF] and define that [figure omitted; refer to PDF] where ∇ is differential operator. From (41), the equations can be expressed as follows: [figure omitted; refer to PDF] We rewrite these equations into the following matrix form: [figure omitted; refer to PDF] Then it can be seen that the coefficient matrix [QdH-[straight epsilon]d-[straight epsilon]dH 0] is the Lagrange matrix. Assuming that this matrix is the nonsingular matrix, the inverse matrix can be expressed as [figure omitted; refer to PDF]
Therefore, according to the definition of the inverse matrix, we can get [figure omitted; refer to PDF] The detailed expression can be written as [figure omitted; refer to PDF] Then, the expressions of matrix JH , VH , and S can be obtained by [figure omitted; refer to PDF] where Qd-H denotes the inverse matrix of QdH . In the following formulas, Qd-1 is the inverse matrix of Qd .
Based on (43), it can be inferred that [figure omitted; refer to PDF] Consequently, η~^ can be obtained by [figure omitted; refer to PDF]
Now, substituting (49) into (29) and (33), Ψ^ and its eigendecomposition can be obtained. According to (31), we can get the DOAs of signals θ^1 ,θ^2 ,...,θ^N .
Since η~^ has been calculated by the process above, based on the definition of (30), gain and phase coefficients can be obtained easily: [figure omitted; refer to PDF] Such a solution process is a joint estimation of DOA and gain and phase coefficients.
3.3. Mutual Coupling Coefficient Estimation
After obtaining DOAs of incident signals and gain and phase coefficients, we will discuss mutual coupling coefficients estimation. Here, the output of the whole array is used. For the whole array with all sensors Mw , the output data is x-(t) ; the corresponding covariance matrix is expressed as R- ; then the noise subspace can be obtained by the eigendecomposition of R- . In addition, we notice that the estimation of gain and phase coefficients has been carried out based on (51); hence, all elements of η- are well known. Due to the banded symmetric Toeplitz structure of MCM C- , the transformation matrix of MCM corresponding to each mutual coupling coefficient can be represented individually. In this case, the C- corresponding to the whole array can be expressed as [figure omitted; refer to PDF] Now, the array manifold taking the new form of C- can be expressed by [figure omitted; refer to PDF] where c-=[c0 ,c1 ,...,cp-1]T , W(θ) is M×p transform matrix, and each element is defined by [figure omitted; refer to PDF] where [figure omitted; refer to PDF] Based on the subspace theory, MCM can be obtained by minimizing fc , which is shown by [figure omitted; refer to PDF] By substituting (56) into (59), fc can be rewritten as [figure omitted; refer to PDF] where Qc =WH (θ)η-HU-^VU-^VH η-W(θ) . Please note Qc is totally independent of the mutual coupling coefficient vector c- .
Since sample data is finite, Qd is a nonsingular matrix with inverse matrix. Without loss of generality, the first element of c- is constrained to 1; namely, c-[straight epsilon] =1 and [straight epsilon] =[1,0,...,0]T . Thus the problem of minimizing fc is transformed into a quadratic programming problem with linear constraints, which is given by [figure omitted; refer to PDF] Through the same method as that in Section 3.2, mutual coupling coefficients can be obtained by [figure omitted; refer to PDF]
Thus mutual coupling coefficients have been estimated. This mutual coupling coefficients estimation method is applicable to any form of antenna array; just the transform matrix W(θ) is different. From the estimation of DOAs, gain and phase coefficients, and mutual coupling coefficients we can see that DOAs as well as gain and phase coefficients can be jointly estimated independently of mutual coupling coefficients firstly; with the known DOAs as well as gain and phase coefficients, mutual coupling coefficients can be estimated. So the three estimations do not affect mutually.
3.4. Algorithm Steps
The estimation method of DOAs, sensor gain and phase coefficients, and mutual coupling coefficients can be summarized as follows.
(1) Assuming a ULA with M+q+2(p+1) elements, where q+2(p+1) elements are auxiliary sensors, we use output data vector x-^(t) of the whole array to compute the covariance R-^ [figure omitted; refer to PDF] and utilize the output data of the first subarray x1 (t) as well as the second subarray x2 (t) of the middle array to get the output vector x~^(t) according to (22). Then covariance R~^ can be obtained: [figure omitted; refer to PDF]
(2) Perform eigendecomposition of R~^ and obtain the signal subspaces U^S1 and U^S2 ; perform eigendecomposition of R-^ and obtain the noise subspaces U-^V .
(3) Based on the extending LS-ESPRIT algorithm of this paper and (53), η~^ can be calculated. Based on the relationship between η~^ and Γ~^ , Γ~^ can be estimated by (29). Substituting Γ~^ into (33), we can obtain Ψ^ . Then by performing eigendecomposition, according to (31), we can obtain the DOAs of signals θ^1 ,θ^2 ,...,θ^N .
(4) Gain and phase coefficients ρ^1 ,...,ρ^M and [straight phi]^1 ,...[straight phi]^M can be estimated by formula (50).
(5) Mutual coupling coefficients 1,c^1 ,...,c^p-1 can be estimated by utilizing (62).
3.5. Performance Analysis and Comparison
In this paper, we focus on the simple solution to tackle sensor gain and phase errors and mutual coupling errors. Based on derivation in previous sections, it can be seen that the proposed method is computationally efficient because no spectral search and iteration are required. The contention WF method in [24] needs iteration and circular process, and the method in [10, 17] makes use of spectral-RARE estimator and thus requires spectral search.
With detailed analysis, there are some necessary conditions worthy of attention. First, in order to make the LS-ESPRIT algorithm work, we must use auxiliary sensors, and the number of q must be no less than two. The more the number of q is, the better the proposed algorithm performs. On each side of the array, p-1 auxiliary sensors are added. The number of p-1 is decided by p nonzero complex coefficients of MCM.
4. Simulation Results
In this section, we present several simulation results to illustrate the performance of the proposed method. In the following examples, we consider a ULA of 14 sensors with half wavelength spacing, and 9 of them are auxiliary sensors, where q=5 , p=3 , and M=5 , namely, Mw =14 and Mm =10 . The mutual coupling vector is set as c-=[1,0.4499+j0.5362,0.2598-j0.1500]T , and the gain and phase vector of uncertain sensors is given by η =[0.8ejπ/5 ,1.25e-jπ/3 ,1.53e-jπ/5 ,0.75ejπ/4 ,1.36e-jπ/10]T . The background noise is assumed to be AWGN and the signal-to-noise ratio (SNR) of the input signals is defined as 10log10 (σs2 /σn2 ) . The following measures are used to evaluate the performance: RMSED (root mean squared error of DOA estimation), RMSEG and RMSEP (RMSE of gain and phase estimation), and RMSEMCR and RMSEMCI (RMSE of real part and imaginary part of mutual coupling coefficient estimation): [figure omitted; refer to PDF] where K is the number of Monte Carlo experiments and N is the number of signals θ^n,i ; [straight phi]^n,i and θn , [straight phi]n are the estimated value and true value by the n th signal and the i th experiment, respectively; G^mi , P^mi and Gm , Pm are the estimated value and true value of gain and phase by the m th sensor and the i th experiment, respectively. In addition, m=1,...,q ; here q=5 ; R^vi , I^vi and Rv , Iv are the estimated value and true value of real part and imaginary part by the v th mutual coupling coefficient and the i th experiment, respectively; v=1,...,p-1 ; here p-1=2 .
4.1. DOA Estimation
In the first example, we present numerical simulation results to evaluate the superior performance of the proposed method in DOA estimation, in comparison with the gain and phase calibration method in [17, 20] and the mutual coupling calibration method in [10]. These methods are chosen because they are based on a similar data model where the auxiliary sensors are required, but these methods cannot calibrate gain and phase errors and mutual coupling errors simultaneously. Moreover, the results of MUSIC using the whole array with known uncertainties are also obtained. Assuming two signals from -10° and 10°, RMSEs of the DOA estimation are obtained by the above methods at different SNRs. 100 Monte Carlo experiments are performed at each SNR, and in each experiment the snapshot number is 200. The experiment results are shown in Figure 2. Figure 3 demonstrates the RMSE of the DOA estimation versus snapshot number with 100 Monte Carlo experiments, and the SNR is 10 dB.
Figure 2: RMSE of two signals DOA estimation versus SNR.
[figure omitted; refer to PDF]
Figure 3: RMSE of two signals DOA estimation versus snapshot number.
[figure omitted; refer to PDF]
According to Figures 2 and 3, we can see that all methods work better as SNR and snapshot number increase. Generally speaking, the estimation performance will stabilize when SNR and snapshot number reach a certain threshold. Moreover, it can be observed that the performance of the proposed method is close to the MUSIC algorithm with known uncertainties, which is much better than other methods. This is mainly because the proposed method can calibrate gain and phase errors as well as mutual coupling errors simultaneously while other methods can only calibrate one error, which seriously affects the performance of direction finding.
Next, performance comparison is made between the proposed method and the WF method which can calibrate these two types of errors at the same time through several experiments. The experiment results of MUSIC using the whole array with known uncertainties are also obtained, and CRB is displayed. Assuming three signals from -10°, 10°, and 20°, Figure 4 shows RMSE of DOA estimation versus SNR; 100 Monte Carlo experiments are performed at each SNR, and the snapshot number is 200; Figure 5 displays RMSE of the DOA estimation versus snapshot number with 100 Monte Carlo experiments, and the SNR is 10 dB.
Figure 4: RMSE of three signals DOA estimation versus SNR.
[figure omitted; refer to PDF]
Figure 5: RMSE of three signals DOA estimation versus snapshot number.
[figure omitted; refer to PDF]
According to Figures 4 and 5, we can see that both methods work better as SNR and snapshot number increase, and the performance of the proposed method is better than that of the WF method. Moreover, the satisfactory performance of the proposed method can be achieved even at low SNRs and few snapshots. When SNR is less than 10 dB, the performance of WF method deteriorates seriously while the proposed method is still effective. As snapshot number increases, the performance of WF method does not improve significantly. One possible explanation is that WF method is limited to the case where array errors are small, and its performance will be poor at low SNR. We can also note that the RMSE of DOA estimated by all methods except MUSIC with known uncertainties can not reach the CRB even at large SNRs and snapshot number.
4.2. Estimation of Gain and Phase Coefficients
Experimental results above have proved good performance of the proposed method for DOA estimation. The proposed method also can be used to estimate gain and phase coefficients. Assuming three signals from -10°, 10°, and 20° the snapshot number is 200. In Figures 6 and 7, RMSEG and RMSEP at different SNRs from 100 Monte Carlo experiments are demonstrated, where Gm , Pm , m=1,...,5 represent gain and phase of each uncertain sensor. Table 1 shows mean and variance of gain and phase coefficients for all uncertain sensors at SNR = 0 dB and SNR = 20 dB.
Table 1: Mean and variance of gain and phase coefficients for all uncertain sensors at SNR = 0 dB and SNR = 20 dB. True values are G1 =0.80 , G2 =1.25 , G3 =1.53 , G4 =0.75 , and G5 =1.36 ; P1 =0.6283 , P2 =-1.0472 , P3 =-0.6283 , P4 =0.7854 , and P5 =-0.3142 (rad).
SNR | 0 dB | 20 dB | ||||||||
m | 1 | 2 | 3 | 4 | 5 | 1 | 2 | 3 | 4 | 5 |
G ^ - m | 0.8087 | 1.2690 | 1.5651 | 0.7702 | 1.3699 | 0.7998 | 1.2504 | 1.5308 | 0.7504 | 1.3600 |
σ G ^ m | 8.2329 e - 4 | 3.8839 e - 3 | 1.0102 e - 2 | 3.2510 e - 3 | 1.2677 e - 2 | 8.0452 e - 5 | 1.5937 e - 4 | 3.0094 e - 4 | 5.0529 e - 4 | 4.8451 e - 4 |
P ^ - m | 0.6316 | -1.0510 | -0.6257 | 0.7878 | -0.3139 | 0.6284 | -1.0470 | -0.6282 | 0.7854 | -0.3137 |
σ P ^ m | 7.2978 e - 5 | 3.2913 e - 5 | 8.4479 e - 5 | 2.8760 e - 5 | 3.0745 e - 5 | 9.0441 e - 6 | 1.5761 e - 5 | 2.4466 e - 5 | 4.0881 e - 5 | 4.2990 e - 6 |
Figure 6: RMSE of gain coefficient estimation versus SNR.
[figure omitted; refer to PDF]
Figure 7: RMSE of phase coefficient estimation versus SNR.
[figure omitted; refer to PDF]
The experimental results show that the proposed method has a better performance in estimation of gain and phase coefficients as SNR increases. When SNR is higher than 5 dB, RMSEs of gain and phase can be less than 10%, which proves that when SNR is high enough, the coefficients could be estimated, with higher accuracy and stability.
4.3. Estimation of Mutual Coupling Coefficients
In the following experiments, the statistical efficiency of the proposed method in mutual coupling coefficient estimation is studied. Assuming three signals from -10°, 10°, and 20° the snapshot number is 200. Figure 8 shows RMSE of real and imaginary part of mutual coupling coefficients RMSEMCR and RMSEMCI versus SNR. 100 Monte Carlo experiments are performed at each SNR. Rv , Iv , v=1,2 represent the real and imaginary part of each mutual coupling coefficient. Table 2 shows mean and variance of the real and imaginary part of all mutual coupling coefficients at different SNRs.
Table 2: Mean and variance of the real and imaginary part of all mutual coupling coefficients at different SNRs. True values are R1 =0.4499 , R2 =0.2598 , I1 =0.5362 , and I2 =-0.1500 .
SNR | -5 dB | 5 dB | 15 dB | 20 dB | ||||
v | 1 | 2 | 1 | 2 | 1 | 2 | 1 | 2 |
R ^ - v | 0.1387 | 0.2266 | 0.4189 | 0.2599 | 0.4459 | 0.2603 | 0.4497 | 0.2598 |
σ R ^ - v | 7.4330 e - 3 | 2.7175 e - 3 | 2.2405 e - 3 | 4.7373 e - 4 | 5.8136 e - 4 | 5.4347 e - 5 | 1.6664 e - 5 | 3.8597 e - 6 |
I ^ - v | 0.3669 | -0.1058 | 0.5185 | -0.1432 | 0.5344 | -0.1492 | 0.5359 | -0.1500 |
σ I ^ - v | 8.6657 e - 3 | 1.2841 e - 3 | 1.1336 e - 3 | 1.3086 e - 4 | 7.3215 e - 5 | 1.6873 e - 5 | 1.0036 e - 5 | 1.8215 e - 6 |
Figure 8: RMSE of mutual coupling coefficients versus SNR.
[figure omitted; refer to PDF]
Figure 8 and Table 2 clearly illustrate that the proposed method has a better performance in estimation of mutual coupling coefficients as SNR increases. It is statistically efficient when SNR is higher than 5 dB, and RMSEs are all less than 10%.
In summary, simulation results have proved the effectiveness and excellent performance of the proposed method. As discussed in Section 4.1, better DOA estimation performance can be achieved when gain and phase errors and mutual coupling errors are simultaneously calibrated. Compared with the WF method, the satisfactory performance of the proposed method can be achieved even at low SNRs and few snapshots. As discussed in Sections 4.2 and 4.3, the proposed method can be used to accurately estimate gain and phase coefficients and mutual coupling coefficients, with even higher accuracy when SNR reaches a certain threshold.
5. Conclusion
In this paper, a simple and efficient method is proposed for direction finding in the presence of gain and phase errors as well as mutual coupling errors. By applying a group of auxiliary sensors, DOAs as well as gain and phase coefficients can be simultaneously estimated, and mutual coupling coefficients can also be estimated by utilizing the proposed decoupling method. Iterative operation, calibration sources, and spectral peak searching are not required for this method. However, auxiliary sensors are necessary to implement this algorithm, whose number affects algorithm performance. Different from other methods, the proposed method can achieve good performance with only a few auxiliary sensors. As shown in simulation results, the effectiveness and excellent performance of the proposed method are compared with other popular methods. DOA estimation can achieve good performance even at low SNRs and few snapshots. Moreover, gain and phase coefficients and mutual coupling coefficients can be accurately estimated when SNR reaches a certain threshold.
Acknowledgments
The authors thank the reviewers and editors for their corresponding contributions in making the paper more presentable. This work was supported in part by the National Science Foundation of China under Grant 61201410 and in part by the Fundamental Research Funds for the Central Universities under Grant no. HEUCF130804.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
[1] H. Krim, M. Viberg, "Two decades of array signal processing research: the parametric approach," IEEE Signal Processing Magazine , vol. 13, no. 4, pp. 67-94, 1996.
[2] D. H. Johnson, D. E. Dudgeon Array Signal Processing-Concepts and Techniques , Prentice-Hall, Englewood Cliffs, NJ, USA, 1993.
[3] P. Stocia, R. Moses Introduction to Spectral Analysis , Prentice-Hall, Upper Saddle River, NJ, USA, 2004., 2nd.
[4] S. Yin, S. Ding, A. Haghani, H. Hao, "Data-driven monitoring for stochastic systems and its application on batch process," International Journal of Systems Science , vol. 44, no. 7, pp. 1366-11376, 2013.
[5] S. Yin, H. Luo, S. Ding, "Real-time implementation of fault-tolerant control systems with performance optimization," IEEE Transactions on Industrial Electronics , vol. 64, no. 5, pp. 2402-2411, 2014.
[6] R. O. Schmidt, "Multiple emitter location and signal parameter estimation," IEEE Transactions on Antennas and Propagation , vol. 34, no. 3, pp. 276-280, 1986.
[7] R. Roy, A. Paulraj, T. Kailath, "ESPRIT--a subspace rotation approach to estimation of parameters of cisoids in noise," IEEE Transactions on Acoustics, Speech, and Signal Processing , vol. 34, no. 5, pp. 1340-1342, 1986.
[8] A. J. Weiss, B. Friedlander, "Effects of modeling errors on the resolution threshold of the music algorithm," IEEE Transactions on Signal Processing , vol. 42, no. 6, pp. 1519-1526, 1994.
[9] Y. Li, M. H. Er, "Theoretical analyses of gain and phase error calibration with optimal implementation for linear equispaced array," IEEE Transactions on Signal Processing , vol. 54, no. 2, pp. 712-723, 2006.
[10] K. C. Ho, L. Yang, "On the use of a calibration emitter for source localization in the presence of sensor position uncertainty," IEEE Transactions on Signal Processing , vol. 56, no. 12, pp. 5758-5772, 2008.
[11] B. C. Ng, C. M. Samson, "Sensor-array calibration using a maximum-likelihood approach," IEEE Transactions on Antennas and Propagation , vol. 44, no. 6, pp. 827-835, 1996.
[12] F. Sellone, A. Serra, "A novel online mutual coupling compensation algorithm for uniform and linear arrays," IEEE Transactions on Signal Processing , vol. 55, no. 2, pp. 560-573, 2007.
[13] C. Qi, Z. Chen, Y. Zhang, Y. Wang, "DOA estimation and self-calibration algorithm for multiple subarrays in the presence of mutual coupling," IEE Proceedings--Radar, Sonar and Navigation , vol. 153, no. 4, pp. 333-337, 2006.
[14] L. Bin, G. C. Zhi, C. Shing, "DOA estimation and tracking of ULAs with mutual coupling," IEEE Transactions on Aerospace and Electronic Systems , vol. 48, no. 1, pp. 891-905, 2012.
[15] Z. Ye, C. Liu, "On the resiliency of MUSIC direction finding against antenna sensor coupling," IEEE Transactions on Antennas and Propagation , vol. 56, no. 2, pp. 371-380, 2008.
[16] Z. Ye, C. Liu, "2-D DOA estimation in the presence of mutual coupling," IEEE Transactions on Antennas and Propagation , vol. 56, no. 10, pp. 3150-3158, 2008.
[17] L. Hong-Bing, G. Yi-Duo, G. Jian, J. Jun, "Mutual coupling self-calibration algorithm for uniform linear array based on ESPRIT," in Proceedings of the 2nd International Conference on Consumer Electronics, Communications and Networks (CECNet '12), pp. 3323-3326, Yichang, China, 2012.
[18] Z. Ye, J. Dai, X. Xu, X. Wu, "DOA estimation for uniform linear array with mutual coupling," IEEE Transactions on Aerospace and Electronic Systems , vol. 45, no. 1, pp. 280-288, 2009.
[19] C. Soon, L. Tong, Y. F. Huang, R. Liu, "Subspace method for estimating sensor gains and phases," IEEE Transactions on Signal Processing , vol. 42, no. 4, pp. 973-976, 1994.
[20] A. Liu, G. Liao, C. Zeng, Z. Yang, Q. Xu, "An eigenstructure method for estimating doa and sensor gain-phase errors," IEEE Transactions on Signal Processing , vol. 59, no. 12, pp. 5944-5956, 2011.
[21] C. Shenghong, Y. Zhongfu, X. Dongyang, X. Xu, "A Hadamard product based method for DOA estimation and gain-phase error calibration," IEEE Transactions on Aerospace and Electronic Systems , vol. 49, no. 2, pp. 1224-11233, 2013.
[22] M. Pesavento, A. B. Gershman, K. M. Wong, "Direction finding in partly calibrated sensor arrays composed of multiple subarrays," IEEE Transactions on Signal Processing , vol. 50, no. 9, pp. 2103-2115, 2002.
[23] C. M. S. See, A. B. Gershman, "Direction-of-arrival estimation in partly calibrated subarray-based sensor arrays," IEEE Transactions on Signal Processing , vol. 52, no. 2, pp. 329-338, 2004.
[24] S. A. Elkader, A. B. Gershman, K. M. Wong, "Rank reduction direction-of-arrival estimators with an improved robustness against subarray orientation errors," IEEE Transactions on Signal Processing , vol. 54, no. 5, pp. 1951-1955, 2006.
[25] L. Bin, C. Shing-Chow, "Direction finding with partly calibrated uniform linear arrays," IEEE Transactions on Antennas and Propagation , vol. 60, no. 2, pp. 922-929, 2012.
[26] C. M. S. See, "Sensor array calibration in the presence of mutual coupling and unknown sensor gains and phases," Electronics Letters , vol. 30, no. 5, pp. 373-374, 1994.
[27] K. V. Stavropoulos, A. Manikas, "Array calibration in the presence of unknown sensor characteristics and mutual coupling," in Proceedings of the European Signal Processing Conference (EUSIPCO '00), vol. 3, pp. 1417-1420, 2000.
[28] W. Ding, "Sensor array calibration in presence of mutual coupling and cain/phase errors by combining the spatial-domain and time-domain waveform information of the calibration sources," Circuits, Systems, and Signal Processing , vol. 32, no. 3, pp. 1257-1292, 2013.
[29] S. Yin, S. Ding, A. Haghani, H. Hao, P. Zhang, "A comparison study of basic datadriven fault diagnosis and process monitoring methods on the benchmark Tennessee Eastman process," Journal of Process Control , vol. 22, no. 9, pp. 1567-1581, 2012.
[30] B. Friedlander, A. J. Weiss, "Direction finding in the presence of mutual coupling," IEEE Transactions on Antennas and Propagation , vol. 39, no. 3, pp. 273-284, 1991.
[31] S. D. Solomon, D. A. Gray, Y. I. Abramovich, S. J. Andersoni, "Receiver array calibration using disparate sources," IEEE Transactions on Antennas and Propagation , vol. 47, no. 3, pp. 496-505, 1999.
[32] S. Yin, G. Wang, H. Karimi, "Data-driven design of robust fault detection system for wind turbines," Mechatronics , 2013.
[33] S. Yin, S. Ding, X. Xie, H. Luo, "A review on basic data-driven approaches for industrial process monitoring," IEEE Transactions on Industrial Electronics , no. 99, 2014.
[34] J. R. Magnus, H. Neudecker, "Mathematical preliminaries," Matrix Differential Calculus with Applications in Statistics and Econometrics , 1999., 2nd.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2014 Weijian Si et al. Weijian Si et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract
Many classical direction of arrival (DOA) estimation algorithms suffer from sensitivity to array errors. A simple but efficient method is presented for direction finding in the presence of gain and phase errors as well as mutual coupling errors. By applying a group of auxiliary sensors, DOAs and gain and phase coefficients can be simultaneously estimated, and mutual coupling coefficients can also be estimated by utilizing a novel decoupling method. The proposed algorithm does not require iterative operation or any calibration sources or spectral peak searching. Simulation results demonstrate the effectiveness of the proposed method.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer