Hybrid Beamforming with Spatial Modulation
in Multi-user Massive MIMO mmWave
Networks
Merve Yüzgeçcioglu and Eduard Jorswieck

arXiv:1709.04826v1 [cs.IT] 14 Sep 2017

Communications Theory, Communications Laboratory
Dresden University of Technology, D-01062 Dresden, Germany
Email: {merve.yuzgeccioglu, eduard.jorswieck}@tu-dresden.de

Abstract—The cost of radio frequency (RF) chains is
the biggest drawback of massive MIMO millimeter wave
networks. By employing spatial modulation (SM), it is
possible to implement lower number of RF chains than
transmit antennas but still achieve high spectral efficiency.
In this work, we propose a system model of the SM scheme
together with hybrid beamforming at the transmitter
and digital combining at the receiver. In the proposed
model, spatially-modulated bits are mapped onto indices
of antenna arrays. It is shown that the proposed model
achieves approximately 5dB gain over classical multi-user
SM scheme with only 8 transmit antennas at each antenna
array. This gain can be improved further by increasing
the number of transmit antennas at each array without
increasing the number of RF chains.

I. I NTRODUCTION
Recent technological developments also bring a challenge to wireless communications since the number
of network devices increases drastically. Even though
there are several works on achieving high data rates,
the bandwidth shortage of today’s commercial networks
also limits the capacity of the channel. In order to
overcome this bottleneck, millimeter wave (mmWave)
spectrum became tempting due to the availability of large
bandwidths. However, the signal in mmWave spectrum
experiences a severe propagation loss and the resulting
channel is poorly scattered. In spite of the unfavorable
characteristics of the channel, beamforming techniques
can be employed at massive multiple-input multipleoutput (mMIMO) networks to direct the beam with high
array gain. Thanks to the small wave length at high
frequencies, it is possible to pack a large amount of
antennas in small areas and asymptotically achieve the
capacity of the channel [1]. However, packing high number of antennas comes with a price: energy consumption.
It is not practical to implement dedicated RF chains
for each antenna at a transceiver due to the power and
space restrictions. Nevertheless, by employing hybrid
beamforming techniques, it is possible to benefit from
The work of Merve Yüzgeçcioglu has received partly funding
from the European Union’s Horizon 2020 research and innovation
programme under the Marie Sklodowzka-Curie grant agreement No
641985
978-1-5386-3531-5/17/$31.00 c 2017 European Union

the gain of a large amount of antennas while consuming
less energy by using lower number of RF chains than
antennas.
There are several works on optimal precoding in hybrid beamforming [2]–[8]. In [4], a method to construct a
near-optimal hybrid beamforming structure is presented
instead of an exhaustive search to maximize the spectral
efficiency. The work is extended to multi-user scenario in
[5] where users employ only analog combiner while the
base station (BS) has a hybrid structure and constructs
baseband precoder in order to mitigate the inter-userinterference. In [6], together with a two-stage hybrid
beamformer/combiner algorithm, the minimum number
of RF chains to realize the performance with a fullyconnected structure is given. Another approach is presented in [7] that exploits the singular vectors of the
channel to generate analog beamformer and combiner
which provides a lower complexity. Furthermore, an
asymptotic rate expression is given for the proposed
scheme.
Another popular research topic recently is SM that
enables transmitting additional bits to conventional modulated symbol without requiring extra power [9]–[11].
There have been extensive work on SM in different
aspects for MIMO networks. Furthermore, implementation of SM at multi-user MIMO networks is studied
at [12] and a precoding method is proposed in order
to cancel inter-user-interference. In this work, SM is
employed at sub-groups which are dedicated for each
user. Combination of SM and beamforming is introduced in [13]. In this work, analog beamforming is
employed to a generalized SM (GSM) scheme for singleuser scenario at Rician fading channels. A closer look
into the transmitter design of SM together with analog
beamforming is given in [14]. Furthermore, the extension
of SM with hybrid beamforming for a mmWave railway
communication system is studied in [15]. Here, the
receive antenna arrays (AA) at the front and end of
the train are assumed as virtual users served by two
different data streams transmitted from the BS. Analog
beamformer and combiner are employed at the transmitter and receiver, respectively. However, the beamformer
and combiner are not designed specifically for multi-

