Next Article in Journal
Characterizing Dynamic Walking Patterns and Detecting Falls with Wearable Sensors Using Gaussian Process Methods
Next Article in Special Issue
DOA Finding with Support Vector Regression Based Forward–Backward Linear Prediction
Previous Article in Journal
Bodily Expression Support for Creative Dance Education by Grasping-Type Musical Interface with Embedded Motion and Grasp Sensors
Previous Article in Special Issue
An Improved DOA Estimation Approach Using Coarray Interpolation and Matrix Denoising
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Low-Complexity DOA and Polarization Method of Polarization-Sensitive Array

College of Information and Communication Engineering, Harbin Engineering University, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(5), 1170; https://doi.org/10.3390/s17051170
Submission received: 20 March 2017 / Revised: 3 May 2017 / Accepted: 12 May 2017 / Published: 20 May 2017

Abstract

:
This paper proposes a low-complexity method to estimate the direction of arrival and polarization based on the polarization sensitive array (PSA) which is composed of cross-dipoles. We built a half-quaternions model through the Cayley–Dickson form to remove the redundant information. Then, the directions of arrival (DOAs) were estimated via the root-MUSIC algorithm. Finally, the polarizations were estimated by generalized eigenvalue method. Unlike some existing searching algorithms, such as multiple signal classification (MUSIC), this method can avoid the peak searching and maintains high estimation accuracy. Moreover, we use the oblique projection operators to filter out the interference signals which are decoys of the target signal. Simulation results demonstrate the effectiveness and favorable performance of the proposed method.

1. Introduction

Polarization-sensitive arrays have received considerable attention in many disciplines, including multiple-input and multiple-output (MIMO) radar, sonar, and mobile communications [1,2,3,4,5,6]. Various array configurations have been proposed, including uniform linear array (ULA) [7] and uniform circular array (UCA) [8]. Moreover, the co-prime array [9,10,11] is also a significant configuration which consists of two uniform sparse linear subarrays. The direction of arrival (DOA) plays an important role in the processing of signals from these vector sensor arrays, which can make full use of the polarization diversity of the impinging sources during the last decades [12,13,14,15]. The most representative joint DOA and polarization estimation methods are based on the subspace algorithms, including the polarized multiple signal classification (MUSIC) method [16] and polarized signal parameters via rotational invariance technique (ESPRIT) method [17,18]. However, the polarized MUSIC method needs the peak searching which makes the computation complexity higher. Especially in polarization sensitive arrays (PSAs), the estimation of the 4-D parameter is a complex process. To solve this problem, an improved method called root-MUSIC has been proposed [19].
To reduce the computation complexity, several quaternion models have been studied in [20], including quaternions, half-quaternion, and biquaternions. In these models, quaternion is a supercomplex number which makes a multidimensional vector into a one-dimensional vector to reduce the computational effort. Miron et al. proposed a quaternion-MUSIC method [21] by using peak searching. A single polarized vector sensor is used in [22] to estimate the DOA and the polarization based on a new ESPRIT algorithm. The maximum likelihood (ML) algorithm [23] and [24] is an optimal algorithm for the DOA estimation that is much more complex than the other algorithms. Besides, various estimated DOA of coherent signal methods have been discussed, including the polarization difference smoothing (PDS) method [25] and polarization angular smoothing (PAS) [26]. In the electronic detection system, the detection of the presence of the decoy is the foundation of countering the radar active decoy. To detect the interference signal, the polarization estimation is an important processes. The decoys also have the same polarization information, which are different from target signal. In [27] , an oblique projection operator has been proposed which can filter out the decoy signals and retain the target signal.
However, it is shown in [16,21] that the MUSIC method can achieve a high accuracy by using peak searching. Hence, we need a method which can achieve a high accuracy and less complexity. To this end, we propose a joint DOA and polarization method for PSA. In the first step, we built a quaternion model to estimate the DOA. In the second step, the polarization information is estimated by the generalized eigenvector. Finally, we filter out the decoy signals via the oblique projection operator. Monte Carlo simulations verify the efficacy of the proposed method.
The remainder of this paper is organized as follows. The mathematical model and array configuration are given in Section 2. Section 3 presents the proposed DOA and polarization method. The computational complexity is discussed in Section 4. The performances of the proposed method in simulations are described in Section 5.

