Next Article in Journal
On the Prediction of the Threshold Voltage Degradation in CMOS Technology Due to Bias-Temperature Instability
Next Article in Special Issue
Compact Antenna for 4G/5G Metal Frame Mobile Phone Applications Using a Tuning Line
Previous Article in Journal
Design of an Highly Efficient AC-DC-AC Three-Phase Converter Using SiC for UPS Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SBL-Based Direction Finding Method with Imperfect Array

1
State Key Laboratory of Millimeter Waves, Southeast Univerity, Nanjing 210096, China
2
School of Electronic and Information Engineering, Shanghai Dianji University, Shanghai 201306, China
3
School of Information Science and Engineering, Southeast University, Nanjing 210096, China
*
Author to whom correspondence should be addressed.
Electronics 2018, 7(12), 426; https://doi.org/10.3390/electronics7120426
Submission received: 27 November 2018 / Revised: 7 December 2018 / Accepted: 10 December 2018 / Published: 11 December 2018
(This article belongs to the Special Issue Intelligent Antennas)

Abstract

:
The imperfect array degrades the direction finding performance. In this paper, we investigate the direction finding problem in uniform linear array (ULA) system with unknown mutual coupling effect between antennas. By exploiting the target sparsity in the spatial domain, the sparse Bayesian learning (SBL)-based model is proposed and converts the direction finding problem into a sparse reconstruction problem. In the sparse-based model, the off-grid errors are introduced by discretizing the direction area into grids. Therefore, an off-grid SBL model with mutual coupling vector is proposed to overcome both the mutual coupling and the off-grid effect. With the distribution assumptions of unknown parameters including the noise variance, the off-grid vector, the received signals and the mutual coupling vector, a novel direction finding method based on SBL with unknown mutual coupling effect named DFSMC is proposed, where an expectation-maximum (EM)-based step is adopted by deriving the estimation expressions for all the unknown parameters theoretically. Simulation results show that the proposed DFSMC method can outperform state-of-the-art direction finding methods significantly in the array system with unknown mutual coupling effect.

1. Introduction

In the direction finding problem, the traditional discrete Fourier transform (DFT)-based method can only find one signal in one beam-width, so the resolution of such a method is too low to estimate multiple signals. Therefore, the super-resolution methods have been proposed including multiple signal classification (MUSIC) method [1,2], Root-MUSIC method [3], and the estimating signal parameters via rotational invariance techniques (ESPRIT) method [4]. Additionally, the subspace methods have also been improved to estimate the correlation signals, such as the spatial smoothing MUSIC method [5]. However, the subspace methods only distinguish the noise and signal subspaces and have not exploited additional characteristics of the received signals.
The compressed sensing (CS)-based methods have been proposed to estimate the directions by exploiting the signal sparsity in the spatial domain [6,7,8,9,10,11,12,13,14]. Notably, the sparse Bayesian learning (SBL) and the relevance vector machine (RVM) proposed in [15] can achieve better estimation performance in the CS-based direction finding methods, where the directions are estimated by reconstructing the sparse signals in the spatial domain with the corresponding distribution assumptions of unknown parameters. Consequently, the SBL-based CS method, named CS-SBL, is developed in [16] to reconstruct the sparse signals. However, in the CS-SBL method, the discrete grids are adopted to formulate the CS-based system model, so the estimation performance is limited by the grid size. To further improve the estimation performance, dense grids can be adopted, but it will improve the computational complexity in the sparse reconstruction algorithm.
Additionally, the dense grids improve the correlation between the grids and decrease the performance of sparse reconstruction. The off-grid CS-based methods were proposed [17,18] to overcome the grid problem in the CS-based model, such as the off-grid sparse Bayesian inference (OGSBI) method proposed in [19]. In the off-grid CS-based system model, the ground-truth directions are approximated by the Taylor expansion, so the performance of direction estimation can be improved in the off-grid methods when the same grids are adopted. Moreover, by solving the roots of a specific polynomial in an off-grid model, the Root-SBL method [20] was also proposed to decrease the computational complexity of the SBL-based method. The grid evolution method was proposed in [21] to refine the grids for the SBL-based method, and a dictionary learning algorithm is proposed in [22].
In a practical direction finding problem, the imperfection of the antenna array will decrease the estimation performance, so the mutual coupling effect between antennas cannot be ignored [23,24]. The direction finding methods were proposed in [25,26,27] to decrease the mutual coupling effect. However, in the existing sparse-based methods, the unknown mutual coupling effect is not considered, especially in a scenario with off-grid effect.
In this paper, a symmetric Toeplitz matrix [28,29,30] is used to describe the mutual coupling effect, and a novel direction estimation method is proposed. With both the off-grid and the mutual coupling effect, the direction finding problem is investigated. A novel system model is formulated to describe both the off-grid and the mutual coupling effect. Then, by exploiting the signal sparsity in the spatial domain, a novel direction finding method based on SBL with unknown mutual coupling effect, named DFSMC, is proposed. Additionally, with the distribution assumptions, we theoretically derive the estimation of all unknown parameters using the expectation-maximum (EM)-based method in DFSMC, where the unknown parameters include the mutual coupling vector, the noise variance, the signals, the off-grid vector, et al. Finally, the proposed DFSMC method is compared with the state-of-art methods in the direction finding performance. To summarize, we make the contributions as follows:
  • The SBL-based system model with mutual coupling effect: With considering both the off-grid and the unknown mutual coupling problems, a novel system model is formulated and transforms the direction finding problems into a sparse reconstruction problem.
  • The DFSM method for direction finding estimation: With the distribution assumptions of all unknown parameters, a novel SBL-based direction finding method with unknown mutual coupling effect, named DFSMC, is proposed. DFSMC method estimates the directions via updating all the unknown parameters alternatively and achieves better estimation performance than the state-of-art methods.
  • The theoretical estimation expressions for all unknown parameters: In the proposed DFSMC method, the EM method is adopted to estimate all the unknown parameters including the noise variance, the received signals, the mutual coupling vector, and the off-grid vectors, et al. With the distribution assumptions, we theoretically derive the expressions for all the unknown parameters.
