Extended version of an invited paper at IEEE CAMSAP 2017.

Position Aided Beam Alignment for Millimeter
Wave Backhaul Systems with Large Phased Arrays
George C. Alexandropoulos

arXiv:1701.03291v3 [] 12 Oct 2017

Mathematical and Algorithmic Sciences Lab, Paris Research Center, Huawei Technologies France SASU,
20 Quai du Point du Jour, 92100 Boulogne-Billancourt, France
email: george.alexandropoulos@huawei.com

Abstract—Wireless backhaul communication has been recently
realized with large antennas operating in the millimeter wave
(mmWave) frequency band and implementing highly directional
beamforming. In this paper, we focus on the alignment problem of
narrow beams between fixed position network nodes in mmWave
backhaul systems that are subject to small displacements due
to wind flow or ground vibration. We consider nodes equipped
with antenna arrays that are capable of performing only analog
processing and communicate through wireless channels including
a line-of-sight component. Aiming at minimizing the time needed
to achieve beam alignment, we present an efficient method that
capitalizes on the exchange of position information between the
nodes to design their beamforming and combining vectors. Some
numerical results on the outage probability with the proposed
beam alignment method offer useful preliminary insights on the
impact of some system and operation parameters.

I. I NTRODUCTION
Millimeter wave (mmWave) communication [1] is a promising technology for addressing the high throughput requirement
for the fifth generation (5G) mobile networks [2], [3]. Shortrange mmWave communication at the unlicensed band of 60
GHz is already standardized in IEEE 802.11ad [4] and initial
theoretical investigations on mmWave cellular systems [1],
[5] have identified their potentials and key challenges. The
mmWave frequencies have been also recently considered for
the wireless backhaul communication of small cells [6], [7],
which are expected to be densely deployed as an efficient
means for increasing the geographic spectrum reusability [8].
Reliable mmWave backhauling depends on very directional
communication, which are implemented in practice either with
antennas having large apertures or with large phased antenna
arrays [6]. By exploiting the fact that the wavelength at very
high frequencies is very small, large phased antennas can be
cheaply packed into small form factors and, thus, have been
effectively used in realizing highly directional BeamForming
(BF) supporting long outdoor links. To achieve the full benefit
from BF in a communication link between multi-antenna
nodes, the entire channel state information needs to be available at both communication ends. However, this information is
hard to acquire in mmWave systems due to the low coherence
time, the Radio Frequency (RF) hardware limitations, and
the small Signal-to-Noise Ratio (SNR) before BF. Although
recent theoretical works [5], [9]–[11] capitalized on the spatial
sparsity of mmWave channels [12] to estimate portions or the
entire channel gain matrix, the presented approaches required

lengthy training phases to estimate the channel in multiple
directions using complex compressed sensing algorithms. Another family of approaches (e.g., [13]–[18]) for efficient BF is
based on beam switching between the communicating nodes in
order to find pair(s) of beams from their available codebooks
meeting a predefined performance threshold. When such beam
pair(s) is(are) found, beam alignment is considered to be
achieved and no further beam searching is needed.
In this paper, we focus on the mmWave backhaul system of
[18] and present a robust beam alignment method for wireless
channels including a Line-Of-Sight (LOS) component. The
proposed method intends at achieving beam alignment in at
most two rounds of control information exchange. The core
idea of the method lies on the control information, which is
considered to be the new positions of the nodes after each
declared beam misalignment event, on its exchange, and on its
utilization in designing the BF technique at all communication
ends. The position information of a node is assumed to be
available to it via an attached position sensor. It is noted that
positioning sensors have been also considered in [19], [20],
however, they were adopted for identifying and circumventing
beam misalignment only at the transmit node.
Notation: Vectors and matrices are denoted by boldface
lowercase letters and boldface capital letters, respectively.
The transpose and Hermitian transpose of a matrix A are
denoted by AT and AH , respectively, diag{a} denotes a
square diagonal matrix with a’s elements in its main diagonal,
whereas In (n ≥ 2) is the n × n identity matrix. The ith
element of a and the (i, j)th element of A are denoted by
[a]i and [A]i,j , respectively, and ||A||F gives the Frobenius
norm of A. C represents the complex number set, card(F )
is the cardinality of set F , | · | denotes the amplitude of
a complex number, and E{·} is the expectation operator.
d(M, N ) denotes the length of the line segment
 connecting
