Next Article in Journal
Theoretical and Experimental Investigations of a Pseudo-Magnetic Levitation System for Energy Harvesting
Next Article in Special Issue
Robust Null Broadening Beamforming Based on Covariance Matrix Reconstruction via Virtual Interference Sources
Previous Article in Journal
Biofeedback Systems for Gait Rehabilitation of Individuals with Lower-Limb Amputation: A Systematic Review
Previous Article in Special Issue
DOA Estimation Using Fourth-Order Cumulants in Nested Arrays with Structured Imperfections
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Array Switching Strategy for Direction of Arrival Estimation with Coprime Linear Array in the Presence of Mutual Coupling

1
College of Electronic and Information Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China
2
Key Laboratory of Dynamic Cognitive System of Electromagnetic Spectrum Space, Ministry of Industry and Information Technology, Nanjing 211106, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(6), 1629; https://doi.org/10.3390/s20061629
Submission received: 16 February 2020 / Revised: 9 March 2020 / Accepted: 10 March 2020 / Published: 14 March 2020

Abstract

:
While the coprime array still suffers from performance degradation due to the mutual coupling dominated by the interleaved subarrays, we propose an array switching strategy for coprime linear array (CLA) by utilizing the large inter-element spacings of the subarrays to mitigate the mutual coupling. Specifically, we first collect the signals by separately activating the two subarrays, where the severe mutual coupling effect is significantly reduced. As a result, well-performed initial direction of arrival (DOA) estimates can be achieved. Subsequently, we establish a quadratic optimization problem by reconstructing the contaminated steering vector of the total CLA elaborately to calculate the mutual coupling coefficients with the initial DOA estimates. Finally, we can obtain refined DOA estimates by an iteration procedure based on the estimated mutual coupling matrix. In addition, numerical simulations are provided to demonstrate the merits of the proposed scheme.

1. Introduction

As a fundamental part of array signal processing, direction of arrival (DOA) estimation plays a critical role in various fields, e.g., radar, acoustics, cognitive radio and wireless communication [1,2,3,4,5,6,7]. Numerous high-resolution DOA estimation algorithms are developed for uniform linear arrays (ULAs), such as MUSIC [8], ESPRIT [9] and the variants [10,11]. However, the compact array structure will lead to limited array aperture and serious mutual coupling effect.
In recent years, coprime array [12,13,14,15,16] has received considerable attention due to the promising merits, such as the enhanced resolution, the increased degrees of freedom (DOFs) and the reduced mutual coupling. Specifically, the DOA estimation methods proposed for coprime linear array (CLA) can be divided into two categories in general. The first type is to decompose the CLA into two uniform subarrays and then eliminate the phase ambiguity problem based on the arrangements with coprime property [17,18,19]. The other one is to construct an equivalent received signal collected from an augmented virtual array to avoid spatial aliasing and enhance the DOF [20,21,22,23]. However, the aforementioned researches critically depend on the array manifold that is not distorted by impairments such as mutual coupling [24]. Since the two subarrays in the prototype CLA are interleaved, the mutual coupling cannot be neglected, which can significantly degrade the performance. To tackle this problem, an unfolded CLA was designed in [25], where the two subarrays are configured along the positive and negative axis respectively and thereby the minimum adjacent distance is enlarged. The coprime array with displaced subarrays (CADiS) in [26] introduces a displacement between the two subarrays to enlarge the smallest inter-element spacing, which can prominently alleviate the mutual coupling. Particularly, an extended coprime array structure in [27] relocates a proper number of sensors in the CADiS, where the mutual coupling is further reduced. In [28], the thinned coprime array exploits the redundancy in the difference co-array to achieve increased DOFs, where the mutual coupling is decreased simultaneously, owing to the reduced number of sensor pairs with small separations. However, since at least one sensor pair is needed to generate the virtual sensors, the sensor pairs with small separations still lead to severe mutual coupling effect. Although significant efforts have been put into mutual coupling reduction techniques for CLA, the investigations about mutual coupling calibration for CLA are neglected.
As alternative approaches to deal with mutual coupling, a series of mutual coupling calibration algorithms are proposed [29,30,31,32,33,34,35,36,37], and a few of them are applicable to nonuniform linear arrays. Specifically, the active calibration [35] utilizes auxiliary sources with exactly known locations to perform calibration, whereas it is time-consuming and impractical in real systems. Moreover, as the mutual coupling effect varies with the environment, repeated off-line calibrations are not conducive to the engineering applications. An iteration-based calibration method was proposed in [36] for nonuniform linear arrays, but the solution may fall into local optimum and the convergence cannot be guaranteed theoretically. Besides, two mutual coupling compensation methods for nonuniform linear arrays were proposed in [37]. The first one is an iteration method which is based on the precisely known mutual impedance matrix obtained by electromagnetic simulation while the second method is proposed by employing a global optimization algorithm, but suffering from heavy computational burden and the loss of the structural characteristics of CLA.
In this paper, the issue of DOA estimation for the prototype CLA with unknown mutual coupling is investigated and an array switching strategy is designed based on the structural features of CLA. Specifically, the two subarrays are separately activated as received arrays first, where the mutual coupling effect can be mitigated to a great extent since the coupling between subarrays is avoided. Then we exploit the uncalibrated polynomial rooting to the two subarrays respectively and obtain the initial DOA estimates by exploiting the coprime property, which is less susceptible to mutual coupling benefiting from sparse arrangement of CLA. Subsequently, the total CLA is activated to collect the signals and the contaminated steering vector is remodeled for decoupling, which is carefully extended from the reconstruction in ULA utilizing two selection matrices. Furthermore, the closed-form solutions to the mutual coupling coefficients (MCCs) are provided by employing a quadratic optimization problem. Finally, re-estimation with few iterations is performed to refine the parameter estimates. Numerical simulations corroborate that the proposed array switching-based scheme for CLA can effectively deal with the mutual coupling and achieve improved estimation accuracy.
We summarize the main contributions of our work as follows.
(1)
The unknown mutual coupling in CLA potentially degrades the estimation performance, whereas the conventional calibration methods for uniform arrays are difficult to apply to CLA due to its nonuniform structure. To tackle this issue, we comprehensively investigate the characteristics of mutual coupling in CLA and significantly mitigate the mutual coupling by exploiting the inherent sparse structural features of CLA.
(2)
We propose an array switching strategy, which can be developed for online calibration, by separately activating the subarrays of CLA to considerably alleviate the severe mutual coupling caused by the two interleaved subarrays in CLA and calculate the well-performed DOA estimates with the signals from subarrays based on the coprime property.
(3)
We reconstruct the contaminated received signal of the total CLA to directly solve the mutual coupling coefficients by utilizing the initial DOA estimates and, in turn, calculate the refined DOA estimates via an iteration procedure. In particular, the reconstruction of the steering vector of CLA for decoupling can be extended to nonuniform linear arrays of arbitrary structure.
The remainder of the paper is given as follows. We provide the preliminary for our work in Section 2, including the data model of CLA and the analysis of mutual coupling characteristics in CLA. Section 3 elaborates the proposed scheme via the array switching strategy and the performance analysis is provided in Section 4. Section 5 exhibits numerical simulations and Section 6 concludes this paper.
Notations: Matrices and vectors are denoted by upper-case and lower-case bold characters, respectively. ( ) T and ( ) H are transpose and conjugate transpose operators, respectively. x a , b represents the set of integers which satisfies a x b . I M denotes an M × M identity matrix. 0 M × N stands for an M × N matrix in which all elements are zero. d i a g ( C ) symbolizes a matrix composed of the principal diagonal elements of C . Toeplitz { c } generates a symmetric Toeplitz matrix having c as its first column and first row. · F is Frobenius norm. angle ( ) represents the phase operator. t r ( · ) is the trace of matrix.

