Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 969042, 11 pages
http://dx.doi.org/10.1155/2015/969042

Research Article

CC-MUSIC: An Optimization Estimator for Mutual Coupling Correction of L-Shaped Nonuniform Array with Single Snapshot

Yuguan Hou,¹ Tongyu Zhang,² and Shaochuan Wu¹

¹School of Electronics and Information Engineering, Harbin Institute of Technology, 92 West Dazhi Street, Nangang, Harbin 150001, China
²Institute of Engineering Mechanics, China Earthquake Administration, No. 29 Xuefu Road, Nangang, Harbin 150080, China

Received 2 September 2014; Revised 6 February 2015; Accepted 6 February 2015

Academic Editor: Massimo Scalia

Copyright © 2015 Yuguan Hou et al. 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.

For the case of the single snapshot, the integrated SNR gain could not be obtained without the multiple snapshots, which degrades the mutual coupling correction performance under the lower SNR case. In this paper, a Convex Chain MUSIC (CC-MUSIC) algorithm is proposed for the mutual coupling correction of the L-shaped nonuniform array with single snapshot. It is an online self-calibration algorithm and does not require the prior knowledge of the correction matrix initialization and the calibration source with the known position. An optimization for the approximation between the no mutual coupling covariance matrix without the interpolated transformation and the covariance matrix with the mutual coupling and the interpolated transformation is derived. A global optimization problem is formed for the mutual coupling correction and the spatial spectrum estimation. Furthermore, the nonconvex optimization problem of this global optimization is transformed as a chain of the convex optimization, which is basically an alternating optimization routine. The simulation results demonstrate the effectiveness of the proposed method, which improve the resolution ability and the estimation accuracy of the multisources with the single snapshot.

1. Introduction

The L-shaped linear array generally consists of two mutually perpendicular configuration uniform linear arrays, which could be horizontal or vertical. It can be used to estimate the target signal’s two-dimensional information of the elevation angle and the azimuth angle. In addition, its structure is simpler than the planar linear array. In practical applications, the antenna performance of the L-shaped linear array is influenced by the mutual coupling between two array elements [1]. Ignoring the impact of the mutual coupling, the direction of arrival (DOA) estimation algorithm makes a serious degradation [2]. Reducing the impact of the mutual coupling of an antenna array, therefore, becomes an important task. The mutual coupling correction approach is a research hotspot in recent years. In the past two decades, many array calibration algorithms have been proposed for the mutual coupling problem. In [2], the authors pointed out that the mutual coupling reduced the eigenstructure decomposition algorithm performance, and a method of compensation and correction was introduced. In [3], a maximum likelihood algorithm for compensating the mutual coupling of the array element gain and phase has been proposed. This scheme considered all the obstacles without the calibration correction; however, a set of calibration sources in known locations were required. In [4], Lin and Yang considered the uniform circular arrays and presented a blind calibration method for the mutual coupling between elements. In [5], a method that links the DOA estimation and correction of the uniform linear array was presented. This method employed the Toeplitz structure of the array covariance matrix to compensate the gain and phase of the array elements. In [6], the self-correction of the array’s mutual coupling was achieved by setting the amount of the auxiliary array element. Furthermore, the problems caused by iteration were also avoided. In [7], an online mutual coupling compensation algorithm for the uniform and linear arrays was presented. It could simultaneously compensate for the mutual coupling and estimate the direction-of-arrivals of signals. An alternating minimization procedure based on the closed-form solutions was performed to estimate the mutual coupling matrix in the field of complex symmetric Toeplitz matrices. In [8], a method for the estimation of the direction of arrival in the presence of multipath propagation and the mutual coupling for a frequency hopping system was proposed. An iterative alternating minimization algorithm for finding the mutual coupling and the DOA parameters in an alternate manner was formulated by using the pilot symbols and assuming perfect time-frequency synchronization for a linear array. In [9], a new mutual coupling compensation method based on the minimum norm solution to an underdetermined system of equations was introduced. The formulation was proved to be independent of the type of the antenna element and provide good results in situations where signal strengths vary considerably. The analysis of mutual coupling was applied in the context of a code division multiple access communication system. In [10], a computationally efficient algorithm for the direction of arrival estimation of uncorrelated sources and for self-calibration of mutual coupling between the sensors consisting of several uniformly spaced subarrays was proposed. In [11], a simple method was presented for the estimation and the compensation of mutual coupling in antenna systems of arbitrary geometries, including antenna arrays in the vicinity of scatters. The method included both theoretical and experimental schemes, whereas it did not resort to assumptions that often encountered in the previous mutual coupling estimation approaches. In [12], a decoupled method for 2D direction of arrival estimation in the presence of the elevation-dependent mutual coupling was proposed for the compact uniform circular arrays (UCAs) based on the rank reduction theory. In [13], a very simple but effective MUSIC DOA estimation algorithm was introduced by decoupling the antenna mutual coupling in the coupled noise component with the assumption that the uncoupled noise power can be determined a priori and can be removed from the array received power. In [14], the mutual-coupling problems in transmitting and receiving antenna arrays were revisited. The differences between the mutual coupling and mutual impedances for transmitting and receiving antenna arrays are explained. In [15], provided the angularly-independent mutual coupling was treated as angularly-dependent complex array gains, the middle subarray was found to have the same complex array gains. Consequently, a way for parameterizing the steering vector was proposed and the corresponding method for joint estimation of DOAs and mutual coupling matrix using the whole array data was derived based on subspace principle. In [16], the system identification method was applied to the wideband mutual coupling compensation of the receiving arrays. Using the receiving mutual impedances of an antenna array which was calculated at different frequencies, a multiport compensation network was identified for the wideband mutual coupling compensation.

