Sub-Nyquist Channel Estimation over IEEE 802.11ad Link

arXiv:1705.03943v2 [cs.IT] 7 Aug 2017

Kumar Vijay Mishra and Yonina C. Eldar

Abstract—Nowadays, millimeter-wave communication centered at the
60 GHz radio frequency band is increasingly the preferred technology
for near-field communication since it provides transmission bandwidth
that is several GHz wide. The IEEE 802.11ad standard has been
developed for commercial wireless local area networks in the 60 GHz
transmission environment. Receivers designed to process IEEE 802.11ad
waveforms employ very high rate analog-to-digital converters, and
therefore, reducing the receiver sampling rate can be useful. In this
work, we study the problem of low-rate channel estimation over the
IEEE 802.11ad 60 GHz communication link by harnessing sparsity in
the channel impulse response. In particular, we focus on single carrier
modulation and exploit the special structure of the 802.11ad waveform
embedded in the channel estimation field of its single carrier physical
layer frame. We examine various sub-Nyquist sampling methods for this
problem and recover the channel using compressed sensing techniques.
Our numerical experiments show feasibility of our procedures up to
one-seventh of the Nyquist rates with minimal performance deterioration.
Keywords—802.11ad, sub-Nyquist, sparse recovery, Golay sequences,
channel estimation.

I.

Introduction

The pervasive presence of wireless networks in indoor data
communication has led to the demand for very high data throughput
links. In this context, recent advances in millimeter-wave technology
provide multi-gigabit connectivity in the 60 GHz unlicensed
radio-frequency (RF) band with 7-9 GHz of contiguous spectrum
[1]. The IEEE 802.11ad standard has been recently developed to
enable throughput of up to 7 Gbps at 60 GHz for short-range wireless
communication [2]. The diminished range is due to heavy attenuation
through physical barriers at 60 GHz. However, in an indoor wireless
setting, such a high-speed link is ideal to enable applications which
require huge data rates such as transmission of uncompressed high
definition (HD) video [3], in-room gaming [4], intra-large-vehicle
communications [5, 6] and indoor positioning systems [7].
The IEEE 802.11ad transmits both single carrier (SC) and
orthogonal frequency-division multiplexing (OFDM) modulation
schemes at operating frequencies of 1.76 GHz and 2.64 GHz,
respectively. The SC physical layer (SCPHY) frame of 802.11ad
protocol encapsulates Golay complementary sequences [8, 9]. These
sequences have the property of perfect aperiodic autocorrelation (zero
sidelobes) making them useful for communication channel estimation
[10, 11] and radar remote sensing [12, 13]. In this paper, we focus
on the channel estimation problem using 802.11ad SCPHY frames.
Environments such as indoor wireless networks and vehicular
communications are highly variable with typical channel coherence
times of nanoseconds. As a result, 802.11ad chip sequences
are of extremely short duration requiring the receivers to use
analog-to-digital converters (ADCs) that operate at very high (∼GHz)
sampling rates. It is desirable to avoid these high-rate ADCs for
the purposes of saving power, circuit area and cost. In this context,
several methods have been proposed recently to address the problem
The authors are with the Andrew and Erna Viterbi Faculty of Electrical
Engineering, Technion - Israel Institute of Technology, Haifa, Israel. E-mail:
{mishra, yonina}@ee.technion.ac.il. This project is funded by the European
Union’s Horizon 2020 research and innovation programme under grant
agreement ERC-BNYQ. K.V.M. acknowledges partial support via Andrew and
Erna Finci Viterbi Fellowship and Lady Davis Fellowship.