The remainder of this paper is organized as follows. The system model for direction finding with unknown mutual coupling effect is formulated in Section 2. The direction finding method based on SBL is given in Section 3. The simulation results are given in Section 4. Finally, Section 5 concludes the paper.
Notations: I N denotes an N × N identity matrix. E · denotes the expectation operation. CN a , B denotes the complex Gaussian distribution with the mean being a and the covariance matrix being B . · 2 , ⊗, Tr · , ( · ) , ( · ) T and ( · ) H denote the 2 norm, the Kronecker product, the trace of a matrix, the conjugate, the matrix transpose and the Hermitian transpose, respectively. R { a } denotes the real part of complex value a. Additionally, for a vector a , a n denotes the n-th entry of a , and diag { a } denotes a diagonal matrix with the diagonal entries from a . For a matrix A , A : , n denotes the n-th column of A , and diag { A } denotes a vector with the entries from the diagonal entries of A .

2. ULA System for Direction Finding

We consider the direction finding problem in the uniform linear array (ULA) system, where N antennas are adopted and the inter-antenna element spacing is d. As shown in Figure 1, K unknown signals with different directions ( θ k , K = 0 , 1 , , K 1 ) are received by the ULA. Thus, the received signals in the N antennas can be expressed as
y ( t ) = k = 0 K 1 C a ( θ k ) s k ( t ) + n ( t ) = CA s ( t ) + n ( t ) ,
where the matrix C C N × N denotes the mutual coupling matrix, and the signals are collected into a vector s ( t ) s 0 ( t ) , s 1 ( t ) , , s K 1 ( t ) T with the k-th signal being s k ( t ) . Then, the received signals in the ULA antennas can be expressed as y ( t ) y 0 ( t ) , y 1 ( t ) , , y N 1 ( t ) T , and the zero-mean additive white Gaussian noise (AWGN) with the variance being σ n 2 is n ( t ) n 0 ( t ) , n 1 ( t ) , , n N 1 ( t ) T . In this paper, we suppose that the noise variance σ n 2 is unknown. A C N × K denotes the steering matrix for the K signals, and can be expressed as
A a ( θ 0 ) , a ( θ 1 ) , , a ( θ K 1 ) ,
where the steering vector for the k-th signal can be written as a ( θ k ) a 0 ( θ k ) , a 1 ( θ k ) , , a N 1 ( θ k ) T , a n ( θ k ) = e j 2 π n d λ sin θ k , and λ denotes the wavelength.
In this paper, we consider the direction finding problem with unknown mutual coupling effect between antennas, and the mutual coupling effect can be described usually by a symmetric Toeplitz matrix [29]. As expressed in (1), the mutual coupling matrix can be represented as
C = 1 c 1 c N 1 c 1 1 c N 2 c N 1 c 1 1 ,
where c n ( n = 1 , 2 , , N 1 ) denotes the mutual coupling coefficient between the n 1 -th antenna and the n 2 -th antenna, and | n 1 n 2 | = n .
The signal model in (1) is a contiiuous domain model, and after the uniform sampling, a discrete model can be obtained in a matrix form as
Y = CAS + N ,
where the sampling interval is T s , the number of the samples is M. Y C N × M , S C K × M and N C N × M are expressed as Y = y ( 0 ) , y ( T s ) , , y ( ( M 1 ) T s ) , S = s ( 0 ) , s ( T s ) , , s ( ( M 1 ) T s ) , N = n ( 0 ) , n ( T s ) , , n ( ( M 1 ) T s ) . To simplify the notations, we define y m y ( m T s ) , n m n ( m T s ) and s m s ( m T s ) , so we have Y = y 0 , y 1 , , y M 1 , S = s 0 , s 1 , , s M 1 , N = n 0 , n 1 , , n M 1 .
However, the system model in (4) is hard to solve directly with the unknown mutual coupling matrix C , so we try to express the matrix C in a vector form. The mutual coupling matrix in (3) can be described alternatively by a vector c as C = Toeplitz { c } , where c 1 , c 1 , , c N 1 T is the first column of C , and Toeplitz { · } denotes the Toeplitz transformation. Therefore, after the simplification, the received signals during the m-th sampling interval in (4) can be rewritten as
y m = CAs m + n m = Q ( s m c ) + n m ,
where the mutual coupling effect is expressed by a vector c , and we use a matrix Q C N × K N to rearrange the steering matrix A .
According to the lemma in [29,31,32], the matrix Q in (5) can be obtained as
Q Q ( θ 0 ) , Q ( θ 1 ) , , Q ( θ K 1 ) .
The k-th sub-matrix Q ( θ k ) C N × N is Q ( θ k ) = Q 1 ( θ k ) + Q 2 ( θ k ) , where Q 1 ( θ k ) and Q 2 ( θ k ) respectively are
Q 1 ( θ k ) a 0 ( θ k ) a 1 ( θ k ) a N 1 ( θ k ) a 1 ( θ k ) a 2 ( θ k ) 0 a N 2 ( θ k ) a N 1 ( θ k ) 0 a N 1 ( θ k ) 0 0 ,
Q 2 ( θ k ) 0 0 0 0 0 a 0 ( θ k ) 0 0 0 a N 3 ( θ k ) a 0 ( θ k ) 0 0 a N 2 ( θ k ) a 1 ( θ k ) a 0 ( θ k ) .
Therefore, by collecting the M samples into a matrix, the received signals in (4) can be finally rewritten as
Y = Q ( S c ) + N .
In this paper, we will propose a high resolution method to estimate the directions ( θ 0 , θ 1 , …, θ K 1 ) from the received signal matrix Y , where the signal matrix S , the mutual coupling vector c and the noise variance σ n 2 are all unknown.

3. Direction Finding Method Based on Sparse Bayesian Learning

In this section, we propose a novel SBL-based method to estimate the directions (named Direction Finding based on SBL with Mutual Coupling effect, DFSMC). The sparse model will be established first, and the DFSMC will be proposed with the distribution assumptions of unknown parameters.

3.1. Sparse-Based Signal Model

