IEEE COMMUNICATIONS LETTERS

1

Joint Fractional Time Allocation and Beamforming
for Downlink Multiuser MISO Systems

arXiv:1706.01771v2 [] 28 Aug 2017

Van-Dinh Nguyen, Hoang Duong Tuan, Trung Q. Duong, Oh-Soon Shin, and H. Vincent Poor

Abstract—It is well-known that the traditional transmit beamforming at a base station (BS) to manage interference in serving
multiple users is effective only when the number of users is less
than the number of transmit antennas at the BS. Non-orthogonal
multiple access (NOMA) can improve the throughput of users with
poorer channel conditions by compromising their own privacy
because other users with better channel conditions can decode the
information of users in poorer channel state. NOMA still prefers
that the number of users is less than the number of antennas at
the BS transmitter. This paper resolves such issues by allocating
separate fractional time slots for serving the users with similar
channel conditions. This enables the BS to serve more users within
the time unit while the privacy of each user is preserved. The
fractional times and beamforming vectors are jointly optimized
to maximize the system’s throughput. An efficient path-following
algorithm, which invokes a simple convex quadratic program at
each iteration, is proposed for the solution of this challenging
optimization problem. Numerical results confirm its versatility.
Index Terms—Beamforming, fractional time allocation, nonconvex optimization, path-following method.

I. I NTRODUCTION
In multi-user communication, the signal received by any user
(UE) is a superposition of the desired information and the
information intended for other UEs. Transmit beamforming at a
base station (BS) plays a pivotal role in focusing the energy of
the desired signal and suppressing the multi-user interference
to achieve better throughput at UEs [1], [2]. For effective
beamforming, the number of UEs usually does not exceed the
number of transmit antennas. Using more transmit antennas
and thus increasing the dimensionality of beamforming vectors
can improve the UEs’ throughput with reduced transmit power.
However, allocating more transmit power under fixed number
of transmit antennas still does not necessarily improve the UEs’
throughput.
This work was supported in part by the National Research Foundation of
Korea (NRF) grant funded by the Korean government (MSIT) (No. NRF2017R1A5A1015596), in part by Basic Science Research Program through
the NRF funded by the Ministry of Education (No. 2017R1D1A1B03030436),
in part by the U.K. Royal Academy of Engineering Research Fellowship under
Grant RF1415/14/22 and by the U.K. Engineering and Physical Sciences
Research Council (EPSRC) under Grant EP/P019374/1, and in part by the
KAUST Grant No. OSR-2016-CRG5-2958-0 and U.S. National Science Foundation Grant ECCS-1647198. (Corresponding author: Oh-Soon Shin.)
V.-D. Nguyen and O.-S. Shin are with the School of Electronic Engineering
and the Department of ICMC Convergence Technology, Soongsil University,
Seoul 06978, Korea (e-mail: {nguyenvandinh, osshin}@ssu.ac.kr).
H. D. Tuan is with the Faculty of Engineering and Information Technology,
University of Technology Sydney, Broadway, NSW 2007, Australia (email:
tuan.hoang@uts.edu.au).
T. Q. Duong is with the School of Electronics, Electrical Engineering and
Computer Science, Queen’s University Belfast, Belfast BT7 1NN, United
Kingdom (e-mail: trung.q.duong@qub.ac.uk).
H. V. Poor is with the Department of Electrical Engineering, Princeton
University, Princeton, NJ 08544 USA (e-mail: poor@princeton.edu).

