Adaptive Beamforming Based on Complex Quaternion Processes

Motivated by the benefits of array signal processing in quaternion domain, we investigate the problem of adaptive beamforming based on complex quaternion processes in this paper. First, a complex quaternion least-mean squares (CQLMS) algorithm is proposed and its performance is analyzed.The CQLMS algorithm is suitable for adaptive beamforming of vector-sensor array. The weight vector update of CQLMS algorithm is derived based on the complex gradient, leading to lower computational complexity. Because the complex quaternion can exhibit the orthogonal structure of an electromagnetic vector-sensor in a natural way, a complex quaternion model in time domain is provided for a 3-component vector-sensor array. And the normalized adaptive beamformer using CQLMS is presented. Finally, simulation results are given to validate the performance of the proposed adaptive beamformer.

An electromagnetic (EM) vector-sensor belongs to a polarization-sensitive sensor.It can measure both horizontal and vertical polarization components of the received waveform.By employing the polarization diversity of polarized signals, the performance of beamformers can be enhanced [16,17].Thus, spatiotemporal-polarization processing based on EM vector-sensor arrays had come to some researcher's notice.In [18], a minimum variance distortionless response (MVDR) beamformer with an EM vector-sensor was presented and its performance was analyzed.In [19], a blind linearly constrained minimum variance (LCMV) beamformer was investigated for multiple wideband fast frequency-hop signals.But, these researches are based on the so-called "long vector" model.In conventional "long vector" based approaches, six components of an EM vector-sensor are treated as six independent and peer-to-peer elements, so that the output data of each component is arranged in a long vector.The main advantage of "long vector" based approaches is that it can use the well-known matrix algebra over the real or the complex field.The "long vector" based approaches, however, have the drawback of destroying locally the orthogonal structure of an EM vector-sensor [2,3].Instead of "long vector" based approach, the complex quaternion domain facilitates modelling and processing of four-dimensional complex signals.The complex quaternion model can exhibit the orthogonal structure of an EM vector-sensor in a natural way.Motivated by the benefits of array signal processing in quaternion domain, the quaternion beamformers [20][21][22][23][24][25][26] were recently 2 Mathematical Problems in Engineering developed.In [20], a quaternion minimum mean square error (MMSE) algorithm was proposed and applied to the beamforming of a two-component vector-sensor array.In [22], a quaternion-capon beamformer using a crossed dipole array was proposed to improve the robustness of Capon beamformer.In [23,24], a quaternion MVDR beamformer with an interference and noise canceller (INC) was proposed to increase the output signal to interference-plus-noise ratio, and its performance was analyzed.In [25], a quaternion spatial matched filter (QSMF) is presented and its characteristics are analyzed.Based on quaternion widely linear processing, a generalized sidelobe canceller (GSC) was proposed in [26].But, these beamformers are only suitable for 2-component vector-sensor array.To make use of 3-component vectorsensor array, the complex quaternions may be introduced in adaptive beamforming.The adaptive beamformers can adapt to the incoming data.They can be divided into three general categories [27]: (1) beamformers with the sample matrix inversion (SMI) or the direct matrix inversion (DMI) technique, implementation of which requires the inversion of the sample covariance matrix.This category includes the LCMV-SMI and MVDR-SMI beamformers, (2) beamformers with the recursive least squares (RLS) algorithm, which implement the inversion recursively.This category includes the MVDR-RLS and LCMV-GSC-RLS beamformers, and (3) beamformers with the least-mean squares (LMS) algorithm, which adapt classical steepest descent algorithms to the optimization problem in order to find the optimal weight vector.This category includes the MMSE-LMS and GSC-LMS beamformers.Unfortunately, adaptive beamforming based on complex quaternion processes has received a little attention.
In this paper, we propose a complex quaternion leastmean squares (CQLMS) algorithm and investigate its performances.The CQLMS algorithm is suitable for adaptive beamforming of vector-sensor array.Unlike the quaternion least-mean squares (QLMS) algorithm [8], the weight vector update of the CQLMS algorithm is derived based on the complex gradient.Thus, the computational complexity of the CQLMS algorithm is less than that of the QLMS algorithm.Further, a complex quaternion model in time domain is provided for a 3-component vector-sensor array.And the normalized adaptive beamformer using the proposed CQLMS algorithm is presented.From the incoming data, the weight vector of this adaptive beamformer can be adaptively updated using the proposed CQLMS algorithm.