Since the received signals are sparse in the spatial domain, we propose a sparse-based model to estimate the directions with unknown mutual coupling effect. In the sparse-based model, the dictionary matrix must be established first, so an over-complete dictionary matrix can be formulated by discretizing the signal direction uniformly in the spatial domain
D Q ( ζ 0 ) , Q ( ζ 1 ) , , Q ( ζ U 1 ) C N × U N ,
where ζ u denotes the u-th discretized direction ( u = 0 , 1 , , U 1 ), U denotes the number of discretized directions, the grid size is defined as δ | ζ u + 1 ζ u | , and we use a vector to contain all the discretized directions ζ ζ 0 , ζ 1 , , ζ U 1 .
With the discretized directions and the system model in (9), a sparse-based on-grid direction finding model can be expressed as
Y = D ( X c ) + n ,
where X is a sparse matrix
X x 0 , x 1 , , x M 1 C U × M .
The structure of sparse matrix X is shown in Figure 2, and the sparse vectors ( x 0 , x 1 , …, x M 1 ) have the same support sets. When the direction of the k-th received signal θ k is equal to the u k -th discretized direction ζ u k , we have X u k , m = S k , m , so the u-th row and m-th column of X is
X u , m = S k , m , u = u k 0 , otherwise .
The sparse-based model in (12) assumes that the directions of received signals are exactly on the discretized grids. However, in the practical direction finding system, when the direction θ k is not on the discretized grids, the direction θ k can be represented by ζ u k , which is a grid nearest to θ k . Thus, the corresponding matrix Q ( θ k ) in (2) can be approximated by
Q ( θ k ) Q ( ζ u k ) + ( θ k ζ u k ) Ω ( ζ u k ) ,
where the first-order derivative is defined as Ω ( ζ u k ) Q ( ζ ) ζ ζ = ζ u k . For example, as shown in Figure 3, the direction of signal s k ( t ) is θ k , and the nearest grid is ζ 3 . Thus, the corresponding matrix Q ( θ k ) in (14) can be written as Q ( θ k ) Q ( ζ 3 ) + ( θ k ζ 3 ) Ω ( ζ 3 ) .
Therefore, with the approximation in (14), the received signal in (9) can be approximated by a sparse-based off-grid model
Y Ψ ( ν ) ( X c ) + N ,
where Ψ ( ν ) D + Ξ diag ν I N , Ξ Ξ 0 , Ξ 1 , , Ξ U 1 , and the u-th submatrix of Ξ is denoted as Ξ u Ω ( ζ u ) . Additionally, a vector ν R U × 1 is used to represent the off-grid directions, and the u-th entry is
ν u = θ k ζ u k , u = u k 0 , otherwise .
Finally, an off-grid sparse-based model is formulated for the direction finding problem in (15). We will estimate the directions by reconstructing the sparse matrix X . The positions of non-zero entries in X indicate the directions of received signals. Simultaneously, the unknown parameters including the mutual coupling vector c , the noise variance σ n 2 and the off-grid vector ν will also be estimated.

3.2. Distribution Assumptions

In the proposed DFSMC method, the sparse Bayesian learning theory is adopted, and the method is established based on the distribution assumptions of all the unknown parameters. We assume that the unknown parameters follow the following distributions:
  • Noise N : Gaussian distribution;
  • The precision of noise variance α n : Gamma distribution;
  • Sparse matrix X : Gaussian distribution;
  • The precision of signal variance ι : Gamma distribution;
  • Mutual coupling vector c : Gaussian distribution;
  • The precision of mutual coupling variance ϑ : Gamma distribution;
  • Off-grid vector ν : Uniform distribution.
The relationships between all the unknown parameters are given in Figure 4, and we will describe the distributions more clear in the following contents.

3.2.1. The Distribution of Noise

When the received signals are independent between different samples, with the assumption of circular symmetric white Gaussian noise, the distribution of noise can be expressed as
p ( N | σ n 2 ) = m = 0 M 1 CN ( n m | 0 N × 1 , σ n 2 I N ) ,
where σ n 2 denotes the noise variance, and the complex Gaussian distribution with the mean being μ and the covariance matrix being Σ is expressed as
CN ( x | μ , Σ ) = 1 π N det ( Σ ) e ( x μ ) H Σ 1 ( x μ ) .

3.2.2. The Distribution of Noise Variance σ n 2

In this paper, the noise variance is unknown. Since the Gamma distribution is a conjugate prior of Gaussian distribution, the posterior distribution also follows a Gamma distribution. Therefore, using the Gamma distribution can simplify the following analysis. With the unknown noise variance σ n 2 , we use a Gamma distribution to describe the precision of noise variance α n σ n 2 , and we have the following Gamma distribution
p ( α n ) = g ( α n ; a , b ) ,
where a and b are the hyperparameters for α n , and g ( α n ; a , b ) Γ 1 ( a ) b a α n a 1 e b α n , Γ ( a ) 0 x a 1 e x d x .

3.2.3. The Distribution of Sparse Matrix

With the independent received signals S among samples, we can assume that the sparse matrix X follows a zero-mean Gaussian distribution
p ( X | Λ x ) = m = 0 M 1 CN ( x m | 0 U × 1 , Λ x ) ,
where the covariance matrix Λ x R U × U is a diagonal matrix with the u-th diagonal entry being σ x , u 2 .

3.2.4. The Distribution of Signal Variance

Similarity, with the unknown signal variance Λ x , we define the precision vector ι ι 0 , ι 1 , , ι U 1 T , where ι u σ x , u 2 , so ι can be expressed by a Gamma prior
p ( ι ; c , d ) = u = 0 U 1 g ( ι u ; c , d ) ,
where c and d are the hyperparmaters for ι .

3.2.5. The Distribution of Mutual Coupling Vector

With the unknown mutual coupling vector, when the mutual coupling coefficients are independent between antennas, the distribution of mutual coupling vector c can be expressed as a Gaussian distribution
p ( c | Λ c ) = n = 0 N 1 CN ( c n | 0 , Λ c ) ,
where the covariance matrix Λ c R N × N is a diagonal matrix with the n-th diagonal entry being σ c , n 2 .

3.2.6. The Distribution of Mutual Coupling Variance

By defining the precisions ϑ ϑ 0 , ϑ 1 , , ϑ N 1 T ( ϑ n σ c , n 2 ), we use a Gamma distribution to describe the distribution of ϑ
p ( ϑ ; e , f ) = n = 0 N 1 g ( ϑ n ; e , f ) ,
where both e and f are the hyperparameters of ϑ .

3.2.7. The Distribution of Off-Grid Vector

We can assume that the off-grid vector ν follows a uniform distribution, and the distribution of ν u can be expressed as
p ( ν u ; δ ) = U ν u 1 2 δ , 1 2 δ ,
where the uniform distribution is defined as
U x a , b 1 b a , a x b 0 , otherwise .

3.3. DFSMC Method

With the distribution assumptions of unknown parameters, a novel direction finding method based on the SBL is proposed with the unknown mutual coupling effect, named DFSMC. In the SBL-based method, the posterior probabilities for all the unknown parameters are theoretically derived.
To estimate the directions, we formulate the following problem to maximize the posterior probability
ˆ = arg max p ( | Y ) ,
where a set X , ν , c , σ n 2 , ι , ϑ is used to contain all the unknown parameters. However, the problem (26) is too complex and cannot be solved directly. The expectation maximum (EM)-based method is used to realize the proposed DFSMC method. Additionally, with the received signal Y , the joint distribution with unknown parameters can be expressed as
p ( Y , ) = p ( Y | ) p ( X | ι ) p ( c | ϑ ) p ( α n ) p ( ι ) p ( ϑ ) p ( ν ) .
The details to estimate all unknown parameters are given as follows.