the points M and N . Notation x ∼ CN 0, σ 2 indicates that
x is a circularly-symmetric complex Gaussian random variable
with zero mean and variance σ 2 , while x ∼ U (α, β) represents
a uniformly distributed random variable in [α, β].
II. S YSTEM

AND

C HANNEL M ODELS

Suppose a wireless backhaul communication system operating in the mmWave frequency band and consisting of two half
duplex multi-antenna transceiver nodes A and B (see Fig. 1).
The fixed positions of the nodes defines a Cartesian coordinate

y

(B)

Mt
dt

θ1(t+1)

(B)

dt-yt+1

(B)

Mt+1

(A)

Mt+1

(A)
yt+1

φ1(t+1)

Ot+1,t+1

(A)

Mt
(A)
-xt+1

0

(B)
xt+1

x

Fig. 1. The considered point-to-point wireless backhaul communication
system. The positions of the nodes A and B at the time instant t are denoted
(A)
(B)
by the points Mt and Mt , respectively, and define a Cartesian coordinate
system. At the time instant t + 1 both nodes A and B move to new positions
(A)
(B)
represented by the points Mt+1 and Mt+1 , respectively. These new positions
determine the coordinates of the point Ot+1,t+1 ; its first subscript refers to
(A)
(B)
Mt+1 and the second to Mt+1 .

system, according to which the positions of nodes A and B
(A)
at the time instant t are represented by the points Mt and
(B)
Mt , respectively. Their respective coordinates are (0, 0) and
(0, dt ), where dt denotes the physical distance between the
nodes at this time instant. Regarding the mmWave access,
each node is assumed to be equipped with one RF chain and
is capable of realizing only analog transmit BF when being in
transmit mode, and only analog receive combining when being
in receive mode. Node A is assumed to have a large uniform
linear antenna array (ULA) with NA elements, whereas node
B is equipped with a large NB -element ULA. We hereinafter
assume for simplicity that the positions of the nodes coincide
with the positions of the centers of their antenna arrays.
To establish wireless communication between nodes A and
B, i.e., to meet a minimum performance threshold for reliable
information exchange, a control phase for channel access is
adopted. We assume that communication is realized in discrete
time instants each including a control and a data phase. The
control phase constitutes of several time slots during which
data communication is being set up. We assume that during
each discrete time instant the channel remains constant, but it
may change between different instances. During the control
phase for channel access, control information needs to be
exchanged between nodes A and B. This control information
exchange requires in general a much less stringent performance threshold than data exchange. This means that BF
gain is not necessary for control signaling, hence, its received
SNR can be afforded to be low. As such, we assume that
control signaling exchange can be in principle handled by the
considered communication systems.
Before signal transmission from node A, the unit power
data stream s ∈ C (chosen from a discrete modulation set)
is processed by a BF vector vA ∈ CNA ×1 , and upon signal
reception at node B, the combining vector uB ∈ CNB ×1 is