user spatial modulation system. Rather the design in [3]
is considered. As performance metric, an upper bound
rate for AWGN channel with Gaussian input is considered and maximum-likelihood detection is employed
for decoding. However, it is not clearly stated how the
inter-user-interference is eliminated and data decoding is
performed at each user.
In this paper, implementation of SM to a multi-user
mmWave network is introduced. In this scheme, BS
employs analog beamformer to direct the beam to the
intended user and digital precoder to cancel inter-userinterference while users employ the digital combiners.
The novelty of the paper is summarized as follows
• Analog beamforming is designed in order to maximize the achievable rate of the link. Two different codebooks are considered to choose analog
beamformer vector such that one provides optimal
beamformer while the other allows practical implementation. The quantization error is characterized
between these two codebooks.
• Digital precoder at the transmitter is designed in
order to eliminate inter-user-interference. Digital
combiner at the receiver is designed to successfully
reconstruct transmitted symbol when it is matched
with the correct analog beamformer.
• A low-complexity maximum-likelihood detector is
introduced.
• Analytical expression of the achievable rate of the
system is derived and the tightness of the expression
is shown.
The rest of the paper is organized as follows. In Sec. II,
the system model of the proposed scheme is described.
In Sec. III, design of the analog beamformer and digital
precoder for transmitter side and the digital combiner
for receiver side is explained in detail, and the quantization error between two codebooks is characterized.
In Sec. IV, analytical expression of achievable rate for
the proposed scheme is derived. The performance is
analyzed numerically in Sec. V and finally the paper
is concluded in Sec. VI.
II. S YSTEM M ODEL
The block diagram of the system is shown in Fig. 1.
There, the BS is equipped with NA uniform linear arrays
(ULA), each with NT transmit antennas connected to
NRF RF chains. We assume without loss of generality
that NT is fixed for all AAs and there are no common
transmit antennas between arrays. The BS serves K users
with NR receive antennas and it is assumed that NRF ≥
K. The number of AAs may vary independently from
the number of RF chains at the transmitter.
Data is transmitted to mobile stations (MS) in two
parts: first part is modulated according to spatial modulation principles by choosing one of the AAs which allows
to transmit log2 NA bits and the second part is modulated
according to a conventional M -ary modulation scheme
which allows to transmit log2 M bits. The index of
the active AA conveys information without consuming

additional power. According to incoming bits for each
bit stream, it is possible that one AA is used for more
than one user at a channel use. This is only possible
by the help of analog beamformer and allows to exploit SM for multi-user scenario without any limitation.
After determining the AA indices for every user, the
transmitter constructs the digital precoder and the analog
beamformer. The received signal at i-th user can be
expressed as follows
ri =

K
√ X
ρ
Haj ,i faj ,j paj ,j s + ni .

(1)

j=1

