Joint Space and Time Processing for Unknown Mutual Coupling Blind Calibration and Mixed Sources Identification Using Uniform Circular Array

In classification and localization of mixed far-field and near-field sources, the unknown mutual coupling degrades the performance of most high-resolution algorithms. In practice, the assumption of an ideal receiving sensor array is rarely satisfied. This paper proposes an effective algorithm of mixed sources identification using uniform circular array under unknown mutual coupling. Firstly, according to rank reduction and joint space–time processing, the directions of arrival of far-field sources is estimated directly without mutual coupling elimination. Addition, the joint space–time processing can improve the estimation results in the case of low signal noise ratio of incoming signal sources and small number of snapshots. Then, these estimates are adopted to reconstruct the mutual coupling matrix. Finally, both direction and range parameters of near-field sources are obtained through spatial search after mutual coupling effects and far-field components elimination. The proposed algorithm is described in detail, and its behavior is illustrated by numerical examples.


Introduction
Source localization using the sensor array techniques has received considerable attention over the past decades.For far-field (FF) sources whose wave fronts is plane wave, only the direction of arrival (DOA) parameter is needed to be estimated.A large number of high-resolution algorithms have been proposed to deal with the DOA estimation problem of FF sources in the past decades, such as the estimation of signal parameters via rotational invariance technique (ESPRIT) [1,2], multiple signal classification (MUSIC) method [3,4], and so on.All aforementioned algorithms generally work based on the ideal receiving sensor array assumption without steering vectors mismatch, such as the unknown mutual coupling [5] and non-plane wave front effect [6].
However, in many practical applications, some incoming sources may locate in the Fresnel Region defined as the near-field (NF) of the array, and these sources would be defined as NF sources [7].For an arbitrary NF source, both the direction and range are required to be estimated since the plane wave front assumption is no longer valid.Therefore, the traditional DOA estimation algorithms cannot deal with the localization problem of NF sources.In order to solve the problem of NF sources localization, various algorithms have been developed in the recent decades, such as the reduced rank (RERA) type methods [8][9][10][11], the covariance approximation (CA) type algorithms [12,13], the two dimensional (2D) MUSIC algorithm [5], and the weighted linear prediction method [14].Although all the aforementioned algorithms focus on the pure FF or NF sources scenario, it is more realistic in many applications that FF and NF sources coexist [15], such as electronic surveillance, seismic exploration and speaker localization using microphone arrays.In the mixed NF and FF sources scenario, the above-mentioned algorithms may fail to deal with the mixed sources problem.
Recently, a large number of algorithms have been developed to deal with the problem of mixed sources classification and localization [16][17][18][19][20][21][22].A two-stage MUSIC (TSMUSIC) algorithm has been proposed to solve the mixed sources issue by Liang, which has used the fourth order cumulant (FOC) technique [16].TSMUSIC constructs a special FOC matrix to eliminate the range parameter in the steering vectors with high computational cost.In [17], an oblique projection MUSIC (OPMUSIC) algorithm has been presented by Zhi.Although this algorithm has low computational complexity, it produces a great loss of array aperture.According to the generalized ESPRIT algorithm in [11], a GESPRIT-like algorithm has been presented by Liu [21].However, the GESPRIT algorithm would produce DOA ambiguity when the distance between adjacent sensors is greater than a quarter wavelength.
It must be noticed that all the mixed source algorithms, the inter-space is set as no greater than a quarter wavelength.In the case of small inter-space, the mutual coupling effect is no longer to be ignored.Hence, the performance of all above algorithms would degrade without array calibration.In order to deal with the problem, many mutual coupling modeling methods and DOA estimation algorithms of FF sources have been presented [23][24][25][26].In [23][24][25], various middle sub-array methods are presented to estimate DOAs by setting auxiliary sensors, without calibration.However, these methods would suffer from a great aperture loss.In [26], a method has been presented to deal with the problem of aperture loss and it provided a reasonable simulation result.However, the above-mentioned algorithms only deal with FF sources DOAs.Fortunately, in [27], Xie proposes a method to deal with the mixed sources problem under mutual coupling.
It is worth noting that uniform circular array (UCA) is preferable over uniform linear array (ULA) because of its 360 • azimuthal coverage and almost unchanged directional pattern [28].In [29][30][31][32], various FF sources estimation algorithms with mutual coupling are presented.However, the algorithm in [29] is based on the root-MUSIC with a high computational cost.In [30], this algorithm produces many pseudo peaks in spectrum searching.And in [33,34], two mixed sources localization algorithms are presented.In [33], a TSMUSIC-like method is proposed, which is based on FOC.In [34], a covariance differencing-like (CD-like) algorithm is presented to reduce computational cost.However, both of the two algorithms deal with the problem without the mutual coupling compensation.Unfortunately, there is no algorithm dealing with the mixed sources under mutual coupling in UCA.
In this paper, an effective algorithm is presented to deal with mixed sources classification and localization under unknown mutual coupling.Firstly, according to space-time processing (STP) [35,36], we construct the space-time covariance matrix.And then, we model the mutual coupling matrix (MCM) based on its symmetric Toeplitz property, and a RERA estimating function is constructed to estimate the FF sources DOAs.The MCM can be reconstructed with the FF sources DOA estimates.After the mutual coupling effect elimination, another differencing covariance is formed to eliminate the FF components and noise.Based on the UCA rotational symmetry property, we construct a GESPRITE like function to estimate the NF sources DOAs.Finally, with estimating the DOAs of both NF and FF sources, a multiple MUSIC spectrum search in the DOA and range joint domain would be reduced to the one-dimensional (1D) spectrum search only in the range domain.