The above-mentioned methods for the mutual coupling correction are based on the case of the multiple snapshots. For the case of the single snapshot, the integrated SNR gain could not be obtained without the multiple snapshots, and the performance of those methods degrades significantly. Therefore, it is a challengeable work to calibrate the mutual coupling in the case of the single snapshot with the lower SNR level. Moreover, for the correlated signals and the nonuniform array geometry cases, an additional optimization for the interpolated matrix is required. The main contribution of the paper is a joint optimization method for the online calibration of the mutual coupling in the cases of the single snapshot, the correlated signals, and the L-shaped nonuniform array geometry. A Convex Chain MUSIC (CC-MUSIC) for the L-shaped nonuniform array with single snapshot is proposed. It uses a novel iterative approach to improve the self-calibration algorithm for the mutual coupling correction. Neither the prior knowledge of the correction matrix initialization nor the calibration source with the known position is required for the proposed algorithm. The virtual array interpolation method is employed to transform the L-shaped linear nonuniform array to the virtual L-shaped linear uniform array for the spatial smoothing technique, which achieves a good DOA resolution and estimation for the coherent source signals. Subsequently an optimum solution for the approximation between the no mutual coupling covariance matrix without the interpolated transformation and the covariance matrix with the coupling and the interpolated transformation is derived. Moreover, a global optimization problem is formed for the mutual coupling correction and the spatial spectrum estimation. To avoid solving the nonconvex optimization problem, a chain of the convex optimization is adopted.

2. Signal Model of Horizontal L-Shaped Nonuniform Linear Array for Coherent Sources

Without loss of generality, we have the following assumptions to simplify the problem: the number of sources is smaller than that of the array elements, thus assuming that the number of the rows of the array bearing matrix is greater than that of the columns. The target signals are the far-field narrowband signals. For the horizontal L-shaped linear array [1], we assume that the reference element is at the origin of the antenna array coordinate axis as shown in Figure 1. The horizontal L-shaped nonuniform linear array is formed based on the horizontal L-shaped uniform linear array and some elements of which are failed or do not exist due to the restriction of the location place.

Figure 1: Horizontal L-shaped nonuniform linear array.

Denote the azimuth and elevation angles as and , respectively. Let the minimum unit of the adjacent distance between the elements be , , where is the radar wavelength. Therefore, the adjacent distances between the elements are , , , and so on. Let the number of the array elements on the -axis be , and let the number of the array elements on the -axis be . Then the total number of the array elements is . Let the array element coordinate be ; then the phase difference between the th element and the reference element is

where

where indicates the distance between the th array element and the 1th element (i.e., the reference element).

We assume that the far-field narrowband signals enter into the antenna array. The th signal of the reference element with the single snapshot could be denoted as

where denotes the complex envelope of the signal and represents the phase of the received signal. The signal on the th array element is denoted as

where is the radar working frequency, is the delay of the th signal reaching the th array element, and is the additive Gaussian colored noise at the snapshot of the th array element, which is correlated with other array elements.

