International Journal of Antennas and Propagation

International Journal of Antennas and Propagation / 2012 / Article
Special Issue

Mutual Coupling in Antenna Arrays 2011

View this Special Issue

Research Article | Open Access

Volume 2012 |Article ID 712318 |

Mohammed Jainul Abedin, Ananda Sanagavarapu Mohan, "A Subspace-Based Compensation Method for the Mutual Coupling in Concentric Circular Ring Arrays for Near-Field Source Localisation", International Journal of Antennas and Propagation, vol. 2012, Article ID 712318, 13 pages, 2012.

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

Academic Editor: Hoi Shun Lui
Received14 Jun 2011
Accepted05 Sep 2011
Published12 Jan 2012


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 [1โ€“4]. 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๐œŒ๐‘™=๐œ†min๎€ท๎€ท4sin๐œ‹/๐‘€๐‘ƒ๎€ธ๎€ธ,(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,โ€ฆ,โ„Ž.

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 2โ€“4, 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.

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 byฬƒ๐š๐‘“(๐‘Ÿ,๐œƒ,๐œ™)=๐‚๐š๐‘“(๐‘Ÿ,๐œƒ,๐œ™).(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๐‘11โŽคโŽฅโŽฅโŽฅโŽฅโŽฅโŽฅโŽฆโ‹ฏ0โ‹ฎโ‹ฑโ‹ฎ00โ‹ฏ1.(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๎€ป๐‘‡+๐‘21๎€บ0,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๐‘‘+๐‘—๐‘ฃ๐‘‘๎€ธ,๐œƒ๐‘‘=sinโˆ’1๎€ท๎€ท๐‘ข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ร—10โˆ’3๐œ† 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.

CCRA configurationsRadii of ringsNo. of elements on each ring


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 2โ€“4, respectively. It can be observed, from Figures 2โ€“4, 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.

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.


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)


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


  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: Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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: Google Scholar
  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: Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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: Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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 Site | 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 Site | Google Scholar
  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: Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  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: Google Scholar
  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: Google Scholar
  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 Site | Google Scholar
  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 Site | Google Scholar
  27. FEKO EM Software and Systems S.A. (Pty) Ltd, 32 Techno Lane, Technopark, Stellenbosch, 7600, South Africa.

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.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.