used for processing the received signal. Similarly, node B
makes use of the BF vector vB ∈ CNB ×1 when transmitting,
and node A utilizes the combining vector uA ∈ CNA ×1
when being in receive mode. Due to practical limitations with
the considered analog antenna arrays [5], [18], the analog
antenna weights at both nodes are supposed to take values
from discrete sets. In particular, we assume that vA and uA
at node A are chosen from the finite set FA of predefined
NA -element vectors. The same holds for node B; vB and uB
belong to the finite set FB of predefined NB -element vectors.
Without loss of generality, it is assumed that for any f ∈ FA
with f ∈ CNA ×1 and any z ∈ FB with z ∈ CNB ×1 holds
|[f ]i |2 , NA−1 with i = 1, 2, . . . , NA and |[z]j |2 , NB−1 with
j = 1, 2, . . . , NB . When node A transmits information to node
B, the output signal after the combiner at node B can be
mathematically expressed as
√
H
yB = puH
(1)
B HBA vA s + uB nB ,
where p is the transmit power, HBA ∈ CNB ×NA denotes
the channel gain matrix between nodes B and A, and nB ∈
CNB ×1 represents the zero-mean Additive White Gaussian
2
Noise (AWGN) vector with covariance matrix σB
INB .
We adopt a geometric channel model with L scatterers
similar to [5], [12]. We assume that during each discrete
time instant the wireless channel remains constant, but it may
change between different instances. According to this channel
model, HBA included in (1) is expressed as
HBA , AB (θ) diag{a}AH
A (φ) ,

(2)

where AA (φ) ∈ CNA ×L , with φ , [φ1 φ2 · · · φL ], and
AB (θ) ∈ CNB ×L , with θ , [θ1 θ2 · · · θL ], are defined as
AA (φ) , [aA (φ1 ) aA (φ2 ) · · · aA (φL )] ,

(3a)

AB (θ) , [aB (θ1 ) aB (θ2 ) · · · aB (θL )] .

(3b)

In (3), variable φℓ ∈ [0, 2π] with ℓ = 1, 2, . . . , L denotes the
ℓth path’s Angle of Departure (AoD) from node A and variable
θℓ ∈ [0, 2π] represents the ℓth path’s Angle of Arrival (AoA) at
node B. Following the investigations in [21], we assume that
the 1st channel path is a LOS one with energy much larger
than each of the rest L − 1 paths. In addition, aA (φℓ ) ∈
CNA ×1 and aB (θℓ ) ∈ CNB ×1 are the array response vectors
at nodes A and B, respectively (these vectors are given by [5,
eq. (5)] for ULAs). In (2), a ∈ CL×1 includes the channel
path gains αℓ ∀ℓ = 1, 2, . . . , L. We further assume that each
path’s amplitude is Rayleigh distributed and, in particular, that
each αℓ ∼ CN (0, NA NB PL ), where PL denotes the average
pathloss between nodes B and A.
In practical deployments of wireless backhaul systems,
the antenna arrays of the communicating nodes are usually
mounted on outdoor structures that are exposed to wind flow
and gusts. These structures are susceptible to movement (or
sway) due to wind or ground vibration, which might cause
unacceptable Outage Probability (OP) if beam alignment is not
frequently performed [18]. To capture the random movements

of the antenna arrays at both nodes A and B, we model the
displacements at the x-axis and y-axis for both of them as
(A)
(A)
x(A) ∼ U(−x(A)
∼ U(−ys(A) , yn(A) ),
w , xe ), y

(4a)

(B)
x(B) ∼ U(−xw
, xe(B) ), y (B) ∼ U(ys(B) , yn(B) ).

(4b)

(A)

(A)

(B)

(B)

In the latter expressions, xw , xe , xw , and xe reveal the
(A)
(A)
position limits in the x-axis for both nodes, while ys , yn ,
(B)
(B)
ys , and yn represent their position limits in the y-axis.
III. P OSITION A IDED B EAM A LIGNMENT
Suppose that at the time instant t + 1 both nodes A
and B move to new positions represented by the points
(B)
(A)
Mt+1 and Mt+1 , respectively, with respective coordinates
(A)
(A)
(B)
(B)
(xt+1 , yt+1 ) and (xt+1 , dt − yt+1 ). According to the presented channel model, the channel between node B and node
A at the discrete time instant t + 1 can be expressed as
HBA (t + 1) , AB (θ(t + 1)) diag{a(t + 1)}AH
A (φ(t + 1)),
with φ(t + 1) , [φ1 (t + 1) φ2 (t + 1) · · · φL (t + 1)] and
θ(t+1) , [θ1 (t+1) θ2 (t+1) · · · θL (t+1)] denoting the AoDs
and AoAs, respectively, resulting from the new node positions,
while a(t+1) includes the channel gains. The beam alignment
objective in this time instant is to design the BF vector at node
A and the combining vector at node B as
{vA (t+1), uB (t+1)} =