Without the mutual coupling error, the steering vector of the L-shaped non-uniform linear array could be written as

where represents the transpose of matrix/vector and and are the steering vectors of -axis and -axis, respectively. We have

In order to facilitate the analysis on the mutual coupling case, we suppose . The total number of array elements is . Denoting as the vector of the single snapshot of the array output data, we have

where

where   ( and ) is the phase difference between the th element and the reference element for the th signal source. is the bearing matrix of the array. is the signal vector. is the noise vector. is the dimensional mutual coupling matrix, which is defined as

For the L-shaped linear array which consists of two mutually perpendicular nonuniform linear arrays, is the dimension mutual coupling matrix within -axis or -axis linear array and is the dimension mutual coupling matrix between the two linear arrays.

The mutual coupling matrix , which has the symmetric structure, could be described by the vector , where is the mutual coupling coefficient between the th array element and the th array element and is the distance. The elements of the matrix are expressed as

where , and is the element in the th row and th column of the mutual coupling matrix . When the array element spacing is greater than , the mutual coupling coefficient between the antenna elements is approximately zero, where is the mutual coupling linear degree of freedom. Therefore, the mutual coupling matrix within the ULA could be represented by a unique characterization of the -dimensional vector. Similarly, the symmetric matrix could be described by

where , and is the element in the th row and -column, is the mutual coupling linear degree of freedom between the linear arrays (similar with ). is the mutual coupling coefficient between the 1st element and the other elements within the same linear array. is the mutual coupling coefficient between the 1st element and the elements on the another linear array (the coupling coefficients between the two linear arrays). The matrix could be divided into two separate parts. Some identical elements with those of are included in the 1st part , and the elements of the 2nd part are different from those of ; that is, . could be described by the vector , where is the mutual coupling coefficient between the th array element and the th array element of the another linear array, and the upper triangle part of the symmetric matrix is filled with the elements of the vector. Thus, the mutual coupling matrix is composed of the parameters including and . And we have

The array covariance matrix of the single snapshot could be expressed as

where indicates the covariance matrix of the noise.

In the case of the single snapshot, the multisources need to be treated as the coherent sources for a better resolution. Therefore, the virtual array technique [17, 18] of the nonuniform linear array is used for the resolution of the coherent sources signals. To design the interpolated array, the field of view of the array is divided into some sectors. For example, two sectors are defined by the interval and . Denote the sets of angles and as

where and are the azimuth and elevation angle intervals of the sectors, respectively. These angles are used only in the design of the interpolation matrices. The steering vectors associated with the sets and for the given array are computed and arranged in the matrices form as follows:

In other words, and are the sections of the array manifold of the real array. Denote and as the sections of the interpolated virtual array manifold computed for the sets of angles and :

where and are the responses of the interpolated virtual array of -axis and -axis. The interpolated virtual array is an L-shaped uniform linear array, which has the same number of elements as the real array. However, the adjacent distances of the elements are different. Assuming that there exists the interpolated transformation matrices and , we have

In fact, the interpolations are not exact; therefore, the equations above are not really hold. The “best” interpolation matrices are those giving the best fit between the interpolated response and the desired response and between the interpolated response and the desired response . Denote the L-shaped interpolated transformation matrix as

Given the interpolation matrix we can calculate the array data vector of the virtual array

The interpolation matrix is equal to

where represents the conjugate transpose of . As and are known and the number of elements of the virtual array is equal to that of the real array, the matrix is determined by the value of the adjacent distances of the elements of the virtual array. With the transformation as shown in (19), the L-shaped nonuniform linear array is transformed into the virtual L-shaped uniform linear array. Therefore, the forward-backward smoothing technique [19] could be employed for the resolution of the coherent sources.

Dividing the virtual array into virtual subarrays within the -axis and the -axis, respectively, the number of the elements of the subarray is . The data vector of the th () subarray within the -axis and -axis is denoted as

where is the th element of the data vector of the virtual array within the -axis and () is the th element of the data vector of the virtual array within the -axis. Thus, the covariance matrix of the th () subarray with the single snapshot is expressed as

The forward-backward smoothing covariance matrix is denoted as