Complex Quaternion Algebra
The complex quaternion algebra (or biquaternion algebra) is well known as a four-dimensional vector space over the complex number field C. In the complex quaternion number field Q, 1, , , and  denote four basis elements and they satisfy the following multiplication laws: and 1 acts as a unity element.

Complex Quaternions. A complex quaternion variable
∈ Q can be written in its Cartesian form as  =  0 +  1  +  2  +  3 , where  0 ,  1 ,  2 ,  3 ∈ C. According to this definition, real numbers, complex numbers, and real quaternions can be regarded as the special cases of complex quaternions.Some known terminologies on complex quaternions are listed below.For  =  0 + 1 + 2 + 3  ∈ Q, it can be written as  = S()+V(), where S() =  0 is the scalar component of  and V() =  1  +  2  +  3  is the vector component of .Also,  can be written as  = R() + I(), where are, respectively, the real and imaginary components of  and R() and I() are real quaternions.The relations between complex imaginary unit  = √ −1 and quaternion imaginary units are The complex conjugate of  is defined as 3 , where  *  ( = 0, 1, 2, 3) are the complex conjugate of complex coefficients of .The quaternion conjugate of  is defined as  (4) Note that the complex quaternions are not a normed algebra under this norm, so in general || ̸ = |||| for ,  ∈ Q.Some known basic properties on complex quaternions are listed below: let ,  ∈ Q be given.Then (1)

Vectors and Matrices of Complex
where ‡ denotes Hermitian conjugate-transposition.With this definition, the norm of a complex quaternion-valued vector a ∈ Q  is given by

Complex Quaternion-Valued Matrices. For any complex quaternion-valued matrix
Several basic properties on complex quaternion-valued matrices are listed below: given two matrices

Complex Representation of a Complex Quaternion-Valued Matrices. For any complex quaternion-valued matrix
Various operation properties on complex representation of complex quaternion-valued matrices are listed below: given three complex quaternion-valued matrices A, B ∈ Q × , D ∈ Q × , and a complex number  ∈ C, then , where The above complex quaternions and their properties are found in [3,28,29].

The Complex Quaternion LMS Algorithm
3.1.The Derivation of CQLMS Algorithm.If complex quaternion-valued series x() and complex quaternion-valued series () denote, respectively, the output signal vector and the desired signal, we assume that x() is a random vector and employ a linear model to estimate ().It is given by d() = w ‡ ()x(), where w() is the complex quaternion-valued weight vector.The estimated error is defined as The real-valued quadratic cost function is defined as Based on the cost function, within the steepest descent optimization, the following complex quaternion gradients need to be calculated: where w() = w 0 () + w 1 () + w 2 () + w 3 ().Based on the operation properties of complex gradients in [30], the gradients can be expressed as (the full derivation can be found in Appendix) Subsequently, the update of the adaptive weight vector can be written as where  denotes the iterative stepsize.Equation ( 12) refers to the update of the adaptive weight vector.When x() and () belong to real quaternion domain, (12) can be rewritten as It is noted that the update (13) is not identical to the update of QLMS algorithm proposed in [8].The update of QLMS is given by Comparing with (13), the computation of ( 14) increases  quaternion multiplications and 2 quaternion additions, where  is the column dimension of x().A quaternion multiplication consists of 16 real multiplications and 12 real additions, and a quaternion addition consists of 4 real additions.At an update, thus, the CQLMS algorithm decreases 16 real multiplications and 20 real additions than the QLMS algorithm.Further, when the complex quaternion-valued series x() and () degenerate to the complex-valued series, (12) can be rewritten as Equation ( 15) is identical to the update of the adaptive weight vector for complex LMS (CLMS) algorithm [31].Thus, the CLMS algorithm can be regarded as the special case of CQLMS algorithm.

The Property of CQLMS Algorithm.
To analyze the convergence of CQLMS algorithm in the mean, we assume that () = w ‡ opt x(), where w opt is the optimal weight vector.Defining the weight error vector as V() = w() − w opt , then the error () is given by Inserting ( 16) into ( 12) and subtracting w opt from both sides of ( 12), we have where R  () = x()x ‡ ().Based on the properties on complex representation of complex quaternion-valued matrices, the complex representation of ( 17) can be written as Since R  () is a complex quaternion Hermitian matrix, Ψ(R  ()) and {Ψ(R  ())} are complex Hermitian matrix, where {⋅} denotes the mathematical expectation operation.Then the maximum eigenvalue  max of complex matrix {Ψ(R  ())} is also real, since the eigenvalues of complex Hermitian matrix are real.According to the convergent analysis of CLMS in [31], we have the following consequence: if the stepsize  satisfies 0 <  < 2/ max , the complex matrix Ψ(V()) converges to zero along with an increase of iterative number .Because V() = (1/4)E 2 Ψ(V())E ‡ 2 , where  denotes the dimension of the column vector V(), the convergence of Ψ(V()) results in the convergence of V().We come to the conclusion that if 0 <  < 2/ max , V() converges to zero along with an increase of iterative number .The convergence of weight error vector V() indicates that the weight vector w() tends to w opt along with an increase of iterative number .