2. Problem Formulation

2.1. Quaternions

Quaternions are a four-dimensional hypercomplex number system, and they are widely used for the estimation of DOA and polarization parameters [21]. As an extension of complex numbers to four-dimensional (4-D) space, a quaternion q can be expressed as
q = a + b · i + c · j + d · k
where a , b , c , d R ,
i 2 = j 2 = k 2 = i j k = 1 i j = j i = k k i = i k = j j k = k j = i
Several properties of quaternions are as follows: The conjugate of quaternion q, denoted q * , is given
q * = a b · i c · j d · k
The norm of a quaternion is given by
q = q q * = q * q = a 2 + b 2 + c 2 + d 2
and its inverse is:
q 1 = q * q 2
By definition, some necessary equations can be obtained:
h · j = j · h *
q * q = q q *
q 1 q 2 q 2 q 1
( q 1 q 2 ) * = q 2 * q 1 *
where h C , q 1 , q 2 Q ( Q is the quaternion field). The Cayley–Dickson form can be expressed as:
q = α + β · j
where α = a + b · i , β = c + d · i . Through the Cayley-Dickson form, we can find that the quaternion is a supercomplex number with its real part and imaginary part all being complex numbers. This notation will be used in the following section to build the polarized signal model.

2.2. Array Configuration and Mathematical Model

PSA can acquire the polarization information besides the spatial information of the electromagnetic signals, which consists of a group of polarization-sensitive elements. Then, consider a vector sensor containing three electric and three magnetic orthogonal dipoles as discussed in [13].
p = e x e y e z h x h y h z = sin φ cos θ cos φ cos φ cos θ sin φ 0 sin θ cos θ cos φ sin φ cos θ sin φ cos φ sin θ 0 cos γ sin γ e j η
where 0 θ < π / 2 denotes the signal’s elevation angle, 0 φ < 2 π denotes the azimuth angle, 0 γ < π / 2 represents the auxiliary polarization angle, and 0 η < 2 π is the polarization phase difference. As demonstrated in [28], to estimate the DOA and polarization parameters, we need to obtain at least two of the six-component vector. In general, the vector sensor is composed of three orthogonally-oriented dipoles or loops, as depicted in Figure 1.
Note that the following analysis is adapted to all the polarized vector-sensor which can receive at least a two-dimensional electromagnetic vector. For convenience, we denote the orthogonal dipole parallel to the x-axis and y-axis, respectively, as depicted in Figure 2.
Consider K completely polarized narrow-band transverse electromagnetic (TEM) waves impinging on this array with M sensors. Then, each sensor can receive two electric field components ( e x and e y ); thus, we can get the output of the mth array element at time t as follows:
x m ( t ) = x 1 m ( t ) x 2 m ( t ) = k = 1 K e x , k · e j τ m k s k ( t ) + n m ( t ) k = 1 K e y , k · e j τ m k s k ( t ) + n m ( t )
where x 1 m and x 2 m are the output of antennas parallel to the x-axis and y-axis, respectively. e x , k and e y , k are given by Equation (11) with respect to the kth source. τ m k is the phase shift associated with the kth signal’s propagation time delay between the mth element and the phase reference point, which is given by
τ m k = 2 π d λ ( m 1 ) sin θ k cos φ k
where θ k is the elevation angles and φ k is the azimuth angles. d is the inter-sensor spacing, as depicted in Figure 1. Then, Equation (12) can be further rewritten as
x ( t ) = x 1 ( t ) x 2 ( t ) = A P 1 s ( t ) + n 1 ( t ) A P 2 s ( t ) + n 2 ( t )
where s ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s K ( t ) ] T and n ( t ) are the source and noise vectors, respectively. A is the M × K steering vector, and P 1 , P 2 are given by
P 1 = e x , 1 e x , K P 2 = e y , 1 e y , K
The basic assumptions utilized throughout this paper are listed as follows.
(1)
The K incoherent arriving signals s ( t ) are narrow band and circular signals, which means E { s s T } = 0 .
(2)
The entries of n ( t ) are white Gaussian noise and uncorrelated with each other. Noise from different sensors are independent, which means E { n 1 n 2 T } = E { s n T } = E { n s T } = E { n n T } = 0 .

