#### Abstract

A short review of the receiving-mutual-impedance method (RMIM) for mutual coupling compensation in direction finding applications using linear array is conducted. The differences between the conventional-mutual-impedance method (CMIM) and RMIM, as well as the three different determination methods for receiving mutual impedance (RMI), will be discussed in details. As an example, direction finding with better accuracies is used for demonstrating the superiority of mutual coupling compensation using RMIM.

#### 1. Introduction

Signal processing algorithms of most antenna array applications usually assume that interactions of signals and antenna elements are independent to each other. Such assumption is not necessary to be true especially when the antenna elements are closely spaced. For transmitting applications, the transmitted signals from an antenna element would be received by the other elements in the array and reradiate. The entire transmitted signal is not just the sum of the individual signals from each antenna element, but also the reradiate field from each antenna element. Similarly in receiving applications, electromagnetic wave is first received and current is then induced on an antenna element. The induced current on the antenna element reradiates electromagnetic field which would be received by other elements on the array. The wavefront of the incoming signal is distorted by the reradiated field. In both cases, interactions between the antenna elements distorted the wavefront of the signals. Such effect is known as mutual coupling and it is usually considered as a defect which degrades the performance of array signal processing algorithms.

Existence of the mutual coupling effect in antenna array has been well known throughout the years and extensive work has already been contributed to either reducing or compensating such undesirable effects. Gupta and Ksienski [1] used a circuit theory approach based on the use of the conventional mutual impedance (CMI) [2] to analyze the effect of antenna mutual coupling. This method, namely the conventional-mutual-impedance method (CMIM), has become the most popular method for mutual coupling analysis as CMIs can easily be measured directly or obtained indirectly from the measurement of the S-parameters. Later, Yeh et al. [3] applied the CMIM to decouple the voltages measured from the antenna terminals of a receiving antenna array in a Direction-of-Arrival (DOA) estimation application. The results with better accuracies of DOA estimation demonstrate that the undesirable mutual coupling effect can be partially compensated via a decoupling process with CMI.

Throughout the last two decades, significant amount of interest has been contributed to mutual coupling compensation in DOA estimations and adaptive nulling. Various methods such as calibration methods [4–6], full-wave method [7–9], and receiving-mutual-impedance method (RMIM) [10–15] have been proposed and demonstrate better mutual coupling compensation performance over the CMIM method. A recent review of various mutual coupling methods can be found in [16]. Direct comparisons of these methods for the same DOA estimation problem can be found in [17]. The superiority of these methods over the CMIM relies on the fact that the mutual coupling of the receiving array is considered, as compared to the CMIM which the CMIs are determined when the array is in transmitting mode. The differences of the mutual coupling problems for transmitting and receiving arrays have been covered in Balanis [18] and Daniel [19]. The differences between the CMIM, RMIM, full-wave method and calibration method have recently been discussed in Lui et al. [20].

In this paper, a review of the RMIM in the context of DOA estimation using linear array is presented. The concept of receiving mutual impedance (RMI) was first introduced by Hui [10, 11] and applied on uniform linear array (ULA) for DOA estimations and adaptive nulling. A review of the RMIM will be given in details in the next section. In particular, the formation of RMIM, the “single-mode” approximation [10] as well as the difference between RMIM and CMIM will be discussed in details. Throughout the years, there have been three methods to determine the RMI for linear arrays [12, 13, 15]. Discussions of these determination methods will be given in details in Section 3. A DOA estimation example using Matrix Pencil Method (MPM) [9, 21, 22] with better accuracies demonstrate the superiority of mutual coupling compensation of RMIM.

#### 2. Receiving Mutual Impedance

Consider a simple case of a two-element dipole array. According to the standard MoM analysis, the entire electromagnetic problem of the antenna array with each antenna discretized into segments can be written as [10, 23]

where the matrix blocks are

respectively. and ranges from 1 and 2 for the case of two antenna elements. The elements in the matrix, , correspond to the impedance between the segment of antenna and segment of antenna , corresponds to the load impedance of the antenna terminal. The elements and of the and matrices, respectively are the induced current and the excitation voltage to the corresponding *i*th segment of the antenna , where . Expanding from (1), this gives

As we can see, the mutual coupling effect is due to the off-diagonal matrix blocks, that is, with . If the entire current distribution of the antenna and the full knowledge of the all elements in the impedance matrices are known, the mutual coupling effect can be accurately quantified and removed [7–9]. However, this can only be done in the modeling domain. In practice, it is not possible to measure the entire current distribution of the antenna, that is, all the values in the matrix. Only the terminal voltages of the antenna terminals, that is, can be measured. As a result, Hui [10] has introduced the RMI based on the “single-mode” assumption, that is, we treat the current induced on the each antenna element of the array consist of a single-mode current based on the current at the antenna terminal. For a two-element receiving array, this can be written as

