About this Journal Submit a Manuscript Table of Contents
International Journal of Antennas and Propagation
Volume 2012 (2012), Article ID 712318, 13 pages
http://dx.doi.org/10.1155/2012/712318
Research Article

A Subspace-Based Compensation Method for the Mutual Coupling in Concentric Circular Ring Arrays for Near-Field Source Localisation

Centre for Health Technologies, School of Electrical, Mechanical and Mechatronic Systems, Faculty of Engineering and Information Technology, University of Technology, Sydney (UTS), P.O. Box 123, Broadway, NSW 2007, Australia

Received 14 June 2011; Accepted 5 September 2011

Academic Editor: Hoi Shun Lui

Copyright © 2012 Mohammed Jainul Abedin and Ananda Sanagavarapu Mohan. 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

We propose a technique for compensating the effect of mutual coupling on parameter estimation that is suitable with any subspace-based super-resolution algorithms. A Concentric circular ring array (CCRA) formed using thin dipole antennas in the receiving mode is employed to estimate the parameters of electromagnetic sources located in the radiating near field of the array. A CCRA geometry that obtains a lowest Cramer-Rao lower bound (CRLB) in the presence of array mutual coupling is chosen for investigation. The mutual coupling among antenna elements of the array would affect the orthogonality of subspaces when MUSIC or ESPRIT algorithms are used for parameter estimation. The proposed method obtains a compensation matrix that restores the orthogonality between the subspaces there by improving the accuracy of estimation. To avoid three-dimensional searches, the range parameter is estimated using a cross-correlation-based method. Numerical simulation using a full-wave electromagnetic (EM) solver is employed to demonstrate the effectiveness of the proposed compensation approach.

1. Introduction

A common practical issue in parameter estimation using array of antennas is the effect of mutual coupling. This effect has long been recognised as one of the leading causes of degradation of the performance of parameter estimation algorithms [1]. This degradation occurs due to the presence of mutual coupling between antenna elements which causes deviations to the array manifold. Many methods to combat the effect of mutual coupling on parameter estimation have been proposed in the literature [14]. The concept of coupling matrix for compensating the effect of mutual coupling was used in [5]; however, it did not properly account for the platform effects. Full-wave EM solvers were also employed by many authors. Adve and Sarkar [6] used the method of moments for computing the mutual coupling effect precisely on wire antenna arrays, and Rogier and Zutter [7] used full-wave EM techniques for computing the mutual coupling effects in a planar array. A minimum norm mutual coupling compensation for the application of direction of arrival (DOA) estimation was also reported in [8]. Yuan et al. [9] proposed a method for DOA estimation using MUSIC algorithm by considering universal steering vectors which does not require any additional mutual coupling compensation method. In addition to the above-mentioned methods, calibration techniques have also been proposed for compensating the effect of mutual coupling in DOA estimation algorithms. For instance, the calibration procedures developed by Weiss and Friedlander [10] require some perfectly calibrated sensors to compensate for the effect of mutual coupling as well as perturbations in gain and phase. An iterative least mean square approach was proposed by Hung [11] to estimate the calibration matrix, which requires some initial calibration. Self-calibration [12] and autocalibration [13] methods were also proposed which seek to iteratively minimise a function with respect to both direction of arrival (DOA) and mutual coupling. All the above-mentioned techniques mainly aim to compensate for the error due to lack of calibration, but all of them require a set of known calibration sources at some known locations. In practice, obtaining more than a single calibration source can be challenging for most of the localisation scenarios. Thus, compensation techniques that do not overly depend on multiple calibration sources are desirable. It should also be noted that all the above-mentioned works mainly aimed at either ULA or UCA geometries for compensating the mutual coupling. Mutual coupling matrix has a banded Toplitz structure for ULA and symmetric circulant structure with three bands for UCA [14], both of which can help to simplify the computation of the compensation matrix. However, such well-defined matrix structure may not be available for most of the complex arbitrary array geometries such as concentric circular ring arrays (CCRA) for which the effect of mutual coupling can be more severe. It would, therefore, be of interest to develop alternate mutual coupling compensation methods useful for complex planar arrays, which is the focus of this paper.

The effect of mutual coupling varies in arrays since the mutual impedances in a transmitting array can be different from that of a receiving array. Lui et al. [15, 16] have analysed this effect especially for the receiving mode of antenna arrays. They also propose a mutual coupling compensation technique based on a least square approach for the receiving mode of antenna arrays [15, 17]. However, their compensation technique can suffer from performance degradation when the number of emitting sources is less than the number of elements in the receiver antenna array. This is mainly because the compensation technique based on least squares method does not have a unique solution when either a single source or a smaller number of sources are present. Further, the investigations reported in [15, 17] on mutual coupling compensation are carried out only for ULA and UCA geometries for the estimation of azimuth angles of far-field sources. However, in many applications, estimation of the parameters of near-field sources when complex arbitrary antenna arrays are employed is of concern and particularly the techniques to reduce the effects of mutual coupling on the parameter estimation are of interest.

For near-field parameter estimation, the use of planar arrays can be more appropriate to overcome the limitations of ULAs, namely, the problems in resolving the noncoplanar emitters and the presence of angular ambiguities within the azimuthal plane [18]. However, the use of complex planar array geometries can pose some extra challenges in the form of the array steering vector imperfections due to mutual coupling which can lead to severe detrimental effects. Further, when parameters need to be estimated, the effect of mutual coupling varies on the type of estimation algorithms used. For subspace-based super-resolution estimation algorithms, the error due to mutual coupling can rotate the subspaces thereby disrupting the orthogonality between signal and noise subspaces. As a result, subspace-based super-resolution algorithms such as MUSIC and ESPRIT would lead to significant degradation of their estimation performance. These errors must be compensated in order to recover their estimation performance. Hence, there is a need for mutual coupling compensation methods to improve the estimation performance of MUSIC and ESPRIT algorithms for localising either single or multiple near-field sources when arbitrary planar antenna arrays are employed.