3. Proposed Algorithm

In this section, a DOA and polarization estimation method is proposed based on the quaternion theory by using a uniform linear array.

3.1. Half-Quaternions Model

Through the Cayley-Dickson form, the half-quaternions model can be writen as
x = x 1 + x 2 · j
Then we can get the covariance matrix of X as follow
R = E { x x H } = E { ( x 1 + x 2 · j ) ( x 1 H j · x 2 H ) } = E { x 1 x 1 H x 1 · j · x 2 H + x 2 · j · x 1 H x 2 · j · j · x 2 H } = E { Ω + Ψ · j }
where Ω = x 1 x 1 H + x 2 x 2 H , Ψ = x 2 x 1 T x 1 x 2 T . ( · ) H is the conjugate transpose matrix, and the quaternion mathematical operations can be found in Equations (6) to (9).
Through Equation (14), we can get
E { Ψ } · j = E { ( A P 2 s + n 2 ) ( A P 1 s + n 1 ) T ( A P 1 s + n 1 ) ( A P 2 s + n 2 ) T } · j = E { A P 2 E { s s T } P 1 T A T + A P 2 E { s n 1 T } + E { n 2 s T } P 1 T A T + E { n 1 n 2 T } A P 1 E { s s T } P 2 T A T A P 1 E { s n 2 T } E { n 1 s T } P 2 T A T E { n 1 n 2 T } } · j
where ( · ) T denotes the transposed matrix. Due to the assumptions of signal as discussed in the Section 2.2, E { Ψ } · j = 0 .
Hence Equation (17) becomes
R = E { Ω } = E { ( A P 1 s + n 1 ) ( A P 1 s + n 1 ) H + ( A P 2 s + n 2 ) ( A P 2 s + n 2 ) H } = A P 1 E { s s H } P 1 H A H + A P 1 E { s n 1 H } + E { n 1 s H } P 1 H A H + E { n 1 n 1 H } + A P 2 E { s s H } P 2 H A H + A P 2 E { s n 2 H } + E { n 2 s H } P 2 H A H + E { n 2 n 2 H } = A P 1 E { s s H } P 1 H A H + A P 2 E { s s H } P 2 H A H + σ 1 2 I + σ 2 2 I
where σ 1 2 and σ 2 2 denote the noise variance. Then we define E { s s H } = d i a g { ρ 1 2 , ρ 2 2 , , ρ K 2 } with ρ k 2 is the power of the kth signal. Thus, we have
P 1 E { s s H } P 1 H = e x , 1 2 ρ 1 2 e x , 2 2 ρ 2 2 e x , K 2 ρ K 2 P 2 E { s s H } P 2 H = e y , 1 2 ρ 1 2 e y , 2 2 ρ 2 2 e y , K 2 ρ K 2
where · denotes the norm. Hence the covariance matrix R is written as
R = A ( e x , 1 2 + e y , 1 2 ) ρ 1 2 ( e x , K 2 + e y , K 2 ) ρ K 2 A H + ( σ 1 2 + σ 2 2 ) I
Note that the noise is a linear superposition of two polarization dimensions, and the signal power is equivalent to the weightings. Due to E { Ψ } 0 when the number of snapshot is limited, the information contained in E { Ψ } is redundant. Through Equation (19), the redundant information has been removed, hence this method can improve utilization of measurement data.

3.2. DOA Estimation

