Next Article in Journal
A Soft-Switching SEPIC with Multi-Output Sources
Next Article in Special Issue
Full-Diversity QO-STBC Technique for Large-Antenna MIMO Systems
Previous Article in Journal / Special Issue
MIMO Channel-State Estimation in the Presence of Partial Data and/or Intermittent Measurements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fully Quaternion-Valued Adaptive Beamforming Based on Crossed-Dipole Arrays

Communications Research Group, Department of Electronic and Electrical Engineering, University of Sheffield, Sheffield S1 4ET, UK
*
Author to whom correspondence should be addressed.
Submission received: 4 March 2017 / Revised: 15 April 2017 / Accepted: 25 April 2017 / Published: 27 April 2017
(This article belongs to the Special Issue Smart Antennas and MIMO Communications)

Abstract

:
Based on crossed-dipole antenna arrays, quaternion-valued data models have been developed for both direction of arrival estimation and beamforming in the past. However, for almost all the models, and especially for adaptive beamforming, the desired signal is still complex-valued as in the quaternion-valued Capon beamformer. Since the complex-valued desired signal only has two components, while there are four components in a quaternion, only two components of the quaternion-valued beamformer output are used and the remaining two are simply discarded, leading to significant redundancy in its implementation. In this work, we consider a quaternion-valued desired signal and develop a fully quaternion-valued Capon beamformer which has a better performance and a much lower complexity. Furthermore, based on this full quaternion model, the robust beamforming problem is also studied in the presence of steering vector errors and a worst-case-based robust beamformer is developed. The performance of the proposed methods is verified by computer simulations.

1. Introduction

Electromagnetic (EM) vector sensor arrays can track the direction of arrival (DOA) of impinging signals as well as their polarization. A crossed-dipole sensor array—firstly introduced in [1] for adaptive beamforming—works by processing the received signals with a long polarization vector. Based on such a model, the beamforming problem was studied in detail in terms of output signal-to-interference-plus-noise ratio (SINR) [2]. In [3,4], further detailed analysis was performed showing that the output SINR is affected by DOA and polarization differences.
Since there are four components for each vector sensor output in a crossed-dipole array, a quaternion model instead of long vectors has been adopted in the past for both adaptive beamforming and direction of arrival estimation [5,6,7,8,9,10,11,12]. In [13], the well-known Capon beamformer was extended to the quaternion domain and a quaternion-valued Capon (Q-Capon) beamformer was proposed with the corresponding optimum solution derived.
However, in most of the beamforming studies, the signal of interest (SOI) is still complex-valued; i.e., with only two components: in-phase (I) and quadrature (Q). Since the output of a quaternion-valued beamformer is also quaternion-valued, only two components of the quaternion are used to recover the SOI, which leads to redundancy in both calculation and data storage. However, with the development of quaternion-valued communications [14,15,16], it is very likely that in the future we will have quaternion-valued signals as the SOI, where two traditional complex-valued signals with different polarisations arrive at the antenna array with the same DOA. In such a case, a full quaternion-valued array model is needed to compactly represent the four-component desired signal and also make sure the four components of the quaternion-valued output of the beamformer are fully utilised. In this work, we develop such a model and propose a new quaternion-valued Capon beamformer, where both its input and output are quaternion-valued.
Based on the proposed full quaternion model, we further study the robust adaptive beamforming problem in the presence of steering vector errors. In the past, many methods have been proposed to improve the robustness of an adaptive beamformer, such as diagonal loading [17,18] and those based on the optimization of worst-case performance [19,20]. In [11], the worst-case-based method is extended to the quaternion-valued case for crossed-dipole arrays; however, it is not a full quaternion model, since the desired signal is still complex-valued. In this work, we extend the worst-case optimisation approach to the full quaternion model and a worst-case-based full-quaternion robust adaptive beamforming method is proposed.
This paper is structured as follows. The full quaternion-valued array model is introduced in Section 2, and the proposed full quaternion-valued Capon beamformer is developed in Section 3. The robust beamforming problem is studied in Section 4 based on the worst-case optimisation approach. Simulation results are presented in Section 5, and conclusions are drawn in Section 6.

2. Quaternion Model for Crossed-Dipole Array

2.1. Basics of Quaternion

Quaternion was introduced by W. R. Hamilton [21]. A quaternion is constructed by four components, with one real part and three imaginary parts, defined as
q = q a + i q b + j q c + k q d
where i, j, k are three basic imaginary units. The multiplication principle among such units is
i 2 = j 2 = k 2 = i j k = 1 ,
and
i j = j i = k , k i = i k = j , j k = k j = i
The conjugate q * of q is q * = q a i q b j q c k q d .
Note that the multiplication of quaternions is generally non-commutative: for q 1 , q 2 H , q 1 q 2 q 2 q 1 , where H represents the quaternion domain. However, if one of the factors (e.g., q 0 ) is real, then we have q 0 q 1 = q 1 q 0 .

2.1.1. Quaternion Vector and Matrix