Signal Model
Consider that K (FF and NF) narrowband and independent sources impinge on a symmetric UCA, as shown in Figure 1.This uniform circular array is composed of L = 2M (M ∈ N + and M > 1) omni-directional sensors with the radius being R. The authors assume that there are K 1 NF incident sources in the Fresnel Region and the rest K 2 incident sources are FF sources.
origin of the coordinate.At the same time, the FF source comes from (θ k , ϕ k ) and the NF source comes from (θ k , ϕ k , r k ), where θ k ∈[0,2π) denotes the azimuth angle measured counterclockwise from the x-axis and ϕ k ∈[0,π/2] denotes the elevation angle measured downward from the z-axis: where, r k ∈[0.62(D 3 /λ) 1/2 , 2D 2 /λ) denotes the range of NF source, and r k ∈[2D 2 /λ, ∞) is the range of FF source, respectively.Therefore, r k ∈[0.62(D 3 /λ) 1/2 , ∞) denotes the range of mixed sources measured from the UCA center, where D = 2R represents the aperture of the array and λ symbolizes the incident sources wavelength.Consequently, the signal received by the l-th sensor can be modeled as: x t s t e s t e n t l L )sin( ) where s k (t) is the k-th narrowband source, n l (t) is the additive Gaussian noise, and noises in different sensors are independent with the same variance σ n 2 .Here, r l ,k is the distance between the k-th source and the l-th sensor, which has the following form: According to the Taylor series expansion, formula (3) can be approximately given as follows Consequently, formula (1) can be expressed as As a result, it is obvious that the FF source can be seen as a generalized NF source, the range parameter of which increases towards infinite.
According to the discussion in [5], the amplitude of mutual coupling coefficient (MCC) is in inverse proportion to the physical distance between each pair of sensors.In another word, the MCCs between neighboring sensors with a same adjacent space are almost equal to each other.Besides, the MCC would be approximately equal to zero when any two sensors are located enough far away.Consequently, we can model an L × L MCM of the UCA containing P + 1 (P ≤ M−1) nonzero MCCs, and the symmetric Toeplitz matrix C is expressed as: We start with modeling the array received signal vector without unknown mutual coupling.Without loss of generality, all the sensors are located on the xy-plane and the UCA centre is set as the origin of the coordinate.At the same time, the FF source comes from (θ k , ϕ k ) and the NF source comes from (θ k ,ϕ k ,r k ), where θ k ∈[0,2π) denotes the azimuth angle measured counterclockwise from the x-axis and ϕ k ∈[0,π/2] denotes the elevation angle measured downward from the z-axis: where, r k ∈[0.62(D 3 /λ) 1/2 ,2D 2 /λ) denotes the range of NF source, and r k ∈[2D 2 /λ,∞) is the range of FF source, respectively.Therefore, r k ∈[0.62(D 3 /λ) 1/2 ,∞) denotes the range of mixed sources measured from the UCA center, where D = 2R represents the aperture of the array and λ symbolizes the incident sources wavelength.Consequently, the signal received by the l-th sensor can be modeled as: where s k (t) is the k-th narrowband source, n l (t) is the additive Gaussian noise, and noises in different sensors are independent with the same variance σ n 2 .Here, r l,k is the distance between the k-th source and the l-th sensor, which has the following form: According to the Taylor series expansion, Formula (3) can be approximately given as follows Consequently, Formula (1) can be expressed as As a result, it is obvious that the FF source can be seen as a generalized NF source, the range parameter of which increases towards infinite.
According to the discussion in [5], the amplitude of mutual coupling coefficient (MCC) is in inverse proportion to the physical distance between each pair of sensors.In another word, the MCCs between neighboring sensors with a same adjacent space are almost equal to each other.Besides, the MCC would be approximately equal to zero when any two sensors are located enough far away.Consequently, we can model an L × L MCM of the UCA containing P + 1 (P ≤ M − 1) nonzero MCCs, and the symmetric Toeplitz matrix C is expressed as: Therefore, the received signal vector under the unknown mutual coupling can be expressed in a matrix form: where S N (t) and S F (t) denote the NF and FF signal vectors, N(t) signifies L × 1 dimensional complex noise vector, and A N and A F denote the steering vectors of NF source and FF source, respectively as Throughout the paper, the following hypotheses are assumed to hold: 1.
The incoming source signals are statistically independent and zero-mean stationary random process; 2.
The sensor noise is the additive Gaussian one, which is independent from the source signals; 3.
The number of sources K 1 and K 2 are known as prior, and it satisfies K ≤ M.