where and are defined as the RMI between the two antennas, and are the induced current in the antenna terminals, and are the induced voltages in the antenna terminal 1 and 2, respectively, and and are the terminal voltages solely based on the incoming signals without any mutual coupling effect. The subscript denotes that the RMI is defined based on the current and voltage parameters at the antenna terminal.

It is apparent that under the “single-mode” assumption, the RMIs defined in (4a) and (4b) are no longer the elements of the matrix in the original integral equation formulation in (1). As shown in (4a) and (4b), the RMI varies as the induced terminal current and voltages of the antennas, as well as the load impedance (if in (3) changes, the entire electromagnetic problem changes and so as the RMI). For omnidirectional antenna elements, the current distribution remains to be the same for any angle of azimuth angle in the elevation plane of ° (horizontal plane) [13].

Now we extend the concept to antenna array with antenna elements and each of them are terminated with the same load impedance . The measured voltage at antenna terminal can be written as

where Once again is the measured terminal voltage due to the incoming signal alone, is the voltage due to the mutual coupling effect from other array elements, is the RMI between antenna elements and , is the current induced at the antenna terminal, given by for . Rearranging (5a) and (5b) the relationship between s and s can be written as

With a *priori* knowledge of the RMIs together with the terminal voltages s from the DOA problem, the uncoupled voltages s can thus be determined using (6) and the undesirable mutual coupling effect is thus compensated.

#### 3. Comparision with the Conventional-Mutual-Impedence Method

Other than the RMIM, another well establishing method, the CMIM has been well known in the literature. The differences between the two methods have not been well discussed and we take this opportunity to give a clear explanation. According to [1], the relationship between the voltages at terminal with other antenna terminals can be written as

where is the CMI between antenna terminals and , corresponds to the self impedance of the antenna element , corresponds to the induced current at the antenna terminal and is defined as the open-circuit voltage at terminal when all other antenna elements under open-circuit condition, that is, for . The relationship between the terminal voltage and current can be given by

where is the load impedance terminated at the feed. With the open-circuit condition and substituting (8) into (7), the relationship between the open-circuit voltages and terminal voltages can be written as Comparing CMIM and RMIM, there are three main differences between the two. They are:

(1)The CMI is determined when the antenna array is in transmitting mode while the RMI is determined when the array is in receiving mode. As discussed in [18–20], the mutual coupling problems of transmitting and receiving arrays are two different problems in general and they should be treated separately.(2)The CMI is defined when the antenna element is open circuit while the RMI is defined when the antenna is loaded. The current distribution of the antenna varies with different load impedance. The antenna is usually terminated with load impedance in practice. Strictly speaking, an open-circuit assumption provides only an approximate solution to the mutual coupling characterization. (3)The two mutual impedances are defined differently and thus the decoupling matrices are different ((6) and (9)). In particular, notice that the knowledge of self-impedance is not required for the RMIM.#### 4. Determination Methods of the Receiving Mutual Impedance

A discussion of the three determination methods of the RMI is covered. The motivations, differences and limitations between these methods are discussed in details.

##### 4.1. Method 1

In this method, the RMI is determined by first taking the current distribution of an antenna element under plane-wave excitation. Here, a two-element dipole array is considered and the incident wave is coming from the elevation plane of ° (perpendicular to the dipole element). Analytical expression of the current distribution can be found in the literature [10]. In the computational domain, such current distribution is then used to excite an antenna element in the array and the coupled induced terminal currents and voltages in other antenna elements can then be obtained. The terminal voltages that are solely under plane wave excitation are determined by considering the current distribution of an antenna element with antenna *k* removed. The RMIs of the two antenna elements *k *and* i *can be calculated by rearranging (4a) or (4b) as

The limitation of this method is that the accuracy of the determined RMI is strongly dependant on the accuracies of the approximated current distribution. In practice, the current distribution varies when are more than one signal or when the incoming elevation angle deviates from °. Furthermore, another limitation is that the RMI can only be determined via numerical modeling. Alternative solution is thus required if one would like to measure the RMI.

##### 4.2. Method 2

In this method [12, 13], the two antenna elements *k *and* i *are first illuminated by plane-wave excitation and the terminal voltages and are obtained. Next, antenna element *i* is removed from the array and the terminal voltage of antenna *k* that is solely under plane wave excitation, , are thus obtained. The RMI, , can then be determined using (10). The study in [12] has shown that the RMI can be measured experimentally and the results are comparable to the numerical modeling using MoM. This method provides an accurate solution as it does not require any assumption of the current distribution. Furthermore, the RMI can be obtained experimentally.

