#### Abstract

This paper proposes a new method for cross array to estimate two-dimensional direction of arrival (2-D DOA) in the presence of mutual coupling. In this method, the array elements which are affected by the same mutual coupling are chosen on -axis and -axis, respectively. Then a new matrix is constructed with the proper entries of cross covariance matrix of the chosen elements outputs on -axis and -axis. Propagation method (PM) and rotational invariance techniques for uniform linear array (ULA) are utilized in the constructed matrix to obtain two parameters correlated with elevations and azimuths. While calculating and pairing the two parameters, only once eigendecomposing and several division operations are required with the relationship among the matrix, its eigenvalues, and corresponding eigenvectors. Simulations are presented to validate the performance of the proposed method.

#### 1. Introduction

Mutual coupling between elements has a serious effect on precision and resolution of the direction of arrival estimation algorithms [1, 2]. Therefore, there are lots of works addressing this issue of how to restrain the effects of mutual coupling. Many calibration algorithms have been proposed in the past decades. Friedlander and Weiss presented an iterative procedure to compensate the mutual coupling and perturbation of gain and phase [3]. In [4], Sellone and Serra constructed an objective function based on covariance matrix matching, dealt with mutual coupling matrix and its conjugate transpose as unrelated ones, and searched the optimal solution with iterative procedure. However, the multidimensional search needed for the associated nonlinear optimal search is computationally complicated in [3, 4]. Wang et al. made use of banded complex symmetric Toeplitz characteristic of the mutual coupling matrix of ULA and rank reduction (RARE) algorithm to obtain azimuths and mutual coupling coefficients [5], while Lin and Yang utilized banded complex symmetric circular Toeplitz characteristic of the mutual coupling matrix of uniform circle array (UCA) [6]. These two methods are only suitable for ULA and UCA, respectively. The idea from [5] was applied to cross array, but the transformed matrix related to the ideal steering vector has no universal form for reasons of the array configuration and the variable mutual coupling degree of freedom [7]. The literature [8] presented a decoupling algorithm for ULA. It chooses the proper elements whose ideal steering vector can be separated from the mutual coupling coefficients and then employs the outputs of the chosen elements to construct MUSIC spatial spectral estimator. Nevertheless this method can only estimate one-dimensional direction and needs the spectral peak search. In [9, 10], this method [8] was applied to uniform rectangular array (URA) and uniform hexagon array (UHA), respectively, to estimate two-dimensional directions, but the configurations of planar arrays are complicated and lots of auxiliary elements are needed.

In order to estimate 2-D DOA in the presence of mutual coupling, a method for cross array which has a simpler configuration than those for planar arrays is presented in this paper. The method is based on PM algorithm utilizing signal subspace eigenvectors algorithm for the array consisting of the chosen elements which are affected by the same mutual coupling. It is not an iterative algorithm, and the spectral peak search is not required. It can estimate the 2-D DOA accurately without estimating the mutual coupling coefficients.

#### 2. Problem Formulation

Consider a cross array consisting of two ULAs located on -axis and -axis, respectively. The ULAs share the middle element as reference at origin of coordinates. Each ULA has sensors with element spacing of , which means that there exist ones in the whole array because of the sharing element. Assume thatnarrowband, far-field, noncoherent signals with the same wave length () impinge on this array, and the th source’s elevation and azimuth are and , respectively. The array configuration is shown in Figure 1.

Let and denote the mutual coupling matrix (MCM) and the time index, respectively. The output of the array can be expressed as where , , , and represent the received signal vector, the ideal steering matrix, the source signal vector, and the stochastic noise (additive white Gauss noise whose mean is zero and variance is ) vector, respectively. They can be described as (the time indexis omitted for brevity) where superscript denotes transpose.

As shown in [3], the MCM of ULA can be expressed as a banded complex symmetric Toeplitz matrix where integer denotes the mutual coupling degree of freedom, which means the coefficient of mutual coupling is assumed to be zero if the distance between the two sensors exceeds.

Based on the MCM model of ULA by (3), the MCM of cross array has the following form: where is the MCM of ULA on - or -axis, and describes the mutual coupling between the two ULAs. The matrix can be presented by where ,denotes the coefficient of mutual coupling between the two elements spacing , and satisfies

#### 3. Proposed Method

The ideal steering vector of ULA elements affected by the same mutual coupling can be separated from the mutual coupling coefficients, which means they can be decoupled easily. Let and be the outputs of elements affected by the same mutual coupling on-axis and-axis, respectively.

can be written as where ,

There are no entries of in (8) as the chosen elements on -axis are not affected by the sensors on-axis, because the distances between them are out of .

Similarly, can be expressed as where ,

In order to introduce rotational invariance techniques for ULA to this array, is partitioned into two parts, so is . They are denoted by , , , and , respectively:

and are the submatrixes of and , respectively, so the cross covariance matrix of and can be built with the proper entries of cross covariance matrix of and : where and are the ideal steering matrixes of arrays corresponding to and , respectively,

Similarly, let , , and , which can be constructed with the proper entries chosen from as in (16), be the cross covariance matrixes of and , , and , and and , respectively, where and are the ideal steering matrixes of arrays corresponding to and , respectively,

The new matrix consisting of the four cross covariance matrixes includes the rotational invariance factors,

Similar to the conventional propagator method, can be partitioned as where and are the first and the last rows, respectively.