Non-orthogonal multiple access (NOMA) [3], [4] has been
introduced to improve the UEs’ throughput by allowing UEs
with better channel conditions to access and decode the signals,
which are intended for the UEs with poorer channel conditions.
In other words, the UEs with poorer channel conditions can
achieve higher throughput by compromising their information
privacy in NOMA [5]. This privacy compromise for the UEs
with poor channel conditions is unavoidable in NOMA. Moreover, by restricting all beamforming vectors in the same space,
beamforming in NOMA still needs that the dimension of this
space, which is equal to the number of transmit antennas,
should not be less than the number of UEs to allow the
suppression of multiuser interference [6].
It is noteworthy that the UEs with good channel conditions
may need only a fraction of the time unit to achieve their needed
throughput. Therefore, by serving them only per fractional
time, the BS still has the remaining time room to serve the
UEs with poor channel conditions. In this way, the information
privacy for each UE is preserved because all UEs are allowed
to decode their own information only. More importantly, the
number of UEs supported at the same fractional time is
effectively reduced. Thus, the BS will not need more transmit
antennas to suppress the multi-user interference. In this letter,
the fractional time allocation to UEs with similar channel
conditions and beamforming are enhanced for improving the
system’s throughput and meeting the UEs’ quality-of-service
(QoS) in terms of the throughput requirement. This problem is
mathematically modelled by a highly nonconvex optimization
problem, for which a path-following computational procedure
of low complexity is then developed for its computation.
Finally, the numerical examples are provided to demonstrate
the advantage of the proposed optimization scheme.
Notation. We use bold upper-case letters for matrices, bold
lower-case letters for column vectors, lower-case letters for
scalars. ℜ{x} denotes the real part of a complex number
x. The notation (·)H stands for the Hermitian transpose.
x ∼ CN (η, Z) means that x is a random vector following
a complex circular Gaussian distribution with mean η and
covariance matrix Z.
II. S YSTEM M ODEL AND P ROBLEM F ORMULATION
Consider a multiuser downlink system over a given frequency
band with a BS equipped with Nt > 1 antennas in serving
2K single-antenna UEs as illustrated by Fig. 1. There are K
UEs (1, k), k = 1, . . . K, which are located in a zone nearer
to the BS, called by zone-1, and K UEs (2, k), k = 1, . . . , K,
which are located in a zone farer from the BS, called by zone-2.
Denote by K , {1, 2, . . . , K} and M = {1, 2}×K. Within one
time unit, BS uses the fraction time (FT) τ1 := τ (0 < τ < 1)

2

IEEE COMMUNICATIONS LETTERS

500

to serve UEs (1, k) and uses the remaining FT τ2 := (1 − τ )
to serve UEs (2, k).
The BS deploys a transmit beamformer wi,k ∈ CNt ×1 to
deliver the information signal xi,k with E{|xi,k |2 } = 1 to UE
(i, k). Let hi,k ∈ CNt ×1 be the channel vector from the BS
to UE (i, k), which is assumed to follow frequency flat fading
with the effects of both large-scale pathloss and small-scale
fading counted. The complex baseband signal received by UE
(i, k) can be expressed as
X

yi,k = hH
i,k wi,k xi,k +

Ri,k wi , τi

2
|hH
i,k wi,k |

2
2
|hH
i,k wi,j | + σi,k

j∈K\{k}

2
2
|hH
i,k wi,j | + σi,k



(1)

under the additional condition
ℜ{hH
i,k wi,k }

≥ 0, (i, k) ∈ M.

2
(ℜ{hH
i,k wi,k })
(i,j)∈M\{(i,k)}

2
2
|hH
i,k wi,j | + σi,k

w,τ

X

(i,k)∈M

y-coordinate (m)

(2,4)

Ri,k (wi , τi ) s.t. (2),


Ri,k wi , τi ≥ R̄i,k , ∀(i, k) ∈ M,

max
τ1 kw1 k2 + τ2 kw2 k2 ≤ Pbs
,
τ1 ≥ 0, τ2 ≥ 0, τ1 + τ2 ≤ 1.

100

300

500

Scenario with K = 4.

III. C ONVEX Q UADRATIC - BASED PATH - FOLLOWING
METHOD

1/α1 + 1/α2 ≤ 1, α1 > 0, α2 > 0.

(i,k)∈M

Ri,k (wi , 1/αi ) ≥ R̄i,k , ∀(i, k) ∈ M,
2



(3)

(4a)
(4b)
(4c)
(4d)

Here R̄i,k sets a minimum throughput requirement for UE (i, k)
max
and Pbs
is a given power budget. Since Ri,k (wi , τi ) is a
nonconcave function, the optimization problem (4) is regarded
as a highly nonconvex optimization problem, for which finding
a feasible point is already computationally difficult. The next
section is devoted to a computational path-following procedure
for its solution.

(5)

The problem (4) can be equivalently expressed as
X
Ri,k (wi , 1/αi ) s.t. (2), (5), (6a)
max Φ(w, α ) ,

(2)