The main drawback of both method 1 and 2 is that only two antenna elements are considered at a time and the scattering effect from other elements (if there are more than two elements in the array) are not taken into account. According to (5a) and (5b), the mutual coupling problem for a *N*-element receiving array can be written as

Consider all the voltages and current values are found as described above, there are unknowns in (11) ( RMIs) which cannot be solved mathematically. As a result, only two antennas can be considered in one time in the previous two methods and removal of “other” antenna elements is needed when . Strictly speaking, this method only provides an accurate solution for the case of two antenna elements. When there are more than two antennas, the “extra” antennas have to be removed and thus the electromagnetic problem has changed. The scattering effect from “other elements” has not been taken into consideration and thus alternative solution is required to handle such drawback.

##### 4.3. Method 3

For antenna arrays with omnidirectional antenna elements, the RMI remains to be a constant as the azimuth angle varies. Consider the terminal voltages are measured at and repeat the measurement at incoming aspects, (11) can be written as where the superscripts on the voltages and currents terms correspond to the voltages and currents measured under the plane-wave illumination from the azimuth angle of , respectively. The RMI s can be determined by solving (11), where .

As previously mentioned, the first advantage this method does not require the removal of antenna elements such that the scattering effect from all antenna elements is taken into consideration in the calculation when we calculate/measure . Preliminary results in [15] have first demonstrated that better accuracies can be obtained for DOA estimations using the RMI determined based on method 3 as compared to method 2. More importantly, as the element separation of the antenna elements becomes smaller, the scattering effect from other elements becomes significant and RMI determined using method 2 may fail to compensate the mutual coupling effect.

To demonstrate the superiority of method 3, a DOA estimation example using two 8-element uniform linear arrays with element separations of and is considered. Decoupling using CMIM is not included here as it has previously been demonstrated [10, 11, 14, 15] that RMIM is better than CMIM and we are not going to repeat it here again. Normal-Mode Helical Antenna (NMHA) element shown in Figure 2 is considered. Three signals coming from the azimuth directions of °, °, and ° (all from the elevation angle of °) shown in Figure 1 are considered. The incoming signals are contaminated with Additive Gaussian White Noise (AGWN) with a Signal-to-Noise Ratio (SNR) of 10 dB. The entire simulations are done in full-wave electromagnetic solver FEKO [24] so that the scattering effect from all antenna elements is taken into account. The signals are first compensated using the RMIM based on the RMI determined using method 2 and a Matrix Pencil Method (MPM) algorithm is used for DOA estimations. The results out of 1000 simulations are tabulated in Tables 1 and 2. For the case with an element separation of , there is not significant differences for the estimated incoming directions when the RMI is determined using the two methods. However, as the element separation is reduced to , as shown in Table 2, it is found that the estimated directions of signal 1 and signal 2 becomes ° and °, respectively, which are not correct. Using the RMI determined using method 3, three incoming signals can be clearly resolved with less than ° of bias, standard deviations, root-mean-square errors (RMSEs).

Further investigations about the solution of the RMI determined using different number of incidence, can be found in our recent publication [25]. The only limitation of this method is that it is only applicable to omnidirectional antenna elements. Alternative solutions are thus required for nonomnidirectional antenna elements [26, 27].

#### 5. Conclusions

A review of the RMIM for mutual coupling compensation in the scope of DOA estimations using ULA is presented. In particular, the definition of RMI, the “single-mode” approximation, as well as the differences between CMIM and RMIM are discussed in details. The advantages and limitations of the three determination methods of RMI have also covered followed by a DOA estimation example which demonstrates the feasibilities and superiorities of mutual coupling compensation using RMI. In addition to theoretical studies, we have recently demonstrated experimentally that the RMIM can be used for mutual coupling compensation [28].

The RMIM has also been applied to mutual coupling compensation of DOA estimations using circular array [29] and phase array in Magnetic Resonance Imaging (MRI) [30–32]. Recently, the RMI has been adopted to model the mutual coupling effect of Multiple-Input Multiple-Output (MIMO) system and new models for channel capacities of the MIMO systems have been developed [33, 34]. It is no doubt that the RMIM is an important model for mutual coupling characterization and it is not difficult to see that such model can be used in other radar, biomedical, and communication applications.

#### Acknowledgments

This work was supported in part by the National University of Singapore (NUS) under Grant no. R-263-000-469-112 and in part by the US ONR research fund under the Project no. PR. NO. 09PR03332-01.