Proposed Solution to Classification and Localization
In [27], the algorithm shows that the DOA estimation of FF sources is the precondition of MCM reconstruction.Therefore, improving the estimation accuracy of FF sources is helpful to improve the global estimation performance, including NF sources DOA and range estimations.In many practical applications, the algorithm must work against the situations of low signal noise ratio (SNR) and small number of snapshots.Fortunately, for a narrowband incoming source, we can use STP to improve the SNR of incoming sources based on its-self temporal coherence.

Space-Time Covariance Matrix
In practice, it is almost impossible to count an ideal covariance matrix.Hence, in algorithm execution process, covariance matrix must be estimated by a limited N snapshots sample, which is given as: R = We assume that all the incoming signals are narrowband sources, therefore the following holds: where m∈N, and z = exp(j2πf c /f s ) is the time delay factor.Moreover, f s denotes the sampling pulse repetition frequency, and f c denotes the target carrier frequency.Here, we set Q time delay taps and only one sampling pulses between adjacent two time delay taps.The received data vector operated by STP can be written as: where ⊗ is Kronecker product, I Q×Q denotes the Q × Q identity matrix and Z denotes the Q × 1matrix expressed as According to the aforementioned assumptions, the received signal covariance matrix by STP can be calculated by: where I QL×QL is the QL × QL identity matrix, and N and F are diagonal matrices,