Without loss of generality, we assume that the azimuth angle of signals φ = 0 . Through Equation (19), the eigenvalue decomposition of covariance matrix can be obtained
R = U S Λ S U S H + U N Λ N U N H
where U S and U N denotes the signal subspace and noise subspace, respectively. The diagonal matrix Λ S and Λ n are composed of K larger eigenvalues and M K smaller eigenvalues, respectively. Invoking the MUSIC algorithm, we know
θ ^ = arg min θ a H ( θ ) U N U N H a ( θ )
where a ( θ ) is the steering vector
a ( θ ) = [ 1 , e j 2 π d sin θ / λ , , e j 2 π M d sin θ / λ ] T
Substituting Equation (24) in Equation (23), we have
a H ( θ ) U N U N H a ( θ ) = m = 0 M 1 n = 0 M 1 e 2 π m d sin θ / λ R N ( m , n ) e j 2 π n d sin θ / λ = l = ( M 1 ) M 1 e j 2 π l d sin θ / λ · r l
where R N = U N U N H is a M × M matrix, and r l is given by
r l = m n = l R N ( m , n )
Let z = e j 2 π d sin θ / λ , then Equation (25) can be written as
D ( z ) = m n = l r l z l
Thus the estimation of Equation (23) has changed into the solution of D ( z ) = 0 . As z = e j 2 π d sin θ / λ , the roots are supposed to lie on the unit circle. Then, the estimation of elevation angle is given by
θ ^ = arcsin ( λ 2 π d arg ( z ) )

3.3. Polarization Parameter Estimation

As can be seen from Equation (21), the polarization parameter become a weighting of the signal power, and hence we need reconstructing the covariance matrix as follow
X ˜ = [ X 1 ; X 2 ] R ˜ = E { X ˜ X ˜ H } = U ˜ S Λ ˜ S U ˜ S H + U ˜ N Λ ˜ N U ˜ N H
Equation (23) can be written as follow
{ γ ^ , η ^ } = arg min γ , η D θ ^ h γ , η D θ ^ h γ , η H U ˜ N U ˜ N H D θ ^ h γ , η D θ ^ h γ , η = arg min γ , η h γ , η H H ( θ ^ ) h γ , η h γ , η H D θ ^ H D θ ^ h γ , η = arg min h γ , η 0 h γ , η H H ( θ ^ ) h γ , η h γ , η H D θ ^ H D θ ^ h γ , η = arg min h γ , η H D θ ^ H D θ ^ h γ , η = 1 h γ , η H H ( θ ^ ) h γ , η h γ , η H D θ ^ H D θ ^ h γ , η
where D θ ^ is the 2 M × 2 matrix only with respect to θ .
D θ ^ = a ( θ ^ ) Ξ ( θ ^ ) , H ( θ ^ ) = D θ ^ H U ˜ N U ˜ N H D θ ^ a ( θ ^ ) = [ 1 , e j 2 π d sin θ ^ / λ , , e j 2 π M d sin θ ^ / λ ] T Ξ ( θ ^ ) = 1 0 0 cos θ ^ h γ , η = cos γ sin γ e j η
where a ( θ ^ ) is the M × 1 steering vector of the ULA. Ξ ( θ ^ ) and h γ , η can be obtained from Equation (11). Note that φ = 0 in the ULA and θ ^ has been estimated in Section 2.2.
Equation (30) can be further rewritten as
min h γ , η h γ , η H H ( θ ) h γ , η s . t . h γ , η H D θ H D θ h γ , η = 1
This optimization problem can be solved by Lagrange multiplier method as follow
Υ θ ( h γ , η , μ ) = h γ , η H H ( θ ) h γ , η + μ ( 1 h γ , η H D θ H D θ h γ , η )
Taking a derivative by h * and μ , then we can get
H ( θ ) h γ , η = μ D θ H D θ h γ , η h γ , η H D θ H D θ h γ , η = 1 h γ , η H H ( θ ) h γ , η = μ = h γ , η H H ( θ ) h γ , η h γ , η H D θ H D θ h γ , η 0
From Equation (34), we can know that min ( h γ , η H H ( θ ) h γ , η ) = min ( μ ) . μ is the generalized eigenvalue of matrix pencil { H ( θ ) , D θ H D θ } , and h γ , η is the generalized eigenvector. Hence the eigenvectors corresponding to the smallest eigenvalue is the estimation of polarization parameters.
h ^ γ , η = min { H ( θ ^ ) , D θ ^ H D θ ^ } γ ^ = arctan { h ^ γ , η ( 2 ) / h ^ γ , η ( 1 ) } η ^ = arg { h ^ γ , η ( 2 ) / h ^ γ , η ( 1 ) }
where min { · } denotes the eigenvectors corresponding to the smallest eigenvalue.