Herein, Haj ,i ∈ CNR ×NT is the L-path channel between
the aj -th AA and i-th user where aj is the selected AA
index to transmit data to j-th user. The channel follows
a geometry-based model shown below [16]
r
L
NT NR X
Haj ,i =
αaj ,il aRaj ,i (θl )aH
Taj ,i (φl ), (2)
L
l=1
q
Here, aRaj ,i (θl ) = N1R[1,ejπ sin θl,. . .,ejπ(NR −1) sin θl]T
q
1
jπ sin φl
,. . .,ejπ(NT −1) sin φl]T
and aTaj ,i (φl ) =
NT [1,e
are the receive and transmit AA responses of the lth path where l = 1, . . . , L. φl and θl are the angle
of departure (AoD) and angle of arrival (AoA) of the
path and drawn from the uniform distribution U(0, 2π].
Finally, αaj ,il , CN (0, 1), represents the channel gain
and the path loss of the path. Throughout the paper,
(.)T and (.)H denote transpose and hermitian operators,
respectively. It is assumed that the antennas are placed on
the array with half-wavelength separation. paj ,j ∈ C1×K
and faj ,j ∈ CNT ×1 are the digital precoder and analog
beamformer vectors for j-th user when aj -th AA array
is chosen for transmission. ρ represents the total transmit
power and with equal power allocation, the transmit
power of i-th user is ρi = ρ/K. The selection procedure
for precoder and beamformer is explained in detail in
Sec. III. s ∈ CK×1 is the data vector contains symbols
for K users and finally ni ∈ CNR ×1 is the noise vector
with elements distributed as CN (0, σ 2 ).
A. Receiver Structure
The receiver structure of the users is shown in Fig. 1.
In this system model, MS has NR receive antennas
connected to the same number of RF chains. The received signal is passed on to the combiner with the
help of RF chains and digital combining is performed
to successfully reconstruct the transmitted symbol. In
order to decode the received symbol, a low complexity maximum-likelihood (ML) detector is employed to
jointly estimate the AA index and M -ary modulated
symbol
2

[âi , ŝi ] = arg min
a,m

1 H
w ri − sm ,
β a,i

(3)

where wa,i ∈ CNR ×1 is the digital combiner from a-th
AA to i-th user, a = 1, . . . , A, β is the normalization

RF chain

Digital
Combiner

RF chain

W1

NR
log2 NA

log2 M
.
.
.
K
.
.
.

log2 NA

RF chain
Digital
Precoder
P

log2 M

.
.
.
NRF
.
.
.

.
.
.
NA
.
.
.

Analog
Beamformer
F

NA

Decoder

â1 , ŝ1

NA

Decoder

âK , ŝK

.
.
.

RF chain

RF chain

Digital
Combiner

RF chain

WK

NR

Fig. 1: Block diagram of hybrid beamforming with spatial modulation in multi-user downlink transmission
coefficient and sm is the m-th symbol from M -ary
constellation diagram. This detector structure allows the
user to decode without having Haj ,i , faj ,j and paj ,j
knowledge but only wa,i vectors and β. The construction
of wa,i is also explained in detail in Sec. III.
III. P RECODER , B EAMFORMER AND C OMBINER
D ESIGN
We consider two different approaches to design the
analog beamformer vector. The first method is choosing the vector among the transmit AA response vectors aT (φ). The second method is using a predefined
beamsteering codebook with quantized angles between
(0, 2π].
Let us define the codebook F = {fn ∈ CNT ×1 :
H
fn fn = 1, n = 1, . . . , N } and name the two different
codebooks as FA for transmit AA response vectors
collection and FB for the beamsteering codebook. For
both cases,
q the individual codewords are constructed
1
T
as f =
NT [1, exp (jπφ), . . . , exp (jπ(NT − 1)φ)]
where φ is the AoD of each path for codebook FA
and quantized angle 2πn
N where n = 0, . . . , N − 1 for
codebook FB .
Furthermore, the analog beamformer f is chosen in
order to maximize the signal-to-noise ratio (SNR) of the
link between the selected AA and the intended user.
For each user, the BS computes the RF beamformer
f1,i , . . . , fNA ,i as follows
fa,i = arg max ||Ha,i fn ||2 ,
∀fn ∈F

(4)

where fa,i is the optimum RF beamformer from the ath AA to i-th user. The resulting analog beamforming
matrix for all users is F = [fa1 ,1 , . . . , faK ,K ].
After determining the optimum RF beamformers for
each AA-user pair, the BS constructs the matrix Hi =
[H1,i f1,i , . . . , HNA ,i fNA ,i ]. Hi is used to calculate the
digital combiner for each user to successfully reconstruct
the transmitted symbol


Wi = w1,i . . . wNA ,i = (H†i )H ,
(5)
where (.)† denotes the pseudo inverse of a matrix. With
this design, it is possible to feedback only Hi matrices
to the users instead of having a training period of the
channel matrices for each AA-user pair. In this case, it is
sufficient for users to have Hi with dimension NA ×NR

instead of full channel matrices H1,i , . . . , HNA ,i with
dimension NA NT × NR .
After determining the optimum RF beamformer and
receive combiner, the transmitter generates the effective
channel in Eq. [6] order to calculate the digital precoder
to eliminate the inter-user-interference
 H

wa1 ,1 Ha1 ,1 fa1 ,1 . . . waH1 ,1 HaK ,1 faK ,K


..
..
..
Hef f = 
.
.
.
.
waHK ,K Ha1 ,K fa1 ,1 . . . waHK ,K HaK ,K faK ,K
(6)
Finally the precoding vector for each user is determined
as P = [pTa1 ,1 . . . pTaK ,K ]T = βH†ef f ,where pak ,k ∈
C1×K is the precoder vector of k-th user when ak -th AA
is chosen by spatially-modulated bits. The coefficient β
satisfies the average power constraint which is calculated
as
(s
)
K
β=E
,
(7)
tr(H†ef f (H†ef f )H )

where the expectation is over channel realizations. Since
the complex channel coefficients are the key parameters
for spatially-modulated bits, instantaneous power constraint is not applicable for such systems [12].
A. Characterization of the Quantization Error
For large scale AAs, the right and left singular vectors
of the channel are the transmit and receive AA response
vectors, respectively [17]. Let us define the optimal
analog beamformer for the channel given in Eq. [2] in
accordance with this statement as follows
fA = arg

max

∀fA,n ∈FA

||HfA,n ||2 ,

(8)

where N = L is the number of paths of the channel
and the total number of elements in the codebook.
H ∈ CNR ×NT is the channel matrix, fA,n ∈ CNT ×1
is the n-th element of the codebook FA . In order to
achieve the performance of the optimal beamformer,
infinite resolution phase shifters should be implemented
at the system. Since this is not possible for practical
systems, beamsteering codebooks with quantized angle
are commonly considered. However, selecting the analog
beamformer among beamsteering codebooks leads to a
quantization error.

Information theoretical methods can be used to find the
achievable rate of the proposed scheme. The mutual
information between transmitted and received symbols
is written as

UB
Exact
LB

2

bits/s/Hz

1.5

I(yi ; xi ) = h(yi ) − h(zi ).
1

Therein, xi is considered as the symbol that contains
the AA index and M -ary modulated symbol to ease
the understanding. zi is the noise term after receive
combining which follows the same distribution with ni .

0.5

0
-10

-5

0

5

10

15

SNR [dB]

Fig. 2: Rate approximation of a user for NA = 2, NT =
8, NR = 1 with BPSK modulation
Proposition 1. Assume that fA and fB are the solution
to Eq. [4] for codebooks FA and FB , respectively. Then,
the quantization error can be characterized as follows
dc,min ≤

max

∀fA ∈G(NT ,1)

d2c (fA , fB ) ≤ dc,max ,

(9)

H
}, dc,max =
where dc,min = 1 − tr{fA fAH fB,n fB,n
−

1

(13)

−

B

c NT −1 2 NT −1 and c is the coefficient of the metric
ball volume defined in [18].
Proof: Let us assume the codebook FB with B bits
resolution phase shifters where codebook size is N =
2B . Eq. [4] can be reformulated as follows
fB , arg
= arg

min

d2c (fA , fB,n )

min

H
1 − tr{fA fAH fB,n fB,n
},

∀fB,n ∈FB
∀fB,n ∈FB

(10)

where dc (., .) is the chordal distance. Eq. [10] is an
example of Grassmannian quantization on the Grassmann manifold G(NT , 1). The codebook FB can be
considered as a Grassmannian subspace sphere-packing
codebook since the elements of FB are distributed over
H
a sphere with radius ||fB,n
fB,n ||2 = 1, ∀n = 1, . . . , N .
Let us define the distance between the optimal beamformer and the selected beamformer from the quantized
codebook as dc (fA , fB ). By using the results in [18]–[20]
we reach the Eq. [9].
It is clearly seen from Eq. [9] that the
quantization error vanishes asymptotically such as
− 1
− B
limB→∞ c NT −1 2 NT −1 = 0.
IV. ACHIEVABLE R ATE
The received signal at i-th user after receive combining is
K

√ X
yi = waHi ,i ρ
Haj ,i faj ,j paj ,j s + waHi ,i ni .

(11)

j=1

Since the interference from the other users is eliminated
by employing transmit precoder, the resulting signal at
the receiver can be reformulated as follows
√
yi = ρi waHi ,i Hai ,i fai ,j si + waHi ,i ni .
(12)

Proposition 2. The received signal yi follows Gaussian
mixture distribution. There is no closed-form solution
for the entropy of a Gaussian mixture but using a tight
approximation the achievable rate can be bounded as
hl (yi ) − h(zi ) ≤ Ri ≤ hu (yi ) − h(zi ),

(14)

where hl (yi ) = (γ + αN,N 0 ) log2 e + log2 σ, hu (yi ) =
(γ + βN,N 0 ) log2 e + log2 σ and h(zi ) = log2 (πeσ 2 ).
Proof: We begin the proof by deriving the probability density function of the received signal that can be
written as follows
fYi(yi) =

NA
M X
X
fYi(yi |a = j,s = sm)P(a = j,s = sm). (15)
m=1 j=1

Since the selection of the AA indices and M -ary symbols
are independent events, joint probability in Eq. [15] can
be separated as P (a = j, s = sm ) = P (a = j)P (s =
1
sm ) = N1A M
. Furthermore, it is clearly seen that
fYi (yi |a = j, s = sm ) follows complex Gaussian distri√
H
bution with different mean values ρi wj,i
Hj,i fj,i sm according to the j and m indices. The resulting probability
density function follows Gaussian mixture distribution
fYi (yi ) =

M NA
1 XX
M NA m=1 j=1
(
) (16)
√
H
|yi − ρi wj,i
Hj,i fj,i sm |2
1
exp −
,
π
σ2

where Hj,i , fj,i and wj,i are the channel matrix, analog
beamformer and receive combiner for j, i AA-user pair,
respectively. Herein, sm denotes the m-th symbol from
M -ary constellation diagram and ρi is the transmit power
for i-th user. By using the tight approximation derived
for Gaussian mixtures in [21], entropy of yi can be
bounded as (γ + αN,N 0 ) log2 e + log2 σ ≤ h(yi ) ≤
(γ + βN,N 0 ) log2 e + log2 σ. Due to space constraints
the detailed explanation on the parameters γ, αN,N 0 and
βN,N 0 are not given here and can be found in [21].
The tightness of the introduced approximation is
shown in Fig. 2. In this figure, the lower (LB) and upper
(UB) bounds are calculated by using the approximation
and the exact rate is found Rby calculating the exact
∞
integral for entropy h(y) = − −∞ fY (y) log fY (y)dy.

100

100

FB - B = 6
FB - B = 7

10-1

FB - B = 8

10-1

FB - B = 9

BER

BER

10-2
-2

10

10-3
10-3

SM - NT = 4
HBFSM - NA
HBFSM - NA
HBFSM - NA
HBFSM - NA

10-4

10-4

10-5
0

5

10

15

20

25

30

SNR [dB]

0

5

= 4,
= 4,
= 4,
= 4,

10

NT
NT
NT
NT

= 8, FB
= 8, FA
= 64, FB
= 64, FA

15

20

25

30

SNR [dB]

Fig. 3: BER performance for different codebook sizes
(K = 2, NA = 4, NT = 8, NR = 1, 4 − QAM )

Fig. 4: Comparison of classical SM and HBF-SM (K =
2, NR = 1, B = 8, 4 − QAM )

V. N UMERICAL R ESULTS

floor when codebook FB is used, the error performance
is better than the classical SM scheme up to high SNR
regime. When the optimal codebook is used, HBF-SM
outperforms classical SM regardless the SNR value. It
is possible to have 5dB gain over classical SM scheme
with only 8 transmit antennas at an AA. If we consider
64 transmit antennas at each array, the gain will increase
up to 12dB. Although there is an intersection point
between HBFM-SM and classical SM, by increasing the
codebook size and transmit antenna numbers at an array,
the performance of HBF-SM can be further increased.
Note that, each AA is served by a single RF chain.
Hence, increasing the number of transmit antennas at an
array does not come with a high cost while significantly
increasing the performance. As it is seen from the
simulation results, the number of antennas at an array is
one of the key parameters that affects the performance.
The larger the number of transmit antennas the better the
beam can be directed to intended direction.

In this section, the performance of the hybrid beamforming with spatial modulation (HBF-SM) is analyzed
numerically. Firstly, the effect of the quantization error
on the uncoded bit error rate (BER) performance is
investigated for different codebook sizes. Then, it is
compared with the conventional SM scheme for multiuser scenario [12].
In Fig. 3, the BER performance of the system is
shown for the cases that the analog beamformer vector
is chosen among the transmit AA response vectors
(FA ) and the beamsteering codebook (FB ). It is seen
that the quantization error of the codebook FB leads
to an error floor for the low codebook sizes. Since
the digital combiner at the receiver side is generated
by using the analog beamforming vector also, when
the beamforming vector is erroneous, it cumulatively
increases the resulting error probability. As it is expected,
when the analog beamformer vector is chosen from
codebook FA , the performance is better and there is
no error floor in the system. Although the codebook
FA provides the optimal beamformer, there is no phase
shifter with infinite resolution. Hence, the codebook
FB with B-bit resolution would be used in practical
systems. When the phase shifters are changed from 6bits to 9-bits the performance is increased significantly.
As it is shown in Sec. III-A the performance of the
beamsteering codebook will eventually converge to the
optimal beamformer.
In Fig. 4, the BER performance of the HBF-SM is
compared with the conventional SM scheme in multiuser scenario. Average power constraint is employed for
both schemes. The considered SM scheme has NT = 4
transmit antennas. For HBF-SM scheme, NA = 4 antenna arrays are considered with NT = 8 and NT = 64
transmit antennas at each array, respectively. In SM
scheme, total log2 NT = 2 bits are transmitted at a
channel use by the selected transmit antenna index.
In HBFM-SM scheme, also log2 NA = 2 bits are
transmitted by the selected AA index. It is seen from
the simulation results, although HBF-SM has an error

VI. C ONCLUSION
In this work, we have introduced a system model
of HBF-SM scheme. Novel design of the beamformer
and digital combiners for the proposed model have been
derived. The effect of the codebook selection for analog
beamformer on performance has been shown and the
quantization error has been studied. Additionally, a tight
approximation for the rate has been introduced. Uncoded
BER performance of the system has been shown for
different parameters and also compared with a classical
SM implementation for multi-user networks. Simulation
results showed that even with low number of transmit
antennas at the AAs, performance of classical SM is
outperformed. For future work, the receiver design can
be improved by employing hybrid beamforming that also
allows reducing the number of RF chains. Since the
most important objective of the SM schemes is providing
higher spectral efficiency while consuming lower energy
compared to conventional schemes, optimizing the energy efficiency of the system is an important direction
for future research.

R EFERENCES
[1] T. L. Marzetta, “Noncooperative cellular wireless with unlimited
numbers of base station antennas,” IEEE Transactions on Wireless Communications, vol. 9, no. 11, pp. 3590–3600, Nov. 2010.
[2] Z. Pi, “Optimal transmitter beamforming with per-antenna power
constraints,” in IEEE Signal Processing for Communications
Symposium, Jun. 2012, pp. 3779–3784.
[3] T. Kim, J. Park, J.-Y. Seol, S. Jeong, J. Cho, and W. Roh, “Tens
of Gbps support with mmwave beamforming systems for next
generation communications,” in IEEE Global Communications
Conference, Dec. 2013, pp. 3685–3690.
[4] O. E. Ayach, S. Rajagopal, S. Abu-Surra, Z. Pi, and R. W. Heath
Jr., “Spatially sparse precoding in millimeter wave MIMO systems,” IEEE Transactions on Wireless Communications, vol. 13,
no. 3, pp. 1499–1513, Mar. 2014.
[5] A. Alkhateeb, G. Leus, and R. W. Heath Jr., “Limited feedback
hybrid precoding for multi-user millimeter wave systems,” IEEE
Transactions on Wireless Communications, vol. 14, no. 11, pp.
6481–6494, Nov. 2015.
[6] F. Sohrabi and W. Yu, “Hybrid digital and analog beamforming
design for large-scale MIMO systems,” in IEEE International
Conference on Acoustics, Speech and Signal Processing, Apr.
2015, pp. 2929–2933.
[7] S. Payami, M. Ghoraishi, and M. Dianati, “Hybrid beamforming
for large antenna arrays with phase shifter selection,” IEEE
Transactions on Wireless Communications, vol. 15, no. 11, pp.
7258–7271, Nov. 2016.
[8] T. E. Bogale, L. B. Le, A. Haghighat, and L. Vandendorpe, “On
the number of RF chains and phase shifters, and scheduling design with hybrid analog-digital beamforming,” IEEE Transactions
on Wireless Communications, vol. 15, no. 5, pp. 3311–3326, May
2015.
[9] R. Y. Mesleh, H. Haas, S. Sinanovic, C. Wook, and S. Yun,
“Spatial modulation,” IEEE Transactions on Vehicular Technology, vol. 57, no. 4, pp. 2228–2241, Jul. 2008.
[10] M. D. Renzo, H. Haas, A. Ghrayeb, S. Sugiura, and L. Hanzo,
“Spatial modulation for generalized MIMO: Challenges, opportunities and implementation,” IEEE in Proceedings, vol. 102, no. 1,
pp. 56–103, Jan. 2014.
[11] M.-C. Lee, W.-H. Chung, and T.-S. Lee, “Generalized precoder
design formulation and iterative algorithm for spatial modulation
in MIMO systems with CSIT,” IEEE Transactions on Communications, vol. 63, no. 4, pp. 1230–1244, Apr. 2015.

[12] S. Narayanan, M. J. Chaudhry, A. Stavridis, M. D. Renzo,
F. Graziosi, and H. Haas, “Multi-user spatial modulation MIMO,”
in IEEE Wireless Communications and Networking Conference,
Apr. 2014, pp. 671–676.
[13] N. Ishikawa, R. Rajashekar, S. Sugiura, and L. Hanzo,
“Generalized-spatial-modulation-based
reduced-RF-chain
millimeter-wave communications,” IEEE Transactions on
Vehicular Technology, vol. 66, no. 1, pp. 879–883, Jan. 2017.
[14] M.-C. Lee and W.-H. Chun, “Transmitter design for analog
beamforming aided spatial modulation in millimeter wave MIMO
systems,” in IEEE International Symposium on Personal, Indoor
and Mobile Radio Communications: Fundamentals and PHY,
Sep. 2016, pp. 1–6.
[15] Y. Cui, X. Fang, and L. Yan, “Hybrid spatial modulation beamforming for mmwave railway communication systems,” IEEE
Transactions on Vehicular Technology, vol. 65, no. 12, pp. 9597–
9606, Dec. 2016.
[16] A. A. M. Saleh and R. A. Valenzuela, “A statistical model for
indoor multipath propogation,” IEEE Journal on Selected Areas
in Communications, vol. SAC-5, no. 2, pp. 128–137, Feb. 1987.
[17] O. E. Ayach, R. W. Heath Jr., S. Abu-Surra, S. Rajagopal, and
Z. Pi, “The capacity optimality of beamsteering in large millimeter wave MIMO systems,” in IEEE International Workshop on
Signal Processing Advances in Wireless Communications, Jun.
2012, pp. 100 – 104.
[18] R. T. Krishnamachari and M. K. Varanasi, “On the geometry
and quantization of manifolds of positive semi-definite matrices,”
IEEE Transactions on Signal Processing, vol. 61, no. 18, pp.
4587–4599, Sep. 2013.
[19] W. Dai, Y. Liu, and B. Rider, “Quantization bounds on Grassmann manifolds and applications to MIMO communications,”
IEEE Transactions on Information Theory, vol. 54, no. 3, pp.
1108–1123, Mar. 2008.
[20] P. Cao, A. Zappone, and E. A. Jorswieck, “Grouping-based
interference alignment with IA-cell assignment in multi-cell
MIMO MAC under limited feedback,” IEEE Transactions on
Signal Processing, vol. 64, no. 5, pp. 1336–1351, Mar. 2016.
[21] K. Moshksar and A. K. Khandani, “Arbitrarily tight bounds on
differential entropy of Gaussian mixtures,” IEEE Transactions on
Information Theory, vol. 62, no. 6, pp. 3340–3354, Jun. 2016.

