Next Article in Journal
Robust Small Target Co-Detection from Airborne Infrared Image Sequences
Previous Article in Journal
3D Power Line Extraction from Multiple Aerial Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Valued 2D MUSIC Algorithm Based on Modified Forward/Backward Averaging Using an Arbitrary Centrosymmetric Polarization Sensitive Array

1
College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China
2
Beijing Institute of Computer Technology and Application, Beijing 100854, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(10), 2241; https://doi.org/10.3390/s17102241
Submission received: 9 July 2017 / Revised: 19 September 2017 / Accepted: 24 September 2017 / Published: 29 September 2017
(This article belongs to the Section Physical Sensors)

Abstract

:
Two-dimensional multiple signal classification (MUSIC) algorithm based on polarization sensitive array (PSA) has excellent performance. However, it suffers a high computational complexity due to a multitude of complex operations. In this paper, we propose a real-valued two-dimensional MUSIC algorithm based on conjugate centrosymmetric signal model, which is applicable to arbitrary centrosymmetric polarization sensitive array. The modified forward/backward averaging, which can be applied to the PSA, is presented. Hence, the eigen-decomposition analysis process and spectrum function computation are converted into real domain, prominently reducing the computational complexity. Then, the direction-of-arrival (DOA) estimation is decoupled from the polarization parameter estimation so that the four-dimensional spectral peak search process is avoided. The theoretical computational complexity is discussed and the Cramer-Rao bound (CRB) of DOA estimation is derived in this paper. The simulation results indicate that the proposed algorithm achieves superior accuracy in DOA estimation and has low computational complexity.

1. Introduction

Parameter estimation is an important area in array signal processing, such as adaptive beamforming [1,2] and direction-of-arrival (DOA) estimation [3,4,5,6,7]. Especially, the two-dimensional (2D) DOA estimation based on polarization sensitive array (PSA) shows great potential [8,9,10], which has gained considerable interests. Compared with traditional scalar array (TSA), the PSA composed of electromagnetic vector sensors (EMVSs) has more inherent merits. Besides spatial domain information, polarization domain information of the electromagnetic signals can also be measured by the PSA [11,12], which can bring greater potential capacity for parameter estimation.
During the past decades, many effective algorithms based on TSA have been extended to PSA, such as multiple signal classification (MUSIC) [13,14], estimation of signal parameters via rotational invariance techniques (ESPRIT) [15,16,17], Root-MUSIC [18,19,20], etc. With the development of research, some specific joint DOA and polarization parameter estimation algorithms for PSA have been presented, including vector cross-product based algorithms [21,22] and hyper complex MUSIC algorithms [23,24]. In almost all of the literature considering PSAs, the data model makes use of the concatenated vector of the received signal data with multi-components, resulting in the so-called “long-vector” approach [24]. An enhanced ESPRIT algorithm which improves estimation accuracy memorably and can solve the ambiguity of DOA estimation was proposed in [22]. However, this algorithm can only be applied to uniform rectangular array (URA) composed of six component electromagnetic vector sensors, leading to great restriction on the application. Moreover, the six component electromagnetic vector sensors have serious mutual coupling effect [25] and bring great economic burden. In addition, literature [26] presented an improved MUSIC algorithm based on PSA without ambiguity. It is quite applicable to the arrays with arbitrary geometries. However, the spectrum function of this algorithm is four-dimensional (4D) with respect to DOA and polarization parameters, leading to a 4D spectral peak search process. High computational load severely decelerates the parameter estimation, therefore, the reduction of computational complexity of MUSIC-based algorithm is a problem worth studying.
Root-MUSIC algorithm is a modified MUSIC algorithm which replaces the spectral peak searching by polynomial rooting, so as to reduce the computational complexity of DOA estimation [19,20]. Nevertheless, Root-MUSIC algorithm is generally used only for uniform linear array (ULA), since it requires a complicated preprocessing for other array structures. Reducing the number of searching points can also improve computation efficiency. For MUSIC-based algorithms using PSAs, decoupling the DOA estimation from polarization parameter estimation can reduce the dimensionality of the spectral peak searching and drastically reduce the number of iterations. A MUSIC-based solution using one-dimensional spectral peak searching has been proposed in [27] base on the ULA composed of cross-dipoles. On the other hand, it is well-known that processing the data in real domain is much faster than in complex domain. In subspace-based algorithms based on TSA, by utilizing the forward/backward (FB) averaging technique, the complex covariance matrix is converted into a real-valued matrix, thus the eigenvalue decomposition (EVD) can be calculated in real domain [28,29]. However, the FB averaging cannot be applied directly to the PSAs consist of EMVSs with multiple vibrators.
In this paper, a novel real-valued MUSIC algorithm based on arbitrary centrosymmetric PSA is proposed to achieve two-dimensional DOA estimation with low computational complexity and high accuracy. The merits are reflected in the following four aspects. First of all, the proposed algorithm is applicable to an arbitrary centrosymmetric PSA, even three-dimensional array geometry. Then, the modified FB averaging which can be applied to PSA is presented in the proposed algorithm to process the data in real domain, thus achieving a large reduction of computational cost. Furthermore, through decoupling the DOA estimation from the polarization parameter estimation, the 4D MUSIC spectrum function is dimension-reduced to a 2D function depending only on the DOA parameters. Hence, a joint search of both DOA and polarization parameters is avoided, which significantly reduces the number of iterations in searching the spectral peaks. The azimuth and elevation simultaneous searches (AESS) [30] can be applied to the proposed algorithm to perform a more efficient 2D DOA peak search process. Finally, electric co-centered crossed-dipole (ECCD) pairs with lower mutual coupling effect are used in the array and it will cause less compromise on the estimation performance than most other form of vector sensors. The simulation results demonstrate that the proposed algorithm can get accurate parameter estimation with a low computational complexity.
The remainder of this paper is organized as follows. Section 2 gives the array signal model used throughout this paper. Then, we explicitly describe the proposed algorithm in Section 3, and discuss the estimation performance and the complexity in Section 4. The simulation results are shown in Section 5 to verify the effectiveness of the proposed algorithm. Finally, the conclusion of this paper is given in Section 6.
Notations: We use ( · ) T , ( · ) * and ( · ) H to represent transpose, conjugate, and conjugate transpose of matrices, respectively. Denote operators and as kronecker product and Khatri-Rao (KR) product, respectively. Re ( · ) and Im ( · ) are used to indicate real part and imaginary part of complex variates, respectively. E { · } denotes the expectation operator, and det ( · ) denotes the determinant of matrices. The non-bold letters, lower-case bold letters, and upper-case bold letters denote scalars, vectors and matrices, respectively. Especially, I N denotes an N × N identity matrix, and J M denotes an M × M anti-identity matrix with ones on its anti-diagonal and zeros elsewhere.

