An Efficient Algorithm for Direction Finding against Unknown Mutual Coupling
Abstract
: In this paper, an algorithm of direction finding is proposed in the presence of unknown mutual coupling. The preliminary direction of arrival (DOA) is estimated using the whole array for high resolution. Further refinement can then be conducted by estimating the angularly dependent coefficients (ADCs) with the subspace theory. The mutual coupling coefficients are finally determined by solving the least squares problem with all of the ADCs utilized without discarding any. Simulation results show that the proposed method can achieve better performance at a low signal-to-noise ratio (SNR) with a small-sized array and is more robust, compared with the similar processes employing the initial DOA estimation and further improvement iteratively.1. Introduction
Direction finding of multiple sources has received tremendous attention in the field of radar, sonar, mobile communication, and so on. High resolution algorithms, such as Multiple Signal Classification (MUSIC) [1] and Estimation of Signal Parameters via Rotational Invariance Techniques (ESPRIT) [2], can distinguish closely-spaced sources by the use of subspace theory and, therefore, have been widely used in the past few decades [3]. In spite of the potential advantages of the eigenstructure methods, their direct application to real systems is difficult, due to the critically required precise knowledge of the array manifold. In other words, the performance of the super-resolution techniques is strongly dependent on the array manifold accuracy. In practice, the array manifold is inevitably affected by mutual coupling, position perturbation and array gain or phase uncertainties. This results in significant distortion of the amplitude and phase of the signals received from the array. The direct employment of the eigenstructure-based methods leads to a serious degradation of direction finding [4–7]. Array calibration is an effective way to alleviate the deviation of direction of arrival (DOA) estimation.
Several calibration algorithms have been discussed in the last few decades. The earliest investigation was made by Schmidt [8] and Weiss [9]. They measured and stored the array manifold directly by interpolation. However, a large amount of memory is required, and the size and cost of the system could obviously be increased. In order to overcome the drawbacks of the above scheme, a kind of calibration method was proposed using a set of calibration sources in known locations. The maximum-likelihood approach [10] can be used to estimate the calibration matrix for array error compensation. Similarly, the algorithm in [11] was developed utilizing an iterative least mean squares approach. Although these techniques can obtain high accuracy and a large scope of calibration, it is impractical to set a collection of calibration sources and get their DOAs exactly as prior knowledge.
An approach to mitigate the influence of array errors is to calibrate the array by the use of the received signals. Such methods for estimating the DOA, unknown coupling, gain and phase simultaneously are called self-calibration. Friedlander and Weiss proposed to use an iterative process to acquire the parameters and DOA [12,13]. Svantesson formulized it as an optimization problem and solved the problem iteratively to estimate the mutual coupling coefficients for coupling compensation in the linear array of dipoles [14]. However, the result will converge to the local optimum if the initial values deviate far from the real ones. Alternatively, another kind of method was proposed conducting multidimensional search based on subspace fitting. The maximization of a posteriori (MAP) estimator proposed in [15] and [16] is one of these algorithms. Although it does not have the problem of convergence, non-linear multidimensional optimization in these methods is computationally consuming, and the convergence rate is relatively slow.
In recent years, several algorithms have been developed based on the characterization of the mutual coupling matrix (MCM). In [17], Ye obtained the initial estimation of DOA by setting the instrumental sensors on each side of the array. As a result, a one-dimensional search of the spatial spectrum can be performed directly using the original array data. Moreover, the result can be refined iteratively by estimating the mutual coupling coefficients for compensation. In [18], the signals received from the middle subarray were directly exploited for the traditional MUSIC, whereas the MCM is assumed to be a complex symmetric Toeplitz matrix. In order to increase the performance of DOA estimation, another method [19] proposed recently takes advantage of the special structure of MCM to parameterize the steering vector. It achieves the estimation of DOA using the whole array and improves the result by mutual coupling compensation.
In this paper, we assume that the MCM is a band-symmetric Toeplitz matrix with finite non-zero elements. By using the parameterized steering vector in [19], the preliminary DOA estimates can be obtained with the whole array. Further compensation of mutual coupling is then made by estimating the ADC first with the orthogonality of the subspace and the mutual coupling coefficients following with the full information of the ADCs. Simulation results show that the proposed algorithm has satisfactory performance compared with the methods in [17] and [19], especially when the SNR is low and small number of sensors exist.
The remainder of the paper is organized as follows. Section 2 is devoted to the problem formulation. Section 3 addresses the direction finding with initial estimation of DOA and coupling compensation by using the subspace theory and all of the parameters estimated. Section 4 gives the concluding remarks.
2. Problem Formulation
We consider a uniform linear array (ULA) consisting of M sensors with K narrowband far-field sources s1(t), s2(t),…, sK(t) received. The sources come from directions θ1, θ2, …, θK with respect to the normal line of the array. Assume that the distance between adjacent sensors is d and the wavelength of the carrier is λ. With the interactions between sensors, the mutual coupling effect cannot be ignored. The M × 1 output vector of the array is then given by:
Several studies of the coupling model [12,21] have shown that the coupling between a pair of sensors is nearly the same. Therefore, the MCM is a banded symmetric Toeplitz matrix for the ULA. Furthermore, based on the fact that the mutual coupling between two sensors is inversely proportional to their distance, the coefficient will be zero if two sensors are several wavelength apart. Let cij = cji = c|i−j| denote the mutual coupling coefficient between the i-th and the j-th element of the ULA, and assume that there are P distinct non-zero elements in the MCM with the self coupling c0 normalized as one. Then, the M × M matrix C can be expressed as:
If the MCM is known, the directions of sources θ1, θ2, …, θK can be estimated based on the spectrum function . However, the vector c is usually unknown, under which circumstance the traditional DOA estimation approach cannot be used, and a new method should be investigated for direction finding.
3. Direction Finding and Mutual Coupling Compensation
In this section, we first reformulate the equivalent steering vector by parameterizing the MCM. Then, we solve the DOA estimation problem in the presence of unknown mutual coupling by applying the whole array data for better performance [19,22]. For further refinement, a new mutual coupling compensation algorithm is proposed by utilizing the special structure of the reformulated MCM with no information lost.
3.1. DOA Estimation Using the Whole Array
Let rk represents the k-th element of the M × 1 equivalent steering vector am(θ). By combining Equations (2) and (3), am(θ) can be expressed as:
For notational clarity, we set , , and in the above equations. In the case of g0 ≠ 0, we can extract g0 out of rk and describe am(θ) as [19]:
Considering the subspace principle given in Equation (6) and am(θ) in Equation (11), we have the following equation on condition that g0 ≠ 0:
It is worth noting that Q(θ) does not contain any information of c. Therefore, the spectrum function given above can be employed even if the mutual coupling is unknown. Compared with the algorithm in [17], this spectrum function takes advantage of the whole array. People do not need to extract the middle array for mutual coupling eliminating. As a result, no information is lost in the spectrum estimation. This method is available only if:
3.2. Mutual Coupling Compensation
Once we get the DOA estimates, further refinement can be conducted based on the angularly dependent expression of am(θ) as shown in Equation (11). We assume θ̂ is the DOA estimated from Equation (14) and that no blind angles exists (i.e., g0 ≠ 0). Based on the fact that the noise subspace is orthogonal to the column subspace of am(θ̂), we have:
Notice that the P-th element of v(θ̂) is one, and the others are ADCs to be estimated. Denote , where qj, j = 1, …, 2P − 1 is the j-th column of Q(θ̂). By replacing the columns of Q(θ̂) and the associated elements of v(θ̂), Equation (16) can then be reformed as:
Consequently, we can get the estimates of μk and αk by solving the above equation as:
Now, we estimate the mutual coupling coefficients by utilizing the estimated ADCs μ̂k, α̂k. Based on the observation of μk, αk in Equations (9) and (10), it is not difficult to find that they are the linear functions of ck, k = 1, …, P − 1. Let B(θ̂) be the coefficient matrix between c′ = [c1, …, cP−1]T and [μ̂1,…, μ̂P−1, α̂1,…, α̂P−1]T. Then, we have:
From Equations (20) and (22), we can see that g0 must be determined before the estimation of Ck. Although containing the unknown mutual coefficient, g0 can be easily obtained by employing the particular composition of μ̂k, α̂k. Notice that μ̂1 and α̂P−1 have the complementary elements of g0. We can therefore get g0 by:
Combining Equations (23) and (22), the mutual coupling coefficient vector c′ can be determined by solving Equation (20) as:
For better performance, all of the DOAs estimated in the above subsection can be used to form the extended coefficient matrix B̃. Let B̃. = [BT(θ̂1),…, BT(θ̂K)]T and with B(θ̂i) and vgi as the matrix and vector evaluated at the i-th estimated DOA θ̂i. Then, the extension of Equation (20) will be:
Solving Equation (25) by the least squares, we can get a more precise estimation of c′ as:
The above approach provides us a means of mutual coupling compensation. That is to say, once the vector c′ is determined, the matrix C can be formed by locating its element on the corresponding sub-diagonal. Therefore, DOA estimation can be further obtained by searching the peak of:
The performance can be further improved by repeating the above procedure. The proposed algorithm for DOA estimation and mutual coupling compensation can be summarized as follows.
- (1)
Get L snapshots of the received signal x(t) at t = t1, …, tL, and form the following matrix as:
- (2)
Generate the covariance matrix using the above data matrix by:
- (3)
Conduct the EVD of R̂x, and get the noise subspace Ên.
- (4)
Scan the direction from −90° to 90° with 1° as the step size. Calculate the special spectrum using Equation (14), and obtain the initial estimation of DOA θ̂1, …, θ̂K.
- (5)
For each θ̂i, estimate the ADCs μ̂k, α̂k for k = 1, …, P − 1 based on Equation (19). Then, calculate Equations (23) and (22) to obtain the values of g0 and vg, respectively.
- (6)
Form the matrix B̃ and ṽg by B(θ̃i) and ṽgi, i = 1, …, K, and solve Equation (26) to get the parameters in the coupling matrix C.
- (7)
Enhance the DOA estimation with the estimated C̃ and Equation (27).
- (8)
Repeat Step (5) to Step (7) to get a more precise estimation of directions.
4. Simulation Results
In this section, simulations will be conducted to validate the performance of the proposed method. Consider two independent sources from the far-field incident on the ULA from θ1 = −10° and θ2 = 20°. Sensors are located in the array with equal spacing d = λ/2. The number of effective mutual coupling coefficients is P = 3. Here, we set c = [1, 0.43301 − 0.25i, 0.14142 − 0.14142i], which is used in the second simulation of [17], to guarantee g0 ≠ 0 at any direction θ.
In the first simulation, we evaluate the performance of DOA estimation in Step (4) without mutual coupling compensation. Assume that the array number is M = 7 and that the snapshot number is 500. The root mean squared error (RMSE) is used to compare the DOA accuracy of different algorithms. It can be calculated, in general, by:
The RMSE as a function of SNR is illustrated in Figure 1a. The method used in our proposed algorithm is superior at low SNRs compared with the method in [17], since the whole array is utilized. As the SNR increases, the RMSE of DOA estimation decreases gradually for all of the methods. When the SNR is greater than 10 dB, the accuracy of the two methods with unknown C is almost the same. Figure 1b shows the effect of the array size on RMSE when SNR = −5 dB. Notice that the choice of M should satisfy Equation (15). It can be seen that our method slightly outperforms the method in [17] when M ≤ 10.
Using the initial DOA estimates, we now proceed to get the mutual coupling vector c. In the second simulation, we consider the same scenario. Define the RMSE of c as , where Ĉk,n represents the estimated coupling coefficient in the n-th Monte Carlo experiment. From Figure 2a, we can see that the proposed method can obtain significant improvement of the coupling coefficient estimation when the SNR is lower than 3 dB. Besides, it is robust compared with the method in [17] and [19].
Now, we keep the SNR at −5 dB and vary the snapshots from 10 to 960. Figure 2b shows that the proposed method can achieve higher accuracy compared with the other two algorithms. From the estimation of the coupling matrix in Steps (5) and (6), it is not difficult to find that the full use of μ̂k, α̂k leads to the superior performance of the mutual coupling estimation.
In the third simulation, we investigate the DOA estimation obtained by mutual coupling compensation. We consider two sources incoming from θ1 and θ2 = θ1 + Δθ with SNR = 0 dB arriving at the ULA for M = 7 and L = 500. Define:
In the fourth simulation, we will investigate the effect of mutual coupling on the DOA estimation. First, consider a signal impinging on the ULA from −10°. Mutual coupling between the adjacent two sensors exists, i.e., P = 2. Then, define rsi = C(i, :)a(θ),i = 1, …, M as the response of the i-th element to the received signal. rs1 is presented in Figure 4a with the amplitude and phase of the coupling coefficient varying. The figure illustrates that when the phase stays near zero, the response becomes higher as the amplitude increases. When the phase grows into the range of [π/2, π], the result will be the inverse to that of [0, π/2). The effect of coupling is to some extent similar to the beamforming. Any change of the amplitude and phase could make the response different.
The RMSE as a function of the amplitude of c1 is presented in Figure 4b with the phase ϕ fixed at π/3 and 5π/6, respectively. From this, we can see that the results get better with the increase of |c1| on the condition that ϕ = π/3. The situation gets worse for ϕ = 5π/6, since the response becomes gradually lower. Figure 4c demonstrates the RMSE as a function of the phase with |c1| fixed at 0.15 and 0.35, respectively. The error increases slightly as the phase grows from zero to π, which coincides with the response shown in Figure 4a. Figure 4d presents the RMSE of the proposed method versus SNR with different coupling coefficients. In this experiment, we assume that there are two signals from −10° and 20° impinging on the ULA and set the four coefficients as c1 = [1, −0.2801 − 0.254i, −0.14 − 0.14i], c2 = [1, −0.125 + 0.108i, −0.066 + 0.858i], c3 = [1, 0.2801 + 0.254i, 0.14 + 0.14i] and c4 = [1, 0.125 + 0.108i, 0.066 + 0.858i]. From Figure 4d we can conclude that the proposed algorithm could achieve better performance when the coupling coefficient has a smaller phase and greater amplitude.
In the fifth simulation, we access the performance of DOA estimation when signals impinge on the ULA with different angle intervals. Figure 5a presents the response of every element to DOA with c = [1, 0.43301 − 0.25i, 0.14142 − 0.14142i]. Figure 5b illustrates the influence of the angle interval to the initial estimation in Step (4). It is shown that not only the angle interval, but also the initial DOA will affect the accuracy of estimation seriously. As long as the DOAs locate in the main lobe of the middle element, the gain of the array will not decline rapidly. Figure 5c gives the RMSE of DOA estimation versus SNR with different angle intervals. With the same initial DOA, the performance obviously decreases when the interval becomes bigger.
In summary, the proposed method can achieve more precise estimates of DOA as opposed to Ye's method [17] and Liao's method [19], especially when the SNR is low and the array size is close to the minimum available value of M = K+2P−1. The reservation of the information obtained guarantees the good performance of the proposed algorithm. The experiments of the mutual coupling effect show that the coupling coefficient, initial DOA and the angle interval could have an influence on the performance of DOA estimation at the same time.
5. Conclusions
This paper addresses the DOA estimation in the presence of unknown mutual coupling. Based on the subspace theory, the initial DOA is first estimated using the whole array without calibration sources and auxiliary sensors, which leads to high accuracy. With the assumption that no blind angles exist in the space, mutual coupling compensation is further conducted by estimating the coupling coefficients indirectly from the angular-dependent coefficients. Finally, with all of the ADCs utilized without discarding any, the mutual coupling coefficients are determined by solving the least squares problem. Simulations show that the proposed method can achieve better performance at low SNR with a small-sized array. The robustness of the method can be verified, as well.
Author Contributions
Weijiang Wang and Yingtao Ding provided the idea of this work. Weijiang Wang and Shiwei Ren conceived of and designed the experiments. Shiwei Ren performed the experiments and provided all of the figures and data for the paper. Haoyu Wang prepared the literature and analyzed the data. Weijiang Wang and Shiwei Ren wrote the paper. Correspondence and requests for the paper should be addressed to Yingtao Ding.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar]
- Roy, R.; Paulraj, A.; Kailath, T. ESPRIT–A subspace rotation approach to estimation of parameters of cisoids in noise. IEEE Trans. Acoust. Speech Signal Process. 1986, 34, 1340–1342. [Google Scholar]
- Krim, H.; Viberg, M. Two decades of array signal processing research: The parametric approach. IEEE Signal Process. Mag. 1996, 13, 67–94. [Google Scholar]
- Roller, C.; Wasylkiwskyj, W. Effects of mutual coupling on superresolution DF in linear arrays. Proceedings of the International Conference on Acoustics, Speech and Signal Processing (ICASSP), San Francisco, CA, USA, 23–26 March 1992; pp. E257–E260.
- Weiss, A.; Friedlander, B. Mutual coupling effects on phase-only direction finding. IEEE Trans. Antennas Propag. 1992, 40, 535–541. [Google Scholar]
- Weiss, A.; Friedlander, B. Effects of modeling errors on the resolution threshold of the MUSIC algorithm. IEEE Trans. Signal Process. 1994, 42, 1519–1526. [Google Scholar]
- Li, F.; Vaccaro, R. Sensitivity analysis of DOA estimation algorithms to sensor errors. IEEE Trans. Aerosp. Electron. Syst. 1992, 28, 708–717. [Google Scholar]
- Schmidt, R. Multilinear array manifold interpolation. IEEE Trans. Signal Process. 1992, 40, 857–866. [Google Scholar]
- Weiss, A.; Friedlander, B. Manifold interpolation for diversely polarized arrays. IEE Proc. Radar Sonar Navig. 1994, 141, 19–24. [Google Scholar]
- Ng, B.C.; See, C.M.S. Sensor-array calibration using a maximum-likelihood approach. IEEE Trans. Antennas Propag. 1996, 44, 827–835. [Google Scholar]
- Hung, E. Matrix-construction calibration method for antenna arrays. IEEE Trans. Aerosp. Electron. Syst. 2000, 36, 819–828. [Google Scholar]
- Friedlander, B.; Weiss, A. Direction finding in the presence of mutual coupling. IEEE Trans. Antennas Propag. 1991, 39, 273–284. [Google Scholar]
- Weiss, A.; Friedlander, B. Eigenstructure methods for direction finding with sensor gain and phase uncertainties. Circuits Syst. Signal Process. 1990, 9, 271–300. [Google Scholar]
- Svantesson, T. Modeling and estimation of mutual coupling in a uniform linear array of dipoles. Proceedings of the International Conference on Acoustics, Speech, and Signal Processing (ICASSP), Phoenix, AZ, USA, 15–19 March 1999; pp. 2961–2964.
- Wahlberg, B.; Ottersten, B.; Viberg, M. Robust signal parameter-estimation in the presence of array perturbations. Proceedings of the International Conference on Acoustics, Speech, and Signal Processing (ICASSP), Toronto, ON, Canada, 14–17 May 1991; pp. 3277–3280.
- Viberg, M.; Swindlehurst, A. A Bayesian approach to auto-calibration for parametric array signal processing. IEEE Trans. Signal Process. 1994, 42, 3495–3507. [Google Scholar]
- Ye, Z.; Liu, C. On the Resiliency of MUSIC Direction Finding against Antenna Sensor Coupling. IEEE Trans. Antennas Propag. 2008, 56, 371–380. [Google Scholar]
- Ye, Z.; Dai, J.; Xu, X.; Wu, X. DOA Estimation for Uniform Linear Array with Mutual Coupling. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 280–288. [Google Scholar]
- Liao, B.; Zhang, Z.G.; Chan, S.C. DOA Estimation and Tracking of ULAs with Mutual Coupling. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 891–905. [Google Scholar]
- Steyskal, H.; Herd, J. Mutual coupling compensation in small array antennas. IEEE Trans. Antennas Propag. 1990, 38, 1971–1975. [Google Scholar]
- Svantesson, T. Mutual coupling compensation using subspace fitting. Proceedings of the 1st IEEE Sensor Array and Multichannel Signal Processing Workshop (SAM 2000), Cambridge, MA, USA, 16–17 March 2000; pp. 494–498.
- Wang, B.; Wang, Y.; Chen, H.; Chen, X. Robust DOA estimation and array calibration in the presence of mutual coupling for uniform linear array. Sci. China Series Inf. Sci. 2004, 47, 348–361. [Google Scholar]
© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license ( http://creativecommons.org/licenses/by/4.0/).
Share and Cite
Wang, W.; Ren, S.; Ding, Y.; Wang, H. An Efficient Algorithm for Direction Finding against Unknown Mutual Coupling. Sensors 2014, 14, 20064-20077. https://doi.org/10.3390/s141120064
Wang W, Ren S, Ding Y, Wang H. An Efficient Algorithm for Direction Finding against Unknown Mutual Coupling. Sensors. 2014; 14(11):20064-20077. https://doi.org/10.3390/s141120064
Chicago/Turabian StyleWang, Weijiang, Shiwei Ren, Yingtao Ding, and Haoyu Wang. 2014. "An Efficient Algorithm for Direction Finding against Unknown Mutual Coupling" Sensors 14, no. 11: 20064-20077. https://doi.org/10.3390/s141120064
APA StyleWang, W., Ren, S., Ding, Y., & Wang, H. (2014). An Efficient Algorithm for Direction Finding against Unknown Mutual Coupling. Sensors, 14(11), 20064-20077. https://doi.org/10.3390/s141120064