where denotes conjugation and is exchange matrix with the antidiagonal elements being 1 and the other elements being 0. The characteristics of the covariance matrix could be obtained through the eigendecomposition. With the noise subspace matrixes and , the MUSIC [20] spatial spectrum estimation equations of the -axis and the -axis are shown in the following equation:

where is a given elevation angle and is a given azimuth angle.

3. A Convex Chain Approach

For the covariance matrix of the th () subarray of the -axis or the -axis, we define the no mutual coupling covariance matrix without the interpolated transformation as . The covariance matrix with the coupling and the interpolated transformation as , and its approximation is denoted as

where or , . The solution of could be obtained by solving the following problem:

In the single snapshot case and the SNR is not high enough. The noise component could not be ignored and the noise covariance matrix could not be expressed with the identical matrix of (26). Equation (26) is rewritten as

In order to alleviate the influence of the noise component, the Toeplitz approximation is introduced for the which is transformed into the virtual uniform linear array. Therefore, we denote the Toeplitz matrix as

Equation (27) is, therefore, rewritten as

The correlation matrix is composed of , and is obtained from the matrix . In order to emphasize the characteristic parameters of , it is denoted as

As is the complex conjugate Toeplitz matrix, it can be written as

where represents the th column of the matrix , and for any column vector matrix can be written as

where the matrixes and are defined as

where , . An equivalent and more compact form could be obtained by defining the following equation:

Thus we can get

Substituting (35) into (29), we have

The complex gradient of the function to be minimized is equal to zero; we have

where is an column vector made by all zeros. Firstly, the cost function is expanded as

According to the theory of complex gradient, the quantity should be treated as a constant, while the only variable term is . Therefore, it is derived that the complex gradients of the addends in (38) are given by

where is a sparse selection matrix denoted as

For any complex matrix and any complex vector , the following equation is true:

Therefore, (37) can be rewritten as

By exploiting the fact that

for any complex matrix and any complex matrix , (42) can be rewritten as

As in [6], the operator vecb(·) is denoted as

for any complex block matrix ,

where is the th block. Equation (44) can be rewritten as

The optimum solution of the optimization problem (37) is given by

When and have different values, the results of the spatial spectrum estimation are different. The optimal selection of and could lead to the highest spatial spectrum amplitudes. Therefore, the global optimization problem is formed as

where (·) is the MUSIC [20] spatial spectrum, which is constructed by and

For this nonconvex problem, we use the iterative alternating chain, which minimizes 1/ alternately over the two matrices and . Suppose we have obtained the two matrices and at the th iteration. Then and can be solved by alternately optimizing the following two subproblems:

Both (51) and (52) are convex. This observation motivates us to use the convex chain approach to solve the original nonconvex optimization problem in (49). We now consider the iterative alternating optimization routine

We can start from and . Here, and are the initial vectors, and we define the covariance matrix with the coupling of   and of the th and ()th iterative routine as

while, practically, “strong” convergence is achieved when

where is the Frobenius norm. An iterative procedure should be adopted, which is shown in Algorithm 1. In each iteration of the alternating optimization routine, a convex optimization is performed. Therefore we refer to this algorithm as Convex Chain MUSIC (CC-MUSIC).

Initialize the virtual L-shaped uniform linear array with the same elements number as the real array, the distance of the

adjacent elements as a random number between 0.01λ to 0.5 λ. We obtain the with (20) and initialize by the estimation

method of the mutual coupling matrix in [6] without auxiliary elements.

For do

        (1) Compute the array data vector of the virtual array by (19).

        (2) Compute the forward-backward smoothing covariance matrix by (23).

        (3) Perform the local optimization for by (51) and the local optimization for by (52).

        (4) If “strong” convergence is achieved by (55), the loop is break and the optimization of the spatial estimation is achieved

             by application of MUSIC algorithm with the covariance matrix obtained by (54). Otherwise, the loop continues.

end for

Algorithm 1: CC-MUSIC.

Remark 1. It should be pointed out that the CC-MUSIC algorithm does not always converge to the real DOAs in the case of low SNR. Although its convergence value is not accurate, it is a simple and effective approach to achieve satisfactory subspace estimation. If it oscillates around a lower level value, this means that the objective function has attained a lower value. Then the algorithm is terminated and the subspace estimate can be obtained from the minimum point among all the iterations. The estimate is an approximation of the real DOAs. The number of the loops is related to the convergence coefficient in (55). If is much smaller, the convergence is “strong” and the number of the loops would be over thirty. To compromise between the convergence speed and the convergence accuracy, we could set a maximum of the number of the loops, such as 30. As the elements of the vectors and are related to the distance of two antenna elements, we could define a linear factor and we have for each search. Hence the number of the search parameter of is reduced to one. Moreover, is determined by one parameter which is the value of the adjacent distances of the elements of the virtual array. Hence the computational complexity of two searches in each loop is low.