3.3.1. The Sparse Matrix

Given the received signal Y and the parameters ( X ) excepting X , the the posterior of X can be expressed as
p ( X | Y , ν , c , α n , ι , ϑ ) = p ( Y | ) p ( X | ι ) p ( Y | ν , c , α n , ι , ϑ ) p ( Y | ) p ( X | ι ) ,
where both p ( Y | ) and p ( X | ι ) follow Gaussian distributions, and can be calculated as
p ( Y | ) = m = 0 M 1 α n N π N e α n y m Ψ ( ν ) ( x m c ) 2 2 ,
p ( X | ι ) = m = 0 M 1 u = 0 U 1 ι u 1 π U e x m H diag { ι } x m .
Therefore, the posterior of X is also a Gaussian function
p ( X | Y , ν , c , α n , ι , ϑ ) = m = 0 M 1 CN ( x m | μ m , Σ X ) ,
where the mean μ m and covariance matrix Σ X are obtained from (29) and (30) as
μ m = α n H Σ X H T H ( ν , c ) y m ,
Σ X = α n T H ( ν , c ) T ( ν , c ) + diag { ι } 1
and we define the following function
T ( ν , c ) Ψ ( ν ) ( I U c ) .
Additionally, to simplify the notations, the u-th entry of μ m is denoted as μ u , m , and we can collect all the mean μ m as a matrix μ μ 0 , μ 1 , , μ M 1 .
To estimate the other unknown parameters X , with (27), we can formulate the following likelihood function
L ( ν , c , α n , ι , ϑ ) = E { ln p ( Y | ) p ( X | ι ) p ( c | ϑ ) p ( α n ) p ( ι ) p ( ϑ ) p ( ν ) } .
where we just use E { · } to represent E X | Y , ν , α n , ι , ϑ { · } . Thus, the EM-based method can be used to estimate X , and the details are given in the following contents. Additionally, the derivatives for the complex vector and matrix are given as the following lemma.
Lemma 1.
With both the complex vectors ( u C P × 1 , v C P × 1 ) and the complex matrix A C M × P being the function of a complex vector x C N × 1 , the following derivations can be obtained
u H v x = v T ( u ) x + u H v x ,
Au x = A x 0 u + A u x 0 , , A x n u + A u x n , .
Proof. 
See: Appendix A. □

3.3.2. The Mutual Coupling Vector

Ignoring terms independent thereof in L ( ν , c , α n , ι , ϑ ) , we can obtain the following likelihood function for the mutual coupling vector c
L ( c ) = E ln p ( Y | X , ν , c , α n ) p ( c | ϑ )   = E ln m = 0 M 1 CN ( y m | Ψ ( ν ) ( x m c ) , α n 1 I N ) + ln n = 0 N 1 CN ( c n | 0 , ϑ n 1 )   α n M G 1 ( c , ν ) m = 0 M 1 α n G 2 , m ( c , ν ) G 3 ( c ) ,
where we have
E ln p ( Y | X , ν , c , α n ) = M N ln α n π α n M G 1 ( c , ν ) m = 0 M 1 α n G 2 , m ( c , ν ) ,
and G 1 ( c , ν ) Tr T H ( ν , c ) T ( ν , c ) Σ X , G 2 , m ( c , ν ) y m Ψ ( ν ) ( μ m c ) 2 2 , G 3 ( c ) n = 0 N 1 ϑ n | c n | 2 .
To estimate the mutual coupling vector c , we can maximize the likelihood function L ( c ) , and we have
c ˆ = arg max c L ( c ) .
Therefore, by setting L ( c ) c = 0 , the mutual coupling vector can be obtained. We can calculate
L ( c ) c = α n M G 1 ( c , ν ) c m = 0 M 1 α n G 2 , m ( c , ν ) c G 3 ( c ) c .
In (41), G 1 ( c , ν ) c , G 2 , m ( c , ν ) c and G 3 ( c ) c can be calculated as follows.
  • For G 1 ( c , ν ) c : With the derivations of complex vector and matrix in Appendix A, G 1 ( c , ν ) c is a row vector, and the n-th entry can be calculated as
    G 1 ( c , ν ) c n = Tr T H ( ν , c ) T ( ν , c ) Σ X c n .
    Additionally, we can calculate
    T H ( ν , c ) T ( ν , c ) Σ X c n = ( I U c ) H c n Ψ H ( ν ) T ( ν , c ) Σ X ] + T H ( ν , c ) Ψ ( ν ) ( I U c ) c n Σ X = T H ( ν , c ) Ψ ( ν ) I U c c n Σ X = T H ( ν , c ) T ( ν , e n N ) Σ X ,
    where e n N is a N × 1 vector with the n-th entry being 1 and other entries being 0. Therefore, the n-th entry in (42) can be simplified as
    G 1 ( c , ν ) c n = c H p = 0 U 1 k = 0 U 1 Ψ p H ( ν ) Ψ k ( ν ) Σ X , k , p e n N ,
    and we finally have the derivation of G 1 ( c , ν ) as
    G 1 ( c , ν ) c = c H p = 0 U 1 k = 0 U 1 Ψ p H ( ν ) Ψ k ( ν ) Σ X , k , p .
  • G 2 , m ( c , ν ) c can be simplified as
    G 2 , m ( c , ν ) c = [ y m Ψ ( ν ) ( μ m c ) ] H Ψ ( ν ) μ m c c   = [ y m Ψ ( ν ) ( μ m c ) ] H Ψ ( ν ) ( μ m I N ) .
  • G 3 ( c ) c can be simplified as G 3 ( c ) c = c H diag { ϑ } .
Therefore, with (41), the mutual coupling vector can be finally estimated as
c ˆ = H 1 z ,
where
H = m = 0 M 1 α n P H ( ν , μ m ) P ( ν , μ m ) + α n M p = 0 U 1 k = 0 U 1 Ψ p H ( ν ) Ψ k ( ν ) Σ X , k , p H + diag { ϑ } ,
z = m = 0 M 1 α n P H ( ν , μ m ) y m ,
and we define P ( ν , μ m ) Ψ ( ν ) ( μ m I N ) .