FF Sources DOA and MCCs Estimation
By implementing eigenvalue decomposition (EVD) of R ST , the following equation holds: where, Λ s is the diagonal matrix which contains the K largest eigenvalues.At the same time, U n is the QL × (QL − K) eigenvectors matrix which spans the noise subspace of R ST , and U s is the QL × K signal eigenvectors matrix of R ST which spans the signal subspace.As the fact that the MCM is a column full rank Toeplitz matrix, the authors can construct a MUSIC spectrum search function to estimate DOAs and ranges parameters, which is expressed as: From Formula (23), it is obvious that the computational cost of this equation is unbearable with the unknown C.Even if C is accurately estimated as a prior, it requires a 3D spectrum search to estimate and match the DOA and range parameters pair as well.In order to decrease the huge computational cost, the DOA estimation of FF sources must be decoupled from mixed sources estimates and MCM.
Referring to the discussion of mutual coupling problem with (ULA) in [5], Ca(θ, ϕ, r) in UCA could be reformulated as: where c is the vector of P + 1 nonzero MCCs, which expressed as: ) matrices, and they are defined as: where, {•} p,q represents the element corresponding to the p-th row and q-th column of the matrix, [•] p represents the element corresponding to p-th element of steering vector, and mod(p,q) denotes the modulus after dividing p by q.
It is obvious that U s and the combination of (I Q×Q ⊗C)(Z⊗A N ) and (I Q×Q ⊗C)(Z⊗A F ) can span the same signal subspace, and the signal subspace is orthogonal to the noise subspace spanned by U n .Therefore, the following equations hold: Therefore, based on Formula ( 24)-( 29), the FF sources' DOAs can be estimated by the following spectrum search function: where W(θ, ϕ) is defined as: Note that c 0 and W(θ, ϕ) is a Hermite and nonnegative definite matrix.Based on the principle of RERA [29], c H W(θ, ϕ)c would be zero only when W(θ, ϕ) is a singular matrix.In other words, the determinant of W(θ, ϕ) would be equal to zero only when both parameters ϕ and θ are equal to any FF source' DOA expressed by ϕ k and θ k (k = 1 + K 1 , . . ., K).Consequently, the FF sources' DOAs could be estimated accurately by searching the K 2 highest spectrum peaks through the following function: where, det[•] signifies the determinant of a matrix.From equation (32), it is easy seen that p F (θ, ϕ) is independently separated from equation (30).On the other hand, since the smallest eigenvalue of W(θ,ϕ) is also equal to zero when parameters ϕ and θ are equal to ϕ k and θ k (k = 1 + K 1 , . . ., K), the DOAs can also be estimated from the following spatial spectrum searching: where, λ min [•] denotes the smallest eigenvalue of a matrix.Therefore, the computational cost of estimating FF sources' DOAs through (32) or ( 33) is effectively reduced compared with the multi-dimensional MUSIC spectrum search function.It is because of that, that only a 2D spectrum search process is required.
It is noteworthy that the proposed algorithm works in a similar way as that defined in [28,29].However, the algorithm in [28] only solves the problem of pure FF signals, whereas the proposed algorithm aims to deal with the problem of mixed FF and NF signals.Moreover, the algorithm in [29] solves the mixed sources problem using a ULA.However, our work explicitly addresses the mixed sources problem under mutual coupling effect, using a UCA.
Substituting the estimated ϕ k and θ k (k = 1 + K 1 , . . .,K) and then perform eigendecomposition of W(θ k ,ϕ k ), we can obtain the smallest eigenvalue λ min and the corresponding eigenvector v min , where v min has been normalized by its first element.Finally, we can obtain the non-zero MCCs by: According to the symmetric Toeplitz structure modeled as Formula ( 6), the MCM can be reconstructed after all MCCs have been calculated.Therefore, the reconstructed MCM could be used to eliminate the mutual coupling effects in the following signal processing process.

NF Sources DOA Estimation
This section begins with the discussion of MCM.It is noteworthy that the MCM of an UCA is different to that in an ULA, and it has some special properties based on its special array structure.

Remark 1:
The MCM C is a symmetric Toeplitz matrix, which holds: and its inverse matrix is also a symmetric Toeplitz matrix, which holds: Remark 2: With C being a complex symmetric Toeplitz matrix, CC H is a real symmetric Toeplitz matrix (see Appendix A), which holds: and the inverse matrix of CC H is also a real symmetric Toeplitz matrix, which holds: As known, the estimation accuracy of NF sources will be improved as the SNR increases.If we can eliminate the noise, the SNR of NF will be approximately equal to infinity.Fortunately, the special properties of the MCM and array structure in UCA can be used to eliminate the noise and FF components effectively.Consequently, we can deal with the NF sources estimation problem only in a space domain with a high enough SNR.
According to Formula ( 17) and ( 19), the received signal covariance matrix only in space domain is expressed as: After estimating and reconstructing mutual coupling matrix, the mutual coupling effects can be eliminated effectively.Therefore, we can get: Based on the symmetric property of UCA configuration with an even number of sensors, the following holds: Due to the fact that R F is a Hermitian matrix but R N only holds a Hermitian structure, we can obtain the following: where J is the L × L particular transformational matrix which can be defined as: where, I M×M denotes the M × M identity matrix.Consequently, the property differences between R F and R N can be utilized to implement the FF sources components elimination in the covariance matrix.
Then, the differencing matrix can be expressed as: where, A D is the virtual steering vectors expressed as: Based on ( 13), ( 41) and (43), the following holds (see Appendix B) where denotes Hadamard product, and D is a rotation factor matrix defined as: And then, by implementing the EVD of R D , the following holds: where, ∆ s is a diagonal matrix which contains 2K 1 non-zero eigenvalues, and γ 0 is the zero eigenvalue, V s is the L × 2K 1 eigenvectors matrix spanning the signal subspace of R D , and V n is the L × (L − 2K 1 ) eigenvectors matrix spanning the noise subspace.It is obvious that the following holds: Consequently, we can construct a diagonal matrix: and the NF sources' DOAs could be estimated accurately by searching the K 1 highest spectrum peaks through the following function: where Q denotes an arbitrary L × 2K 1 full column rank matrix.The typical value of Q is set as that composed of arbitrary 2K 1 columns of the L × L identity matrix.Conversely, since Q H {JV s −Φ(θ,ϕ)V s } contains zero eigenvalue only when ϕ and θ are equal to ϕ k and θ k (k = 1, . . ., K 1 ), the DOAs can also be estimated from the following spatial spectrum searching: In this subsection, the covariance differencing operation in UCA eliminate the FF components and noise effectively based on the structure differences between the NF and FF sources covariance matrices.As a conclusion, the proposed algorithm successfully avoids the 3D spectrum searching and sources classification processing with the help of covariance differencing.