Quaternion vectors and matrices are quite different from real and complex ones. For quaternion vectors, the space can be denoted by two different bases—left-spanned and right-spanned—given by
s p a n L { u 1 , u 2 , . . . u M } , s p a n R { u 1 , u 2 , . . . u M } .
Here we use q L to denote the left linear combination case and q R to denote the right linear combination case. Note that q = q L = q R ,
q L = m = 1 M τ m u m q R = m = 1 M u m γ m
where τ m , γ m H , for m = 1 , 2 , . . . , M .
Similarly, the eigendecomposition operation for quaternion matrices also differs from real and complex ones. We can also define the left eigendecomposition Q L and right eigendecomposition Q R for a Hermitian matrix Q H N × N , ( Q = Q H , where { } H is the Hermitian transpose, a combination of the quaternion-valued conjugate and transpose operations)
Q L = n = 1 N α n u n u n H Q R = n = 1 N u n u n H β n
where α and β are the left and right eigenvalues [5,22,23]. In this paper, only right-spanned spaces and left eigendecomposition are employed.

2.1.2. The Gradient for a Quaternion Function

Let l ( w ) be a scalar function of the quaternion-valued vector w H N × 1 , with w = [ w 1 . . . w n . . . w N ] T , and w n = w n a + w n b i + w n c j + w n d k .
Its gradient with respect to w is defined by [24]
w l = [ w 1 l . . . w n l . . . w N l ] T
where
w n l = 1 4 ( w n a l + w n b l i + w n c l j + w n d l k )
w n * l = 1 4 ( w n a l w n b l i w n c l j w n d l k )

2.2. Model for Crossed-Dipole Array

A quaternion number q = q a + i q b + j q c + k q d can be conveniently expressed as a combination of two complex numbers c 1 = q a + j q c and c 2 = q b + j q d , as follows
q = c 1 + i c 2 = ( q a + j q c ) + i ( q b + j q d ) = q a + i q b + j q c + k q d
We will use this form later to represent our quaternion-valued signal of interest.
Consider a uniform linear array (ULA) with N crossed-dipole sensors, as shown in Figure 1, where the adjacent vector sensor spacing d equals half wavelength, and the two components of each crossed-dipole are parallel to x- and y-axes, respectively. A quaternion-valued narrowband signal s 0 ( t ) impinges upon the vector sensor array among other M uncorrelated quaternion-valued interfering signals { s m ( t ) } m = 1 M , with background noise n ( t ) . s 0 ( t ) can be decomposed into
s 0 ( t ) = s 01 ( t ) + i s 02 ( t )
where s 01 ( t ) and s 02 ( t ) are two complex-valued sub-signals with the same DOA but different polarizations.
Assume that all signals are ellipse-polarized. Here we use θ to denote the azimuth angle, ϕ to denote the elevation angle, and γ and η the polarization parameters of the signal [11]. The parameters, including DOA and polarization of the m-th signal, are denoted by ( θ m , ϕ m , γ m 1 , η m 1 ) for the first sub-signal and ( θ m , ϕ m , γ m 2 , η m 2 ) for the second sub-signal. Each crossed-dipole sensor receives signals both in the x and y sub-arrays.
For signal s m ( t ) , the corresponding received signals in the x and y sub-arrays are, respectively, given by [11]:
x ( t ) = a m 1 p x m 1 s m 1 ( t ) + a m 2 p x m 2 s m 2 ( t ) y ( t ) = a m 1 p y m 1 s m 1 ( t ) + a m 2 p y m 2 s m 2 ( t )
where x ( t ) represents the received part in the x-sub-array and y ( t ) represents the part in the y-sub-array, and ( p x m 1 , p y m 1 ) and ( p x m 2 , p y m 2 ) are the polarizations of the two complex sub-signals in x and y directions, respectively, given by,
p x m 1 = cos γ m 1 p y m 1 = cos ϕ m sin γ m 1 e j η m 1 p x m 2 = cos γ m 2 p y m 2 = cos ϕ m sin γ m 2 e j η m 2 , when θ m = π 2
Note that a m 1 and a m 2 are the steering vectors for both sub-signals, which are the same since the two sub-signals share the same DOA ( θ m , ϕ m ) .
a m 1 = [ 1 , e j 2 π sin θ m sin ϕ m λ , . . . , e j ( N 1 ) 2 π sin θ m sin ϕ m λ ] T a m 2 = [ 1 , e j 2 π sin θ m sin ϕ m λ , . . . , e j ( N 1 ) 2 π sin θ m sin ϕ m λ ] T
A quaternion model can be constructed by combining the two parts as below:
q m ( t ) = x ( t ) + i y ( t ) = a m 1 ( p x m 1 + i p y m 1 ) s m 1 ( t ) + a m 2 ( p x m 2 + i p y m 2 ) s m 2 ( t )
= b m 1 s m 1 ( t ) + b m 2 s m 2 ( t )
where { b m 1 , b m 2 } H N × 1 can be considered as the composite quaternion-valued steering vector. Combining all source signals and the noise together, the result is given by:
q ( t ) = m = 0 M ( b m 1 s m 1 ( t ) + b m 2 s m 2 ( t ) ) + n q ( t )
where n q ( t ) = n x ( t ) + i n y ( t ) is the quaternion-valued noise vector consisting of the two sub-array noise vectors n x ( t ) and n y ( t ) .

3. The Full Quaternion-Valued Capon Beamformer

