# DOA Estimation of Cylindrical Conformal Array Based on Geometric Algebra

Authors: Minjie Wu; Xiaofa Zhang; Jingjian Huang; Naichang Yuan
Journal: International Journal of Antennas and Propagation
 (2016)

---

## Abstract

Due to the variable curvature of the conformal carrier, the pattern of each element has a different direction. The traditional method of analyzing the conformal array is to use the Euler rotation angle and its matrix representation. However, it is computationally demanding especially for irregular array structures. In this paper, we present a novel algorithm by combining the geometric algebra with Multiple Signal Classification (MUSIC), termed as GA-MUSIC, to solve the direction of arrival (DOA) for cylindrical conformal array. And on this basis, we derive the pattern and array manifold. Compared with the existing algorithms, our proposed one avoids the cumbersome matrix transformations and largely decreases the computational complexity. The simulation results verify the effectiveness of the proposed method.


---

## Body

## 1. Introduction

A conformal antenna is an antenna that conforms to a prescribed shape. The shape can be some part of an airplane, high-speed missile, or other vehicle [1]. Their benefits include reduction of aerodynamic drag, wide angle coverage, and space-saving [2]. Nevertheless, due to the complex curved surface structure, the pattern of each antenna is inconsistent. Thus, the conformal array can no longer be regarded as a simple isotropic one. The pattern multiplication theorem is not available as well. Most classical DOA estimation algorithms cannot be directly transplanted to such scene.In recent years, there has been a considerable interest in estimating DOAs for conformal array. Milligan used Euler rotation angles to find the patterns with elements in a conformal array that requires one to rotate not only the direction but also the polarization [3]. In [2], Wang et al. proposed a uniform method for the element polarized pattern transformation of arbitrary 3D conformal arrays based on Euler rotation. Yang et al. introduced a conformal array DOA algorithm with an unknown source number; the method was realized by virtue of the pseudo expected signal [4]. However, the root mean square error (RMSE) deteriorated severely when the number of snapshots was small. Up to the present, we have observed that most of DOA estimation algorithms for conformal array are based on the Euler rotation transformation which converts the local coordinate system to the global coordinate system. Though the Euler rotation angle is a useful tool for spatial rotation transformation [5], a huge amount of computation is incurred.Geometric algebra is the largest possible associative algebra that integrates all algebraic systems (algebra of complex numbers, matrix algebra, quaternion algebra, etc.) into a coherent mathematical language [6]. Three-dimensional pattern analysis of arbitrary conformal arrays using the mathematical framework of the geometric algebra was introduced [7]. Nevertheless, this mathematical language was not transplanted to the DOA estimation. In [8], Zou et al. took several elements as a new one and transformed the original array into another regular array to estimate the DOA. However, this method was only suitable for some particular array structures. Combining the MUSIC with geometric algebra to solve the DOA estimation has not been addressed in the literature. In this paper, we fill this gap and study the problem based on the cylindrical conformal array. Compared with the existing methods, the proposed one has three main advantages. Firstly, it does not need to calculate the rotation matrices and therefore has a much lower computational complexity. Subsequently, it is not limited to the cylindrical conformal array due to its strong commonality. Finally, it can still work effectively even when the number of polarized signals is larger than that of the array elements.The structure of this paper is as follows. In Section2, the rotors in geometric algebra which establishes the mathematical knowledge of transformation is briefly introduced. In Section 3, we derive the cylindrical conformal array manifold using rotors and present the GA-MUSIC algorithm. In addition, to better explain the superiority of GA-MUSIC in reducing the computational complexity, we briefly introduce the Euler angle and compare it with the proposed algorithm. Simulations using the proposed method for cylindrical conformal array are given in Section 4. Finally, the conclusions are drawn.


## 2. Rotors in Geometric Algebra