NF Sources Range Estimation
After estimating the NF sources' DOAs and replacing the unknown DOA variable pairs in the following equation by the estimated ones, the k-th NF source' range parameter can be estimated by: From the above equation, it shows that the estimation results of both DOA and range parameters would be paired together automatically without any parameters matching processing.

The Proposed Algorithm Summary
In order to show the specific execution process, the proposed algorithm can be summarized as following.STEP 1. Estimate the STP covariance matrix through Formula (19).STEP 2. Eigen decompose the STP covariance matrix to generate its noise subspace.STEP 3. Estimate the FF sources' DOAs through the RERA estimator by Formula (32) or (33).STEP 4. Obtain the MCCs by Formula (34), and construct the MCM by Formula (6).STEP 5. Compensate the MCM and by Formula (40).STEP 6. Construct the differencing matrix by Formula (45).STEP 7. Eigen decompose the differencing matrix to generate its noise subspace.STEP 8. Estimate the DOAs of NF sources through Formula (52) or (53).STEP 9. Estimate the range parameter of NF sources through Formula (54).

Number of Required Array Elements
In order to get the unique solution of (32) or (33), W(θ,ϕ) must be invertible at non-incident direction.Due to the fact that all the incident sources have the same carrier frequency, the extended freedom degree by time delay is invalid.Furthermore, the maximum rank of Z⊗B(θ, ϕ, ∞) maybe M, and the space domain freedom degree of U n U n H is L−K, they must satisfy: consequently, the parameter identifiability condition of the algorithm is K ≤ M.Moreover, in order to estimate the MCCs, [W(θ,ϕ)] must be equal to zero at least once.In the other word, at least one FF source is required to be used for MCCs estimation.

Contribution of STP
By eigenvalue decomposition, the covariance matrix R ST can be decomposed as shown in ( 22), here Λ s is the diagonal matrix which contains the K largest eigenvalues: and the eigenvalues hold the following comparison expression: In fact, eigenvalue composition is the unitized orthogonalization processing of steering vectors.Therefore, each eigenvalue is a linear combination of all signals and noise power, and the i-th eigenvalue can be expressed as: where α i,k ∈(0,1) is the power factor whose value depends on the projection of steering vectors and relative power ratios of all sources, and s k 2 denotes the power of the k-th source.Compared with pure space domain processing, the eigenvalues are almost enlarged Q times.Thus, we define the power gain by STP as: Consequently, it is equivalent to improving the SNR of incoming signals.

Computational Complexity
In order to compare computational complexity intuitively, TSMUSIC-like and CD-like and algorithms are set as comparison objects of the proposed algorithm.For all three algorithms, the major computation processes are calculating statistical matrices, eigenvalue decomposing and spectrum searching.It is defined that the search step of azimuth θ∈(0 • ,360 • ] is θ ∆ , the search step of elevation ϕ ∈(0 • ,90 • ] is ϕ ∆ , and the search step of r∈[0.62(D 3 /λ) 1/2 ,2D 2 /λ] is r ∆ for K 1 NF sources.We assume that N is the number of snapshots.
For TSMUSIC-like algorithm, the major computations are to compute an L × L FOC matrix and an L × L second-order cumulant one, and to implement the EVDs of the two matrices for spatial searching twice.For CD-like algorithm the major computations are to form an L × L covariance matrix and an L × L noise subspace matrix, and to implement once EVD of the covariance matrix.After that, three times spatial searching are required to estimate the NF sources' DOAs, the ranges of NF sources and the FF sources' DOAs.The proposed algorithm constructs a second-order cumulant space-time matrix.The major computations are to implement the EVDs of the QL × QL covariance matrix and the L × L differencing matrix, and to perform three times spatial search for estimating DOA and range.Consequently, the computational complexity of the three algorithms is shown in Table 1.