max

f ∈FA ,z∈FB

2

zH HBA (t+1)f , (5)

which implies BF gain maximization. If perfect knowledge of
HBA (t+1) was available at both nodes, nodes A and B should
use the principal right and left singular vectors of HBA (t + 1)
as the BF and combining vectors, respectively, to maximize
the BF gain [22]. However, neither HBA (t + 1) for any t can
be estimated at any of the nodes due to the assumed hardware
limitations nor the nodes are capable of realizing any arbitrary
vector (they possess the finite vector sets FA and FB ).
A. Proposed Method
Let us consider that, upon installation of the considered
mmWave backhaul system at the initial discrete time instant
t, both nodes A and B are aware of the coordinate system
(A)
(B)
defined by the coordinates of the points Mt and Mt . In
the next discrete time instant t + 1, suppose that both nodes
(A)
(B)
move to the new position points Mt+1 and Mt+1 . Each node’s
displacement in the x and y axes is assumed to be available
to the node. The latter indicates that each node is aware of
the coordinates of its new position, i.e., node A learns the
(A)
(A)
(B)
(B)
coordinates (xt+1 , yt+1 ) and B obtains (xt+1 , dt − yt+1 ).
If also the coordinates of the position of node B become
available to A (through, for example, a dedicated control
channel), then the latter node may estimate the AoD of the
LOS channel path for its transmission at this time instant as
(
(A)
(B)
arcsin (gt+1 ) ,
xt+1 ≥ xt+1
(6)
φ̂1 (t + 1) =
(A) ,
(B)
π − arcsin (gt+1 ) , xt+1 < xt+1

where the positive real gt+1 is given by


(B)
d Ot+1,t+1 , Mt+1


gt+1 =
(A)
(B)
d Mt+1 , Mt+1

(7)

(B)

(A)

with point Ot+1,t+1 having the coordinates (xt+1 , yt+1 ). From
the specific nodes’ positions at the time instant t + 1 yields
(A)

gt+1 = r

(B)

dt − yt+1 − yt+1
2 .
2 
(A)
(B)
(A)
(B)
xt+1 + xt+1 + dt − yt+1 − yt+1

(8)

In a similar way, if node A shares its coordinates at the same
time instant with B, then node B can estimate the AoA of the
LOS path for A’s transmission as
(
(A)
(B)
π + φ1 (t + 1), xt+1 ≥ xt+1
(9)
θ̂1 (t + 1) =
(B)
(A) .
2π − φ1 (t + 1), xt+1 < xt+1
The AoD and AoA of the LOS path can be obtained similarly
if node B transmits and node A operates in receive mode.
Capitalizing on (6) and on our channel model for each time
instant t+1, we propose that transmit node A uses its estimate
φ̂1 (t + 1) to realize a beam steering towards the direction of
receive node B. To accomplish this, it searches inside its beam
codebook FA for the vector that is closest to aA (φ̂1 (t + 1)).
In mathematical terms and without loss of generality, node A
designs its BF vector at each time instant t + 1 as


vA (t + 1) = min f − aA φ̂1 (t + 1)
f ∈FA

2
F

.

(10)

Similarly, receive node B uses its estimate θ̂1 (t + 1) at each
t + 1 to realize a beam as close as possible to the direction of
node A. It, therefore, constructs its combining vector as


uB (t + 1) = min z − aB θ̂1 (t + 1)
z∈FB

2
F

.

(11)

By using the vectors given by (10) and (11) at nodes A and
B, respectively, the instantaneous received SNR at node B is
2
2
given by γt+1 = |µt+1 | /σB
, where µt+1 ∈ C is defined as
µt+1 =

L
X
ℓ=1

×



αℓ (t + 1)uH
B (t + 1)aB θ̂ℓ (t + 1)

aH
A



φ̂ℓ (t + 1) vA (t + 1).

(12)

It is noted that for the special case of L = 1 LOS channel
path, perfectly estimated AoD and AoD for this path, and
infinite resolution beam codebooks at both nodes A and B
yields µt+1 = α1 (t + 1). This indicates that for this ideal
case the vectors given by (10) and (11) maximize the BF gain
2
described in (5) and given by |µt+1 | . The proposed beam
alignment method is summarized in Algorithm 1.

IV. N UMERICAL R ESULTS

AND

D ISCUSSION

In this section we evaluate the performance of the proposed
beam alignment method over the considered mmWave channel
model. We particularly evaluate the OP performance defined
as the probability that the instantaneous SNR falls below a
minimum SNR threshold γo . To carry out this evaluation, we
have simulated 104 channel samples according to (2) with
2
normalized σB
and PL = d−3.75
, where each channel sample
1

100

Outage Probability (OP)

Algorithm 1 Position Aided Beam Alignment
Initialization: Construct a coordinate system upon installation of nodes A and B at the initial time instant
t. Determine the minimum required SNR γo for data
communication.
1: for n = t + 1, t + 2, . . . do
Node A Recovery Phase 1:
2:
if γn ≥ γo , then
3:
Set vA (n) = vA (n − 1), uB (n) = uB (n − 1),
and transmit data.
4:
else
(A)
5:
Obtain the coordinates of the new position Mn .
6:
Compute the AOD φ̂1 (n) using (6) and the
(A)
(B)
coordinates of the positions Mn and Mn−1 , and
the point On,n−1 .
7:
Compute aA (φ̂1 (n)) and design vA (n) using (10).
8:
end if
9:
if γn ≥ γo then
10:
Set vA (n) as in step 7, uB (n) = uB (n − 1),
and transmit data.
11:
else
(A)
12:
Send the coordinates of Mn to Node B.
13:
end if
Node B Recovery Phase:
14:
Upon reception of a control signal with the coordinates
(A)
(B)
of Mn , obtain the coordinates of Mn .
15:
Compute the AOA θ̂1 (n) using (9) and the coordinates
(A)
(B)
of the positions Mn and Mn , and On,n .
16:
Compute aB (θ̂1 (n)) and design uB (n) using (11).
17:
if γn ≥ γo , then
18:
Trigger node A to transmit data.
19:
Set vA (n) as in step 7, uB (n) as in step 16,
and transmit data.
20:
else
(B)
21:
Send the coordinates of Mn to Node A and use
uB (n) designed in step 16.
22:
end if
Node A Recovery Phase 2:
23:
Upon reception of a control signal with the coordinates
(B)
of Mn , compute the AOD φ̂1 (n) using (6) and the
(A)
(B)
coordinates of Mn and Mn , and On,n .
24:
Compute aA (φ̂1 (n)) and design vA (n) using (10).
25:
Set vA (n) as in step 23, uB (n) as in step 16,
and transmit data.
26: end for

10-1

No Action
A1
A1 and B
A1, B and A2
Optimum BF
10

-2

1

1.5

2

2.5

3

3.5

4

SNR Threshold γo in dB

Fig. 2. OP versus the SNR threshold γo in dB for NA = NB = 16,
2 = 5dB, card(F ) = card(F ) = 17, L = 3, κ = 13.2dB, and
p/σB
A
B
different phases of the proposed method. Curves for the unfeasible optimum
BF case and the case of no action for beam alignment are included.

appears at one discrete time instant. We have considered
Ricean fading channels with the κ-factor denoting the ratio
of the energy in the LOS channel path to the sum of the
energies in the other non LOS paths [21]. The distance of
the nodes A and B at the initial time instant t = 1 is
assumed to be d1 = 10m and the random displacements of the
nodes in the x-axis and y-axis due to wind gusts are obtained
(B)
(B)
(A)
(A)
from (4) with xw = xe = xw = xe = 1.5m and
(A)
(A)
(B)
(B)
yw = ye = yw = ye = 1.5m. The beam codebooks
FA and FB are constructed by quantizing the feasible sets
of departure and arrival angles, respectively. Specifically, for
the ith BF vector belonging in FA with i ∈ card(FA )
the departure angle χi takes discrete values with step size
21−q (χmax − χmin ) within [χmin , χmax ], where q represents
the number of angle quantization bits. Similarly, for each jth
combining vector inside FB with j ∈ card(FB ) the arrival
angle ψj ∈ [ψmin , ψmax ] has been quantized with q bits. We
have set χmin = 60o , χmax = 120o, ψmin = 240o , and
ψmax = 300o with respect to the nodes’ orientation.
In Fig. 2, we plot the OP with the proposed beam alignment
method as a function of the SNR threshold γo in dB for
2
the transmit SNR p/σB
= 5dB and NA = NB = 16,
card(FA ) = card(FB ) = 17 resulting from q = 5, L = 3,
and the Ricean factor κ = 13.2dB. Within this figures, the OP
curves for the case of no action for beam alignment and the
case of optimum BF [22] when perfect channel knowledge
is available are also sketched. As for the proposed method,
the performance for different sequences of phases is also
demonstrated. In particular, we provide the performance for
the case where only Node A Recovery Phase 1 is used, denoted
by A1; the case where Node A Recovery Phase 1 followed
by Node B Recovery Phase are used, denoted by A1 and
B; and the case where all phases are utilized, denoted by
A1, B and A2. As seen and as expected, OP degrades with
increasing γo . It is also shown that the exchange of position
information improves this probability for any γo value. In fact,

the availability of nodes’ position at both nodes, with only 2
control signals according to the proposed method, results in the
best performance. Evidently, actions from only the transmitter
or exchange of only the latter’s position information result in
poor OP. We have also compared the proposed beam alignment
method with an exhaustive beam search similar to [15] that
seeks to find the first beam pair meeting γo whenever beam
misalignment occurs. For the parameter settings of Fig. 2
we quote the following representative results: i) the proposed
method achieves 98% of the average received SNR of optimum
BF irrespective of γo , while the exhaustive search only the
28% for γo = 1dB after 30 control signals on average; and ii)
for γo = 3dB the exhaustive search reaches the 72% of the
optimum BF SNR requiring on average 153 control signals.
V. C ONCLUSION
In this paper, we investigated the problem of beam alignment between network nodes in wireless backhaul communication systems operating in the mmWave frequency band.
We considered the practical case of random movement of the
large phased antenna arrays of the nodes due to wind flow or
ground vibration, and presented a simple model that captures
their small displacements. A robust beam alignment method
for wireless environments including a LOS component was
presented that capitalizes on the exchange of position information between the nodes to design their BF and combining
vectors. Through representative OP performance evaluation
results the impact of some key parameters on the performance
of the proposed method was highlighted.
R EFERENCES
[1] T. S. Rappaport et al., “Millimeter wave mobile communications for 5G
cellular: It will work!” IEEE Access, vol. 1, pp. 335–349, May 2013.
[2] F. Boccardi et al., “Five disruptive technology directions for 5G,” IEEE
Commun. Mag., vol. 52, no. 2, pp. 74–80, Feb. 2014.
[3] J. G. Andrews et al., “What will 5G be?” IEEE J. Sel. Areas Commun.,
vol. 32, no. 6, pp. 1065–1082, Jun. 2014.
[4] “IEEE wireless LAN MAC and PHY specifications– Amendment 3:
Enhancements for very high throughput in the 60 GHz band,” IEEE
Std. 802.11ad, 2012.
[5] A. Alkhateeb et al., “Channel estimation and hybrid precoding for
millimeter wave cellular systems,” IEEE J. Sel. Topics Signal Process.,
vol. 8, no. 5, pp. 831–846, Oct. 2014.
[6] S. Chia et al., “The next challenge for cellular networks: Backhaul,”
IEEE Microwave Mag., vol. 10, no. 5, pp. 54–66, Aug. 2009.
[7] X. Ge et al., “5G wireless backhaul networks: Challenges and research
advances,” IEEE Netw., vol. 28, no. 6, pp. 6–11, Dec. 2014.
[8] N. Bhushan et al., “Network densification: The dominant theme for
wireless evolution into 5G,” IEEE Commun. Mag., vol. 52, no. 2, pp.
82–89, Feb. 2014.
[9] Z. Marzi et al., “Compressive channel estimation and tracking for large
arrays in mm-wave picocells,” IEEE Trans. Signal Process., vol. 11,
no. 3, pp. 514–527, Apr. 2016.
[10] A. Bazzi et al., “A comparative study of sparse recovery and compressed
sensing algorithms with application to AoA estimation,” in Proc. IEEE
SPAWC, Edinburgh, UK, 3-6 Jul. 2016, pp. 1–5.
[11] G. C. Alexandropoulos and S. Chouvardas, “Low complexity channel
estimation for millimeter wave systems with hybrid A/D antenna processing,” in Proc. IEEE GLOBECOM, Washington D.C., USA, 4-8 Dec.
2016, pp. 1–6.
[12] H. Zhang et al., “Channel modeling and MIMO capacity for outdoor
millimeter wave links,” in Proc. IEEE WCNC, Sydney, Australia, 18-21
Apr. 2010, pp. 1–6.