2. Array Signal Model for DOA Estimation

2.1. Received Signal Model

Assume that K uncorrelated transverse electromagnetic plane waves, which have traveled through a homogeneous isotropic medium, impinge on the electromagnetic vector sensors. Then, the k-th unit-power completely polarized incident source has the electric-field vector e k and the magnetic-filed vector h k , which can expressed in Cartesian coordinates as [31,32,33]
a p ( θ k , ϕ k , γ k , η k ) = [ e k h k ] = [ e x , k e y , k e z , k h x , k h y , k h z , k ] = [ sin θ k cos ϕ k cos θ k cos θ k cos ϕ k sin θ k 0    sin ϕ k cos ϕ k cos θ k sin θ k cos ϕ k sin θ k cos θ k sin ϕ k    0 ] [ cos γ k sin γ k e j η k ] ,
where ϕ k [ 0 , π ] and θ k [ 0 , 2 π ) denote the elevation angle measured from the vertical z-axis and the azimuth angle measured from the positive x-axis, respectively, see in Figure 1; γ k [ 0 , π / 2 ] and η k [ π , π ) denote the auxiliary polarization angle and the polarization phase difference, respectively.
Assume that a centrosymmetric polarization sensitive array consists of M ECCD pairs is arranged in the spatial plane of the coordinate system which includes x-axis, y-axis and the original point. For brevity, we abbreviate this plane as xoy-plane. The two vibrators of each ECCD pair are deployed parallel to x-axis and y-axis, respectively. Since ECCD pairs only receive components e x , k and e y , k of electric-field vector, (1) is simplified as
a p ( θ k , ϕ k , γ k , η k ) = [ e x , k e y , k ] = [ sin θ k cos ϕ k cos θ k cos θ k cos ϕ k sin θ k ] v p ( θ k , ϕ k ) = v p , k [ cos γ k sin γ k e j η k ] h r ( γ k , η k ) = h r , k .
Note that v p , k and h r , k depend only on the source spatial angular locations and the incent signal polarization states, respectively.
The spatial steering vector of the k-th narrow-band incent signal equals
a s ( θ k , ϕ k ) = [ e j Δ k 1 , e j Δ k 2 , , e j Δ k M ] T ,
where Δ k m is the inter-vector-sensor spatial phase factor relating the k-th source to the m-th vector sensor at location ( x m , y m , z m ) and is defined as
Δ k m = 2 π ( x m sin ϕ k cos θ k + y m sin ϕ k sin θ k + z m cos ϕ k ) λ , m = 1 , , M ,
where λ means wavelength. For simplicity, we use   a p , k and a s , k to respectively represent a p ( θ k , ϕ k , γ k , η k ) and a s ( θ k , ϕ k ) in the following parts. As such, the received signal model (RSM) of array can be expressed as a long-vector form
x ( t ) = ( A s A P ) s ( t ) + n ( t ) = A s ( t ) + n ( t ) ,
A = [ a 1 , , a K ] = [ a s , 1 a p , 1 , , a s , K a p , K ] ,
A s = [ a s , 1 , , a s , K ] ,   A p = [ a p , 1 , , a p , K ] ,
s ( t ) = [ s 1 ( t ) , , s K ( t ) ] T ,   n ( t ) = [ n 1 ( t ) , , n M ( t ) ] T ,
where n m ( t ) symbolizes the 2 × 1 zero-mean additive complex white Gaussian noise vector at the m-th vector sensor.

2.2. Conjugate Centrosymmetric Signal Model