In this paper, we propose a mutual coupling compensation technique applicable for MUSIC and ESPRIT algorithms and investigate its performance when it is applied on a CCRA for the near-field parameter estimation. CCRA geometry can be formed by having rings that have either uniform or nonuniform radii. Similarly, the number of antenna elements over the circumference of each ring could also be variable. All these variations have an impact on the array mutual coupling of CCRA which in turn affects the parameter estimation. Thus, in order to choose an optimised CCRA geometry that can lead to lowest estimation bias; we derive CRLB for four different cases each with different ring radii and number of antenna elements so as to obtain a particular combination which provides the lowest root mean square error. The optimised array geometry is then employed for further processing, that is, the mutual coupling compensation and near-field parameter estimation using both MUSIC and ESPRIT algorithms. Since the parameter set required to be estimated for near-field localisation is larger than for the corresponding far-field case, it would be useful to reduce the computational load of conventional MUSIC and ESPRIT algorithms. To achieve this objective, here the range parameter is estimated using a cross-correlation-based method. Further, beamspace processing [19] is applied on CCRA so as to make it suitable for use with MUSIC and ESPRIT algorithms. We assume that the antennas receive correlated signals. Since the beamspace processing maps the array steering manifold of a planar array onto that of a virtual linear array, forward-backward smoothing can also be applied for decorrelation [18]. We employ a three-ring CCRA with thin wire dipole antennas as elements and obtain the mutual coupling matrix using a full-wave electromagnetic simulator. Simulation results are provided on the 3D localization of radiating sources placed in the near field of the CCRA by incorporating the proposed mutual coupling compensation method on MUSIC and ESPRIT algorithms.

2. Problem Formulation

2.1. Concentric Circular Ring Array (CCRA)

A CCRA is designed with 𝑃 concentric rings, and each ring is populated with the same number of 𝑀 omnidirectional dipole antennas on its circumference, as shown in Figure 1. The location of 𝑘th element on 𝑙th ring is denoted by its Cartesian coordinate {𝑥𝑘𝑙,𝑦𝑘𝑙}={𝜌𝑙cos𝜑𝑘𝑙,𝜌𝑙sin𝜑𝑘𝑙}, where 𝑘𝑙=11,21,,𝑀1,12,22,,𝑀2,,𝑀𝑃, with respect to the origin chosen at the centre of the circular ring structure, where 𝜌𝑙 is the radius of the 𝑙th ring and 𝜑𝑘𝑙 is the angular position of the 𝑘th element of the 𝑙th ring. For the uniform placement of array elements on the circumference of each ring, 𝜑𝑘𝑙=2𝜋(𝑘𝑙1)/𝑀𝑃. The radius of the 𝑙th ring is given by𝜌𝑙=𝜆min4sin𝜋/𝑀𝑃,(1) where 𝜆 is the wavelength. Now, the elements of the array steering vector for the 𝑙th ring can be obtained as𝑎𝑘𝑙=𝑒𝑗𝜔𝜏𝑘𝑙(𝑟,𝜃,𝜙),(2) where 𝜏𝑘𝑙 is the time delay of incoming signal coming from a source located at (𝑟,𝜃,𝜙) to 𝑘th element of 𝑙th ring, and time delay is determined as𝜏𝑘𝑙1(𝑟,𝜃,𝜙)=𝑐𝑟𝑘𝑙𝑟,(3) where𝑟𝑘𝑙𝑟=sqrt2+𝑥2𝑘𝑙+𝑦2𝑘𝑙2𝑥2sin𝜃𝑘𝑙cos𝜙+𝑦𝑘𝑙sin𝜙.(4) The array steering vector for the 𝑙th ring is given by𝐚𝑙=𝑒𝑗𝜔𝜏1𝑙(𝑟,𝜃,𝜙),𝑒𝑗𝜔𝜏2𝑙(𝑟,𝜃,𝜙),,𝑒𝑗𝜔𝜏𝑀𝑙(𝑟,𝜃,𝜙)𝑇,(5) where []𝑇 denotes transpose operation. Considering all the rings in the CCRA, the array steering vector can be written as𝐚𝐚(𝑟,𝜃,𝜙)=1(𝑟,𝜃,𝜙),𝐚2(𝑟,𝜃,𝜙),,𝐚𝑃(𝑟,𝜃,𝜙)𝑇.(6) A CCRA can be formed either by having same or varying number of elements on each ring [20]. Hence, the array steering vector of every ring may not be equal in length and thus can be represented as a group of separate single array steering vectors. The dimension of the total array steering vector is 𝑀𝑃×1, given by𝑒𝐚(𝑟,𝜃,𝜙)=𝑗𝜔𝜏11(𝑟,𝜃,𝜙),𝑒𝑗𝜔𝜏21(𝑟,𝜃,𝜙),,𝑒𝑗𝜔𝜏𝑀1(𝑟,𝜃,𝜙),𝑒𝑗𝜔𝜏12(𝑟,𝜃,𝜙),𝑒𝑗𝜔𝜏22(𝑟,𝜃,𝜙),,𝑒𝑗𝜔𝜏𝑀2(𝑟,𝜃,𝜙),𝑒𝑗𝜔𝜏1𝑙(𝑟,𝜃,𝜙),𝑒𝑗𝜔𝜏2𝑙(𝑟,𝜃,𝜙),,𝑒𝑗𝜔𝜏𝑀𝑃(𝑟,𝜃,𝜙)𝑇.(7)With an aim to apply beamspace transformation, we express the array steering vector of CCRA into a phase mode as𝐚𝑓=𝑓𝐻𝑝𝐚(𝑟,𝜃,𝜙).(8) This transformation maps the array steering vector of CCRA to that of an equivalent virtual linear array. Defining “𝐁” as the discrete Fourier transform matrix given by𝐁=1𝑒𝑗𝜑𝑒𝑗(𝑀1)𝜑1𝑒𝑗𝜑𝑒𝑗(𝑀1)𝜑,𝑓𝐁=𝐻𝑓𝐻,(9) where ()𝐻 denotes complex conjugate transpose and the elements of (9) can be expressed as𝑓𝑝=1𝑀1,𝑒𝑗𝑝𝜑,,𝑒𝑗𝑝(𝑀1)𝜑𝑇,(10) where 𝑝=,+1,,.

712318.fig.001
Figure 1: Concentric circular ring array (CCRA).
2.2. Selection of Optimal CCRA Configuration Based on CRLB