4. Simulation and Analysis

In this section, we conduct several simulations to assess the proposed approach. We assume that the received array is an L-shaped nonuniform linear array with   -axis and   -axis receiving antenna elements, respectively. The antenna elements are of the equal gains and are spaced with the base distance of a half wavelength. We compare the performance of CC-MUSIC algorithm with three other approaches. They are the approaches which combine the MUSIC algorithm based on the spatial smoothing technique with the virtual array technique in the presence of the unknown coupling, the known coupling, and the estimated coupling coefficients by the estimation method in [6] without auxiliary elements. We consider targets with the same ranges and azimuths but with the different elevations.

In the first simulation, we set . Both of the distance aggregates of the adjacent elements are {1 1 1 1 1 2 1 1 3 1 1 1 1 1 3 1 1 1 1 2 1 1 1 1 1 2 1 1 1 2 1 1 1 3 1 1 2 1 1 1}, where the “1” is a half wavelength, the “2” is two times of a half wavelength, and “3” is three times of a half wavelength. The vector , and . Moreover, and the vector . We consider three targets with the different elevations which are 3°, 6°, and 20°. The spatial spectrum of the first simulation with 5 dB SNR in each receiving antenna element is shown in Figure 2. It is shown that the proposed approach is effective for the estimation of the three targets. However, the three other approaches fail to resolve the targets. The spatial spectrum of the first simulation with 10 dB SNR in each receiving antenna element is shown in Figure 3. Although all of the four approaches are effective for the estimation of the three targets, the proposed approach achieves higher spatial amplitude than the other approaches.

Figure 2: Spatial spectrum of the first simulation with 5 dB SNR.

Figure 3: Spatial spectrum of the first simulation with 10 dB SNR.

We consider the estimation bias of the four approaches in the case of SNR = 0 dB~30 dB. Fifty simulations are carried on for each different SNR. The estimation bias in each Monte Carlo simulation is measured by the root-mean-squared error (RMSE) as , where is the estimate of in the th Monte Carlo simulation and is the number of the signals. The estimation bias is the mean of the fifty simulations bias, which is equal to , where is the repetition number of the Monte Carlo simulation. If the approaches fail to resolve the three signals in some simulation repetitions, the of the repetitions are set to be 15°. The RMSE of the first simulation versus SNR is shown in Figure 4. It is shown that the proposed approach has the smallest RMSE. The variance of the first simulation versus SNR is also calculated. The result is shown in Figure 5. It is shown that the proposed approach has smallest variance.

Figure 4: RMSE of the first simulation versus SNR.

Figure 5: Variance of the first simulation versus SNR.

The second simulation is carried out to evaluate the estimation bias and variance of the four approaches in different simulation conditions. Keep the other conditions of the array the same as the first simulation and consider the different vector and the different target information. We consider , and . Two targets have the different elevations which are 3° and 5°. The RMSE of the second simulation versus SNR is shown in Figure 6. It is shown that the proposed approach also has smallest RMSE. The variance of the second simulation versus SNR is shown in Figure 7. It is shown that the proposed approach has smallest variance.

Figure 6: RMSE of the second simulation versus SNR.

Figure 7: Variance of the second simulation versus SNR.

In the third simulation, we keep the other conditions of the array the same as the second simulation and consider the different array antenna number and the different target elevations. The number of array antennas is 32 and two targets have the different elevations which are −6° and −1°. The RMSE of the second simulation versus SNR is shown in Figure 8. It is shown that the proposed approach has the smallest RMSE. The variance of the second simulation versus SNR is shown in Figure 9. It is shown that the proposed approach has smallest variance.

Figure 8: RMSE of the third simulation versus SNR.

Figure 9: Variance of the third simulation versus SNR.

5. Conclusion