Geometric algebra was first introduced by the British mathematician, named Clifford, in the nineteenth century. He constructed the geometric product by combining the inner product with the outer product. The main advantage of the geometric algebra is embodied in processing the rotation transformation [9]. Various rotations can be described by an element called the rotor. A rotor is more general than an Euler rotation angle because a rotor can be used in an arbitrary dimensional space.We begin by introducing a new product between vectors that we call the outer product. Let us use the wedge symbol “∧” to denote outer product with the properties listed in Table 1.Table 1
Properties of the outer product.
Property
Meaning
Antisymmetry
(
x
∧
y
)
=
-
(
y
∧
x
)
Scaling
x
∧
(
γ
y
)
=
γ
(
x
∧
y
)
Distributivity
x
∧
(
y
+
z
)
=
(
x
∧
y
)
+
(
x
∧
z
)
Associativity
x
∧
(
y
∧
z
)
=
(
x
∧
y
)
∧
zThe outer product is regarded as the “addition operator” of subspaces, in that the outer productx
∧
y spans the subspace that x and y span together, as long as x and y are independent. The geometry is illustrated in Figure 1.Figure 1
The geometry of outer product (a)x
∧
y and (b) y
∧
x.
(a)
(b)Next, we will introduce the fundamental product of the geometric algebra, namely, geometric product. It is simply the sum of the outer and inner product:(1)
x
y
=
x
•
y
+
x
∧
y
.Reversing the order ofx and y in (1), by means of the symmetry of the inner product and the antisymmetry of the outer product, it follows that(2)
y
x
=
x
•
y
-
x
∧
y
.Thus, combining (1) with (2), the inner product and the outer product can be uniformly represented by the geometric product:(3)
x
•
y
=
x
y
+
y
x
2
x
∧
y
=
x
y
-
y
x
2
.Generally, we call an outer product ofk vectors a k-blade. The value of k is called the grade of the blade. Specifically, the top-grade blades E
n in an n-dimensional space are called pseudoscalars. In principle, blades are just elements of the geometric algebra. To be useful for doing geometry, blades can be interpreted as subspaces. x
∧
y is a 2-blade as shown in Figure 1(a). By introducing vector z, we can construct the 3-blade x
∧
y
∧
z, which coincides with the Cartesian coordinate system. As shown in Figure 2, the unit vectors e
x, e
y, and e
z, in that order, form a right-handed Cartesian coordinate system. E
3 is the pseudoscalar, relative to the origin denoted by O. The 3-blade is drawn as a parallelepiped. The volume depicts the weight of the 3-blade, but, in principle, blades have no specific shape.Figure 2
The geometry of 3-blade.As shown in Figure3, vector y is obtained by rotating vector x with θ. We can regard the rotation as two consecutive reflections, first in a and then in b. The expression that reflects x in the line with direction a is (4)
x
′
=
a
x
a
=
2
a
•
x
a
-
x
.Figure 3
Rotation of vectorx.In (4), a must be unit length. We give the detailed derivation process in Appendix. The regularity condition can be dropped by using an inverse geometric product, as in(5)
x
′
=
a
x
a
-
1
=
2
a
•
x
a
-
1
-
x
,where(6)
a
-
1
=
a
a
•
awith ·
-
1 representing the inverse operator.Then,y can be obtained by reflecting x
′ in the line with direction b:(7)
y
=
b
x
′
b
-
1
=
b
a
x
a
-
1
b
-
1
=
b
a
x
b
a
-
1
=
R
x
R
-
1
.Thus, we can identifyR as the rotor. To proceed further, we rewrite R according to the definition of the geometric product:(8)
R
=
b
a
=
b
•
a
+
b
∧
a
.Here, we consider the case that the vectors are unit length. This assumption is reasonable, because the basic vectors of the Cartesian coordinate system satisfy it as well. The geometric product ofb
∧
a itself is(9)
b
∧
a
b
∧
a
=
b
a
-
b
•
a
b
•
a
-
a
b
=
b
•
a
b
a
+
a
b
-
b
•
a
2
-
b
a
a
b
=
b
•
a
2
b
•
a
-
b
•
a
2
-
b
a
a
b
=
b
•
a
2
-
b
b
=
cos
2
⁡
θ
-
1
=
-
sin
2
⁡
θ
.Thus, we define the 2-bladeE
2:(10)
E
2
=
b
∧
a
sin
⁡
θ
.R can be further simplified by substituting (10) into (8):(11)
R
=
cos
⁡
θ
-
E
2
sin
⁡
θ
.The expression is similar to the polar decomposition of a complex number with the unit imaginary replaced by the 2-bladeE
2. It can also be written as the exponentials of E
2:(12)
R
=
e
-
E
2
θ
.This formalism is more useful for the log-space of rotors is linear. Up to this point, the angle of rotation in theb
∧
a-plane remains to be determined. We split x into part (x
p) parallel to b
∧
a-plane and part (x
o) orthogonal to b
∧
a-plane. Then, x
o is not affected by application R. And we infer that the rotation must be in b
∧
a-plane. As stated above, the rotation consists of two successive reflections which are orthogonal (angle-preserving) transformations. Thus, it allows us to pick any vector in b
∧
a-plane to determine the angle. Without loss of generality, we choose vector a and construct the “sandwich product” R
a
R
-
1 as shown in (7):(13)
R
a
R
-
1
=
b
a
a
a
-
1
b
-
1
=
b
a
b
-
1
,where b
a
b
-
1 is the reflection of a in b. From this, it is clear that the rotation must be over twice the angle between a and b, since the angle between a and b
a
b
-
1 is twice the angle between a and b. The negative signature in (12) represents the rotation direction.Consequently, if we want to rotate a vector counterclockwise by a specific angle, we only need to apply the rotor to the vector.


## 3. GA-MUSIC Algorithm

### 3.1. Array Manifold Modeling Based on GA-MUSIC