Here, our aim is to select an optimal CCRA configuration that provides lowest estimation bias using Cramer-Rao lower bound (CRLB) and employ that array for further processing in this paper. CRLB is derived in this section for CCRA in the presence of mutual coupling. Assuming the received signal vector 𝐬 to be complex valued and Gaussian distributed, the covariance matrix can be formed as 𝐑=𝐀𝐑𝑠𝐀𝐻+𝜎2𝑛𝐈,(11) where 𝐀=𝐂𝐀 and 𝐑𝑠=𝐸[𝐬𝐬𝐻]. The unknown parameters which need to be estimated can be written in a vector form 𝛼Θ=𝐫,𝜽,𝝓,𝑘1𝑀𝑃,𝛽𝑘1𝑀𝑃𝑇,(12) where 𝜽=[𝜃1,,𝜃𝑑]𝑇,𝝓=[𝜙1,,𝜙𝑑]𝑇,𝛼𝑘={𝐶re𝑘}, and 𝛽𝑘={𝐶im𝑘}. For convenience, we are using notations 𝛼𝑘 and 𝛽𝑘 instead of 𝛼𝑘𝑙 and 𝛽𝑘𝑙, for 𝑘=1,2,,𝑀𝑃. Considering 𝑁 snapshots of the received signal, the Fisher Information Matrix (FIM) can be formed as𝐹𝑚𝑛𝑅=𝑁𝑡𝑟1𝜕𝑅𝜕Θ𝑚𝑅1𝜕𝑅𝜕Θ𝑛.(13) The elements of FIM are derived [21, 22] and presented in the appendix. Referring to the appendix, the FIM which is a (3𝐷+𝑀𝑃)×(3𝐷+𝑀𝑃) matrix can be written as||||||||||||||𝐹FIM𝐅=𝜃𝜃𝐹𝜃𝜙𝐹𝜃𝑟𝐹𝜃𝜶𝐹𝜃𝜷𝐹𝑇𝜃𝜙𝐹𝜙𝜙𝐹𝜙𝑟𝐹𝜙𝜶𝐹𝜙𝜷𝐹𝑇𝜃𝑟𝐹𝑇𝑟𝜙𝐹𝑟𝑟𝐹𝑟𝜶𝐹𝑟𝜷𝐹𝑇𝜃𝜶𝐹𝑇𝜶𝜙𝐹𝑇𝑟𝜶𝐹𝜶𝜶𝐹𝜶𝜷𝐹𝑇𝜃𝜷𝐹𝑇𝜙𝜷𝐹𝑇𝑟𝜷𝐹𝜷𝜶𝐹𝜷𝜷||||||||||||||,(14) where each term represents a block matrix, for example, 𝐹𝜃𝜶=[𝐹𝜃𝛼𝑘]𝑘=1𝑀𝑃,𝐹𝜶𝜶=[𝐹𝛼𝑘𝛼𝑘]𝑘=1𝑀𝑃, and 𝐹𝜷𝜷=[𝐹𝛽𝑘𝛽𝑘]𝑘=1𝑀𝑃 and the other block matrices follow the same notational convention. CRLB is the inverse of FIM given by 𝐉=𝐅1, and for all the desired parameters it can be obtained asCRLBΘ1=sqrt2𝐷3𝐷𝑑=1𝐽𝑑𝑑,(15)CRLB𝑐1=sqrt𝑐23𝐷+𝑀𝑃𝑑=3𝐷+1𝐽𝑑𝑑×100%.(16) In (16), CRLB with mutual coupling is calculated as a percentage as it is a relative value [21]. The CRLB for CCRA is computed by varying either the number of antenna elements or the ring radii, so that we can chose an optimum array configuration that offers the lowest bias for estimating the near-field parameters in the presence of array mutual coupling. In this paper, we consider four different variations of CCRA for computing the CRLB and the results are indicated in Figures 2, 3, and 4. It can be observed, from Figures 24, that lowest CRLB results when the number of antenna elements on each ring is not the same and also when the incremental variation of ring radii is not uniform. Thus, the optimal CCRA configuration for performing the near-field parameter estimation is chosen.

712318.fig.002
Figure 2: CRLB of estimated azimuth for four different arrangements: (a), (b), (c), and (d) of CCRA.
712318.fig.003
Figure 3: CRLB of estimated elevation angle for four different arrangements: (a), (b), (c), and (d) of CCRA.
712318.fig.004
Figure 4: CRLB of estimated range for four different arrangements: (a), (b), (c), and (d) of CCRA.
2.3. Signal Modelling in the Presence of Mutual Coupling

The array response of an 𝑀-element CCRA due to 𝑑th near-field source located at (𝑟,𝜃,𝜙) is given by𝑥(𝑡)=𝐷𝑑=1𝐂𝐚𝑓𝑟𝑑,𝜃𝑑,𝜙𝑑𝑠𝑑(𝑡)+𝑛(𝑡),(17) where 𝐂 denotes a matrix that contains the mutual coupling among the array elements. The coupling is inversely proportional to the distance between the elements. After normalising the main diagonal elements to unity, the mutual coupling matrix is expressed as follows:𝑍𝐂=𝐴+𝑍𝑇𝐙+𝑍𝑇𝐈𝑁1,(18) where 𝑍𝐴 is the impedance of an isolated individual element, 𝑍𝑇 is the impedance of each receiver element, and 𝐈𝑀 is the identity matrix. The impedance 𝑍𝑇 is considered to be the complex conjugate of 𝑍𝐴 ensuring impedance matching for maximum power transfer. The elements of 𝐙 matrix are||||||||||||𝑍𝐙=𝐴+𝑍𝑇𝑍1121𝑍11𝑀𝑃𝑍2111𝑍𝐴+𝑍𝑇𝑍21𝑀𝑃𝑍𝑀𝑃11𝑍𝑀𝑃21𝑍𝐴+𝑍𝑇||||||||||||.(19) Now, 𝑍 is defined as𝑧𝑍=𝑘𝑙𝑘𝑙=𝑅𝑘𝑙𝑘𝑙+𝑗𝑋𝑘𝑙𝑘𝑙,(20) where 𝑅𝑘𝑙𝑘𝑙 and 𝑋𝑘𝑙𝑘𝑙 are the resistive and reactive components of the impedance 𝑍. Using a number of snapshots of the received signal in (17), the data covariance matrix can be expressed as𝐑=𝐷𝑑=1𝜎2𝑑𝐂𝐚𝑓𝑟𝑑,𝜃𝑑,𝜙𝑑𝐚𝐻𝑓𝑟𝑑,𝜃𝑑,𝜙𝑑𝐂𝐻+𝜎2𝑛𝐈.(21) In order to apply subspace-based algorithms such as MUSIC or ESPRIT for parameter estimation, an eigenvalue decomposition has to be performed which is given by𝐑=𝐔𝑠Λ𝑠𝐔𝐻𝑠+𝜎2𝑛𝐔𝑛𝐔𝐻𝑛,(22) where 𝐔𝑠 and 𝐔𝑛 contain the signal and noise eigenvectors, Λ𝑠 is the diagonal matrix that contains signal power, and 𝜎2𝑛 denotes noise power. In order to estimate parameters accurately, the signal and noise subspaces must be mutually orthogonal. But, in the presence of mutual coupling 𝐂, the orthogonality no longer exists between signal and noise subspaces. Hence, we need to obtain a matrix 𝐂 that can compensate the effects of mutual coupling and thereby restoring orthogonality between the subspaces. Thus, incorporation of 𝐂 makes 𝐔𝐻𝑛𝐂𝐚(𝑟𝑑,𝜃𝑑,𝜙𝑑)=0 for 𝑑=1,,𝐷 which will lead to proper estimation of the parameters using subspace-based algorithms. However, in practice, 𝐂 is not known. Hence, methods are required to accurately estimate 𝐂 prior to applying any subspace-based parameter estimation procedure when the mutual coupling is present among the elements of the array.