Algorithms
Statistical Matrices EVDs Spectral Search However, in order to ensure sufficient measuring precision, the search steps are usually set as very small values.Therefore, the computational burden of spectral search is huge, which is usually more than 90 percent of the total computational cost.In this case, the computational cost of proposed algorithm is approximately equal to that of CD-like, and almost twice as that of TSMUSIC-like.

Capacity for Mixed Sources Classification
From the above analysis, it is obvious that TSMUSIC-like algorithm cannot classify the mixed sources autonomously.Especially if the FF sources and the NF sources have the same DOA, the TSMUSIC-like algorithm would be invalid.Fortunately, the CD-like algorithm and proposed algorithm can eliminate FF component when estimating the NF sources.However, the CD-like algorithm ignores the mutual coupling effect.

DOA Ambiguity
In order to avoid DOA ambiguity in UCA, we must construct Φ(θ,ϕ) without ambiguity, i.e., the steering vector d(θ,ϕ) must be unique at different incident-directions.Therefore, the aperture of UCA dealing with mixed sources problem must be set as a half of that, solving pure FF sources estimation without ambiguity.For the similar question in ULA, the inter-space is usually set as a quarter of wavelength, half of the normal inter-space.

Simulation Results
Some simulations are conducted in this section to evaluate the proposed algorithm.A UCA of 10 elements with radius R = λ is taken into consideration.The input signal to noise ratio (SNR) of the k-th source is defined as 10 × log 10 (σ k 2 /σ n 2 ), where σ k 2 denotes the power of the k-th source, and σ n 2 denotes the noise power.We assume that all the sources are with equal power.In the following experiments, the performance is measured by the root mean square error (RMSE) of 200 independent Monte Carlo experiments.It is noted that the estimation performance of the range parameters is only for the NF sources experiment.Hence, we do not give the estimation performance of the range parameters of FF sources.