To recover the SOI among interfering signals and noise, the basic idea is to keep a unity response to the SOI at the beamformer output and then reduce the power/variance of the output as much as possible [25,26]. The key to constructing such a Capon beamformer in the quaternion domain is to design an appropriate constraint to make sure the quaternion-valued SOI can pass through the beamformer with the desired unity response.
Again note that the quaternion-valued SOI can be expressed as a combination of two complex sub-signals. To construct such a constraint, one choice is to make sure the first complex sub-signal of the SOI passes through the beamformer and appears in the real and j components of the beamformer output, while the second complex sub-signal appears in the i and k components of the beamformer output. Then, with a quaternion-valued weight vector w, the constraint can be formulated as
w H C = f
where C = [ b 01 b 02 ] , f = [ 1 i ] .
With this constraint, the beamformer output z ( t ) is given by
z ( t ) = w H q ( t ) = s 01 ( t ) + i s 02 ( t ) s 0 ( t ) + w H n q ( t ) + m = 1 M w H [ b m 1 s m 1 ( t ) + b m 2 s m 2 ( t ) ]
Clearly, the quaternion-valued SOI has been preserved at the output with the desired unity response.
Now, the full-quaternion Capon (full Q-Capon) beamformer can be formulated as
min w H R w subject to w H C = f
where
R = E { q ( t ) q H ( t ) }
Applying the Lagrange multipliers method, we have
l ( w , λ ) = w H R w + ( w H C f ) λ H + λ ( C H w f H )
where λ is a quaternion-valued vector of the Lagrange multipliers.
The minimum can be obtained by setting the gradient of (22) with respect to w * equal to a zero vector [24]. It is given by
w * l ( w , λ ) = 1 2 R w + 1 2 C λ H = 0
Considering all the constraints above, we obtain the optimum weight vector w o p t as follows
w o p t = R 1 C ( C H R 1 C ) 1 f H
A detailed derivation for the quaternion-valued optimum weight vector can be found at the Appendix A.
In the next, we give a brief analysis to show that by this optimum weight vector, the interference part at the beamformer output z ( t ) in (19) has been suppressed effectively.
Expanding the covariance matrix, we have
R = E { q ( t ) q H ( t ) } = R i + n + σ 1 2 b 01 b 01 H + σ 2 2 b 02 b 02 H
where σ 1 2 , σ 2 2 are the power of the two sub-signals of SOI and R i + n denotes the covariance matrix of interferences plus noise. Using (24), we have
w o p t = R i + n 1 C β
where β = ( C H R i + n C ) 1 f H H 2 × 1 is a quaternion vector.
Applying left eigendecomposition for quaternion matrices,
R i + n = n = 1 N α n u n u n H
with α 1 . . . α M 2 > α M 1 = . . . = α N = 2 σ 0 2 R , where 2 σ 0 2 denotes the noise power.
With sufficiently high interference-to-noise ratio (INR), the inverse of R i + n can be approximated by
R i + n 1 n = M + 1 N 1 2 σ 0 2 u n u n H
Then we have
w o p t = n = M + 1 N 1 2 σ 0 2 u n u n H C β = n = M + 1 N u n ρ n
where ρ n = 1 2 σ 0 2 u n H C β . Since u n H H 1 × N , C H N × 2 , and β H 2 × 1 , ρ n is a quaternion-valued scalar. Clearly, w o p t is the right linear combination of { u M + 1 , u M + 2 , . . . , u N } , and w s p a n R { u M + 1 , u M + 2 , . . . , u N } .
For those M interfering signals, their quaternion steering vectors belong to the space right-spanned by the related M eigenvectors; i.e., b m 1 , b m 2 s p a n R { u 1 , u 2 , . . . , u M } . As a result,
w o p t H b m 1 0 , w o p t H b m 2 0 , m = 1 , 2 , . . . , M
which shows that the beamformer has eliminated the interferences effectively.

4. Worst-Case-Based Robust Adaptive Beamforming

4.1. Worst-Case Constrained Algorithm

The proposed full quaternion beamformer is based on accurate steering vectors of incident array signals. However, steering vector mismatch usually exists in array processing as a result of DOA or polarization mismatch and various model errors. With such a mismatch, there will be an error vector e ( | e | ε , where ε is the upper bound of its norm) between the assumed quaternion-valued steering vector b m and the actual quaternion-valued steering vector b ¯ m , which can be expressed as [11]
b ¯ m = b m + e
According to (31), the actual quaternion-valued steering vector could be any vector within the multi-dimensional sphere, centered at b m with a radius ε .
In order to achieve a robust response, the beamformer can be constrained to have a response greater than unity for all steering vectors within such a sphere. For the proposed quaternion-valued beamformer in the last section, the constraint can be formulated as
min w H R w s . t . min | w H b ¯ 01 | 1 , min | w H b ¯ 02 | 1
According to (31), for the first constraint (corresponding to b ¯ 01 ), we have
| w H b ¯ 01 | = | w H b 01 + w H e | | w H b 01 | | w H e | .
Since | e | ε , we further have
| w H b ¯ 01 | | w H b 01 | | w H e | | w H b 01 | ε | | w | |
From (34), it can be derived that
min | w H b ¯ 01 | = | w H b 01 | ε | | w | | 1 | w H b 01 | 1 + ε | | w | |
Now the constraint becomes
min w H R w s . t . | w H b 01 | 1 + ε | | w | |
We should ensure that the response for the first sub-signal s 01 ( t ) is constrained in the real and j domain and the response for s 02 ( t ) in the i and k domain. For a given quaternion-valued vector w , the final output power w H R w keeps unchanged if w undergoes any phase shift. That is to say, we can always find a quaternion-valued vector w to make the product w H b 01 real-valued. Then, the constraint can be changed to
min w H R w s . t . R e { w H b 01 } 1 + ε | | w | | I m ( i ) { w H b 01 } = 0 I m ( j ) { w H b 01 } = 0 I m ( k ) { w H b 01 } = 0
We can transform the constraint of the b ¯ 02 part in the same way. Then, the overall formulation for the worst-case based method can be described as
min w H R w s . t . R e { w H b 01 } 1 + ε | | w | | I m ( i ) { w H b 01 } = 0 I m ( j ) { w H b 01 } = 0 I m ( k ) { w H b 01 } = 0 R e { w H b 02 } = 0 I m ( i ) { w H b 02 } 1 + ε | | w | | I m ( j ) { w H b 02 } = 0 I m ( k ) { w H b 02 } = 0
This is a convex problem, and can be solved by the second-order cone (SOC) programming method. We refer to this method as the full quaternion worst-case constraint beamformer (FQWCCB).