2.4. Subspace Rotation due to Array Mutual Coupling

For the CCRA, in an ideal case, that is, when the mutual coupling is ignored, the array steering vector due to an incident signal from 𝑑th source is given by (8). Under the effect of mutual coupling, the modified array steering vector ̃𝐚𝑓(𝑟,𝜃,𝜙) is given bỹ𝐚𝑓(𝑟,𝜃,𝜙)=𝐂𝐚𝑓(𝑟,𝜃,𝜙).(23) As noted earlier, the mutual coupling disrupts the orthogonality of subspaces since the signal subspace that is spanned by the steering vector is rotated from its true position. Therefore, it can be seen that𝐴span𝑓=span𝐂𝐴𝑓,(24) where 𝐴𝑓 is the array manifold matrix of size 𝑀𝑃×𝐷 that is created from the array steering vector (7) in the presence of 𝐷 sources. Referring to (17), the mutual coupling matrix 𝐂 for CCRA is given by|||||||||||||||||||||||𝑐𝐂=1111𝑐1121𝑐1131𝑐11𝑀1𝑐1112𝑐1122𝑐11𝑘𝑙𝑐11𝑀𝑃𝑐2111𝑐2121𝑐2131𝑐21𝑀1𝑐2112𝑐2122𝑐21𝑘𝑙𝑐21𝑀𝑃𝑐1211𝑐1𝑝11𝑐1𝑝𝑀1𝑐1𝑝𝑘𝑙𝑐𝑀𝑝11𝑐𝑀𝑝21𝑐𝑀𝑝𝑀1𝑐𝑘𝑃𝑘𝑙𝑐𝑀𝑝𝑀𝑝|||||||||||||||||||||||.(25)Now, the rotated subspace due to 𝐂 can be expressed as̃𝐮(𝑟,𝜃,𝜙)=𝐂𝐮(𝑟,𝜃,𝜙),(26) where ̃𝐮 are the elements of the rotated signal and noise subspaces due to the presence of mutual coupling. The effect of mutual coupling 𝐂 must be compensated so as to make subspaces orthogonal in order to accurately estimate parameters.

2.5. Subspace Rotation due to Simplified Mutual Coupling Model in a CCRA

Our aim here is to derive the rotated subspaces for CCRA under simplified assumptions to demonstrate our method. The mutual coupling in an array depends on the separation between two adjacent antenna elements, and hence nonadjacent array elements have weaker coupling. Realising this, a single coupling coefficient was considered for computing the compensation matrix in the literature [14, 23]. Also, the coupling effect was ignored for an ULA for interelement separation greater than 0.707𝜆 [21]. In this section, in order to derive the subspace rotation due to mutual coupling, we simplify (25) by assuming large interring separation for CCRA so as to make the interring coupling insignificant and thus can be ignored. Also, in this case, coupling among elements of outer ring can also be ignored due to larger interelement separation since the elements are distributed uniformly on rings with larger circumferences thus leading to larger interelement separation. Under these assumptions, only the mutual coupling among the adjacent elements of the innermost rings will be considered to simplify the problem. Thus, the simplified mutual coupling matrix 𝐂 can be obtained from (25) as 𝐂=1𝑐11𝑐21𝑐0111𝑐110001.(27) Now, simplifying (25) for the above case, the expression for the rotated subspace can be expressed as ̃𝐮(𝑟,𝜃,𝜙)=1+𝑐11𝑢,𝑐11+𝑎+𝑐21𝑢2,,𝑐21𝑢𝑀3+𝑢𝑀2+𝑐21𝑢𝑀1,𝑐21𝑢𝑀2+𝑢𝑀1𝑇.(28) After further manipulation, (28) can be represented as ̃𝐮(𝑟,𝜃,𝜙)=1+𝑐21𝑢1,𝑢,,𝑢𝑀2,𝑢𝑀1𝑇+𝑐210,1,,𝑢𝑀3,𝑢𝑀2𝑢𝑀𝑇=𝐮1𝑚𝑐(𝑟,𝜃,𝜙)+𝐮2𝑚𝑐(𝑟,𝜃,𝜙),(29) where 𝐮𝑚𝑐 denotes the modified eigenvector in the presence of mutual coupling. Thus, the expression for the estimated covariance matrix modified due to the presence of mutual coupling becomes 𝐑=𝐔𝟏𝐻𝑚𝑐𝑟𝑑,𝜃𝑑,𝜙𝑑Λs𝐔1𝑚𝑐𝑟𝑑,𝜃𝑑,𝜙𝑑+𝐔2𝐻𝑚𝑐𝑟𝑑,𝜃𝑑,𝜙𝑑Λ𝑠𝐔2𝑚𝑐𝑟𝑑,𝜃𝑑,𝜙𝑑+𝜎2𝑛𝐈.(30) Comparing (30) and (21), for the simplified case considered here, one can observe that the variation in covariance matrix is due to the contribution of the factor 1+𝑐21𝑢 which causes subspaces to rotate. As a result, it would lead to performance degradation of the subspace-based super-resolution algorithms.