3.4. Oblique Projection Operators

In a practical application, we need to enhance signals while nulling interferences. Hence the oblique projection operator was proposed in [27]. We assume that all the doa and polarization parameters of signals have been estimated.
a s = 1 e j 2 π d sin θ 0 / λ e j 2 π M d sin θ 0 / λ ; w i = 1 e j 2 π d sin θ i / λ e j 2 π M d sin θ i / λ , i = 1 , , K
where a s is the steering vector of the target signal, and w i is the the steering vector of the ith interference ( 1 + K = K ). Then the oblique projection operators are given by
P w i = I w i ( w i H w i ) 1 w i H E s w i = s ( a s H P w i a s ) 1 a s H P w i
where P w and E s w are the orthogonal projection operators and oblique projection operators, respectively. The array output is written as
X ( t ) = E s w 1 E s w K X 1 ( t ) E s w 1 E s w K X 2 ( t )
Now the received data X ( t ) that only contains the information of the target signal has been gotten. The main steps of the proposed method are summarized in Algorithm 1 as follow.
Algorithm 1 Steps in the Proposed Method
Input: X ( 1 ) , X ( 2 ) , , X ( N )
1. obtain X according to Equation (11)
DOA Estimation:
2. Calculate the covariance matrix R = E { Ω } via Equation (8)
3. Divide R into U S and U N according to Equation (19)
4. Calculate r l according to Equation (23)
5. Calculate the roots of r l which lie on the unit circle
6. The estimates of DOA ( θ ^ ) are obtained from Equation (25)
Polarization Parameter Estimation:
7. Calculate the covariance matrix X ˜ = [ X 1 ; X 2 ] and the noise subspace U ˜ N via Equation (26)
8. Calculate D θ ^ and H ( θ ^ ) according to Equation (28)
9. Obtain the generalized eigenvectors corresponding to the smallest eigenvalue from Equation (32)
10. Estimate Polarization Parameters ( η ^ and γ ^ ) via Equation (32)
Oblique Projecting Filter
11. Find out the target signal A S and the interferences W through Equation (33)
12. Compute the oblique projection operators via Equation (34)
13. Filter out interfering signals using Equation (35)

4. Computational Complexity

The complexity of the algorithm mainly depends on two aspects: First, the the maximum likelihood estimation of the covariance matrix. Second, the process of peak searching. To demonstrate the advantage of the proposed method, we discuss the computational complexities of the proposed method and the traditional polarization MUSIC algorithm.
Table 1 presents the comparison of computational complexities of the two methods. N denotes the snapshot number, and Δ s denotes the number of searching points.
As can be seen from Table 1, the proposed method based on quaternions has a lower computational complexity. That is because the covariance matrix R ( X ) C 2 N × 2 N has become R ( Ω ) C N × N via Equation (16). As discussed in Section 3.2 and Section 3.3, we use the root algorithm and the generalized eigenvector algorithm to estimate the DOA and polarization parameters without peak searching. Furthermore, the advantage of the proposed method becomes increasingly obvious on the array which can estimate 2-D DOA information.

5. Simulation