4.2. SOC Implementation of FQWCCB

By Cholesky decomposition of R , we have [27]
R = Q Q H
The output power w H R w can be rewritten as w H Q Q H w , which is the l-2 norm of the vector w H Q . To simplify the calculation, the constraint can be transformed into
min | | w H Q | | s . t . R e { w H b 01 } 1 + ε | | w | | I m ( i ) { w H b 01 } = 0 I m ( j ) { w H b 01 } = 0 I m ( k ) { w H b 01 } = 0 R e { w H b 02 } = 0 I m ( i ) { w H b 02 } 1 + ε | | w | | I m ( j ) { w H b 02 } = 0 I m ( k ) { w H b 02 } = 0
Note that Q , w , b 01 , and b 02 are quaternion-valued, and we cannot use quaternion-valued variables directly in SOC programming. As a result, it is necessary to convert the quaternion-valued matrices or vectors into real-valued ones.
We can rewrite Q , w , b 01 , and b 02 as follows
Q = Q 1 + Q 2 i + Q 3 j + Q 4 k w = w 1 + w 2 i + w 3 j + w 4 k b 01 = b 01 , 1 + b 01 , 2 i + b 01 , 3 j + b 01 , 4 k b 02 = b 02 , 1 + b 02 , 2 i + b 02 , 3 j + b 02 , 4 k
where Q l , w l , b 01 , l , and b 02 , l ( 1 l 4 ) are all real-valued.
We can also define a new form for such quaternion-valued matrices and vectors
Q ¯ Q 1 Q 2 Q 3 Q 4 Q 2 Q 1 Q 4 Q 3 Q 3 Q 4 Q 1 Q 2 Q 4 Q 3 Q 2 Q 1 w ¯ H [ w 1 T w 2 T w 3 T w 4 T ] b 01 1 [ b 01 , 1 T b 01 , 2 T b 01 , 3 T b 01 , 4 T ] T b 01 2 [ b 01 , 2 T b 01 , 1 T b 01 , 4 T b 01 , 3 T ] T b 01 3 [ b 01 , 3 T b 01 , 4 T b 01 , 1 T b 01 , 2 T ] T b 01 4 [ b 01 , 4 T b 01 , 3 T b 01 , 2 T b 01 , 1 T ] T b 02 1 [ b 02 , 1 T b 02 , 2 T b 02 , 3 T b 02 , 4 T ] T b 02 2 [ b 02 , 2 T b 02 , 1 T b 02 , 4 T b 02 , 3 T ] T b 02 3 [ b 02 , 3 T b 02 , 4 T b 02 , 1 T b 02 , 2 T ] T b 02 4 [ b 02 , 4 T b 02 , 3 T b 02 , 2 T b 02 , 1 T ] T
It can be verified that
| | w H Q | | = | | w ¯ H Q ¯ | | R e { w H b 01 } = w ¯ H b 01 1 , R e { w H b 02 } = w ¯ H b 02 1 I m ( i ) { w H b 01 } = w ¯ H b 01 2 , I m ( i ) { w H b 02 } = w ¯ H b 02 2 I m ( j ) { w H b 01 } = w ¯ H b 01 3 , I m ( j ) { w H b 02 } = w ¯ H b 02 3 I m ( k ) { w H b 01 } = w ¯ H b 01 4 , I m ( k ) { w H b 02 } = w ¯ H b 02 4
Based on the above real-valued matrices and vectors, the formulation changes to
min | | w ¯ H Q ¯ H | | s . t . w ¯ H b 01 1 1 + ε | | w | | , w ¯ H b 02 1 = 0 w ¯ H b 01 2 = 0 , w ¯ H b 02 2 1 + ε | | w | | w ¯ H b 01 3 = 0 , w ¯ H b 02 3 = 0 w ¯ H b 01 4 = 0 , w ¯ H b 02 4 = 0
Solving the above convex optimization problem, we can obtain a real-valued weight vector w ¯ H R 1 × 4 N . The quaternion-valued weight vector can be recovered from the four corresponding elements of w ¯ H by (41) and (42).