3. Compensation for the Mutual Coupling

The effect of mutual coupling must be compensated for achieving a desired performance. This can be achieved by formulating a compensation matrix that minimises the error by which the array manifold gets deviated from that of actual one [14]. For subspace-based estimation algorithms, the deviation in array manifold due to mutual coupling does appear as an error that disrupts the orthogonality between signal and noise subspaces of the covariance matrix. Similar concept was also used in [24] where genetic algorithm (GA) was used to minimise the error to obtain a compensation matrix in order to recover the estimation performance of the MUSIC algorithm. However, the method proposed in [24] requires MUSIC pseudospectrum for a known calibration source to correlate with the pseudospectra of unknown sources at each search point in the space which can be computationally demanding. Alternatively, in this paper, we propose a method that uses an ideal array response without mutual coupling to estimate a compensation matrix by minimising the error between the rotated and ideal subspaces, so that the introduction of the compensation matrix restores the orthogonality between the subspaces. The advantage of the method proposed here as compared to that used in [24] is that, in the present case, the compensation matrix can be determined with a lower computational load as it does not require performing correlation of MUSIC pseudospectra at every search point within the parameter space. As a result, the proposed compensation method achieves quicker convergence. Also, the proposed method is applicable for any subspace-based parameter estimation algorithm, as will be demonstrated below where we also provide a detailed description of the proposed technique.

3.1. Role of the Compensation Matrix

In order to remove the effect of mutual coupling completely, one needs to estimate the compensation matrix as 𝐂=𝐂1. However, in actual practice, the array output model cannot account for all the sources of signal impairments in an antenna array, hence it is difficult to achieve 𝐂 equal to 𝐂1. Considering (17), even in a noise-free situation, that is, when 𝐧(𝑡) becomes a null vector, 𝐂1 can only be considered to be proportional to 𝐂. Therefore, only an estimated compensation matrix can be obtained. Assuming that ideal array response from reference sources is known, the estimation procedure for 𝐂 can be expressed as a minimisation of the function which is given by𝐂𝐚min𝐻𝑘𝑙(𝑟,𝜃,𝜙)𝐂𝐑𝑥𝐂𝐻𝐚𝑘𝑙(𝑟,𝜃,𝜙)2𝐹,(31) where 𝐚𝑘𝑙(𝑟,𝜃,𝜙) is the orthogonal null steering vector that spans the equivalent null space of the reference source. Incorporating 𝐂 obtained from (31) into (26), the subspace 𝐔 will approximately be equal to 𝐔 since 𝐂𝐂𝐈. Also, in (28), all coupling terms will almost vanish once the compensation terms are incorporated and thereby the subspace rotation due to the effect of mutual coupling will almost be nullified resulting in the effective restoration of orthogonality between subspaces.

Comparing the proposed technique with the other published compensation techniques in the literature shows that the method used in [11, 15] requires at least 𝑀𝑃 different near-field calibration sources around CCRA to achieve 𝐚𝑘𝑙(𝑟,𝜃,𝜙) of (31), where 𝑘𝑙=11,21,,𝑀1,12,22,,𝑀2,,𝑀𝑃. This requires estimation of a compensation matrix 𝐂 of the size 𝑀𝑃×𝑀𝑃 by minimising (31) for all 𝑀𝑃 calibration sources. For the same case, the least square approach would require at least more than 𝑀𝑃 calibration sources to achieve a unique solution for the estimation of compensation matrix 𝐂. Therefore, methods that do not require multiple calibration sources for estimating the compensation matrix are desirable.

3.2. Proposed Compensation Method

Here, it is hypothesized that the error that destroys the orthogonality between signal and noise subspaces is only due to the mutual coupling and all other possible sources of impairment are ignored. As a result, 𝐔𝐻𝑛𝐴(𝑟,𝜃,𝜙)0, which is due to the presence of mutual coupling among the elements of the receiver antenna array. Here, 𝐔𝑛 is the noise subspace that is obtained from EVD of covariance matrix formed by the measured data at the receiver. Since 𝐔𝑛 is corrupted by the array mutual coupling, we need to estimate 𝐂. For that, instead of minimising the functional given in (31), we reformulate an objective function 𝑄 as𝐂𝐔𝑄=min𝐻𝑛𝐴(𝑟,𝜃,𝜙)𝐔𝐻𝑛𝐂𝐴(𝑟,𝜃,𝜙)2𝐹.(32) Now, the estimated compensation matrix 𝐂 can be computed by performing an optimisation method to minimise (32). Here, the true noise subspace 𝐔𝑛 is obtained from a single known near-field calibration source for an ideal array without mutual coupling with known array geometry and element positions and assuming no gain and phase deviation to occur. We employ genetic algorithm to determine 𝐂 by minimising 𝑄 of (32). Incorporating 𝐂 restores the orthogonality between subspaces, thereby recovering the performance of parameter estimation even in the presence of array mutual coupling.

4. Parameter Estimation after Compensation

In this section, we show as to how the proposed technique of compensation is incorporated into both MUSIC and ESPRIT algorithms for estimating the parameters of near-field sources. In order to reduce the computational burden in MUSIC and ESPRIT, we estimate range parameter by using a method that employs cross-correlation among the eigenvectors of the signal subspace [25]. A brief description of the procedure for parameter estimation of near-field sources is given below.

4.1. MUSIC Algorithm for Estimating Azimuth and Elevation Angles

The azimuth and elevation angles of near-field sources can be obtained from the two-dimensional MUSIC pseudospectrum given by1𝑃=𝐔𝑛𝐂𝐚𝑓(𝜃,𝜙)2,(33) where 𝐂 is the estimated compensation matrix obtained by minimising 𝑄 in (32).

4.2. ESPRIT Algorithm for Estimating Azimuth and Elevation Angles