Spectrums of DOA and Range Estimation
In the first simulation, the authors consider a scenario in where two NF sources and two FF sources coexist, and they are located at {θ 1 = 43 From Figure 2a, it shows that TSMUSIC-like generates four highest sharp spectrum peaks on the directions of NF and FF sources.However, there are other false peaks caused by the unknown mutual coupling, which would lead to false estimation.Compared with the DOA spectrum after MCM elimination, it obviously shows that this algorithm generates serious error in DOA estimation.Assuming that NF and FF sources have been classified, the range spectrum of NF sources by TSMUSIC-like is shown in Figure 2b.Similarly, this algorithm generates unreliable results in range estimation, without MCM elimination.
From Figure 2a, it shows that TSMUSIC-like generates four highest sharp spectrum peaks on the directions of NF and FF sources.However, there are other false peaks caused by the unknown mutual coupling, which would lead to false estimation.Compared with the DOA spectrum after MCM elimination, it obviously shows that this algorithm generates serious error in DOA estimation.Assuming that NF and FF sources have been classified, the range spectrum of NF sources by TSMUSIC-like is shown in Figure 2b.Similarly, this algorithm generates unreliable results in range estimation, without MCM elimination.In Figure 3a, it shows that CD-like generates worse spectrums of NF and FF sources than that in Figure 2a.This is because that the structural properties of array steering vectors are invalid under unknown mutual coupling.In particular, there may be fatal errors in DOA estimation of NF sources, as shown in detail.In Figure 3b, due to the propagation error and model mismatching, this algorithm generates serious error in range estimation without MCM elimination.In Figure 3a, it shows that CD-like generates worse spectrums of NF and FF sources than that in Figure 2a.This is because that the structural properties of array steering vectors are invalid under unknown mutual coupling.In particular, there may be fatal errors in DOA estimation of NF sources, as shown in detail.In Figure 3b, due to the propagation error and model mismatching, this algorithm generates serious error in range estimation without MCM elimination.In Figure 4a, it shows that the proposed algorithm generates two sharp peaks in spectrum of NF sources and another two sharp peaks in that of FF sources, corresponding to the actual DOAs of mixed sources.The range estimation is shown in Figure 4b, it shows that the proposed algorithm can generate reliable range estimation results as well.In Figure 4a, it shows that the proposed algorithm generates two sharp peaks in spectrum of NF sources and another two sharp peaks in that of FF sources, corresponding to the actual DOAs of mixed sources.The range estimation is shown in Figure 4b, it shows that the proposed algorithm can generate reliable range estimation results as well.

CD-like.
In Figure 4a, it shows that the proposed algorithm generates two sharp peaks in spectrum of NF sources and another two sharp peaks in that of FF sources, corresponding to the actual DOAs of mixed sources.The range estimation is shown in Figure 4b, it shows that the proposed algorithm can generate reliable range estimation results as well.

Performance versus Snapshot Number
In the second simulation, the authors consider a scenario in where one NF source and one FF sources coexist, and the location parameters are {θ 5 = 77°, ϕ 5 = 47°, r 5 = 2.2λ} and {θ 6 = 259°, ϕ 6 = 58°, r 6 = ∞}, with the SNR being set as 10 dB.And the time delay tap number of the proposed algorithm is set as Q = 3.The nonzero MCCs are set as the same as those in the first experiment and the number of snapshots varies from 100 to 1000.The RMSEs of the DOA and range estimates, as the snapshot number changes, are shown in Figure 5 and Figure 6.It obviously shows that both the DOA and the range estimation RMSEs of the proposed algorithm decrease monotonically as the snapshot number increases.It is due to the fact that a larger sampling number will produce better estimate of the covariance matrix for stationary data.Compared with the proposed algorithm, the DOA RMSEs of the other two algorithms no longer decrease as the number of snapshots increases, shown in Figure 5a.It is because that steering mismatching caused by mutual coupling effect.Moreover, as shown in Figure 5b, TSMUSIC-like estimates the DOAs of mixed sources in the same search, hence the mixed sources have the similar

Performance versus Snapshot Number
In the second simulation, the authors consider a scenario in where one NF source and one FF sources coexist, and the location parameters are {θ 5 = 77 • , ϕ 5 = 47 • , r 5 = 2.2λ} and {θ 6 = 259 • , ϕ 6 = 58 • , r 6 = ∞}, with the SNR being set as 10 dB.And the time delay tap number of the proposed algorithm is set as Q = 3.The nonzero MCCs are set as the same as those in the first experiment and the number of snapshots varies from 100 to 1000.
The RMSEs of the DOA and range estimates, as the snapshot number changes, are shown in Figures 5 and 6.It obviously shows that both the DOA and the range estimation RMSEs of the proposed algorithm decrease monotonically as the snapshot number increases.It is due to the fact that a larger sampling number will produce better estimate of the covariance matrix for stationary data.Compared with the proposed algorithm, the DOA RMSEs of the other two algorithms no longer decrease as the number of snapshots increases, shown in Figure 5a.It is because that steering mismatching caused by mutual coupling effect.Moreover, as shown in Figure 5b, TSMUSIC-like estimates the DOAs of mixed sources in the same search, hence the mixed sources have the similar RMSEs with known MCM.On the other hand, CD-like performances better results of NF sources than that of FF sources, due to the covariance differencing operation with known MCM.As known, in proposed algorithm, FF source DOA estimation is the precondition of NF source DOA estimation, hence the propagation error will degrade NF DOA estimation.In addition, the proposed algorithm produces the best performance because of STP.A similar result is also obtained in range estimation as shown in Figure 6.RMSEs with known MCM.On the other hand, CD-like performances better results of NF sources than that of FF sources, due to the covariance differencing operation with known MCM.As known, in proposed algorithm, FF source DOA estimation is the precondition of NF source DOA estimation, hence the propagation error will degrade NF DOA estimation.In addition, the proposed algorithm produces the best performance because of STP.A similar result is also obtained in range estimation as shown in Figure 6.

Performance versus SNR
In the third simulation, almost all simulation conditions are adopted as the same as those in the second experiment, except that the number of snapshots is set as 400, and the SNR varies from −5 dB to 20 dB.
The RMSEs of the DOA and range estimates, as the changes of SNR, are shown from Figure 7 and Figure 8.Compared with the proposed algorithm, both the DOA and range RMSEs of the other two algorithms no longer decrease as the SNR increases due to mutual coupling effect.It shows that only the proposed algorithm can provide satisfactory DOA estimation performance of mixed sources in Figure 7a.Moreover, in Figure 7b, FF source DOA RMSEs by TSMUSIC-like are similar to NF ones by itself with known MCM.Respectively, CD-like produces better results of NF sources and worse results of FF sources than those by TSMUSIC-like.Additionally, in the proposed algorithm, the propagation error will degrade NF DOA estimation.As a result, the proposed algorithm outperforms the other two algorithms in both DOA and range estimation shown in Figure 8, because of STP.

Performance versus SNR
In the third simulation, almost all simulation conditions are adopted as the same as those in the second experiment, except that the number of snapshots is set as 400, and the SNR varies from −5 dB to 20 dB.
The RMSEs of the DOA and range estimates, as the changes of SNR, are shown from Figures 7  and 8. Compared with the proposed algorithm, both the DOA and range RMSEs of the other two algorithms no longer decrease as the SNR increases due to mutual coupling effect.It shows that only the proposed algorithm can provide satisfactory DOA estimation performance of mixed sources in Figure 7a.Moreover, in Figure 7b, FF source DOA RMSEs by TSMUSIC-like are similar to NF ones by itself with known MCM.Respectively, CD-like produces better results of NF sources and worse results of FF sources than those by TSMUSIC-like.Additionally, in the proposed algorithm, the propagation error will degrade NF DOA estimation.As a result, the proposed algorithm outperforms the other two algorithms in both DOA and range estimation shown in Figure 8, because of STP.

Conclusions
In this paper, a high performance joint space-time classification and localization algorithm for mixed sources in the presence of unknown mutual coupling in UCA is proposed.Compared with the aforementioned algorithms, the proposed algorithm is effective in mixed sources classification and localization, and it can provide a better DOA and range estimation performance under the mutual coupling effect.In addition, for computational complexity, the proposed algorithm successfully avoids multi-dimensional spectrum search (including MCCs and range), parameter matching and high-order cumulant.Finally, many simulation results give forceful proof that the proposed algorithm is efficient for mixed sources classification and localization under unknown mutual coupling in UCA, especially in the case of low SNR and small snapshot number.

Appendix A
Without loss of generality, we replace the zero coefficients with unknown complex ones, and C is an L × L symmetric Toeplitz matrix.We defined that:

Conclusions
In this paper, a high performance joint space-time classification and localization algorithm for mixed sources in the presence of unknown mutual coupling in UCA is proposed.Compared with the aforementioned algorithms, the proposed algorithm is effective in mixed sources classification and localization, and it can provide a better DOA and range estimation performance under the mutual coupling effect.In addition, for computational complexity, the proposed algorithm successfully avoids multi-dimensional spectrum search (including MCCs and range), parameter matching and high-order cumulant.Finally, many simulation results give forceful proof that the proposed algorithm is efficient for mixed sources classification and localization under unknown mutual coupling in UCA, especially in the case of low SNR and small snapshot number.
it is clear that the elements in R C are expressed as: and it clearly shows that there are P pairs of conjugate coefficients and a real one in each element, and the elements on each minor diagonal are same.When L is even, i.e., L = 2P, the MCM is constructed as: Therefore, the elements in R C are expressed as:

Figure 1 .
Figure 1.The uniform circular array configuration.

Figure 1 .
Figure 1.The uniform circular array configuration.

Figure 5 .
Figure 5. (a) Root mean square error (RMSEs) of the direction of arrival (DOA) estimates versus snapshot number; (b) RMSEs of the DOA estimates versus snapshot number with known mutual coupling matrix (MCM).

Figure 5 .
Figure 5. (a) Root mean square error (RMSEs) of the direction of arrival (DOA) estimates versus snapshot number; (b) RMSEs of the DOA estimates versus snapshot number with known mutual coupling matrix (MCM).

Figure 5 .
Figure 5. (a) Root mean square error (RMSEs) of the direction of arrival (DOA) estimates versus snapshot number; (b) RMSEs of the DOA estimates versus snapshot number with known mutual coupling matrix (MCM).

Figure 6 .
Figure 6.RMSEs of the range estimates versus snapshot number.

Figure 6 .
Figure 6.RMSEs of the range estimates versus snapshot number.

Figure 7 .
Figure 7. (a) RMSEs of the DOA estimates versus signal noise ratio (SNR); (b) RMSEs of the DOA estimates versus SNR with known MCM.

Figure 8 .
Figure 8. RMSEs of the range estimates versus SNR.
DOA of FFsource TSMUSIC-like DOA of NFsource CD-like DOA of FFsource CD-like DOA of NFsource Proposed DOA of FFsource Proposed DOA range of NFsource TSMUSIC-like range of NFsource with known MCM CD-like range of NFsource CD-like range of NFsource with known MCM Proposed range of NFsource

Figure 7 .Figure 7 .
Figure 7. (a) RMSEs of the DOA estimates versus signal noise ratio (SNR); (b) RMSEs of the DOA estimates versus SNR with known MCM.

Figure 8 .
Figure 8. RMSEs of the range estimates versus SNR.

Figure 8 .
Figure 8. RMSEs of the range estimates versus SNR.