4.3. Complexity Analysis

In this section, we make a comparison of the computation complexity between the Q-Capon beamformer in [13] and our proposed full Q-Capon beamformer. To deal with a quaternion-valued signal, the Q-Capon beamformer has to process the two complex sub-signals separately to recover the desired signal completely, which means we need to apply the beamformer twice for a quaternion-valued SOI. However, it is not needed for the full Q-Capon beamformer, and the SOI is recovered directly by applying the beamformer once.
For the Q-Capon beamformer, the weight vector is calculated by w = R 1 a 0 ( a 0 H R 1 a 0 ) 1 , where a 0 is the steering vector for the complex-valued SOI. As an example, we use Gaussian elimination to calculate the matrix inversion R 1 and 1 3 ( N 3 N ) quaternion-valued multiplications are needed, equivalent to 16 3 ( N 3 N ) real-valued multiplications. Additionally, R 1 a 0 requires 16 N 2 real-valued multiplications, while 16 ( N 2 + N ) real multiplications are needed for ( a 0 H R 1 a 0 ) 1 . In total, 16 3 N 3 + 32 N 2 + 80 3 N real multiplications are needed. When processing a quaternion-valued signal, this number will be doubled and the total number of real multiplications becomes 32 3 N 3 + 64 N 2 + 160 3 N .
For the proposed full Q-Capon beamformer, in addition to calculating R 1 , 32 N 2 real multiplications are needed to calculate R 1 C and 32 M 2 + 32 M + 96 real multiplications for ( C H R 1 C 0 ) 1 f . In total, the number of real-valued multiplications is 16 3 M 3 + 64 M 2 + 272 3 M + 96 , which is roughly half that of the Q-Capon beamformer.

5. Simulations Results

In this section, simulations results are provided in terms of the resultant beam pattern, output SINR, and robustness against steering vector mismatch. In the beam pattern part, both 2-D and 3-D beam patterns are presented for the proposed full Q-Capon beamformer to show the unity response to the desired signal and effective suppression to interferences. The beam pattern of the FQWCCB is also provided to show its tolerance around the desired DOA and polarization region. In the output SINR performance part, the full Q-Capon beamformer, FQWCCB, and the Q-Capon beamformer are compared in two scenarios: one without steering vector mismatch, and one with steering vector mismatch. A solid-line is also displayed as the optimal (ideal) beamforming result. The error constant ε for FQWCCB is set to 1 . 3 . In the robustness part, the beamformers are compared in terms of output SINR versus the snapshots number with 1° and 5° DOA and polarization mismatch errors, respectively.

5.1. Beam Pattern

We consider 10 pairs of cross-dipoles with half wavelength spacing. All signals are assumed to arrive from the same plane of θ = 90 °, and all interferences have the same polarization parameter γ = 60 °. For the SOI, the two sub-signals are set to (90°, 1.5°, 90°, 45°) and (90°, 1.5°, 0°, 0°), with interferences coming from (90°, 30°, 60°, −80°), (90°, −70°, 60°, 30°), (90°, −20°, 60°, 70°), (90°, 50°, 60°, −50°), respectively. The background noise is zero-mean quaternion-valued Gaussian. The power of SOI and all interfering signals are set equal, and SNR (INR) is 20 dB.
Figure 2 shows the resultant 3-D beam pattern by the proposed full Q-Capon beamformer, where the interfering signals from ( ϕ , η ) = (30°, −80°), (−70°, 30°), (−20°, 70°), and (50°, −50°) have all been effectively suppressed, while the gain of SOI from ϕ = 1 . 5 ° stays almost constant.
Figure 3 and Figure 4 show the 2-D beam patterns for η = 80 ° , 30 ° , 70 ° , 50 ° , respectively, where deep nulls have been formed at locations ϕ = 30 ° , 70 ° , 20 ° , 50 ° .
The 3-D beam pattern for the proposed FQWCCB is provided in Figure 5. We can still see that the interferences from ( ϕ , η ) = (30°, −80°), (−70°, 30°), (−20°, 70°) and (50°, −50°) are suppressed, but the suppression is not as deep as in the case of the full Q-Capon beamformer; moreover, a large magnitude response are formed around the direction of the desired signal.

5.2. Output SINR Performance