of channel estimation for ultra wideband wireless systems using fewer
samples than obtained at the Nyquist rate (see [14] and references
therein). Most exploit the fact that the channel is sparse facilitating
the use of compressed sensing (CS) techniques [15, 16].
However, sub-Nyquist channel estimation for IEEE 802.11ad
remains relatively unexamined in prior studies. In this paper, we
develop an 802.11ad-based channel estimation procedure that operates
at up to one-seventh of the Nyquist rate without significant loss of
performance. Although a comprehensive statistical characterization of
the 802.11ad channel model is still ongoing, it has been found that it
is similar to the well-documented IEEE 802.15.3c channel model [17,
18], and can be assumed to be sparse. Our recovery algorithm uses
the Xampling framework where Fourier coefficients of the received
signal are acquired from their low-rate samples [16].
Among prior studies, the closest to our work is [14] where analog
sub-Nyquist sampling methods are developed for channel estimation.
The transmit waveform in [14] is an Ipatov sequence that has the
property of perfect periodic autocorrelation. For the IEEE 802.11ad
link, [19, 20] explore channel estimation at Nyquist sampling
rates. Recently, the special structure of 802.11ad SCPHY frame
has also been exploited for parameter estimation in an automotive
radar setup that operates at Nyquist rates and uses Wireless Local
Area Network (WLAN) for vehicle-to-vehicle communications [21].
CS-based channel estimation for 802.11ad link is investigated in
[18] where it is assumed that low rate measurements are already
available. Compared to these studies, our method directly deals with
the sub-Nyquist sampling of the IEEE 802.11ad waveform, leverages
the special construction of the IEEE 802.11ad SCPHY frame, and
suggests feasible modifications in 802.11ad hardware for possible
implementations. As in [14], we observe that foldable sampling leads
to the best performance over all hardware realizable sub-Nyquist
methods examined.
In the next section, we describe the special structure of the
IEEE 802.11ad SCPHY preamble including properties of Golay
complementary sequence pairs. Then, we explain our system model
and conventional Nyquist-based channel estimation. We present
sub-Nyquist channel estimation in Section IV. Various sampling
strategies to obtain few Fourier coefficients are discussed in Section V
where we validate our model and methods via numerical experiments.
We provide concluding remarks in Section VI.
II.

System Model and Nyquist Processing

The IEEE 802.11ad protocol frame has a complex structure. We
rely only on its SCPHY preamble for channel estimation.
A. Transmitter
The IEEE 802.11ad SCPHY frame consists of a short training
field (STF), a channel estimation field (CEF), header, data and
beamforming training field (Fig. 1). The STF and CEF together form
the SCPHY preamble. Within the CEF lie two 512-point sequences
Gu512 [n] and Gv512 [n]. Each of them contains a Golay complementary
pair of length 256, {Gau256 , Gbu256 } and {Gav256 , Gbv256 }, respectively.

hCH (t) =

L
X

αl .δ(t − τl ),

(8)

l=1

where αl ∈ C and τl are the path losses and delays, respectively. The
received signal is then given by
xR (t) = (xT ∗ hCH )(t) + η(t) = (sT ∗ hT ∗ hCH )(t) + η(t),
Fig. 1. Composition of IEEE 802.11ad SCPHY frame.

A Golay complementary pair consists of two sequences GaN and
GbN each of the same length N with entries ±1, such that their
aperiodic autocorrelation functions have sidelobes equal in magnitude
but opposite in sign. The sum of the two autocorrelations, therefore,
has a peak of 2N and a sidelobe level of zero:
GaN [n] ∗ GaN [−n] + GbN [n] ∗ GbN [−n] = 2Nδ[n],

(1)

where ∗ denotes linear convolution. To observe this property in the
frequency domain, consider FaN and FbN as the Discrete Fourier
Transforms (DFTs) of the sequences GaN and GbN , respectively.
Then, the sum of the spectral densities of the two sequences is flat,
i.e., for each DFT index k,
|FaN [k]|2 + |FbN [k]|2 = 2N.

(2)

The convolution in (1) is linear while the DFT has the circular
convolution property. Therefore, for these relations to hold, we must
choose a DFT of length NDFT ≥ 2N − 1. Other perfect autocorrelation
sequences [22] and complementary pairs [23] also exist and have been
used for radar and communication applications.
The 802.11ad CEF transmit signal is a concatenated sequence
rCE S C [n] = Gu512 [n] + Gv512 [n − 512], n = 0, 1, · · · , 1151.

(3)

Here, the sequences Gu512 [n] and Gv512 [n] are defined for 0 ≤
n ≤ 511; for other values of n, they are set to zero. Typical
802.11ad transmitters employ DACs that operate at 2.64 GHz to also
enable OFDM signals. Therefore, the discrete-time signal rCES C [n] is
upsampled by a factor of 1.5 to obtain a sequence sT [n] just before
sending it to a digital-to-analog-converter (DAC) module (Fig. 2):
r1CES C [n] =

