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 x-axis and z-axis, respectively. Then a new matrix is constructed with the proper entries of cross covariance matrix of the chosen elements outputs on x-axis and z-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 x-axis and z-axis, respectively. The ULAs share the middle element as reference at origin of coordinates. Each ULA has 2M+1sensors with element spacing of d, which means that there exist 4M+1ones in the whole array because of the sharing element. Assume thatKnarrowband, far-field, noncoherent signals with the same wave length λ (λ=2d) impinge on this array, and the kth source’s elevation and azimuth are αk and βk, respectively. The array configuration is shown in Figure 1.
The cross array configuration.
Let C and t denote the mutual coupling matrix (MCM) and the time index, respectively. The output of the array can be expressed as
(1)Y(t)=CAS(t)+N(t),
where Y, A, S, and N 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 σ2) vector, respectively. They can be described as (the time indextis omitted for brevity)
(2)Y=[x1,…,xM,…,x2M+1,z1,…,zM,…,z2M+1]T,A=[a(α1,β1),…,a(αk,βk),…,a(αK,βK)],a(αk,βk)=[ej2πdMcosαkcosβk/λ,…,1,…,e-j2πdMcosαkcosβk/λ,ej2πdMsinαk/λ,…,1,…,e-j2πdMsinαk/λ]T,S=[s1,…,sk,…,sK]T,N=[nx1,…,nxM,…,nx2M+1,nz1,…,nzM,…,nz2M+1]T,
where superscript T denotes transpose.
As shown in [3], the MCM of ULA can be expressed as a banded complex symmetric Toeplitz matrixB(3)B=Toeplitz{1,b1,…,bp},
where integer p 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 exceedspd.
Based on the MCM model of ULA by (3), the MCM of cross array has the following form:
(4)C=[BDDB],
where B is the MCM of ULA on x- or z-axis, and D describes the mutual coupling between the two ULAs. The matrix D can be presented by
(5)D=[JD1Jh1TJD1h11h1JD1J(h1J)TD1],
where
(6)h1=[0(M-p)×1bp⋯b1],J=[00⋯100⋰⋮⋮1⋱01⋯00]M×M.D1=(dij)M×M,dijdenotes the coefficient of mutual coupling between the two elements spacing i2+j2d, and dij satisfies
(7)dij=dji,dij={≠0i2+j2<p=0i2+j2>p.
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 Y~′ and Y~′′ be the outputs of elements affected by the same mutual coupling onx-axis andz-axis, respectively.
Y~′ can be written as
(8)Y~′=HY=∑k=1K[bpukM+⋯+ukM-p+⋯bpukM-2p⋮bpuk2p+⋯+ukp+⋯bpbp+⋯+uk-p+⋯bpuk-2p⋮bpuk2p-M+⋯+ukp-M+⋯bpuk-M]sk+N~1=∑k=1Kfkska~k+N~1=A~1T1S+N~1,
where fk=bpukM+⋯+ukM-p+⋯bpukM-2p,
(9)H=[H10︸2(M-2p+1)×(2M+1)],H1=[0︸(M-2p+1)×pI0︸(M-2p+1)×(2p-1)00︸(M-2p+1)×p000I0],a~k=[1,…,uk-(M-2p),uk-M,…,uk-2M+2p]T,uk=exp(j2πdcosαkcosβkλ),A~1=[a~1,…,a~k,…,a~K],T1=diag{f1,…,fk,…,fK}.
There are no entries of D in (8) as the chosen elements on x-axis are not affected by the sensors onz-axis, because the distances between them are out of pd.
Similarly, Y~′′ can be expressed as
(10)Y~′′=A~2T2S+N~2,
where A~2=[b~1,…,b~k,…,b~K],
(11)T2=diag{g1,…,gk,…,gK},(12)b~k=[1,…,vk-(M-2p),vk-M,…,vk-2M+2p]T,(13)gk=bpvkM+⋯+vkM-p+⋯bpvkM-2p,(14)vk=exp(j2πdsinαkλ).
In order to introduce rotational invariance techniques for ULA to this array, Y~′ is partitioned into two parts, so is Y~′′. They are denoted by Y~11, Y~12, Y~21, and Y~22, respectively:
(15)Y~11=[Y~′(1:M-2p)Y~′(M-2p+2:2M-4p+1)],Y~12=[Y~′(2:M-2p+1)Y~′(M-2p+3:2M-4p+2)],Y~21=[Y~′′(1:M-2p)Y~′′(M-2p+2:2M-4p+1)],Y~22=[Y~′′(2:M-2p+1)Y~′′(M-2p+3:2M-4p+2)].
Y~11 and Y~21 are the submatrixes of Y~′ and Y~′′, respectively, so the cross covariance matrix of Y~11 and Y~21 can be built with the proper entries of cross covariance matrix of Y~′ and Y~′′:
(16)R~1=A~11R~ssA~21H=[R~11R~12R~13R~14],
where A~11 and A~21 are the ideal steering matrixes of arrays corresponding to Y~11 and Y~21, respectively,
(17)R~=E{Y~′Y~′′H},(18)R~11=R~(1:M-2p,1:M-2p),(19)R~12=R~(1:M-2p,M-2p+2:2M-4p+1),(20)R~13=R~(M-2p+2:2M-4p+1,1:M-2p),(21)R~14=R~(M-2p+2:2M-4p+1,M-2p+2:2M-4p+1),(22)R~ss=T1E{SSH}T2H.
Similarly, let R~2, R~3, and R~4, which can be constructed with the proper entries chosen from R~ as R~1 in (16), be the cross covariance matrixes of Y~12 and Y~21, Y~11, and Y~22, and Y~12 and Y~22, respectively,
(23)R~2=A~12R~ssA~21H=A~11𝚽1R~ssA~21H,(24)R~3=A~11R~ssA~22H=A~11R~ss(A~21𝚽2)H=A~11𝚽2HR~ssA~21H,(25)R~4=A~12R~ssA~22H=A~11𝚽1𝚽2HR~ssA~21H,
where A~12 and A~22 are the ideal steering matrixes of arrays corresponding to Y~12 and Y~22, respectively,
(26)𝚽1=diag{u1-1,…,uk-1,…,uK-1},𝚽2=diag{v1-1,…,vk-1,…,vK-1}.
The new matrix R~′ consisting of the four cross covariance matrixes includes the rotational invariance factors,
(27)R~′=[R~1R~2R~3R~4]=[A~11A~11𝚽1A~11𝚽2HA~11𝚽1𝚽2H]R~ssA~21H.
Similar to the conventional propagator method, R~′ can be partitioned as
(28)R~′=[GTQT]T,
where G and Q are the first K and the last 8(M-2p)-K rows, respectively.
A matrix P (the propagator parameter) is defined to satisfy
(29)P=minP∥Q-PG∥2.
The solution of (29) can be obtained as follows:
(30)P=QG+=QGH(GGH)-1.
With the relationship between the propagator parameter and ideal steering matrix, we obtain an equation as
(31)P′=[IP]=[P1P2P3P4]=[A~11A~11𝚽1A~11𝚽2HA~11𝚽1𝚽2H]T3,
where Pi (i=1,2,3,4) is the (2M-4p)×K matrix and T3is the K×K matrix with rank K.
Similar to ESPRIT, from (31) we can have the equations correlated with azimuths and elevations as
(32)P1+P2=T3-1𝚽1T3,P3+P4=T3-1𝚽1T3,P1+P3=T3-1𝚽2HT3,P2+P4=T3-1𝚽2HT3.
To solve the above equation sets and pair the azimuths and elevations, four times eigen-decomposing and 3[K(K+1)/2-1] 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 P1+P2, and the eigenvalues and the eigenvectors are denoted by λ1k and ρk (k=1,2,…,K), respectively Then substitute ρk into other equations, as the first one of the right equation sets in (32); the corresponding eigenvalue λ3k can be obtained from
(33)P1+P3ρk=λ3kρk.
Let P1+P3ρk=ωk, and λ3k can be expressed by
(34)λ3k=1K∑i=1Kωkiρki,
where ωki and ρki are the ith element of ωk and ρk, respectively.
Assume that λ2k and λ4k denote the corresponding eigenvalues of P3+P4, and P2+P4 respectively, and they can be calculated by (34).
Using the relationship among λik (i=1,2,3,4), uk, vk, αk, and βk, the elevation αk and the azimuth βk of the kth source can be represented as
(35)αk={arcsin{arg(λ3k)/π}+arcsin{arg(λ4k)/π}}2,βk={arccos{-arg(λ1k)/πcosαk}+arccos{-arg(λ2k)/πcosαk}}2.
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 x-axis and z-axis, eigen-decomposing a matrix with rank K, and pairing the parameters, so the computation load is O((4M+2)3)+N(2M+1)2+O(4K3)+3K2[K(K+1)/2-1](Nis the snapshots number); the computation load on the proposed method is focused on cross matrix of output of elements on x-axis and z-axis, once eigen-decomposing a matrix with rank K, and pairing the parameters, so the computation load is N(2M-4p+2)2+O(K3)+3K(K2+K).
Finally, the necessary conditions for the proposed method are given as follows:
for unique estimation, 2(M-2p)>K is required;
the sources can not be located at the blind angles which make the determinant of T1 or T2 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 (M=13). The mutual coupling degree of freedom p is two and the coefficients of mutual coupling are denoted by b1=0.3065+0.3951j, b2=0.1018-0.1721j, and d11=0.2154-0.2738j.
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.
2-D spatial spectrum of the PM algorithm in the presence of unknown mutual coupling.
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.
2-D spatial spectrum of the PM algorithm in the presence of known mutual coupling.
2-D spatial spectrum of the proposed method.
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:
(36)RMSE=12∑k=12E{(α^k-αk)2+(β^k-βk)2}.
RMSE versus SNR.
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.
RMSE versus number of snapshots.
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.
RollerC.WasylkiwskyjW.Effects of mutual coupling on super resolution DF in linear arrays5Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal ProcessingMarch 1992257260WeissA. J.FriedlanderB.Mutual coupling effects on phase-only direction finding19924055355412-s2.0-002686972310.1109/8.142628FriedlanderB.WeissA. J.Direction finding in the presence of mutual coupling19913932732842-s2.0-0026117417SelloneF.SerraA.A novel online mutual coupling compensation algorithm for uniform and linear arrays20075525605732-s2.0-3384761109610.1109/TSP.2006.885732WangB.-H.WangY.-L.ChenH.A robust DOA estimation algorithm for uniform linear array in the presence of mutual coupling3Proceedings of the IEEE Antennas and Propagation Society International SymposiumJune 20039249272-s2.0-0041941147LinM.YangL.Blind calibration and DOA estimation with uniform circular arrays in the presence of mutual coupling2006513153182-s2.0-3454714494110.1109/LAWP.2006.878898HuX.ChenH.WangY.ChenJ.A self-calibration algorithm for cross array in the presence of mutual coupling20115448368482-s2.0-7995547364110.1007/s11432-010-4144-0YeZ.DaiJ.XuX.WuX.DOA estimation for uniform linear array with mutual coupling20094512802882-s2.0-6544916803410.1109/TAES.2009.4805279YeZ.LiuC.2-D DOA estimation in the presence of mutual coupling20085610315031582-s2.0-5424913868410.1109/TAP.2008.929446LiuC.YeZ.ZhangY.Autocalibration algorithm for mutual coupling of planar array20109037847942-s2.0-7035073348010.1016/j.sigpro.2009.08.014