In the following, the output SINR performance of the two Capon beamformers (full Q-Capon and Q-Capon) is studied with the DOA and polarization (90°, 1.5°, 90°, 45°) and (90°, 1.5°, 0°, 0°) for SOI and (90°, 30°, 60°, −80°), (90°, −70°, 60°, 30°), (90°, −20°, 60°, 70°), (90°, 50°, 60°, −50°) for interferences, respectively. Again, we have set SNR = INR = 20 dB. All results are obtained by averaging 1000 Monte-Carlo trials.
Firstly, we consider a scenario without steering vector mismatch. Figure 6 shows the output SINR performance versus SNR with 100 snapshots, where the solid-line is for the optimal beamformer (i.e., with infinite number of snapshots). Note that with only 100 snapshots, there will be non-negligible data model errors. It can be seen that for most of the input SNR range, the FQWCCB has the best output SINR performance versus the input SNR, while the full Q-Capon has a worse performance, but it still outperforms the Q-Capon beamformer. In the lower SNR range, FQWCCB and the full Q-Capon beamformer have a similar result, but the gap becomes greater from input SNR = 10 dB onwards. Additionally, in the lower range, the proposed full Q-Capon beamformer has a better performance than the Q-Capon beamformer; for very high input SNR values, these two beamformers have a very similar performance.
Next, another set of simulations is performed in a similar setting, which is shown in Figure 7. The only difference is that there is 5° DOA and polarization mismatch for the SOI: the actual parameters of the two desired sub-signals are (95°, 6.5°, 95°, 50°) and (95°, 6.5°, 5°, 5°). For lower input SNR values (from −20 dB to 0 dB), the three beamformers have a similar output SINR, all worse than the example without steering vector mismatch. As the input SNR increases, difference among the three beamformers appears. From SNR = 0 dB onwards, the output SINR of FQWCCB and the full Q-Capon beamformer continue to increase, while the output SINR of the Q-Capon beamformer starts to drop slowly. Compared with the full Q-Capon beamformer, the FQWCCB has roughly achieved a 3 dB higher output SINR than the full Q-Capon beamformer, which means that the FQWCCB has a better performance than the full Q-Capon beamformer.

5.3. Performance with DOA and Polarization Mismatch

Next, we investigate their performance in the presence of both DOA and polarization errors with SNR = SIR = 15 dB. The output SINR with respect to the number of snapshots is shown in Figure 8 in the presence of 1° error for the SOI, where the actual DOA and polarization parameters are (91°, 2.5°, 91°, 46°) and (91°, 2.5°, 1°, 1°).
We can see that the FQWCCB has the best performance in the presence of DOA and polarization mismatch, with output SINR above 20 dB for a varied number of snapshots. Both the full Q-Capon beamformer and the Q-Capon beamformer experience a loss of SINR (below 5 dB at the beginning and increasing with the snapshots number). However, the full Q-Capon beamformer has still achieved a much higher output SINR than the Q-Capon beamformer, and this gap increases with the increase of snapshots number. Figure 9 shows a similar trend in the presence of a 5° error. For the Q-Capon beamformer, the output SINR stays below 5 dB with little fluctuation, while the full Q-Capon beamformer and its worst-case constraint version can achieve a relatively high output SINR level of above 15 dB with 500 snapshots. Overall, we can see that the proposed full Q-Capon beamformer (both with worst-case constraint and without worst-case constraint) is more robust against array steering vector errors than the original Q-Capon beamformer.

6. Conclusions

In this paper, a full quaternion model has been developed for adaptive beamforming based on crossed-dipole vector sensor arrays, with a new full quaternion Capon beamformer derived. Different from previous studies in quaternion-valued adaptive beamforming, we have considered a quaternion-valued desired signal, given the recent research development in quaternion-valued communications. The proposed beamformer has a better performance and a much lower computational complexity than a previously proposed Q-Capon beamformer, and is also shown to be more robust against array pointing errors. Furthermore, based on this full quaternion model, the robust beamforming problem has also been studied in the presence of steering vector errors, and a worst-case-based robust beamformer was developed. The effectiveness of the full Q-Capon beamformer and the robustness of the further developed worst-case-based beamformer was verified by computer simulations.

Appendix A

The gradient of a quaternion vector u = w H C λ H with respect to w * can be calculated as follows:
w * u = [ w 1 * u w 2 * u . . . w n * u ] T
where w n , n = 1 , 2 , , N is the n-th quaternion-valued coefficient of the beamformer. Then,
w 1 * u = 1 4 ( w 1 a u + w 1 b u i + w 1 c u j + w 1 d u k )
where
w 1 * = w 1 a w 1 b i w 1 c j w 1 d k
Since w 1 a is real-valued, with the chain rule [24], we have
w 1 a u = w 1 a ( w H ) C λ H + w H w 1 a ( C λ H ) = [ 1 0 0 . . . 0 ] C λ H
Similarly,
w 1 b u = [ i 0 0 . . . 0 ] C λ H w 1 c u = [ j 0 0 . . . 0 ] C λ H w 1 d u = [ k 0 0 . . . 0 ] C λ H
Hence,
w 1 * u = 1 4 ( 4 R e { C λ H } 1 ) = R e { C λ H } 1
where the subscript { } 1 in the last item means taking the first entry of the vector.
Finally,
w * u = R e { C λ H }
The gradient of the quaternion vector v = λ C H w H with respect to w * can be calculated in the same way:
w 1 a v = λ C H w 1 a ( w H ) + w 1 a ( λ C H ) w H = λ C H [ 1 0 0 . . . 0 ] T
Similarly,
w 1 b v = λ C H [ i 0 0 . . . 0 ] T w 1 c v = λ C H [ j 0 0 . . . 0 ] T w 1 d v = λ C H [ k 0 0 . . . 0 ] T
Thus, the gradient can be expressed as
w 1 * v = 1 2 { C λ H } 1 *
Finally,
w * v = 1 2 ( C λ H ) *
The gradient of c w = w H R w can be calculated as follows.
w * c w = [ w 1 * c w w 2 * c w . . . w n * c w ] T
w 1 * c w = 1 4 ( w 1 a c w + w 1 b c w i + w 1 c c w j + w 1 d c w k )
Now we calculate the gradient of c w with respect to the four components of w 1 .
w 1 a c w = w 1 a ( w H R ) w + w H R w 1 a w = [ 1 0 0 . . . 0 ] R w + w H R [ 1 0 0 . . . 0 ] T
The other three components are,
w 1 b c w = [ i 0 0 . . . 0 ] R w + w H R [ i 0 0 . . . 0 ] T w 1 c c w = [ j 0 0 . . . 0 ] R w + w H R [ j 0 0 . . . 0 ] T w 1 d c w = [ k 0 0 . . . 0 ] R w + w H R [ k 0 0 . . . 0 ] T
Hence,
w 1 * c w = R e { R w } 1 1 2 { R w } 1 * = 1 2 { R w } 1
Finally,
w * c w = 1 2 R w
Combining (A7), (A11) and (A16), with (22), we have
w * l ( w , λ ) = 1 2 ( R w + C λ H ) = 0
Further,
w = 2 R 1 C λ H
Subsituting (A18) into (18),
λ = 1 2 f ( C H R 1 C ) 1
Finally,
w = R 1 C ( C H R 1 C ) 1 f H