In this section, we compare the proposed method with the standard root-MUSIC [13] and the Cramér-Rao lower bound (CRB) is also used as a benchmark [29]. Consider a ULA which contains a total of M = 8 × 2 dipoles and the inter-sensor spacing d = 0 . 5 λ . Two-hundred independent Monte Carlo trials are conducted for the following simulations, and the root mean squared error (RMSE) is chosen as a performance metric with different SNR, which are defined as
R M S E = 1 200 K i = 1 200 k = 1 K ( θ ^ k θ k ) 2
S N R = 10 log 10 ρ 2 σ 2
In the first simulation, we evaluate the DOA estimation performance of the proposed method with different SNR and snapshot number. Assume that the four far-field narrowband completely polarized electromagnetic wave sources impinge on the array as depicted in Figure 2. We set { θ 1 = 5 , γ 1 = 5 , η 1 = 11 } , { θ 2 = 16 , γ 2 = 12 , η 2 = 34 } , { θ 3 = 28 , γ 3 = 26 , η 3 = 74 } , and { θ 4 = 40 , γ 4 = 34 , η 4 = 124 } . The snapshot number and the wavelength are set to be 1000 and λ = 0 . 13 , respectively. Figure 3 depicts the performance versus SNR with snapshot number = 1000, and Figure 4 depicts the performance versus snapshot number SNR = 20 dB. It can be seen from Figure 3 and Figure 4 that the proposed method has much better performance than the standard root-MUSIC.
In the second simulation, we evaluate the polarization estimation performance of the proposed method with different SNR. As discussed in the first simulation, the proposed method has a better DOA estimation performance. The proposed method have a lower RMSE of the polarization estimation, because the polarization estimation of the two methods are all based on their DOA estimation. Hence, we only need to compare the proposed method with CRB as in Figure 5 and Figure 6. The results from Figure 3, Figure 4, Figure 5 and Figure 6 demonstrate that the proposed method yields more accurate DOA and polarization estimates than the standard root-MUSIC.
The third simulation compares the estimation performance of the oblique Projection Operators. For comparison purposes, we set a target signal and three interference signals. The target signal is parameterized by { θ 1 = 20 , γ 1 = 8 , η 1 = 25 } , and the three interference signals are parameterized by { θ 2 = 16 , γ 2 = 20 , η 2 = 60 } , { θ 3 = 30 , γ 1 = 20 , η 1 = 60 } , and { θ 4 = 8 , γ 4 = 20 , η 4 = 60 } . Note that the three interference signals have the same polarization parameters, which are different from the target signal. By comparing Figure 7 and Figure 8, we can find that the interference signals have been filtered out.

6. Conclusions

In this paper, we have proposed a low-complexity DOA and polarization estimation method based on quaternions for PSA. We firstly built a mode of quaternions to remove the redundant information, which can reduce the complexity. Then, we obtain the DOA and polarization parameters by the root-MUSIC and generalized eigenvalue methods. Finally, the interference signals are filtered out through the oblique projection operators, and then we can obtain the target signal only. Simulation results show that the proposed method has a better performance with a low-complexity.

Acknowledgments

This work was supported in part by the National Science Foundation of China under Grant 61571149. The authors would like to thank the anonymous reviewers and the associate editor for their valuable comments and suggestions, which have greatly improved the quality of this paper.

Author Contributions