The Beamforming of Vector-Sensor Array Using CQLMS Algorithm
In this section, we apply CQLMS algorithm to adaptive beamforming of vector-sensor array.First, the complex (or real) quaternion model of vector-sensor is introduced.This quaternion model is different from the quaternion model proposed in [2,3].The quaternion model proposed in this paper is a time domain model of output data, whereas the quaternion model proposed in [2,3] is a frequency-domain model of output data.Second, the normalized adaptive beamformer using CQLMS is presented.

Complex Quaternion Model of 3-Component Vector-
Sensor.Consider a scenario with a vector-sensor mounted on an aircraft, assuming that the Cartesian coordinates of vector-sensor are consistent with the airframe coordinates; that is, the airframe locates on - plane and the -axis points into the airframe.The airframe may form a reflection plane when the wavelength of an incidence EM-wave is less than the size of airframe.The studies in [32,33] had indicated that the mutual coupling between two collocated orthogonally oriented magnetic loops, which are along and -axis, and airframes is minimum.Thus, we construct a 3-component vector-sensor by using two collocated orthogonally oriented magnetic loops along and -axis and one dipole along ; see Figure 1(a).In the Cartesian coordinates, the 3-component vector-sensor has the following responses (ignoring a common constant that is determined by the antenna parameters and the distance to the antenna) [34]: where −/2 ≤  ≤ /2 is the elevation angle of incidence EM-wave and − ≤  <  is the azimuth angle of incidence EM-wave.Consider a narrowband, transverse electromagnetic plane-wave traveling in an isotropic, homogeneous medium; we assume that its electric field can be represented by where  ℎ () and  V () are the time-variant amplitudes of two orthogonal components in the electric field and  ℎ () and  V () are the time-variant phases of two orthogonal components.We assume that  ℎ (),  V (),  ℎ (), and  V () are a stationary stochastic process;  V ()/ ℎ () = tan  = sin / cos  and  V () −  ℎ () = , where  and  are constants.In this scenario, the electromagnetic wave is a completely polarized wave and the waveform polarization can be described by the so-called polarization ellipse.Thus 0 ≤  < /2 represents the auxiliary polarization angle, and − ≤  <  signifies the polarization phase difference.() is the scalar complex envelope of the waveform.When a narrowband plane-wave impinges on 3-component vector-sensor from direction (, ), three temporal signals   (),   (), and   () are recorded on magnetic loop antennae in and -axis and dipole antenna in , respectively.They are given by where () is a complex stationary stochastic process.

Real Quaternion Model of 2-Component Vector-Sensor.
Using only two collocated orthogonally oriented magnetic loops along and -axis, we can also construct a 2-component vector-sensor; see Figure 1(b).When a narrowband plane-wave impinges on 2-component vector-sensor, we have where is the responses of 2-component vector-sensor.Since a quaternion may be expressed by two complex numbers, the output model of 2-component vector-sensor can be compactly written in terms of quaternion: where (, , , ) is a quaternion; that is, In this way, we define a transformation to map a two-dimensional complex-valued signal space to a quaternion-valued signal space.