3.3.3. For the Precision of Signal Variance

Ignoring terms independent thereof in L ( ν , c , α n , ι , ϑ ) , we can obtain the likelihood function of ι as
L ( ι ) = E ln p ( X | ι ) p ( ι ) = E ln m = 0 M 1 CN ( x m | 0 U × 1 , Λ x ) + ln u = 0 U 1 g ( ι u ; c , d ) .
Then, the precision of signal variance can be estimated by ι ˆ = arg max ι L ( ι ) .
By setting L ( ι ) ι = 0 , the u-th entry of ι can be obtained as
ι ˆ u = M + c 1 d + M Σ X , u , u + m = 0 M 1 | μ u , m | 2 .
In the iterative algorithm, (50) can be rewritten as
ι ˆ u i + 1 M 1 ι ˆ u i m = 0 M 1 | μ u , m | 2 d + M Σ X , u , u ,
where ι ˆ u i + 1 and ι ˆ u i are the esitmated results at the ( i + 1 ) -th and the i-th iterations, respectively.

3.3.4. For α n

Ignoring terms independent thereof in L ( ν , c , α n , ι , ϑ ) , we can obtain the likelihood function
L ( α n ) = E ln p ( Y | X , ν , c , α n ) p ( α n )   = E ln m = 0 M 1 CN y m | Ψ ( ν ) ( x m c ) , σ n 2 I + ln g ( α n ; a , b ) .
The precision of noise variance can be estimated by
α ˆ n = arg max α n L ( α n ) .
By setting L ( α n ) α n = 0 , we can obtain
α ˆ n = M N + a 1 M G 1 ( c , ν ) + m = 0 M 1 G 2 , m ( c , ν ) + b .
In the iterative algorithm, (54) can be rewritten as
α ˆ n i + 1 M N 1 α ˆ n i m = 0 M 1 G 2 , m ( c , ν ) M G 1 ( c , ν ) + b . ,
where α ˆ n i + 1 and α ˆ n i are the esitmated results at the ( i + 1 ) -th and the i-th iterations, respectively.

3.3.5. For the Precision of Mutual Coupling Variance

Ignoring terms independent thereof in L ( ν , c , α n , ι , ϑ ) , we can obtain the likelihood function
L ( ϑ ) = E ln p ( c | ϑ ) p ( ϑ ) = E ln n = 0 N 1 CN ( c n | 0 , ϑ n 1 ) + ln n = 0 N 1 g ( ϑ n ; e , f ) .
The precision of mutual coupling variance can be estimated by ϑ ˆ = arg max ϑ L ( ϑ ) .
By setting L ( ϑ ) ϑ = 0 , we can obtain the n-th entry of ϑ as
ϑ ˆ n 1 f + c n H c n .

3.3.6. For the Off-Grid Vector

Ignoring terms independent thereof in L ( ν , c , α n , ι , ϑ ) , we can obtain the likelihood function
L ( ν ) = E ln p ( Y | X , ν , c , α n ) p ( ν ) M G 1 ( c , ν ) m = 0 M 1 G 2 , m ( c , ν ) .
The off-grid vector can be estimated by
ν ˆ = arg max ν L ( ν ) .
Then, G 1 ( c , ν ) ν R 1 × U is a row vector, and the u-th entry is
G 1 ( c , ν ) ν u = Tr T H ( ν , c ) T ( ν , c ) Σ X ν u = Tr 0 , T H ( ν , c ) Ξ u c , 0 Σ X + Tr 0 , T H ( ν , c ) Ξ u c , 0 H Σ X = 2 R m = 0 U 1 c H Ψ m H ( ν ) Ξ u c Σ X , u , m = 2 R T H ( ν , c ) Ξ u c H Σ X , : , u .
G 1 ( c , ν ) ν can be simplified as
G 1 ( c , ν ) ν = 2 R diag Σ X T H ( ν , c ) Ξ ( I U c ) T .
Additionally, G 2 , m ( c , ν ) ν can be obtained as
G 2 , m ( c , ν ) ν = 2 R [ y m Ψ ( ν ) ( μ m c ) ] H Ψ ( ν ) ( μ m c ) ν = 2 R [ y m Ψ ( ν ) ( μ m c ) ] H Ξ ( diag { μ m } c ) .
Therefore, with L ( ν ) ν u = 0 , we can obtain
ν ˆ = G 1 z ,
where the entry of the u-th row in G R U × U is
G u , : = R M c H Ξ u H Ξ ( diag Σ X , : , u c ) + m = 0 M 1 R μ m , u c T Ξ u T Ξ ( I U c ) diag μ m ,
and the u-th entry of z R U × 1 is
z u = m = 0 M 1 R [ y m D ( μ m c ) ] H Ξ u μ m , u c M R c H Ξ u H D ( I U c ) Σ X , : , u .
In Algorithm 1, we show the details of the proposed DFSMC method for the direction finding with the unknown mutual coupling effect. In the proposed DFSMC algorithm, after the iterations, we can obtain the spatial spectrum P X of the sparse matrix X from the received signal Y . Then, by searching all the values of P X , the corresponding peak values can be found. By selecting positions of peak values corresponding to the K maximum values, we can estimate the directions with ζ + ν .
Algorithm 1 DFSMC algorithm for direction finding with the unknown mutual coupling effect.
 1:
Input: received signal Y , the number of samples M, the numbers of iterations N 1 , N 2 and N 3 dictionary matrix D , the first order derivative of dictionary matrix Ξ . Usually, we have N 1 = 10 3 , N 2 = 300 , N 3 = 50 , b = d = f = 10 3 , and a = c = e = 1 + b .
 2:
Initialization: i iter = 1 , s method = 0 , i method = 1 , c ˆ = ϑ ˆ = [ 1 , 0 1 × ( N 1 ) ] T , and ν ˆ = 0 U × 1 .
 3:
Ψ ( ν ˆ ) D + Ξ diag ν ˆ I N .
 4:
while i iter N 1 do
 5:
 Obtain T ( ν , c ) from (34).
 6:
 Obtain the mean μ m ( m = 0 , 1 , , M 1 ) and covariance matrix Σ X from (32) and (33), respectively.
 7:
 Update the precision of noise variance α ˆ n from (55).
 8:
 Update the precision of signal variance ι ˆ from (51).
 9:
 Obtain the spatial spectrum P X = 1 ι ˆ 0 , 1 ι ˆ 1 , , 1 ι ˆ N 1 T .
10:
if i iter N 2 and s method = 1 then
11:
   i method i method + 1 .