with the full inter-zone interference. Here and in the sequence
w , (w1 , w2 ) and τ , (τ1 , τ2 ).
We are interested in the following problem of jointly designing FT (τ1 , τ2 ) and the beamformers (w1 , w2 ) to maximize the
system sum throughput (ST):
max

(2,1)

(1,3)

−500
−500 −300 −100

α
w,α

The main advantage of this FT-based beamforming scheme is
that there is no inter-zone interference in (1) that is in contrast
with the conventional scheme to concurrently serve all UEs,
under which the throughput at UE (i, k) is

′
Ri,k (w) = ln 1 + P

(2,2)

The most important step is to provide an effective lower
bounding approximation for the function Ri,k (wi , τi ) defined
by (1) to facilitate a path-following computational procedure of
the problem (4). We use the variable changes α , (α1 , α2 ) =
(1/τ1 , 1/τ2 ), which satisfy the following convex constraint:



which by [7] can be equivalently reformulated by
2
(ℜ{hH
i,k wi,k })

(1,4)

(1,1)

−300

Fig. 1.

j∈K\{k}



Ri,k wi , τi = τi ln 1 + P

−100

(1,2)

x-coordinate (m)

where the first term is the desired signal, the second term is the
2
multi-user interference, and the third term ni,k ∼ CN (0, σi,k
)
is additive noise. For wi , (wi,k )k∈K , the throughput at UE
(i, k) is

= τi ln 1 + P

(2,3)

100

hH
i,k wi,j xi,j + ni,k

j∈K\{k}



300

BS
UEs in zone-1
UEs in zone-2

2

(1 − 1/α2 )kw1 k + kw2 k /α2 ≤

(6b)

max
Pbs
.

(κ)

(6c)
(κ)

At a feasible point (w(κ) , α (κ) ), define xi,k , ℜ{hH
i,k wi,k } >
P
(κ)
(κ)
(κ) 2
(κ)
(κ) 2
H
2
0, yi,k , j∈K\{k} |hi,k wi,j | + σi,k , di,k , (xi,k ) /yi,k >
(κ)

(κ)

(κ)

(κ)

(κ)

(κ)

0, ai,k , 2Ri,k (wi , 1/αi ) + di,k /αi (di,k + 1) >
(κ)
(κ)
(κ)
(κ) (κ)
0, bi,k , (di,k )2 /αi (di,k + 1) > 0, and ci,k ,
(κ)
(κ)
Ri,k (wi , 1/(αi )2 ) > 0. It follows from the inequality (13)
in the appendix that
(κ)

Ri,k (wi , 1/αi ) ≥ Ri,k (wi , αi )
over the trust region
(κ)

H
2ℜ{hH
i,k wi,k } − ℜ{hi,k wi,k } > 0, ∀(i, k) ∈ M,

(7)

for the concave function
(κ)
Ri,k (wi , αi )

,

(κ)
(κ)
ai,k −bi,k

2
2
H
j∈K\{k} |hi,k wi,j | + σi,k
(κ)
−ci,k αi .
(κ)
(κ)
H
xi,k (2ℜ{hi,k wi,k } − xi,k )

P

Next, due to the convexity of function kw1 k2 /α2 , it is true that
(κ)
(κ)
(κ)
(κ)
kw1 k2 /α2 ≥ 2ℜ{(w1 )H w1 }/α2 − (kw1 k2 /(α2 )2 )α2 .
An inner convex approximation of nonconvex constraint (6c)
is then given by
(κ)

(κ)

kw1 k2 + kw2 k2 /α2 − 2ℜ{(w1 )H w1 }/α2
(κ)

(κ)

+(kw1 k2 /(α2 )2 )α2

max
≤ Pbs
. (8)

Initialized by a feasible point (w(0) , α (0) ) for (6), the following
convex quadratic program (QP) is solved at the κ-th iteration

NGUYEN et al.: JOINT FRACTIONAL TIME ALLOCATION AND BEAMFORMING FOR DOWNLINK MULTIUSER MISO SYSTEMS

(κ+1)

TABLE I
S IMULATION PARAMETERS
Parameters
Noise power density
Path loss from the BS to UE (i, k), σPL
Radius of cell
Coverage of zone-1 UEs
Distance between the BS and nearest user

α(κ+1)