In this subsection, we will combine the geometric algebra with MUSIC to estimate the DOA. To illustrate the versatility of this algorithm, we considerM
×
N cylindrical conformal array as shown in Figure 4. The array contains N uniformly spaced rings on the surface. In addition, there are M dipoles distributed on each ring. We assume that each dipole is a short dipole whose output voltage is proportional to the electric field along the dipole. The angle between two consecutive elements on the same ring is β.Figure 4
The cylindrical conformal array consisting ofM
×
N short dipoles.Assume that there areP far-field narrow band polarization sources. θ
p represents pth signal elevation angle which is measured down from e
z-axis. φ
p indicates pth signal azimuth angle and is measured counterclockwise from e
x-axis. The polarization ellipse of pth signal is depicted by constants γ
p and η
p, representing the auxiliary polarization angle and the polarization phase difference [10], respectively. The array element spatial phase matrix of the pth signal can be described as follows:(14)
Υ
p
=
Υ
θ
p
,
φ
p
=
u
p
,
1
⋱
u
p
,
P
,where Υ
p is the diagonal matrix depicting the output signal spatial coherent structures. kth diagonal element, u
p
,
k
=
e
-
j
2
π
(
ε
s
T
(
θ
p
,
φ
p
)
r
k
)
/
λ
p, represents the space phase factor about pth array element. Among which, ε
s
(
θ
p
,
ϕ
p
)
=
-
[
sin
⁡
θ
p
cos
⁡
φ
p
,
sin
⁡
θ
p
sin
⁡
φ
p
,
cos
⁡
θ
p
]
T and λ
p are pth signal propagation vector and the wavelength, respectively. Symbol r
k is the element location vector and ·
T denotes the transpose operator.It is worthwhile to note that the aforementioned element spatial phase matrix,Υ
p, is derived under the global coordinate system. The azimuth and elevation are defined in Figure 4 as well. However, due to the effects of the curvature of conformal carriers, the local coordinate system is distinct from the global one. As stated above, the rotor can be used to realize the rotation between the two coordinate systems. Thus, we define the local coordinate system of (
m
,
n
)th element as shown in Figure 5.Figure 5
The local coordinate of(
m
,
n
)th element.e
x
m
n-axis is the same as e
x-axis in the global coordinate system, e
z
m
n is normal to the element surface, and e
y
m
n is tangent to the surface so as to form a right-handed coordinate system. Naturally, transforming the global coordinate into the local one is equivalent to rotating the global coordinate around e
x-axis. From the last section, we have known that (12) denotes the rotation in b
∧
a-plane with twice the value of θ. The rotation angle is(15)
ξ
=
m
-
1
β
-
M
-
1
2
β
=
m
-
M
+
1
2
β
.Substitutinge
z and e
y for b and a, respectively, then the rotor is(16)
R
m
n
=
e
-
e
z
∧
e
y
ξ
/
2
.In addition,e
x, e
y, and e
z are orthogonal to each other and the inner products between them are zero. So,(17)
e
z
∧
e
y
=
e
z
e
y
.According to the antisymmetry property of the outer product,(18)
e
z
e
y
=
e
z
∧
e
y
=
-
e
y
∧
e
z
=
-
e
y
e
z
.AdoptingE
3
=
e
x
e
y
e
z as the pseudoscalar in the three-dimensional space, (16) can be further simplified:(19)
R
m
n
=
e
E
3
e
x
ξ
/
2
.Through (7), we obtain the standard orthogonal basis in the local coordinate, and the specific calculation process can refer to [11]. Here, the results are given directly:(20)
e
x
m
n
=
R
m
n
e
x
R
m
n
-
1
=
e
x
(21)
e
y
m
n
=
R
m
n
e
y
R
m
n
-
1
=
cos
⁡
m
-
M
+
1
2
β
e
y
-
sin
⁡
m
-
M
+
1
2
β
e
z
(22)
e
z
m
n
=
R
m
n
e
z
R
m
n
-
1
=
sin
⁡
m
-
M
+
1
2
β
e
y
+
cos
⁡
m
-
M
+
1
2
β
e
z
.Referring back to (14), the remaining unknown variable is r
k. The position vector of (
m
,
n
)th element in the global frame is(23)
r
m
n
=
-
n
δ
e
x
+
R
sin
⁡
m
-
M
+
1
2
β
e
y
+
R
cos
⁡
m
-
M
+
1
2
β
e
z
m
=
1
,
…
,
M
;
n
=
1
,
…
,
N
,where δ means the spacing between adjacent rings. Thus, Υ
p can be obtained. If g is the gain when the signal perfectly matched the antenna polarization, then the element of generalized polarization sensitive matrix Q can be represented as(24)
q
=
g
cos
⁡
φ
′
sin
⁡
θ
′
sin
⁡
φ
′
sin
⁡
θ
′
cos
⁡
θ
′
T
,where θ
′ and φ
′ indicate the elevation and azimuth pointing directions of the short dipole.According to this, we can get the array manifold:(25)
a
p
=
a
θ
p
,
ϕ
p
,
γ
p
,
η
p
=
Υ
p
Q
Ψ
p
h
p
,where h
p is pth signal polarization vector [12] and can be described by γ
p and η
p, that is, h
p
=
cos
⁡
γ
p
sin
⁡
γ
p
e
j
η
p
T. Ψ
p is the steering vector of the angle field [13] and is independent of the space location:(26)
Ψ
p
=
-
sin
⁡
φ
p
cos
⁡
θ
p
cos
⁡
φ
p
cos
⁡
φ
p
cos
⁡
θ
p
sin
⁡
φ
p
0
-
sin
⁡
θ
p
.The received signals of the array are a superposition of the response of each signal, and the output can be expressed as(27)
x
t
=
∑
p
=
1
P
a
p
s
p
t
+
n
t
,where s
p
(
t
) is pth signal and n
(
t
) is assumed to be zero mean, complex Gaussian processes statistically independent of each other, with covariance σ
n
2.Up to the present, we have perfectly applied the geometric algebra to the derivation of the steering vector of the conformal array and acquired the GA-MUSIC algorithm. This is also the focus of our work. As for the content of constructing the spatial spectrum and searching the peak, the readers can refer to literature [14]. As the spatial location of the short dipole is arbitrary when introducing the “rotors,” the GA-MUSIC is not limited to any specific array geometry.


### 3.2. Comparing with the Conventional Modeling Using Euler Angle

For better understanding of the superiority of the geometric algebra in modeling the conformal array, we will briefly introduce the conventional methods of analysis based on Euler angle. In general, the transformations from the element local coordinates to the array global coordinates can be realized by three successive Euler rotations [2]. The corresponding rotation matrix can be written as (28)
R
C
,
D
,
F
=
R
x
C
R
y
D
R
z
F
=
1
0
0
0
cos
⁡
C
-
sin
⁡
C
0
sin
⁡
C
cos
⁡
C
cos
⁡
D
0
-
sin
⁡
D
0
1
0
sin
⁡
D
0
cos
⁡
D
cos
⁡
F
-
sin
⁡
F
0
sin
⁡
F
cos
⁡
F
0
0
0
1
=
cos
⁡
D
cos
⁡
F
-
cos
⁡
D
sin
⁡
F
-
sin
⁡
D
cos
⁡
C
sin
⁡
F
-
cos
⁡
F
sin
⁡
C
sin
⁡
D
cos
⁡
C
cos
⁡
F
+
sin
⁡
C
sin
⁡
D
sin
⁡
F
-
cos
⁡
D
sin
⁡
C
sin
⁡
C
sin
⁡
F
+
cos
⁡
C
cos
⁡
F
sin
⁡
D
cos
⁡
F
sin
⁡
C
-
cos
⁡
C
sin
⁡
D
sin
⁡
F
cos
⁡
C
cos
⁡
D
,where C
,
D
,
F are three successive Euler rotation angles about e
x-axis, e
y-axis, and e
z-axis, respectively. R
x
(
C
), R
y
(
D
), and R
z
(
F
) are the corresponding Euler rotation matrices. Note that, for cylindrical conformal array, two successive Euler rotations are usually sufficient [2]. The third Euler rotation matrix is added here to account for some irregular or complex conformal arrays. Moreover, as shown in (28), the rotation matrix is invertible. Thus, the inversion is taken with respect to R
(
C
,
D
,
F
), resulting in (29)
R
C
,
D
,
F
-
1
=
R
z
-
1
F
R
y
-
1
D
R
x
-
1
C
=
cos
⁡
F
sin
⁡
F
0
-
sin
⁡
F
cos
⁡
F
0
0
0
1
cos
⁡
D
0
sin
⁡
D
0
1
0
-
sin
⁡
D
0
cos
⁡
D
1
0
0
0
cos
⁡
C
sin
⁡
C
0
-
sin
⁡
C
cos
⁡
C
=
cos
⁡
D
cos
⁡
F
cos
⁡
C
sin
⁡
F
-
cos
⁡
F
sin
⁡
C
sin
⁡
D
sin
⁡
C
sin
⁡
F
+
cos
⁡
C
cos
⁡
F
sin
⁡
D
-
cos
⁡
D
sin
⁡
F
cos
⁡
C
cos
⁡
F
+
sin
⁡
C
sin
⁡
D
sin
⁡
F
cos
⁡
F
sin
⁡
C
-
cos
⁡
C
sin
⁡
D
sin
⁡
F
-
sin
⁡
D
-
cos
⁡
D
sin
⁡
C
cos
⁡
C
cos
⁡
D
.From (28) and (29), it is not hard to find that(30)
R
C
,
D
,
F
-
1
=
R
C
,
D
,
F
T
.So,R
(
C
,
D
,
F
) is the orthogonal matrix. Then, taking the transformation from the element local coordinates to the array global coordinates is equivalent to taking the transposition/inversion with respect to the aforementioned rotation matrix. If we model the conformal array adopting the Euler angle, three matrix multiplications and one matrix transposition are involved for each element.In practical applications, the matrix operations are essentially the multiplication and addition operations between elements. To quantify what we mean by this, the amounts of multiplication and addition operations of the two algorithms (i.e., GA-MUSIC and Euler angle) are computed, respectively, as shown in Table2. Suppose that one matrix transposition is regarded as one multiplication or addition operation. As is known to all, the multiplication between two 3 × 3 matrices requires 9 × 3 multiplications and 9 × 2 additions. For convenience, the multiplication operation and the addition operation are collectively referred to as the operation. Then, the calculation of (28) contains 2 × 9 × 3 + 2 × 9 × 2 operations. For the conformal array consisting of M
×
N elements, the transformation between different coordinates involves 91
×
M
N operations. In contrast to the Euler angle, GA-MUSIC avoids the cumbersome matrix transformations.Table 2
The computational complexity of GA-MUSIC and Euler angle.
Multiplications
Additions
Transpositions
Operations
Euler angle
2
×
9
×
3
×
M
N
2
×
9
×
2
×
M
N
M
N
91
×
M
N
GA-MUSIC
4
×
M
N
2
×
M
N
0
6
×
M
NFrom (20)–(22), we know that e
y
m
n and e
z
m
n are independent of e
x. Moreover, e
x
m
n can be obtained directly from (20) without extra computations. Thus, (21)-(22) can be rewritten in 2 × 2 matrix. While applying the rotor to construct the array manifold, the computational process is equivalent to 2 × 2 matrix multiplied by 2 × 1 vector. Then, the operations for each element involve 4 multiplications and 2 additions. The total amount of operations is 6
×
M
N
.Thus, the latter method largely decreases the computational complexity.To sum up, the Euler rotation and its matrix representation cannot intuitively display the whole procedure. Moreover, as the configuration of the conformal array becomes more irregular and complex, the level of complexity involved in the transformations and the number of calculations required increase significantly.


## 3.1. Array Manifold Modeling Based on GA-MUSIC

In this subsection, we will combine the geometric algebra with MUSIC to estimate the DOA. To illustrate the versatility of this algorithm, we considerM
×
N cylindrical conformal array as shown in Figure 4. The array contains N uniformly spaced rings on the surface. In addition, there are M dipoles distributed on each ring. We assume that each dipole is a short dipole whose output voltage is proportional to the electric field along the dipole. The angle between two consecutive elements on the same ring is β.Figure 4
The cylindrical conformal array consisting ofM
×
N short dipoles.Assume that there areP far-field narrow band polarization sources. θ
p represents pth signal elevation angle which is measured down from e
z-axis. φ
p indicates pth signal azimuth angle and is measured counterclockwise from e
x-axis. The polarization ellipse of pth signal is depicted by constants γ
p and η
p, representing the auxiliary polarization angle and the polarization phase difference [10], respectively. The array element spatial phase matrix of the pth signal can be described as follows:(14)
Υ
p
=
Υ
θ
p
,
φ
p
=
u
p
,
1
⋱
u
p
,
P
,where Υ
p is the diagonal matrix depicting the output signal spatial coherent structures. kth diagonal element, u
p
,
k
=
e
-
j
2
π
(
ε
s
T
(
θ
p
,
φ
p
)
r
k
)
/
λ
p, represents the space phase factor about pth array element. Among which, ε
s
(
θ
p
,
ϕ
p
)
=
-
[
sin
⁡
θ
p
cos
⁡
φ
p
,
sin
⁡
θ
p
sin
⁡
φ
p
,
cos
⁡
θ
p
]
T and λ
p are pth signal propagation vector and the wavelength, respectively. Symbol r
k is the element location vector and ·
T denotes the transpose operator.It is worthwhile to note that the aforementioned element spatial phase matrix,Υ
p, is derived under the global coordinate system. The azimuth and elevation are defined in Figure 4 as well. However, due to the effects of the curvature of conformal carriers, the local coordinate system is distinct from the global one. As stated above, the rotor can be used to realize the rotation between the two coordinate systems. Thus, we define the local coordinate system of (
m
,
n
)th element as shown in Figure 5.Figure 5
The local coordinate of(
m
,
n
)th element.e
x
m
n-axis is the same as e
x-axis in the global coordinate system, e
z
m
n is normal to the element surface, and e
y
m
n is tangent to the surface so as to form a right-handed coordinate system. Naturally, transforming the global coordinate into the local one is equivalent to rotating the global coordinate around e
x-axis. From the last section, we have known that (12) denotes the rotation in b
∧
a-plane with twice the value of θ. The rotation angle is(15)
ξ
=
m
-
1
β
-
M
-
1
2
β
=
m
-
M
+
1
2
β
.Substitutinge
z and e
y for b and a, respectively, then the rotor is(16)
R
m
n
=
e
-
e
z
∧
e
y
ξ
/
2
.In addition,e
x, e
y, and e
z are orthogonal to each other and the inner products between them are zero. So,(17)
e
z
∧
e
y
=
e
z
e
y
.According to the antisymmetry property of the outer product,(18)
e
z
e
y
=
e
z
∧
e
y
=
-
e
y
∧
e
z
=
-
e
y
e
z
.AdoptingE
3
=
e
x
e
y
e
z as the pseudoscalar in the three-dimensional space, (16) can be further simplified:(19)
R
m
n
=
e
E
3
e
x
ξ
/
2
.Through (7), we obtain the standard orthogonal basis in the local coordinate, and the specific calculation process can refer to [11]. Here, the results are given directly:(20)
e
x
m
n
=
R
m
n
e
x
R
m
n
-
1
=
e
x
(21)
e
y
m
n
=
R
m
n
e
y
R
m
n
-
1
=
cos
⁡
m
-
M
+
1
2
β
e
y
-
sin
⁡
m
-
M
+
1
2
β
e
z
(22)
e
z
m
n
=
R
m
n
e
z
R
m
n
-
1
=
sin
⁡
m
-
M
+
1
2
β
e
y
+
cos
⁡
m
-
M
+
1
2
β
e
z
.Referring back to (14), the remaining unknown variable is r
k. The position vector of (
m
,
n
)th element in the global frame is(23)
r
m
n
=
-
n
δ
e
x
+
R
sin
⁡
m
-
M
+
1
2
β
e
y
+
R
cos
⁡
m
-
M
+
1
2
β
e
z
m
=
1
,
…
,
M
;
n
=
1
,
…
,
N
,where δ means the spacing between adjacent rings. Thus, Υ
p can be obtained. If g is the gain when the signal perfectly matched the antenna polarization, then the element of generalized polarization sensitive matrix Q can be represented as(24)
q
=
g
cos
⁡
φ
′
sin
⁡
θ
′
sin
⁡
φ
′
sin
⁡
θ
′
cos
⁡
θ
′
T
,where θ
′ and φ
′ indicate the elevation and azimuth pointing directions of the short dipole.According to this, we can get the array manifold:(25)
a
p
=
a
θ
p
,
ϕ
p
,
γ
p
,
η
p
=
Υ
p
Q
Ψ
p
h
p
,where h
p is pth signal polarization vector [12] and can be described by γ
p and η
p, that is, h
p
=
cos
⁡
γ
p
sin
⁡
γ
p
e
j
η
p
T. Ψ
p is the steering vector of the angle field [13] and is independent of the space location:(26)
Ψ
p
=
-
sin
⁡
φ
p
cos
⁡
θ
p
cos
⁡
φ
p
cos
⁡
φ
p
cos
⁡
θ
p
sin
⁡
φ
p
0
-
sin
⁡
θ
p
.The received signals of the array are a superposition of the response of each signal, and the output can be expressed as(27)
x
t
=
∑
p
=
1
P
a
p
s
p
t
+
n
t
,where s
p
(
t
) is pth signal and n
(
t
) is assumed to be zero mean, complex Gaussian processes statistically independent of each other, with covariance σ
n
2.Up to the present, we have perfectly applied the geometric algebra to the derivation of the steering vector of the conformal array and acquired the GA-MUSIC algorithm. This is also the focus of our work. As for the content of constructing the spatial spectrum and searching the peak, the readers can refer to literature [14]. As the spatial location of the short dipole is arbitrary when introducing the “rotors,” the GA-MUSIC is not limited to any specific array geometry.


## 3.2. Comparing with the Conventional Modeling Using Euler Angle

For better understanding of the superiority of the geometric algebra in modeling the conformal array, we will briefly introduce the conventional methods of analysis based on Euler angle. In general, the transformations from the element local coordinates to the array global coordinates can be realized by three successive Euler rotations [2]. The corresponding rotation matrix can be written as (28)
R
C
,
D
,
F
=
R
x
C
R
y
D
R
z
F
=
1
0
0
0
cos
⁡
C
-
sin
⁡
C
0
sin
⁡
C
cos
⁡
C
cos
⁡
D
0
-
sin
⁡
D
0
1
0
sin
⁡
D
0
cos
⁡
D
cos
⁡
F
-
sin
⁡
F
0
sin
⁡
F
cos
⁡
F
0
0
0
1
=
cos
⁡
D
cos
⁡
F
-
cos
⁡
D
sin
⁡
F
-
sin
⁡
D
cos
⁡
C
sin
⁡
F
-
cos
⁡
F
sin
⁡
C
sin
⁡
D
cos
⁡
C
cos
⁡
F
+
sin
⁡
C
sin
⁡
D
sin
⁡
F
-
cos
⁡
D
sin
⁡
C
sin
⁡
C
sin
⁡
F
+
cos
⁡
C
cos
⁡
F
sin
⁡
D
cos
⁡
F
sin
⁡
C
-
cos
⁡
C
sin
⁡
D
sin
⁡
F
cos
⁡
C
cos
⁡
D
,where C
,
D
,
F are three successive Euler rotation angles about e
x-axis, e
y-axis, and e
z-axis, respectively. R
x
(
C
), R
y
(
D
), and R
z
(
F
) are the corresponding Euler rotation matrices. Note that, for cylindrical conformal array, two successive Euler rotations are usually sufficient [2]. The third Euler rotation matrix is added here to account for some irregular or complex conformal arrays. Moreover, as shown in (28), the rotation matrix is invertible. Thus, the inversion is taken with respect to R
(
C
,
D
,
F
), resulting in (29)
R
C
,
D
,
F
-
1
=
R
z
-
1
F
R
y
-
1
D
R
x
-
1
C
=
cos
⁡
F
sin
⁡
F
0
-
sin
⁡
F
cos
⁡
F
0
0
0
1
cos
⁡
D
0
sin
⁡
D
0
1
0
-
sin
⁡
D
0
cos
⁡
D
1
0
0
0
cos
⁡
C
sin
⁡
C
0
-
sin
⁡
C
cos
⁡
C
=
cos
⁡
D
cos
⁡
F
cos
⁡
C
sin
⁡
F
-
cos
⁡
F
sin
⁡
C
sin
⁡
D
sin
⁡
C
sin
⁡
F
+
cos
⁡
C
cos
⁡
F
sin
⁡
D
-
cos
⁡
D
sin
⁡
F
cos
⁡
C
cos
⁡
F
+
sin
⁡
C
sin
⁡
D
sin
⁡
F
cos
⁡
F
sin
⁡
C
-
cos
⁡
C
sin
⁡
D
sin
⁡
F
-
sin
⁡
D
-
cos
⁡
D
sin
⁡
C
cos
⁡
C
cos
⁡
D
.From (28) and (29), it is not hard to find that(30)
R
C
,
D
,
F
-
1
=
R
C
,
D
,
F
T
.So,R
(
C
,
D
,
F
) is the orthogonal matrix. Then, taking the transformation from the element local coordinates to the array global coordinates is equivalent to taking the transposition/inversion with respect to the aforementioned rotation matrix. If we model the conformal array adopting the Euler angle, three matrix multiplications and one matrix transposition are involved for each element.In practical applications, the matrix operations are essentially the multiplication and addition operations between elements. To quantify what we mean by this, the amounts of multiplication and addition operations of the two algorithms (i.e., GA-MUSIC and Euler angle) are computed, respectively, as shown in Table2. Suppose that one matrix transposition is regarded as one multiplication or addition operation. As is known to all, the multiplication between two 3 × 3 matrices requires 9 × 3 multiplications and 9 × 2 additions. For convenience, the multiplication operation and the addition operation are collectively referred to as the operation. Then, the calculation of (28) contains 2 × 9 × 3 + 2 × 9 × 2 operations. For the conformal array consisting of M
×
N elements, the transformation between different coordinates involves 91
×
M
N operations. In contrast to the Euler angle, GA-MUSIC avoids the cumbersome matrix transformations.Table 2
The computational complexity of GA-MUSIC and Euler angle.
Multiplications
Additions
Transpositions
Operations
Euler angle
2
×
9
×
3
×
M
N
2
×
9
×
2
×
M
N
M
N
91
×
M
N
GA-MUSIC
4
×
M
N
2
×
M
N
0
6
×
M
NFrom (20)–(22), we know that e
y
m
n and e
z
m
n are independent of e
x. Moreover, e
x
m
n can be obtained directly from (20) without extra computations. Thus, (21)-(22) can be rewritten in 2 × 2 matrix. While applying the rotor to construct the array manifold, the computational process is equivalent to 2 × 2 matrix multiplied by 2 × 1 vector. Then, the operations for each element involve 4 multiplications and 2 additions. The total amount of operations is 6
×
M
N
.Thus, the latter method largely decreases the computational complexity.To sum up, the Euler rotation and its matrix representation cannot intuitively display the whole procedure. Moreover, as the configuration of the conformal array becomes more irregular and complex, the level of complexity involved in the transformations and the number of calculations required increase significantly.


## 4. Simulation Results

In this section, Monte-Carlo simulation experiments are employed to verify the effectiveness of the GA-MUSIC algorithm. The array structure is shown in Figure4. We select M and N as 4 and 2, respectively, for the validations. The pattern of single short dipole is shown in Figure 6(a), while the cylindrical conformal array pattern is displayed in Figure 6(b). Compared with the existing methods, adopting the geometric algebra will be more simple and intuitive. For the solved array pattern or manifold is the superposition of various elements, the proposed method can be applied to arbitrary array structures. This also proves the superiority of the rotor in solving conformal problems.Figure 6
Pattern analysis: (a) pattern of single short dipole and (b) pattern of cylindrical conformal array.
(a)
(b)Next, the performance of the GA-MUSIC algorithm is to be verified. We make some notations. Firstly, the absolute value of the differences between the estimated mean and the true value is regarded as the deviation. Secondly, the RMSE is utilized as the performance measure. Under these premises, 100 independent simulation experiments are carried out. The RMSE is defined as(31)
RMSE
=
1
100
∑
i
=
1
100
θ
^
i
-
θ
i
2
+
φ
^
i
-
φ
i
2
,where θ
^
i
,
φ
^
i are the estimates of elevation angles and azimuth angles, respectively, at ith run. Wang’s method [2] and Qun’s method [4] are included for comparison.Provided that there are three signals that can be received, the incident angles are (15°, 40°), (35°, 10°), and (60°, 65°), respectively. The corresponding polarization auxiliary angle and the polarization phase difference are (20°, 25°), (50°, 45°), and (65°, 65°). The snapshot,K, is selected as 100. Figure 7(a) shows the simulation results of the GA-MUSIC algorithm. The position of the spectrum peak represents the corresponding signal DOA. Intuitively, the estimation accuracy of the GA-MUSIC algorithm is high.Figure 7
The performance of GA-MUSIC algorithm: (a) DOA estimation using GA-MUSIC and (b) the spatial spectra of multiple signals.
(a)
(b)We assume that there are nine distinct signals impinging on the cylindrical conformal array at the same time. The spatial spectra are shown in Figure7(b). It is obvious that the method still works effectively when the number of incident signals is larger than that of the array elements.Figure8 displays the performance with a varying SNR from 0 dB to 30 dB. Among them, Figure 8(a) displays the relationship between SNR and deviation of DOA estimation, while Figure 8(b) reveals the Cramer-Rao Lower Bound (CRLB) and the RMSE versus the SNR, respectively. It is clear that the deviation varies inversely with SNR. The higher the SNR, the lower the deviation. The trend of the variance of the CRLB with SNRs is the same as the RMSE, which is expected. Obviously, in the engineering design, the higher the SNR, the better the estimation performance that we can obtain. Note that, since statistical data have some randomness, the simulation curves in Figure 8(a) are not smooth and do not decline monotonically.Figure 8
The DOA estimation accuracy withK
=
100: (a) deviation of estimates versus SNR and (b) the RMSE and CRB versus SNR.
(a)
(b)To proceed further, we increase the number of snapshots to 200 and leave the other conditions unchanged. The respective results are shown in Figure9. Compared with Figure 8, both the deviation and CRB were improved. If we choose the point at some SNR, we can find that the CRB of Figure 8 is nearly twice as much as of Figure 9. For example, the CRB of Target 1 is −82.58 dB when the SNR is 20 dB as shown in Figure 8(b), while, in Figure 9(b), the value is −86.08 dB. In fact, these improvements can be predicted from the derivation of CRLB. For details on the specific derivation process, see the work of Stoica and Nehorai [15, 16]. The number of snapshots can be extracted from the Fisher information matrix. Moreover, the CRLB is found as the element of the inverse of that matrix. So, we can conclude that CRLB is inversely proportional to K. Thus, the estimation precision will be higher.Figure 9
The DOA estimation accuracy withK
=
200: (a) deviation of estimates versus SNR and (b) the RMSE and CRB versus SNR.
(a)
(b)Figure10 shows the improvements of GA-MUSIC over existing algorithms, such as Qun’s method [4] and Wang’s method [2]. We study the performance with a varying SNR from 0 dB to 30 dB and the performance with the number of snapshots varying from 100 to 1000, respectively. Without loss of generality, we select the first source (T1) and the second source (T2), respectively, to verify it. As shown in Figure 10(a), the proposed method outperforms Qun’s method by exploiting the polarization information of the received data. Moreover, the elliptic covariance matrix of the GA-MUSIC is nonzero which increases the information utilization rates as well.Figure 10
Improvements of GA-MUSIC: (a) RMSE versus SNR with the snapshots being 100 and (b) RMSE versus the snapshots with the SNR fixed at 20 dB.
(a)
(b)The RMSE of the GA-MUSIC is close to that of the Wang’s method. The reason is that the estimation accuracy mainly depends on the steering vector. Both the Euler rotation angle and the geometric algebra can be used to derive the steering vector. As previously mentioned, the array manifold was obtained by using the rotor in this paper. However, in [2], Wang derived the array manifold by means of the Euler rotation angle. Therefore, in this case, both methods exhibit the same performance. From Section 3.2, we clearly know the GA-MUSIC is superior to the Euler rotation angle in computational complexity.Figure10(b) illustrates the RMSE versus the number of snapshots with the SNR fixed at 20 dB. Compared with Figure 10(a), we can draw similar conclusions.


## 5. Conclusion

In this paper, we proposed a novel algorithm named GA-MUSIC to estimate the DOA for cylindrical conformal array, which combines the geometric algebra with MUSIC. Compared with existing methods based on Euler rotation angle, it avoids complex matrix transformations and largely decreases the computational complexity. In contrast to the method introduced in the literature [8], our presented method has a strong commonality which can be used for arbitrary array structure. In addition, it can also be suited for the case that the number of polarized signals is larger than that of the array elements. At last, the simulated results verify the effectiveness of the GA-MUSIC algorithm.



---

*Source: 7832475-2016-12-04.xml*