1

Robust Secure Transmission of Using
Main-Lobe-Integration Based Leakage Beaforming
in Directional Modulation MU-MIMO Systems

arXiv:1704.07161v2 [cs.IT] 1 Aug 2017

Feng Shu, Wei Zhu, Xiangwei Zhou, Jun Li, and Jinhui Lu .

Abstract—In the paper, we make an investigation of robust
beamforming for secure directional modulation in the multi-user
multiple-input and multiple output (MU-MIMO) systems in the
presence of direction angle measurement errors. When statistical
knowledge of direction angle measurement errors is unavailable,
a novel robust beamforming scheme of combining main-lobeintegration (MLI) and leakage is proposed to simultaneously
transmit multiple different independent parallel confidential message streams to the associated multiple distinct desired users. The
proposed scheme includes two steps: designing the beamforming
vectors of the useful confidential messages and constructing
artificial noise (AN) projection matrix. Here, in its first step,
the beamforming vectors for the useful confidential messages of
desired user k are given by minimizing the useful signal power
leakage from main-lobe of desired user k to the sum of main-lobes
of the remaining desired directions plus main-lobes of all eavesdropper directions. In its second step, the AN projection matrix is
constructed by simultaneously maximizing the AN power leakage
to all eavesdropper directions such that all eavesdroppers are
disrupted seriously, where AN is viewed by the transmitter as a
useful signal for eavedroppers. Due to independent beamforming
vectors for different desired users, a more secure transmission
is achieved. Compared with conventional non-robust methods,
the proposed method can provide a significant improvement in
bit error rate along the desired directions and secrecy-sumrate towards multiple desired users without needing statistical
property or distribution of angle measurement errors.
Index Terms—Robust, directional modulation, leakage, mainlobe-integration, multi-user MIMO, secrecy-sum-rate, artificial
noise.

I. I NTRODUCTION
The privacy and security of information transmission are
extremely important for wireless communications and networking [1]–[6]. Due to the broadcasting nature and lack of
physical boundaries of wireless transmission, the information
is readily intercepted by unauthorized users, and the growing
This work was supported in part by the National Natural Science Foundation
of China (Nos. 61472190, 6147210, 61501238, and 61271230), the Open
Research Fund of National Key Laboratory of Electromagnetic Environment,
China Research Institute of Radiowave Propagation (No. 201500013), and the
open research fund of National Mobile Communications Research Laboratory,
Southeast University, China (No. 2013D02).
Feng Shu, Wei Zhu, Jun Li, and Jinhui Lu are with the School of Electronic
and Optical Engineering, Nanjing University of Science and Technology,
210094, CHINA. Email: {shufeng, wei.zhu, jinhui.lu, Jun.li}@njust.edu.cn
Xiangwei Zhou is with the Division of ECE, Louisiana State University,
102 South Campus Drive, Baton Rouge, LA 70803, USA.
Feng Shu is also with the National Key Laboratory of Electromagnetic
Environment, China Research Institute of Radiowave Propagation, Qingdao
266107, China, and National Mobile Communications Research Laboratory,
210096, Southeast University, China.

cyber criminal events in mobile terminals have exposed the
enormous hidden risks in wireless communications and networking. Traditional encryption techniques can only provide
computational security and rely heavily on the complexity of
their keys, that are very easy to be cracked if an efficient
method to solve the corresponding mathematical problem is
found [1]–[4]. In the emerging social aware network, in order
to take full advantage of wireless network resources, several
incentive mechanisms based on credit or friendship had been
presented to stimulate selfish users to forward data for other
terminals, which will result in more serious security problems
if the useful messages are intercepted by malicious users [5],
[6]. In 1975, Wyner first proposes the remarkable wiretap
channel model and lays the foundations of information theory
of physical-layer security communications in [7]–[9]. In recent
years, the physical-layer security has ever been becoming
a promising research field in wireless communications and
networking by exploiting the physical-layer characteristics of
wireless channels [10]–[15].
As a novel physical-layer security technology in wireless
communications and networking, directional modulation (DM)
has been attracting widespread attention and research activities
from both academia and industry. In [16], the authors introduce
a technique of near-field direct antenna modulation (NFDAM),
changing the antenna boundary conditions at the symbol rate
thereby modulating the phase and amplitude of the antenna
pattern. A simplified structure of DM synthesis relying on
actively excited elements in phased arrays is described in [17],
[18]. Subsequently, the artificial-noise-aided security emerges
and is applied to confidential transmission [19]–[21], which
enables the transmitter to transmit artificial noise (AN), as an
interference signal, and confidential messages together, such
that the AN is only used to interfere with the eavesdroppers
without affecting the legitimate receivers. The algorithm in
[11] allocates that all the remaining available transmit power to
transmit AN in the case that the desired receivers are guaranteed to achieve a target signal-to-interference-plus-noise-ratio
(SINR). Based on the previous work, the authors in [22], [23]
initiate the DM research on the baseband signal processing
and propose to add the orthogonal vector, which can be chosen
and updated in the null space of channel vector at the desired
direction, to the transmitted baseband signal as AN [22], [23].
An orthogonal-vector-approach-based synthesis of multi-beam
DM follows immediately from [24]. The methods in [22] and
[23] are shown to perform very well for the perfect direction
angles but be quite sensitive to the measurement errors of

2