12:
  if i method = N 3 then
13:
    i method 1 .
14:
    s method 0 .
15:
  end if
16:
  Update the off-grid vector ν ˆ from (63).
17:
   Ψ ( ν ˆ ) D + Ξ diag ν ˆ I N .
18:
end if
19:
if i iter N 2 and s method = 0 then
20:
   i method i method + 1 .
21:
  if i method = N 3 then
22:
    i method 1 .
23:
    s method 1 .
24:
  end if
25:
  Update the precision of mutual coupling variance ϑ ˆ from (57).
26:
  Update the mutual coupling vector c ˆ from (46).
27:
end if
28:
i iter i iter + 1 .
29:
end while
30:
Output: the spatial spectrum P X , and the directions ( ζ + ν ) .

4. Simulation Results

Extensive simulation results have been conducted. All experiments are conducted in Matlab R2017b on a PC with a 2.9 GHz Intel Core i5 and 8 GB of RAM, and Matlab codes have been made available online at https://drive.google.com/drive/folders/1XwzbNtHXfjTrN4-wylAhGI-CwIY3u3K1. The mutual coupling effect can be generated by the following expression
c n = ( 1 + ξ c ) e j ϕ c 10 α c ( 1 + 0.5 n ) 20 , n < 5 0 , otherwise ,
where ξ c U ξ c 0.05 , 0.05 , ϕ c U ϕ c 0 , 2 π , and we use the parameter α c in dB to measure the mutual coupling effect between adjacent antennas. Additionally, we use the independent Gaussian distribution to generate the received signals, and for the m-th sample in the n-th antenna, we have s n , m CN 2 e j π 2 , 1 .
In this paper, to compare with the state-of-art direction estimation methods, we compare the proposed DFSMC method with the following algorithms:
With the simulation parameters in Table 1 and the the mutual coupling between adjacent antennas being α c = 8 dB, the spatial spectrum is given in Figure 5, where the proposed DFSMC is compared with MUSIC, CS-SBL and OGSBI methods. The estimated directions for K = 3 signals are given in Table 2. Additionally, the iteration processes of DFSMC, CS-SBL, and OGSBI methods are also given in Figure 6. With both the mutual coupling effect and off-grid, the proposed DFSMC method is advantageous in this scenario.
As shown in Figure 6, in the first 300 iterations, DFSMC method only updates the parameters μ m , Σ X and ι . Then, during the 301 to 350 iterations, the mutual coupling parameters c and ϑ are updated. For the next 50 iterations, the off-grid parameter ν is updated. With repeating the 50 iterations to update the mutual coupling parameters and the off-grid parameter, the direction error can be decreased. Moreover, as shown in Figure 6, when only the mutual coupling parameters are updated, the direction estimation performance can be not improved with the correct estimated directions. However, for the next off-grid estimation, the better performance can be achieved with the updated mutual coupling parameters.
The estimated spatial spectrum is shown in Figure 5. It can be seen that the positions of peak spectrum are closer to the ground-truth directions using the DFSMC method than the OGSBI, CS-SBL and MUSIC methods. The corresponding estimated directions are given in Table 2. When we use the following expression to measure the estimation performance
e 1 = 1 K θ ˆ θ 2 2 ,
where θ ˆ denotes the estimated directions. Then, the estimation errors (in deg) of DFSMC, OGSBI, CS-SBL and MUSIC methods can be obtained as 0.092°, 0.346°, 0.301° and 0.301°, respectively. Therefore, since the mutual coupling effect is estimated in the proposed DFSMC method, the direction estimation performance is much better than the existing methods including OGSBI, CS-SBL and MUSIC.
When the mutual coupling effect increases from α c = 8 dB to α c = 5 dB, the corresponding iteration processes and the spatial spectrums of DFSMC, CS-SBL, OGSBI and MUSIC methods are given in Figure 7 and Figure 8. The estimated directions are given in Table 3, and the estimation errors (in deg) of DFSMC, OGSBI, CS-SBL and MUSIC methods can be obtained as 0.134°, 1.024°, 1.128° and 1.495°. Compared with the direction estimation performance in the scenario α c = 8 dB, the performance in the scenario α c = 5 dB decreases for all the methods, so the mutual coupling effect has a great effect on the direction estimation performance. However, the proposed DFSMC method can also achieve much better performance than existing methods.
With the 100 trails, the direction estimation performance with different SNRs is given in Figure 9, where we use the following root-mean-square error (RMSE) expression to measure the estimation performance
e 2 = 1 K P p = 0 P 1 θ ˆ p θ p 2 2 ,
where P denote the number of trails, θ p denotes the directions in the p-th trail, and θ ˆ p denotes the estimated directions in the p-th trail. As shown in Figure 9, the proposed DFSMC method achieves the best estimation performance when the SNR of received signals is greater than 0 dB. Almost the same estimation performance is achieved by the MUSIC and CS-SBL method. However, with the mutual coupling effect, the direction grids usually cannot be estimated correctly, so the further off-grid optimization in OGSBI cannot improve the estimation performance. Figure 9 indicates that our proposed DFSMC method is very advantageous in the cases when the SNR of received signals is large.
With different mutual coupling effects between antennas, we show the simulation results in Figure 10, where the mutual coupling effect α c between adjacent antennas is from 16 dB to 2 dB. Since the proposed DFSMC method estimates the mutual coupling vector c iteratively, and DFSMC achieves the best estimation performance among the existing methods including CS-SBL, OGSBI, and MUSIC. It can be seen that with optimizing the off-grid and the mutual coupling vector, the performance of direction estimation can be improved by estimating the sparse signals in the continue domain using the DFSMC method. The computational complexity of the proposed algorithm mainly depends on step 6, step 16 and step 26. The computational complexity of step 6, step 16 and step 26 can be obtained as O ( M U 2 N + U N 2 + U 3 ) , O ( U 3 + M N U 2 ) and O ( U 3 + M U 2 N ) . Therefore, the computational complexity of the proposed algorithm can be obtained as O ( U 3 + M U 2 N + U N 2 ) . Additionally, with U N , the computational complexity can be simplified as O ( M U N 2 ) . The computational complexity of the proposed algorithm has the same order of the SBL-based algorithms, such as the OGSBI algorithm and the SBL algorithm.

5. Conclusions