For ESPRIT algorithm, the signal subspace matrices can be written as [26]𝚪𝐔=𝐔1Ψ+𝐃𝐈𝐔1Ψ,(34) where Ψ=𝑇1Φ𝑇, 𝐃=diag{(1)𝑀𝑃2,,(1)1,(1)0,(1)1,,(1)𝑀𝑃} and 𝐔𝐂𝐔= is the matrix of eigenvectors that is obtained after mutual coupling compensation. This system of equations has a unique solution when the number of sources is less than the number of antenna elements in CCRA. The solution can be obtained as Φ=𝑇Ψ𝑇1, and the eigenvalues of Ψ provide the diagonal elements of Φ which therefore yield automatically paired source azimuth and elevation angles as𝜙𝑑𝑢=arg𝑑+𝑗𝑣𝑑,𝜃𝑑=sin1𝑢sqrt2𝑑+𝑣2𝑑.(35)

4.3. Cross-Correlation among Signal Eigenvectors for Time-Delay Estimation

Since the eigenvectors contain the time delay information, the time delay among the signal paths can be estimated via the cross-correlation of eigenvectors. The compensated signal eigenvectors are 𝐮𝑖 and 𝐮𝑗, where 𝐮𝑖=𝐂̃𝐮𝑖 and 𝐮𝑗=𝐂̃𝐮𝑗. A normalised cross-correlation of 𝐮𝑖 and 𝐮𝑗 is given by𝑅𝑖,𝑗𝑟(𝜏)=𝑖,𝑗(𝜏)𝑟𝑖,𝑖(0)𝑟𝑗,𝑗(0),(36) where the 𝑟𝑖,𝑗(𝜏) is the cross-correlation of 𝐮𝑖(𝑘𝑙) and 𝐮𝑗(𝑘𝑙) and it can be expressed as𝑟𝑖,𝑗1(𝜏)=𝑁𝑠𝜏𝑀𝑃𝜏𝑘𝑙=1𝑢𝑖𝑘𝑙𝑢+𝜏𝑗𝑘𝑙,(37) where 𝜏=𝑀𝑃+1,,𝑀𝑃1 and 𝑢𝑖(𝑛) is the 𝑛th element of 𝐮. An estimated time delay can be obtained by maximising the 𝑅𝑖,𝑗(𝜏). Using the time delay, the range parameter is estimated.

5. Simulation Results and Discussion

In order to demonstrate the performance of the proposed compensation method for near-field source localisation, we consider a CCRA with 27 half-wave, thin, dipole antennas which are distributed uniformly along the circumferences of three concentric rings as shown in Figure 1. Antenna element distribution on each ring and the radii of rings are given in Table 1. The radius of each half-wave dipole is considered to be equal to 5×103𝜆 at an operating frequency of 2.4 GHz and a 50 Ω load is assumed to be connected to the terminal of each of the dipole antenna. The sources of EM radiation to be localised are modelled using three vertically polarised thin dipole antennas which are located at (2.2𝜆,120,40), (3𝜆,70,70), and (4𝜆,80,20) which lie within the radiating near field of the receiver antenna array. Each of the source dipole is excited with a voltage source. Induced current at every dipole element of CCRA due to the incident field emanated from three dipole sources is computed by using the full-wave electromagnetic simulation package FEKO [27]. The antenna terminal voltage due to the received field is calculated as 𝑉𝑚=𝑍𝐿𝐼𝑚 where this terminal voltage represents (17). Since terminal voltage measured at every antenna element in the array includes the contribution of the mutual coupling, the proposed compensation method is applied to compensate the effect of mutual coupling. A single (first) row of estimated compensation matrix which includes the elements ̃𝑐1111 to ̃𝑐1193 using proposed method is tabulated in Table 2.

tab1
Table 1: The different CCRA configurations used for obtaining the optimal case.
tab2
Table 2: Values of estimated compensation matrix using the proposed method.
5.1. Performance of CCRA for Different Arrangements of Antenna Elements

In order to analyse the effect of mutual coupling on the near-field parameter estimation using CCRA and chose optimum array configuration, CRLB for different arrangements of CCRA are calculated in the presence of mutual coupling as discussed earlier. Here, we consider four different arrangements of ring radii and element distributions for forming an optimum CCRA as given in Table 1. In every CCRA configuration under consideration, only three concentric rings are considered and the total number of dipole antenna elements in all rings put together is equal to 27. The antenna elements are placed uniformly on the circumference of the each ring. In order to derive CRLB for all the four cases, a single transmitting half-wave dipole is considered to be positioned at a known location of (4𝜆,80,20). The CRLB for estimated azimuth angle, elevation angle, and range for these cases are plotted in Figures 24, respectively. It can be observed, from Figures 24, that the array configuration of case (d) has the lowest CRLB when compared to other cases (a), (b), and (c) and can be considered as an optimal configuration. Hence, we choose the optimal CCRA configuration given in case (d) for demonstrating compensation and parameter estimation.

5.2. Near-Field Parameter Estimation

Here, MUSIC and ESPRIT algorithms are applied for parameter estimation of near-field sources using the optimal CCRA configuration. Zero mean white Gaussian noise is added to the received signal at every element of CCRA. We apply the beamspace transformation on the array manifold of CCRA to map it to an equivalent virtual linear array. The decorrelation technique available in [26] is applied since we consider the received signals to be correlated. For an idealised situation, that is, when the effect of mutual coupling is ignored, the estimated azimuth and elevation angles are shown by their pseudospectra in Figure 5. When the effect of mutual coupling is taken into consideration, without adopting any compensation, MUSIC pseudospectrum fails to create sharp peaks as shown in Figure 6. After applying the proposed compensation method, improved accuracy of parameter estimation can be observed from the peaks of MUSIC pseudospectra as shown in Figure 7. Similarly, in the presence of mutual coupling, that is, when no compensation is applied, ESPRIT algorithm also fails to estimate azimuth and elevation angles properly as shown in histogram plots of Figures 8 and 9, respectively. However, after applying the proposed compensation method, the estimation accuracy of azimuth and elevation angles using ESPRIT improves significantly as shown in histogram plots of Figures 10 and 11, respectively. The estimated ranges for near-field sources are calculated from the estimated time delay which is obtained by using cross-correlation among eigenvectors of signal subspaces. The coefficients of cross-correlations before compensating the effect of mutual coupling are plotted in Figure 12 whereas Figure 13 shows coefficients of cross-correlation after applying the proposed compensation method. Comparing Figures 12 and 13, it is obvious that, after the compensation, the peaks have become distinct and clearer. The optimal CCRA configuration, as explained earlier, is chosen for testing the performance of both MUSIC and ESPRIT algorithms, and their performance is plotted in Figures 14, 15, and 16 along with CRLB. It can be observed from these results that the MUSIC algorithm performs slightly better at lower SNRs when compared to ESPRIT. The RMSE of cross-correlation-based range estimation method is closer CRLB at higher SNRs.