A CC-MUISIC algorithm for the mutual coupling correction of the nonuniform L-shaped array with the single snapshot has been proposed. In applying the convex chain optimization, the computational complexity is reduced. Basically, the optimization routine is a process to find the optimization of the interpolated transformation matrix and the mutual coupling matrix for the spatial spectrum estimation of the multisources which are treated as the coherent signal sources for better resolution and estimation. We compared the performance of the CC-MUSIC algorithm with the approaches which combine the MUSIC algorithm based on the spatial smoothing technique with the virtual array technique in the presence of the unknown coupling, the known coupling, and the estimated coupling. The CC-MUSIC algorithm is more robust and more accurate. The simulation results demonstrate the validity of the proposed method in terms of the statistics of bias and variance. Although the proposed algorithm is for the L-shaped nonuniform linear array, it also suits other types of array, such as the circular array and the semicircular array.

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 anonymous reviewers for their constructive comments, which helped improve the presentation of the paper. And this paper is partially supported by the National Natural Science Foundation of China (Grant no. 61301209).

References

Y. Hua, T. K. Sarkar, and D. D. Weiner, “An L-shaped array for estimating 2-D directions of wave arrival,” IEEE Transactions on Antennas and Propagation, vol. 39, no. 2, pp. 143–146, 1991.

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.

B. C. Ng and C. M. Samson, “Sensor-array calibration using a maximum-likelihood approach,” IEEE Transactions on Antennas and Propagation, vol. 44, no. 6, pp. 827–835, 1996.

M. Lin and L. Yang, “Blind calibration and DOA estimation with uniform circular arrays in the presence of mutual coupling,” IEEE Antennas and Wireless Propagation Letters, vol. 5, no. 1, pp. 315–318, 2006.

A. Paulraj and T. Kailath, “Direction of arrival estimation by eigenstructure methods with unknown sensor gain and phase,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '85), vol. 10, pp. 640–643, April 1985.

Z. Ye and C. Liu, “On the resiliency of MUSIC direction finding against antenna sensor coupling,” IEEE Transactions on Antennas and Propagation, vol. 56, no. 2, pp. 3741–3780, 2008.

F. Sellone and A. Serra, “A novel online mutual coupling compensation algorithm for uniform and linear arrays,” IEEE Transactions on Signal Processing, vol. 55, no. 2, pp. 560–573, 2007.

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.

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.

C. Qi, Z. Chen, Y. Zhang, and Y. Wang, “DOA estimation and self-calibration algorithm for multiple subarrays in the presence of mutual coupling,” IEE Proceedings: Radar, Sonar and Navigation, vol. 153, no. 4, pp. 333–337, 2006.

S. Henault, S. K. Podilchak, S. M. Mikki, and Y. M. Antar, “A methodology for mutual coupling estimation and compensation in antennas,” IEEE Transactions on Antennas and Propagation, vol. 61, no. 3, pp. 1119–1131, 2013.

B. H. Wang, H. T. Hui, and M. S. Leong, “Decoupled 2D direction of arrival estimation using compact uniform circular arrays in the presence of elevation-dependent mutual coupling,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 3, pp. 747–755, 2010.

C. H. Niow and H. T. Hui, “Improved noise modeling with mutual coupling in receiving antenna arrays for direction-of-arrival estimation,” IEEE Transactions on Wireless Communications, vol. 11, no. 4, pp. 1616–1621, 2012.

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, pp. 171–176, 2009.

B. Liao, Z.-G. Zhang, and S.-C. Chan, “DOA estimation and tracking of ULAs with mutual coupling,” IEEE Transactions on Aerospace and Electronic Systems, vol. 48, no. 1, pp. 891–905, 2012.

B. H. Wang and H. T. Hui, “Wideband mutual coupling compensation for receiving antenna arrays using the system identification method,” IET Microwaves, Antennas and Propagation, vol. 5, no. 2, pp. 184–191, 2011.

B. Friedlander and A. J. Weiss, “Direction finding using spatial smoothing with interpolated arrays,” IEEE Transactions on Aerospace and Electronic Systems, vol. 28, no. 2, pp. 574–587, 1992.

A. J. Weiss and B. Friedlander, “Performance analysis of spatial smoothing with interpolated arrays,” IEEE Transactions on Signal Processing, vol. 41, no. 5, pp. 1881–1892, 1993.

S. U. Pillai and B. H. Kwon, “Forward/backward spatial smoothing techniques for coherent signal identification,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 1, pp. 8–15, 1989.

R. O. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986.