(
rCES C [m], n = 3m,
0, otherwise,

sT [n] = r1CES C [2n].

(4)
(5)

The corresponding upsampled constituting sequences of sT [n] are
denoted as Gu512up [n], Gv512up [n], and so on. We represent the
discrete-time sequence sT [n] as a weighted sum of Dirac impulses
sT (t) =

1535
X

where η(t) is additive white Gaussian noise (AWGN). Prior to
sampling at the receiver, an anti-aliasing continuous-time filter hR (t)
(usually a milled analog filter at 60 GHz) is used to remove
high-frequency noise from the signal xR (t) and also to attenuate
adjacent channels. The signal processed by the filter hR (t) is
sR (t) = (xR ∗ hR )(t) = (sT ∗ hT ∗ hCH ∗ hR )(t) + (η ∗ hR )(t)
= (sT ∗ hT R ∗ hCH )(t) + z(t),

C. Channel estimation
In a Nyquist receiver, the received signal samples are
downsampled by a factor of 1.5. Since the constant spectrum property
(2) of Golay sequences is not affected due to upsampling, here we
skip the downsampling stage. The samples sR [n] are passed through
four correlators, each corresponding to one of the Golay sequences of
length 256. For example, if we assume that the overall analog filter
has a flat unity frequency response over the band of interest, then the
signal trail of one of the correlators is
ĥ1au [n] = sR [n] ∗ Gau256 [−n] = (hCH (nT c ) ∗ sT [n]) ∗ Gau256 [−n].

(11)

By correlating the received samples with Gbu256 , Gav256 , and Gbv256 ,
we similarly obtain ĥ1bu [n], ĥ1av [n], and ĥ1bv [n].
Since we intend to exploit the property in (1), we delay the
correlator outputs: ĥ1a [n] = ĥ1au [n],ĥ1b [n] = ĥ1bu [n + 256], ĥ2a [n] =
ĥ2av [n+512], ĥ2b [n] = ĥ2bv [n+768]. The outputs of the complementary
pairs are then summed up. As shown in Fig. 2, this results in two
estimates of the channel, hˆ1 [n] and hˆ2 [n]. For instance, the first
estimate of the channel is
ĥ1 [n] = ĥ1a [n] + ĥ1b [n]
= hCH [n] ∗ (Gau256 [n] ∗ Gau256 [−n] + Gbu256 [n] ∗ Gbu256 [−n])
= 512 · hCH [n].

(12)

1 1
We average to obtain the final estimate ĥ[n] = 512
(ĥ [n] + ĥ2 [n]). In
2 1
the absence of noise, ĥ[n] recovers the true channel perfectly.

Sub-Nyquist Channel Estimation

(6)

m=0

where Fc = 2.64GHz = 1/T c . This signal is passed through the
transmit shaping filter hT (t) (Fig. 2) to obtain:
xT (t) = (sT ∗ hT )(t).

(10)

where hT R (t) = (hT ∗ hR )(t) represents the overall analog filter impulse
response and z(t) = (η∗hR )(t) is the filtered noise. The received signal
sR (t) is sampled at rate Fc = 2.64 GHz to obtain the discrete-time
signal sR [n] = sR (nT c ).

III.
sT [m]δ(t − mT c ),

(9)

(7)

The 802.11ad protocol specifies a spectral mask for the transmit signal
to limit inter-symbol interference [2, section 21.3.2]. We assume that
hT (t) includes a low-pass baseband filter with an equivalent amplitude
characteristic of the spectral mask.
B. Channel and receiver
For the entire duration of the SCPHY packet, the channel is
treated as linear and time-invariant [17]. Assuming the received signal
contains echoes from L paths, the continuous channel response is

In order to reduce the GHz sampling rates of 802.11ad receivers,
we propose sub-Nyquist channel estimation. It has been shown
[24–26] that sparse channels as in (8) can be recovered from low-rate
samples using only K ≥ 2L Fourier coefficients. We use the Xampling
framework to obtain Fourier coefficients from low-rate samples.
A. Fourier series representation
The sequence sT (t) consists of two Golay complementary pairs
and each pair is processed separately. Therefore, hereafter, it suffices
to consider only one pair (say {Gau256up (t), Gbu256up (t)}) for our
discussion. The continuous-time signals Gau256up (t) and Gbu256up (t)
are derived as weighted sums of Diracs from their discrete-time
counterparts. The received signal corresponding to one of the transmit
Golay sequences Gau256up (t) is

Fig. 2.

Nyquist processing for 802.11ad-based channel estimation.

Fig. 3.

Processing flow for 802.11ad protocol channel estimation using sub-Nyquist sampling.

sau (t) = (Gau256up ∗ hT R ∗ hCH )(t) + z(t).

(13)

The channel impulse response is of finite duration and we observe the
signal for time τCH . The Fourier series representation of the received
signal in (13) is
X

sau (t) =

dau [k]e j2πkt/τCH , t ∈ [0, τCH ].

(14)

k∈Z

Here, dau [k] are the Fourier series coefficients:
dau [k] = cau [k] + Z[k],

(15)

where cau [k] and Z[k] are the Fourier series coefficients of the signal
and noise trails. In particular,

cau [k] =

=

1

ZτCH
sau (t)e− j2πkt/τCH dt

τCH
1
τCH

=

τCH

l=0

αl

dau [k] =

1
τCH

=

rau (t − τl )e− j2πk(t−τl )/τCH e− j2πkτl /τCH dt.

(16)

0

=

τCH
1
τCH

1
τCH

|Fau (2πk/τCH )|2 HT R (2πk/τCH )

L−1
X

αl e− j2πkτl /τCH

l=0

+ Z[k]Fau (−2πk/τCH ),

d̃bu [k] =

(20)

1
τCH

|Fbu (2πk/τCH )|2 HT R (2πk/τCH )

L−1
X

αl e− j2πkτl /τCH

l=0

+ Z[k]Fbu (−2πk/τCH ).

−∞

Rau (2πk/τCH )

αl e− j2πkτl /τCH + Z[k]. (19)

l=0

for all k ∈ |κ|. Similarly, for the other part of the complementary pair,

Substituting into (16), we get
1

L−1
X

The time delays are aligned to a grid τl = nl ∆τ where 0 ≤ nl ≤ N and
N = τCH /∆τ is an integer. In Nyquist processing, all N complex
Fourier coefficients are obtained. For sub-Nyquist processing, the
received signal is sampled at low rates. Corresponding to these
low-rate samples, let κ denote a subset of N complex Fourier
τCH
. In Fig. 3, this
coefficients such that κ = {k1 , · · · , k|κ| } and k|κ| < fT2π
process corresponds to the “Sub-Nyquist sampler” block and yields
|κ| = K Fourier coefficients.

=rau (t)
ZτCH

Suppose the continuous-time Fourier Transforms (CTFTs) of
Gau512up (t) and hT R (t) are Fau (ω) and HT R (ω), respectively. Then,
R∞
the CTFT of rau (t) is Rau (ω) =
rau (t)e− jωt dt = Fau (ω)HT R (ω).

cau [k] =

Fau (2πk/τCH )HT R (2πk/τCH )

d̃au [k] = dau [k]Fau (−2πk/τCH )

(Gau256up ∗ hT R ∗hCH )(t)e− j2πkt/τCH dt
|
{z
}

L−1
1 X

Substituting (17) into (15) produces

When the sampled signal is correlated with the Golay sequence
Gau256up , the Fourier coefficients of the resulting signal are

0
ZτCH

0

Following the aforementioned procedure, we similarly obtain dbu [k],
dav [k], and dbv [k] corresponding to continuous-time Golay sequence
signals Gbu256 (t), Gav256 (t), and Gbv256 (t), respectively.
B. Sampling and matched filtering

(21)

The sum of the two complementary Fourier coefficients yields
L−1
X

αl e− j2πkτl /τCH

l=0

Fau (2πk/τCH )HT R (2πk/τCH )

L−1
X

αl e− j2πkτl /τCH .

(17)

d̃u [k] = d̃au [k] + d̃bu [k]
1
=
{|Fau (2πk/τCH )|2 + |Fbu (2πk/τCH )|2 } HT R (2πk/τCH )
{z
}
τCH |
Cu

l=0

If the low-pass filter hT R (t) has a cut-off frequency fT (say), then the
Fourier series coefficients of the filtered signal above fT vanish:

·

L−1
X

αl e− j2πkτl /τCH + Z̃[k]

l=0
L−1

cau [k] = 0, ∀ k ≥

fT τCH
.
2π

(18)

=

X
Cu
HT R (2πk/τCH )
αl e− j2πkτl /τCH + Z̃[k],
τCH
l=0

(22)

where Cu = 512 (i.e. twice the Golay sequence length), and Z̃[k] =
Z[k]Fau (−2πk/τCH ) + Z[k]Fbu (−2πk/τCH ) is the corresponding noise
trail. We assume that only AWGN enters the receiver so that the
filtered, sampled noise is white discrete noise. If the power density
of white noise z(t) is N0 /2 and the bandwidth of the anti-aliasing filter
hR (t) is B, then the power of discrete white noise is N0 B/2. For both
Nyquist and sub-Nyquist processing, the response HT R (2πk/τCH ) is
typically known via prior calibration.
C. Sparse recovery
We can express the coefficients in (22) in matrix form as
d̃u =

Cu
HT R By1 + z̃,
τCH

(23)

iT
where d̃u = d̃u [k1 ] · · · d̃u [k|κ| ] , y1 is an L-sparse vector of
length N whose non-zero values are the path losses αl , HT R =

diag HT R (2πk1 /τCH ), · · · , HT R (2πk|κ| /τCH ) is a K×K diagonal matrix,
z̃ is the noise vector, and B is a K × N matrix with {l, m}th element
Blm = e− j2πkl m/N . Equation (23) models a sparse recovery problem
h

minimize ky1 k0
y1

subject to

d̃u −

Cu
HTR By1
τCH

≤ σ1 ,

(24)

2

where σ1 is some positive constant and k·k0 is the number of non-zero
elements of the vector. This problem can be solved using one of the
several compressed sensing recovery algorithms such as orthogonal
matching pursuit (OMP) and `1 minimization [15]. Let the solution of
(24) be ŷ1 . For the other Golay sequence pair, we similarly obtain d̃v
and the corresponding sparse vector ŷ2 for a different set of sampled
coefficients. Our final channel estimate is then ŷ = (ŷ1 + ŷ2 )/2.
Alternatively, if the sampling coefficients are the same for both
Golay pairs, then we first average to obtain d̃ = (d̃u + d̃v )/2 and then
obtain the channel estimate ŷ by solving a single sparse recovery
problem
minimize kyk0
y

subject to

d̃ −

Cu
HTR By
τCH

where σ is some positive constant.
IV.

≤ σ,

(25)

2

Sub-Nyquist Sampling Methods

There are several ways to implement the sub-Nyquist sampler.
For simplicity, consider |κ| = K such that q = N/K is an integer
defining the sampling reduction factor. The 802.11ad single carrier
Nyquist rate is fNyq = Fc = 1/T c = 2.64 GHz. In direct sampling
(Fig. 4a), the signal sR (t) obtained after the anti-aliasing filter hR (t) is
passed through as many analog chains as the number of sub-Nyquist
coefficients K. Each branch is modulated by a complex exponential,
followed by integration over τCH . This technique gives most flexibility
in choosing the Fourier coefficients, but is expensive in terms of
hardware.
Another approach is to limit the bandwidth of hR (t) such
that only the lowest K frequencies are free of aliasing (Fig. 4b).
Then, we sample these lowest K frequencies with a (Nyquist)
low-rate ADC. This method requires modification in 60 GHz
hardware so that the analog receive filter h̃R (t) has reduced passband.
In multiband sampling, M disjoint randomly-chosen groups of
consecutive Fourier coefficients are obtained such that the total
sampled coefficients are still K. This translates to splitting the
signal across M branches, passing the downconverted signal through

reduced-bandwidth anti-aliasing filters h̃R (t), and then performing
Nyquist sampling on each band with a low-rate ADC (Fig. 4c).
The previous three methods alter only the receiver. We now
discuss two sampling procedures where both transmitter and receiver
are modified. In foldable sampling [14], the received signal is sampled
below the Nyquist rate resulting in folding over or aliasing of the
signal spectrum that exceeds the Nyquist sampling rate of the ADC.
Aliasing is avoided by notching out those coefficients from the
transmit signal that will alias over the low-frequency part of the
sampled signal. The entire signal consisting of N DFT taps can be
assumed to be divided into q groups of K coefficients. If we subsample
this signal at the rate fNyq /q, then K coefficients of each group will
alias over each other in the sampled signal. But if q − 1 coefficients
at each tap are notched out in the transmit signal, then aliasing is
circumvented. We randomly choose each of the K coefficients from
any of the q groups and notch out the remaining q − 1 coefficients
in the transmit signal. The notch filter is a digital filter (Fig. 4d),
and can be realized by zeroing out the undesired FFT coefficients of
the signal sT [n]. Let 1κ be a vector that is one only at coefficients
in the set κ and zero everywhere. Then, the transmit signal s˜T [n]
is obtained by s˜T [n] = IFFT N {FFT N {sT [n]} · 1κ }, where ‘·0 denotes
entry-wise product. This has more flexibility in randomizing the
selected coefficients over the full signal spectrum. Changes in the
60 GHz filter hR (t) are also not necessary. This sampling may suffer
from SNR loss due to aliasing of out-band noise; fortunately, the
stopband attenuation of 60 GHz filters is sufficiently high to avert
this. The foldable multiband sampling (Fig. 4e), is similar to foldable
sampling except that instead of selecting each coefficient randomly
from q groups, we choose M sets of consecutive coefficients from
each group such that the total number of coefficients is still K. The
corresponding transmitter uses digital bandpass filters.
We verified our procedures through numerical experiments. We
used channel impulse responses of length 192 with a sparsity level
K = 10 and a flat response for the filter HT R as is typical in
802.11ad scenarios. We used four bands for multiband and multiband
foldable sampling. The effective SNR for foldable methods after
sampling is kept the same as the Nyquist in these experiments
showing an improvement in RMSE at low SNRs; here, the noise in
unsampled frequencies is completely filtered out. We characterized the
recovery
performance through the root mean square error: RMSE =
q
PN−1
1/N n=0 |hCH [n] − ŷ[n]|2 . Fig. 5 shows the estimation performance
for different algorithms, reduction factors (q = 2 and 7) and SNR
levels. We observe that the foldable and direct sampling are closest
to the Nyquist sampling in performance especially at low SNRs and
high reduction factor q. At high SNRs, Nyquist sampling shows
superior performance over all other methods. The performance of
average-first-OMP-next is slightly better than averaging the individual
optimized solutions.
V. Summary
We presented a sub-Nyquist channel estimation procedure for
the IEEE 802.11ad link. Among various low-rate sampling options,
the foldable sampling method offers minimum degradation in RMSE
and can be realized using a digitally programmable transmitter.
Other methods necessitate changing the analog receive filter which
is technologically difficult at 60 GHz. With some tolerance in
performance loss, sub-Nyquist sampling methods can complement
the promise of very high throughput in IEEE 802.11ad links.
Acknowledgements
We would like to thank Alecsander Eitan of Qualcomm Israel,
Tamir Bendory and Kfir Cohen for insightful conversations.

Fig. 4.

Sub-Nyquist sampling methods: (a) direct sampling (b) low frequencies only (c) multiband sampling (d) foldable sampling (e) foldable multiband.

Fig. 5.

Performance of average-first-OMP-next algorithm with reduction factor (a) q = 2 (b) q = 7. OMP-first-average-next with (c) q = 2 (d) q = 7.

References
[1]

[2]
[3]
[4]
[5]

[6]
[7]
[8]
[9]
[10]
[11]
[12]
[13]

T. S. Rappaport, S. Sun, R. Mayzus, H. Zhao, Y. Azar, K. Wang, G. N.
Wong, J. K. Schulz, M. Samimi, and F. Gutierrez, “Millimeter wave
mobile communications for 5G cellular: It will work!” IEEE Access,
vol. 1, pp. 335–349, 2013.
“Wireless LAN Medium Access Control (MAC) and Physical Layer
(PHY) Specifications Amendment 3: Enhancements for Very High
Throughput in the 60 GHz Band,” IEEE Std. 802.11ad, 2012.
R. C. Daniels and R. W. Heath Jr., “60 GHz wireless communications:
Emerging requirements and design recommendations,” IEEE Veh.
Technol. Mag, vol. 2, no. 3, 2007.
“Wireless LAN at 60 GHz IEEE 802.11ad explained,” Agilent
Technologies Application Note, 2013.
J. Choi, V. Va, N. Gonzalez-Prelcic, R. Daniels, C. R. Bhat, and
R. W. Heath, “Millimeter-wave vehicular communication to support
massive automotive sensing,” IEEE Commun. Mag., vol. 54, no. 12,
pp. 160–167, 2016.
N. González-Prelcic, R. Méndez-Rial, and R. W. Heath Jr., “Radar aided
beam alignment in mmwave V2I communications supporting antenna
diversity,” in IEEE Inf. Theory and Appl. Workshop, 2016.
C. Chen, Y. Han, Y. Chen, and K. R. Liu, “Indoor global positioning
system with centimeter accuracy using Wi-Fi [Applications Corner],”
IEEE Signal Process. Mag., vol. 33, no. 6, pp. 128–134, 2016.
M. J. E. Golay, “Complementary series,” IRE Trans. on Inf. Theory,
vol. 7, no. 2, pp. 82–87, 1961.
——, “Note on complementary series,” Proc. IRE, vol. 50, p. 84, 1962.
J. Mo, P. Schniter, and R. W. Heath Jr, “Channel estimation in
broadband millimeter wave mimo systems with few-bit adcs,” arXiv
preprint arXiv:1610.02735, 2016.
T. Qureshi, M. Zoltowski, and R. Calderbank, “MIMO-OFDM channel
estimation using Golay complementary sequences,” in Int. Waveform
Diversity Des. Conf., 2009, pp. 253–257.
E. C. Farnett and G. H. Stevens, “Pulse compression radar,” in Radar
Handbook, 2nd ed., M. I. Skolnik, Ed. McGraw-Hill Professional,
1990.
Y. Chi, A. Pezeshki, and R. Calderbank, “Complementary waveforms
for sidelobe suppression and radar polarimetry,” in Principles of
Waveform Diversity and Design, M. Wicks, E. Makole, S. Blunt,
R. Schneible, and V. Amuso, Eds. SciTech Publishing, 2011.

[14]

[15]
[16]
[17]
[18]

[19]

[20]
[21]
[22]
[23]
[24]
[25]
[26]

K. M. Cohen, C. Attias, B. Farbman, I. Tselniker, and Y. C. Eldar,
“Channel estimation in UWB channels using compressed sensing,” in
IEEE Int. Conf. on Acoustics, Speech and Signal Process., 2014, pp.
1966–1970.
Y. C. Eldar and G. Kutyniok, Compressed Sensing: Theory and
Applications. Cambridge University Press, 2012.
Y. C. Eldar, Sampling Theory: Beyond Bandlimited Systems.
Cambridge University Press, 2015.
T. S. Rappaport, R. W. Heath Jr., R. C. Daniels, and J. N. Murdock,
Millimeter Wave Wireless Communications. Prentice Hall, 2014.
R. T. Suryaprakash, M. Pajovic, K. J. Kim, and P. Orlik, “Millimeter
wave communications channel estimation via bayesian group sparse
recovery,” in IEEE Int. Conf. on Acoustics, Speech and Signal Process.,
2016, pp. 3406–3410.
W.-C. Liu, T.-C. Wei, Y.-S. Huang, C.-D. Chan, and S.-J. Jou,
“All-digital synchronization for SC/OFDM mode of IEEE 802.15. 3c
and IEEE 802.11 ad,” IEEE Trans. Circuits Syst. I, Reg. Papers, vol. 62,
no. 2, pp. 545–553, 2015.
B. Ye and Z. Zhang, “Improved pilot design and channel estimation for
60 Ghz OFDM based on IEEE 802.11ad,” in IEEE Wireless Commun.
Netw. Conf., 2013, pp. 4129–4133.
P. Kumari, N. González-Prelcic, and R. W. Heath Jr., “Investigating
the IEEE 802.11ad standard for millimeter wave automotive radar,” in
Vehicular Technology Conf., 2015, pp. 1–5.
H. D. Luke, “Sequences and arrays with perfect periodic correlation,”
IEEE Trans. Aerosp. Electron. Syst., vol. 24, no. 3, pp. 287–294, 1988.
E. E. Hollis, “Another type of complementary sequence,” IEEE Trans.
Aerosp. Electron. Syst., vol. 5, pp. 916–920, 1975.
M. Vetterli, P. Marziliano, and T. Blu, “Sampling signals with finite
rate of innovation,” IEEE Trans. Signal Process., vol. 50, no. 6, pp.
1417–1428, 2002.
E. Baransky, G. Itzhak, I. Shmuel, N. Wagner, E. Shoshan, and Y. C.
Eldar, “A sub-Nyquist radar prototype: Hardware and applications,”
IEEE Trans. Aerosp. Electron. Syst., vol. 50, no. 2, pp. 809–822, 2014.
K. V. Mishra, E. Shoshan, M. Namer, M. Meltsin, D. Cohen,
R. Madmoni, S. Dror, R. Ifraimov, and Y. C. Eldar, “Cognitive
sub-Nyquist hardware prototype of a collocated MIMO radar,” in
Compressed Sens. Theory Appl. Radar, Sonar Remote Sens., 2016, pp.
56–60.