The main idea was proposed by Wen Dong and Lutao Liu. Lipeng Gao and Ming Diao performed the experiments and analyzed the simulation results. Wen Dong wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Wan, L.T.; Han, G.J.; Han, G.J. PD Source Diagnosis and Localization in Industrial High-Voltage Insulation System via Multimodal Joint Sparse Representation. IEEE Trans. Ind. Electron. 2016, 63, 2506–2516. [Google Scholar] [CrossRef]
  2. Gu, C.; He, J.; Li, H.; Zhu, X. Target localization using MIMO electromagnetic. Signal Process. 2013, 93, 2103–2107. [Google Scholar] [CrossRef]
  3. Zheng, G.; Wu, B. Polarisation smoothing for coherent source direction finding with multiple-input and multiple-output electromagnetic vector sensor array. IET Signal Process. 2016, 10, 873–879. [Google Scholar] [CrossRef]
  4. 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, in press. [Google Scholar] [CrossRef]
  5. 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]
  6. Gu, Y.; Goodman, N.A.; Hong, S.; Li, Y. Robust adaptive beamforming based on interference covariance matrix sparse reconstruction. Signal Process. 2014, 96, 375–381. [Google Scholar] [CrossRef]
  7. 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, 1–23. [Google Scholar] [CrossRef] [PubMed]
  8. Goossens, R.; Rogier, H. A Hybrid UCA-RARE/Root-MUSIC Approach for 2-D Direction of Arrival Estimation in Uniform Circular Arrays in the Presence of Mutual Coupling. IEEE Trans. Antenna Propag. 2007, 55, 841–849. [Google Scholar] [CrossRef]
  9. Shen, Q.; Zhang, Y.; Amin, M.G. Generalized coprime array configurations for direction-of-arrival estimation. IEEE Trans. Signal Process. 2015, 63, 1377–1390. [Google Scholar]
  10. Sun, F.; Gao, B.; Chen, L.; Lan, P. A Low-Complexity ESPRIT-Based DOA Estimation Method for Co-Prime Linear Arrays. Sensors 2016, 16, 1367. [Google Scholar] [CrossRef] [PubMed]
  11. Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A.; Qu, F. Source Estimation Using Coprime Array: A Sparse Reconstruction Perspective. IEEE Sens. J. 2017, 17, 755–765. [Google Scholar] [CrossRef]
  12. Liu, L.; Jiang, Y.; Wan, L.; Tian, Z. Beamforming of Joint Polarization-Space Matched Filtering for Conformal Array. Sci. World J. 2013, 2013, 1–10. [Google Scholar] [CrossRef] [PubMed]
  13. Wong, K.T.; Li, L.; Zoltowski, M.D. Root-MUSIC-based direction-finding and polarization estimation using diversely polarized possibly collocated antennas. IEEE Antenna Wirel. Propag. Lett. 2004, 3, 129–132. [Google Scholar] [CrossRef]
  14. Rahamim, D.; Tabrikian, J.; Shavit, R. Source localization using vector sensor array in a multipath environment. IEEE Trans. Signal Process. 2004, 52, 3096–3103. [Google Scholar] [CrossRef]
  15. Hajian, M.; Nikookar, H.; Der Zwan, F. Branch correlation measurements and analysis in an indoor Rayleigh fading channel for polarization diversity using a dual polarized patch antenna. IEEE Microw. Wirel. Compon. Lett. 2005, 15, 555–557. [Google Scholar] [CrossRef]
  16. Cheng, Q.; Hua, Y. Performance analysis of the MUSIC and Pencil-MUSIC algorithms for diversely polarized array. IEEE Trans. Signal Process. 1996, 32, 284–299. [Google Scholar]
  17. Zoltowski, M.D.; Wong, K.T. ESPRIT-based 2-D direction finding with a sparse uniform array of electromagnetic vector sensors. IEEE Trans. Signal Process. 2000, 48, 2195–2204. [Google Scholar] [CrossRef]
  18. Li, J.; Compton, R.T. Angle and polarization esrimation using ESPRIT with a polarization sensitive array. IEEE Trans. Antenna Propag. 1991, 39, 1376–1383. [Google Scholar] [CrossRef]
  19. Silverstein, S.D.; Zoltowski, M.D. The mathematical basis for element and fourier beamspace MUSIC and root-MUSIC algorithms. Digit. Signal Process. 1991, 1, 161–175. [Google Scholar] [CrossRef]
  20. Schutte, H.D.; Wenzel, J. Hypercomplex numbers in digital signal processing. In Proceedings of the IEEE International Symposium on Circuits and Systems, New Orleans, LA, USA, 1–3 May 1990; pp. 1557–1560. [Google Scholar]
  21. Miron, S.; Bihan, N.L.; Mars, J.I. Quaternion-MUSIC for vector-sensor array processing. IEEE Trans. Signal Process. 2006, 54, 1218–1229. [Google Scholar] [CrossRef]
  22. Yuan, X. Estimating the DOA and the Polarization of a Polynomial-Phase Signal Using a Single Polarized Vector-Sensor. IEEE Trans. Signal Process. 2012, 60, 1270–1282. [Google Scholar] [CrossRef]
  23. Weiss, A.J.; Friedlander, B. Maximum Likelihood Signal Estimation for Polarization Sensitive Arrays. IEEE Trans. Antennas Propag. 1993, 41, 918–925. [Google Scholar] [CrossRef]
  24. Lee, H.; Stovall, R. Maximum Likelihood Methods for Determining the Direction of Arrival for a Single Electromagnetic Source with Unknown Polarization. IEEE Trans. Signal Process. 1994, 43, 474–479. [Google Scholar] [CrossRef]
  25. He, J.; Jiang, S.; Wang, J. Polarization difference smoothing for direction finding of coherent signals. IEEE Trans. Aerosp. Electron. Syst.. 2010, 46, 469–480. [Google Scholar] [CrossRef]
  26. Xu, Y.; Liu, Z. Polarimetric angular smoothing algorithm for an electromagnetic vector-sensor array. IET Radar Sona. Navig. 2007, 1, 230–240. [Google Scholar] [CrossRef]
  27. Behrens, R.T.; Scharf, L.L. Signal Processing Applications of Oblique Projection Operators. IEEE Trans. Signal Process. 1994, 42, 1413–1424. [Google Scholar] [CrossRef]
  28. Yuan, X.; Wong, K.T.; Xu, Z. Various compositions to form a triad of collocated dipoles/loops, for direction finding and polarization esitmation. IEEE Sens. J. 2012, 12, 1763–1771. [Google Scholar] [CrossRef]
  29. Häge, M.; Oispuu, M. DOA and polarization accuracy study for an imperfect dual-polarized antenna array. In Proceedings of the 19th European Signal Processing Conference, Barcelona, Spain, 29 August–2 September 2011; pp. 599–603. [Google Scholar]