Suppose that there a 2 M × 2 M permute matrix
P ˜ e = P e I 2 ,
where P e is an M × M permute matrix, and it can be obtained readily that P ˜ e P ˜ e T = P ˜ e T P ˜ e = I 2 M . Then, a novel signal model is proposed as follow
y ( t ) = P ˜ e x ( t ) = P ˜ e ( A s ( t ) + N ( t ) ) = A e s ( t ) + n e ( t ) ,
A e = P ˜ e A = P ˜ e ( A s A P ) = ( P ˜ e A s ) A P = A s e A P ,
A s e = P ˜ e A s = [ P ˜ e a s , 1 , , P ˜ e a s , K ] = [ a s e , 1 , , a s e , K ] .
According to (3) and (4), it is obvious that the inter-vector-sensor spatial phase factors a s , k of two electromagnetic vector-sensors which are centrosymmetric around the origin of coordinate system are conjugate. Thus, for a centrosymmetric array, there must be a calculable and unique P ˜ e such that the spatial steering vector of the proposed signal model can be expressed as
A s e = { [ A s e ( T ) A s e ( B ) ] = [ A s e ( T ) ( J M / 2 A s e ( T ) ) * ] , for   even   M , [ A s e ( T ) 1 K T A s e ( B ) ] = [ A s e ( T ) 1 K T ( J ( M 1 ) / 2 A s e ( T ) ) * ] , for   odd   M ,
where 1 K denotes a K × 1 matrix of ones, A s e ( T ) denotes the top M / 2 or ( M 1 ) / 2 rows of A s e and A s e ( B ) denotes the bottom part. Denote the proposed signal model as the conjugate centrosymmetric signal model (CCSM). Furthermore, note that the matrix multiplications involving P ˜ e or J are not counted as floating operations because they merely represent the permutations of the rows or columns of the matrix being multiplied [34].
Without loss of generality, consider a uniform circular array (UCA) with 12 vector sensors and a uniform rectangular array (URA) with 3 rows and 4 columns which are shown in Figure 2. Consequently the permute matrix P ˜ e corresponding to each array is expressed as
P ˜ e UCA = [ I 6 0 6 0 6 J 6 ] I 2 , P ˜ e URA = I 12 ,
where 0 6 denotes a 6 × 6 null matrix.
Theorem 1.
The noise-subspace and spatial steering vectors of CCSM are orthogonal.
Proof. 
The covariance matrices of RSM and CCSM are
R x = E { x ( t ) x H ( t ) } = U x Λ x U x H = U x , S Λ x , S U x , S H + U x , N Λ x , N U x , N H ,
R y = E { y ( t ) y H ( t ) } = U y Λ y U y H = U y , S Λ y , S U y , S H + U y , N Λ y , N U y , N H ,
where U and Λ denote the eigenvector matrix and eigenvalue diagonal matrix, subscript x and y denote the RSM and CCSM, subscript S and N denote the signal- and noise-subspace, respectively.
According to (10) and (15), the covariance matrix of CCSM can be also expressed as
R y = E { ( P ˜ e x ( t ) ) ( P ˜ e x ( t ) ) H } = P ˜ e R x P ˜ e T = ( P ˜ e U x ) Λ x ( P ˜ e U x ) H .
Compared the expressions of R y in (16) and (17), we have
U y , S = P ˜ e U x , S , U y , N = P ˜ e U x , N .
Note that the elements in U y , S and U y , N are just a displacement of the elements in U x , S and U x , N due to the permuted matrix P ˜ e . Based on the MUSIC algorithm, the signal subspace is orthogonal to the noise subspace. The signal subspace is spanned by the steering vectors, which means that the steering vectors of sources are also orthogonal to the noise subspace, i.e.,
a k H U x , N = 0 , k = 1 , , K .
As such, we can obtain
a e , k H U y , N = ( P ˜ e a k ) H ( P ˜ e U x , N ) = a k H U x , N = 0 , k = 1 , , K
which illustrates that the noise-subspace and the steering vectors of CCSM are still orthogonal. ☐

3. Real-valued 2D DOA Estimation Based on Modified FB Averaging

3.1. The Modified FB Averaging

Considering about finite received array data, the covariance matrix of CCSM is expressed as
R ^ y = 1 L l = 1 L y l ( t ) y l H ( t ) = Y Y H ,
where L is the number of snapshots and Y = 1 / L [ y 1 ( t ) , , y L ( t ) ] is the CCSM observation matrix.
In contrast to [34,35], define the modified FB covariance matrix as
R ^ FB = 1 2 ( R ^ y + J ˜ 2 M R ^ y * J ˜ 2 M ) = Y FB Y FB H ,
where Y FB = 1 / 2 [ Y J ˜ 2 M Y *   ] is the modified FB observation matrix in which J ˜ 2 M is a 2 M × 2 M unitary permute matrix denoted as
J ˜ 2 M = J M I 2 .
It is noteworthy that, different from the original FB averaging [34], the modified FB averaging uses the unitary permute matrix J ˜ 2 M conforming to (23) instead of an anti-identity matrix. This significant change avoids the confusion of the different components of the electric-field vector of the incident signals, and allows the modified FB averaging to be applied to PSAs consist of ECCD pairs. In addition, it is straightforward to certify that
R ^ FB = J ˜ 2 M R ^ FB * J ˜ 2 M
Then, let Q be a unitary matrix which is defined as
Q = { 1 2 [ I M j I M J ˜ M j J ˜ M ] , for   even   M , 1 2 [ I M 1 0 ( M 1 ) × 2 j I M 1 0 ( M 1 ) × 2 T 2 I 2 0 ( M 1 ) × 2 T J ˜ M 1 0 ( M 1 ) × 2 j J ˜ M 1 ] , for   odd   M ,
where 0 ( M 1 ) × 2 denotes an ( M 1 ) × 2 null matrix, and J ˜ M 1 is defined similarly as (23). It is obvious that Q = J ˜ 2 M Q * .
Therefore, we can obtain the real-valued covariance matrix as
R RV = Q H R ^ FB Q .
By substituting (22) into (26), R RV can be further simplified as
R RV = 1 2 ( Q H R ^ y Q + Q H J ˜ 2 M R ^ y * J ˜ 2 M Q ) = 1 2 ( Q H R ^ y Q + Q T R ^ y * Q * ) = Re ( Q H R ^ y Q ) ,
The second equation in (27) is derived readily by using (24). (27) indicates that R RV only depends on the covariance of CCSM.
The real-valued covariance matrix R RV is also symmetric due to the following properties
R RV T = ( Q H R ^ FB Q ) T = Q T R ^ FB T Q * = Q T ( J ˜ 2 M R ^ FB H J ˜ 2 M ) Q * = Q H R ^ FB Q = R RV .
Thus, the EVD of R RV require only real-valued computation [36]. Furthermore, the eigenvalues and eigenvectors of R RV are real-valued as well.
The EVDs of R ^ FB and R RV are
R ^ FB = E Μ E H = E S Μ S E S H + E N Μ N E N H ,
R RV = V Π V T = V S Π S V S T + V N Π N V N T ,
where Μ = diag [ μ 1 , μ 2 , , μ M ] and   Π = diag [ π 1 , π 2 , , π M ] are eigenvalue diagonal matrices and E = [ e 1 , e 2 , , e M ] and V = [ v 1 , v 2 , , v M ] are eigenvector matrices. The subscript S and N represent signal- and noise-subspace, respectively.
Substituting (30) into (26), R ^ FB can be expressed as
R ^ FB = Q V Π V T Q H = Q V S Π S V S T Q H + Q V N Π N V N T Q H .
Compared (31) with (29), we have the following equations
Π = Μ ,   E N = Q V N ,   E S = Q V S .
Define Z FB as
Z FB = Q H Y FB P ,
where P is defined as
P = 1 2 [ I L j I L I L j I L ] .
Then, we can simplify R RV by using (26), (33) and (34)
R RV = Q H R ^ FB Q = ( Q H Y FB P ) ( P H Y FB H Q ) = Z FB Z FB T .
Note that Z F B is a real-valued matrix.
Partition CCSM observation matrix as Y = [ Y 1 Y 2 ] } M   rows   } M   rows for even M, or Y = [ Y 1 Y 0 Y 2 ] } ( M 1 )   rows } 2   rows   } ( M 1 )   rows for odd M. Then, substituting Q H , P , and Y FB yields
Z FB = { 1 2 [ Re ( Y 1 + J ˜ M Y 2 ) Im ( Y 1 + J ˜ M Y 2 ) Im ( Y 1 J ˜ M Y 2 ) Re ( Y 1 J ˜ M Y 2 ) ] ,  for even  M , 1 2 [ Re ( Y 1 + J ˜ M 1 Y 2 ) Im ( Y 1 + J ˜ M 1 Y 2 ) 2 Re ( Y 0 )     2 Re ( Y 0 ) Im ( Y 1 J ˜ M 1 Y 2 ) Re ( Y 1 J ˜ M 1 Y 2 ) ] ,  for odd  M .
It is easy to determine that Z FB is always real-valued. According to the previous conclusion, the matrix multiplication involving permute matrix does not contain any floating operations, but only exchanges the corresponding rows or columns. Thus, from (36), we can deduce that calculating Z FB does not involve general matrix multiplication, but is a simple construction based on Y (the CCSM data matrix) involving only additions. Therefore, it has significant consequences for processing modified FB covariance matrix.

3.2. 2D DOA Estimation with Real-Valued Spectrum Function

Applying MUSIC to the modified FB covariance matrix, the joint spatial-polarized spectrum function is
f ( θ , ϕ , γ , η ) = 1 a e H E N E N H a e .
In contrast with the joint spectrum function of the customary sensor-vector MUSIC algorithm, the signal steering vectors of CCSM and the noise-subspace of the modified FB covariance matrix are used in (37). Note that calculating the noise-subspace of modified FB covariance matrix is almost a real operation, thus reducing the substantial computational amount. However, f ( θ , ϕ , γ , η ) is a four-dimensional function with respect to both 2D DOA parameters and 2D polarization parameters. Hence, the process of parameter estimation is actually a four-dimensional iterative peak search process which leads to large computation cost. If the four-dimensional joint spectrum function can be turned into a two-dimensional function only with respect to the 2D DOA parameters, the four-dimensional iterative search will be transformed to a two-dimensional iterative search, namely, the DOA parameters can be estimated without polarization parameters estimation.
Substituting (7), (11) and (12) into (37), the reciprocal of the joint spectrum function equals
g ( θ , ϕ , γ , η ) = 1 f ( θ , ϕ , γ , η ) = a e H E N E N H a e = h r H H ( θ , ϕ ) h r ,
where
H ( θ , ϕ ) = ( a s e v p ) H E N E N H ( a s e v p )
only depends on DOA parameters, while h r only depends on polarization parameters. Therefore, we can obtain that
h r , k H H ( θ k , ϕ k ) h r , k = 0 ,
where ( θ k , ϕ k , γ k , η k ) are the DOA and polarization states of the k-th source, respectively.
According to the MUSIC algorithm, it is known that
E N H a e = E N H ( a s e , k v p ) h r , k = 0 .
Furthermore, by using (39) and (41), we have
H ( θ k , ϕ k ) h r , k = 0 .
It is obvious that the solution to (40) and (42) is unique, which indicates that H ( θ , ϕ ) is rank-deficient only when ( θ , ϕ ) equals the actual DOA, i.e., ( θ , ϕ ) { ( θ k , ϕ k ) } k = 1 K . Thus, the DOA can be estimated by solving the following optimization problem
{ θ ^ k , ϕ ^ k } = argmax θ , ϕ 1 det ( H ( θ , ϕ ) ) .
Therefore, the direction parameter estimation is successfully decoupled from the polarization estimation, thus reducing the number of variables that need to be optimized.
Then, exchange the expression of H ( θ , ϕ ) in order to estimate the DOAs only through real-valued operations. Substituting (13) and (32) into (39), it can be obtained that
H ( θ , ϕ ) = D T ( θ , ϕ ) V N V N T D ( θ , ϕ ) ,
where
D ( θ , ϕ ) = Q H a s e v p = { 2 [ Re ( a s e ( T ) ) Im ( a s e ( B ) ) ] v p , for   even   M , 2 [ Re ( a s e ( T ) ) 1 Im ( a s e ( B ) ) ] v p , for   odd   M .
Obviously, H ( θ , ϕ ) is a real-valued matrix since D ( θ , ϕ ) and V N is real-valued. From (39) and (44), it is clear that these two forms of H ( θ , ϕ ) are exactly identical. However, calculating H ( θ , ϕ ) according to (44) only involves the operations between real-valued matrices. Therefore, (44) is a much better choice to calculate the spectrum function and estimate the azimuth angle and the elevation angle.
The whole process of the proposed algorithm is summarized in Table 1.

4. Discussion

4.1. Cramer-Rao Bound (CRB) Analysis

Cramer-Rao bound (CRB) is an important parameter in discussing the performance of a system [37]. In this subsection, the CRB of 2D DOA estimation is discussed. The array received signal model is shown in Equation (5).
The Fisher information matrix (FIM) with respect to azimuth angle θ = [ θ 1 , , θ K ] , elevation angle ϕ = [ ϕ 1 , , ϕ K ] , auxiliary polarization angle γ = [ γ 1 , , γ K ] and polarization phase difference η = [ η 1 , , η K ] is written as
F = [ F θ θ F θ ϕ F ϕ θ F ϕ ϕ F θ γ F θ η F ϕ γ F ϕ η F γ θ F γ ϕ F η θ F η ϕ F γ γ F γ η F η γ F η η ] .
Here, the CRB of 2D DOA estimation is solely considered. Thus, the sub-matrix of FIM with respect to azimuth and elevation is expressed as
F DOA = [ F θ θ F θ ϕ F ϕ θ F ϕ ϕ ] .
According to literatures [38,39], the (i,j)-th element of F θ θ is given by
F θ θ ( i , j ) = L · t r { R x 1 R x θ i R x 1 R x θ j } = 2 L · t r { A ¯ θ i R s A H R x 1 A R s A ¯ θ j H R x 1 + A ¯ θ i R s A H R x 1 A ¯ θ j H R s A H R x 1 } = 2 L · Re { e i T R s A H R x 1 A R s e j e j T A ¯ θ H R x 1 A ¯ θ e i + e i T R s A H R x 1 A ¯ θ e j e j T R s A H R x 1 A ¯ θ e i } .
Hence, the sub-matrix of F DOA with respect to azimuth is
F θ θ = 2 L · Re { ( R s A H R x 1 A R s ) ( A ¯ θ H R x 1 A ¯ θ ) T + ( R s A H R x 1 A ¯ θ ) ( R s A H R x 1 A ¯ θ ) } .
Similarly, the sub-matrices with respect to elevation and the cross terms of F DOA can be expressed as
F ϕ ϕ = 2 L · Re { ( R s A H R x 1 A R s ) ( A ¯ ϕ H R x 1 A ¯ ϕ ) T + ( R s A H R x 1 A ¯ ϕ ) ( R s A H R x 1 A ¯ ϕ ) } ,
F θ ϕ = 2 L · Re { ( R s A H R x 1 A R s ) ( A ¯ ϕ H R x 1 A ¯ θ ) T + ( R s A H R x 1 A ¯ ϕ ) ( R s A H R x 1 A ¯ θ ) } ,
F ϕ θ = 2 L · Re { ( R s A H R x 1 A R s ) ( A ¯ θ H R x 1 A ¯ ϕ ) T + ( R s A H R x 1 A ¯ θ ) ( R s A H R x 1 A ¯ ϕ ) }
where t r { } denotes the trace of a matrix, e i denotes the i-th column of the identity matrix. R s is the covariance matrix of sources which is defined as R s = E { s ( t ) s H ( t ) } . A ¯ θ = [ A θ 1 , , A θ K ] and A ¯ ϕ = [ A ϕ 1 , , A ϕ K ] are the derivations of A with respect to θ i and ϕ j with i , j = 1 , 2 , , K .
Then, the CRB matrix C DOA with respect to azimuth and elevation is the inverse of F DOA , i.e.,
C DOA = F DOA 1 .
As such, the CRBs of azimuth and elevation of the k-th source are
CRB θ k = C DOA ( k , k ) ,
CRB ϕ k = C DOA ( K + k , K + k ) ,
where C DOA ( k , k ) denotes the (k,k)-th element of C DOA .

4.2. Theoretical Computational Complexity Analysis

To illustrate the computational effectiveness of the proposed algorithm, the computational complexity is explicitly discussed in this subsection comparing with LV-MUSIC algorithm [24]. The LV-MUSIC algorithm is an improved MUSIC algorithm making use of long-vector data form. It has the advantages of high resolution and superior accuracy, and can be applied to the PSAs with arbitrary structure composed of vector sensors. In view of the fairness of comparison, when considering of LV-MUSIC algorithm, assume that the polarization parameters are priori or have been estimated accurately, so that the search dimensions of the two algorithms are consistent. We abbreviate the LV-MUSIC algorithm with 2D DOAs peak searching as 2D LV-MUSIC temporarily.
The computational complexity of the proposed algorithm mainly includes the following three parts. (1) Computing the real-valued covariance matrix R RV from RSM requires about ( 16 L 4 ) M 2 + 6 L M flops; (2) The real-valued EVD (based on symmetric QR algorithm [40]) of R RV requires about 72 M 3 flops; (3) Computing the spectrum function in (44) for each searching point requires about 16 M 3 + ( 28 8 K ) M 2 + 9 M + 8 flops. Assume that there are N S P searching points in the interested search space. Thus, the computational complexity of the proposed algorithm is
C proposed = 72 M 3 + ( 16 L 4 ) M 2 + 6 L M + N S P ( 16 M 3 + ( 12 8 K ) M 2 + 29 M + 4 ) .
Similarly, the computational complexity of 2D LV-MUSIC algorithm mainly consists of following three parts. (1) Computing complex covariance matrix from RSM requires about 8 ( 3 L 2 + L 1 ) M 2 flops; (2) The EVD of complex covariance matrix requires about 288 M 3 flops; (3) Computing the spectrum function for each searching point requires about 64 M 3 + ( 56 43 K ) M 2 + 88 M + 21 flops. Under the assumption that N S P searching points are used in the interested search space, we can obtain the computational complexity of 2D LV-MUSIC algorithm is
C LV MUSIC = 288 M 3 + 8 ( 3 L 2 + L 1 ) M 2 + N S P ( 64 M 3 + ( 56 32 K ) M 2 + 88 M + 21 ) .
Figure 3 shows the complexity of the proposed algorithm and 2D LV-MUSIC algorithm versus the number of snapshots, the number of sensors and the search step of parameters. From Figure 3a–c, we can find that the theoretical computational cost of both the algorithms increases as the number of snapshot and the number of the sensors become larger, while reduces as the search step of the parameters decreases. In addition, it is obvious that the computational complexity of the proposed algorithm is always lower than 2D LV-MUSIC algorithm.

5. Simulation Results

In this section, several simulation results are presented to illustrate the performance of the proposed algorithm. A UCA and a URA as shown in Figure 2 are used in the following simulations. The radius of UCA is 2 λ and the interspace of URA is λ . Two uncorrelated far-field narrowband signals whose DOAs and polarization parameters ( θ , ϕ , γ , η ) are ( 50 ° , 10 ° , 40 ° ,   120 ° ) and ( 100 ° , 20 ° , 45 ° , 70 ° ) , respectively, are assumed to impinge on the array. The first source is left-circularly elliptically polarized, whereas the second source is right-circularly elliptically polarized. Complex additive Gaussian white noise is added into the system.
The root mean square error (RMSE) of the k-th source is defined as
RMSE k = 1 N n = 1 N ( ς ^ k , n ς k ) 2 ,   ς =   θ , ϕ ,
where N is the number of independent Monte Carlo trails. ς k is the parameter of the k-th source, and ς ^ k , n is the estimated values of ς k in the n-th trail.

5.1. The Simulation Results Distribution Scatter Plots

Figure 4 shows the simulation results of the elevation and the azimuth of the proposed algorithm and 2D LV-MUSIC algorithm in 200 trials. The signal-noise ratio (SNR) is 10 dB, and the number of snapshots is 300. The search step is set as 0.1°. The green asterisk symbols represent the actual DOA values, while the blue points represent the estimated DOA values. The small windows inside each figure are the enlarged version of the dotted areas. In Figure 4, the estimated DOA values always cluster around the actual DOA values, which demonstrate that both two algorithms derive a small estimation bias and variance. Note that in Figure 4, all the simulation results are distributed on the searching points, which is a characteristic of the spectral peak search-based algorithms. In addition, from Figure 4a, it can be seen that the maximum evaluated error of elevation and azimuth angles of proposed algorithm are less than 0.2° and 0.5°, respectively; from Figure 4b, we can see that the maximum evaluated error of elevation and azimuth angles of the 2D LV-MUSIC algorithm are less than 0.2° and 0.9°, respectively. That is, the scatter plots of proposed algorithm concentrate on a smaller area than 2D LV-MUSIC algorithm, which manifests that the proposed algorithm has a better estimation performance.

5.2. The Estimation Performance versus SNR and the Number of Snapshots

In this subsection, the RMSEs of the proposed algorithm and 2D LV-MUSIC algorithm versus SNR and the number of snapshots are compared. The azimuth angle and elevation angle are searched in the range of 0° to 360° and 0° to 90° with a step of 0.1°, respectively. 500 Monte Carlo trails are performed to acquire the RMSE.
We first examine the RMSE versus SNR where the SNR of the signals varies from −10 dB to 30 dB with a step of 5 dB and the number of snapshots is set as 300. As shown in Figure 5, it is obvious that for the estimation of the azimuth, the proposed algorithm can achieve a lower RMSE than the 2D LV-MUSIC especially for the UCA. As to the estimation of the elevation, in low SNR region, the proposed algorithm using UCA outperforms the 2D LV-MUSIC. In high SNR region, the performance of the proposed algorithm and the 2D LV-MUSIC is similar for estimating the elevation. However, as analyzed in Section 4.2, the computational complexity of the proposed algorithm is much lower than the 2D LV-MUSIC.
Then, we will examine the DOA estimation performance of the proposed algorithm and 2D LV-MUSIC algorithm versus number of snapshots by varying the number of snapshots from 50 to 500 with a step of 50. The SNR is set as 10 dB. The simulation results are plotted in Figure 6. We can easily find that for the estimation of both azimuth and elevation, the proposed algorithm can get a better estimation result than the 2D LV-MUSIC.
Figure 5 and Figure 6 verified the effectiveness of the proposed algorithm. This is mainly because when the modified FB observation matrix is obtained by the CCSM observation matrix, the size of the data matrix yields an effective doubling, which achieves an expected improvement in the estimator variances. The improvement makes the proposed algorithm have preferable robustness and higher estimation accuracy of direction finding. Furthermore, a better estimation performance of the proposed algorithm can be obtained when UCA is used.

5.3. Running Time Comparison

The runtime of the proposed algorithm and LV-MUSIC algorithm under the same conditions is compared in this subsection. The simulation results are obtained using a PC with an Inter(R) Pentium(R) G2010 2.8 GHz CPU and 8 GB RAM by running the Matlab (Ver. R2016b) codes in the same environment. Two signals are considered here. 200 Monte Carlo trails are carried out with 10 dB SNR and 300 snapshots. The azimuth angle and the elevation angle have been searched in the range of 0° to 360° and 0° to 90°, respectively. The average running time with different numbers of sensors and different search steps are shown in Table 2. It can be seen that the running time of the proposed algorithm is lower than the 2D LV-MUSIC, namely, the proposed algorithm is more efficient than the 2D LV-MUSIC in computational complexity. This is mainly because the proposed algorithm transforms the computation of the spectrum function from the complex space into the real-valued space, so that the computing time is significantly reduced at each searching point. It is worth noting that the number of searching points in the LV-MUSIC algorithm with 4D search is several orders of magnitude larger than that of the algorithm with 2D search. Thus, its running time is unacceptable, and it is not considered as a comparison here.

6. Conclusions

This paper proposes a real-valued 2D MUSIC algorithm based on CCSM data matrix, which ensures high computational efficiency and can be applied to arbitrary centrosymmetric PSA. Due to the modified FB averaging presented, the complex covariance matrix and spectrum function are converted into real domain. As a result, the computational cost is reduced significantly. By utilizing the rank deficiency theorem, the spectrum function which only depends on the DOA parameters is obtained. As such, the 2D azimuth and elevation can be estimated without estimating polarization parameters. Furthermore, in practical implementation, the search space and search step can be set reasonably under different circumstances in order to minimize the time of DOA estimation with the satisfaction of estimation accuracy. Taking UCA and URA for examples, the DOA estimation performance and computational complexity are discussed by several simulations. The simulation results validate that the proposed algorithm outperforms the 2D LV-MUSIC algorithm and has a lower running time.

Acknowledgments

This work was supported by the China NSF Grant (61671168), the Fundamental Research Funds for the Central Universities (HEUCF160807) and the Nature Science Foundation of Heilongjiang Province of China (QC2016085). The authors wish to thank the anonymous reviewers for their valuable comments on improving this paper.

Author Contributions

The main idea was proposed by Weijian Si and Yan Wang. Changbo Hou and Hong Wang conceived the experiments and provided many valuable suggestions. Yan Wang and Changbo Hou conducted experiments and collected data. Yan Wang wrote the manuscript and all the authors participated in amending the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gu, Y.; Leshem, A. Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation. IEEE Trans. Signal Process. 2012, 60, 3881–3885. [Google Scholar]
  2. Zhou, C.; Gu, Y.; He, S.; Shi, Z. A robust and efficient algorithm for coprime array adaptive beamforming. IEEE Trans. Veh. Technol. 2017. [Google Scholar] [CrossRef]
  3. Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A. Source estimation using coprime array: A sparse reconstruction perspective. IEEE Sens. J. 2017, 17, 755–765. [Google Scholar] [CrossRef]
  4. Zhou, C.; Gu, Y.; Zhang, Y.D.; Shi, Z.; Jin, T.; Wu, X. Compressive sensing coprime array signal for high-resolution DOA estimation. IET Commun. 2017, 11, 1719–1724. [Google Scholar] [CrossRef]
  5. Guo, M.; Chen, T.; Wang, B. An improved DOA estimation approach using coarray interpolation and matrix denoising. Sensors 2017, 17, 1140. [Google Scholar] [CrossRef] [PubMed]
  6. Wang, B.; Zhang, Y.D.; Wang, W. Robust DOA estimation in the presence of miscalibrated sensors. IEEE Signal Process. Lett. 2017, 24, 1073–1077. [Google Scholar] [CrossRef]
  7. Wang, B.; Zhang, Y.D.; Wang, W. Robust group compressive sensing for DOA estimation with partially distorted observations. EURASIP J. Adv. Signal Process. 2016, 2016. [Google Scholar] [CrossRef]
  8. Wu, N.; Qu, Z.; Si, W.; Jiao, S. DOA and polarization estimation using an electromagnetic vector sensor uniform circular array based on the ESPRIT algorithm. Sensors 2016, 16, 2109. [Google Scholar] [CrossRef] [PubMed]
  9. Si, W.; Zhao, P.; Qu, Z. Two-dimensional DOA and polarization estimation for a mixture of uncorrelated and coherent sources with sparsely-distributed vector sensor array. Sensors 2016, 16, 789. [Google Scholar] [CrossRef] [PubMed]
  10. Zhao, A.; Ma, L.; Ma, X.; Hui, J. An improved azimuth angle estimation method with a single acoustic vector sensor based on an active sonar detection system. Sensors 2017, 17, 412. [Google Scholar] [CrossRef] [PubMed]
  11. Yuan, X.; Wong, K.T.; Agrawal, K. Polarization estimation with a dipole-dipole pair, a dipole-loop pair, or a loop-loop pair of various orientations. IEEE Trans. Antennas Propag. 2012, 60, 2442–2452. [Google Scholar] [CrossRef]
  12. Yuan, X.; Wong, K.T.; Xu, Z.; Agrawal, K. Various compositions to form a triad of collocated dipoles/loops, for direction finding and polarization estimation. IEEE Sens. J. 2012, 12, 1763–1771. [Google Scholar] [CrossRef]
  13. Wang, L.; Yang, L.; Wang, G.; Chen, Z.; Zou, M. Uni-vector-sensor dimensionality reduction MUSIC algorithm for DOA and polarization estimation. Math. Probl. Eng. 2014, 2014, 682472. [Google Scholar] [CrossRef]
  14. Zoltowski, M.D.; Wong, K.T. Closed-form eigenstructure-based direction finding using arbitrary but identical subarrays on a sparse uniform cartesian array grid. IEEE Trans. Signal Proces. 2000, 48, 2205–2210. [Google Scholar] [CrossRef]
  15. Wong, K.T.; Zoltowski, M.D. Uni-vector-sensor ESPRIT for multisource azimuth, elevation, and polarization estimation. IEEE Trans. Antennas Propag. 1997, 45, 1467–1474. [Google Scholar] [CrossRef]
  16. Wong, K.T.; Zoltowski, M.D. Closed-form direction finding and polarization estimation with arbitrarily spaced electromagnetic vector-sensors at unknown locations. IEEE Trans. Antennas Propag. 2000, 48, 671–681. [Google Scholar] [CrossRef]
  17. Gong, J.; Lou, S.; Guo, Y. ESPRIT-like algorithm for computational-efficient angle estimation in bistatic multiple-input multiple-output radar. J. Appl. Remote Sens. 2016, 10, 025003. [Google Scholar] [CrossRef]
  18. Li, L. Root-MUSIC-based direction-finding and polarization estimation using diversely polarized possibly collocated antennas. IEEE Antennas Wirel. Propag. Lett. 2004, 3, 129–132. [Google Scholar]
  19. Qian, C.; Huang, L.; So, H.-C. Improved unitary root-MUSIC for DOA estimation based on pseudo-noise resampling. IEEE Signal Process Lett. 2014, 21, 140–144. [Google Scholar] [CrossRef]
  20. Yan, F.-G.; Shen, Y.; Jin, M. Fast DOA estimation based on a split subspace decomposition on the array covariance matrix. Signal Process. 2015, 115, 1–8. [Google Scholar] [CrossRef]
  21. Wong, K.T.; Yuan, X. “Vector cross-product direction-finding” with an electromagnetic vector-sensor of six orthogonally oriented but spatially noncollocating dipoles/loops. IEEE Trans. Signal Process. 2011, 59, 160–171. [Google Scholar] [CrossRef]
  22. Zoltowski, M.D.; Wong, K.T. ESPRIT-based 2D direction finding with a sparse uniform array of electromagnetic vector sensors. IEEE Trans. Signal Process. 2000, 48, 2195–2204. [Google Scholar] [CrossRef]
  23. Miron, S.; Le Bihan, N.; Mars, J.I. Quaternion-MUSIC for vector-sensor array processing. IEEE Trans. Signal Process. 2006, 54, 1218–1229. [Google Scholar] [CrossRef]
  24. Le Bihan, N.; Miron, S.; Mars, J.I. MUSIC algorithm for vector-sensors array using biquaternions. IEEE Trans. Signal Process. 2007, 55, 4523–4533. [Google Scholar] [CrossRef]
  25. Wang, G.; Tao, H.; Su, J.; Guo, X.; Zeng, C.; Wang, L. Mutual coupling calibration for electromagnetic vector sensor array. In Proceedings of the International Symposium on Antennas, Propagation and EM Theory (ISAPE), Xi’an, China, 22–26 October 2012; pp. 261–264. [Google Scholar]
  26. Hua, Y. A pencil-MUSIC algorithm for finding two-dimensional angles and polarizations using crossed dipoles. IEEE Trans. Antennas Propag. 1993, 41, 370–376. [Google Scholar] [CrossRef]
  27. Wang, K.; He, J.; Shu, T.; Liu, Z. Localization of mixed completely and partially polarized signals with crossed-dipole sensor arrays. Sensors 2015, 15, 31859–31868. [Google Scholar] [CrossRef] [PubMed]
  28. Chen, T.; Wu, H.; Zhao, Z. The real-valued sparse direction of arrival (DOA) estimation based on the Khatri-Rao product. Sensors 2016, 16, 693. [Google Scholar] [CrossRef] [PubMed]
  29. Akkar, S.; Gharsallah, A. Reactance domains unitary MUSIC algorithms based on real-valued orthogonal decomposition for electronically steerable parasitic array radiator antennas. IET Microw. Antennas Propag. 2012, 6, 223–230. [Google Scholar] [CrossRef]
  30. Nie, W.; Xu, K.; Feng, D.; Wu, C.Q.; Hou, A.; Yin, X. A fast algorithm for 2D DOA estimation using an omnidirectional sensor array. Sensors 2017, 17, 515. [Google Scholar] [CrossRef] [PubMed]
  31. Li, J. Direction and polarization estimation using arrays with small loops and short dipoles. IEEE Trans. Antennas Propag. 1993, 41, 379–387. [Google Scholar] [CrossRef]
  32. Dong, W.; Diao, M.; Gao, L.; Liu, L. A low-complexity DOA and polarization method of polarization-sensitive array. Sensors 2017, 17, 1170. [Google Scholar] [CrossRef] [PubMed]
  33. Gong, X.; Liu, Z.; Xu, Y.; Ahmad, M.I. Direction-of-arrival estimation via twofold mode-projection. Signal Process. 2009, 89, 831–842. [Google Scholar] [CrossRef]
  34. Linebarger, D.A.; DeGroat, R.D.; Dowling, E.M. Efficient direction-finding methods employing forward/backward averaging. IEEE Trans. Signal Process. 1994, 42, 2136–2145. [Google Scholar] [CrossRef]
  35. Pesavento, M.; Gershman, A.B.; Haardt, M. Unitary root-MUSIC with a real-valued eigendecomposition: A theoretical and experimental performance study. IEEE Trans. Signal Process. 2000, 48, 1306–1314. [Google Scholar] [CrossRef]
  36. Yan, F.; Shen, Y.; Jin, M.; Qiao, X. Computationally efficient direction finding using polynomial rooting with reduced-order and real-valued computations. J. Syst. Eng. Electron. 2016, 27, 739–745. [Google Scholar]
  37. Gu, Y.; Goodman, N.A. Information-theoretic compressive sensing kernel optimization and Bayesian Cramér-Rao bound for time delay estimation. IEEE Trans. Signal Process. 2017. [Google Scholar] [CrossRef]
  38. Weiss, A.J.; Friedlander, B. Performance analysis of diversely polarized antenna arrays. IEEE Trans. Signal Process. 1991, 39, 1589–1603. [Google Scholar] [CrossRef]
  39. Wang, Q.; Yang, H.; Chen, H.; Dong, Y.; Wang, L. A low-complexity method for two-dimensional direction-of-arrival estimation using an L-shaped array. Sensors 2017, 17, 190. [Google Scholar] [CrossRef] [PubMed]
  40. Golub, G.H.; Van Loan, C.F. Matrix Computations; JHU Press: Baltimore, MD, USA, 2012; Volume 3. [Google Scholar]
Figure 1. The coordinate system for the k-th source.
Figure 1. The coordinate system for the k-th source.
Sensors 17 02241 g001
Figure 2. Two kinds of centrosymmetric array. (a) A UCA with 12 ECCDs; (b) A URA of 3 rows and 4 columns with totally 12 ECCDs. The numbers represent the order of vector sensors.
Figure 2. Two kinds of centrosymmetric array. (a) A UCA with 12 ECCDs; (b) A URA of 3 rows and 4 columns with totally 12 ECCDs. The numbers represent the order of vector sensors.
Sensors 17 02241 g002
Figure 3. Computational complexity versus (a) the number of snapshots (for 12 sensors and 1° search step); (b) the number of sensors (for 300 snapshots and 1° search step), and (c) the search step of parameters (for 12 sensors and 300 snapshots).
Figure 3. Computational complexity versus (a) the number of snapshots (for 12 sensors and 1° search step); (b) the number of sensors (for 300 snapshots and 1° search step), and (c) the search step of parameters (for 12 sensors and 300 snapshots).
Sensors 17 02241 g003
Figure 4. Scatter plots of 2D DOA estimation for (a) the proposed algorithm and (b) 2D LV-MUSIC algorithm.
Figure 4. Scatter plots of 2D DOA estimation for (a) the proposed algorithm and (b) 2D LV-MUSIC algorithm.
Sensors 17 02241 g004
Figure 5. RMSE and CRB of (a) azimuth and (b) elevation versus SNR (for 300 snapshots).
Figure 5. RMSE and CRB of (a) azimuth and (b) elevation versus SNR (for 300 snapshots).
Sensors 17 02241 g005
Figure 6. RMSE and CRB of (a) azimuth and (b) elevation versus number of snapshots (for SNR = 10 dB).
Figure 6. RMSE and CRB of (a) azimuth and (b) elevation versus number of snapshots (for SNR = 10 dB).
Sensors 17 02241 g006
Table 1. Implementation steps of the entire approach.
Table 1. Implementation steps of the entire approach.
InputThe Received Signal Model (RSM) of Array x ( t ) = ( A s A P ) s ( t ) + n ( t ) .
Output2D DOA Estimation.
Step 1Obtain the CCSM observation matrix Y = 1 / L [ y 1 ( t ) , , y L ( t ) ] .
Step 2Compute the real-valued matrix Z FB by (36).
Step 3Compute the symmetric real-valued covariance matrix R RV = Z FB Z FB T .
Step 4Perform the EVD of R RV by (30) to get V N .
Step 5Compute the spectrum function at each searching point via H ( θ , ϕ ) = D T ( θ , ϕ ) V N V N T D ( θ , ϕ ) , and estimate 2D azimuth angle and elevation angle of K signals by { θ ^ k , ϕ ^ k } = argmax 1 / det ( H ( θ , ϕ ) ) .
Table 2. Comparison of the average running time (s) with 200 Monte Carlo trails.
Table 2. Comparison of the average running time (s) with 200 Monte Carlo trails.
Search StepNumber of SensorsProposed AlgorithmLV-MUSIC (2D Search)
0.251217.609621.4674
0.5124.85175.5413
1121.15391.4175
0.563.97734.5010
0.5186.24518.5407

Share and Cite

MDPI and ACS Style

Si, W.; Wang, Y.; Hou, C.; Wang, H. Real-Valued 2D MUSIC Algorithm Based on Modified Forward/Backward Averaging Using an Arbitrary Centrosymmetric Polarization Sensitive Array. Sensors 2017, 17, 2241. https://doi.org/10.3390/s17102241

AMA Style

Si W, Wang Y, Hou C, Wang H. Real-Valued 2D MUSIC Algorithm Based on Modified Forward/Backward Averaging Using an Arbitrary Centrosymmetric Polarization Sensitive Array. Sensors. 2017; 17(10):2241. https://doi.org/10.3390/s17102241

Chicago/Turabian Style

Si, Weijian, Yan Wang, Changbo Hou, and Hong Wang. 2017. "Real-Valued 2D MUSIC Algorithm Based on Modified Forward/Backward Averaging Using an Arbitrary Centrosymmetric Polarization Sensitive Array" Sensors 17, no. 10: 2241. https://doi.org/10.3390/s17102241

APA Style

Si, W., Wang, Y., Hou, C., & Wang, H. (2017). Real-Valued 2D MUSIC Algorithm Based on Modified Forward/Backward Averaging Using an Arbitrary Centrosymmetric Polarization Sensitive Array. Sensors, 17(10), 2241. https://doi.org/10.3390/s17102241

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