The direction finding problem with the unknown mutual coupling effect has been investigated in this paper. The novel DFSMC method has been proposed to estimate the directions, the means, and variance of received signals, the mutual coupling vector, the noise variance, and the off-grid vector, et al. iteratively. Additionally, the expressions to estimate the unknown parameters have been theoretically derived using the EM method. Simulation results confirm that the proposed DFSMC method outperforms the existing direction finding methods in the ULA system with the unknown mutual coupling effect. Future work will focus on the extension of the proposed DFSMC method in the scenario with correlated signals.

Author Contributions

Conceptualization, P.C. and Z.C.; methodology, P.C. and X.Z.; software, Z.C.; validation, L.L.; formal analysis, Z.C.; investigation, P.C.; resources, P.C.; data curation, L.L.; writing—original draft preparation, P.C.; writing—review and editing, Z.C.; visualization, Z.C.; supervision, X.Z.; project administration, P.C.; funding acquisition, P.C.

Funding

This work was supported in part by the National Natural Science Foundation of China (Grant No. 61801112, 61601281), the Natural Science Foundation of Jiangsu Province (Grant No. BK20180357), the Open Program of State Key Laboratory of Millimeter Waves at Southeast University (Grant No. Z201804).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Lemma 1

When the complex vectors u and v are the functions of x , we can obtain
u H v x = u H v x 0 , u H v x 1 , , u H v x N 1 = m = 0 M 1 u m v m x 0 , , m = 0 M 1 u m v m x n ,   = , m = 0 M 1 u m x n v m + u m v m x n , = , u x n T v + u H v x n ,   = v T u x 0 , , u x n , + u H v x 0 , , v x n ,   = v T ( u ) x + u H v x .
With A and u being the function of x , we can obtain the entry in m-th row and n-th column of Au x as
Au m x n = p = 0 P 1 A m , p u p x n = p = 0 P 1 A m , p x n u p + A m , p u p x n   = u T [ A T ] m x n + [ A T ] m T u x n = A x n u + A u x n m ,
so the n-th column of Au x is
Au x n = A x n u + A u x n ,
and
Au x = A x 0 u + A u x 0 , , A x n u + A u x n , .

References

  1. Schmidt, R.O. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
  2. Schmidt, R. A Signal Subspace Approach to Multiple Emitter Location Spectrum Estimation. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 1981. [Google Scholar]
  3. Zoltowski, M.; Kautz, G.; Silverstein, S. Beamspace Root-MUSIC. IEEE Trans. Signal Process. 1993, 41, 344–364. [Google Scholar] [CrossRef]
  4. 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]
  5. Pham, G.T.; Loubaton, P.; Vallet, P. Performance analysis of spatial smoothing schemes in the context of large arrays. IEEE Trans. Signal Process. 2016, 64, 160–172. [Google Scholar] [CrossRef]
  6. Chen, P.; Cao, Z.; Chen, Z.; Yu, C. Sparse DOD/DOA estimation in a bistatic MIMO radar with mutual coupling effect. Electronics 2018, 7, 341. [Google Scholar] [CrossRef]
  7. Carlin, M.; Rocca, P.; Oliveri, G.; Viani, F.; Massa, A. Directions-of-arrival estimation through Bayesian compressive sensing strategies. IEEE Trans. Antennas Propag. 2013, 61, 3828–3838. [Google Scholar] [CrossRef]
  8. Chen, P.; Qi, C.; Wu, L.; Wang, X. Estimation of Extended Targets Based on Compressed Sensing in Cognitive Radar System. IEEE Trans.Veh. Technol. 2017, 66, 941–951. [Google Scholar] [CrossRef]
  9. Yu, Y.; Petropulu, A.P.; Poor, H.V. Measurement matrix design for compressive sensing-based MIMO radar. IEEE Trans. Signal Process. 2011, 59, 5338–5352. [Google Scholar] [CrossRef]
  10. Carlin, M.; Rocca, P.; Oliveri, G.; Viani, F.; Massa, A. Novel wideband DOA estimation based on sparse Bayesian learning with dirichlet process priors. IEEE Trans. Signal Process. 2016, 64, 275–289. [Google Scholar]
  11. Chen, P.; Qi, C.; Wu, L. Antenna placement optimisation for compressed sensing-based distributed MIMO radar. IET Radar Sonar Navig. 2017, 11, 285–293. [Google Scholar] [CrossRef]
  12. Yang, Z.; Xie, L. Enhancing sparsity and resolution via reweighted atomic norm minimization. IEEE Trans. Signal Process. 2016, 64, 995–1006. [Google Scholar] [CrossRef]
  13. Shen, Q.; Liu, W.; Cui, W.; Wu, S. Underdetermined DOA estimation under the compressive sensing framework: A review. IEEE Access 2016, 4, 8865–8878. [Google Scholar] [CrossRef]
  14. Yang, Z.; Xie, L. Exact joint sparse frequency recovery via optimization methods. IEEE Trans. Signal Process. 2016, 64, 5145–5157. [Google Scholar] [CrossRef]
  15. Tipping, M.E. Sparse Bayesian Learning and the Relevance Vector Machine. J. Mach. Learn. Res. 2001, 1, 211–244. [Google Scholar]
  16. Ji, S.; Xue, Y.; Carin, L. Bayesian compressive sensing. IEEE Trans. Signal Process. 2008, 56, 2346–2356. [Google Scholar] [CrossRef]
  17. Wu, X.; Zhu, W.; Yan, J. Direction of arrival estimation for off-grid signals based on sparse Bayesian learning. IEEE Sens. J. 2016, 16, 2004–2016. [Google Scholar] [CrossRef]
  18. Chen, P.; Zheng, L.; Wang, X.; Li, H.; Wu, L. Moving target detection using colocated MIMO radar on multiple distributed moving platforms. IEEE Trans. Signal Process. 2017, 65, 4670–4683. [Google Scholar] [CrossRef]
  19. Yang, Z.; Lihua, X.; Cishen, Z. Off-grid direction of arrival estimation using sparse Bayesian inference. IEEE Trans. Signal Process. 2013, 61, 38–43. [Google Scholar] [CrossRef]
  20. Dai, J.; Bao, X.; Xu, W.; Chang, C. Root sparse Bayesian learning for off-grid DOA estimation. IEEE Signal Process. Lett. 2017, 24, 46–50. [Google Scholar] [CrossRef]
  21. Wang, Q.; Zhao, Z.; Chen, Z.; Nie, Z. Grid evolution method for DOA estimation. IEEE Trans. Signal Process. 2018, 66, 2374–2383. [Google Scholar] [CrossRef]
  22. Zamani, H.; Zayyani, H.; Marvasti, F. An iterative dictionary learning-based algorithm for DOA estimation. IEEE Commun. Lett. 2016, 20, 1784–1787. [Google Scholar] [CrossRef]
  23. Clerckx, B.; Craeye, C.; Vanhoenacker-Janvier, D.; Oestges, C. Impact of Antenna Coupling on 2 × 2 MIMO Communications. IEEE Trans. Veh. Technol. 2007, 56, 1009–1018. [Google Scholar] [CrossRef]
  24. Zheng, Z.; Zhang, J.; Zhang, J. Joint DOD and DOA estimation of bistatic MIMO radar in the presence of unknown mutual coupling. Signal Process. 2012, 92, 3039–3048. [Google Scholar] [CrossRef]
  25. Rocca, P.; Hannan, M.A.; Salucci, M.; Massa, A. Single-snapshot DoA estimation in array antennas with mutual coupling through a multiscaling BCS strategy. IEEE Trans. Antennas Propag. 2017, 65, 3203–3213. [Google Scholar] [CrossRef]
  26. Liu, J.; Zhang, Y.; Lu, Y.; Ren, S.; Cao, S. Augmented nested arrays with enhanced DOF and reduced mutual coupling. IEEE Trans. Signal Process. 2017, 65, 5549–5563. [Google Scholar] [CrossRef]
  27. Hawes, M.; Mihaylova, L.; Septer, F.; Godsill, S. Bayesian compressive sensing approaches for direction of arrival estimation with mutual coupling effects. IEEE Trans. Antennas Propag. 2017, 65, 1357–1367. [Google Scholar] [CrossRef]
  28. Basikolo, T.; Ichige, K.; Arai, H. A novel mutual coupling compensation method for underdetermined direction of arrival estimation in nested sparse circular arrays. IEEE Trans. Antennas Propag. 2018, 66, 909–917. [Google Scholar] [CrossRef]
  29. Zhang, C.; Huang, H.; Liao, B. Direction finding in MIMO radar with unknown mutual coupling. IEEE Access 2017, 5, 4439–4447. [Google Scholar] [CrossRef]
  30. Liao, B.; Zhang, Z.G.; Chan, S.C. DOA estimation and tracking of ULAs with mutual coupling. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 891–905. [Google Scholar] [CrossRef]
  31. Termos, A.; Hochwald, B.M. Capacity benefits of antenna coupling. In Proceedings of the 2016 Information Theory and Applications (ITA), La Jolla, CA, USA, 31 January–5 February 2016; pp. 1–5. [Google Scholar]
  32. Liu, X.; Liao, G. Direction finding and mutual coupling estimation for bistatic MIMO radar. Signal Process 2012, 92, 517–522. [Google Scholar] [CrossRef]
