Direction Finding with Gain / Phase Errors and Mutual Coupling Errors in the Presence of Auxiliary Sensors

Many classical direction of arrival (DOA) estimation algorithms suffer from sensitivity to array errors. A simple but efficientmethod 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, andmutual coupling coefficients can also be estimated by utilizing a novel decouplingmethod.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.


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][2][3][4][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][12][13][14][15][16][17][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][20][21][22][23][24][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 2 Mathematical Problems in Engineering 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][23][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][27][28][29] and self-calibration methods [30][31][32][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 leastsquares 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.
According to (1), the array covariance matrix of the array output is where R  = {s()s H ()} is signal covariance matrix, {⋅} represents the statistic expectation, and  2 and I denote noise variance and  ×  identity matrix, respectively.The eigendecomposition of R is expressed as where U  is signal subspace matrix containing eigenvectors with the largest eigenvalues and U  is noise subspace matrix containing eigenvectors with the smallest eigenvalues.Σ  and Σ  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 where  is snapshot number.After the eigendecomposition of array covariance matrix is accomplished, signal subspace matrix Û and noise subspace matrix Û can be estimated.
Considering mutual coupling errors and gain and phase uncertainties of the array, the corresponding array output vector is rewritten by 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  nonzero complex coefficients in MCM, with  < .The 1 ×  coefficient vector can be written as where   ,  = 1, . . .,  − 1 is the mutual coupling coefficient.The corresponding  ×  matrix C is given by Actually, mutual coupling errors and gain and phase errors coexist.Assuming that the gain and phase of  sensor are uncertain, gain and phase coefficients can be presented by  1 , . . .,   and  1 , . . .  , respectively. × 1 gain and phase vector is defined as T (8) and  ×  GPM can be written as Here, we are interested in the problem of how to estimate the DOAs {  }  =1 in the presence of sensors gain and phase errors and mutual coupling errors and how to obtain their coefficients, respectively.

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  sensors.First, we add  sensors in front of  sensors ULA and then add  − 1 auxiliary sensors on each side of the  +  sensors array.Thus a new array with a total of  +  + 2( − 1) sensors is formed, and the space between all neighboring sensors is .Here  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  affects algorithm performance.Besides, the array with a total of   =  +  + 2( − 1) sensors is called whole array and that with   =  +  middle sensors is called middle array.
Assuming that  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   − 1 sensors where  of them are well calibrated and the other comprising last   − 1 sensors where  − 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.
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 where 0   ×(−1) is an   × ( − 1) zero matrix and I   ×  is an   ×   unit matrix.Thus the MCM C of the whole array can be expressed as The MCM C of the middle array is written by Taking the unknown gain and phase uncertainties into account, their gain and phase vectors can be denoted by respectively, where 1 is a vector with all elements equal to one.
Here gain and phase matrix can be written as
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 where A 1 and A 2 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 The relationship between A 1 and A 2 is where Φ = diag{e 2 sin  1 / , e 2 sin  2 / , . . ., e 2 sin   / }.
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 where . Then the covariance matrix R = [x()x H ()] can be obtained.By taking the eigendecomposition of R, we arrive at the following equation: where Ũ is the signal subspace spanned by eigenvectors corresponding to the big eigenvalues and Ũ is the noise subspace spanned by eigenvectors corresponding to the small eigenvalues; Σ and Σ 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 As it is known,  ×  unique nonsingular matrix T satisfies the following formula: Obviously, (25) can be applied to the two subarrays equally, so we can get where U  1 and U  2 are signal subspaces of the two subarrays, respectively.Based on (26), the relationship between their signal subspaces can be expressed as where Ψ = T −1 ΦT.For clarity, we rewrite the formula above as where Γ = Γ 1 Γ −1 2 .Combining the expression of Γ 1 and Γ 2 , we can get Thus η can be expressed as As it is known that Ψ is the similarity matrix of  ×  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  eigenvalues  1 ,  1 . . .,   can be obtained with eigendecomposition.Finally, DOAs of incident signals can be estimated by Considering the finite sampling data, the two signal subspace Û 1 and Û 2 are easily estimated.According to LS-ESPRIT algorithm and based on (27), we only need to minimize the following function to estimate the DOAs: where ‖ ⋅ ‖  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 Substituting formula ( 33) into ( 32), ( 32) can be expressed as where the (  − 1) × (  − 1) matrix K is and K H K = K.Substituting ( 29) into (34) and applying the conclusion proved in [34], we get the following equation: where ⊙ is the Hadamard product.Here we assume that In that way, ‖ ΓÛ According to (30), it can be known that front the  − 1 elements of η are defined as 1, so without loss of generality, the ( − 1) × ( − 1) matrix   can be used to constrain η: Thus solving the problem by minimizing the objective function ( 32) is transformed into It should be noted that Q  is a singular matrix for infinite sample data.Actually, sample data is finite, which implies that Q  is a nonsingular matrix; thus the inverse matrix of Q  exists, and the inverse matrix of Q H  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 and define that where ∇ is differential operator.From (41), the equations can be expressed as follows: We rewrite these equations into the following matrix form: Then it can be seen that the coefficient matrix [ is the Lagrange matrix.Assuming that this matrix is the nonsingular matrix, the inverse matrix can be expressed as Therefore, according to the definition of the inverse matrix, we can get The detailed expression can be written as Then, the expressions of matrix J H , V H , and S can be obtained by where Q −H  denotes the inverse matrix of Q H  .In the following formulas, Q −1  is the inverse matrix of Q  .Based on (43), it can be inferred that Consequently, η can be obtained by 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 , . . ., θ .
Since η has been calculated by the process above, based on the definition of (30), gain and phase coefficients can be obtained easily: Such a solution process is a joint estimation of DOA and gain and phase coefficients.

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   , 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 Now, the array manifold taking the new form of C can be expressed by where c = [ 0 ,  1 , . . .,  −1 ] T , W() is  ×  transform matrix, and each element is defined by where Based on the subspace theory, MCM can be obtained by minimizing   , which is shown by By substituting (56) into (59),   can be rewritten as where Please note Q  is totally independent of the mutual coupling coefficient vector c.
Since sample data is finite, Q  is a nonsingular matrix with inverse matrix.Without loss of generality, the first element of c is constrained to 1; namely, c = 1 and  = [1, 0, . . ., 0] T .Thus the problem of minimizing   is transformed into a quadratic programming problem with linear constraints, which is given by Through the same method as that in Section 3.2, mutual coupling coefficients can be obtained by 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.

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 ++2(+1) elements, where  + 2( + 1) elements are auxiliary sensors, we use output data vector x() of the whole array to compute and utilize the output data of the first subarray x 1 () as well as the second subarray x 2 () of the middle array to get the output vector x() according to (22).Then covariance R can be obtained: (2) Perform eigendecomposition of R and obtain the signal subspaces Û 1 and Û 2 ; perform eigendecomposition of R and obtain the noise subspaces Û .

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  must be no less than two.The more the number of  is, the better the proposed algorithm performs.On each side of the array,  − 1 auxiliary sensors are added.The number of  − 1 is decided by  nonzero complex coefficients of MCM.

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  = 5,  = 3, and  = 5, namely,   = 14 and   = 10.The mutual coupling vector is set as c = [1, 0.4499 + 0.5362, 0.2598 − 0.1500]T , and the gain and phase vector of uncertain sensors is given by  = [0.8e/5 , 1.25e −/3 , 1.53e −/5 , 0.75e /4 , 1.36e −/10 ] T .The background noise is assumed to be AWGN and the signal-to-noise ratio (SNR) of the input signals is defined as 10log 10 ( 2  / 2  ).The following measures are used to evaluate the performance: RMSE D (root mean squared error of DOA estimation), RMSE  and RMSE  (RMSE of gain and phase estimation), and RMSE MCR and RMSE MCI (RMSE of real part and imaginary part of mutual coupling coefficient estimation): , Gain and phase calibration in [17] Gain and phase calibration in [20] Mutual coupling calibration in [10] The proposed method MUSIC with known uncertainties where  is the number of Monte Carlo experiments and  is the number of signals θ, ; φ, and   ,   are the estimated value and true value by the th signal and the th experiment, respectively; Ĝ , P and   ,   are the estimated value and true value of gain and phase by the th sensor and the th experiment, respectively.In addition,  = 1, . . ., ; here  = 5; RV , ÎV and  V ,  V are the estimated value and true value of real part and imaginary part by the Vth mutual coupling coefficient and the th experiment, respectively; V = 1, . . ., − 1; here  − 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 Gain and phase calibration in [17] Gain and phase calibration in [20] Mutual coupling calibration in [10] The proposed method MUSIC with known uncertainties The WF method in [24] The proposed method MUSIC with known uncertainties CRB  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.
The WF method in [24] The proposed method MUSIC with known uncertainties CRB 100 200 300 400 500 600 700 800 900 1000 Snapshots DOA RMSE (deg)  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.
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.

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, RMSE  and RMSE  at different SNRs from 100 Monte Carlo experiments are demonstrated, where   ,   ,  = 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.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.

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 RMSE MCR and RMSE MCI versus SNR. 100 Monte Carlo experiments are performed at each SNR. V ,  V , 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.
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.

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.

Figure 1 :
Figure 1: Array geometry.Solid points represent uncertain sensors and circles represent auxiliary sensors.

Figure 2 :
Figure 2: RMSE of two signals DOA estimation versus SNR.

Figure 3 :
Figure 3: RMSE of two signals DOA estimation versus snapshot number.

Figure 4 :
Figure 4: RMSE of three signals DOA estimation versus SNR.

Figure 5 :
Figure 5: RMSE of three signals DOA estimation versus snapshot number.
∘ , RMSEs of the DOA estimation are obtained by the above methods at different SNRs. 100 Monte Carlo experiments