712318.fig.005
Figure 5: MUSIC pseudospectrum for ideal case where the effect of mutual coupling is ignored.
712318.fig.006
Figure 6: MUSIC pseudospectrum when mutual coupling is uncompensated.
712318.fig.007
Figure 7: MUSIC pseudospectrum after the proposed compensation.
712318.fig.008
Figure 8: Histogram plot of estimated azimuth angles using ESPRIT algorithm when mutual coupling is uncompensated.
712318.fig.009
Figure 9: Histogram plot of estimated elevation angles using ESPRIT algorithm when mutual coupling is uncompensated.
712318.fig.0010
Figure 10: Histogram plot of estimated azimuth angles using ESPRIT algorithm after applying proposed compensation.
712318.fig.0011
Figure 11: Histogram plot of estimated elevation angles using ESPRIT algorithm after applying proposed compensation.
712318.fig.0012
Figure 12: The cross-correlation among the eigenvectors of signal subspace when mutual coupling is uncompensated.
712318.fig.0013
Figure 13: The cross-correlation among the eigenvectors of signal subspace after applying proposed compensation.
712318.fig.0014
Figure 14: Performance comparison for the estimation of azimuth angle using MUSIC and ESPRIT algorithms for CCRA design of case (d).
712318.fig.0015
Figure 15: Performance comparison for the estimation of elevation angle using MUSIC and ESPRIT algorithms using CCRA design of case (d).
712318.fig.0016
Figure 16: Performance of range estimation using cross-correlation based time-delay estimation method.

6. Conclusions

A mutual coupling compensation technique that is applicable with subspace-based super-resolution estimation algorithms is proposed in this paper. Concentric circular ring arrays are employed for testing the proposed method of mutual coupling compensation by estimating the near-field parameters of transmitting dipoles using MUSIC and ESPRIT algorithms. In order to achieve a better accuracy on parameter estimation, an optimal CCRA configuration is chosen to obtain lowest CRLB in the presence of mutual coupling. The proposed compensation method is applied with MUSIC and ESPRIT algorithms for estimating azimuth and elevation angles, and the range parameter is estimated from the time delay obtained from the cross-correlation of signal subspace eigenvectors so as to avoid parameter search over 3D space and reduce the associated computational burden. The estimation performance after applying the proposed compensation method of mutual coupling is compared with CRLB for different SNRs. From simulation results, it can be observed that the proposed compensation method can effectively remove the effect of mutual coupling in CCRA and helps to improve the estimation performance of subspace-based algorithms in the presence of array mutual coupling.

Appendix

Elements of Fisher information matrix (FIM) (14) can be obtained as𝐹𝜃𝜃𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑅𝑠𝐴𝐻𝜃𝑅𝑥1𝐴𝜃𝑇,𝐹𝜙𝜙𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑅𝑠𝐴𝐻𝜙𝑅𝑥1𝐴𝜙𝑇,𝐹𝑟𝑟𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑅𝑠𝐴𝐻𝑟𝑅𝑥1𝐴𝑟𝑇,𝐹𝜃𝜃𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑇,𝐹𝜃𝜙𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝐻𝜃𝑅𝑥1𝐴𝜙𝑇,𝐹𝜙𝜃𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝐻𝜙𝑅𝑥1𝐴𝜃𝑇,𝐹𝑟𝜙𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑅𝑠𝐴𝐻𝜙𝑅𝑥1𝐴𝑟𝑇,𝐹𝑟𝜃𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝐻𝑟𝑅𝑥1𝐴𝜃𝑇,𝐹𝜙𝑟𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑅𝑠𝐴𝐻𝜙𝑅𝑥1𝐴𝑟𝑇,𝐹𝜃𝑟𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝐻𝑟𝑅𝑥1𝐴𝜃𝑇,𝐹𝜃𝛼𝑘𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝛼𝑘𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝐻𝛼𝑘𝑅𝑥1𝐴𝜃𝑇,𝐹𝜃𝛽𝑘𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝛽𝑘𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜃𝑅𝑠𝐴𝛽𝐻𝑘𝑅𝑥1𝐴𝜃𝑇,𝐹𝜙𝛼𝑘𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝛼𝑘𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑅𝑠𝐴𝐻𝛼𝑘𝑅𝑥1𝐴𝜙𝑇,𝐹𝜙𝛽𝑘𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝛽𝑘𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝜙𝑅𝑠𝐴𝐻𝛽𝑘𝑅𝑥1𝐴𝜙𝑇,𝐹𝑟𝛼𝑘𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝛼𝑘𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑅𝑠𝐴𝐻𝛼𝑘𝑅𝑥1𝐴𝑟𝑇,𝐹𝑟𝛽𝑘𝑅=2𝑁Re𝑠𝐴𝐻𝑅𝑥1𝐴𝛽𝑘𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑇+𝑅𝑠𝐴𝐻𝑅𝑥1𝐴𝑟𝑅𝑠𝐴𝐻𝛽𝑘𝑅𝑥1𝐴𝑟𝑇.(A.1)

Acknowledgment

The work reported in this paper is supported by the Australian Research Council through a Discovery Project Grant DP 0773234.