to generate the next feasible point (w
,α
):
X
(κ)
Ri,k (wi , αi ) s.t.
max Φ(κ) (w, α ) ,
α
w,α

(i,k)∈M

≥ R̄i,k , ∀(i, k) ∈ M, (2), (5), (7), (8).

(9)

As problem (9) involves m = 2(3K + 1) quadratic and linear
constraints, and n = 2(KNt + 1) real decision variables, its
computational complexity is O(n2 m2.5 + m3.5 ).
Note that Φ(w, α )
≥
Φ(κ) (w, α ) ∀(w, α ),
(κ)
(κ)
(κ)
and Φ(w , α )
=
Φ (w(κ) , α (κ) ). Moreover,
(κ)
(κ+1)
(κ+1)
Φ (w
,α
)
>
Φ(κ) (w(κ) , α (κ) ) whenever
(κ+1)
(κ+1)
(κ)
(w
,α
) 6= (w , α (κ) ) because the former
and the latter, respectively, are the optimal solution and
feasible point for (9). Therefore, Φ(w(κ+1) , α (κ+1) ) ≥
Φ(κ) (w(κ+1) , α (κ+1) ) > Φ(κ) (w(κ) , α (κ) ) = Φ(w(κ) , α (κ) ),
showing that (w(κ+1) , α (κ+1) ) is a better feasible point than
(w(κ) , α (κ) ) for (6). The sequence {(w(κ) , α (κ) )} of improved
feasible points for (6) thus converges at least to a locally
optimal solution satisfying the Karush-Kuh-Tucker conditions
[8]. We summarize the proposed QP-based path-following
procedure in Algorithm 1.
Generation of an initial point.: Initialized from a feasible
point (w(0) , α (0) ) for constraints (5) and (8), we iterate the
convex program
(κ)

max min Ri,k (wi , αi )/R̄i,k
α (i,k)∈M
w,α

(κ)

s.t. (2), (5), (7), (8) (10)
(κ+1)