References

  1. Compton, R., Jr. On the performance of a polarization sensitive adaptive array. IEEE Trans. Antennas Propag. 1981, 29, 718–725. [Google Scholar] [CrossRef]
  2. Nehorai, A.; Ho, K.C.; Tan, B.T.G. Minimum-noise-variance beamformer with an electromagnetic vector sensor. IEEE Trans. Signal Process. 1999, 47, 601–618. [Google Scholar] [CrossRef]
  3. Godara, L.C. Application of antenna arrays to mobile communications. II. Beam-forming and direction-of-arrival considerations. Proc. IEEE 1997, 85, 1195–1245. [Google Scholar] [CrossRef]
  4. Xu, Y.G.; Liu, T.; Liu, Z.W. Output SINR of MV beamformer with one EM vector sensor of and magnetic noise power. In Proceedings of the International Conference on Signal Processing, Montreal, QC, Canada, 17–21 May 2004; pp. 419–422. [Google Scholar]
  5. Le Bihan, N.; Mars, J. Singular value decomposition of quaternion matrices: A new tool for vector-sensor signal processing. Signal Process. 2004, 84, 1177–1199. [Google Scholar] [CrossRef]
  6. Miron, S.; Le Bihan, N.; Mars, J.I. High resolution vector-sensor array processing using quaternions. In Proceedings of the IEEE Workshop on Statistical Signal Processing, Bordeaux, France, 17–20 July 2005; pp. 918–923. [Google Scholar]
  7. 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]
  8. Gong, X.; Xu, Y.; Liu, Z. Quaternion ESPRIT for direction finding with a polarization sentive array. In Proceedings of the International Conference on Signal Processing, Beijing, China, 26–29 October 2008; pp. 378–381. [Google Scholar]
  9. Tao, J.W.; Chang, W.X. A novel combined beamformer based on hypercomplex processes. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 1276–1289. [Google Scholar] [CrossRef]
  10. Tao, J.W. Performance analysis for interference and noise canceller based on hypercomplex and spatio-temporal-polarisation processes. IET Radar Sonar Navig. 2013, 7, 277–286. [Google Scholar] [CrossRef]
  11. Zhang, X.R.; Liu, W.; Xu, Y.G.; Liu, Z.W. Quaternion-valued Robust Adaptive Beamformer for Electromagnetic Vector-sensor Arrays with Worst-case Constraint. Signal Process. 2014, 104, 274–283. [Google Scholar]
  12. Jiang, M.D.; Liu, W.; Li, Y. Adaptive Beamforming for Vector-Sensor Arrays Based on Reweighted Zero-Attracting Quaternion-Valued LMS Algorithm. IEEE Trans. Circuits Syst. II Express Br. 2016, 63, 274–278. [Google Scholar] [CrossRef]
  13. Gou, X.; Xu, Y.; Liu, Z.; Gong, X. Quaternion-capon beamformer using crossed-dipole arrays. In Proceedings of the IEEE International Symposium on Microwave, Antenna, Propagation, and EMC Technologies for Wireless Communications (MAPE), Beijing, China, 1–3 November 2011; pp. 34–37. [Google Scholar]
  14. Isaeva, O.M.; Sarytchev, V.A. Quaternion presentations polarization state. In Proceedings of the 2nd IEEE Topical Symposium of Combined Optical-Microwave Earth and Atmosphere Sensing, Atlanta, GA, USA, 3–6 April 1995; pp. 195–196. [Google Scholar]
  15. Wysocki, B.J.; Wysocki, T.A. On an Orthogonal Space-Time-Polarization Block Code. J. Commun. 2009, 4, 20–25. [Google Scholar]
  16. Liu, W. Channel equalization and beamforming for quaternion-valued wireless communication systems. J. Frankl. Inst. 2016. [Google Scholar] [CrossRef]
  17. Li, J.; Stoica, P. (Eds.) Robust Adaptive Beamforming; John Wiley & Sons: Hoboken, NJ, USA, 2005. [Google Scholar]
  18. Liu, W.; Weiss, S. Wideband Beamforming: Concepts and Techniques; John Wiley & Sons: Chichester, UK, 2010. [Google Scholar]
  19. Vorobyov, S.A.; Gershman, A.B.; Luo, Z.Q. Robust adaptive beamforming using worst-case performance optimization: A solution to the signal mismatch problem. IEEE Trans. Signal Process. 2003, 51, 313–324. [Google Scholar] [CrossRef]
  20. Yu, L.; Liu, W.; Langley, R.J. Novel robust beamformers for coherent interference suppression with DOA estimation errors. IET Microw. Antennas Propag. 2010, 4, 1310–1319. [Google Scholar] [CrossRef]
  21. Hamilton, W.R., II. On quaternions; or on a new system of imaginaries in algebra. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1844, 25, 10–13. [Google Scholar]
  22. Huang, L.; So, W. On left eigenvalues of a quaternionic matrix. Linear Algebra Appl. 2001, 323, 105–116. [Google Scholar] [CrossRef]
  23. Zhang, F. Quaternions and matrices of quaternions. Linear Algebra Appl. 1997, 251, 21–57. [Google Scholar] [CrossRef]
  24. Jiang, M.D.; Li, Y.; Liu, W. Properties of a General Quaternion-Valued Gradient Operator and Its Application to Signal Processing. Front. Inf. Technol. Electron. Eng. 2016, 17, 83–95. [Google Scholar]
  25. Capon, J. High-resolution Frequency-wavenumber Spectrum Analysis. Proc. IEEE 1969, 57, 1408–1418. [Google Scholar] [CrossRef]
  26. Frost, O.L., III. An algorithm for linearly constrained adaptive array processing. Proc. IEEE 1972, 60, 926–935. [Google Scholar] [CrossRef]
  27. Wang, M.; Ma, W. A structure-preserving algorithm for the quaternion Cholesky decomposition. Appl. Math. Comput. 2013, 223, 354–361. [Google Scholar] [CrossRef]