2. Preliminaries

In this section, we first introduce the data model of CLA in the absence of mutual coupling. Subsequently, we modify the data model with mutual coupling and explore the characteristics of mutual coupling in CLA.

2.1. Data Model without Mutual Coupling

The CLA under study in this paper is composed of two uniform linear subarrays with M 1 and M 2 omnidirectional sensors, as illustrated in Figure 1. The inter-element spacing of the subarray 1 with M1 sensors is d 1 = M 2 d while that of the subarray 2 with M 2 sensors is d 2 = M 1 d , where M 1 and M 2 are coprime integers, d = λ / 2 and λ represents the wavelength. Without loss of generality, assume that M 1 < M 2 . The CLA has M = M 1 + M 2 1  sensors in total and the location set of elements can be given by L CLA = S CLA d , where S CLA  is expressed as
S CLA = { m 1 M 2 ,   m 1 0 , M 1 1 } { m 2 M 1 ,   m 2 0 , M 2 1 } .
Assume that K far-field uncorrelated narrowband signals impinging on the CLA from distinct directions θ k ( k 1 , K ) . The received data of the CLA can be represented by
x ( t ) = [ A 1 A 2 ] s ( t ) + [ n 1 ( t ) n 2 ( t ) ] = A s ( t ) + n ( t ) ,
where s ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s K ( t ) ] T represents the signal vector, n i ( t ) is the white Gaussian noise, t 1 , L and L is the number of snapshots. A = [ A 1 T , A 2 T ] T is the total steering matrix of the CLA and the corresponding steering vector is a ( θ k ) = [ a 1 T ( θ k ) , a 2 T ( θ k ) ] T . A i = [ a i ( θ 1 ) , a i ( θ 2 ) , , a i ( θ K ) ] ( i = 1 , 2 ) is the steering matrix of the i-th subarray, where a i ( θ k ) = [ 1 , e j 2 π d i sin θ k / λ , , e j 2 π ( M i 1 ) d i sin θ k / λ ] T . The total steering matrix of the CLA can be also constructed by A c = [ a c ( θ 1 ) , a c ( θ 2 ) , , a c ( θ K ) ] after row exchange, where a c ( θ k ) = [ e j π l 1 sin θ k , e j π l 2 sin θ k , , e j π l M sin θ k ] T and l m S C L A ( m 1 , M ) .
In particular, the CLA can be extracted from a filled ULA (FULA) with inter-element spacing d = λ / 2 and M F = M 1 ( M 2 1 ) + 1 sensors. As a result, we have A c = G 1 A F , where A F = [ a F ( θ 1 ) , a F ( θ 2 ) , , a F ( θ K ) ] denotes the steering matrix of the FULA and a F ( θ k ) = [ 1 , e j π sin θ k , , e j π ( M F 1 ) sin θ k ] T . G 1 is an M × M F selection matrix composed of the t m - th ( t m = l m + 1 ,   m 1 , M ) rows of I M F .
In practice, the covariance of the received signal x ( t ) is usually estimated by using L snapshots
R ^ = 1 L t = 1 L x ( t ) x H ( t ) .  
Performing eigenvalue decomposition of R ^ and we have
R ^ = U ^ s Λ ^ s U ^ s H + U ^ n Λ ^ n U ^ n H ,  
where Λ ^ s and Λ ^ n denotes the diagonal matrices composed of the K largest and the remaining eigenvalues of R ^ , respectively. U ^ s and U ^ n represent the signal subspace and noise subspace composed of the corresponding eigenvectors, respectively.

2.2. Data Model with Mutual Coupling