Figure 1. The ULA system for direction finding.
Figure 1. The ULA system for direction finding.
Electronics 07 00426 g001
Figure 2. The structure of sparse matrix X .
Figure 2. The structure of sparse matrix X .
Electronics 07 00426 g002
Figure 3. The off-grid approximation for direction.
Figure 3. The off-grid approximation for direction.
Electronics 07 00426 g003
Figure 4. Graphical model of sparse bayesian learning for direction estimation.
Figure 4. Graphical model of sparse bayesian learning for direction estimation.
Electronics 07 00426 g004
Figure 5. The spatial spectrum for direction estimation ( α c = 8 dB).
Figure 5. The spatial spectrum for direction estimation ( α c = 8 dB).
Electronics 07 00426 g005
Figure 6. The estimation error with iterations ( α c = 8 dB).
Figure 6. The estimation error with iterations ( α c = 8 dB).
Electronics 07 00426 g006
Figure 7. The estimation error with iterations ( α c = 5 dB).
Figure 7. The estimation error with iterations ( α c = 5 dB).
Electronics 07 00426 g007
Figure 8. The spatial spectrum for direction estimation ( α c = 5 dB).
Figure 8. The spatial spectrum for direction estimation ( α c = 5 dB).
Electronics 07 00426 g008
Figure 9. The direction estimation performance with different SNRs.
Figure 9. The direction estimation performance with different SNRs.
Electronics 07 00426 g009
Figure 10. The direction estimation performance with different mutual coupling effect.
Figure 10. The direction estimation performance with different mutual coupling effect.
Electronics 07 00426 g010
Table 1. Simulation Parameters.
Table 1. Simulation Parameters.
ParameterValue
The signal-to-noise ratio (SNR)20 dB
The number of samples M100
The number of antennas N20
The number of signals K3
The space between antennas d0.5 wavelength
The grid space δ
The direction range[−60° 60°]
The hyperparameters b , d , f 10 3
N 1 in Algorithm 1 N 1 = 10 3
N 2 in Algorithm 1 N 2 = 300
N 3 in Algorithm 1 N 3 = 50
Table 2. Estimated directions ( α c = 8 dB).
Table 2. Estimated directions ( α c = 8 dB).
MethodsSignal 1Signal 2Signal 3
Ground-truth directions−8.268°18.128°30.428°
OGSBI−8.267°17.69°30.02°
CS-SBL−8°18°30°
MUSIC−8°18°30°
DFSMC−8.254°18.13°30.27°
Table 3. Estimated directions ( α c = 5 dB).
Table 3. Estimated directions ( α c = 5 dB).
MethodsSignal 1Signal 2Signal 3
Ground-truth direction−8.268°18.128°30.428°
OGSBI−8.222°17.29°31.99°
CS-SBL−8°17°32°
MUSIC−8°18°33°
DFSMC−8.260°18.11°30.66°

Share and Cite

MDPI and ACS Style

Chen, P.; Chen, Z.; Zhang, X.; Liu, L. SBL-Based Direction Finding Method with Imperfect Array. Electronics 2018, 7, 426. https://doi.org/10.3390/electronics7120426

AMA Style

Chen P, Chen Z, Zhang X, Liu L. SBL-Based Direction Finding Method with Imperfect Array. Electronics. 2018; 7(12):426. https://doi.org/10.3390/electronics7120426

Chicago/Turabian Style

Chen, Peng, Zhimin Chen, Xuan Zhang, and Linxi Liu. 2018. "SBL-Based Direction Finding Method with Imperfect Array" Electronics 7, no. 12: 426. https://doi.org/10.3390/electronics7120426

APA Style

Chen, P., Chen, Z., Zhang, X., & Liu, L. (2018). SBL-Based Direction Finding Method with Imperfect Array. Electronics, 7(12), 426. https://doi.org/10.3390/electronics7120426

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