direction angles. To reduce the impact of the measurement
errors, the authors in [25] propose a robust DM synthesis
method that is capable of achieving an excellent performance
of bit error rate under imperfect direction angle in a singledesired-user scenario. In [25], a closed-form projection matrix
is derived to force the AN to the null space of steering vector
of the desired direction by utilizing the uniform distribution
of measurement angle errors. In [26], the authors propose a
robust beamforming scheme for multi-beam DM broadcasting
systems and derive the expressions of different beamformers
for various scenarios, in which case only one confidential
message stream is broadcasted to multiple desired users.
However, most existing work focuses on robust and nonrobust synthesis schemes of DM. Such schemes require perfect
knowledge of direction angles or imperfect one with distribution of angle measurement errors at the transmitters. In a
practical DM system, it is impossible to obtain perfect knowledge of direction angles or precisely model the distribution
of angle measurement errors. This motivates us to develop a
low-complexity and robust synthesis scheme needing only the
estimated direction angles, without requiring perfect direction
angles or the distribution of angle measurement errors under
the imperfect case.
In this paper, we consider multiple confidential message
stream transmission of using directional modulation in MUMIMO scenario. We propose a main-lobe-integration (MLI)
based leakage synthesis scheme. This method is based on
leakage idea in [27]–[30] and MLI, and shown robust to angle
measurement errors in our simulation. Our main contributions
are as follows:
• We propose a totally distinct robust beamfroming method
of using the concept of leakage in multi-user MIMO (MUMIMO) situation, which does not require the statistical knowledge of direction angle measurement error, i.e., to estimate
the variances of direction error. However, in [25], [26], the
proposed robust methods need to predict the variances of
direction estimation errors or even know the distributions of
measurement errors in advance.
• The robust method proposed by us simultaneously sends
multiple distinct independent parallel confidential message
streams to the associated multi-desired-receivers. However,
those methods in [25] and [26] transmit only single confidential message stream towards one or more desired users.
• The beamforming vectors for different desired users are
individually designed and distinct in this paper. This guarantees that any desired receiver can not intercept confidential
messages of others desired user, which further improves the
transmit security. However, in [26], all desired-user channels
are combined into one total large virtual channel and only
single identical desired beamforming vector is required and
devised for all desired users due to broadcasting scenario.
• In the following, AN is viewed as a virtual useful signal
for eavesdroppers and we construct the AN projecting matrix
by minimizing the leakage of AN power to the main-lobes
of the desired directions. This will reduce the effect of AN
on desired receivers and maximize the effect of AN on
eavesdroppers. In [26], the AN projection matrix is designed

by using the criterion of maximizing receiving-AN-to-signalnoise ratio at undesired receivers.
The rest of the paper is organized as follows. System model
is described in Section II. Section III presents the proposed
robust synthesis method with desired angle uncertainty under
two cases: imperfect and unknown eavesdropper directions.
The numerical results are shown and discussed in section IV,
and Section V concludes this paper.
Notations: throughout the paper, matrices, vectors, and
scalars are denoted by letters of bold upper case, bold lower
case, and lower case, respectively. Signs (·)T , (·)H , (·)−1
and tr(·) denote matrix transpose, conjugate transpose, MoorePenrose inverse and trace, respectively. Operation (x)+ returns
zero if x is negative, otherwise x is returned. The notation E{·}
refers to the expectation operation. The symbol IN denotes the
N × N identity matrix.
Since there are a lot of variables in this paper, we have
summarized main variables in Table I.
TABLE I
S UMMARY OF MAIN VARIABLES
N
K
M
θ dk
θem
s
dk
vk
TAN
z
Ps
α1
α2
β1
β2
h(θ)
ϕθ (n)
y(θ)
ω
Ck (θ)
Csec
θ̂dk
θ̂em
θBW
S
RS