In the presence of mutual coupling, the model in Equation (2) should be reconstructed by
x ˜ ( t ) = C A s ( t ) + n ( t ) ,  
where C represents the mutual coupling matrix (MCM) and depends on practical factors, e.g., the type of antennas. For instance, the MCM of a linear dipole array with M elements can be written as [24,33,38,39]
C = ( Z A + Z L ) ( Z + Z L I M ) 1 ,  
where Z A and Z L denote the antenna/load impedance, respectively. Z is the mutual impedance matrix and can be expressed as
{ Z = η 0 4 π ( 0.5772 + ln ( 2 β l ) C i ( 2 β l ) + j S i ( 2 β l ) ) , p = q Z = η 0 4 π ( p , q + j p , q ) , p q ,  
where β = 2 π / λ , l denotes the length of dipole antennas and
p , q = sin ( β l ) ( S i ( u 0 ) + S i ( v 0 ) + 2 S i ( u 1 ) 2 S i ( v 1 ) ) + c o s ( β l ) ( C i ( u 0 ) + C i ( v 0 ) 2 C i ( u 1 ) 2 C i ( v 1 ) + 2 C i ( β d p , q ) ) ( 2 C i ( u 1 ) + 2 C i ( v 1 ) 4 C i ( β d p , q ) ) ,
p , q = sin ( β l ) ( C i ( u 0 ) + C i ( v 0 ) + 2 C i ( u 1 ) 2 C i ( v 1 ) ) + c o s ( β l ) ( S i ( u 0 ) S i ( v 0 ) + 2 S i ( u 1 ) + 2 S i ( v 1 ) 2 S i ( β d p , q ) ) + ( 2 S i ( u 1 ) + 2 S i ( v 1 ) 4 S i ( β d p , q ) ) ,
where d p , q is the distance between the p-th element and the q-th element, u 0 = β ( d p , q 2 + l 2 l ) , u 1 = β ( d p , q 2 + 0.25 l 2 0.5 l ) , v 0 = β ( d p , q 2 + l 2 + l ) and v 1 = β ( d p , q 2 + 0.25 l 2 + 0.5 l ) . S i ( u ) and C i ( u ) are defined as
S i ( u ) = u sin x x d x ,   C i ( u ) = u cos x x d x .
In order to facilitate the theoretical study of parameter estimation, a series of simplified models of mutual coupling are established based on certain assumptions [24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40]. In this paper, we assume that c p , q can be characterized by [24,33,38,39]
c p , q = { 0 , | d p d q | B c | d p d q | , | d p d q | < B ,
where d p , d q S CLA , c p , q represents the coupling effect between the q-th sensor and the p-th sensor in CLA. c b = c | d p d q | are defined as the MCCs which satisfy 1 = c 0 > | c 1 | > | c 2 | > > | c B | = 0 . The threshold of mutual coupling range B = 3 and B = 5 are commonly selected assumptions in the existing work [32,40,41]. Specifically, as the mutual coupling depends on practical factors, the value of B is generally determined by the actual scenario. To explain it more intuitively and provide a reliable foundation for the proposed scheme, we provide the S-parameter characteristics of two kinds of arrays below as an example. In practical applications, S-parameters are typically used to characterize the mutual coupling and S p , q captures the coupling degree between the p-th element and the q-th element [42]. Figure 2 exhibits the S-parameter characteristics from an electromagnetic simulation software, i.e., HFSS, where two ten-element ULAs composed of microstrip antennas and dipole antennas are studied. The operating frequencies corresponding to the two arrays are 0.8 GHz and 2.95 GHz. It is generally believed that the mutual coupling is relatively weak when S p , q is lower than -40dB [43,44], and by assuming that the mutual coupling effect can be ignored when S p , q < 55 dB , the value of B in the above two arrays can take 5 and 6, respectively.
According to the labels of the sensors in Figure 3, C can be partitioned by
C = [ 1 c 1 , 2 c 1 , 3 c 1 , M 1 c 1 , M 1 + 1 c 1 , M 1 + 2 c 1 , M c 2 , 1 1 c 2 , 3 c 2 , M 1 c 2 , M 1 + 1 c 2 , M 1 + 2 c 2 , M c 3 , 1 c 3 , 2 1 c 3 , M 1 c M 1 , 1 c M 1 , 2 c M 1 , 3 1 c M 1 , M 1 + 1 c M 1 , M 1 + 2 c M 1 , M c M 1 + 1 , 1 c M 1 + 1 , 2 c M 1 + 1 , M 1 1 c M 1 + 1 , M 1 + 2 c M 1 + 1 , M c M 1 + 2 , 1 c M 1 + 2 , 2 c M 1 + 2 , M 1 c M 1 + 2 , M 1 + 1 1 c M 1 + 2 , M c M , 1 c M , 2 c M , M 1 c M , M 1 + 1 c M , M 1 + 2 1 ] , = [ D 1 B 1 B 2 D 2 ]
where D i M i × M i denotes the MCM in the i-th ( i = 1 , 2 ) subarray while B 1 M 1 × M 2 and B 2 M 2 × M 1 capture the coupling effect between the two subarrays. In particular, C = C T .
Since the inter-element spacings in the two subarrays of CLA are large, we can conclude that the mutual coupling in CLA is mainly caused by the interaction between the two subarrays. It is noteworthy that most of the elements in D 1 and D 2 are zero, and especially when M i B , no mutual coupling exists in the j-th subarray ( i = 1 , 2 ; j = 3 i ) , which motivates us to propose the array switching scheme to mitigate the mutual coupling in CLA by exploiting the inter-element spacings of the two subarrays.

3. The Proposed Parameter Estimation Scheme

In this section, we introduce the proposed scheme to deal with mutual coupling for improved parameter estimation performance through the array switching strategy.
Figure 4 presents the implementation flow of the proposed scheme. Based on the inherent sparse structural characteristic of CLA, we first estimate the initial DOAs by separately activating the two subarrays, which can avoid the mutual coupling resulting from the interlaced subarrays. Then the two interactive subarrays receive signals at the same time and the MCCs of the total CLA can be estimated by utilizing the initial DOA estimates. Furthermore, re-estimation is performed by an iteration procedure based on the estimated MCM to refine parameter estimates.

3.1. Initial DOA Estimation

To begin with, the two subarrays are separately activated to collect signals X 1 i n i and X 2 i n i , which are inherently less susceptible to mutual coupling benefiting from the large inter-element spacings. As a result, we directly apply the uncalibrated polynomial rooting technique to the two subarrays respectively, then the initial unambiguous DOAs can be determined by utilizing the coprime property [17].
In this situation, the MCMs of the two subarrays are assumed to be identity matrices and then the MUSIC spectrum function of the i-th subarray can be represented by [17]
f i ( θ ) = 1 a i H ( θ ) U i , n i n i ( U i , n i n i ) H a i H ( θ ) ,  
where U i , n i n i is the noise subspace corresponding to X i i n i . Let z i = e j 2 π d i sin θ / λ , we have a i ( θ ) = a ( z i ) = [ 1 , z i , , z i M i 1 ] T and the DOA estimation is transformed into solving the roots of the polynomial [11]
a i T ( z 1 ) U i , n i n i ( U i , n i n i ) H a i ( z ) = 0 ,  
Consequently, the initial ambiguous DOA estimates can be achieved by
sin ( θ ^ i , k i n i ) = a n g l e ( z ^ i , k ) λ / 2 π d i   ( k 1 , K ) ,  
The phase ambiguity, which stems from the inter-element spacing larger than λ / 2 , can be eliminated based on the coprime property [17]. As a result, the true estimates θ ^ k i n i can be uniquely determined from the K pairs of closest estimates of the two subarrays. Once the DOA estimates θ ^ k i n i free from severe mutual coupling is obtained, we can further use it to estimate the MCCs by receiving the signal X of the total CLA.

3.2. Mutual Coupling Estimation

Herein, we first reconstruct the contaminated steering vector of the whole CLA to decouple the DOA and mutual coupling. Subsequently, a quadratic optimization problem is established to achieve the mutual coupling estimation. Finally, refined DOA estimates can be obtained via an iteration procedure.
Due to the banded symmetric Toeplitz characteristic of the MCM in ULA [24], the contaminated steering vector of the FULA can be transformed into C F a F ( θ ) = T ( a F ( θ ) ) c F , where C F M F × M F and c F = [ 1 , c 1 , c 2 , , c B 2 ] T are the MCM and mutual coupling vector of the FULA, respectively. The steering vector transformation matrix T ( a F ( θ ) ) M F × B can be obtained by
T ( a F ( θ ) ) = T 1 ( a F ( θ ) ) + T 2 ( a F ( θ ) ) ,  
[ T 1 ( a F ( θ ) ) ] p , q = { [ a F ( θ ) ] p + q 1 , p + q M F + 1 0 , otherwise ,
[ T 2 ( a F ( θ ) ) ] p , q = { [ a F ( θ ) ] p q 1 , p q 2 0 , otherwise .
Due to the nonuniform array structure, the MCM of CLA no longer has such special banded symmetric Toeplitz form. As a result, the reconstruction of the steering vector in ULA cannot be directly applied for further decoupling. To address this, we utilize two matrices G 0 and G 1 derived from the relation between FULA and CLA to reconstruct the contaminated steering vector of CLA as
a ˜ ( θ ) = C a ( θ ) = G 1 C F ( G 0 a F ( θ ) ) = G 1 T ( a F ( θ ) ) c F = H ( θ ) c F ,
where a F ( θ ) = G 0 a F ( θ ) and G 0 M F × 1 is a column vector whose t m - th ( t m = d m + 1 ,   m 1 , M ) elements are 1 and the rest entries are 0. It is noteworthy that this reconstruction method of the contaminated steering vector is applicable to other nonuniform linear arrays.
Based on the transformation in Equation (19), the mutual coupling estimation can be achieved from the following quadratic optimization problem
min c F c F H H H ( θ ^ i n i ) U n U n H H ( θ ^ i n i ) c F , s . t . e H c F = 1 ,  
where U n is the noise subspace of the signal X received by the whole CLA and e 1 = T B × 1 . Construct the Lagrange function as
L ( θ , c F ) = c F H Q ( θ ) c F ε ( e 1 H c F 1 ) ,  
where Q ( θ ) = H H ( θ ) U n U n H H ( θ ) , ε stands for a Lagrange multiplier. By taking the partial derivative of L ( θ , c F ) , i.e., L ( θ , c F ) / c F = 2 Q ( θ ) c F ε e 1 = 0 , we can obtain the estimates of c F by
c ^ F , k = Q 1 ( θ ^ k i n i ) e 1 e 1 H Q 1 ( θ ^ k i n i ) e 1 ,  
c ^ F = 1 K k = 1 K c ^ F , k .  
Subsequently, the MCM of CLA can be constructed by
C ^ = G 1 C ^ F G 1 T ,  
where C ^ F = Toeplitz { c ^ F , 0 1 × ( M F B ) } .

3.3. Iteration Procedure for Refined Estimation

Since there exists a certain deviation between the spectral function in Equation (13) and the actual one contaminated by mutual coupling, the performance of the initial estimation is still potentially degraded. For further improvement in estimation performance, an iteration approach based on the initial MCM estimate C ^ is designed in this part. The detailed steps are outlined in Table 1.
Note that the conventional iteration-based mutual coupling calibration algorithm [30,36] usually sets C ( 0 ) = I , whereas the convergence speed and estimation accuracy heavily depends on C ( 0 ) . The proposed scheme can provide well-performed initial DOA estimates which contributes to the estimation of MCM, thus it requires few iterations and can gain superior estimation performance simultaneously.

3.4. Procedure of the Proposed Scheme

The detailed steps for implementation of the proposed scheme is summarized below.
1. Activate the two subarrays of CLA separately and compute R i i n i with X i i n i to obtain U i , n i n i ( i = 1 , 2 ) .
2. Employ the polynomial rooting according to Equations (14)–(15) and distinguish the unambiguous angles θ ^ k i n i based on the coprime property.
3. Activate the whole CLA and calculate R with the received signal X to obtain U n .
4. Construct Q ( θ ) and achieve the mutual coupling estimation according to Equations (20)–(24).
5. Perform the re-estimation to further improve the estimation accuracy.
Remark 1.
It is assumed that the DOAs and statistical characteristics of the signals are constant during the array switching process.
Remark 2.
Based on the steering vector reconstruction in Equation (19), the well-known rank-reduction (RARE) estimator [33,34] can be extended from ULA to CLA to achieve DOA and mutual coupling estimation, which is provided in the comparison and termed as RARE-based calibration in this paper.

4. Performance Analysis

4.1. Complexity Analysis

The computational complexity of the proposed scheme is analyzed in this part and it is represented by the number of complex multiplications. In the process of initial DOA estimation, calculating the covariance R i i n i ( i = 1 , 2 ) with L snapshots requires O ( L ( M 1 2 + M 2 2 ) ) and performing eigenvalue decomposition costs O ( M 1 3 + M 2 3 ) . Constructing the polynomial needs O ( 2 M 1 ( M 1 K ) + M 1 + 2 M 2 ( M 2 K ) + M 2 ) and the complexity of root-finding operation is O ( ( M 1 1 ) 3 + ( M 2 1 ) 3 ) . As for the mutual coupling estimation, constructing Q ( θ i n i ) takes O ( M 3 + L M 2 + B 2 M + M M F B + 2 M ( M K ) B ) operations and calculating Q 1 ( θ i n i ) costs O ( B 3 ) . The complexity of the recovery of MCM is O ( M M F 2 + M 2 M F ) . The complexity calculation of the iteration process is similar to that of the previous part. Table 2 lists the complexity of the proposed scheme in steps, where n i t e r stands for the iteration times.
In addition, we compare the complexity of the proposed scheme with the RARE-based calibration scheme. Table 3 gives the total complexities of the two schemes, where n 1 = 180 / d s is the peak search times and d s is the searching grid. The complexity comparisons of the two schemes are illustrated in Figure 5 and Figure 6, where B = 6 , d s = 0.001 and n i t e r = 3 . Figure 5 exhibits the complexity comparison versus M 2 , where M 1 = 5 , K = 3 and L = 500 , while the comparison versus snapshots is presented in Figure 6, where M 1 = 5 , M 2 = 6 and K = 3 . It can be seen clearly in Figure 5 and Figure 6 that the proposed scheme requires significantly lower computational cost than the RARE-based calibration method, which is attractive for online calibration.

4.2. Mutual Coupling Analysis

In this part, we employ the coupling leakage [38] to measure the inhibition of the proposed scheme on mutual coupling which is defined as γ = C d i a g ( C ) F / C F . Conceptually, the smaller γ , the weaker the mutual coupling is.
Table 4 exhibits the coupling leakage comparison of the ULA, the general CLA [12] and the CLA based on array switching with the number of sensors, where the MCCs are provided as in Equation (11) with c 1 = 0.3 e j π / 3 and c b = c 1 e - j π ( b 1 ) / 8 / b for 2 b < B . As depicted in Table 4, compared with the general CLA, the mutual coupling effect is significantly alleviated in the CLA by exploiting the proposed array switching strategy, as the severe mutual coupling caused by the interleaved subarrays is avoided.

4.3. Cramer-Rao Bound

To provide a benchmark of parameter estimation performance, the Cramer-Rao Bound (CRB) [24] with CLA under unknown mutual coupling is derived. Define the vector of unknown parameters as
p T = [ θ 1 , θ 2 , , θ K , ρ 1 , ρ 2 , , ρ B 1 , η 1 , η 2 , , η B 1 ] ,
where θ k ( k 1 , K ) represents the DOAs of the signals. ρ b and η b ( b 1 , B 1 ) are the real and imaginary parts of the MCCs. The CRB corresponding to p can be determined by
CRB = d i a g ( F - 1 ) ,  
where F is the Fisher information matrix (FIM) and can be represented by
F = [ F θ θ F θ ρ F θ η F ρ θ F ρ ρ F ρ η F η θ F η ρ F η η ] .
The (u,v)-the element of the FIM can be specified as
F u , v = L × tr [ R 1 R p u R 1 R p v ] = 2 L × Re { tr [ D u R s A ˜ H R 1 A ˜ R s D v R 1 ] + tr [ D u R s A ˜ H R 1 D v H R s A ˜ H R 1 ] }
where A ˜ = C A , D u = A ˜ / p u and D v = A ˜ / p v .

4.4. Advantages

The proposed scheme for CLA against mutual coupling has the following advantages:
(1)
The proposed scheme can be employed as an online calibration technique, which requires no extra auxiliary sources or auxiliary sensors.
(2)
The proposed scheme can significantly alleviate the mutual coupling by exploiting the structural characteristics of CLA. In particular, it outperforms the RARE-based and iterative calibration methods in parameter estimation, which is illustrated in Section 5.
(3)
The proposed scheme is computationally efficient since no spectral search is required, which is attractive in practical applications.

5. Simulation Results

In this section, extensive simulations are carried out to corroborate the superiority of the proposed scheme. The root mean square error (RMSE) with 1000 Monte Carlo trials is defined as
RMSE = 1 1000 K ν = 1 1000 k = 1 K ( θ k θ ^ k , ν ) 2 ,  
where θ k and θ ^ k , ν denote the theoretical DOA and the estimate for the v-th trial of the k-th signal, respectively.

5.1. Verification of the Parameter Estimation

We provide the DOA estimation results of the proposed scheme in Figure 7, where M 1 = 5 , M 2 = 6 , SNR = 10 dB , L = 500 and B = 6 . Consider K = 3 narrowband uncorrelated signals with [ θ 1 , θ 2 , θ 3 ] = [ 30 , 40 , 50 ] . The MCCs are set as c 1 = 0.7 e j π / 3 and c b = c 1 e - j π ( b 1 ) / 4 / b , 2 b < B . It is illustrated in Figure 7 that all signals can be resolved correctly. In addition, Figure 8 presents the estimation results of MCCs in polar coordinates and Table 5 lists the mean values and estimation biases of the estimated MCCs, where the estimation bias is defined as ε r = ν = 1 1000 c b c ^ b , ν F / c b F / 1000 , c b and c ^ b , ν represent the theoretical value and the estimate for the v-th trial of the b-th MCC, respectively. The results demonstrate that the proposed scheme can obtain well-performed MCC estimates, which, in turn, contribute to the refined DOA estimates.

5.2. RMSE Performance of Different Schemes

In this example, we compare the RMSE results of different schemes, including the proposed scheme without iterations (marked as ‘‘proposed scheme-initial’’), the proposed scheme with iterations (marked as ‘‘proposed scheme-refined’’), the mutual coupling calibration schemes (the RARE-based, auxiliary source-based [35] and iterative calibration [36] methods), the uncalibrated root-MUSIC algorithm with unknown mutual coupling, where K = 2 , [ θ 1 , θ 2 ] = [ 30 , 40 ] , M 1 = 6 , M 2 = 7 , B = 8 , c 1 = 0.7 e j π / 3 , c b = c 1 e - j π ( b 1 ) / 8 / b and the searching grid of spectral search is set to d s = 0.0001 . The estimation performance of the root-MUSIC [45] in the absence of mutual coupling and the CRB are also plotted for reference.
Figure 9 exhibits the RMSE results versus signal-to-noise ratio (SNR), where L = 500 . It is indicated explicitly that the uncalibrated root-MUSIC method [45] performs the worst and achieves little improvement in performance with the increase of SNR, which signifies that even though the mutual coupling effect in CLA is relatively weak, it still degrades the performance. By comparison, since the subarrays are less susceptible to mutual coupling benefiting from the large inter-element spacing, the proposed scheme can obtain well-performed initial DOA estimates, which is even superior to the RARE-based algorithm at low SNRs and contributes to the estimation of MCM. As the mutual coupling is not compensated or calibrated in the initial estimation process, the estimation performance of the initial estimates is inevitably restricted. By comparison, it is observed that refined DOA estimates, especially at high SNRs, can be achieved via an iteration procedure by utilizing the MCM estimates. In particular, the proposed scheme can obtain similar DOA estimates to the method with auxiliary source. The estimation performance of the auxiliary source-based method generally improves with the increase of the number of auxiliary sources to some extent, whereas it is exceedingly difficult to implement in practical systems.
Meanwhile, the RMSE comparison in terms of snapshots is captured in Figure 10, where SNR = 10 dB . It is illustrated that the RMSE results improve with the number of snapshots increasing due to the more accurate estimated covariance matrix. We can also conclude from Figure 10 that the proposed scheme can provide as good DOA estimation performance comparable to the calibration methods, indicating that the array switching strategy can effectively alleviate the mutual coupling effect.

5.3. RMSE Performance of Different Mutual Coupling

In this example, we provide the RMSE results of the proposed scheme with different B versus SNR and snapshots in Figure 11 and Figure 12, where K = 2 , [ θ 1 , θ 2 ] = [ 30 , 40 ] , M 1 = 5 , M 2 = 7 and the MCCs are the same as those in Section 5.2. In the two examples, we set L = 500 and SNR = 10 dB , respectively. The RMSE results of the RARE-based calibration are also given for comparison. It is indicated by Figure 11 and Figure 12 that the DOA estimation performance of the two methods deteriorates with strong mutual coupling effect. Besides, the proposed scheme yields more accurate DOA estimates than the RARE-based calibration in all three cases and has better robustness of the range of mutual coupling.

6. Conclusions

In this paper, we have presented an array switching strategy for DOA estimation with CLA in the presence of mutual coupling. Due to the nonuniform array configuration of CLA, the conventional calibration methods for ULAs become difficult to apply. Motivated by the sparse arrangement and the characteristics of mutual coupling in CLA, we design an array switching-based scheme to achieve the DOA estimation against mutual coupling. In the proposed scheme, the two subarrays are separately activated first to collect the signals and estimate the initial DOAs, then the total CLA is used to achieve the mutual coupling estimation and further re-estimate for refined estimates. As a consequence, the well-performed DOA estimates free from the severe mutual coupling effect can be obtained. Numerical simulations corroborate the superiority of the proposed scheme for CLA with regard to computational complexity, mutual coupling and parameter estimation performance. Considering that the mutual coupling is more complicated in practice, the calibration technique with modified coupling model, e.g., direction-dependent model, based on engineering applications will constitute our future work.

Author Contributions

J.S. designed the original idea, performed the simulations and wrote the manuscript; Y.H. contributed to the background and performance analysis; J.L. provided guidance of the study. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by China NSF Grants (61371169, 61631020, 61971218, 61601167), the Fundamental Research Funds for the Central Universities (NT2019013).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, X.; Xu, L.; Xu, L.; Xu, D. Direction of departure (DOD) and direction of arrival (DOA) estimation in MIMO radar with reduced-dimension MUSIC. IEEE Commun. Lett. 2010, 14, 1161–1163. [Google Scholar] [CrossRef]
  2. Li, J.; Zhang, X.; Chen, H. Improved two-dimensional DOA estimation algorithm for two-parallel uniform linear arrays using propagator method. Signal Process. 2012, 92, 3032–3038. [Google Scholar] [CrossRef]
  3. Wen, F. Computationally efficient DOA estimation algorithm for MIMO radar with imperfect waveforms. IEEE Commun. Lett. 2019, 23, 1037–1040. [Google Scholar] [CrossRef]
  4. Wang, H.; Wan, L.; Dong, M.; Ota, K.; Wang, X. Assistant vehicle localization based on three collaborative base stations via SBL-based robust DOA estimation. IEEE Internet Things J. 2019, 6, 5766–5777. [Google Scholar] [CrossRef]
  5. Zhou, C.; Gu, Y.; He, S.; Shi, Z. A robust and efficient algorithm for coprime array adaptive beamforming. IEEE Trans. Veh. Technol. 2018, 67, 1099–1112. [Google Scholar] [CrossRef]
  6. Wen, F.; Shi, J.; Zhang, Z. Joint 2D-DOD, 2D-DOA and polarization angles estimation for bistatic EMVS-MIMO radar via PARAFAC analysis. IEEE Trans. Veh. Technol. 2019, 69, 1626–1638. [Google Scholar] [CrossRef]
  7. Wang, X.; Wan, L.; Huang, M.; Shen, C.; Zhang, K. Polarization channel estimation for circular and non-circular signals in massive MIMO systems. IEEE J. Sel. Top. Signal Process. 2019, 13, 1001–1016. [Google Scholar] [CrossRef]
  8. Schmidt, R.O. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef] [Green Version]
  9. Roy, R.; Kailath, T. ESPRIT-estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef] [Green Version]
  10. Viberg, M.; Ottersten, B.; Kailath, T. Detection and estimation in sensor arrays using weighted subspace fitting. IEEE Trans. Signal Process. 1991, 39, 2436–2449. [Google Scholar] [CrossRef]
  11. Barabell, A. Improving the resolution performance of eigenstructure-based direction-finding algorithms. In Proceedings of the 1983 IEEE International Conference on Acoustics, Speech, and Signal Processing, Boston, MA, USA, 14–16 April 1983; pp. 336–339. [Google Scholar]
  12. Vaidyanathan, P.P.; Pal, P. Theory of sparse coprime sensing in multiple dimensions. IEEE Trans. Signal Process. 2011, 59, 3592–3608. [Google Scholar] [CrossRef]
  13. Pal, P.; Vaidyanathan, P.P. A novel array structure for directions-of-arrival estimation with increased degrees of freedom. In Proceedings of the 2010 IEEE International Conference on Acoustics, Speech and Signal Processing, Dallas, TX, USA, 14–19 March 2010; pp. 2606–2609. [Google Scholar]
  14. Zhou, C.; Gu, Y.; Fan, X.; Shi, Z.; Mao, G.; Zhang, Y.D. Direction-of-arrival estimation for coprime array via virtual array interpolation. IEEE Trans. Signal Process. 2018, 66, 5956–5971. [Google Scholar] [CrossRef]
  15. Tan, Z.; Eldar, Y.C.; Nehorai, A. Direction of arrival estimation using co-prime arrays: A super reso-lution viewpoint. IEEE Trans. Signal Process. 2014, 62, 5565–5576. [Google Scholar] [CrossRef] [Green Version]
  16. Li, J.; Jiang, D.F.; Zhang, X. DOA estimation based on combined unitary ESPRIT for coprime MIMO radar. IEEE Commun. Lett. 2016, 21, 96–99. [Google Scholar] [CrossRef]
  17. Zhou, C.; Shi, Z.; Gu, Y.; Shen, X. DECOM: DOA Estimation with Combined MUSIC for Coprime Array. In Proceedings of the 2013 International Conference on Wireless Communications and Signal Processing (WCSP), Hangzhou, China, 24–26 October 2013; pp. 1–5. [Google Scholar]
  18. Sun, F.; Gao, B.; Chen, L.; Lan, P. A low-complexity ESPRIT-based DOA estimation method for co-prime linear arrays. Sensors 2016, 16, 1367. [Google Scholar] [CrossRef] [Green Version]
  19. Sun, F.; Gao, B.; Lan, P. Partial spectral search-based DOA estimation method for co-prime linear arrays. Electron. Lett. 2015, 51, 2053–2055. [Google Scholar] [CrossRef]
  20. Pal, P.; Vaidyanathan, P.P. Coprime sampling and the MUSIC algorithm. In Proceedings of the 2011 Digital Signal Processing and Signal Processing Education Meeting (DSP/SPE), Sedona, AZ, USA, 4–7 January 2011; pp. 289–294. [Google Scholar]
  21. Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A.; Qu, F. Source estimation using coprime array: A sparse reconstruction perspective. IEEE Sens. J. 2016, 17, 755–765. [Google Scholar] [CrossRef]
  22. Zhou, C.; Zhou, J. Direction-of-Arrival Estimation with Coarray ESPRIT for Coprime Array. Sensors 2017, 17, 1779. [Google Scholar] [CrossRef] [Green Version]
  23. Cheng, Z.; Zhao, Y.; Li, H.; Shui, P. Two-dimensional DOA estimation algorithm with co-prime array via sparse representation. Electron. Lett. 2015, 51, 2084–2086. [Google Scholar] [CrossRef]
  24. Friedlander, B.; Weiss, A.J. Direction finding in the presence of mutual coupling. IEEE Trans. Antennas Propag. 1991, 39, 277–284. [Google Scholar] [CrossRef]
  25. Li, J.; Zhang, X. Direction of Arrival Estimation of Quasi-Stationary Signals Using Unfolded Coprime Array. IEEE Access 2017, 5, 6538–6545. [Google Scholar] [CrossRef]
  26. Qin, S.; Zhang, Y.; Amin, M. Generalized coprime array configurations for direction-of-arrival estimation. IEEE Trans. Signal Process. 2015, 63, 1377–1390. [Google Scholar] [CrossRef]
  27. Zheng, W.; Zhang, X.; Wang, Y.; Zhou, M.; Wu, Q. Extended Coprime Array Configuration Generating Large-Scale Antenna Co-Array in Massive MIMO System. IEEE Trans. Veh. Technol. 2019, 68, 7841–7853. [Google Scholar] [CrossRef]
  28. Raza, A.; Liu, W.; Shen, Q. Thinned Coprime Array for Second-Order Difference Co-Array Generation with Reduced Mutual Coupling. IEEE Trans. Signal Process. 2019, 67, 2052–2065. [Google Scholar] [CrossRef] [Green Version]
  29. Ye, Z.; Dai, J.; Xu, X. DOA estimation for uniform linear array with mutual coupling. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 280–288. [Google Scholar]
  30. Zhang, Y.; Liu, C.; Xu, X.; Ye, Z. Direction of arrival estimation in the presence of array perturbations. In Proceedings of the 2010 IEEE International Conference on Education Technology & Computer, Shanghai, China, 22–24 June 2010; pp. 12–16. [Google Scholar]
  31. Wang, X.; Meng, D.; Huang, M.; Wan, L. Reweighted regularized sparse recovery for DOA estimation with unknown mutual coupling. IEEE Commun. Lett. 2019, 23, 290–293. [Google Scholar] [CrossRef]
  32. Liao, B.; Zhang, Z.; Chan, S. DOA estimation and tracking of ULAs with mutual coupling. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 891–905. [Google Scholar] [CrossRef] [Green Version]
  33. Lin, M.; Yang, L. Blind calibration and DOA estimation with uniform circular arrays in the presence of mutual coupling. IEEE Antennas Wirel. Propag. Lett. 2006, 5, 315–318. [Google Scholar] [CrossRef]
  34. Qi, C.; Wang, Y.; Zhang, Y.; Chen, H. DOA estimation and self-calibration algorithm for uniform circular array. Electron. Lett. 2005, 41, 1092–1094. [Google Scholar] [CrossRef]
  35. See, C.M.S. Sensor array calibration in the presence of mutual coupling and unknown sensor gains and phases. Electron. Lett. 1994, 30, 373–374. [Google Scholar] [CrossRef]
  36. Dai, J.; Zhao, D.; Ye, Z. DOA estimation and self-calibration algorithm for nonuniform linear array. In Proceedings of the 2010 IEEE International Symposium on Intelligent Signal Processing and Communication Systems, Chengdu, China, 6–8 December 2010; pp. 1–4. [Google Scholar]
  37. BouDaher, E.; Ahmad, F.; Amin, M.G.; Ahmad, H. Mutual coupling effect and compensation in non-uniform arrays for direction-of-arrival estimation. Digit. Signal Process. 2017, 61, 3–14. [Google Scholar] [CrossRef] [Green Version]
  38. Liu, C.L.; Vaidyanathan, P.P. Super Nested Arrays: Linear Sparse Arrays with Reduced Mutual Coupling—Part I: Fundamentals. IEEE Trans. Signal Process. 2016, 64, 3997–4012. [Google Scholar] [CrossRef]
  39. Svantesson, T. Modeling and estimation of mutual coupling in a uniform linear array of dipoles. In Proceedings of the 1999 IEEE International Conference on Acoustics, Speech, and Signal Processing, Phoenix, AZ, USA, 15–19 March 1999; pp. 2961–2964. [Google Scholar]
  40. Shi, J.; Hu, G.; Zhang, X.; Zhou, H. Generalized nested array: Optimization for degrees of freedom and mutual coupling. IEEE Commun. Lett. 2018, 22, 1208–1211. [Google Scholar] [CrossRef]
  41. Liu, C.L.; Vaidyanathan, P.P. Hourglass arrays and other novel 2-D sparse arrays with reduced mutual coupling. IEEE Trans. Signal Process. 2017, 65, 3369–3383. [Google Scholar] [CrossRef]
  42. Henault, S.; Antar, Y. Unifying the Theory of Mutual Coupling Compensation in Antenna Arrays. IEEE Antennas Propag. Mag. 2015, 57, 104–122. [Google Scholar] [CrossRef]
  43. Chen, X.; Zhang, S.; Li, Q. A review of mutual coupling in MIMO systems. IEEE Access 2018, 6, 24706–24719. [Google Scholar] [CrossRef]
  44. Lee, J.Y.; Kim, S.H.; Jang, J.H. Reduction of Mutual Coupling in Planar Multiple Antenna by Using 1D EBG and SRR Structures. IEEE Trans. Antennas Propag. 2015, 63, 4194–4198. [Google Scholar] [CrossRef]
  45. Zhang, D.; Zhang, Y.; Zheng, G.; Zheng, G.; Feng, C.; Tang, J. Improved DOA estimation algorithm for co-prime linear arrays using root-MUSIC algorithm. Electron. Lett. 2017, 53, 1277–1279. [Google Scholar] [CrossRef]
Figure 1. Array configurations. (a) the two subarrays; (b) the coprime linear array (CLA).
Figure 1. Array configurations. (a) the two subarrays; (b) the coprime linear array (CLA).
Sensors 20 01629 g001
Figure 2. S-parameter characteristics of two kinds of antenna array. (a) Microstrip array; (b) Dipole array.
Figure 2. S-parameter characteristics of two kinds of antenna array. (a) Microstrip array; (b) Dipole array.
Sensors 20 01629 g002
Figure 3. The labels of the sensors in CLA.
Figure 3. The labels of the sensors in CLA.
Sensors 20 01629 g003
Figure 4. The implementation flow of the proposed scheme.
Figure 4. The implementation flow of the proposed scheme.
Sensors 20 01629 g004
Figure 5. Complexities of different schemes versus different M 2 .
Figure 5. Complexities of different schemes versus different M 2 .
Sensors 20 01629 g005
Figure 6. Complexities of different schemes versus different snapshots.
Figure 6. Complexities of different schemes versus different snapshots.
Sensors 20 01629 g006
Figure 7. Direction of arrival (DOA) estimation results of the proposed scheme.
Figure 7. Direction of arrival (DOA) estimation results of the proposed scheme.
Sensors 20 01629 g007
Figure 8. Mutual coupling coefficients (MCCs) estimation results of the proposed scheme.
Figure 8. Mutual coupling coefficients (MCCs) estimation results of the proposed scheme.
Sensors 20 01629 g008
Figure 9. Root mean square error (RMSE) of different schemes versus signal-to-noise ratio (SNR) (L=500).
Figure 9. Root mean square error (RMSE) of different schemes versus signal-to-noise ratio (SNR) (L=500).
Sensors 20 01629 g009
Figure 10. RMSE of different schemes versus snapshots (SNR=10dB).
Figure 10. RMSE of different schemes versus snapshots (SNR=10dB).
Sensors 20 01629 g010
Figure 11. RMSE with different mutual coupling versus SNR (L=500).
Figure 11. RMSE with different mutual coupling versus SNR (L=500).
Sensors 20 01629 g011
Figure 12. RMSE with different mutual coupling versus snapshots (SNR=10dB).
Figure 12. RMSE with different mutual coupling versus snapshots (SNR=10dB).
Sensors 20 01629 g012
Table 1. Steps of the iteration procedure.
Table 1. Steps of the iteration procedure.
Procedure: Iterations for refined estimates
Step 1: Initialization: Initialize i t e r = 0 and C ( 0 ) = C ^ .
Step 2: DOA Estimation: Obtain the DOA estimates θ ( i t e r ) by solving the roots of polynomial f ( i t e r ) ( z ) = a T ( z 1 ) ( C ( i t e r ) ) H U n U n H ( C ( i t e r ) a ( z ) ) = 0 , then let i t e r = i t e r + 1 .
Step 3: Mutual Coupling Estimation: Calculate c F ( i t e r ) = Q 1 ( θ ( i t e r ) ) e 1 / ( e 1 H Q 1 ( θ ( i t e r ) ) e 1 ) and further obtain C ( i t e r ) = G 1 C F ( i t e r ) G 1 T .
Step 4: Convergence Determination: Go to step 2 until C ( i t e r ) A ( i t e r 1 ) C ( i t e r 1 ) A ( i t e r 2 ) 2 < δ , where δ is a given threshold, such as 10 - 2 .
Table 2. Complexity of the proposed scheme.
Table 2. Complexity of the proposed scheme.
StepComputational complexity
Initial DOA Estimation O ( L ( M 1 2 + M 2 2 ) + M 1 3 + M 2 3 + 2 M 1 ( M 1 K ) + 2 M 2 ( M 2 K )
+ M 1 + M 2 + ( M 1 1 ) 3 + ( M 2 1 ) 3 )
Mutual Coupling Estimation O ( M 3 + ( L + M F ) M 2 + M ( B 2 + M F 2 + M F B ) + 2 M ( M K ) B + B 3 )
Iteration Process O ( n i t e r ( 2 M ( M K ) + M + ( M F 1 ) 3 ) + ( n i t e r 1 ) ( M F M 2
+ M ( B 2 + M F 2 + M F B ) + 2 M ( M K ) B + B 3 ) )
Total O ( L ( M 1 2 + M 2 2 + M 2 ) + M 1 3 + M 2 3 + M 3 + ( M 1 1 ) 3 + ( M 2 1 ) 3
+ 2 M 1 ( M 1 K ) + 2 M 2 ( M 2 K ) + M 1 + M 2 + n i t e r ( 2 M ( M K ) ( 1 + B )
+ M + ( M F 1 ) 3 + M F M 2 + M ( B 2 + M F 2 + M F B ) + B 3 )
Table 3. Complexities of different schemes.
Table 3. Complexities of different schemes.
SchemeComputational complexity
Proposed O ( L ( M 1 2 + M 2 2 + M 2 ) + M 1 3 + M 2 3 + M 3 + ( M 1 1 ) 3 + ( M 2 1 ) 3
+ 2 M 1 ( M 1 K ) + 2 M 2 ( M 2 K ) + M 1 + M 2 + n i t e r ( 2 M ( M K ) ( 1 + B )
+ M + ( M F 1 ) 3 + M F M 2 + M ( B 2 + M F 2 + M F B ) + B 3 )
RARE-based calibration O ( 2 L M 2 + 2 M 3 + n 1 ( 2 M M F B + 2 M ( M K ) B + M B 2 ) + M M F B
+ M 2 B + 2 M ( M K ) B + B 3 + M M F 2 + M 2 M F )
Table 4. Comparison of coupling leakage γ among uniform linear arrays (ULA) and CLA in two cases.
Table 4. Comparison of coupling leakage γ among uniform linear arrays (ULA) and CLA in two cases.
ConditionULACLA (General)CLA (Switching-Based)
M 1 = 5 , M 2 = 6 , M = 10 , B = 5 0.61770.35300
M 1 = 5 , M 2 = 6 , M = 10 , B = 10 0.62510.37410.1156
M 1 = 8 , M 2 = 13 , M = 20 , B = 5 0.63120.25780
M 1 = 8 , M 2 = 13 , M = 20 , B = 10 0.64240.27360.0654
Table 5. Estimation results of mutual coupling coefficients.
Table 5. Estimation results of mutual coupling coefficients.
Theoretical Value of c b Mean Value of c ^ b Estimation Biases ε r
c 1 0.3500 + 0.6062i0.3484 + 0.6081i0.0139
c 2 0.3381 + 0.0906i0.3391 + 0.0893i0.0241
c 3 0.2021 − 0.1167i0.2000 − 0.1206i0.0555
c 4 0.0453 − 0.1690i0.0529 − 0.1640i0.0805

Share and Cite

MDPI and ACS Style

Shen, J.; He, Y.; Li, J. An Array Switching Strategy for Direction of Arrival Estimation with Coprime Linear Array in the Presence of Mutual Coupling. Sensors 2020, 20, 1629. https://doi.org/10.3390/s20061629

AMA Style

Shen J, He Y, Li J. An Array Switching Strategy for Direction of Arrival Estimation with Coprime Linear Array in the Presence of Mutual Coupling. Sensors. 2020; 20(6):1629. https://doi.org/10.3390/s20061629

Chicago/Turabian Style

Shen, Jinqing, Yi He, and Jianfeng Li. 2020. "An Array Switching Strategy for Direction of Arrival Estimation with Coprime Linear Array in the Presence of Mutual Coupling" Sensors 20, no. 6: 1629. https://doi.org/10.3390/s20061629

APA Style

Shen, J., He, Y., & Li, J. (2020). An Array Switching Strategy for Direction of Arrival Estimation with Coprime Linear Array in the Presence of Mutual Coupling. Sensors, 20(6), 1629. https://doi.org/10.3390/s20061629

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop