Professional Documents
Culture Documents
Manikas A. - Differential Geometry in Array Processing-ICP (2004)
Manikas A. - Differential Geometry in Array Processing-ICP (2004)
ARRAY
PROCESSING
This page intentionally left blank
Athanassios Manikas
Imperial College London, UK
DIFFERENTIAL GEOMETRY IN
ARRAY
PROCESSING
Distributed by
World Scientific Publishing Co. Pte. Ltd.
5 Toh Tuck Link, Singapore 596224
USA office: 27 Warren Street, Suite 401–402, Hackensack, NJ 07601
UK office: 57 Shelton Street, Covent Garden, London WC2H 9HE
For photocopying of material in this volume, please pay a copying fee through the Copyright
Clearance Center, Inc., 222 Rosewood Drive, Danvers, MA 01923, USA. In this case permission to
photocopy is not required from the publisher.
ISBN 1-86094-422-1
ISBN 1-86094-423-X (pbk)
Printed in Singapore.
To my wife, Eleni
v
This page intentionally left blank
July 6, 2004 9:31 WSPC/Book Trim Size for 9in x 6in fm
Preface
During the past few decades, there has been significant research into
sensor array signal processing, culminating in the development of super-
resolution array processing, which asymptotically exhibits infinite resolu-
tion capabilities.
Array processing has an enormous set of applications and has recently
experienced an explosive interest due to the realization that arrays have a
major role to play in the development of future communication systems,
wireless computing, biomedicine (bio-array processing) and environmental
monitoring.
However, the “heart” of any application is the structure of the employed
array of sensors and this is completely characterized [1] by the array mani-
fold. The array manifold is a fundamental concept and is defined as the locus
of all the response vectors of the array over the feasible set of source/signal
parameters. In view of the nature of the array manifold and its signifi-
cance in the area of array processing and array communications, the role
of differential geometry as the most particularly appropriate analysis tool,
cannot be over-emphasized.
Differential geometry is a branch of mathematics concerned with the
application of differential calculus for the investigation of the properties of
geometric objects (curves, surfaces, etc.) referred to, collectively, as “mani-
folds”. This is a vast subject area with numerous abstract definitions,
theorems, notations and rigorous formal proofs [2,3] and is mainly confined
to the investigation of the geometrical properties of manifolds in three-
dimensional Euclidean space R3 and in real spaces of higher dimension.
However, the array manifolds are embedded not in real, but in
N -dimensional complex space (where N is the number of sensors). There-
fore, by extending the theoretical framework of R3 to complex spaces, the
vii
July 6, 2004 9:31 WSPC/Book Trim Size for 9in x 6in fm
Acknowledgments
Keywords
Contents
Preface vii
1. Introduction 1
1.1 Nomenclature . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 Main Abbreviations . . . . . . . . . . . . . . . . . . . . . . 4
1.3 Array of Sensors — Environment . . . . . . . . . . . . . . . 5
1.4 Pictorial Notation . . . . . . . . . . . . . . . . . . . . . . . 8
1.4.1 Spaces/Subspaces . . . . . . . . . . . . . . . . . . . 8
1.4.2 Projection Operator . . . . . . . . . . . . . . . . . . 9
1.5 Principal Symbols . . . . . . . . . . . . . . . . . . . . . . . 11
1.6 Modelling the Array Signal-Vector and Array Manifold . . 11
1.7 Significance of Array Manifolds . . . . . . . . . . . . . . . . 17
1.8 An Outline of the Book . . . . . . . . . . . . . . . . . . . . 19
ix
July 6, 2004 9:31 WSPC/Book Trim Size for 9in x 6in fm
Contents xi
Bibliography 215
Index 217
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Chapter 1
Introduction
1
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
The four types of problem described above are inter-related and the
solution to one problem may result in a partial or complete solution to
another. For example, the successful operation of all parametric parameter
estimation algorithms requires solving firstly the detection problem (i.e.
a priori knowledge of the number of emitters present). Furthermore, once
the number and directions-of-arrival (DOAs) of signals received at the array
site are estimated by solving the detection and direction-finding problems,
nulls may be readily placed along the directions of the unwanted signals,
hence achieving interference cancellation.
The applications of arrays in various scientific disciplines (such as the
ones already mentioned) are extensive and suffice to reveal the multi-
dimensional significance of the array concept. For instance, although array
processing has been extensively used in high frequency communications in
the past, the explosive growth in demand for cellular services in recent years
has placed it at the centre of interest. Spatial diversity is considered to be
one of the most promising solutions for increasing capacity and spectral
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Introduction 3
1.1 Nomenclature
It is assumed that the reader is familiar with the fundamentals of vector and
matrix algebra. In this book, for typographical convenience, matrices will be
denoted by blackboard bold symbols (e.g. A, T, I) or, in the absence of a cor-
responding blackboard bold symbol, by boldface (e.g. r, k, Γ). Any under-
lined symbol will represent a column vector, e.g. A, a, a. Derivatives with
respect to a general parameter p will be denoted with a “dot” (e.g. ȧ), while
the “prime” symbol (e.g. a ) will be reserved for differentiation with respect
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Introduction 5
FOV Field-of-View
SNR Signal-to-Noise Ratio
RMS Root Mean Square
ULA Uniform Linear Array
UCA Uniform Circular Array
Fig. 1.1 Relative geometry between a far-field emitting source and an array of sensors.
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Introduction 7
Fig. 1.2 Examples of array geometries: (a) linear array, (b) uniform circular array,
(c) “cross” array.
It can be easily deduced from the geometry of the problem that a lin-
ear array is incapable of distinguishing directions which are symmetric
with respect to the array line or which have the same elevation. If a
360◦ FOV is instead required, then a planar or 2D array should be
employed.
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
1.4.1 Spaces/Subspaces
Because the visualization of a space greater than a 3-dimensional real
space is impossible, the notation shown in Fig. 1.3 is employed to pro-
vide some illustrative aid. However, this pictorial representation should be
used with care.
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Introduction 9
Fig. 1.3 Illustrative Notation: (a) denotes a N -dimensional complex (or real) obser-
vation space. Note that any vector in this space has N elements; (b) denotes a one-
dimensional subspace/space spanned by the vector a; (c) denotes an M -dimensional
subspace/space (with M ≥ 2) spanned by the columns of the matrix A.
Then:
P⊥
A = the projection operator onto L[A]
⊥
= IN − PA (1.12)
(N × N ) matrix
N.B.: Properties of PA PA PA = PA (1.13)
P = PH
A A
Fig. 1.4 Illustrative representation of the subspaces L[A] and L[A]⊥ lying in an
N -dimensional complex space.
Fig. 1.5 Illustrative representation of the projection of any vector x ∈ H onto L[A] and
onto L[A]⊥ .
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Introduction 11
(3) For any vector x ∈ L[A] the following expressions are valid
PA x = x
(1.14)
P⊥
A x = 0N
p, q generic parameters
θ azimuth angle
φ elevation angle
α alpha — cone angle
β beta — cone angle
s arc length
N number of sensors
r (3 × N ) real matrix with columns the sensor locations
a array manifold vector
lm array manifold length
la array aperture
A array manifold curve
M array manifold surface
H space/subspace
la = max
ri − rj
(1.15)
∀i,j
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
where mi (t) represents the complex envelope (message) of the signal and
exp(j2πFc t) represents the carrier. If τik represents the propagation delay
of the ith signal between the phase reference location and the kth sensor
then the signal arriving at the kth sensor will be a delayed version of the
signal given in Eq. (1.16), expressed as shown below
This propagation delay τik is a function of the DOA of the ith signal and
the position of the kth sensor with respect to the reference point. Indeed,
with reference to Fig. 1.6, τik can be derived as
ik
Pu r
i k
τik = =
c c
rTk PTui Pui rk
=
c
rk Pui rk
T
= (using Eq. (1.13))
c −1 T
rTk ui uTi ui ui r k
=
c
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Introduction 13
Fig. 1.6 Geometry of a travelling plane wave, relative to the kth sensor and the array
reference point.
2
rTk ui uTi rk rTk ui rTk ui
= = =
c c c
i.e.
rTk ui
τik = (1.18)
c
It is reasonable to assume that the envelope (baseband signal mi (t)) does
not change significantly as it traverses the array. Hence τik τmax , ∀i, k
and mi (t − τmax ) mi (t) has been used in Eq. (1.17), where τmax is the
maximum possible time for a signal to traverse the array. This is due to the
narrowband assumption, i.e. the highest frequency in the message signal
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
is much less than the carrier frequency. By taking into account the above
assumption, as well as the complex response (gain and phase) γk (θi , φi ) of
the kth sensor at frequency Fc in the azimuth-elevation direction (θi , φi ),
the received signal at the kth sensor due to the ith emitting source can be
written as
rT · u
γk (θi , φi ) · mi (t) · exp −j2πFc k i exp(j2πFc t)
c
= γk (θi , φi ) · mi (t) · exp(−jrTk · ki ) exp(j2πFc t) (1.19)
Using vector notation, the output from all the N sensors, can be
expressed as
T
x(t) = x1 (t), x1 (t), . . . , xN (t)
M
= mi (t) ai + n(t) (1.21)
i=1
T
where n(t) = n1 (t) n2 (t) · · · nN (t) ∈ CN ×1 denotes the baseband
additive white Gaussian noise of power σ . The (N × 1) complex vector
2
ai is known as the array manifold vector associated with the ith emitting
source (signal), representing the complex array response to a unit amplitude
plane wave impinging from direction (θi , φi ), compactly written as
ai a(θi , φi ) = γ(θi , φi ) exp −j [rx , ry , rz ]k(θi , φi ) (1.22)
=rT
where
γ(θi , φi ) = [γ1 (θi ,φi ), γ2 (θi , φi ), . . . , γN (θi , φi )]T ∈ CN ×1
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Introduction 15
As can be seen from Eq. (1.22), a shift in the reference point of the array
coordinate system results in a manifold vector ar (θ, φ) which is related to
the initial manifold vector:
where g(θ, φ) is a complex number with |g(θ, φ)| = 1, ∀(θ, φ). In fact,
g(θ, φ) = exp(−jrT0 k(θ, φ)) where r0 is the translation vector of the ref-
erence point from the centroid of the array. For simplicity, all arrays
henceforth will have their centroid as the coordinate reference point, i.e.
g(θ, φ) = 1, ∀(θ, φ).
Note that in the case of a linear array of isotropic sensors Eq. (1.22)
simplifies to
It can easily be seen that r(θ) in Eq. (1.27) represents the positions of
the projections of the sensor locations of the planar array onto the line of
azimuth θ (see Fig. 1.7), and is known as the “Equivalent Linear Array”
(ELA) along the direction θ. Notice the similarity between the response
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Fig. 1.7 Equivalent Linear Array (ELA) associated with the azimuth angle θ.
vectors of the linear and planar arrays in Eqs. (1.26) and (1.27) respectively.
Using Eq. (1.22) the observed array signal-vector x(t) ∈ CN ×1 of Eq. (1.21)
can be written concisely as
A = exp(−jrT k) (1.31)
Introduction 17
It is clear from Eqs. (1.22), (1.26) or (1.27) that the manifold vector for
a particular direction contains all the information about the geometry
involved when a wave is incident on the array from that direction. By
recording the locus of the manifold vectors as a function of direction, a
“continuum” (i.e. a geometrical object such as a curve or surface) is formed
lying in an N -dimensional space. This geometrical object (locus of manifold,
or response, vectors) is known as the array manifold. The array manifold
can be calculated (and stored) from only the knowledge of the locations and
directional characteristics of the sensors. Thus, according to Schmidt [1],
“the array manifold completely characterizes any array and provides a rep-
resentation of the real array into N -dimensional complex space.”
The significance of the array manifold concept becomes apparent when it
is realized that all subspace-based parameter estimation algorithms involve
searching over the array manifold for response vectors which satisfy a
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
given criterion. For example, in the case of the MUSIC algorithm [7],
the manifold is searched for the array manifold vectors which are (nearly)
orthogonal to the estimated noise subspace. In other words the manifold
is searched for vectors belonging to the subspace L[A] i.e. the subspace
spanned by the columns of the matrix A = [a1 , a2 , . . . , aM ] (the signal sub-
space). Figure 1.8 illustrates a single-parameter manifold curve embedded
in N -dimensional complex space together with the subspace L[A] associated
with two sources. These two sources can be asymptotically resolved, irre-
spective of their angular separation — hence the title “superresolution.” It
is clear that the accuracy of the DOA estimates is dependent on how accu-
rately the noise subspace is estimated; which in turn is a function of the
observation interval and the signal-to-noise ratio (SNR). Thus the DOAs
can be exactly estimated from the array output only on the basis of either
an infinite number of snapshots (i.e. L → ∞) or, an infinite signal-to-noise
ratio scenario (i.e. n(t) = 0) — asymptotic conditions.
Another determining factor which might not be initially apparent is
the shape of the array manifold. A little thought reveals that a poten-
tially unresolvable situation might arise if two response (manifold) vectors
corresponding to different DOAs happen to be identical. For instance, this
“abnormality” arises when a single-parameter array manifold crosses upon
itself, as shown in Fig. 1.9. Such phenomena are commonly referred to as
“ambiguities” and are a direct consequence of the array geometry dictating
Introduction 19
the behavior of the array manifold. As will be discussed below, the behavior
of the array manifold, in particular its local geometry, plays an important
role in handling “abnormalities” and also in defining the capabilities (e.g.
the resolving power) of an array system.
are examined for non-linear arrays (2D and 3D geometries). The significance
of these parametrizations is demonstrated by a number of representative
examples. Furthermore, properties such as Gaussian and geodesic curva-
tures are defined and their implications with regards to isometric mappings
are discussed.
In the next two chapters of the book (Chapters 6 and 7) the funda-
mental effects of the array geometry behavior on the performance of the
system irrespective of the type of algorithm used, are studied. The array
structure is incorporated into the observed array signal-vectors x(t) (and
therefore to the estimation problem) through the array manifold. Thus the
geometry of an array plays a crucial role in dictating the shape, properties
July 9, 2004 13:51 WSPC/Book Trim Size for 9in x 6in chap01
Introduction 21
Chapter 2
ȧ(p) = 0N , ∀p ∈ Ωp (2.2)
and
ds
ṡ(p) = ȧ(p) (2.4)
dp
22
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
da(s)
da(p)/dp
ȧ(p) ṡ(p)
a (s) =
=
ds
ds/dp
= ṡ(p) = ṡ(p) = 1, ∀s (2.5)
which is a result used next in defining the unit tangent vector u1 (s) (see
Fig. 2.3)
u1 (s) a (s) (2.6)
Fig. 2.3 Arc length s is an “invariant” parameter — tangent vector u1 (s) has unit
norm ∀s.
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
Fig. 2.4 The parameter p is not an invariant parameter — i.e. speed = ȧ(p) = ṡ(p) =
variable = 1.
moving object while the total distance travelled, s(p), can be obtained by
integrating speed with respect to “time” p. Parametrization in terms of s,
i.e. a(s), indicates that the “velocity” vector (tangent vector) a (s) always
has constant length equal to 1 (moving object with unity speed).
Note that, in general, it is not necessary to mention the origin of
the arc length (s = 0) since most concepts are defined only in terms
of the derivatives of a(s). Furthermore, the rate of change of arc length
ṡ(p) is a local property of the curve and, as will be demonstrated in
Chapter 8, plays a crucial role in dictating the resolution/detection capa-
bilities offered by the array manifold and consequently by the geometry of
the array.
As was mentioned before, the tangent vector a (s) is of unity length
∀s, indicating that s is an “invariant” parameter. In this case the norm
of the second derivative of the manifold vector (parametrized in terms of
s) a (s) measures how fast the curve pulls away from the tangent line
at s, in the neighborhood of s (see Fig. 2.5). In other words, it measures
the rate of change of the angle which neighboring tangents make with the
tangent at s. The number a (s) = κ1 (s) is called the first curvature of
the manifold curve at s. As we will see next, although many curvatures can
be defined at a point on the manifold curve A ∈CN , the first curvature is
the most important one and for this reason this is also known as “principal”
curvature.
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
Fig. 2.5 How fast the curve pulls away from the tangent line, at point s, is measured
by a (s) with a (s) = κ1 (s).
We have seen that the arc length s is the invariant parameter of a manifold
curve of an array of N sensors. With the array manifold curve embedded in
an N -dimensional complex space (or equivalently in a 2N real space) and
parametrized in terms of the arc length s, it is essential to attach/define
to the running point a = a(s) on the curve a number of curvatures as well
as a continuous differential and orthonormal system of coordinates u1 (s),
u2 (s), . . . , u2N (s).
(1) Symmetric: m = N , i.e. all the sensors occur in symmetric pairs about
the origin.
(2) Non-Symmetric: m < N
(a) Partially symmetric: 0 < m < N , i.e. at least one sensor has a
symmetric counterpart about the origin;
(b) Fully asymmetric: m = 0, i.e. no sensor has a symmetrical coun-
terpart about the origin.
and is derived from a fixed known frame U(0), (i.e. at s = 0 say) by rotation,
using the transformation matrix F(s) ∈ Cd×d , i.e.
To summarize, the d×d real matrix F(s) is known as the “frame matrix”
while the set of coordinates represented by the N × d complex matrix U(s)
is the “moving frame.”
where U (s) = [u1 (s), u2 (s), . . . , ud (s)] ∈ CN ×d , d is the dimensionality of
the subspace in which the curve is embedded. The matrix C(s) ∈ Cd×d is
termed the Cartan matrix, which is a real skew-symmetric matrix of the
curvatures defined as follows:
0, −κ1 (s), 0, ..., 0, 0
κ (s), −κ2 (s), . . . ,
1 0, 0, 0
0, κ2 (s), 0, ..., 0, 0
C(s) .. .. .. .. .. .. (2.14)
. . . . . .
0, 0, 0, ..., 0, −κd−1 (s)
0, 0, 0, . . . , κd−1 (s), 0
Starting with Eq. (2.13) and then using Eq. (2.10), we have
This is a first order differential equation of the frame matrix F(s) with
initial condition F(0) = Id . Therefore, its solution is
where expm(·) denotes the matrix exponential. Equation (2.16) provides the
relationship between the frame matrix F(s) at the running point s and the
curvatures (Cartan matrix) of the manifold attached to this point, which
can also be used to rewrite Eq. (2.10) as follows:
Finally, Eq. (2.15) implies that the Cartan matrix can always be written
as a function of the frame matrix F(s) as follows
where the 3rd property of F(s) given in Table 2.1 has been used. Remember
that the Cartan matrix, as a formation of the curvatures, contains all the
information about the local behavior and shape of the manifold curve.
This, however, is only true for the broad class of arrays which are symmetric
with respect to their centroid, in which case m = N . In general, if the array
is nonsymmetric (m < N ), then the orthogonality is not valid although the
complex vectors ui (s) have constant (unity) magnitude. In particular, the
inner product of a complex vector of constant norm with its derivative is
purely imaginary, i.e.
uH
i (s)ui (s) = imaginary (2.19)
as follows:
Definition 2.1 Wide sense and narrow sense orthogonality of two com-
plex vectors of unity magnitude:
Thus, for a general linear array, that is, for a non-symmetric array, the
wide-sense orthogonality should be used, i.e.
Indeed
C(s) = FT (s) · F (s) = FT (s) · Id · F (s) = FT (s) · Re(UH (0) · U(0)) · F (s)
=Id
= Re FT (s) · UH (0) · U(0)F (s) = Re(UH (s)U (s))
=UH (s) =U (s)
to the curvatures by
(
N −1
−
= κN
(s) (2.25)
=1
A {a(p) ∈ CN , ∀p : p ∈ Ωp } (2.26)
where
a(p) = exp(−j(πr cos p + v)) (2.27)
calculated as:
This explains the better resolution of sources at broadside and the better
resolving power of the 7-sensor standard ULA, as evidenced by Fig. 2.10.
It can therefore be concluded that the rate of change of arc length of the
array manifold plays a prominent role in determining the resolving power
of the array, as will be further detailed in Chapter 8.
Another important parameter of the array manifold curve is its total
length,
lm s(π) = 2π r (2.33)
In general, it can be shown that the ith coordinate vector and curvature
of the manifold of a linear array of isotropic sensors can be calculated
according to the following theorem.
fix( i−1
(j)
i
2 )+1
fix( 2 )+1
1
κi =
(−1)n−1
b )
r i−2n+3
(2.39)
κ1 κ2 · · · κi−1
i,n
n=1
where
) r
r = r (normalized sensor positions)
κ = ||)2
r ||
1
κi−1 = 0
sum(r) = 0 (i.e. phase reference = array centroid)
i−2n+5
i−2n+3
i−1
bi,n = ··· κ2m1 κ2m2 · · · κ2mn−1 ; i, n > 2
m1 =1 m2 =2+m1 mn−1 =2+mn−2
(2.40)
with
bi,1 = 1; i 1
* 2
i−1
bi,n = κm i > 1
m=1
or recursively by
bi,n = bi−1,n + κ2i−1 bi−2,n−1 , i > 2, n > 1 (2.41)
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
where the non-zero entry is in the same position as the sensor at the centroid
in the vector of sensor locations )
r.
are mutually orthogonal and that the same holds for the even-indexed coor-
dinate vectors forming the matrix
Ueven (s) = u2 (s), u4 (s), u6 (s), . . . , udeven (s) (2.46)
where deven = fix( d2 ) and dodd = fix( d−1 2 ) + 1. This can be expressed in a
compact form as
.
UH
odd (s)Uodd (s) = Idodd
(2.47)
UHeven (s)Ueven (s) = Ideven
with .
UHodd (s)Ueven (s) = Ododd ×deven
(2.48)
Re{UH odd (s)Ueven (s)} = Ododd ×deven
Note also that, at broadside (i.e. p = 90◦ ), Uodd (s) is an imaginary matrix
while Ueven (s) is a real matrix. In addition, in the special case of “symmetric
arrays”, the coordinate vectors form a mutually orthogonal frame
since the sum of odd powers of the sensor locations is equal to zero.
1
κ1 > √ (2.50)
N
Table 2.3 The first eight curvatures of the manifold of uniform linear arrays.
N =3 N =4 N =5 N =6 N =7 N =8 N =9 N = 10
Fig. 2.11 Curvatures of uniform linear arrays as functions of the number of sensors N .
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
The above relation also stems from the following inner product aH u2
in conjunction with the expression
such arrays depends both on the relative spacing and the first two deriva-
tives of each elemental pattern γk (p). In the case of the manifold of a
4-element uniform linear array with non-isotropic, but identical, sensors
with directional pattern γk (p) = cos p, the variation of the first curvature
with the azimuth direction is shown in Fig. 2.13. The constant first curva-
ture of the corresponding isotropic array is also shown in the same figure.
It is seen that there exist two bearings where the first curvature becomes
equal to zero. From the variation of the first curvature, it is possible to
deduce the geometrical object to which the manifold best fits. In the case
considered, it is apparent that the effect of the sinusoidal elemental pattern
is to deform the hyperhelix to a geometrical figure resembling an “eight”
with a double point at the origin of the coordinates of CN .
It is noted that the manifolds of uniform linear arrays with an odd number
of sensors and half-wavelength spacing consist of one winding (round) since
the manifold vectors corresponding to the manifold boundaries coincide
(i.e. a(0◦ ) = a(180◦ )). In addition, the manifolds of uniform linear arrays
with an even number of sensors and half wavelength spacing, consist of one
half winding (half round) since the manifold vectors corresponding to the
manifold boundaries are opposite (i.e. a(0◦ ) = −a(180◦ )). This situation
has forced the separation of a hyperhelix corresponding to a general array
with an odd number of sensors into consecutive equal arcs of one winding
and a hyperhelix corresponding to an array with an even number of sensors
into consecutive equal arcs of half-winding.
The following theorem is concerned with the estimation of the number
of windings (N odd) and the number of half windings (N even) existing in
the manifold of an arbitrary linear array.
manifold length=lm
2π r
nw = (2.51)
lw
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
Fig. 2.13 The principal curvature as a function of azimuth (top) and its polar repre-
sentation (κ, p) (bottom) for a 4-element ULA of directional sensors — γk (p) = cos p.
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
where lw represents the arc length of half winding. For an array with an
odd number of sensors, the function ξw (Eq. (2.52)) does not generally have
a root at the point on the manifold that corresponds to half the length of
one manifold winding.
Similarly, the definition of one winding is restricted to linear arrays with
an odd number of sensors and
a(0) = −a(lw )
U(0) = −U(lw )
F(0) = Id = −F(lw )
⇒ Tr(F(0)) = −Tr(F(lw )) = d
where lw represents the arc length of one winding. For an array with an
even number of sensors, the function ξw does not generally have a root
at the point on the manifold that corresponds to twice the length of a
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
Fig. 2.14 Tr(F(s)) and ξw (s), given by Eq. (2.52), for a 5-element and a 6-element ULA
of half wavelength spacing. Note that A1 = 5 = A2 and B1 = 6 = −B2.
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
In the previous section it has been shown that the shape of the manifold
of a general linear array of omnidirectional sensors is a complex circular
√
hyperhelix lying on the complex N -dimensional sphere with radius N . In
this section, an attempt is made to identify the orientation of the array
manifold (hyperhelix) on the complex N -dimensional space. Therefore, it is
necessary to establish a measure of the orientation of the manifold curve in
CN using the concept of array inclination angle. This is intuitively defined
as follows:
Having defined the matrix Usub (s) as a submatrix of Ueven (s) and
because of the orthonormality of columns of the matrix Ueven (s) (see
Eq. (2.47)) it is apparent that Eq. (2.54) can be rewritten as follows:
0
1
1 1
N
=2
2
cos ζinc |aH (s) ui (s)| (2.55)
N i=even
Fig. 2.15 Inclination angle for a number of linear arrays with 4 and 6 sensors.
different arrays resulting from uniform linear arrays with 4 and 6 elements,
of half-wavelength spacing, by varying the rightmost (leftmost) spacing
from 0.5 half-wavelengths to 3 half-wavelengths is shown.
Equation (2.55) shows that the inclination angle is directly related to the
inner products of a manifold vector aH (s) with the coordinate vectors u2 ,
u4 , u6 , . . . , etc. This demonstrates that the relationship between aH (s) and
its associated coordinate vectors ui (s), ∀i is very important.
For example, the relationships between aH (s) and the first four coordi-
nate vectors, at a point, are:
aH (s)u1 (s) = j1TN )
r
1
a (s)u2 (s) = −
H
κ1
j & ' (2.56)
aH (s)u3 (s) = r
r3 − κ21 1TN )
1TN )
κ1 κ2
κ2
a (s)u4 (s) = −
H
κ1 κ3
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
If the array is also symmetric, then at every point along the manifold, the
corresponding response vector is orthogonal to the odd-indexed coordinate
vectors.
To provide a more general expression of the inner product
aH (s)ui (s); ∀i
Lemma 2.1 The inner product aH (s)ui (s), ∀i: For a linear array of
N sensors the inner products of a manifold vector a(s) with its d coor-
dinate vectors ui (s) are constant and given by the following expressions:
Note that the inner product of a manifold vector and the corresponding
normal vector (or second coordinate vector) is the negative of the radius
of curvature (i.e. −1/κ1 ) and hence the other interior products can be
considered as higher order “manifold radii.” Using the above lemma the
inclination angle of Eq. (2.55) can be rewritten as
Rsub
ζinc = arccos √ = constant (2.62)
N
Definition 2.4 Dual arrays: The dual array of a linear array, non-
symmetrical with respect to the array centroid, is the array with sensors
having positions given by rdual
i = −rN -i+1 for i = 1, . . . , N .
Using the above definition it can be seen that, although the manifolds
of two dual arrays are not identical, they have the same curvatures and
length. Thus, apart from the fact that they have opposite inclinations, they
have the same differential geometry properties.
Based on the above concepts the array manifold vector of Eq. (2.27) or
(2.35) can be written as
real vector
It is clear that Eq. (2.69) is the equation of a closed circle (assuming half-
wavelength spacing) with radius 1/κ1 for −rπ ≤ s ≤ rπ. If the two
spacings of the above array become unequal, κ2 does not vanish and the
manifold becomes a helix.
T
v(s) =v (s)
H
= Tr
U (0)P U(0) F(s)R RT F(s)
n
=Pn =V(s)
= Tr Pn V(s) (2.71)
Fig. 2.16 MUSIC spectrum based on Eq. (2.71) for a ULA of 7 sensors operating in the
presence of 3 emitting sources.
(2) If Rsub is the vector formed by the first N components of the manifold
radii vector, then
√
Rsub ≤ R = N (2.74)
From the above it can be deduced that, for any symmetric sequence of
N real numbers ri , Eq. (2.65) gives an expression
√ for a real hyperhelix
(on a real N -dimensional sphere with radius N ) of latitude
π − arccos √1 N odd
latitude = 2 N
0 N even
Finally the Frobenius norm of the Cartan matrix (i.e. CF ) is related
to the array symmetricity as follows:
(2.77)
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
2.7 Appendices
Assuming that
Re a(i−1) (s)H ui−1 (s) = κ1 (s)κ2 (s) · · · κi−2 (s) (2.79)
For convenience let us drop the parameter s. The argument of the first part
of Eq. (2.80) can be written as follows:
(a(i) )H ui
= (a(i−1) )H ui
& '
= (a(i−1) )H ui − (a(i−1) )H ui
& '
= (a(i−1) )H ui − (a(i−1) )H (κi ui+1 − κi−1 ui−1 )
=ui (see Eq. (2.12))
& '
= (a(i−1) )H ui − κi (a(i−1) )H ui+1 + κi−1 (a(i−1) )H ui−1
=κ1 κ2 ···κi−2
(2.81)
However, by using the fact that the vectors a(i−1) can be expressed as a
linear combination of coordinate vectors with no higher order than (i − 1),
Eq. (2.81) can be rewritten as follows:
H H
i−1
i−1
(2.81) = c j u j u i − κi cj uj ui+1 + κ1 κ2 · · · κi−2 κi−1
j=1 j=1
imaginary imaginary
(2.82)
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
where cj are scalar coefficients. Then, by taking the real part of both sides
and using the property that
j ui } = 0
Re{uH for j = i (2.83)
we have & 'H .
Re a(i) ui = 0 + 0 + κ1 κ2 · · · κi−2 κi−1 (2.84)
= j r̃ a (2.88)
It can be seen that, since the origin of the coordinates is at the array
centroid,
aH u1 (s) = 0 (2.89)
1
= j r̃ (jπ sin pr a) ·
π sin pr
ȧ
= −r̃ a
2
(2.90)
u1 (s)
u2 (s) =
u1 (s)
1 2
=− r̃ a (2.92)
κ1
Furthermore,
du2 dp
u2 (s) =
dp ds
1 −1
= − r̃2 ȧ (π sin p r)
κ1
=ṡ(p)
=u̇2
1 2 1
=− r̃ (jπ sin pr a) ·
κ1
π sin pr
ȧ
1 3
= −j r̃ a (2.93)
κ1
and
1 3
u2 (s) + κ1 u1 (s) = −j r̃ − κ1 r̃ a (2.94)
κ1
Hence, the second curvature (see Table 2.2) and the third coordinate
vector are
and
1
u3 (s) = (u (s) + κ1 u1 (s))
κ2 2
1 1 3
= −j · r̃ − κ1 r̃ a
κ2 κ1
1
= −j · (r̃3 − κ21 r̃) a (2.96)
κ1 κ2
In order to estimate the third curvature it is necessary to estimate the
terms u3 (s) and κ2 u2 (s). Indeed,
du3 dp
u3 (s) =
dp ds
1 −1
= −j · r̃3 − κ21 r̃ ȧ · (π sin pr)
κ1 κ2
=ṡ(p)
=u̇3
1
= −j · (r̃3 − κ21 r̃) (jπ sin pr a)(π sin pr)−1
κ1 κ2
=ȧ
1
= · (r̃4 − κ21 r̃2 ) a (2.97)
κ1 κ2
and
κ2 2
κ2 u2 (s) = − r̃ a (2.98)
κ1
which implies
1 4
u3 (s) + κ2 u2 (s) = r̃ − (κ21 + κ22 )r̃2 a (2.99)
κ1 κ2
By taking the magnitude of the above equation the third curvature is
estimated as
r̃5
=
2 2 2 3
−(κ1 + κ2 + κ3 )r̃
(2.102)
κ1 κ2 κ3
+(κ1 κ3 )2 r̃
r̃6
−(κ
1 2 2 2 2
1 + κ2 + κ3 + κ4
)r̃
= (2.104)
κ1 κ2 κ3 κ4
2
+ (κ1 κ3 ) + (κ1 κ4 )
2
2 r̃2
+(κ2 κ4 )
The above equations can be generalized for any index i, to give the
curvature κi and the vectors ui by the following expressions:
i−1
i−3
i−1
1
i+1 i−1 i−3
κi =
)
r −)
r · κn + )
2
r · (κn κ ) − · · ·
2
κ1 κ2 · · · κi−1
n=1 n=1 =n+2
(2.106)
July 6, 2004 9:29 WSPC/Book Trim Size for 9in x 6in chap02
or, equivalently,
i
fix( 2 )+1
1
κi =
(−1)n−1
bi,n )
r i−2n+3
(2.107)
κ1 κ2 · · · κi−1
n=1
and
fix( i−1
(j)i 2 )+1
Chapter 3
The locus of all array manifold vectors a(p, q) ∀(p, q) forms a surface known
as the array manifold. This is a two-parameter manifold embedded in an
N -dimensional complex space and is formally defined as
where Ω denotes the parameter space. The shape of the surface is very
important. In the same way that the shape of a curve A is uniquely
defined by its curvatures, the shape of a two-parameter manifold (surface)
may be quantitatively expressed in terms of intrinsic geometrical parame-
ters such as the Gaussian curvature KG (p, q) of the surface (see Fig. 3.1)
and the geodesic curvature κg (p, q) of the curves lying on the surface (see
Fig. 3.2).
More specifically, for a point a a(p, q) on the manifold surface, the
parameters of interest, intrinsic to the manifold surface, together with neces-
sary building blocks of surface differential geometry, will be identified in this
chapter and extended from three-dimensional real space to N -dimensional
complex space. In order to achieve this all the ramifications and subtleties
involved in this process have been taken into consideration.
Thus, it is first necessary to introduce the building blocks of surface
differential geometry such as the manifold metric, first fundamental coeffi-
cients and the Christoffel symbols.
However, to proceed, it is essential to start with two definitions which
will aid in the development of the basic theory. These are
59
July 9, 2004 13:52 WSPC/Book Trim Size for 9in x 6in chap03
Fig. 3.1 The Gaussian curvature KG (p, q) provides information about the local shape
of the surface M in the neighbourhood of a(p, q).
Fig. 3.2 The geodesic curvature κg provides information about a curve on a surface.
a = a(p, q) (3.2)
The vectors ȧp and ȧq at a specific point A on the manifold represent
the tangent vectors to the p- and q-parameter curves respectively passing
through A and also form a basis for the tangent plane to the surface at that
point, i.e.
Tangent plane = L [T] where T = [ȧp , ȧq ] ∈ CN ×2 (3.4)
The regularity condition of Eq. (3.3) ensures that a tangent plane exists at
all points on the surface. Note that the basis created by T is, in general,
not orthonormal.
is said to be the manifold metric. The elements gpp , gqq , and gpq (note
gqp = gpq ) of G are known, in Differential Geometry terms, as the first
fundamental coefficients (or metric coefficients) expressing the magnitudes
and inner products of the parameter-curve tangent vectors and entirely
describing the manifold properties of the surface. For instance, the angle
between the two parameter curves of the surface at a point (p, q) can be
expressed as
gpq
ȧp , ȧq = arccos √ (3.6)
gpp gqq
July 9, 2004 13:52 WSPC/Book Trim Size for 9in x 6in chap03
xT Gx 0, ∀x ∈ R2 (3.12)
which implies that at some points on the surface M this matrix can possibly
be singular in order to account for the fact that there may not exist a
unique tangent space (L[T]) at each point. That is, at a point (p, q) with
a non-unique tangent space, the matrix G is singular. These points are the
boundaries of the manifold surface.
2
(3.13)
4 5
dp
= dp dq G = dpT Gdp (3.14)
dq
with T
dp = dp, dq
is known as the first fundamental form and essentially represents the dis-
tance between two neighboring points a(p, q) and a(p + dp, q + dq) on the
manifold surface.
Two properties of I which may be readily established are
the Christoffel symbols of first kind represent the inner products between
the tangent vectors ȧp , ȧq (columns of T) and their derivatives (columns of
Ṫp and Ṫq ), and are defined as
Γi,jk Re ȧH i äjk with i, j, k = p or q (3.19)
Thus there are eight Christoffel symbols of first kind forming the two
Christoffel matrices of the first kind Γ1p and Γ1q defined as follows:
4Γ Γp,qζ
5
p,pζ
Γ1ζ Re T Ṫζ =
H
with ζ = p or q (3.20)
Γq,pζ Γq,qζ
It is obvious that all Christoffel symbols are functions only of the first
fundamental coefficients gij (i, j = p, q) and their derivatives. Thus, for
this reason they play a central role in expressing the intrinsic geometry of
a surface, which is described in the next section.
where T T (p, q) is the unique normal to the surface at point (p, q).
As the array manifold is embedded in CN (and not in R3 ) a unique
normal does not exist to each point and, therefore, Eq. (3.28) cannot be
defined for array manifolds. Here we will only consider the intrinsic for-
mula for the Gaussian curvature which is independent of the normal space.
For instance, an intrinsic formula for the Gaussian curvature, KG (p, q),
of a surface, can be shown to be of the form [17] (also known as the
July 9, 2004 13:52 WSPC/Book Trim Size for 9in x 6in chap03
curvature scalar ):
2
1 d gpq d2 gpp d2 gqq
2 2 − 2
−
1 dp dq dq
7 dp2 8
KG (p, q) = * ij Γi,pq Γi,pp (3.29)
det(G)
+ g det
i,j=p,q Γj,qq Γj,qp
The above expression may be more compactly represented in terms of the
Christoffel symbols:
d d
Γp,qq − Γp,qp
dp dq
1
KG (p, q) = + col2 {Γ2p }T col1 {Γ1q } (3.30)
det (G)
T
− col2 {Γ2q } col1 {Γ1p }
where coli {matrix} denotes the ith column of the matrix.
Here the following intrinsic formula is used, adapted from [19], and
expressed as a function of the metric G as follows:
√ √
det(G) q det(G) q
d gpp Γpq d gpp Γ pp
1 (3.31)
KG (p, q) = − 6 −
det(G) dp dq
where Γqpq , Γqpp are elements of the 2 × 2 real matrices Γ2p and Γ2q .
The sign of the Gaussian curvature of Eq. (3.31) or (3.29), provides
an indication of the local shape of the surface in that neighborhood.
For instance, the surface around a point (p, q) is locally:
• elliptic, if KG (p, q) > 0 (the whole neighborhood of the surface at the
point considered lies on one side of the tangent hyperplane),
• hyperbolic, if KG (p, q) < 0 (one part of the surface at the point considered
lies on one side of the tangent hyperplane and the other part on the other
side),
• either parabolic or planar (i.e. flat), if KG (p, q) = 0 (there is a straight
line of the surface, lying totally on the tangent space).
Example 3.1 The following lists a number of surfaces in R3 and their
respective Gaussian curvatures:
The converse of this is in general not true. However if two surfaces have
the same constant Gaussian curvature, then any two sufficiently small
neighborhoods of the surfaces are related by an isometric mapping. In
particular
At {a(t) ∈ CN , ∀t : t ∈ Ωt } (3.32)
where
a(t) a(p(t), q(t))
on the array manifold surface M, the arc length s is the total distance
travelled along the curve At from an “initial” point (p(t1 ), q(t1 )) to the
“current” point (p(t2 ), q(t2 )) under consideration. Its expression as a func-
tion of the manifold metric is, based on Eq. (3.16):
9
t2
dp T dp
arc length : s= G dt (3.33)
(curve on a surface) t1 dt dt
where
4 5T
T
dp dp dq
p = [p(t), q(t)] and, hence, = , (3.34)
dt dt dt
where s is the arc length of the curve under consideration. From Eq. (3.36),
a geodesic curve, i.e. a curve for which the geodesic curvature is invariably
July 9, 2004 13:52 WSPC/Book Trim Size for 9in x 6in chap03
equal to
9
det (G) qo
κg,p κg (p, qo ) = 3
Γpp (3.42)
gpp
In Section 3.4.1, we have seen that by examining the sign of the Gaussian
curvature KG of the array manifold it is possible to study its shape. It is
recognized that apart from the fact that the Gaussian curvature offers an
indication about the manifold shape, its importance is enhanced by the fact
that if it satisfies certain stringent conditions (valid for a very broad class
of arrays) the geodesic mappings of the individual manifold curves result in
a consistent mapping of the whole manifold on the parameter plane, called
the development of the manifold (see Fig. 3.3), which is defined as follows:
Definition 3.6 Development: Let At be a differentiable curve lying on
the manifold M expressed in terms of a parameter t as follows:
Furthermore, let s(t) be the arc length along At given by Eq. (3.33). The
development of At on the plane R2 is the plane curve
surface is mapped to another surface and the mapping is such that geodesic
curves in one surface are mapped onto geodesic curves in the other sur-
face (geodesic mapping) then if one surface has a constant Gaussian cur-
vature the other must also have a constant Gaussian curvature, in order
to have a distortionless mapping. Surfaces of constant KG are known as
“developable” surfaces [17]. Thus the condition for the development1 to
exist is KG to be constant ∀(p, q), since a manifold surface M with a vary-
ing Gaussian curvature KG cannot be adequately mapped onto the plane.
3.6 Summary
3.7 Appendices
A {a(s) ∈ CN , ∀s : s ∈ Ωs } (3.47)
1 Note that the “development” is on a plane and planes have constant (equal to zero)
Gaussian curvature.
July 9, 2004 13:52 WSPC/Book Trim Size for 9in x 6in chap03
where
a(s) a(p(s), q(s))
lying on M and let us define p(s) [p(s), q(s)]T where s is the arc length
of the curve. Then the tangent vector at a point (p(s), q(s)) on the curve
A is given as follows:
da(s)
a ≡ = T p
ds
where T = [ȧp , ȧq ] ∈ CN ×2 such that L[T] is the tangent space to
the manifold surface at (p(s), q(s)). Furthermore, let us define the matrix
B = [b1 , b2 ] ∈ CN ×2 as an orthonormal basis for the space L[T] with b1
chosen to be equal to a /a . This gives the expression:
where
9
det (G) gpq
1 −√
Ep = 6 gpp gpp ∈ R2×2 (3.49)
det (G) √
0 gpp
where ψ is the rotation angle between the tangent vector to the p-curve ap
and the tangent vector to the curve under consideration a :
: ;
Re{ȧH p a }
ψ = ȧp ,a = arccos
ȧp
a (3.51)
d (a ) dT dp d2 p
a = +T 2 (3.52)
ds ds ds ds
of the manifold vector at (p(s), q(s)) onto the subspace L[B] (note that
L[B] = L[T]), we have
PB a = κ1 b1 + κg b2 (3.53)
July 9, 2004 13:52 WSPC/Book Trim Size for 9in x 6in chap03
−1 H
where PB = B BH B B is the projection operator. The coefficient of b2 ,
that is κg , is defined as the geodesic curvature of the curve at (p(s), q(s)) and
represents the component of a in the direction of b2 . Hence its expression is
2 a }
κg = Re{bH (3.54)
The vector b = [0, 1]T will be used for notational purposes such that
b2 = Bb. Substituting the expressions of b2 and a into Eq. (3.54) gives,
after some manipulation:
<
dp dq dp d2 p
κg (s) = b H Ep G
T T T
Γ2p + Γ2q + 2 . (3.55)
ds ds ds ds
6
In Eq. (3.55) the term bT HT ETp G can be simplified to det(G) · [−p , q ],
providing
4 5T
6 −p
κg (s) = det (G) (Γ2p p + Γ2q q ) p + p (3.56)
q
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap04
Chapter 4
Non-Linear Arrays:
(θ, φ)-Parametrization of Array
Manifold Surfaces
where r is the matrix with columns the sensor location and k(θ, φ) is the
wavenumber vector given by Eq. (1.3). Thus by using the directional param-
eters (θ, φ), instead of the generic parameters (p, q) of the previous chapter,
the two-parameter array manifold M of Eq. (3.1) becomes
77
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap04
dk (θ, φ) d2 k(θ, φ)
k̇p = and k̈pq = (4.5)
dp dp dq
with
(p and q) = (θ or φ)
=Γθφθ
Γθθθ , − tan φ
Γ2θ =
Γφ ,
θθ 0
=Γφ
φθ
It should be noted from Eq. (4.9) that Γφθφ = 0. This simplifies the
Gaussian curvature of Eq. (3.31), for a 3D array of N -omnidirectional
sensors, to
:6 ;
1 d det(G) φ
KG (θ, φ) = 6 Γθθ (4.10)
det(G) dφ gθθ
Thus, in 3D-grid arrays the vectors rx , ry and rz are not only orthogonal
but also have the same magnitude. The following theorem is concerned
with the shape of the manifolds of 3D-grid arrays. To prove this theorem
it suffices to show that these arrays have constant and positive Gaussian
curvature.
Theorem 4.2 The manifold surface of a 3D-grid array of N -
omnidirectional sensors is spherical with radius πρ (and hence developable)
embedded in an N -dimensional complex space.
Example 4.1 The 8-element cube array with all sides equal to one half-
wavelength is a 3D-grid array
√ and has a spherical manifold with Gaussian
curvature (2π)−1 , or radius 2π half-wavelengths.
Table 4.1 summarizes the results of the array manifold parameters for
3D-grid arrays.
Table 4.1 Manifold parameters of 3D-grid arrays.
det(G) ρ4 π 4 cos2 φ
0, −1
Γ1θ ρ2 π 2 cos φ sin φ
1, 0
(Eq. (3.20), (ζ=θ))
2 2 −1, 0
Γ1φ ρ π cos φ sin φ
0, 0
(Eq. (3.20), (ζ=φ))
0, − tan φ
Γ2θ
cos φ sin φ, 0
(Eq. (3.22), (ζ=θ))
− tan φ, 0
Γ2φ
(Eq. (3.22), (ζ=φ))
0, 0
1
KG
(Eq. (3.31)) ρ2 π 2
A planar (2D) array can be seen as a special case of a 3D-array where all
the elements of the last (3rd) row of the matrix r are zeros (i.e. rz = 0N ).
In this case, the term rrT has the following form
L, 02 c1 , c2
where L = ∈ R2×2 (4.12)
0T2 , 0 c2 , c3
Fig. 4.2 Manifold Surface M of a Planar Array parametrized in terms of (θ,φ). This is
a developable manifold-conoid generated by rotating a φ-curve (geodesic curve) around
the apex.
π4
det(G) det(L) sin2 (2φ) ρ4 π 4 cos2 φ sin2 φ
4
g
θφ
, −gθθ tan φ
tan φ 0, −1
Γ1θ gφφ ρ2 π 2 cos φ sin φ
(Eq. (3.20), (ζ=θ)) , −gθφ tan φ 1, 0
tan φ
gθφ
−gθθ tan φ,
tan φ −1, 0
Γ1φ gφφ ρ2 π 2 cos φ sin φ
(Eq. (3.20), (ζ=φ)) −gθφ tan φ, 0, 1
tan φ
0, − tan φ 0, − tan φ
Γ2θ 1 1
(Eq. (3.22), (ζ=θ)) tan φ
, 0 tan φ
, 0
− tan φ, 0 − tan φ, 0
Γ2φ 1 1
(Eq. (3.22), (ζ=φ))
0, tan φ
0, tan φ
KG 0 0
(Eq. (3.31))
Based on the results presented in Tables 4.1 and 4.2 a number of com-
ments can be made. For instance, the Gaussian curvature of 3D-grid arrays
is always positive and constant whereas for planar arrays it is always zero
(constant). This implies that the array manifold of a 3D-grid array is iso-
metric with a sphere of radius ρπ while the manifold of a planar array is
isometric with a plane. In addition, for all 2D-grid and 3D-grid arrays, the
manifold metric G is a diagonal matrix. The determinant of this matrix is
given by the 2nd row of Tables 4.1 and 4.2 and can be easily proven for grid
4
arrays. For planar arrays, the proof of det(G) = π4 det(L) sin2 (2φ) is given
in Appendix 4.7.2. All these expressions indicate that the determinant of
the manifold metric G is independent of the azimuth angle θ.
The array manifold can also be described and analyzed by treating this
surface as one of the following two families of curves:
Fig. 4.3 An illustrative example of a θ-curve and a φ-curve on the manifold surface M.
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap04
It can be proven that all 3D and all planar arrays have zero geodesic
curvature for the φ-curves. This implies that
Furthermore, the geodesic curvature κg,θ for the θ-curves for 3D-grid and
2D-grid arrays is constant (depending on φo ) but is not zero. Hence
while constant κg,θ implies that the development of the θ-curves of grid
arrays are circles whose radius is a function of φo .
Using Eq. (3.42) it can be proven that for a general planar array the
geodesic curvature κg,θ is given as follows
9
det (G) 1
κg,θ = 3 (4.20)
gθθ tan φo
and this may not even be constant. The above results are summarized in
Table 4.3.
Table 4.3 Geodesic curvature.
Aφ|θo κg,φ 0 0 0
(hyperhelix)
For the 2D-grid and 3D-grid arrays, the off-diagonal elements of the
manifold metric G (i.e. gθφ , gφθ ) are zero (see Tables 4.1, 4.2) and hence
the θ-curves and φ-curves are orthogonal and constitute geodesic coordi-
nates i.e.
.
(θ-curves) ⊥ (φ-curves)
=⇒ geodesic coordinates (4.21)
(φ-curves) = geodesic curve
Figure 4.4 shows an illustrative representation of a planar array manifold
surface and its families of θ- and φ-curves as geodesic coordinates.
Theorem 4.5 The geodesic curves (φ-curves) on the manifold surface
are complex hyperhelices embedded in CN .
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap04
Fig. 4.4 Families of θ-curves and φ-curves on the manifold M of a planar array.
Hence hyperhelical curves on the manifold surface M are not only ana-
lytically “convenient” but also geodesic curves. Overall it may be concluded
that, although both families of curves can be used to describe the manifold
surface M, the θ-curves are the more complicated of the two families.
Corollary 4.1 All members of the family of φ-curves of planar arrays are
identical apart from their length lm and position in complex N-dimensional
space.
which implies that the two φ-curves at θ = θo and θ = θo + 180◦ have the
same length and curvatures. As illustrated in Fig. 4.5, these two φ-curves
can be considered as a continuation of one another, forming a composite
φ-curve having a hyperhelical shape of double the length of the original
curves. That is
where κg,θ (θ, φ) is the geodesic curvature of the θ-curves on the manifold
and x0 (φ) and y0 (φ) ∈ R are the integration constants such that:
=π =
x0 = − 12 0 κ−1 1 π −1
g,θ (ϑ, φ) cos(ϑ)dϑ and y0 = − 2 0 κg,θ (ϑ, φ) sin(ϑ)dϑ
The first condition of the theorem (i.e. KG = constant) is valid for 3D-grid
and all planar arrays, which implies that these arrays are “developable,”
with their “development” expressions summarized in Table 4.4.
π
θ
π cos φ f (ϑ) cos ϑdϑ − 1
f (ϑ) cos ϑdϑ
D (θ,φ) = 0θ 2 0π
det (L) 0 f (ϑ) sin ϑdϑ − 1
2 0 f (ϑ) sin ϑdϑ
c1 c2
Planar r 2 , rT r
x x y2
with L =
rT r ,
x y r y
c2 c3
3
and f (ϑ) c3 cos2 ϑ − 2c2 cos ϑ sin ϑ + c1 sin2 ϑ
sin θ
2D-grid D (θ,φ) = ρπ cos φ
− cos θ
Example 4.2 Figures 4.6, 4.7 and 4.8 provide the development of rep-
resentative 3D-grid, planar and 2D-grid array geometries. For instance the
development of a 3D-grid array with sensors located on the eight vertices
of a cube of side one half-wavelength and with reference point the centre
of the cube (array centroid) is shown in Fig. 4.6. For the planar array with
the positions of the sensors, in half-wavelengths, as given below:
−2.1 −1.1 0.4 0.9 1.9
r = −2.4 1.1 0.6 −0.4 1.1 (4.31)
0 0 0 0 0
the development is shown Fig. 4.7, while the development for a 5-sensor
uniform circular array (2D-grid array) of radius one half wavelength is pre-
sented in Fig. 4.8.
Fig. 4.6 θ-curves and φ-curves development of the manifold of a 3D-grid array with
sensors located on the eight vertices of a cube of side one half-wavelength and with
reference point the centre of the cube (array centroid).
have an ellipsoidal shape with the two main axes at angles ψ0 and ψ0 + π/2
with respect to the x-axis where:
1 2c2
ψ0 = tan−1 (4.32)
2 c1 − c3
! "
with ci denoting the elements of the matrix L, i.e. L = cc12 ,, cc23 . Indeed the
angle ψ0 can be found by examining the points of extremum curvatures.
Since the curvature of the development is equal to κgθ , the geodesic cur-
vature of the θ-curves of the planar array manifold is differentiated with
respect to θ and equated to zero:
dκgθ
=0 (4.33)
dθ
d(c3 cos2 ψ0 − 2c2 cos ψ0 sin ψ0 + c1 sin2 ψ0 )−3/2
⇔ =0
dψ0
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap04
Fig. 4.7 θ-curves and φ-curves development of the manifold of a planar array of 5-sensor
with positions given by Eq. (4.31).
⇔ or (4.34)
(c1 − c3 ) sin(2ψ0 ) − 2c2 cos(2ψ0 ) =0
the discriminant gives −4 det(L) < 0. Hence the first equation has no solu-
tion. The second equation gives
2c2 1 2c2
tan(2ψ0 ) = and hence ψ0 = tan−1 (4.35)
c1 − c3 2 c1 − c3
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap04
Fig. 4.8 θ-curves and φ-curves development of the manifold of a 5-sensor uniform cir-
cular array (2D-grid array) of radius one half-wavelength.
Fig. 4.9 Orientation of the “development” of θ-curves — axes at angles ψ0 and ψ0 +90◦ .
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap04
4.6 Summary
4.7 Appendices
& '
dgφφ dgθφ
2 gθθ dφ + gφφ dφ − 2gθφ dφ
1 dgθθ
6
×
6 ×(& det(G))−1 gθθ tan φ
'
− det(G) dφ tan φ + gθθ (1 + tan φ)
dgθθ 2
dgφφ dgθφ
gθθ dφ + gφφ dgdφθθ − 2gθφ dφ
=
2 det(G)gθθ tan φ
dgθθ
dφ tan φ + gθθ (1 + tan2 φ)
− 2 tan2 φ
(4.36)
gθθ
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap04
By using the Christoffel symbols of the second kind, Γθθφ and Γφφφ , which
can be expressed in terms of the elements of the manifold metric G as
& '
Γθθφ = 2 det(G)
1 dg
gφφ dgdφθθ − gθφ dφφφ
& ' (4.37)
Γφ = 1
g
dgφφ
+ g
dgφφ
− 2g
dgθφ
φφ 2 det(G) θθ dφ θφ dφ θφ dφ
KG is simplified to
−2gθθ tan φ + 2gθθ tan φ
KG = − 2 tan φ =0 (4.41)
gθθ
g = π 2 sin2 φ c1 cos
2
θ + 2c2 cos θ sin θ + c3 sin2 θ
(4.42)
φφ
gθφ = π cos φ sin φ c2 sin θ − cos θ + (c1 − c3 ) cos θ sin θ
2 2 2
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap04
Let the first term inside the braces be denoted by X and the second by
Y , i.e.
dD(ϑ, φo )
sd (θ, φo ) =
dϑ (4.44)
dϑ
o)
and hence
dD(θ,φ
dθ
= κ−1
gθ (θ, φo ) proving Eq. (4.29).
(iii) The curvature of a curve is the magnitude of the derivative of the
normalized tangent vector with respect to its arc length. Hence, con-
sidering the development of the manifold θ-curves for a constant φo ,
let u(θ) be its normalized tangent vector given by
dD(θ,φo ) 4 5
cos θ
u(θ, φo ) =
dθ
=±
dD(θ,φo )
(4.45)
dθ
sin θ
du (θ, φo )
du (θ, φo ) dθ
dsd
dθ dsd
4 5
− sin θ
=
κg (θ, φo )
cos θ
= κg (θ, φo ) (4.48)
Chapter 5
Non-Linear Arrays:
(α, β)-Parametrization
In the previous chapter the array manifold surface M was defined in terms
of azimuth θ and elevation φ angles and its intrinsic geometrical properties
were highlighted. Furthermore, the surface M was also considered based on
the families of θ-curves and φ-curves, whose characteristics were examined
using the differential geometry framework presented in Chapters 2 and 3.
In Chapter 4 it was observed that while the φ-curves were hyperhelical and
geodesic, the θ-curves possessed neither property. It was also pointed out
that since the θ-curves are not hyperhelical, their curvatures depend on θ
and so analytical evaluation of curvatures of order greater than two can
become exceedingly laborious and impractical.
In spite of the intuitive appeal of using azimuth and elevation angles
as the medium of array parametrization, this choice of angles is by no
means unique. In this chapter an alternative parametrization, known as
“cone-angle”, is presented. This parametrization is particularly suitable for
planar arrays resulting in two families of hyperhelical curves which can
provide a great deal of additional insight into the nature of planar array
behavior.
97
July 9, 2004 14:41 WSPC/Book Trim Size for 9in x 6in chap05
where
Ω = {(θ, φ) : θ ∈ [0◦ , 360◦ ) and φ ∈ [0◦ , 90◦ )}
with
a(θ, φ) = exp(−jrT k(θ, φ)) (5.2)
and is illustrated in Fig. 5.1. These angles α and β are called cone-angles
because the loci of wavenumber vectors of constant α = α0 (or equivalently
constant β = β0 ) form cones about the x̌ (equivalently the y̌) axis. In other
words, all possible incident signals from a fixed angle α, or β, form a cone
in space. Fig. 5.2 shows this for a random choice of cone-angles α0 and β0 .
The mapping
Fig. 5.3 Mapping from the (θ, φ) parameter space on to the (α, β) parameter space.
From the above equations it is obvious that, unlike in the θ-φ parametriza-
tion, α and β angles are not completely independent.
Furthermore, not all combinations of α ∈ [0, 180◦ ] and β ∈ [0, 180◦ ] are
acceptable/possible. The acceptable limits for α, when β = βo , and the
limits of β, when α = αo , are independent of the value of the rotation of
the x-y frame and may easily be derived from the limitations imposed by
Eqs. (5.6) and (5.7). These limits are given below:
if α = αo then
βmin βmax
90◦ − α β 90◦ + α for αo 90◦
o o
(5.11)
αo − 90◦ β 270◦ − αo for 90◦ < αo 180◦
βmin βmax
if β = βo then
αmin αmax
90◦ − βo α 90◦ + βo for βo 90◦
(5.12)
βo − 90◦ α 270◦ − βo
for 90◦ < βo 180◦
αmin αmax
angle Θ, then the sensor positions ř ∈ R3×N with respect to the Cartesian
frame (x̌-y̌-z) can be expressed as follows:
=r
cos Θ sin Θ 0
ř = cos(Θ + 90◦ ) sin(Θ + 90◦ ) 0 [rx , ry , rz ]T
0 0 1
with
r(Θ) rx cos Θ + ry sin Θ (5.14)
Based on the above, the wavenumber vector expressed in terms of the cone
angles α and β, and for a rotation angle Θ, can be written as
> 6 ?T
ǩ(α, β) = π cos α, cos β, 1 − cos2 α − cos2 β (5.16)
and thus, for isotropic sensors, the manifold vector of Eq. (1.22) can be
expressed as a function of α and β parameters as follows:
with a(α, β) given by Eq. (5.17). As might be expected, the intrinsic shape
of a manifold surface, as defined by its Gaussian curvature KG , is invariant
under any parameter transformation. For instance, using the general frame-
work presented in Chapter 3 and replacing p = α and q = β, the Gaussian
curvature of the manifold surface for a 3D grid array can be found to be
KG = (ρπ)−2 (5.19)
confirming (see Table 4.1) that the manifold is “developable” and that the
shape of the manifold is independent of parametrization.
Next we will focus on the differential geometry properties of the manifold
surfaces of planar arrays of omnidirectional sensors. Results for the 2D-grid
arrays can easily be deduced from the following results and those for 3D-grid
arrays can be analyzed with no difficulty.
Using Eq. (5.14) the manifold vector of Eq. (5.17) for a planar array is
simplified to
α β β
7 8
π 2 r(Θ)T r(Θ + 90◦ ) sin α sin β
2
π 2 r(Θ) sin2 α,
=
π 2 r(Θ)T r(Θ + 90◦ ) sin α sin β, π 2 r(Θ + 90◦ ) sin2 β
2
where, as usual,
∂a
ȧα = (5.22)
∂α
and
∂a
ȧβ = (5.23)
∂β
Whenever α or β equal 0 or π, the matrix G of the first fundamental coef-
ficients becomes singular. At such singularities, one or both of the tangent
vectors ȧα or ȧβ vanish and so the tangent plane is no longer defined. Fortu-
nately, the requirement φ = 0 (for the independence of α and β) established
July 9, 2004 14:41 WSPC/Book Trim Size for 9in x 6in chap05
at the outset (see Eq. (5.10)) avoids these singularities. Also note that the
cone parametrization avoids the singularity at the apex (φ = π2 ) which
exists with θ-φ parametrization and hence is more suitable for the study of
the manifold in that region.
The other variables required for the study of intrinsic geometry are the
Christoffel symbols. From the definitions of Chapter 3, it can be shown
that if
Note that the Christoffel matrices are quite sparse (even more so in the case
of balanced-symmetric arrays where gαβ = gβα = 0), which is a desirable
feature since it simplifies the differential geometry considerably.
We have seen that the manifold surface of any planar array (see
Table 4.2) has a zero Gaussian curvature. This is a property which is
independent of the parametrization. Indeed following any of the Eqs.
(3.31), (3.29) or (3.30), in conjunction with the Christoffel symbol matri-
ces defined above (see Eqs. (5.25) and (5.26)), the Gaussian curvature of
July 9, 2004 14:41 WSPC/Book Trim Size for 9in x 6in chap05
the array manifold can be found to be zero. Indeed, by using for instance
Eq. (3.30), we have:
∂
∂α Γα,ββ − ∂β Γα,βα
∂
1
Kgauss = +col2 {Γ2α }T col1 {Γ1β }
det[G]
−col2 {Γ2β }T col1 {Γ1α }
1 ∂gαβ
= cot β − 0 + 0 − gβα cot β cot α
det[G] ∂α
1
= (gαβ cot α cot β − gβα cot α cot β) (5.27)
det[G]
=0 (5.28)
Fig. 5.4 An example of an α-curve and a β-curve lying on the array manifold M.
Fig. 5.4 shows an illustrative representation of one “α-curve” and one “β-
curve” lying on the array manifold surface M.
The differential geometry of the α- and β-curves on the array manifold
surface are next investigated. Naturally, due to the symmetry of Eq. (5.20),
the parameter curves are expected to exhibit similar characteristics.
5.5.1 Geodecity
Following the definitions of Chapter 3, the geodesic curvatures of the α- and
β-curves (see Eqs. (3.42) and (3.43)) are given by
Γβ
κg,α = det[G] (g αα
)3/2
αα
(5.34)
κg,β Γα
= − det[G] (g ββ
)3/2
ββ
July 9, 2004 14:41 WSPC/Book Trim Size for 9in x 6in chap05
which is a very interesting result, implying that the two parameter curves
are of minimum length on the array manifold surface, and leads to the
following theorem.
Theorem 5.1 Both the α- and β-curves are geodesic curves for all val-
ues of Θ.
To check if, in addition to being geodesic curves, the α- and β-curves
are also “orthogonal” (constituting a geodesic coordinate system), the angle
between the tangent vectors ȧα and ȧβ should be found. Indeed
and so the tangent vectors ȧα and ȧβ are strictly orthogonal. For such arrays
matrix G is diagonal and the hyperhelical α-β curves form an orthogonal
coordinate system ∀Θ.
1• balanced-symmetric arrays: r x = r y and r x ⊥ ry
• unbalanced-symmetric arrays: r x = r y and r x ⊥ ry
July 9, 2004 14:41 WSPC/Book Trim Size for 9in x 6in chap05
Then, by integrating Eq. (5.42) over its limits αmin and αmax , the length
lm of an α-curve can be found as
Clearly the α- (or β-) curve degenerates into a single point (i.e. a curve
of zero length) whenever β (or α) approaches either 0 or π radians.
is apparent from Eq. (2.39) which indicates that curvatures are indepen-
dent of α and β respectively and are only functions of the array geometry
projected along the lines of azimuth at Θ and Θ + π2 . Consequently, the
curvatures and hence the shapes of the α, β-curves depend strongly on Θ,
reminiscent of the drastic variations in array performance which can be
observed along different lines of azimuth. Thus, Theorem 2.1 is the centre-
piece of the investigation of the α/β-curves and, as will be demonstrated
next, has some far-reaching implications.
From the above it is clear that all α-curves have the same ELA r(Θo )
but different “visible” areas specified by the minimum and maximum per-
missible values of α, i.e. Ωα|βo = [αmin , αmax ]. The same is true for the
β-curves but with r(Θo + 90◦ ) as their ELA with Ωβ|αo = [βmin , βmax ].
Fig. 5.5 An illustrative diagram of the length of the hyperhelical α-curve Aα|βo .
July 9, 2004 14:41 WSPC/Book Trim Size for 9in x 6in chap05
Corollary 5.2 For a given frame rotation angle Θo , all members of the
family of α- (or β-) curves of a planar array of N sensors are identi-
cal (congruent) in shape apart from their length and position in complex
N -dimensional space.
The above result means that the whole of the manifold surface may be
fully covered by the simple translation of either of the two individual curves
with displacement values given by v = r(Θo + 90◦ ) cos β0 for α-curves and
v = r(Θo ) cos α0 for β-curves. Note that α-curves are the same as β-curves
for rotation Θo + 90◦ .
Corollary 5.3 The α- and β-curves are identical to the θ-parameter
manifold of a linear array of isotropic sensors with locations r(Θ) and
r(Θ + π2 ) respectively.
Figure 5.6 illustrates the α-β parametrization of the array manifold in
the special case of Θ = 0. Notice that for Θ = 0, cone-angles are mea-
sured with respect to the x and y axes and consequently, according to
Corollary 5.3, the α- and β-curves are identical to the θ-parameter man-
ifolds of linear arrays with sensor locations r(0) = rx and r( π2 ) = ry
respectively.
Although the hyperhelical α (or β) family of curves shown above covers
the whole of the manifold surface and is sufficient to fully define the man-
ifold shape, it is not unique. In fact, to understand the behavior of the
array for all possible directions of arrival, it is necessary to consider the
hyperhelical coordinates corresponding to all values of Θ from 0 to π2 .
Fig. 5.6 (θ, φ) and (α, β) parametrizations of the same manifold surface M with the
frame rotation for the cone-angles Θo = 0.
July 9, 2004 14:41 WSPC/Book Trim Size for 9in x 6in chap05
Fig. 5.7 Families of α and β curves for a frame rotation Θo = 0◦ . and for Θo < 90◦ .
Corollary 5.5 The α-curve with β = 90◦ for a frame rotation Θo and
the composite φ curve Aφ|θo + Aφ|(θo +π) are one and the same curve if and
only if θo = Θo .
Similarly the β-curve with α = 90◦ for a frame rotation Θo is one and
the same curve with the composite φ curve Aφ|θo + Aφ|(θo +π) if and only if
θo = Θo + π2 .
The most significant implication of the property that the array manifold
M ∈CN of a planar array has a zero Gaussian curvature, is that M is a
July 9, 2004 14:41 WSPC/Book Trim Size for 9in x 6in chap05
• Length of A(d)α|βo
∈ R2 = Length of Aα|βo ∈ CN ,
• Principal curvature of A(d) α|βo
∈ R2 = Geodesic curvature of Aα|βo ∈ CN ,
• (A(d)
α|βo
, Aβ|α
(d)
o
) in R2 = (Aα|βo , Aβ|αo ) in CN .
Parameter 3D Planar
1
KG 0
(Eq. (3.31)) ρ2 π 2
cot β
κg,α 0
ρπ
cot α
κg,β 0
ρπ
opment
are
orthogonal straight lines (geodesics) and of length 2π
rx
sinβ
and 2π
ry
Chapter 6
Array Ambiguities
The parameter Ψ0 , shown in Fig. 6.1 and used in Eq. (6.3), is defined as
◦
0 for any linear array on the x-axis
θo for one φ-curve Aφ|θo
Ψ0 (6.4)
Θo
for all α-curves
Θo + 90◦ for all β-curves
113
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
p (generic) θ φ α β
The array manifold A of Eq. (6.1) is essentially a mapping from the param-
eter space Ωp ∈ R1 to the complex N -dimensional space CN
array geometry
p ∈ Ωp −→ a(p) ∈ CN (6.5)
However, in this study we will consider the problem which arises when
the mapping of Eq. (6.5) is not one-to-one. In this case the array cannot dis-
tinguish between two (or more) different signal environments. For instance
two different sets of signals impinging on the array can provide identical
responses at the array output, i.e. the same measurements. In such a case
any estimation algorithm may be unable to distinguish the true parameter
values (e.g. directions) from the false ones. This is termed as the ambiguity
problem.
from at least one of its subsets (or one subset from another).
array geometry
p1 −→ a(p1 ) ∈ CN
array geometry
(6.7)
p2 −→ a(p2 ) ∈ CN
with
a(p1 ) = ka(p2 ) (6.8)
and
k = scalar
Trivial ambiguities are the simplest type of ambiguity and can be easily
identified/detected. For instance it is impossible to distinguish whether a
signal is impinging on a linear array from direction-of-arrival p1 or from the
direction p2 = 360◦ − p1 , i.e. the mirror image of p1 with respect to 180◦ .
In this case
and for this reason, in order to avoid this ambiguity problem, the parameter
space in the case of linear arrays is confined to Ωp = [0◦ , 180◦ ).
However, Eq. (6.9) is a special case of Eq. (6.8) for k = 1. Eq. (6.8) indi-
cates that there exists at least one manifold vector a(p2 ) which is colinear
to another manifold vector a(p1 ). Then the two waves with bearings p1 , p2
are indistinguishable by the array even if p1 , p2 ∈ Ωp = [0◦ , 180◦ ).
Non-trivial ambiguities are a more complicated type of ambiguity arising
when a manifold vector can be written as a linear combination of two or
more different manifold vectors.
array geometry
pi ∈ Ωp −→ a(pi ) ∈ CN ; ∀i = 1, 2, . . . , M (6.10)
with
a(pM ) = c1 a(p1 ) + c2 a(p2 ) + · · · + cM −1 a(pM −1 ) (6.11)
Then this set (p1 , p2 , . . . , pM ) and all its subsets of M − 1 elements are
indistinguishable by the array.
In such a case the array will have identical responses for the sets of
bearings (p1 , p2 , . . ., pM −1 ), (p2 , p3 , . . ., pM ), (p1 , p3 , p4 , . . ., pM ), etc.
Fig. 6.2 illustrates a non-trivial ambiguous situation where the manifold
vector a(p3 ) belongs to the subspace L [A] where A = [a(p1 ), a(p2 )].
This type of ambiguity is much more difficult to identify since,
unlike a trivial ambiguity, it cannot be detected by a simple search
of the manifold. For example consider the non-uniform linear array
of four sensors with sensor locations (measured in half wavelengths)
given by
(0◦ , 60◦ , 90◦ ), (0◦ , 60◦ , 120◦ ), (60◦ , 90◦ , 120◦ ), (0◦ , 60◦ , 90◦ , 120◦ )
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Fig. 6.2 Illustrative representation of array ambiguity (the vectors a(p1 ), a(p2 ) and
a(p3 ) are linearly dependent).
For instance, for this array, the above four different signal environments
provide identical MUSIC spectra. Note that for this specific array,
from Fig. 6.3 where in the MUSIC spectrum, four nulls appear at directions
5◦ , 35◦ , 109.70◦ , 120.94◦ rather than two at 5◦ and 35◦ .
Note that in this study the array centroid will be taken as the reference
point, implying that sum(r)= rT 1N = 0. However, it can be proved that
the ambiguities are independent of the choice of the array reference point.
with si ∈ RN , ∀i.
is an ambiguous set of arc lengths with rank of ambiguity ρa , then any set
The proof of this theorem can be found in Appendix 6.7.1 while its
essential features are illustrated in the following example.
which is also an ambiguous set. This is also true for any ∆s so long as the
last element of š is smaller than, or equal to, the manifold length lm . Note
that the rotation should be carried out in the “arc length” domain and not
in the “directions-of-arrival” domain (parameter p domain). For instance,
if the directions of the original sources are rotated by ∆p = 20◦ so that the
new set of directions
(a) s1 = 0 and sj = 0 with 2 ≤ j ≤ c. That is, all the elements of the set
but the first element are non-zero.
(b) the rank of the N × c matrix A(s) with columns the manifold vectors
associated with the elements of the set is less than c, i.e.
(i.e. third condition is not satisfied). On the other hand, the subset
[s1 , s2 , s3 ]T is an ambiguous generator set since it satisfies all the three
conditions of Definition 6.5.
The next section is concerned with the identification of ambiguous gen-
erator sets (AGS) existing in the manifold of a linear array of arbitrary
geometry.
4 5T
lm lm lm
s∆rij = 0, ,2 , . . . , (c − 1) (6.23)
∆rij ∆rij ∆rij
where
1 + fix(∆rij ) if ∆rij ∈
/ N+
c=
∆rij if ∆rij ∈ N +
rs − πr + v))
a(s) = exp(j() (6.24)
in conjunction with
r
lm = 2π r and )
r= (6.25)
r
by AN
> & & ''
AN = exp (j (−πr + v)) , exp j ∆r 2π
ij
r − πr + v , . . . ,
& & '' ?
exp j (c − 1) ∆r
2π
ij
r − πr + v
)N
= diag(exp (j (−πr + v)) A (6.27)
where
4 5
) N = 1 , exp j2π 1 r , . . . , exp j(2π(N − 1) 1 r
A (6.28)
N
∆rij ∆rij
Consider now the ith and jth rows of A) N with i > j. By using the
property that
∆rij
exp −j2πk = 1, ∀k ∈ N + (6.29)
∆rij
) N ) = IN · rowT (A
rowTi (A )N )
i
1 0 0 ... 0
0 exp −j2π ∆rij
0 ... 0
∆rij
= .
∆rij
0 0 exp −j2π2 ∆rij . . . 0
0
0 0 ...
0
∆rij
0 0 0 . . . exp −j2π(N − 1) ∆rij
> & ' & '?T
1N , exp j2πri ∆r1ij , . . . , exp j2πri (N − 1) ∆r1ij
> & ' & '?T
r −∆r r −∆r
= 1N , exp j2π i∆rij ij , . . . , exp j2π(N − 1) i∆rij ij
> & ' & '?T
r r
= 1N , exp j2π ∆rjij , . . . , exp j2π(N − 1) ∆rjij
)N )
= rowTj (A
i.e.
) N ) = rowj (A
rowi (A )N )
& '
=⇒ det A ) N = 0 =⇒ det (AN ) = 0 (6.30)
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
It is easy to see that the fact that AN was chosen to consist of the first N
columns of A, and not any N columns, is not restrictive in the least since
the ith and the jth rows of any N × N submatrix of A will be equal. Since
all submatrices AN of A are singular, it results that A is rank deficient and
therefore, any subset of s∆ri,j with exactly N elements is an ambiguous set.
Note that
1 + fix(∆rij ) if ∆rij ∈
/ N+ lm
c= =⇒ c < 1 + ∆rij =⇒ (c − 1) < lm
∆rij if ∆rij ∈ N + ∆r ij
(6.31)
indicating that the last element of the vector s∆rij of Eq. (6.23), is smaller
that the manifold length lm . The vector s∆rij with elements c ∈ N + man-
ifold points will be known as the “uniform basic set” (UBS) of arc lengths
associated with the intersensor spacing ∆rij as it is constructed by uni-
formly partitioning the manifold length by the spacing ∆rij between the
ith and jth sensors.
According to Theorem 6.2, all subsets of N elements from a set s∆rij ,
of the form of Eq. (6.23) are ambiguous sets but not necessarily ambiguous
generator sets. In order for such a subset to be an ambiguous generator
set it must have its first element equal to zero and also satisfy the third
condition of the ambiguous generator set in Definition 6.5.
From the investigation of the precise cause of the ambiguous nature
of the set of arc lengths s∆rij constructed by uniformly partitioning the
manifold length by the spacing ∆rij , it is clear that the most common
AGSs are of rank N − 1. However, AGSs of rank less than (N − 1) can exist
if more pairs of rows of the determinant matrix A )N s
∆rij are identical.
This is possible if the array geometrical structure satisfies certain special
conditions and for this reason we need the following definition.
For instance, if ∆rij and ∆rjk are members of a distinct set with i =
j = k then ∆rik is not a member.
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Based on the above definition, the following theorem states that if there
are distinct sets of intersensor spacings then AGSs of rank less than (N − 1)
may be found.
Theorem 6.3 — Number of AGSs of a given rank: Consider an N -
sensor linear array. If k∆rij denotes the elements of a distrinct set of inter-
sensor spacings which are integer multiples of the smallest (∆rij ) amongst
them, and
∆rij N − k∆rij (6.32)
then the Uniform Basic Set s∆rij of c elements is an ambiguous set from
which
min(N −1,c−1)
c−1
(6.33)
=N −k∆rij
ambiguous generator sets of rank- N − k∆rij can be constructed.
The implication of the above theorem to the number of ambiguous gen-
erator sets of a linear array is twofold. Firstly, it reveals the existence of
additional ambiguous generator sets corresponding to intersensor spacings
which are less than (N − 1), provided that they occur in integer multiples of
each other, which includes the special case of repeated occurrences. In par-
ticular, if a spacing,
say ∆r
ij , occurs k∆rij times then the size of this spacing
need only be N − k∆rij to give rise to ambiguous generator sets. Secondly,
the larger the k∆rij is, the more rank- N − k∆rij ambiguous generator sets
can be constructed. This point can be illustrated from the intersensor spac-
T
ings ∆r13 = ∆r24 = 4.1 of the linear array r = [ −3 , −2.2 , 1.1 , 1.9 , 2.2 ]
where k∆r13 = 2 and ∆r13 (N − k∆r13 ) = 3. Thus from s∆r1,3 , five
rank-3 ambiguous generator sets can be constructed (instead of a single
rank-4 set).
From the array’s
geometrical
perspective, the following statement with
respect to rank- N − k∆rij ambiguity can be made:
Corollary 6.1 If k intersensor spacings of an N -sensor linear array
are integer multiples of the smallest among them, with the latter being at
least (N − k) half-wavelengths, then the array suffers from rank-(N − k)
ambiguity.
In summary the reason for the existence of these lower than N − 1 rank
ambiguities can easily be deduced from the definition of the ambiguous
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
generator set and especially its third condition (see Definition 6.5). More
analytically, it can be easily derived from the proof of Theorem 2.1, that
if two or more intersensor spacings of a linear array (or ELA) are integer
multiples of each other (which is the case for either of the above conditions)
and if s is an ambiguous generator set associated with the smallest of these
intersensor spacings, then the manifold matrix A(s) may have more than
two identical rows. This implies that its rank will be smaller than N − 1.
As a direct result of Theorem 6.3 and of the previous corollary it can
be said that the uniform linear array (or ELA) with intersensor spacing
∆r1,2 1 suffers from trivial ambiguities. Indeed, the intersensor spacing
∆r = ∆r1,2 occurs N −1 times. That is ∆r1,2 = ∆r2,3 = · · · = ∆rN −1,N , i.e.
k∆rij = N − 1 =⇒ ρa = N − N + 1 = 1 (6.34)
This well known result is rediscovered by observing that in such a case the
set s∆ri,j will contain at least two elements and the matrix A will be rank
deficient since all its rows will be equal.
By setting the value of ∆rij = |ri −rj | in the above equation to be equal
to the biggest possible difference between two array sensor locations, which
is obviously the aperture la = ∆r1,N , and by combining the above with the
Theorem 6.2(a), the following sufficient condition can be obtained.
Corollary 6.2 A SUFFICIENT condition for the presence of ambiguities
in ALL linear arrays of N sensors, is
since an array with sensor locations given by ř = br, with b > 1, might
have many intersensor spacings ∆rij which result in ambiguous sets of the
form given in Eq. (6.23).
So far it has been shown that Theorem 6.2 can be used in order to
identify ambiguous sets inherent in linear arrays of any geometry. When
this theorem is focused on specific array geometries it can produce some
more useful results as the following corollaries indicate.
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Corollary 6.3 All the ambiguous sets that exist in a two-element array
can be calculated from Eq. (6.23).
Proof. The proof is as follows. Consider the two-element array r =
[−r1 , r1 ]. Let s = [s1 , s2 ]T with s1 = 0 being an ambiguous set of arc
lengths. The matrix A with columns the manifold vectors corresponding
to s is
A = [a(0), a(s2 )]
s2
exp(jπr 1 ) exp jπr 1 1 −
= π r (6.36)
s2
exp(−jπr1 ) exp −jπr1 1 −
π r
Since it is assumed that s is an ambiguous set, it results that A is singular
and therefore
s2
det(A) = exp(jπr1 ). exp −jπr1 1 −
π r
s2
− exp(−jπr1 ). exp jπr1 1 − =0
πr
s2 s2
=⇒ exp jπr1 − exp −jπr1 =0
πr πr
s2
=⇒ sin r1 =0
r
s2
=⇒ r1 = kπ with k = 1, 2, 3, . . . (since s2 r1 = 0)
r
πr 2πr 2πr lm
=⇒ s2 = k =k =k =k
r1 2r1 ∆r1,2 ∆r1,2
That is
lm
s1 = 0 and s2 = k (6.37)
∆r1,2
with
k = 1, 2, 3, . . .
Therefore the only ambiguous set of arc lengths that can possibly exist for
a two-element array are of the form of Eq. (6.23). As a result, a necessary
and sufficient condition for the presence of ambiguities in a two element
linear array is
lm
< lm =⇒ ∆r12 > 1 (6.38)
∆r12
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
√ √ √
Corollary 6.4 The set of arc lengths s = [0, π 2, 2π 2, . . . , (c−1)π 2]T
is an ambiguous Uniform Basic Set for all three-element symmetric linear
arrays, as long as ∆r1,3 > 2.
Proof. For all three-element symmetric linear arrays with sensor pos-
itions given by r = [−r1 , 0, r1 ]T , the manifold length is equal to
6 √
lm = 2πr = 2π (−r1 )2 + (r1 )2 = 2πr1 2 (6.39)
The Uniform Basic Set s∆r1,3 associated with the difference between the
first and the third sensors is
> ?T
s = 0, ∆rlm1,3 , 2 ∆rlm1,3 , . . . , (c − 1) ∆rlm1,3
> √ √ √ ?T
= 0, 2πr1 2
2r1 , 2 2πr2r11 2 , . . . , (c − 1) 2πr2r112
√ √ √ T
= 0, π 2, 2π 2, . . . , (c − 1)π 2 (6.40)
Corollary 6.5 Let two arrays of N sensors have a common difference
∆r between two sensor locations. The Uniform Basic Set of arc lengths
associated with the common ∆r is different for these two arrays but the
corresponding ambiguous sets of DOA’s are the same.
Proof. Consider two different linear arrays (or ELAs) each having N
sensors with sensor locations r1 and r2 . Let there be a common dif-
ference ∆r between two sensor locations, which, for both arrays, satisfy
the conditions of Theorem 6.2. In this case, the two Uniform Basic Sets
will be:
First array:
> ?T
s1 = 0, 2πr 1 2πr 1 2πr 1
∆r , 2 ∆r , . . . , (c − 1) ∆r
(6.41)
Second array:
> ?T
s2 = 0, 2πr 2 2πr 2 2πr 2
∆r , 2 ∆r , . . . , (c − 1) ∆r
(6.42)
with c N , which implies that the above two sets are ambiguous. Using
Eq. (2.28) the corresponding set of DOA’s for the first array is
4 2πr 1
2πr 1
5T
θ1 = 0, arccos 1 − πr , . . . , (c − 1) arccos 1 − πr
∆r ∆r
1 1
T
= 0, arccos 1 − 2
∆r , . . . , (c − 1) arccos 1 − 2
∆r
(6.43)
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
T
= 0, arccos 1 − 2
∆r , . . . , (c − 1) arccos 1 − 2
∆r
(6.44)
c = 1 + fix(∆rij ) = 1 + 0 = 1 (6.45)
rank(A(s∆rij )) = c < N
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Form the “Table of Ambiguous Uniform Basic Sets” with rows all the
remaining vectors s∆rij . All these vectors are ambiguous sets and pro-
duce ambiguous generator sets. It is obvious that the set with the max-
imum number of elements is the set s∆r1,N .
(6) For each s∆rij (i.e. for each row of the “Table of Uniform Basic Sets”)
construct the AGSs (and associated rank of ambiguity), according to
the following rules:
Rule (a) If the non-zero elements of s∆rij cannot be found in any
other “Uniform Basic Set” vector, then
c−1
ambiguous generator sets
N −1
can be produced by this s∆rij . These ambiguous genera-
tor sets are all the possible subsets of N elements of s∆ri,j
with their first element zero and N − 1 non-zero elements
of s∆rij .
(3) Eliminating those entries of r r which are smaller than unity results in
∆r1,3 = 4.1
∆r = 4.9
1,4
∆r = 5.2
1,5
∆r = ∆r2,3 = 3.3 (6.46)
∆r2,4 = 4.1
∆r2,5 = 4.4
∆r3,5 = 1.1
(4) The seven Uniform Basic Sets (vectors) s∆ri,j are given in Table 6.2.
(5) The Uniform Basic Sets s∆rij which have less than N = 5 elements (i.e.
c < N ) are examined. These are the sets s∆r2,3 and s∆r3,5 .
• The number of elements in s∆r2,3 is 4 < N = 5 and, furthermore,
rank (A(s∆r2,3 )) = 4 = c and hence this row is eliminated.
• Also, s∆r3,5 is eliminated since c = 2 < N = 5 and
rank(A(s∆r2,3 )) = 2 = c. (6.47)
Thus the remaining sets are s∆r1,3 , s∆r1,4 , s∆r1,5 , s∆r2,4 and s∆r2,5 .
(6) In this step, the AGSs arising from the remaining Uniform Basic Sets
are calculated based on the rules (6a) and (6b).
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
s1 s2 s3 s4 s5 s6
sT
∆r1,3 0 7.4290 14.8580 22.2870 29.7160 —
sT
∆r1,4 0 6.2161 12.4322 18.6483 24.8644 —
sT
∆r1,5 0 5.8575 11.7150 17.5725 23.4300 29.2875
sT
∆r2,3 0 9.2300 18.4600 27.6900 — —
sT
∆r2,4 0 7.4290 14.8580 22.2870 29.7160 —
sT
∆r2,5 0 6.9225 13.8450 20.7675 27.6900 —
sT
∆r3,5 0 27.6900 — — — —
c ρa rule
sT
∆r1,3 5 3 (6b)
sT
∆r1,4 5 4 (6a)
sT
∆r1,5 6 4 (6a)
sT
∆r2,3 4 4 removed
sT
∆r2,4 5 3 (6b)
sT
∆r2,5 5 4 (6a)
sT
∆r3,5 2 2 removed
4
• s∆r2,5 : Finally, s∆r2,5 has c = 5 and it generates Nc−1
−1 = 4 = 1
ambiguous generator set. This AGS is the same as set s∆r2,5 and has
ρa = 4.
Based on Rule (6b)
However, any one of the four subsets š with four elements of s∆r1,3
involving one zero and three non-zero elements, is also an ambiguous
set. The matrix with columns the corresponding manifold vectors has
rank(A(š)) = 3 < c = 5. Therefore, the set s∆r1,3 satisfies the third
condition of the definition of the ambiguous generator set. Similarly,
the four subsets š are ambiguous generator sets since all their sub-
sets with three elements are unambiguous. Hence, the following five
ambiguous generator sets can be defined from s∆r1,3 , all of which
have ρa = 3. Indeed
4
4 4 4
= = + = 5 (6.49)
3 4
=3
• s∆r2,4 : Since s∆r2,4 is the same as s∆r1,3 , it follows that the ambigu-
ous generator sets obtained from this row are the ones which were
already obtained.
(7) The ambiguous generator table, Table 6.4, with rows the AGSs pro-
duced in Step-6 is now constructed. Their associated rank of ambiguity
is also shown.
AGS UBS
s1 s2 s3 s4 s5 ρa
The material presented in this chapter is valid for any hyperhelical curve
(of the form of Eq. (6.2)) embedded in CN . However, for p-curves Ap|qo
(p = φ, α or β) some extra steps should be taken to accommodate the
property that these curves of hyperhelical shape are lying on the manifold
surface of a planar array. These curves can be seen as the array manifolds
of an ELA with location given by r(Ψ0 ) represented by Eqs. (6.3) and (6.4).
Therefore, the Uniform Basic Sets s∆rij (see Eq. (6.23)), and consequently
the AGSs, can be constructed according to the procedure described in the
previous section with
where ∆rij represents the intersensor spacing between the ith and jth sen-
sors of the ELA r(Ψ0 ).
Before we proceed it is important to point out that
implies that the two φ-curves Aφ|θo and Aφ|θo +180◦ have the same length
and curvatures and can be considered as a continuation of one another
forming a composite φ-curve having a hyperhelical shape. The family of
these composite curves (from now on known as ϕ-curves) can be described
with an alternative but equivalent parameter space as
is essential before using the technique described in Section 6.3.1 and is given
as follows:
ϕ=φ
if θo < 180◦ then (6.60)
ϑ o = θo
ϕ = 180◦ − φ
if θo 180◦ then (6.61)
ϑo = θo − 180◦
AGS UBS
1st s∆r1,6
2nd s∆r1,5 ;s∆r2,6
3rd s∆r1,4 ;s∆r3,6
4th s∆r2,4 ;s∆r3,5
Note that the directions in degrees associated with the above ambiguous
generator sets of arc lengths can be estimated using Eq. (2.28) of Chapter 2
and are given in Table 6.6.
s1 s2 s3 s4 s5 s6 ρa
Table 6.6 Ambiguous generator sets (in degrees) corresponding to Table 6.5.
(1) (5◦ , 0◦ ) (5◦ , 51.79◦ ) (5◦ , 76.29◦ ) (185◦ , 81.69◦ ) (185◦ , 58.26◦ ) (185◦ , 24.84◦ ) 5
(2) (5◦ , 0◦ ) (5◦ , 53.24◦ ) (5◦ , 78.64◦ ) (185◦ , 78.20◦ ) (185◦ , 52.69◦ ) — 4
(3) (5◦ , 0◦ ) (5◦ , 55.49◦ ) (5◦ , 82.36◦ ) (185◦ , 72.51◦ ) (185◦ , 42.78◦ ) — 4
(4) (5◦ , 0◦ ) (5◦ , 57.29◦ ) (5◦ , 85.36◦ ) (185◦ , 67.75◦ ) (185◦ , 33.05◦ ) — 4
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
with αmin and αmax given by Eq. (5.12). As we have seen in Chapter 5
this family covers and represents the whole of the manifold surface M of a
planar array while all its members are associated with the same ELA
but with different “visible” areas specified by the minimum and maximum
permissible values of α, i.e. the parameter space is Ωa = [αmin , αmax ].
The same is true for the β-curves but with r(Θo + 90◦ ) as their ELA with
parameter space Ωβ = [βmin , βmax ].
A few other points must be emphasized here. Firstly, from Eq. (6.3),
the following relation can be derived:
which implies that for a given rotation Θo of the x–y frame and a given
βo , the resulting α-curve Aα|βo (Θo ) has the same length as the α-curve
Aα|βo (Θo + 180◦ ) resulting from rotation Θo + 180◦ . This is the reason
why only frame rotations Θo in the region [0◦ , 180◦ ] should be considered.
Secondly, as already mentioned, the roles of α and β can be interchanged
by simply replacing the rotation Θo of the x–y frame with Θo + 90◦ and
therefore the β-curves are effectively transposed versions of the α-curves.
This means that the α-curves for some rotation Θo , are the same as the
β-curves for rotation Θo +90◦ . This is the reason the β-curves are neglected
throughout this analysis and only the α-curves are considered.
Based on the above discussion the ambiguous generator sets of arc
lengths of the α-curve of βo = 90◦ , i.e. Aα|90◦ , can be estimated by par-
titioning this curve according to the elements (arc lengths) of Eq. (6.50)
using Ψ0 Θo .
Note that the α cone-angle parameter, corresponding to the arc lengths
sk of the α-curve, can be calculated by the following equation:
sk
αk = arccos sin βo − (6.72)
πr(Θo )
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
and that the α-curve with βo = 90◦ for a frame rotation Θo has identical
AGSs with the composite φ-curve associated with the composite azimuth
ϑo = Θo .
Indeed if βo = 90◦ then
Θo if α ≤ 90◦
θo = (6.73)
Θo + 180 if α > 90◦
◦
α if α ≤ 90◦
φ= (6.74)
180◦ − α if α > 90◦
Now, based on Eq. (6.50), let us assume that all AGSs of arc lengths of
the α-curve with βo = 90◦ , have been constructed. Then consider another
α-curve of βo different to 90◦ . For this second curve let smin and smax
be the arc lengths corresponding to the parameter angles αmin and αmax
(given by Eq. (5.12)), respectively. This curve has ambiguous sets but not
ambiguous generator sets, since smin = 0 for any β = 90◦ , which implies
that the first condition of the definition of the ambiguous generator set (see
Definition 6.5) is not satisfied. Let us define the ambiguous sets with their
first element equal to smin as the first permissible ambiguous sets. These can
be found by rotating all the ambiguous generator sets of the α-curve with
βo = 90◦ by smin , subject to the condition that the maximum element of
each set cannot exceed smax . Thus it is clear that the ambiguous generator
sets for the whole family of α-curves, for a given frame rotation Θo , can be
provided by examining only the α-curve of βo = 90◦ . The ambiguities of
any other α-curve can be generated by a simple rotation (Theorem 6.1) of
those ambiguous generator sets, making sure that its maximum element sc
is smaller than the manifold length lm of Aα|βo .
Based on the above discussion it is clear that the Ambiguous Generator
Table in arc lengths should be constructed for the α-curve with βo = 90◦ ,
using the technique of Section 6.3.1. Then, for a specific α-curve with
βo < 90◦ only the AGSs should be kept whose last element sc is smaller
than the length of the α-curve given by lmα|βo = 2πr(Θo ) sin βo . These
sets can be expressed from arc lengths to degrees using Eq. (6.72) but using
an additional term smin which takes into account the fact that for βo = 90◦
the minimum value of α is not 0◦ . That is
smin
sk + πr(Θo )(1 − sin βo )
αk = arccos 1 −
(6.75)
πr(Θo )
s1 s2 s3 s4 ρa
(1) 0 19.8692 — — 1
(2) 0 6.6231 13.2462 — 2
(3) 0 4.9673 24.8365 — 2
(4) 0 9.9346 14.9019 24.8365 3
(5) 0 4.9673 9.9346 14.9019 3
AGS UBS
The sets of Table 6.7 are the ambiguous generator sets for the whole
family of α-curves (since ambiguities of any α-curve with βo < 90◦ will
relate to a subset of the above AGSs). To see this property let us consider
βo = 90◦ . In this case, all the AGSs of Table 6.7 are certainly included
since the length lm of this α-curve Aα|90◦ is equal to the manifold length
of the ELA lm = 25.8299. The values in Table 6.7, transformed in degrees,
are shown in Table 6.9. If a different α-curve of βo < 90◦ is chosen then its
length is
In this case only the rows in Table 6.7 with the last element less than lm
should be kept. Indeed for the α-curve with βo = 74◦ , i.e. Aα|74◦ , then
lm = 24.8293 and therefore the AGSs will be as shown in Table 6.10.
The α-curve Aα|50◦ , having a length lm = 19.7869 will be associated
with only two out of the five AGS of the family. These are the ones with
their last element smaller than, or equal to, the length of the curve 19.7869
(i.e. s∆r1,3 and s∆r1,3 ) and are given in Table 6.11.
Finally, Table 6.12 gives the only AGS associated with the α-curve
Aα|35◦ (curve length lm = 14.8154), while no ambiguity can be identified
for any α-curve with βo < 30◦ or βo > 150◦ .
All these are shown in Fig. 6.4 which plots the length of the α-curves
versus the value of βo where the darker the shade of the area, the larger
the number of ambiguous sets that can be identified.
Note that in order to have a set of M (azimuth, elevation) directions
belonging to the same α-curve, there should be a frame rotation Θ = Θo
such that
where k is a constant. Then these directions which satisfy Eq. (6.79) belong
to the same α-curve Aα|βo , i.e correspond to a common βo = k for this
rotation Θo . For instance this can be seen for the following planar array:
T
−1, 2, −4, 3
rx , ry , 0 = 0, 0.5, 1.5, −2 (6.80)
0, 0, 0, 0
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Fig. 6.4 Manifold length for the family of α-curves. In the region of the darker shade
the number of AGS is 5 while no AGS can be identifed in the white region (0 AGSs).
Indeed in this case Eq. (6.79) is satisfied for Θ = Θo = 39◦ with constant
k = 82◦ (= βo ) as can be seen from Fig. 6.5.
This indicates that these directions are on the same α-curve which is
the manifold of the ELA r(Ψ0 ) with Ψ0 = Θo , i.e.
This α-curve (manifold) can be found using Eqs. (5.3) and (5.4) to trans-
form the set of (θ, φ) to cone-angles for the frame rotation Θo . Thus for
Θo = 39◦ the directions (θ, φ) are transformed to cone-angles (α, β) given
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Fig. 6.5 Plot of the four equations given by Eq. (6.79) for different rotation Θ.
This implies that the common βo is 82◦ for this rotation Θo (i.e. belong to
the same α-curve Aα|βo ). This curve is the manifold of the following ELA
r(Ψ0 ) with Ψ0 = Θo :
T
r(39◦ ) = rx cos 39◦ + ry sin 39◦ = −2.1646, −0.7771, 1.0728, 1.8690
(6.83)
Note that if there exists an intersensor difference ∆rij = |ri (Θo ) − rj (Θo )|
for which all elements of α satisfy
2
α = arccos sin β0 1M − [1, 2, . . . , M ] T
(6.84)
∆rij
array. Based on similar arguments to the linear array case the following
statements about ambiguities inherent in planar arrays can be made:
Proposition 6.1
r = [rx , ry , 0N ]
= [r1 , r2 , . . . , rN ]T (6.89)
Geometrically, this is equivalent to the line joining the pth and qth sensors
and that joining the ith and jth sensors, being parallel, with the length of
the former being an integer multiple of the latter. A special case is when
∆rpq = n∆rqj ; n ∈ Z + , resulting in
rp − rq = ± n rq − rj , n ∈ Z + (6.91)
In geometrical terms, this means that the three sensors are colinear, with
the spacing between the pth and qth sensors being an integer multiple of
that between the ith and jth sensors. Therefore, in general (k 2).
Proposition 6.2 A planar array will suffer from rank-(N − k) ambiguity
if either
(a) the spacings between k pairs of sensors are integer multiples of the
smallest among them, the latter being at least (N − k) half-wavelengths,
and the lines joining the sensor pairs are parallel,
or
(b) (k+1) sensor locations (vectors) are colinear and the k intersensor spac-
ings between adjacent sensor pairs are integer multiples of the smallest
among them, with the latter being at least (N − k) half-wavelengths.
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Note that the same factors which contribute to the growth of the number
of ambiguous generator sets in linear arrays are applicable to planar arrays
via the equivalent linear array.
Corollary 6.6 The causes of ambiguities in planar arrays are namely
the colinearity of the sensor locations, uniformity in the sensor spacings
and the size of the common spacing.
while an infinite number of ELAs are associated with the family of φ-curves
set of ELAs of the family of φ-curves = r(θo ) ∈ RN , ∀θo ∈ Ωθ
(6.93)
while its elements are shown as dots in Fig. 6.7. Furthermore, in the same
figure the locus of the manifold lengths of all ELAs, i.e. {lm (Ψ0 ), ∀Ψ0 } is
also shown.
By varying Ψ0 over its parameter space ΩΨ0 = [0◦ , 180◦ ) we will get
the set of ambiguous generator lines of the φ-curves to which the elements
of the ambiguous generator sets of arc lengths belong.
In these lines the AGS (in arc length) are associated with all ELA
r(Ψ0 ), ∀Ψ0 with Ψ0 = Θo or ϑo . This implies that, for variable Θo in
the region 0◦ to 180◦ , and for the array of Eq. (6.64) the ambiguous
generator lines of all families of α-curves are numerically, in arc lengths,
equal to those of the family of φ-curves. Thus for the array of Eq. (6.64),
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Fig. 6.7 The locus of the manifold lengths of ELA of Example 6.4. The dots represent
the 1st AGS (rank 5) given in Table 6.5.
Fig. 6.8 The set of ambiguous generator lines of rank 5 in which the elements of the
1st AGS of Table 6.5 (Example 6.4) belong.
Figs. 6.8–6.11 not only represent the set of ambiguous generator lines of
the φ-curves but also of the α and β-curves by generalizing the polar axes
to (Ψ0 , s).
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Fig. 6.9 The set of ambiguous generator lines of rank 4 in which the elements of the
2nd AGS of Table 6.5 (Example 6.4) belong.
Fig. 6.10 The set of ambiguous generator lines of rank 4 in which the elements of the
3rd AGS of Table 6.5 (Example 6.4) belong.
Thus, Fig. 6.8 shows the set of ambiguous generator lines of rank 5 in
which the first row of Table 6.5 belongs. As it can be seen in this figure,
the set of ambiguous generator lines ceases to exist at θo = 0 and for values
of θo in the region 61◦ to 119◦ while for θo in the region 119◦ to 180◦ its
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
Fig. 6.11 The set of ambiguous generator lines of rank 4 in which the elements of the
4th AGS of Table 6.5 (Example 6.4) belong.
values are the mirror image of the values from 0◦ to 61◦ . Note that the set
of ambiguous generator lines is defined only in those areas of the parameter
space Ωθ, at which the last line (e.g. the 6th line in Fig. 6.8) in which the
largest element of the associated AGS is located, is below the locus of the
manifold length of all ELAs.
Fig. 6.12 “Discrete” AGSs of corresponding to ELA of 45◦ , 90◦ and 135◦ for the array
of the Example 6.4.
Fig. 6.13 All the sets of ambiguous generator lines of the array of Example 6.4 shown
together with the “discrete” AGSs.
In the case of a symmetric linear array, the normalized sensor locations are
simply given by j times the eigenvalues of the corresponding Cartan matrix
(since the array and its mirror image about the centroid are identical).
Therefore, in the symmetric linear array case, the denominator ∆r̃1,N in
Eq. (6.97) is equal to the eigenvalue spread of the Cartan matrix-spread(C).
Note that the spread of a matrix A with eigenvalues denoted by λ, is
defined as:
spread (A) max |λi − λj | (6.98)
i,j
Recall that the Cartan matrix is a real normal matrix (CCT = CT C), which
in the case of a symmetric array, is given by
0 −κ1 0 ... 0
κ1 0 −κ2 ... 0
0 κ2 0 ... 0
C= (6.101)
. .. .. .. ..
.. . . . .
0 0 0 κN −1 0
July 7, 2004 14:4 WSPC/Book Trim Size for 9in x 6in chap06
(N − 1) π
lm (6.104)
κ1
then, rank-(N −1) ambiguity, resulting from uniform partitions of the mani-
fold, is bound to exist. If instead,
√
lm < 2 (N − 1) π (6.105)
It is however important to note that the condition set out in Eq. (6.105)
does not guarantee the absence of other types of rank-(N − 1) ambiguity,
for instance, those corresponding to non-uniform partitions of the array
manifold.
6.7 Appendices
s = [s1 , s2 , . . . , sc ]T (6.106)
=0 (6.111)
and is hence singular. This implies that Ǎ is rank deficient and therefore
the set š is also an ambiguous set. Furthermore the ranks of ambiguity of
s and of š are the same. This is because the submatrices of A(s) with less
than c rows are non-singular which implies that the submatrices of Ǎ with
less than c rows are also non-singular.
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap07
Chapter 7
More on Ambiguities:
Symmetrical Arrays
1 Every sensor has a symmetric counterpart with respect to the the array centroid.
157
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap07
The matrix AN (s) is a square matrix and, therefore, its determinant can
be used to examine the linear dependence of its columns. An expression for
the determinant of AN (s) can be provided as follows:
det(AN (s)) = det(U(0)) · det([ F(s1 )R, F(s2 )R, . . . , F(sN )R]) (7.6)
Noting that for a symmetric linear array, the columns of U(s) are mutually
orthogonal, that is
det(AN (s)) = ±det[ F(s1 )R, F(s2 )R, . . . , F(sN )R] (7.8)
Because of the special structures of the radii vector R and of the Frame
matrix F, the det(A(s)) can be further written as
RT F(s1 )T
RT F(s2 )T
det(AN (s)) = ± det
...
R F(sN )
T T
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap07
RT F(−s1 )
RT F(−s )
2
= ±det
...
RT F(−sN )
RT F(s1 )
RT F(s )
2
= ±det (7.9)
...
RT F(sN )
From Eq. (2.61) evaluated at s = 0 and for a symmetric linear array:
RT = aH (0)U(0) (7.10)
and therefore,
a(0)H U(0)F(s1 )
a(0)H U(0)F(s )
2
det(AN (s)) = ±det
..
.
a(0)H U(0)F(sN )
a(0)H U(s1 )
a(0)H U(s )
2
= ±det
..
(using Eq. (7.5))
.
a(0) U(sN )
H
a(0)H u1 (s1 ) a(0)H u2 (s1 ) ... a(0)H uN (s1 )
a(0)H u1 (s2 ) a(0)H u2 (s2 ) ... a(0)H uN (s2 )
= ±det
.. .. .. ..
. . . .
H
a(0) u1 (sN ) a(0)H u2 (sN ) . . . a(0)H uN (sN )
(7.11)
0 ≤ s ≤ lm = 2πr (7.12)
where u1 (s) denotes the tangent vector at the manifold point s, i.e.
u1 (s) = a (s).
This is the starting point for the identification of a new class of ambiguous
generator set having a rank of ambiguity ρa equal to (N − 1).
The existence, however, of at least N characteristic points within the
manifold length lm is very important for the identification of ambiguous sit-
uations and this, in turn, relates to the number of windings. In particular we
have seen in Chapter 2 (Theorem 2.2, page 41) that the manifold of a linear
array of N sensor has a hyperhelical shape with the number of windings
(N = odd) or half-windings (N = even) nw given by the following equation
lm
nw = (7.15)
lw
where lw is the length of one winding or one half-winding and is given by
the (N − 1)th positive root of the following expression:
This expression and its derivatives with respect to s are the pivots for
identifying ambiguities which are only present in symmetric linear arrays
(or symmetric ELAs).
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap07
Equation (7.15) indicates that when nw > 1, then there is more than
one winding (or half-winding) contained in the manifold length of the curve.
It must be noted that uniform linear arrays with intersensor spacing equal
to one half-wavelength have lm = lw =⇒ nw = 1. A spiral of a spring could
be used to provide a visualization of the concept of a winding.
lw 2
nw = (7.17)
2πr | cos p1 − cos pN |
Figure 7.1 is the plot of Eq. (7.16) with respect to the arc length. From
this plot it can be seen that the manifold curve consists of more than one
half-winding. Actually, the length of one half-winding is equal to lw =
14.0476, i.e. the arc length at the (N − 1)th positive characteristic point.
The number of half-windings of the hyperhelical curve is also calculated
using Eq. (7.15) and is equal to nw = 1.2603. The first four roots (arc
lengths) of Eq. (7.16) are
Fig. 7.1 Plot of Tr(C expm(sC)) with respect to s for the array of Eq. (7.18).
Furthermore the set s has its first element equal to zero and also satisfies
all the other conditions of the definition of the ambiguous generator set
(see Definition 6.5, page 120). Hence s is also an ambiguous generator set
of parameters.
AGSs and this complicates the representation and identification of this type
of ambiguity.
where the expression F(s) = Cexpm(sC) has been used. Then by differen-
tiating Eq. (7.23) we get
(n−1)
a(0)H u1 (s) = Tr(Cn expm(sC)), n = 1, 2, . . . (7.24)
(n)
where u1 denotes the nth derivative of the tangent vector at the manifold
point s. This leads to the following definition:
0 ≤ s ≤ lm = 2πr (7.25)
the following table gives the AGSs produced from the Uniform Basic Sets
s1 s2 s3 s4 ρa
0 6.7672 13.5344 — 2
0 4.8337 9.6674 14.5011 3
0 4.8337 9.6674 19.3349 3
0 4.8337 14.5011 19.3349 3
0 9.6674 14.5011 19.3349 3
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap07
and since it is a symmetric array, the roots of Eq. (7.29), (for n = 1) i.e.
the set of characteristic points,
s(1)
cp = [0, 5.5147, 9.3228, 14.0360, 19.8000]
T
s1 s2 s3 s4 ρa
However for n = 1.9 the first four roots of Eq. (7.29) for s ∈ [0, lm ] forming
the vector
s(1.9)
cp = [2.9041, 7.6098, 12.0167, 17.1851]T (7.32)
is also an ambiguous set and of rank N − 1 = 3. Eq. (7.29) has been plotted
in Fig. 7.2 for n = 1.9.
From Example 7.2 it can be observed that the first root of Eq. (7.29) is
not equal to zero for all values of n. This means that although in the case
of n = 1 the AGSs can be simply identified as all the combinations of three,
(1)
or in general N − 1, non-zero elements of scp with a zero prepended, the
same cannot be done for n = 1.9. In order to find the AGSs in this case, the
first element of the ambiguous set must be subtracted from all its elements.
Note that the first root of Eq. (7.29) is equal to zero only if n = 1, 3, 5,
etc. This can be seen from Eq. (7.29) which for s = 0 becomes equal to
Tr(Cn ) that is zero only if n is odd, since in that case its structure is
identical with the structure of C (zero diagonal, more precisely only its
two secondary diagonals are non-zero). This result is confirmed by Fig. 7.3
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap07
Fig. 7.3 Plot of smallest root of Tr(Cn expm(sC)) versus n for the array in Example 7.2.
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap07
where the first root of Tr(Cn expm(sC)) is plotted for n = 0.1 to 10 in steps
of 0.1 for the Cartan matrix of the array in Example 7.2.
From the previous discussion it is clear that an infinite number of AGSs
exist in symmetric linear arrays because n takes infinite values. However
for large values of n there exists a relationship between Theorems 6.2 and
7.1. This is stated in the following lemma:
Lemma 7.1 For n → ∞ (large), the roots of Tr(Cn expm(sC)) = 0,
n > 0 are a shifted version of the set of arc lengths created by partitioning
the manifold in segments equal to lm /la , where la = ∆r1,N is the aperture
of the array.
Example 7.3 Consider the array
Using the aperture intersensor spacing ∆r1,5 = 5.3740 the UBS s∆r1,5 can
be found as
s(13)
cp = [0, 4.9162, 9.8322, 14.7483, 19.6645, 24.5807]
T
(7.35)
cp s∆r1,5
s(13) (7.36)
s(1)
cp = [0, 5.9142, 10.5499, 14.2941] (7.39)
Furthermore, as can be seen from Fig. 7.5, the gain pattern of the array
steered towards endfire exhibits stationary points at
Fig. 7.4 The function Tr(Cexpm(sC)) and the array pattern, plotted as a function of
the arc length, for a symmetric array with locations r = [−1.45, −0.8, 0.8, 1.45 ]T . Note
that the arc length space s ∈ [0, lm = 14.7153) is equivalent to the azimuth space
θ ∈ [0◦ , 180◦ ), while the black dots constitute a rank-3 AGS of arc lengths.
is also ambiguous.
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap07
Corollary 7.2 All equivalent linear arrays r(Ψ0 ),∀Ψ0 (see Eq. 6.4) asso-
ciated with the families of p-curves (p = ϕ, α, β) of a symmetric planar
array are symmetric and therefore can suffer from both uniform and non-
uniform types of ambiguity.
with two collocated sensors at −1.45. Discarding one of the collocated sen-
sors, the resulting array is found to be symmetric and given by
7.6 Conclusions
Although all “uniform” ambiguities can be found using the Uniform Basic
Sets of arc lengths of Eq. (6.23), there is no single coherent framework
July 6, 2004 9:30 WSPC/Book Trim Size for 9in x 6in chap07
Chapter 8
Array Bounds
174
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
with the vector u1 (s̆) as its tangent, and has radius R which is equal to the
inverse of the curvature of the curve at s̆. This implies that
∆si ∆si
u1 (s̆ ± ∆si ) = cos u1 (s̆) ± sin u2 (s̆) (8.2)
R R
∆s2i ∆si
1−
2 R
2R
A = {a(p) ∈ CN , ∀p : p ∈ Ωp } (8.3)
Combining Eqs. (8.4) and (8.5), and using the special skew symmetric form
of the Cartan matrix (Eq. (2.14)), the following approximation to the mani-
fold tangent vector at s̆ ± ∆si (i.e. the first column of U(s̆ ± ∆si )) may be
obtained:
(∆si )2 2 (∆si )2
u1 (s̆ ± ∆si ) 1 − κ1 u1 (s̆) ± ∆si κ1 u2 (s̆) + κ1 κ2 u3 (s̆)
2 2
(8.6)
However, the 3rd term of Eq. (8.6) can be ignored so long as a tolerance
factor T OL 1 is defined such that:
@ @
@ 3rd term @ 1
@ @
@ 2nd term @ = 2 ∆si κ2 ≤ T OL (8.7)
2κ1
∆si ≤ 3 T OL (8.9)
r − κ21 r
(∆si )2 2
u1 (s̆ ± ∆si ) 1− κ1 u1 (s̆) ± ∆si κ1 u2 (s̆) (8.10)
2
Fig. 8.1 H12 subspace used to approximate the array manifold at a local level.
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
∆si = κ̂−1
1 ∆ψi (8.16)
Fig. 8.2 Two alternative ways of visualizing the ‘circular approximation’ of array man-
ifold curve at an arbitrary point s̆: (a) view showing H12 plane, (b) view orthogonal to
L ([H12 , a(s̆)]).
or equivalently, as
µyi
sin(∆ψi ) = (8.19)
κ̂−1
1
Hence, substituting ∆ψi given by Eq. (8.19) back into Eq. (8.16) we have
∆si = κ̂−1
1 ∆ψi
arccos(1 − κ̂1 µxi )
= (8.21)
κ̂1
arcsin(κ̂1 µyi )
(or, equivalently) =
κ̂1
It must be noted that if the inclination angle ζinc (see Eq. (2.55)) of the
manifold is 0◦ then ζ = 90◦ and κ̂1 = κ1 , which is an important feature of
symmetric linear arrays.
σ2 & '−1
CRB[p] = Re(H RTm ) ∈ RM ×M (8.22)
2L
where
H
H = Ȧ P⊥
A Ȧ ∈ C
M ×M
(8.23)
RTm = E{m(t)m(t)H } ∈ CM ×M = source covariance matrix
with
A = [a1 , a2 , . . . , aM ]
Ȧ = [ȧ1 , ȧ2 , . . . , ȧM ] (8.24)
⊥ H −1
PA = IN − A(A A) AH
E{n(t) nH (t)} = σ 2 IN ,
In this section the nature of the CRB and its relationship with the
differential geometry of the array manifold (and hence the array geometry)
is clarified by focusing on the special cases of one emitter as well as two
emitters located close together.
P = E{m(t)m∗ (t)}
from bearing p. Then, since Rm = P and A = a(p), Eq. (8.22) implies that
the CRB may be expressed as
σ2 1
CRB[p] = (8.25)
2LP ȧH (p)P⊥a ȧ(p)
Recalling that ȧ(p) = u1 (s) ṡ(p) and u1 (s) = 1 the above equation
becomes
σ2 1
CRB[p] = (8.26)
2LP ṡ(p) u1 (s)P⊥
2 H
a u1 (s)
σ2 1
CRB[p] = =
2LP ṡ(p) 2 2(SNR × L)ṡ(p)2 (8.27)
2
where SNR = P/σ
Note that the above expressions are valid for any manifold curve, not
necessarily of hyperhelical shape. For instance Eq. (8.27) can be written so
as to provide information about the accuracy of p estimates in terms of the
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
Table 8.1 Rate-of-change of arc length for linear and planar arrays. The
planar arrays are parametrized in terms of (θ, φ) and (α, β).
CRB[p|Ap|qo ] (8.28)
is adopted to represent the CRB for the p-estimates and show its depen-
dence on the properties of the Ap|qo curve, with the subscript ‘o’ indicating
that q is assumed constant and known. Similarly, the q-estimates are asso-
ciated with the curve Aq|p on M as the parameter p is assumed constant
and known, with the CRB of q-estimates represented by
CRB[q|Aq|po ] (8.29)
Using this notation in conjunction with Eq. (8.27) and Table 8.1 for an
emitter at a direction (θ, φ), the CRB for the φ-estimates and θ-estimates
are expressed as follows
1
CRB[φ|Aφ|θo ] =
2 (SNR × L) (π r(θo ) sin φ)
2
(8.30)
1
CRB[θ|A θ|φo ] = 2
2 (SNR × L) (π ṙ(θ) cos φo )
dr(θ)
where as usual r(θ) = rx cos θ + ry sin θ and ṙ(θ) = dθ .
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
therefore:
1
CRB[φ|Aφ|θo ] =
2
2L × SNR (π rx sin(φ))
⇒ φ accuracy is independent of θ
(8.34)
1
CRB[θ|Aθ|φo ] =
2L × SNR (π rx cos (φo ))
2
⇒ θ accuracy is independent of θ
σ2 1
CRB[s] = = (8.35)
2LP 2(SNR × L)
In other words,
1 1 1
CRB[p1 ] = (8.36)
2(SNR1 × L) ȧ1 P⊥
H
A ȧ1 1 −
Re2 [ρȧH ⊥
2 PA ȧ1 ]
⊥ ȧ )(ȧH P⊥ ȧ )
P1 P2 (ȧH P
1 A 1 2 A 2
matrix with columns the M manifold vectors. To make Eq. (8.36) more
tractable, consider the scenario where all the M emitters are uncorrelated
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
The dependence of the CRB on the rate of change of manifold arc length,
ṡ(p), become apparent by considering a1 and a2 as two points on the
same manifold curve A. Then using the expression ȧ(p1 ) = u1 (s1 ) ṡ(p1 ),
Eq. (8.37) can be written as follows:
1 1
CRB[p1 |A] = (8.38)
2(SNR1 × L) ṡ(p1 ) u1 (s1 )P⊥
2 H
A u1 (s1 )
1 2
CRB[p1 |A] = (8.39)
(SNR1 × L) ṡ(p1 ) (∆s) κ̂21 (p̆) −
2 2 1
N
∆s = ṡ(p̆) ∆p and
where (8.40)
κ̂1 (p) = κ1 (p) sin(ζ(p))
Note that the bearing p̆ corresponds to the point with arc length s̆, and
to a first-order approximation, equals (p1 + p2 )/2 or equivalently p̆ =
p1 + ∆p2 .
It is important to point out that, since ∆s = ṡ(p̆)∆p, Eq. (8.39) indicates
that the CRB is inversely proportional to the square of the emitters’ angular
separation,
1
CRB[p1 |A] ∝ 2 (8.41)
(∆p)
i.e. the accuracy of the system deteriorates rapidly as the emitters approach
one another.
The technique deployed for the derivation of Eq. (8.39) can also be
directly applied to scenarios involving more than two closely spaced emit-
ters. However, due to the proliferating number of manifold vector inner
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
Table 8.2 Principal curvature for linear and planar arrays parametrized —
planar arrays are in terms of (θ, φ) and (α, β).
= R and Ṙ = dR
where R
R dp
products, the expression for the CRB can become very cumbersome.
Nevertheless, it can be shown that so long as the additional manifold points
s3 , s4 , . . . , sM are not in the neighborhood of s̆ (i.e. in the neighborhood of
s1 and s2 ), the value of CRB[p1 |A] is primarily dominated by the presence
of the source at bearing p2 .
Example 8.2 In Fig. 8.4 the value of CRBθ [θ1 ], as defined in [24], is
compared to that given by Eq. (8.39) for two unit-power emitters located
at bearings θ1 = 20◦ and θ2 = θ1 + ∆θ, and both at a common elevation of
φ = 20◦ . The expressions are evaluated for the 24-element Y-shaped array,
assuming the availability of L = 100 snapshots and a signal-to-noise ratio
of 10 dB (for each emitter).
As can be seen, for small emitter separations, Eq. (8.39) provides an
excellent estimate of the CRB. However for increasing separations the
circular approximation breaks down and Eq. (8.39) is no longer valid. Note
that for large emitter separations, the exact CRB settles down to the value
corresponding to a single emitter.
Fig. 8.4 Plots of the CRB Eqs. (8.37) (exact) and (8.39) (appr.) for two emitters with
angular separation ∆θ (Y-shaped array of 24 sensors). The CRB for a single source is
also shown for comparison.
where the subscripts θ and φ have been added for reference to θ-curves and
φ-curves respectively.
Also since the φ-curves are hyperhelical, κ̂1,φ (φ) is independent of φ.
Furthermore, for large arrays, κ̂1,θ (θ) is independent of φ so long as elevation
is not close to 90◦ . The combination of these results, along with Eq. (8.39)
reveal the following dependencies on elevation φ:
• For emitters corresponding to two points on the same θ-curve Aθ|φo , that
is having the same elevation angle, i.e. (θ1 ,φo ) and (θ2 , φo ),
1
CRB[θ1 |Aθ|φo ] ∝ (8.43)
ṡ2θ (θ1 )ṡ2θ (θ̆)
1
∝ (8.44)
cos4 φo
θ1 + θ 2
where θ̆ = (8.45)
2
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
• For emitters with the same azimuth angle, i.e. (θo ,φ1 ) and (θo , φ2 ), hence
on the same φ-curve Aφ|θo
1
CRB[φ1 |Aφ|θo ] ∝ (8.46)
ṡ2φ (φ1 )ṡ2φ (φ̆)
1
∝ (8.47)
sin φ1 sin2 φ̆
2
φ1 + φ2
where φ̆ = (8.48)
2
Consequently, the variations of CRB[θ1 |Aθ|φo ] and CRB[φ1 |Aφ|θo ] with
respect to φ are independent of the array geometry and 90◦ out of phase.
Incorporating these results into Eq. (8.39), and on the grounds that for
closely spaced emitters
ṡθ (θ1 ) ṡθ (θ̆) and ṡφ (φ1 ) ṡφ (φ̆), (8.51)
Example 8.3 Figures 8.5(a) and (b) show values of CRB[θ1 |Aθ|φ ]
(Eq. (8.39)) as a function of θ and φ, in the case of two unit-power
emitters of common elevation φ but separated by ∆θ = 1◦ about azimuth θ.
The bound is evaluated for the 24-element Y-shaped and circular arrays
assuming the availability of L = 100 snapshots and signal-to-noise ratio
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
of 10 dB (for each emitter). In Fig. 8.5(c) and (d) the exercise is repeated
for CRB[φ1 |Aφ|θ ] with emitters of common azimuth θ but separated by
∆φ = 1◦ about elevation φ. While the circular array again exhibits uni-
form accuracy for all azimuths, the Y-shaped array shows fluctuations with
a period of 60◦ due to its special shape. Again, the performances of the
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
Fig. 8.6 View of the H12 subspace (plane): An illustrative representation of the geom-
etry at ‘resolution threshold’ (circular approximation)
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
These two subspace are illustrated in Figs. 8.7 and 8.8 together with
the two manifold vectors a(s1 ) and a(s2 ), and their associated uncertainty
spheres. Note that in these figures the manifold curve may be the manifold
of a linear array or a p-curve lying on the manifold surface of a planar
array.
∆s ∆sdet-thr (8.55)
Based on this definition and once again using the circular approxima-
tion/representation of the array manifold, the scenario of Fig. 8.2, at the
detection threshold, pertains:
∆sdet-thr = κ̂−1
1 ∆ψ
∆ψ1 ∆ψ2
arcsin(κ̂1 σe1,d ) + arcsin(κ̂1 σe2,d )
= (8.57)
κ̂1
Noting that arcsin(x) x for x 1, one may write
1
∆pdet-thr = (σe + σe2,d ) (8.61)
π r sin p̆ 1,d
where
σei,d =
P⊥
H ⊥
Hdet a(si ) = a (si )PHdet a(si ) for i = 1, 2 (8.62)
∆s ∆sdet-thr (8.63)
Fig. 8.10 View of the H12 subspace (plane): an illustrative representation of the geom-
etry at “resolution threshold” on the H12 subspace (circular approximation).
at the point s̆. Once again using the circular representation of the array
manifold the situation of Definition 8.2 may be illustrated as in Figs. 8.10
and 8.11. With reference to these two figures, it can be seen that σe2,r (see
Fig. 8.11) and µx2 (see Fig. 8.10) as well as σe1,r and µx1 may be given as
follows:
σe1,r = P⊥ ⊥
Hres a(s1 ) and σe2,r = PHres a(s2 ) (8.64)
σe σe
µx1 = 1,r and µx2 = 2,r (8.65)
sin γ sin γ
Fig. 8.11 View orthogonal to Hres and H12 subspaces (planes): an illustrative view of
the geometry at “resolution threshold”.
σe
∆ψ2 = arccos 1 − κ̂1 2
sin γ
σ e2
2κ̂1 (8.67b)
sin γ
P⊥
H1,2 a(s̆) 1/κ̂1 1
cos γ = = √ ⇒ sin2 γ = 1 − cos2 γ = 1 −
a(s̆) N N κ̂21
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
where
H1,2 = L{[u1 (s̆), u2 (s̆)]}
Therefore,
9
4 √ √
∆sres-thr = 4 σe1,r + σe2,r (8.69)
κ̂21 − 1
N
where
σei,r =
P⊥
H ⊥
Hres a(si ) = a (si )PHres a(si ), for i = 1, 2 (8.71)
4 1
∆sres-thr = 4 4 ∆pres-thr = 4
1 πr sin p̆
2 1 2
κ 1− N
κ1− N
√ √ √ √
× σe1,d + σe2,d × σe1,r + σe2,r
σei,d = PHdet a(si ), i = 1, 2
with
σe = ⊥ a(s ),
PH i i = 1, 2
i,r res
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
where the fourth factor is dependent on the available SNR and the number
of array output snapshots L, and will be modelled next.
1
σe2 ∝ (8.73)
SNR
This implies that the uncertainty sphere shrinks, or equivalently, that there
would be a gradual decrease in its radius σe , thus providing an effective
radius σe according to the following model:
/ 9
σ2 1
σe = = ; (in arc lengths) (8.74)
2LP 2(SNR × L)
It should be pointed out that the factor “2” appearing in the denominator
is a direct result of using complex numbers to represent signal envelopes in
the array model.
The shrinking rate of the uncertainty sphere depends also on the esti-
mation algorithm deployed. Some algorithms exhibit better properties than
others thereby “masking” the overall array performance and this makes it
difficult to identify the full detection-resolution capabilities inherent in the
array structure. Thus we may introduce the factor C (with 0 < C ≤ 1) to
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
The equality of the single source Cramer Rao lower bound with the square
of the uncertainty sphere radius can be interpreted as follows:
Based on the previous discussion and the model of Eq. (8.74), at the detec-
tion and resolution “thresholds” the uncertainty spheres can be expressed as
9 @
@
1 @
σei,d = aH (si )P⊥
Hdet a(si ) = @ (8.76a)
2(SNR × L)C @
det
9 @
@
1 @
H
σei,r = a (si )P⊥
Hres a(si ) = @ (8.76b)
2(SNR × L)C @
res
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
Substituting Eq. (8.76a) back into Eqs. (8.60) and (8.69), with C = 1, it
follows that:
(SNR1 × L) κ̂1 − N1 P2
Since ṡ(p̆) = πr sin p̆, Eq. (8.78) indicates that the resolution capabilities
of a linear array (or ELA) are maximum at broadside (p = 90◦ ) and zero
along the endfire directions (p = 0◦ or 180◦ ). Zero resolution along end-fire
directions is expected since, due to cylindrical symmetry, two signals at
directions which are symmetrical to the array axis cannot be distinguished,
irrespective of their angular separation.
For the special case of equi-powered sources (P1 = P2 = P ) the ultimate
product (signal-to-noise ratio times L) required to detect and resolve two
sources at bearings p1 and p2 = p1 + ∆p may be written as
√
2 1
∆pdet-thr ≈ √
ṡ(p̆) SNR × L
9 (8.79)
1 4 32 1
∆pres-thr ≈ 2 √
ṡ(p̆) κ̂1 − N1 4 SNR × L
falls more rapidly than the resolution threshold as the number of sensors
N increases for a constant SNR × L.
It is interesting to note that the resolution threshold ∆pres-thr of a linear
array is a function of sin(p). This indicates that an algorithm which per-
forms an exhaustive search of the manifold (e.g. MUSIC) should evaluate
its cost function along a non-uniform grid of azimuths whose values follow
a sinusoidal distribution. Such a non-uniform angular grid corresponds to a
uniform arc length grid along the manifold and would allow full use of the
array’s resolution capabilities.
Finally Eq. (8.77) can be used to derive other useful system require-
ments. For example, the number of snapshots, Ldet-thr , required to detect
two sources at bearings p1 and p2 = p1 + ∆p or Lres-thr to resolve them,
may be written as shown below.
1 P2
∆s4 κ̂21 −
N
∆s = πr| cos p2 − cos p1 |
where
∆s1 κ̂1 1, ∆s2 κ̂1 1
(8.83)
Above, the expressions are also provided for the minimum product
(SNR1 × L) required for a source of power P1 to be detected and then
resolved in the presence of a source of power P2 after L snapshots when the
two sources are separated by ∆s.
Finally Table 8.4 summarizes the main expressions for the special case
of two equi-power sources.
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
Detection Resolution
√
∆pdet-thr ≈ 2 √ 1
∆pres-thr ≈ 1 4 32
1 √
4
1
ṡ(p̆) SNR×L ṡ(p̆) (κ21 − N ) SNR×L
(SNR × L)det-thr = 2
∆s2
(SNR × L)res,thr = 32
1
1− N )
∆s4 (κ2
2 1 32 1
Ldet-thr = ∆s2 SNR
Lres-thr = 1
1− N )
∆s4 (κ2 SNR
8.6 Comments
for large N and inclination angle ζ = 90◦ . That is for symmetric arrays
with a large number of sensors.
It is important to note that for two equi-powered sources the above expres-
sion is simplified to
∆p ṡ(p̆)
CRB[p1 ]res .
4 ṡ(p1 )
∆p
→ (8.88)
4
[as (p̆−p1 )→0]
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
Fig. 8.13 Simulation results of the CRB (exact and approximate expressions) as a
function of azimuth separation, averaged over 50 trials. The resolution threshold and the
single-source CRB are also shown for comparison.
1 −1
CRB[p1 |A] = 1 P ȧ1 )
(ȧH
2 SNR1 × L|g(p1 )|2
1 2
(SNR1 × L)|g(p1 )|2 (∆s)22 (p1 )(
κ21 (p) − 1/N )
!2
1 |g(p1 )|2 P1
(SNR1 × L)det = 1+
2(∆p × (p))2 |g(p1 )|2 |g(p2 )|2 P2
!4
2 |g(p1 )|2 P1
(SNR1 × L)res = 1+ 4
(∆p × (p))4 (
κ21 (p) − 1
N
)|g(p1 )|2 |g(p2 )|2 P2
or
◦
(θ = Θo , φ = 180◦ −α1 ) and (θ = Θo , φ = 180◦ −α2 ) if α1 , α2 > 90
(8.90)
Proof.
(i) According to Corollary 5.2 and Theorem 5.1, all members of the family
of α-curves are identical, i.e. their differential geometry is independent
of β.
(ii) From Corollary 5.5 and Theorem 5.1, each member of the α-curve
family is identical to the combination of the φ-curves corresponding
to θo = Θo and θo = Θo + 180◦ . This implies that α-estimation and
φ-estimation performance must by necessity be equivalent.
Fig. 8.14 The loci of two DOAs separated by ∆α = 2◦ for all values of β and the
corresponding DOA loci on the (θ, φ) parameter plane revealing the “conic” nature of
the array performance.
8.8 Summary
8.9 Appendices
However,
2
κ−1
1
= (8.93)
sin ζ
κ̂−1
1 = PH12 a(s̆)
κ−1
1
= (8.94)
sin ζ
“Circular” “Y”
Sensor
number x y x y
⊥ H ⊥ H ⊥ −1 H ⊥
11 Pa1 u11 − u11 Pa1 Ar (Ar Pa1 Ar )
= uH Ar Pa1 u11
H ⊥ −1 H
= 1 − uH
11 Ar (Ar Pa1 Ar ) Ar u11
−1
1 H
= 1 − u11 Ar Ar Ar − Ar a1 a1 Ar
H H H
AH
r u11 (8.96)
N
⊥ N |uH11 a2 |
2
11 PA u11 = 1 −
uH (8.97)
N 2 − |aH1 a2 |
2
The physical proximity of the arriving signals allows the use of local
differential geometry for the evaluation of Eq. (8.38) (or, of Eq. (8.36) for
the general case of M > 2).
or
/
1 2
u1 (s1 )H a(s2 ) ∆a 1− ∆a κ̂21 (s̆) (8.100)
4
(b) Evaluation of the term aH H
1 a2 = a(s1 ) a(s2 ). Clearly:
The array manifold has constant norm and hence lies on the surface
of a hypersphere. Consequently, for sufficiently small ∆s, vectors a(s̆)
and ∆a are strictly orthogonal as indicated in Fig. 8.16. Using simple
trigonometry:
√
a(s1 )H ∆a N ∆a cos(90◦ + γ) (8.102)
√
− N ∆a sin(γ) (8.103)
√ ∆a 1
− N ∆a √ = − ∆a2 (8.104)
2 N 2
Therefore
Substituting Eqs. (8.104) and (8.105) back into Eq. (8.97) and then
using Eq. (8.38):
: ;−1
1 N ∆a2 1 − 14 ∆a2 κ̂21 (s̆)
CRB[p|A] = 1−
2(SNR1 × L)ṡ(p1 )2 N 2 − (N − 12 ∆a2 )2
Fig. 8.16 For sufficiently small ∆s, the vectors a(s̆) and ∆a are strictly orthogonal.
July 27, 2004 13:9 WSPC/Book Trim Size for 9in x 6in chap08
−1
1 N (1 − 14 ∆a2 κ̂21 (s̆))
= 1−
2(SNR1 × L)ṡ(p1 )2 N − 14 ∆a2
1 4N − ∆a2
=
2(SNR1 × L)ṡ(p1 )2 ∆a2 (N κ̂21 (s̆) − 1
So finally:
1 2
CRB[p|A] = (8.106)
(SNR1 × L) ṡ(p1 )2 (∆s)2 κ̂21 (p̆) − 1
N
Bibliography
215
July 6, 2004 9:31 WSPC/Book Trim Size for 9in x 6in bib
Index
217
July 6, 2004 9:32 WSPC/Book Trim Size for 9in x 6in INDEX