Figure 1. (a) Three orthogonally-oriented dipoles; (b) Three orthogonally-oriented dipoles.
Figure 1. (a) Three orthogonally-oriented dipoles; (b) Three orthogonally-oriented dipoles.
Sensors 17 01170 g001
Figure 2. Uniform linear array configuration.
Figure 2. Uniform linear array configuration.
Sensors 17 01170 g002
Figure 3. RMSE versus of θ estimates signal-to-noise ratio (SNR) for two methods with fixed snapshot number 1000.
Figure 3. RMSE versus of θ estimates signal-to-noise ratio (SNR) for two methods with fixed snapshot number 1000.
Sensors 17 01170 g003
Figure 4. RMSE versus θ estimates snapshot number for two methods with fixed SNR of 20 dB.
Figure 4. RMSE versus θ estimates snapshot number for two methods with fixed SNR of 20 dB.
Sensors 17 01170 g004
Figure 5. RMSE of γ estimates versus SNR with snapshot number 1000.
Figure 5. RMSE of γ estimates versus SNR with snapshot number 1000.
Sensors 17 01170 g005
Figure 6. RMSE of η estimates versus SNR with snapshot number 1000.
Figure 6. RMSE of η estimates versus SNR with snapshot number 1000.
Sensors 17 01170 g006
Figure 7. The DOA and polarization estimation before the filtering with fixed SNR 20.
Figure 7. The DOA and polarization estimation before the filtering with fixed SNR 20.
Sensors 17 01170 g007
Figure 8. The DOA and polarization estimation after the filtering with fixed SNR 20.
Figure 8. The DOA and polarization estimation after the filtering with fixed SNR 20.
Sensors 17 01170 g008
Table 1. Comparison of computational complexity of two methods.
Table 1. Comparison of computational complexity of two methods.
MethodsCovariance MatrixPeak Search
Proposed M 2 N without
MUSIC ( 2 M ) 2 N ( 4 M + 4 ) × 2 M × Δ s

Share and Cite

MDPI and ACS Style

Dong, W.; Diao, M.; Gao, L.; Liu, L. A Low-Complexity DOA and Polarization Method of Polarization-Sensitive Array. Sensors 2017, 17, 1170. https://doi.org/10.3390/s17051170

AMA Style

Dong W, Diao M, Gao L, Liu L. A Low-Complexity DOA and Polarization Method of Polarization-Sensitive Array. Sensors. 2017; 17(5):1170. https://doi.org/10.3390/s17051170

Chicago/Turabian Style

Dong, Wen, Ming Diao, Lipeng Gao, and Lutao Liu. 2017. "A Low-Complexity DOA and Polarization Method of Polarization-Sensitive Array" Sensors 17, no. 5: 1170. https://doi.org/10.3390/s17051170

APA Style

Dong, W., Diao, M., Gao, L., & Liu, L. (2017). A Low-Complexity DOA and Polarization Method of Polarization-Sensitive Array. Sensors, 17(5), 1170. https://doi.org/10.3390/s17051170

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