Figure 1. A crossed-dipole linear array with N vector sensors.
Figure 1. A crossed-dipole linear array with N vector sensors.
Electronics 06 00034 g001
Figure 2. Beam pattern obtained from the full Q-Capon beamformer with θ = 90 ° and γ = 60 °.
Figure 2. Beam pattern obtained from the full Q-Capon beamformer with θ = 90 ° and γ = 60 °.
Electronics 06 00034 g002
Figure 3. The resultant beam pattern for (1) θ = 90 ° , γ = 60 ° , η = 80 ° ; (2) θ = 90 ° , γ = 60 ° , η = 30 ° .
Figure 3. The resultant beam pattern for (1) θ = 90 ° , γ = 60 ° , η = 80 ° ; (2) θ = 90 ° , γ = 60 ° , η = 30 ° .
Electronics 06 00034 g003
Figure 4. The resultant beam pattern for (1) θ = 90 ° , γ = 60 ° , η = 70 ° ; (2) θ = 90 ° , γ = 60 ° , η = 50 ° .
Figure 4. The resultant beam pattern for (1) θ = 90 ° , γ = 60 ° , η = 70 ° ; (2) θ = 90 ° , γ = 60 ° , η = 50 ° .
Electronics 06 00034 g004
Figure 5. The resultant full quaternion worst-case constraint beamformer (FQWCCB) beam pattern with θ = 90 ° and γ = 60 ° .
Figure 5. The resultant full quaternion worst-case constraint beamformer (FQWCCB) beam pattern with θ = 90 ° and γ = 60 ° .
Electronics 06 00034 g005
Figure 6. Output signal-to-interference-plus-noise ratio (SINR) versus input signal-to-noise ratio (SNR) for snapshots number 100, without steering vector mismatch.
Figure 6. Output signal-to-interference-plus-noise ratio (SINR) versus input signal-to-noise ratio (SNR) for snapshots number 100, without steering vector mismatch.
Electronics 06 00034 g006
Figure 7. Output SINR versus input SNR for snapshots number 100, with 5° mismatch.
Figure 7. Output SINR versus input SNR for snapshots number 100, with 5° mismatch.
Electronics 06 00034 g007
Figure 8. Output SINR versus snapshot number with SNR = SIR = 15 dB and 1° error.
Figure 8. Output SINR versus snapshot number with SNR = SIR = 15 dB and 1° error.
Electronics 06 00034 g008
Figure 9. Output SINR versus snapshot number with SNR = SIR = 15 dB and 5° error.
Figure 9. Output SINR versus snapshot number with SNR = SIR = 15 dB and 5° error.
Electronics 06 00034 g009

Share and Cite

MDPI and ACS Style

Lan, X.; Liu, W. Fully Quaternion-Valued Adaptive Beamforming Based on Crossed-Dipole Arrays. Electronics 2017, 6, 34. https://0-doi-org.brum.beds.ac.uk/10.3390/electronics6020034

AMA Style

Lan X, Liu W. Fully Quaternion-Valued Adaptive Beamforming Based on Crossed-Dipole Arrays. Electronics. 2017; 6(2):34. https://0-doi-org.brum.beds.ac.uk/10.3390/electronics6020034

Chicago/Turabian Style

Lan, Xiang, and Wei Liu. 2017. "Fully Quaternion-Valued Adaptive Beamforming Based on Crossed-Dipole Arrays" Electronics 6, no. 2: 34. https://0-doi-org.brum.beds.ac.uk/10.3390/electronics6020034

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