A matrix (the propagator parameter) is defined to satisfy

The solution of (29) can be obtained as follows:

With the relationship between the propagator parameter and ideal steering matrix, we obtain an equation as where () is the matrix and is the matrix with rank .

Similar to ESPRIT, from (31) we can have the equations correlated with azimuths and elevations as

To solve the above equation sets and pair the azimuths and elevations, four times eigen-decomposing and times calculating inner product of vectors are necessary in a general way. Now a method, which needs only once eigen-decomposing and some division operations to obtain and pair the azimuths and the elevations, is introduced with the relationship among the matrix, its eigenvalues, and corresponding eigenvectors.

First eigendecompose , and the eigenvalues and the eigenvectors are denoted by and (), respectively Then substitute into other equations, as the first one of the right equation sets in (32); the corresponding eigenvalue can be obtained from

Let , and can be expressed by where and are the th element of and , respectively.

Assume that and denote the corresponding eigenvalues of , and respectively, and they can be calculated by (34).

Using the relationship among (), , , , and , the elevation and the azimuth of the th source can be represented as

Now give the computation load on the proposed method and PM algorithm with known mutual coupling. The computation load on PM algorithm with known mutual coupling is focused on the inverse of MCM, cross matrix of output of elements on -axis and -axis, eigen-decomposing a matrix with rank , and pairing the parameters, so the computation load is (is the snapshots number); the computation load on the proposed method is focused on cross matrix of output of elements on -axis and -axis, once eigen-decomposing a matrix with rank , and pairing the parameters, so the computation load is .

Finally, the necessary conditions for the proposed method are given as follows:(1)for unique estimation, is required;(2)the sources can not be located at the blind angles which make the determinant of or zero.

#### 4. Simulation Experiments

In this section, we carry out some representative simulation experiments to demonstrate the validity of the proposed calibration method. In all the experiments, we assume there are two uncorrelated equivalent-power sources at (20°, 10°) and (80°, 60°) that impinge on the cross array which consists of 53 elements (). The mutual coupling degree of freedom is two and the coefficients of mutual coupling are denoted by , , and .

In the first experiment, the 2-D spatial spectrums of the PM algorithm with unknown mutual coupling, the PM algorithm with known mutual coupling, and the proposed method (estimate the 2-D DOAs to get the coefficients of mutual coupling with the method in [7]) are shown in Figures 2 to 4 in the condition that the signal-to-noise ration (SNR) is 5 dB and the number of the snapshots is 300.

From Figure 2 we can see that the PM algorithm without calibration has very poor resolution in the presence of mutual coupling. Though the peaks appear at the signal directions, there exist lots of pseudopeaks especially at (80°, 10°) whose spatial spectrum is greater than the sources’, which may result in the inaccurate estimation.

Comparing with Figure 4 and Figure 3, the spatial spectrum of proposed method has two sharp peaks at the signal directions, which is also shown in Figure 3. In other words, the proposed method can restrain the effects of mutual coupling and has a high resolution. Though the performance of proposed method is worse than that in Figure 3 for reasons of the array elements consumption, it eliminates pseudo peaks caused by mutual coupling, which can distinguish the two sources visibly.

In the second experiment, the effect of SNR on the performance of proposed method is explored. The number of snapshots is 300 and the SNR varies from −10 to 10 dB with the interval 1 dB. The root mean square error (RMSE) which is defined as (36) is shown in Figure 5 with 100 independent experiments:

The simulation results show the proposed method has an obvious performance improvement over the PM algorithm with unknown coupling, though it is a little inferior to the PM algorithm with known coupling as we think. As the SNR increases, the difference of the RMSE between the proposed method and the PM algorithm with known coupling gets closer; however, the accuracy of DOA estimates of the PM algorithm with unknown coupling stays at a low level as the main error is caused by the manifold error between the ideal and the actual one.

In the last experiment, the influence of the number of the snapshots on the performance of proposed method is investigated. SNR is 0 dB and the number of the snapshots varies from 100 to 800 with the interval 100. 100 independent experiments are carried out, and the RMSE curves are shown in Figure 6.

From Figure 6 we can see that the result is similar to that of the second experiment: the capability of proposed method is close to the one of the PM algorithm with known coupling, and the performance of the PM algorithm with unknown coupling is still poor when the number of the snapshots increases.

From the above three experiments, it is seen that the proposed method has compensated the effect caused by mutual coupling. It provides a significant improvement of DOA estimation accuracy and can achieve almost the same performance as that of PM algorithm with known mutual coupling.

#### 5. Conclusion

A method compensating the effect of the mutual coupling is proposed for cross array to estimate the 2-D DOA. This method is based on PM algorithm and rotational invariance techniques for ULA. It can compensate the effect of the mutual coupling without any calibration sources, and it does not require any iterative procedure or spectral peak search. Simulations have shown that the proposed method can improve the DOA estimation accuracy well.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors would like to thank the editor and reviewers for their helpful comments and suggestions. In addition, this work is supported by the National High Technology Research and Development Program of China (863 Program) under Grant 2012AA0601101, the Key Laboratory of Intelligent Perception and Systems for High-Dimensional Information (Nanjing University of Science and Technology), Ministry of Education under Grant 3092013012205, and Research Fund for the Doctoral Program of Higher Education of China under Grant 20133219110027.