The Beamforming of Vector-Sensor Array.
Consider that a three-dimensional array consists of  arbitrarily spaced but identically oriented 3-component (or 2-component) vectorsensors.In this three-dimensional array, the th vector-sensor has the location r  = [  ,   ,   ]  (1 ≤  ≤ ).Taking the first vector-sensor as reference, the inter-vector-sensor spatial phase-factor vector of array is where   = (2/)u  r  .u = [sin  cos , sin  sin , cos ]  denotes direction-cosine vector and  is the wavelength of incidence EM-waves.Assume a scenario of  incoherent and complete polarization waves with 2D direction angles (  ,   ) and polarization parameters (  ,   ) impinging on this three-dimensional array; the output of array is given by a column vector x() ∈ Q  : where a  = (  ,   ,   ,   )q(  ,   ) denotes the steering vector.
The series x() obtained by sampling x() is multiplied by a weight vector to produce the array output series; it is given by When all of vector-sensors in array are the 3-component vector-sensor, a complex quaternion-valued output model ( 24) can be employed.Thus, the output series x() and weight vector series w() are the complex quaternion-valued vectors.Using the CQLMS algorithm, the update of the adaptive weight vector can be given by (12).In order to speed the convergence and reduce the MSE of estimates, the normalized CQLMS algorithm is adopted.Thus, the update of the adaptive weight vector becomes w ( + 1) = w () +  () x ()  † () , where () is the time-variant stepsize.In most cases, a sample-dependent stepsize is adopted [27]; that is, with  > 0 and 0 <  < 1.Note that  is a weight value.It can avoid the fact that () becomes too large when x ‡ ()x() is too small.Normally,  will be close to one.The constant  can be selected in advance.If  is too small, the stepsize leads to slow convergence.If  is too large, the algorithm has a large excess error or has stability problems.When all of vector-sensors in array are the 2-component vector-sensor, a quaternion-valued output model ( 29) can be used.Thus, the output series x() and weight vector series w() are the real quaternion-valued vectors.Thus, the update of the adaptive weight vector is given by (13).

Monte Carlo Simulations
In this section, we evaluate the performance of the proposed algorithm by two experiments.Consider a linear array consisting of two vector-sensors; see Figure 1(c).Let  = /2.Assume that two narrowband, completely polarized planewaves impinge on this array.One is the desired signal characterized by arrival angles (  ,   ) and polarization parameters (  ,   ), whereas the other is the interference characterized by (  ,   ) and (  ,   ).The desired signal is uncorrelated with the interference, and the interference is also uncorrelated with added noise.We define, respectively, the signal to noise ratio SNR = 10 log(  /  ), the interference to noise ratio INR = 10 log(  /  ), and the signal to interference-plus-noise ratio SINR = 10 log(  /(  +   )), where   ,   , and   are the power of desired signal, the power of interference, and the power of added noise, respectively.All simulation results are based on 50 times Monte Carlo runs.
Experiment 1 (beamforming using a 3-component vector-sensor array).In this experiment, two 3-component vectorsensors are adopted in this linear array.Assume that   = 45 ∘ ,   = 45 ∘ ,   = 60 ∘ ,   = 60 ∘ ,   = 60 ∘ ,   = 45 ∘ ,   = 30 ∘ ,   = 60 ∘ , SNR = 10 dB, and INR = 20 dB.Employing the CQLMS algorithm, Figure 2 shows the norm of weight vector ‖w()‖ and the cost function () (i.e., MSE) as a function of the iteration number , with  = 0.9 and various .To compare, Figure 3 shows the simulation results of the complex LMS (CLMS) algorithm based on the "long vector" model in [16].From Figure 2, it is shown that the CQLMS algorithm is convergent.When  increases, the convergence is speeded up.For example, 20 iteration steps are only requested to reach stable value at  = 0.8.After reaching stability, the average MSE of CQLMS algorithm is 0.15.The convergent speed of CLMS algorithm is identical to the CQLMS algorithm, but its average MSE is 0.6.Thus, the excess error of the CLMS algorithm is larger than that of the CQLMS algorithm.Further, we compare the performance of the proposed beamformer with the MVDR-SMI beamformer in [18].Figure 4 shows the output SINR of two beamformers versus   , where the adaptive beamformer with CQLMS algorithm ( = 0.4 and  = 0.9) is shown in Figure 4(a) and the MVDR-SMI beamformer is shown in Figure 4(b).From these curves, it is seen that the output SINR increases along with an increase of the iterative number (or snapshots)  and Mathematical Problems in Engineering   the worst output SINR is obtained if   is close to   .This implies that when the interference is nullified, the desired signal is also nullified.When   deviates from   , the output SINR increases.The maximum output SINR approaches up to 10 dB.Compared with the proposed beamformer, the output SINR of the MVDR-SMI beamformer is larger when   is close to   .This implies that the interference is not nullified effectively.In other words, the performance of the beamformer with the CQLMS algorithm is superior to that of the MVDR-SMI beamformer.

Figure 4 :
Figure 4: The output SINR of two beamformers versus   .
The norm of weight vector versus iteration number  The cost function () versus iteration number