Number of antennas at the base station(BS)
Number of desired users
Number of eavesdroppers
The kth desired direction angle
The mth eavesdropping direction angle
Transmit signal vector at the BS
The kth confidential message
Beamforming vector of the kth confidential message
Projection matrix of AN
Random AN vector
Total transmit power constraint at the BS
Normalized power factors for confidential messages
Normalized power factors for AN
Power allocation of confidential messages
Power allocation of AN
Normalized steering vector along the direction θ
Phase difference between the nth element in h(θ) and the
array phase center
Received signal along the direction θ
Additive white Gaussian noise(AWGN) with distribution
2)
CN (0, σω
Achievable rate of receiving the kth useful data along the
direction θ
Secrecy sum-rate
The kth estimated desired direction angle
The mth estimated eavesdropping direction angle
Beam width between first nulls (BWFN) for a long broadside
array
Integral interval with single continuous interval or the union
of several subintervals
H
Integral result
R of the matrix h(θ)h(θ) within the interval S,
i.e., RS = S h(θ)hH (θ) dθ

II. S YSTEM M ODEL
A schematic diagram of directional modulation MU-MIMO
system is shown in Fig. 1, where N is the number of elements
at the base station (BS). In such a system, we assume that
K independent confidential messages {dk }K
k=1 are transmitted towards K distinct desired users with direction angles
{θd1 , θd2 , . . . , θdK }. Additionally, there are M eavesdroppers
with direction angles {θe1 , θe2 , . . . , θeM }. Due to the irregular
reflection and refraction in the multipath fading channel, the

3

Desired user K
Eavesdropper 1

Desired user 1



If E dk dH
= 1, vkH vk = 1, and E zzH =
k
(3) and (4) can be further simplified as


α22 tr TAN TH
AN = N − K,
and

α21 =

Eavesdropper M

Ă

1
N −K IN −K ,

(5)

1
.
K

(6)

The N × 1 vector s passes through the LoS channel, the
received signal along direction θ is given by
y(θ) = hH (θ)s + ω,

Ă
Element 1

Element 2

Element N

Antenna array at base station

where ω is the additive white Gaussian noise (AWGN) with
distribution CN (0, σω2 ), and

T
1
j2πϕθ (n)
j2πϕθ (N ) 
θ (1)
h(θ) = √ e|j2πϕ
{z }, · · · , e| {z }, · · · , e| {z }
N
h1 (θ)

Fig. 1. Directional modulation MU-MIMO system, where different gray-scale
reds denote different confidential message streams for different desired users.

propagation direction of the signal is unpredictable, while
the directional modulation is very sensitive to the arrival
direction at desired receiver in the wireless channel. Until
now, the extension of directional modulation to multi-path
fading channel is a challenging problem. For example, if the
blocking object locates at the eavesdropper direction, then it
can reflect the AN to the desired receiver. We call the result
as the effect of gathering AN. This effect will dramatically
degrade the performance of the desired receiver by collecting
a large amount of AN from reflecting objects at the undesired
directions. Therefore, in this paper, we only consider the lineof-propagation channels.

= 1 is sent
The confidential message dk with E dH
k dk
to the desired receiver k, and multiplied by an N × 1
beam-forming vector vk before being transmitted through the
channel, where vk is called the confidential useful vector for
desired user k below. Thus, the transmit signal vector is written
as
K
p
p X
(1)
s = α1 β1 Ps
vk dk + α2 β2 Ps TAN z,
{z
}
|
k=1
AN
{z
}
|
Confidential messages

where Ps is the total transmit power constraint at BS, β1 and
β2 are the power allocation between confidential messages and
AN such that
β12 + β22 = 1,
(2)

α1 and α2 are the normalized power factors for confidential
messages and AN such that

 
= 1,
(3)
α22 E tr TAN zzH TH
AN
and TAN is the projection matrix of forcing AN to the
eavesdropping directions. The transmit signal vector in (1)
satisfies the following power constraint
(K K
)
XX
α21 E
= 1.
(4)
vkH′ vk dk dH
k′
k=1 k′ =1

(7)

hn (θ)

hN (θ)

(8)
is the normalized steering vector along the direction θ with
function ϕθ (n) defined by
ϕθ (n) ,

(n−(N +1)/2)d cos θ
, (n
λ

= 1, 2, . . . , N ),

(9)

where d denotes the spacing between two adjacent elements
of transmit antenna array and λ is the wavelength of transmit
carrier. According to (7), the received signal at the kth desired
user is given by
y(θdk ) = hH (θdk )s + ωdk
K
X
p
p
Ps hH (θdk )vk dk + α1 β1 Ps hH (θdk )
vi di
{z
}
i=1,i6=k
Useful data
{z
}
|
Interference from other users
p
+ α2 β2 Ps hH (θdk )TAN z + ωdk ,
(10)
|
{z
} |{z}

= α1 β1
|

AN

AWGN

where the first term of the above expression is the useful
received signal for user k, the second one is the multi-user
interference from other users, the third one is the AN, and the
last one is the AWGN with distribution CN (0, σd2k ) at receiver.
Similarly, the received signal at eavesdropper m is
y(θem ) = hH (θem )s + ωem
= α1 β1

K
X
p
vk dk
Ps hH (θem )
k=1

|

+ α2 β2
|

{z

Confidential messages

}

p
Ps hH (θem )TAN z + ωem ,
{z
} |{z}
AN

(11)

AWGN

where the first term of the above expression is composed of
the confidential messages intercepted by eavesdropper m, the
second one is the AN to disturb eavesdropper m, and the last
one is the AWGN with distribution CN (0, σe2m ).
To evaluate the security performance for multi-user scenario
in this paper, the secrecy sum-rate is adopted and can be defined as the sum of difference in available rate receiving useful

4

Ck (θdk ) , I(y(θdk ); [dk , θdk ])
α21 β12 Ps hH (θdk )vk vkH h(θdk )

= log2 (1 +
σd2k +

K
P

i=1,i6=k

).

(13)

α21 β12 Ps hH (θdk )vi viH h(θdk ) + α22 β22 Ps hH (θdk )TAN TH
AN h(θdk )

Ck (θem ) , I(y(θem ); [dk , θem ])

α21 β12 Ps hH (θem )vk vkH h(θem )

= log2 (1 +
σe2m +

K
P

i=1,i6=k

data between secure transmission channel and eavesdropper
channel, [31], [32],
Csec =

K h
i+
X
Ck (θdk ) − max Ck (θem ) ,

k=1

m

).

(14)

α21 β12 Ps hH (θem )vi viH h(θem ) + α22 β22 Ps hH (θem )TAN TH
AN h(θem )

(12)

where Ck (θdk ) and Ck (θem ) shown in (13) and (14) are
the achievable rates of receiving useful data dk along the
kth desired direction θdk and the mth eavesdropper direction
θem , respectively. Function I(y; [d, θ]) denotes the mutual
information along direction θ between the input d and the
output y.
The secrecy sum-rate in (12) is one of the most important
metrics for assessing the performance of the DM system.
Under the multi-user scenario, the expression of the secrecy
sum-rate in (12) contains multiple variables, i.e., Rk (θdk ) and
Rk (θem ). In addition, Rk (θdk ) and Rk (θem ) are both K + 1
coupled variables with {vdk }K
k=1 and TAN as shown in (13)
and (14), respectively. Apparently, it is challenging to solve
the K + 1 coupled optimization problems. Maximizing the
secrecy sum-rate in (12) directly is an NP-hard problem with
exponential complexity. Maybe, with a high-computational
amount, only suboptimal solution is achieved. On the other
hand, minimizing the bit error rate (BER), which is another
important metric for evaluating the DM system, of a selected
user will depend heavily on other users’ beamforming vectors and direction angles. The coupled property requires the
iteration operation. Additionally, it is not easy to guarantee
its convergence. To address the above, we propose a mainlobe-integration based leakage beamforming method, which
can provide an approximate closed-form solution with lowcomplexity and high performance. More importantly, due
to main-lobe-integration, it is also robust to direction angle
measurement errors.
III. P ROPOSED ROBUST S YNTHESIS M ETHOD WITH
D ESIRED A NGLE U NCERTAINTY
In Fig. 2, we sketch the basic idea of the proposed robust
precoding method based on MLI and leakage. As shown In
Fig. 2, if desired user k is chosen as the current desired user,
the corresponding confidential useful beamforming vector is
given by minimizing its useful signal power leakage to the
remaining desired users and all eavesdroppers. Confidential
beamforming vector corresponding to each desired user is
designed individually in order to safeguard each desired user

privacy. The AN projection matrix TAN is optimized by
maximizing the AN power leakage to all eavesdroppers. In
other words, the influence of AN on all desired users is
minimized. Here, TAN is constructed in the all-in-one way not
individually. In the following, we mainly consider the imperfect desired direction, where imperfect means that there exists
with measurement errors on the measured desired direction
angles. Eavesdropper direction angles fall into two categories:
imperfect (See Subsection III.A) and unknown (See Subsection III.B). In the first scenario, both desired and eavesdropper
directions are imperfect. In other words, there usually exists
errors in the measured desired and eavesdropper directions.
In the second scenario, eavesdropper directions are unknown
while the measured desired directions are imperfect. In the
two situations, we accordingly show how to design the AN
projection matrix TAN and confidential useful beamforming
vectors vdk .
BS

AN
ĂĂ

ĂĂ
dk

dk

d1 d k -1 d k +1 d K

dk

Unselected
desired users

Selected
desired users

Ă

Ă

e1

ĂĂ

eM

Eavesdroppers

Fig. 2. Schematic diagram of the proposed scheme.

A. Imperfect desired and eavesdropper directions
In practice, the desired and eavesdropper directions are
unknown. In this case, the BS can estimate both the values of
desired and eavesdropper direction angles with traditional spatial spectrum estimation such as MUSIC, Capon and ESPRIT
[33]. Due to the effect of channel noise, there always exist
errors in the estimated direction angles. Given the estimated
direction angles θ̂dk and θ̂em associated with desired user k
and eavesdropper m, we define their main-lobe intervals as


θBW
θBW
,
(15)
, θ̂dk +
Sdk = θ̂dk −
2
2

5

and



θBW
θBW
Sem = θ̂em −
,
, θ̂em +
2
2

(16)

2λ
(17)
Nd
is the beam width between first nulls (BWFN) for a long
broadside array [34]. Therefore, the overall direction angle
intervals of all main-lobes of desired users and eavesdroppers
are

K 
[
θBW
θBW
, S̄d = [0, π] \Sd ,
Sd =
, θ̂dk +
θ̂dk −
2
2
k=1
(18)
and

M 
[
θBW
θBW
,
(19)
θ̂em −
, θ̂em +
Se =
2
2
m=1
θBW =

respectively. Based on the above discussion, the average power
of confidential message stream dk sent to desired user k is
(Z
)
2 2
H H
H
Pdk ,U = E
α1 β1 Ps dk vdk h(θ)h (θ)vdk dk dθ
=

Sdk

Sdk

α21 β12 Ps vdHk h(θ)hH (θ)vdk dθ.

(20)

The remaining power leakage of confidential message stream
dk to other desired users and all eavesdroppers is represented
as
K
X

E

i=1,i6=k

+

M
X

m=1

where S is the integral interval with single continuous interval
or the union of several subintervals and each element of RS
is the definite integration of the corresponding element in the
N × N matrix h(θ)hH (θ) , which is a function of θ, over
the interval S. The detailed derivation of matrix RS is given
in Appendix B. Maximizing the MLI-SLNR in (24) by using
the generalized Rayleigh-Ritz theorem in [35] yields that the
optimal vdk is the normalized eigenvector corresponding to
the largest eigenvalue of
"
#−1
σd2k θBW
IN + RSd \Sdk + RSe
R Sdk .
(26)
α21 β12 Ps
Until now we complete the design of vdk . Below, we
similarly construct the projection matrix TAN of AN. The
basic idea is to project less power of AN onto the subspace
spanned by all desired steering vectors and more onto its null
space. Here, we view the AN as a useful signal. The average
AN power sent to all main-lobes of all eavesdropper directions
is as follows
PAN,U
Z

 H

2 2
H H
=E
α2 β2 Ps tr h (θ)TAN zz TAN h(θ) dθ
Se
Z



α22 β22 Ps tr hH (θ)TAN E zzH TH
=
AN h(θ) dθ, (27)
Se

Pdk ,L
=

(25)

S

where

Z

where matrix RS is defined as
Z
RS =
h(θ)hH (θ) dθ,

E

(Z

Sdi

(Z

Sem

H
H
α21 β12 Ps dH
k vdk h(θ)h (θ)vdk dk dθ

H
H
α21 β12 Ps dH
k vdk h(θ)h (θ)vdk dk dθ

)

)
,
(21)

+

M
X

m=1

Z

Sd

Sdi

Sem

α21 β12 Ps vdHk h(θ)hH (θ)vdk dθ,


α22 β22 Ps  H
(28)
tr TAN RSe TAN .
N −K
The average leakage power of AN to main-lobes of all desired
directions is given by
PAN,U =

PAN,L
Z
=E

which can be simplified as
Z
K
X
Pdk ,L =
α21 β12 Ps vdHk h(θ)hH (θ)vdk dθ
i=1,i6=k


1
IN −K and tr (AB) = tr (BA), the
Given E zzH = N −K
above equation is reduced to

(22)

MLI-SLNR(vdk ) = R

Sdk

Pdk ,U
.
σd2k dθ + Pdk ,L

(23)

which, according to Appendix A, can be simplified as
vdHk RSdk vdk
σd2 θBW
k
α21 β12 Ps

z



H
TH
AN h(θ)h (θ)TAN z



dθ .
(29)

Similar to (28), we have

α22 β22 Ps  H
tr TAN RSd TAN .
(30)
N −K
Combining the above two expressions and using the definition
of SLNR, we have the MLI-SLNR of AN as
PAN,U
dθ + PAN,L


tr TH RSe TAN
= (N −K)Mθ σAN
,
 H
2
BW e
+
tr
T
R
T
2
2
S
AN
d
AN
α2 β2 Ps
(31)

MLI-SLNR(TAN ) = R

σ2
Se e

MLI-SLNR(vdk ) =
vdHk

H

PAN,L =

where the first and second terms are the leakage powers to
other desired users and all eavesdroppers, respectively. We
define the corresponding MLI-SLNR as follows





α22 β22 Ps tr



IN + RSd \Sdk + RSe vdk

, (24)

where
σe2 =

M
1 X 2
σ .
M m=1 em

(32)

6

As a result, the optimal TAN when maximizing the MLISLNR in (31) is composed of the N − K normalized eigenvectors corresponding to the N − K largest eigenvalues of
matrix

−1
(N − K)M θBW σe2
(33)
I
+
R
R Se .
N
Sd
α22 β22 Ps
B. Unknown eavesdropper directions
In the following, we consider a more practical scenario,
where the BS doesn’t know the exact or estimated direction
values of eavesdroppers. In such a situation, all the remaining
angle region excluding the union of main-lobes of all desired
directions, which is actually the complementary set S̄d of
Sd , will be viewed as the potential eavesdropper directions.
For desired user k, the main-lobe profile of all eavesdropper
and the remaining desired directions are virtually viewed
as a potential set of intercepting directions, defined as the
complementary set S̄dk of its main-lobe region Sdk , i.e.,
S̄dk = [0, π] \Sdk

 

θBW [
θBW
= 0, θ̂dk −
θ̂dk +
, π .
2
2

(34)

Note that the useful part Pd′ k ,U of transmit signal power of
desired user k has the same expression as Pdk ,U in (20). The
power leakage of confidential message stream dk transmitted
by BS to the set S̄dk of all potential eavesdropper directions
is represented as
(Z
)
′
2 2
H H
H
Pdk ,L = E
α1 β1 Ps dk vdk h(θ)h (θ)vdk dk dθ
Z

=

S̄dk

S̄dk

α21 β12 Ps vdHk h(θ)hH (θ)vdk dθ.

(35)

Similar to (23), we can readily obtain the MLI-SLNR expression corresponding to the kth desired user as follows
MLI-SLNR′ (vdk ) = R

Pd′ k ,U

Sdk

=

vdHk

σd2k



dθ +

Pd′ k ,L

vdHk RSdk vdk
σd2 θBW
k
α21 β12 Ps

IN + RS̄d

k



. (36)
vdk

The optimal vdk to maximize the MLI-SLNR in (36) is the
generalized eigenvector corresponding to the largest normalized eigenvalue of
"
#−1
σd2k θBW
IN + RS̄d
R Sdk .
(37)
k
α21 β12 Ps
For the design of TAN with unknown directions of eavesdroppers, the angle range of all potential eavesdropping directions is S̄d , the average transmit AN power in S̄d is
′
PAN,U
Z
=E

S̄d



α22 β22 Ps tr hH (θ)TAN zzH TH
AN h(θ) dθ


α2 β 2 Ps 
= 2 2 tr TH
AN RS̄d TAN ,
N −K


(38)

′
to main-lobes of all
and the average leakage AN power PAN,L
desired directions is the same as (30). Based on the definition
of SLNR, we have the MLI-SLNR expression of AN as

MLI-SLNR′ (TAN )
′
PAN,U
= R
′
σ 2 dθ + PAN,L
S̄d e
 H

tr TAN RS̄d TAN
= (N −K)( π−Kθ )σ2
.

BW
e
+ tr TH
AN RSd TAN
α2 β 2 Ps
2

(39)

2

Similar to (33), via maximizing the MLI-SLNR in (39),
the optimal TAN is composed of the N − K normalized
eigenvectors corresponding to the N − K largest eigenvalues
of matrix

−1
(N − K)( π − KθBW )σe2
RS̄d .
(40)
I
+
R
N
S
d
α22 β22 Ps
IV. S IMULATION

AND

D ISCUSSION

To evaluate the bit error rate (BER) and secrecy-sum-rate
(SSR) performance of the proposed robust method, quadrature
phase shift keying (QPSK) is chosen, and main simulation
parameters are listed in Table II.
TABLE II
M AIN SIMULATION PARAMETERS AND THEIR VALUES
Parameter
d
N
K
{θdk }K
k=1
M
M
{θem }m=1
Ps
β1
β2
∆θmax

Value
λ/2
16
2
{60◦ , 120◦ }
3
{30◦ , 90◦ , 150◦ }
√1
√0.9
0.1
5◦

Here, signal-to-noise ratio is defined as SNR =
10 log10 (α21 β12 Ps /σω2 ), where σd2k = σe2m = σω2 , ∀k ∈
{1, 2}, ∀m ∈ {1, 2, 3}. All measurement errors of desired and
eavesdropper direction angles are assumed to be independently
uniform distributed over the interval [−∆θmax , ∆θmax ]. The orthogonal projection (OP) method in [24] and the conventional
leakage-based method in [28] are adopted as performance references, respectively. All the remaining angle region excluding
the union of main-lobes of all desired directions will be viewed
as the potential eavesdropper directions when eavesdroppers’
information are unable to be available.
Given the estimated desired and eavesdropper directions,
Fig. 3 (a) and (b) illustrate the curves of BER versus direction
angle of the proposed robust method in Section III-A, the OP
method in [24] and the conventional leakage-based method in
[28] for the first and second desired receivers, respectively. In
the Fig. 3 (a), it can be observed that our proposed method at
the desired direction 60◦ shows about an order-of-magnitude
improvement over that of the OP method and the conventional
leakage-based method. The lowest BER values of the OP
method and the conventional leakage-based method are approximately at the same level. As the receiver direction angle

7

100

100

10-1

10-1
10-3

10-3

BER

10-2

BER

10-2

10-3

10-3

10-4

10-4
60

60
10-4

10-4

Proposed method
Orthogonal method in [24]
Leakage-based method in [28]

0

30

60

90

120

150

180

Proposed method
Orthogonal method in [24]
Leakage-based method in [28]

0

30

Direction Angle (degree)

60

90

120

150

180

150

180

Direction Angle (degree)

(a) θd1 = 60◦

(a) θd1 = 60◦

100

100

10-1

10-1
-3

10

10-3

BER

10-2

BER

10-2

10-3 10-4

10-3 10-4
120

120
10-4

10-4

Proposed method
Orthogonal method in [24]
Leakage-based method in [28]

0

30

60

90

120

150

180

Direction Angle (degree)
(b) θd2 = 120◦

Proposed method
Orthogonal method in [24]
Leakage-based method in [28]

0

30

60

90

120

Direction Angle (degree)
(b) θd2 = 120◦

Fig. 3. BER versus direction angle with imperfect desired and eavesdropper
directions (SNR=14dB).

Fig. 4. BER versus direction angle with unknown eavesdropper directions
(SNR=14dB).

steers away from the desired direction, the BER performance
of both methods degrade rapidly. Compared to the remaining
two methods, the BER curve of our proposed method fluctuates less outside the main-lobe of the desired direction. For the
second desired user with direction θd2 = 120◦ in Fig. 3 (b),
the performance trend is similar to the first desired user.
Given the estimated desired directions and under unknown
eavesdropper directions, Fig. 4 (a) and (b) plot the curves of
BER versus direction angle of the proposed robust method in
Section III-B, the OP method in [24] and the conventional
leakage-based method in [28], respectively, for the first and
second desired receivers. From Fig. 4 (a), similar to Fig. 3 (a),
we can see that the proposed method still outperforms the
OP by approximately an order-of-magnitude in the desired
direction 60◦ , where the BER value of the proposed method
is about one-third of the conventional leakage-based method.

More importantly, the BER of the proposed method is about
40 percent in the remaining region outside the main-lobe of
the desired direction, and the curve fluctuates less than that of
the other two methods. This provides an effective barrier to
prevent potential eavesdroppers recovering confidential information. Fig. 4 (b) presents a similar BER performance trend
for the second desired direction 120◦ as Fig. 4 (a).
In the following, we will evaluate the performance of the
proposed method from the SSR aspect. Fig. 5 demonstrates
the curves of SSR versus SNR for the proposed method with
imperfect desired and eavesdropper directions in Section IIIA, the OP method in [24] and the conventional leakage-based
method in [28]. As can be seen from Fig. 5, with the increase
in SNR, the SSRs of all three methods increase continuously
and monotonously. The SSR of the proposed method is always
larger than that of the remaining two methods. Their SSR

18

16

16

14

Secrecy Sum Rate (bits/s/Hz)

Secrecy Sum Rate (bits/s/Hz)

8

14
12
10
8
6
4
Proposed method
Orthogonal method in [24]
Leakage-based method in [28]

2
0

0

5

10

15

20

25

30

12
10
8
6
4
Proposed method
Orthogonal method in [24]
Leakage-based method in [28]

2
35

0

SNR (dB)
Fig. 5. Secrecy sum rate versus SNR with imperfect desired and eavesdropper
directions.

values are up to 16.4 bit/s/Hz, 12.2 bit/s/Hz and 11.1 bit/s/Hz
at SNR = 35dB, respectively. Hence, the SSR improvement of
the proposed method over the OP method and the conventional
leakage-based method is significant. Additionally, the SSRs of
the three methods approach the same value as SNR goes to 0
dB. The rate gap between the proposed method and the other
two methods grows gradually with the increase in SNR.
In the scenario of unknown eavesdropper directions, Fig. 6
demonstrates the curves of SSR versus SNR for the proposed
method in Section III-B, the OP method in [24] and the
conventional leakage-based method in [28]. As indicated in
Fig. 6, the proposed method performs better than the OP
method and the conventional leakage-based method in terms
of SSR, in particular, in the medium and high SNR regions.
Their SSR values are up to 14.0 bit/s/Hz, 12.2 bit/s/Hz and
12.8 bit/s/Hz at SNR = 35 dB, respectively, in which case that
the proposed method shows an approximate 15 and 9 percent
SSR improvement over the OP and the conventional leakagebased method, respectively. Additionally, the SSRs of all three
methods tend to be the same in the low SNR region. The rate
gap between the proposed method and the other two methods,
similar to Fig. 5, grows gradually with the increase in SNR.
Note that the SSRs of the proposed method in Fig. 6 are less
than the corresponding SSRs of the proposed method in Fig. 5
in the middle and high regions. This implies that the measured
values of eavesdropper direction angles play an important role
in improving the SSR performance.
V. C ONCLUSION
In this paper, we study robust precoding in MU-MIMO
directional modulation systems under imperfect desired directions and eavesdropper directions or unknown eavesdropper directions. To realize the high-performance robust secure multimessage-stream simultaneous transmission in MU-MIMO systems, we propose an MLI-based precoding method. In the first
stage, the precoding vector per desired user is optimized to

0

5

10

15

20

25

30

35

SNR (dB)
Fig. 6. Secrecy sum rate versus SNR with unknown eavesdropper directions.

maximize the transmit power of useful confidential messages
along the main-lobe of the corresponding desired direction and
correspondingly minimize the transmit power of AN along all
eavesdropper directions. In the second stage, the AN projection
matrix is designed to force AN to all eavesdropper directions
with only a small amount of residual AN along the mainlobes of the desired directions. Due to the use of MLI, the
proposed method requires no perfect direction knowledge or
the distribution of measurement errors of direction angle, and
at the same time provides a robust performance. Simulation
results verify the performance benefits of our method. Compared with OP, the proposed method can achieve an oneorder-magnitude improvement on BER and at the same time
its secrecy-sum-rate along the desired directions is shown to
have a substantial enhancement over that of OP in the high
SNR region. The proposed scheme can be applied to the future
satellite communications, mobile communications, D2D, V2V,
unmanned-aerial-vehicles networks, and internet of things.
A PPENDIX A
S IMPLIFICATION OF MLI-SLNR(vdk )
Proof: The simplification of MLI-SLNR(vdk ) is shown
a
in (41), where = is achieved by extracting integral term
R
b
h(θ)hH (θ), = is achieved by replacing Sd h(θ)hH (θ) dθ,
k
R
R
H
H
Sdi h(θ)h (θ) dθ and Sem h(θ)h (θ) dθ by RSdk , RSdi

and RSem , respectively.
A PPENDIX B
D ERIVATION OF RS
Proof: The integral matrix RS in (25) is an N × N matrix.
Here we assume that set S is a union of I separate subintervals
as follows
I
[
Si ,
(42)
S=
i=1

where

i
i
Si = [θmin
, θmax
].

(43)

9

R

Sdk

MLI-SLNR(vdk ) =

R

Sdk

σd2k dθ +

Sdk

σd2k

K
P

i=1,i6=k

R

Sdi

a

α21 β12 Ps vdHk h(θ)hH (θ)vdk dθ

M R
P
α2 β 2 P vH h(θ)hH (θ)vdk dθ
α21 β12 Ps vdHk h(θ)hH (θ)vdk dθ +
Sem 1 1 s dk
m=1
R
α21 β12 Ps vdHk Sd h(θ)hH (θ) dθvdk
k

=
R

b

dθ +

=
σd2k θBW +
=

σd2k θBW +

=

σd2 θBW
k

α21 β12 Ps

=



vdHk

α21 β12 Ps vdHk

K
P

R

Sdi

RSi (p, q) =
=

i
θmax

i
θmin

Z

i
θmax

i
θmin

vdHk RSdk vdk
σd2 θBW
k
α21 β12 Ps



Sem

h(θ)hH (θ) dθvdk

.

(41)

Using the above definition or transformation, Eq. (44) is
rewritten as

hp (θ)hH
q (θ) dθ
1 j2π(p−(N +1)/2)d cos θ
λ
√ e
N

i
θi + θmax
= min
,
2

(44)

(45)

and
i
i
∆θi =θmax
− θmin
,

RSi (p, q)
Z
∆θ
1 π j2π(p−q)d cos( 2πi x+θ0i ) ∆θi
λ
·
dx
e
=
N −π
2π
Z π j2π(p−q)d cos( ∆θi x+θi )
0
∆θi
2π
λ
e
=
dx
2πN −π
Z π j2π(p−q)d[cos(θi ) cos( ∆θi x)−sin(θi ) sin( ∆θi x)]
0
0
2π
2π
∆θi
λ
dx
e
=
2πN −π
Z π
i
i
∆θi
e[apq cos(ci x)+bpq sin(ci x)] dx
=
2πN −π
∆θi
=
g (ai , bi , ci ),
(49)
N B pq pq
where gB (·) is the extension of the modified Bessel function
of the first kind with the integer order 0 [36], i.e.,
Z π
i
i
1
gB (aipq , bipq , ci ) =
e[apq cos(ci x)+bpq sin(ci x)] dx, (50)
2π −π
aipq ,

(46)

respectively, and the new integral variable is
2π
(θ − θ0i ),
∆θi

∆θi
x + θ0i .
2π

j2π(p − q)d cos(θ0i )
,
λ

(51)

j2π(p − q)d sin(θ0i )
,
λ

(52)

bipq , −
(47)

and
ci ,

i.e.,
θ=

m=1

IN + RSd \Sdk + RSe vdk

Let us define the center point and length of integral interval
as

x=

M R
P

+ vdHk RSd \Sdk vdk + vdHk RSe vdk

j2π(q−(N +1)/2)d cos θ
1
λ
dθ
· √ e−
N
Z i
1 θmax j2π(p−q)d cos θ
λ
dθ.
e
=
i
N θmin

θ0i

+ α21 β12 Ps vdHk

i=1,i6=k
2 2
α1 β1 Ps vdHk RSdk vdk
M
K
P
P
RSem vdk
α21 β12 Ps vdHk
RSdi vdk + α21 β12 Ps vdHk
m=1
i=1,i6=k
α21 β12 Ps vdHk RSdk vdk
α21 β12 Ps vdHk RSd \Sdk vdk + α21 β12 Ps vdHk RSe vdk
vdHk RSdk vdk

Then the (p, q) entry of RSi has the following form
Z

h(θ)hH (θ) dθvdk

(48)

∆θi
.
2π

(53)

This completes the derivation of the (p, q) entry RSi (p, q)
of matrix RSi . Note that set S is a union of I separated
subintervals as shown in (41). Then the (p, q) element of

10

matrix RS is expanded as the summation of integrations over
all subintervals
RS (p, q) =

I
X
i=1

RSi (p, q) =

I
X
∆θi
i=1

N

gB (aipq , bipq , ci ). (54)

This completes the derivation of matrix RS .



R EFERENCES
[1] Y.-W. P. Hong, P.-C. Lan, and C.-C. J. Kuo, Signal Processing Approaches to Secure Physical Layer Communications in multi-antenna
wireless systems. New York, NY: Springer, SpringerBriefs in Electrical
and Computer Engineering, 2013.
[2] Y. Zou and J. Zhu, Physical-Layer Security for Cooperative Relay
Networks. Springer International Publishing, 2016.
[3] C. Huang, M. Ma, X. Liu, A. Liu, and Z. Zuo, “Unequal probability
marking approach to enhance security of traceback scheme in tree-based
wsns,” Sensors, vol. 17, no. 6, Jun. 2017.
[4] Y. Liu, M. Dong, K. Ota, and A. Liu, “Activetrust: Secure and trustable
routing in wireless sensor networks,” IEEE Trans. Inf. Forensics Security,
vol. 11, no. 9, pp. 2013–2027, Sept. 2016.
[5] Z. Ning, F. Xia, X. Hu, Z. Chen, and M. S. Obaidat, “Social-oriented
adaptive transmission in opportunistic internet of smartphones,” IEEE
Trans. Ind. Informat., vol. 13, no. 2, pp. 810–820, Apr. 2017.
[6] Z. Ning, L. Liu, F. Xia, B. Jedari, I. Lee, and W. Zhang, “CAIS: A
copy adjustable incentive scheme in community-based socially aware
networking,” IEEE Trans. Veh. Technol., vol. 66, no. 4, pp. 3406–3419,
April 2017.
[7] A. D. Wyner, “The wire-tap channel,” Bell. Syst. Tech. J., vol. 54, no. 8,
pp. 1355–1387, Oct. 1975.
[8] S. Leung-Yan-Cheong and M. Hellman, “The gaussian wire-tap channel,” IEEE Trans. Inf. Theory, vol. 24, no. 4, pp. 451–456, Jul. 1978.
[9] I. Csiszar and J. Korner, “Broadcast channels with confidential messages,” IEEE Trans. Inf. Theory, vol. 24, no. 3, pp. 339–348, May. 1978.
[10] M. Bloch, J. Barros, M. R. D. Rodrigues, and S. W. McLaughlin, “Wireless information-theoretic security,” IEEE Trans. Inf. Theory, vol. 54,
no. 6, pp. 2515–2534, Jun. 2008.
[11] A. Mukherjee and A. L. Swindlehurst, “Robust beamforming for security
in MIMO wiretap channels with imperfect CSI,” IEEE Trans. Signal
Process., vol. 59, no. 1, pp. 351–361, Jan. 2011.
[12] Y. L. Zou, J. Zhu, X. Wang, and V. Leung, “Improving physical-layer
security in wireless communications through diversity techniques,” IEEE
Netw., vol. 29, no. 1, pp. 42–48, Jan. 2015.
[13] Y. L. Zou, B. Champagne, W.-P. Zhu, and L. Hanzo, “Relay-selection
improves the security-reliability trade-off in cognitive radio systems,”
IEEE Trans. Commun., vol. 63, no. 1, pp. 215–228, Jan. 2015.
[14] X. Chen, D. W. K. Ng, and H. H. Chen, “Secrecy wireless information and power transfer: challenges and opportunities,” IEEE Wireless
Communications, vol. 23, no. 2, pp. 54–61, Apr. 2016.
[15] N. Zhao, F. R. Yu, M. Li, Q. Yan, and V. C. Leung, “Physical layer
security issues in interference-alignment-based wireless networks,” IEEE
Commun. Mag., vol. 54, no. 8, pp. 162–168, Aug. 2016.
[16] A. Babakhani, D. B. Rutledge, and A. Hajimiri, “Transmitter architectures based on near-field direct antenna modulation,” IEEE J. Solid-State
Circuits, vol. 43, no. 12, pp. 2674–2692, Dec. 2008.
[17] M. P. Daly, E. L. Daly, and J. T. Bernhard, “Directional modulation
technique for phased arrays,” IEEE Trans. Antennas Propag., vol. 57,
no. 9, pp. 2633–2640, Sept. 2009.
[18] ——, “Demonstration of directional modulation using a phased array,”
IEEE Trans. Antennas Propag., vol. 58, no. 5, pp. 1545–1550, May.
2010.
[19] S. Goel and R. Negi, “Guaranteeing secrecy using artificial noise,” IEEE
Trans. Wireless Commun., vol. 7, no. 6, pp. 2180–2189, Jun. 2008.
[20] X. Zhou and M. R. McKay, “Physical layer security with artificial
noise: Secrecy capacity and optimal power allocation,” in International
Conference on Signal Processing and Communication Systems, Sept.
2009, pp. 1–5.
[21] N. Zhao, F. R. Yu, M. Li, and V. C. Leung, “Anti-eavesdropping schemes
for interference alignment (IA)-based wireless networks,” IEEE Trans.
Wireless Commun., vol. 15, no. 8, pp. 5719–5732, Aug. 2016.
[22] Y. Ding and V. Fusco, “A vector approach for the analysis and synthesis
of directional modulation transmitters,” IEEE Trans. Antennas Propag.,
vol. 62, no. 1, pp. 361–370, Jan. 2014.

[23] ——, “Directional modulation far-field pattern separation synthesis
approach,” IET Microwaves Antennas and Propagation, vol. 9, no. 1,
pp. 41–48, Jan. 2015.
[24] ——, “Orthogonal vector approach for synthesis of multi-beam directional modulation transmitters,” IEEE Antennas Wireless Propag. Lett.,
vol. 14, pp. 1330–1333, Feb. 2015.
[25] J. Hu, F. Shu, and J. Li, “Robust synthesis method for secure directional modulation with imperfect direction angle,” IEEE Commun. Lett.,
vol. 20, no. 6, pp. 1084–1087, Jun. 2016.
[26] F. Shu, X. Wu, J. Li, R. Chen, and B. Vucetic, “Robust synthesis scheme
for secure multi-beam directional modulation in broadcasting systems,”
IEEE Access, vol. 4, pp. 6614–6623, 2016.
[27] A. Tarighat, M. Sadek, and A. H. Sayed, “A multi-user beamforming
scheme for downlink MIMO channels based on maximizing signalto-leakage ratios,” in Proceedings. (ICASSP ’05). IEEE International
Conference on Acoustics, Speech, and Signal Processing, 2005., vol. 3,
Philadelphia, PA, Mar. 2005, pp. 1129–1132.
[28] M. Sadek, A. Tarighat, and A. H. Sayed, “A leakage-based procoding
scheme for downlink multi-user MIMO channels,” IEEE Trans. Wireless
Commun., vol. 6, no. 5, pp. 1711–1721, Jun. 2007.
[29] F. Shu, M. M. Wang, Y. X. Wang, H. Q. Fan, and J. H. Lu, “An
efficient power allocation scheme for leakage-based precoding in multicell multiuser MIMO downlink,” IEEE Commun. Lett., vol. 6, no. 5, pp.
1053–1055, Oct. 2011.
[30] F. Shu, J. J. Tong, X. H. You, C. Gu, and J. J. Wu, “Adaptive robust MaxSLNR precoder for MU-MIMO-OFDM systems with imperfect CSI,”
Science China Information Sciences, vol. 59, no. 6, pp. 1–14, Jun. 2016.
[31] F. Oggier and B. Hassibi, “The secrecy capacity of the MIMO wiretap
channel,” IEEE Trans. Inf. Theory, vol. 57, no. 8, pp. 4961–4972, Aug.
2011.
[32] N. Li, X. Tao, and J. Xu, “Ergodic secrecy sum-rate for downlink
multiuser MIMO systems with limited csi feedback,” IEEE Commun.
Lett., vol. 18, no. 6, pp. 969–972, Jun. 2014.
[33] F. Gross, Smart Antennas for Wireless Communications: With MATLAB.
McGraw-Hill Education, 2005.
[34] J. D. Kraus and R. J. Marhefka, Antennas for All Applications, Third
Edition. McGraw-Hill, 2006.
[35] R. A. Horn and C. R. Johnson, Matrix Analysis. Cambridge, U.K.:
Cambridge Univ. Press, 1987.
[36] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and
Products, Seventh Edition. San Diego, CA, USA: Academic, 2007.