[13] K. Hosoya et al., “Multiple sector ID capture (MIDC): A novel beamforming technique for 60 GHz band multi-Gbps WLAN/PAN systems,”
IEEE Trans. Antennas Propag., vol. 63, no. 1, pp. 81–96, Jan. 2015.
[14] J. Wang et al., “Beam codebook based beamforming protocol for multiGbps millimeter-wave WPAN systems,” IEEE J. Sel. Areas Commun.,
vol. 27, no. 8, pp. 1390–1399, Oct. 2009.
[15] C. Jeong et al., “Random access in millimeter-wave beamforming
cellular networks: Issues and approaches,” IEEE Commun. Mag., vol. 53,
no. 1, pp. 180–185, Jan. 2015.
[16] C. Barati et al., “Directional cell discovery in millimeter wave cellular
networks,” IEEE Trans. Wireless Commun., vol. 14, no. 12, pp. 6664–
6678, Dec. 2015.
[17] V. Raghavan et al., “Beamforming tradeoffs for initial UE discovery in
millimeter-wave MIMO systems,” IEEE J. Sel. Topics Signal Process.,
vol. 10, no. 3, pp. 543–559, Apr. 2016.
[18] S. Hur et al., “Millimeter wave beamforming for wireless backhaul and
access in small cell networks,” IEEE Trans. Commun., vol. 61, no. 10,
pp. 4391–4403, Oct. 2013.
[19] R. Maiberger et al., “Location based beamforming,” in Proc. IEEE
IEEEI, Eliat, Israel, 17-20 Nov. 2010, pp. 184–187.
[20] A. W. Doff, K. Chandra, and R. V. Prasad, “Sensor assisted movement
identification and prediction for beamformed 60 GHz links,” Available
at https://arxiv.org/pdf/1502.04054.pdf, 2015.
[21] Z. Muhi-Eldeen et al., “Modelling and measurements of millimetre
wavelength propagation in urban environments,” IET Microwaves, Ant.
Propag., vol. 4, no. 9, pp. 1300–1309, May 2010.
[22] D. J. Love and R. W. Heath, Jr., “Equal gain transmission in multipleinput multiple-output wireless systems,” IEEE Trans. Commun., vol. 51,
no. 7, pp. 1102–1110, Jul. 2003.