(κ+1)

)/R̄i,k ≥ 1 to
, αi
till reaching min(i,k)∈M Ri,k (wi
make (w(κ+1) , α (κ+1) ) feasible for (6) and thus usable as an
initial feasible point for implementing Algorithm 1.
IV. N UMERICAL R ESULTS
Monte Carlo simulations have been implemented to evaluate
the performance of the proposed algorithm for K = 4 (8 UEs)
and Nt = 5 per the scenario in Fig. 1. The channel vector hi,k
between the BS and UE√(i, k) at a distance di,k (in kilometres)
is generated as hi,k = 10−σPL /10 h̃i,k , where σPL is the path
loss (PL) in dB and h̃i,k ∼ CN (0, INt ) represents small-scale
effects. The other parameters are given by Table I. Without
loss of generality, R̄i,k ≡ R̄ is set. The numerical results are
obtained using the parser YALMIP [9].
We compare the performance of the proposed FTbased beamforming scheme with five other beamforming schemes: (i) “Conventional DL,” under which the
problem
is formulated similarly as:
P of ST maximization
′
′
maxw (i,k)∈M Ri,k (w) s.t. Ri,k (w) ≥ R̄i,k , (i, k) ∈

50

Proposed FT
NOMA
FT + NOMA in both zones
FT + NOMA in zone-1
FT + NOMA in zone-2
Conventional DL

40

30

20
20

25
30
35
max
Transmit power at the BS, Pbs
[dBm]

40

(a) R̄ = 0 bps/Hz.
60
Average sum throughput (bps/Hz)

(κ)
Ri,k (wi , αi )

Value
-174 [dBm/Hz]
128.1 + 37.6log10 (di,k ) [dB]
500 [m]
200 [m]
> 10 [m]

60
Average sum throughput (bps/Hz)

Algorithm 1 QP-based path-following algorithm for ST maximization problem (4)
Initialization: Iterate (10) for an initial feasible point
(w(0) , α (0) ). Set κ := 0
1: repeat
2:
Solve convex quadratic program (9) to obtain the optimal
solution: (w(κ+1) , α (κ+1) ).
3:
Set κ := κ + 1.
4: until Convergence

3

50
40

Proposed FT
NOMA
FT + NOMA in both zones
FT + NOMA in zone-1
FT + NOMA in zone-2
Conventional DL

30
20
10
20

25
30
35
max
Transmit power at the BS, Pbs
[dBm]

40

(b) R̄ = 1 bps/Hz.
Fig. 2.

max .
Average sum throughput versus Pbs

max
M, kw1 k2 + kw2 k2 ≤ Pbs
under the definition (3); (ii)
“NOMA”: each UE in zone-1 is paired with an UE in zone2 according to the clustering algorithm in [10] to create a
virtual cluster. In each cluster, both UEs decode the signal
intended for the UE in zone-2 and then the UE in zone-1
processes successive interference cancellation (SIC) to cancel
the interference of the UE in zone-2 in decoding its own signal;
(iii) “FT + NOMA in both zones”: under FT, NOMA is zonewide adopted; (iv) “FT + NOMA in zone-1”: under FT, NOMA
is adopted only in zone-1; and (v) “FT + NOMA in zone-2”:
under FT, NOMA is adopted only in zone-2. The reader is
referred to [6, Sec. V] for beamforming under NOMA, which
is used in these five schemes. The computational complexity of
each iteration in NOMA is similar to that of (9). Note that the
performance of NOMA-based beamforming can be improved
by involving more UEs in virtual clusters [6, Sec VI] but the
UEs’ privacy is more compromised. On average, Algorithm 1
requires about 10 iterations for convergence.
Fig. 2 plots the average achievable ST versus the transmit
max
power Pbs
for Nt = 5. For R̄ = 0 bps/Hz shown in Fig. 2(a),
one can see that the ST of the proposed FT-based beamforming
is higher than that achieved by the other schemes in the high
transmit power region. On the other hand, the conventional DL
outperforms NOMA and FT+NOMA schemes for high transmit

4

IEEE COMMUNICATIONS LETTERS

Average sum throughput (bps/Hz)

60

40

20

0.4
0.6
0.8
1
Throughput threshold, R̄ [bps/Hz]

1.2

max = 30 dBm.
Average sum throughput versus R̄ for Pbs

Average max-min throughput (bps/Hz)

Fig. 4.

V. C ONCLUSIONS
The paper has proposed a fractional time-based beamforming
scheme at a base station to serve two groups of users, which
is able to improve the network throughput while preserving
the information privacy for the users. Accordingly, a pathfollowing computational procedure for the joint design of
fractional times and beamforming vectors to maximize the
network throughput has been developed. Extensive simulations
have been provided to demonstrate the superior performance
of the proposed scheme over the exiting schemes.

30

10
0.2

Fig. 3.

others schemes.

Proposed FT
NOMA
FT + NOMA in both zones
FT + NOMA in zone-1
FT + NOMA in zone-2
Conventional DL

50

5
4
3

Proposed FT
NOMA
FT + NOMA in both zones
FT + NOMA in zone-1
FT + NOMA in zone-2
Conventional DL

A PPENDIX
For all x > 0, x̄ > 0, y > 0, ȳ > 0, t > 0, and t̄ > 0, it is
true that
ln(1 + x2 /y)
y
(12)
≥ a − b 2 − ct
t
x
y
− ct
(13)
≥ a−b
x̄(2x − x̄)

2
1
0
20

25
30
35
max
Transmit power at the BS, Pbs
[dBm]

over the trust region

40

2x − x̄ > 0,

max .
Average max-min throughput versus Pbs

power by making the ST concentrated at those UEs with good
channel conditions. Apparently, NOMA does not look efficient
when there is no UEs’ QoS requirement. In Fig. 2(b) with R̄ =
1 bps/Hz, the conventional DL performs worse than NOMA
max
at high Pbs
. The low STs of FT+NOMA-based schemes are
probably attributed to the fact that NOMA is more efficient by
max
exploiting their channel differentiation. Increasing Pbs
leads
to a remarkable gain in ST by the proposed FT compared with
the other schemes. In addition, the proposed FT is feasible in
max
all range of Pbs
while the other schemes cannot offer such
max
high QoS to UEs at low Pbs
.
The plot of the ST versus QoS requirement threshold R̄ ∈
[0.2, 1.2] bps/Hz is illustrated in Fig. 3. We can observe
that the proposed FT-based beamforming performs quite well
and only slightly degrades when R̄ increases. The performance
gap between the proposed FT-based beamforming and other
schemes substantially increases by increasing R̄. It is expected
because the proposed FT-based beamforming can tune the
power allocation in meeting zone-2 UEs’ QoS requirements
without causing interference to the zone-1’s UEs.
Next, we look for the max-min UE throughput optimization
problem
max min Ri,k (wi , τi )
w,τ (i,k)∈M

s.t. (4c), (4d)

(11)

which can be addressed similarly by solving the convex program (10) (with R̄i,k ≡ 1) at each iteration.
Fig. 4 plots the average UEs’ worst throughput as a function
max
of Pbs
. It shows that the worst throughput achieved by the
conventional DL and FT + NOMA-based schemes is saturated
max
once Pbs
is beyond a threshold. It also reveals that the
proposed FT-based beamforming and FT + NOMA in zone1 achieve worst throughput that is higher than that achieved by

(14)
2

where a = 2[ln(1+d)]/t̄+d/t̄(d+1) > 0, b = d /t̄(d+1) > 0,
c = [ln(1 + d)]/t̄2 > 0, d = x̄2 /ȳ. Inequality (12) follows
from [11] while inequality (13) is obtained by using x2 ≥
x̄(2x − x̄) > 0 over the trust region (14).
R EFERENCES
[1] X. Shang, B. Chen, and H. V. Poor, “Multiuser MISO interference
channels with single-user detection: Optimality of beamforming and the
achievable rate region,” IEEE Trans. Inf. Theory, vol. 57, no. 7, pp. 4255–
4273, July 2011.
[2] G. Caire and S. Shamai, “On the achievable throughput of a multiantenna
Gaussian broadcast channel,” IEEE Trans. Inf. Theory., vol. 49, no. 7, pp.
1691 – 1706, July 2003.
[3] Z. Ding, R. Schober, and H. V. Poor, “A general MIMO framework for
NOMA downlink and uplink transmission based on signal alignment,”
IEEE Trans. Wireless Commun., vol. 15, no. 6, pp. 4483–4454, June
2016.
[4] Z. Ding, Y. Liu, J. Choi, Q. Sun, M. Elkashlan, C. L. I, and H. V. Poor,
“Application of non-orthogonal multiple access in LTE and 5G networks,”
IEEE Commun. Mag., vol. 55, no. 2, pp. 185–191, Feb. 2017.
[5] Z. Ding, F. Adachi, and H. V. Poor, “The application of MIMO to nonorthogonal multiple access,” IEEE Trans. Wireless Commun., vol. 15,
no. 1, pp. 537–552, Jan. 2016.
[6] V.-D. Nguyen, H. D. Tuan, T. Q. Duong, H. V. Poor, and O.-S. Shin,
“Precoder design for signal superposition in MIMO-NOMA multicell
networks,” IEEE J. Select. Areas Commun., 2017, to appear, available
at http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=7974739.
[7] A. Wiesel, Y. Eldar, and S. Shamai, “Linear precoding via conic optimization for fixed MIMO receivers,” IEEE Trans. Signal Process., vol. 54,
no. 1, pp. 161 – 176, Jan. 2006.
[8] B. R. Marks and G. P. Wright, “A general inner approximation algorithm
for nonconvex mathematical programms,” Operations Research, vol. 26,
no. 4, pp. 681–683, July 1978.
[9] J. Löfberg, “YALMIP: A toolbox for modeling and optimization in
MATLAB,” in Proc. 2004 IEEE Inter. Symposium on Computer Aided
Control Systems Design, Sept. 2004, pp. 284–289.
[10] B. Kimy, S. Lim, H. Kim, S. Suh, J. Kwun, S. Choi, C. Lee, S. Lee, and
D. Hong, “Non-orthogonal multiple access in a downlink multiuser beamforming system,” in Proc. IEEE Military Commun. Conf. (MILCOM),
Nov. 2013, pp. 1278–1283.
[11] V.-D. Nguyen, T. Q. Duong, H. D. Tuan, O.-S. Shin, and H. V. Poor,
“Spectral and energy efficiencies in full-duplex wireless information and
power transfer,” IEEE Trans. Commun., vol. 65, no. 5, pp. 2220–2233,
May 2017.