References

  1. I. J. Gupta and A. A. Ksienski, “Effect of mutual coupling on the performance of adaptive arrays,” IEEE Transactions on Antennas and Propagation, vol. 31, no. 5, pp. 785–791, 1983. View at Scopus
  2. H. T. Hui, “Improved compensation for the mutual coupling effect in a dipole array for direction finding,” IEEE Transactions on Antennas and Propagation, vol. 51, no. 9, pp. 2498–2503, 2003. View at Publisher · View at Google Scholar · View at Scopus
  3. T. T. Zhang, H. T. Hui, and Y. L. Lu, “Compensation for the mutual coupling effect in the ESPRIT direction finding algorithm by using a more effective method,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 4, pp. 1552–1555, 2005. View at Publisher · View at Google Scholar · View at Scopus
  4. T. T. Zhang, Y. L. Lu, and H. T. Hui, “Compensation for the mutual coupling effect in uniform circular arrays for 2D DOA estimations employing the maximum likelihood technique,” IEEE Transactions on Aerospace and Electronic Systems, vol. 44, no. 3, pp. 1215–1221, 2008. View at Publisher · View at Google Scholar · View at Scopus
  5. K. R. Dandekar, H. Ling, and G. Xu, “Experimental study of mutual coupling compensation in smart antenna applications,” IEEE Transactions on Wireless Communications, vol. 1, no. 3, pp. 480–487, 2002. View at Publisher · View at Google Scholar · View at Scopus
  6. R. S. Adve and T. K. Sarkar, “Compensation for the effects of mutual coupling on direct data domain adaptive algorithms,” IEEE Transactions on Antennas and Propagation, vol. 48, no. 1, pp. 86–94, 2000. View at Publisher · View at Google Scholar · View at Scopus
  7. H. Rogier and D. De Zutter, “Beamforming strategies for compact arrays in mobile terminals using the exact active element pattern method,” Microwave and Optical Technology Letters, vol. 35, no. 3, pp. 201–203, 2002. View at Publisher · View at Google Scholar · View at Scopus
  8. C. K. E. Lau, R. S. Adve, and T. K. Sarkar, “Minimum norm mutual coupling compensation with applications in direction of arrival estimation,” IEEE Transactions on Antennas and Propagation, vol. 52, no. 8, pp. 2034–2041, 2004. View at Publisher · View at Google Scholar · View at Scopus
  9. Q. Yuan, Q. Chen, and K. Sawaya, “Accurate DOA estimation using array antenna with arbitrary geometry,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 4, pp. 1352–1357, 2005. View at Publisher · View at Google Scholar · View at Scopus
  10. A. J. Weiss and B. Friedlander, “DOA and steering vector estimation using a partially calibrated array,” IEEE Transactions on Aerospace and Electronic Systems, vol. 32, no. 3, pp. 1047–1057, 1996. View at Scopus
  11. E. K. L. Hung, “Matrix-construction calibration method for antenna arrays,” IEEE Transactions on Aerospace and Electronic Systems, vol. 36, no. 3, pp. 819–828, 2000. View at Scopus
  12. Q. Bao, C. C. Ko, and W. Zhi, “DOA estimation under unknown mutual coupling and multipath,” IEEE Transactions on Aerospace and Electronic Systems, vol. 41, no. 2, pp. 565–573, 2005. View at Publisher · View at Google Scholar · View at Scopus
  13. S. Kikuchi, H. Tsuji, and A. Sano, “Autocalibration algorithm for robust Capon beamforming,” IEEE Antennas and Wireless Propagation Letters, vol. 5, no. 1, Article ID 874070, pp. 251–255, 2006. View at Publisher · View at Google Scholar · View at Scopus
  14. B. Friedlander and 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. View at Scopus
  15. H. S. Lui and H. T. Hui, “Improved mutual coupling compensation in compact antenna arrays,” IET Microwaves, Antennas and Propagation, vol. 4, no. 10, pp. 1506–1516, 2010. View at Publisher · View at Google Scholar · View at Scopus
  16. H. S. Lui, H. T. Hui, and M. S. Leong, “A note on the mutual-coupling problems in transmitting and receiving antenna arrays,” IEEE Antennas and Propagation Magazine, vol. 51, no. 5, Article ID 5432083, pp. 171–176, 2009. View at Publisher · View at Google Scholar · View at Scopus
  17. H. S. Lui and H. T. Hui, “Mutual coupling compensation for direction-of-arrival estimations using the receiving-mutual-impedance method,” International Journal of Antennas and Propagation, vol. 2010, Article ID 373061, 7 pages, 2010. View at Publisher · View at Google Scholar
  18. C. P. Mathews and M. D. Zoltowski, “Eigenstructure techniques for 2-D angle estimation with uniform circular arrays,” IEEE Transactions on Signal Processing, vol. 42, no. 9, pp. 2395–2407, 1994. View at Publisher · View at Google Scholar · View at Scopus
  19. M. D. Zoltowski, M. Haardt, and C. P. Mathews, “Closed-form 2-D angle estimation with rectangular arrays in element space or beamspace via unitary ESPRIT,” IEEE Transactions on Signal Processing, vol. 44, no. 2, pp. 316–328, 1996. View at Scopus
  20. R. L. Haupt, “Optimized element spacing for low sidelobe concentric ring arrays,” IEEE Transactions on Antennas and Propagation, vol. 56, no. 1, pp. 266–268, 2008. View at Publisher · View at Google Scholar · View at Scopus
  21. Z. Ye and 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. View at Publisher · View at Google Scholar · View at Scopus
  22. P. Stoica and A. Nehorai, “MUSIC, maximum likelihood, and Cramer-Rao bound,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 5, pp. 720–741, 1989. View at Publisher · View at Google Scholar · View at Scopus
  23. T. Svantesson, “Modeling and estimation of mutual coupling in a uniform linear array of dipoles,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '99), pp. 2961–2964, Phoenix, Ariz, USA, March 1999. View at Scopus
  24. T. Huang and A. S. Mohan, “Effects of array mutual coupling on near-field DOA estimation,” in Proceedings of the Canadian Conference on Electrical and Computer Engineering: Toward a Caring and Humane Technology (CCECE '03), pp. 1881–1884, May 2003. View at Scopus
  25. S. H. Lee, C. S. Ryu, and K. K. Lee, “Near-field source localisation using bottom-mounted linear sensor array in multipath environment,” IEE Proceedings: Radar, Sonar and Navigation, vol. 149, no. 4, pp. 202–206, 2002. View at Publisher · View at Google Scholar · View at Scopus
  26. C. P. Mathews and M. D. Zoltowski, “Performance analysis of the UCA-ESPRIT algorithm for circular ring arrays,” IEEE Transactions on Signal Processing, vol. 42, no. 9, pp. 2535–2539, 1994. View at Publisher · View at Google Scholar · View at Scopus
  27. FEKO EM Software and Systems S.A. (Pty) Ltd, 32 Techno Lane, Technopark, Stellenbosch, 7600, South Africa.