1

On Non-Orthogonal Multiple Access with
Finite-Alphabet Inputs in Z-Channels

arXiv:1707.08305v1 [cs.IT] 26 Jul 2017

Zheng Dong, He Chen, Jian-Kang Zhang, and Lei Huang

Abstract—This paper focuses on the design of non-orthogonal
multiple access (NOMA) in a classical two-transmitter tworeceiver Z-channel, wherein one transmitter sends information
to its intended receiver from the direct link while the other
transmitter sends information to both receivers from the direct
and cross links. Unlike most existing designs using (continuous)
Gaussian input distribution, we consider the practical finitealphabet (i.e., discrete) inputs by assuming that the widelyused quadrature amplitude modulation (QAM) constellations
are adopted by both transmitters. To balance the error performance of two receivers, we apply the max-min fairness design
criterion in this paper. More specifically, we propose to jointly
optimize the scaling factors at both transmitters, which control
the minimum Euclidean distance of transmitting constellations,
to maximize the smaller minimum Euclidean distance of two
resulting constellations at the receivers, subject to an individual
average power constraint at each transmitter. The formulated
problem is a mixed continuous-discrete optimization problem and
is thus intractable in general. By resorting to the Farey sequence,
we manage to attain the closed-form expression for the optimal
solution to the formulated problem. This is achieved by dividing
the overall feasible region of the original optimization problem
into a finite number of sub-intervals and deriving the optimal
solution in each sub-interval. Through carefully observing the
structure of the optimal solutions in all sub-intervals, we obtain
compact and closed-form expressions for the optimal solutions
to the original problem in three possible scenarios defined by
the relative strength of the cross link. Simulation studies are
provided to validate our analysis and demonstrate the merits of
the proposed design over existing orthogonal or non-orthogonal
schemes.
Index Terms—Non-orthogonal multiple access (NOMA), Zchannel, finite-alphabet inputs, quadrature amplitude modulation, max-min fairness, Farey sequence.

I. I NTRODUCTION
Multiple access technologies have been playing an important role in determining the performance of each generation of
mobile communication systems. Based on how the resources
are allocated to users, multiple access technologies can generally be categorized into two types: orthogonal multiple access
(OMA) and non-orthogonal multiple access (NOMA) [1,
Ch. 6]. The current generation of cellular networks, known as
Z. Dong is with Shenzhen University, Shenzhen 518060, China, and
also with McMaster University, Hamilton, ON L8S 4K1, Canada (e-mail:
dongz3@mcmaster.ca).
H. Chen is with The University of Sydney, Sydney, NSW 2006, Australia
(e-mail: he.chen@sydney.edu.au).
J.-K. Zhang is with McMaster University, Hamilton, ON L8S 4K1, Canada
(e-mail: jkzhang@mail.ece.mcmaster.ca).
L. Huang is with Shenzhen University, Shenzhen 518060, China (e-mail:
lhuang@szu.edu.cn).

4G, and all previous generations have primarily adopted the
OMA technologies, which include frequency-division multiple
access (FDMA) for 1G, time-division multiple access (TDMA)
for 2G, code-division multiple access (CDMA) for 3G, and
orthogonal frequency-division multiple access (OFDMA) for
4G [2]. In these OMA schemes, the resource is partitioned into
orthogonal blocks in time/frequency/code domain, and each
resource block is then assigned to one single user exclusively.
In this sense, there is no inter-user interference in OMA,
leading to low-complexity receiver and scheduling algorithms.
Moreover, after the resource allocation, the multiple-user problem is divided into several point-to-point problems such that
the well-established single-user encoder/decoder techniques
can be directly applied. However, early information-theoretic
studies showed that compared with NOMA, OMA has lower
spectral efficiency as it normally cannot achieve the multi-user
channel capacity region [3]. Besides, OMA is not scalable as
the total number of orthogonal resources and their granularity
strictly limit the maximum number of served users.
Different from OMA, NOMA exploits the power domain
to multiplex multiple users together such that they can be
served in the same time/frequency/code resources [2]–[6]. As
such, with proper multi-user detection techniques to deal with
the inter-user interference at the receiver side (e.g., successive interference cancellation (SIC) [7]), NOMA is capable
of achieving improved spectral efficiency and serving much
more users simultaneously. In fact, the uplink and downlink
versions of NOMA, well-known as multiple access channel
(MAC) and broadcast channel (BC) respectively, have been
intensively investigated for several decades in the information
theory community, see, e.g., [8]–[11]. However, due to the
high complexity of interference cancellation, these studies
mainly lied in the theoretical aspects and their results were
not implemented in practical communication systems. With
the fast advances of hardware, the implementation of NOMA
with interference cancellation becomes more affordable and
feasible. Actually, NOMA has been regarded as a key enabling
technology to meet the unprecedented requirements of 5G
wireless networks due to its significant network throughput
gain and great potential to support massive connectivity, low
latency and user fairness [2], [5], [7], [12]–[17]. Furthermore,
a two-user downlink scenario of NOMA, termed multiuser
superposition transmission (MUST), has been incorporated in
the 3rd Generation Partnership Project (3GPP) Long Term
Evolution-Advanced (LTE-A) [18], [19].
Most conventional information-theoretic and recent studies
on NOMA adopted the assumption of Gaussian input distribution [5], [8], [9], [20]–[25]. Although the designs with

2

Gaussian signaling can approach most of the known capacity
inner bounds, such as in [8]–[11], their direct implementation
in practical communication systems may lead to significant
performance loss [26]. Moreover, Gaussian signaling will
require unaffordable encoding and decoding efforts, which
could lead to extremely high hardware cost, huge storage
capability, high computational complexity, and long delay.
Therefore, Gaussian inputs could arguably be infeasible for
current hardware and it acts mostly as the theoretical benchmark. The inputs of practical wireless systems are actually
drawn from finite constellations, such as phase shift keying (PSK) modulation or quadrature amplitude modulation
(QAM), which are essentially different from the continuous
Gaussian inputs. When it comes to a NOMA system with finite
input constellations, the key design challenge is to guarantee
that each user’s codeword can be uniquely decoded from their
sum signal at the receiver side [27]–[29]. For the two-user
MAC with finite-alphabet inputs, a constellation rotation (CR)
scheme and a constellation power allocation (CPA) scheme
were proposed in [30] and [31] to construct an unambiguous
sum constellation at the receiver, respectively. This is achieved
by strategically introducing certain angle of rotation between
the input constellations in the CR scheme and appropriately
controlling the transmit power of each user in the CPA scheme.
The results in [30] and [31] have been extended to various multiple-antenna scenarios, see, e.g., [32] and references
therein. The aforementioned NOMA designs were primarily
for the PSK modulations by utilizing its circular structure. The
studies on NOMA with QAM, another practical modulation
scheme that has been widely adopted in cellular systems due
to its higher spectral efficiency, are quite limited. Very recently,
the mutual information were used as the performance metric to
optimize NOMA systems with QAM in [33], [34]. However,
the optimal NOMA designs in [33], [34] were achieved by
numerical approaches with high computational complexities.
We also notice that the existing studies on NOMA mainly
focused on MAC and BC, due to their wide applications
in centralized systems like cellular networks. Also there are
some initial efforts considering the NOMA design for the
interference channel (IC) [3]. With recent advances in noncentralized networks (e.g., wireless sensor networks and ad
hoc networks), Z-channel (ZC) was proposed in [35] and
attracted considerable attention in the past decade [36]–[42].
As a special case of the classical two-user IC [43], a two-user
ZC consists of two transmitters and two receivers, wherein
one transmitter sends information to its intended receiver from
the direct link without causing interference to the unintended
receiver, while the other transmitter sends information to both
receivers from the direct and cross links. In this sense, the
ZC is a general channel model that includes MAC and BC
as special cases. The ZC is also closely related to the Zinterference channel (ZIC) [44]–[47], wherein each transmitter
transmits information to its corresponding receiver only from
the direct link and the received signal from the cross link carries no desired information, and thus is treated as interference
at the receiver side. It is also worth emphasizing that there are
two messages transmitted from the two direct links in ZIC,
while three messages are sent via the two direct links and the

cross link in ZC. Both ZC and ZIC are proper models for the
multi-cell downlink transmission, where one user is located
near to the cell edge and thus can receive signals from both
transmitters, while the other user is near the cell center and
suffers from no interference. Another example corresponds to
the two-user IC, where one cross link is blocked by obstacles
with large pathloss such as tall buildings or thick walls, while
the other user is still exposed to interference [45], [48]. Despite
of the existing great efforts on ZC, to our best knowledge, the
design of NOMA with practical finite-alphabet inputs in ZC
is still an open problem in the literature.
Motivated by this gap, in this paper we concentrate on the
practical design of NOMA with QAM in a two-transmitter
two-receiver ZC. It is worth emphasizing that the design of
NOMA with QAM is much more challenging than that for
PSK modulation. This is mainly because the unambiguity of
sum QAM constellations is much harder to maintain since its
signal points are distributed more evenly and there is a higher
probability that more than one signal points coincide or close
to each other on the sum constellation. The main contributions
of this paper can be summarized as follows:
1) We, for the first time, develop a practical NOMA
framework with QAM and max-min fairness in ZC.
In our framework, we optimize the scaling factors of
both transmitters, which adjust the minimum Euclidean
distance of the transmitting constellations, to maximize
the smaller minimum Euclidean distance among the
resulting constellations at both receivers subject to an
individual power constraint on each transmitter. Through
our design, the average error performance of both transmitters in the considered ZC can be minimized with
good user fairness, which is fundamentally different
from the existing designs that mainly focused on the
channel capacity maximization.
2) The formulated optimization problem is shown as a
mixed continuous-discrete optimization problem, which
is challenging to solve in general. By carefully observing
the features of the formulated problem, we realize that
the Farey sequence (also known as Farey series) [49] can
be applied to resolve the problem. More specifically, by
taking the advantage of Farey sequence associated with
the finite-alphabet, we strategically partition the entire
feasible region of the original optimization problem into
a finite number of sub-intervals and attain the closedform solution in each sub-interval. Then, by a careful
observation on the structure of the solutions in all subintervals, the overall solution is obtained in a compact
closed-form for three complementary scenarios divided
by the relative strength of the cross link.
3) We verify the correctness of the analytical results by
conducting simulations in both deterministic and random
fading channels. Simulation results show that the sum
constellation at the receiver side is still a regular QAM
constellation with a larger size for most scenarios, but
could be a hierarchical QAM with two-resolution [50],
e.g., when the cross link is very strong relative to the
direct link. We adopt the bit error rate (BER) as the

3

(a)
(b)
Fig. 1. (a) Two-user complex Gaussian ZC. (b) Two-user real Gaussian ZC.

performance metric to compare the proposed NOMA
design with the existing OMA and NOMA schemes
under random fading channels. The comparison illustrates that our scheme can achieve a significant lower
BER performance than the benchmark schemes, which
validates the effectiveness of our design.
II. S YSTEM M ODEL OF C OMPLEX G AUSSIAN ZC WITH
QAM C ONSTELLATIONS
We consider a two-user complex Gaussian ZC consisting of
two transmitters S1 and S2 , and two receivers D1 and D2 , as
depicted in Fig. 1(a). We consider that each node is equipped
with a single antenna and works in a half-duplex mode. As
per the ZC, only the cross link between S2 and D1 is assumed
to be available. Moreover, S1 sends one unicasting message
x1 to D1 , while S2 transmits one multicasting message x2
to both D1 and D2 . S1 and S2 transmit their messages
simultaneously using the same frequency band. The equivalent
complex baseband signals observed at D1 and D2 can be
given, respectively, by the following equations:
z1 = h11 x1 + h21 x2 + ξ1 ,
z2 = h22 x2 + ξ2 ,

(1a)
(1b)

where hkℓ ∈ C, k, ℓ ∈ {1, 2} denotes the complex channel
coefficients from transmitter Sk to receiver Dℓ . Hereafter, we
call h11 and h22 the direct links, while h21 is referred to as
the cross link. In line with [35]–[39], [41], [42], [45], all the
channel links are assumed to be known perfectly at all the
terminals. The additive noise processes ξ1 , ξ2 ∼ CN (0, 2σ 2 )
are independent and identically distributed (i.i.d.) over time
and are assumed to be circularly symmetric complex Gaussian
(CSCG). Note that the case with different noise levels at
receivers can be incorporated into our model by scaling
operations. We suppose that QAM constellations are used
by both transmitters since it is more spectrally efficient than
other frequently-used modulation schemes such as phase-shift
keying (PSK), and is also relatively easy to implement [51,
Ch. 5.3.3]. A predefined average transmitted power constraint
is imposed to both transmitters1 , i.e., E[|x1 |2 ] ≤ P1 and
E[|x2 |2 ] ≤ P2 . In this paper, the system signal-to-noise ratio
(SNR) is defined by ρ = 2σ1 2 .
For analytical simplicity, we decompose the considered
complex Gaussian ZC given in (1) into two parallel real scalar
Gaussian ZCs [35], [36], [45], which are called the in-phase
and quadrature components, respectively. We note that this
method was commonly used in the study of IC and the real
1 Our design can also be generalized to the case with peak power constraint
straightforwardly.

IC was studied directly in [48], [52]–[54]. Actually, designing
two-dimensional QAM constellations is an extremely challenging problem even for two-user MAC, see e.g., [30], [31]
and references therein. In this paper, instead of designing the
two-dimensional QAM constellations directly, we propose a
practical design that decomposes the complex Gaussian ZC
into two parallel real scalar Gaussian ZCs such that we can
split the two-dimensional QAM constellation into two onedimensional PAM constellations. In fact, although we use
complex baseband representation in (1), the actual modulated
and demodulated signals are all real since the oscillator at
the transmitter can only generate real sinusoids rather than
complex exponentials, and the channel just introduces amplitude and phase distortion to the transmitted signals [51,
Ch. 2.2]. By this means, the original two-dimensional QAM
constellation can be split into two one-dimensional pulse amplitude modulation (PAM) constellations for both the in-phase
and quadrature components. Mathematically, for the complex
Gaussian ZC described in (1), the in-phase and quadrature
components can be attained by rotating x1 and x2 according
to the instantaneous channel coefficients to compensate for
the phase offset, and then taking the real and imaginary parts,
respectively. First of all, we note that (1) is equivalent to
z1 = |h11 | exp(j arg(h11 ))x1 + |h21 | exp(j arg(h21 ))x2 + ξ1

h21 
exp j arg
z2 = |h22 | exp(j arg(h21 ))x2
h22

h21 
ξ2 .
(2)
+ exp j arg
h22
Now, we set

h21  
y1 = Re(z1 ), y2 = Re exp(j arg
)z2 ,
h22

w1 s1 = Re exp(j arg(h11 ))x1 ,

w2 s2 = Re exp(j arg(h21 ))x2 ,
h21  
)ξ2 ;
n1 = Re(ξ1 ), n2 = Re exp(j arg
h22

h21  
y1′ = Im(z1 ), y2′ = Im exp(j arg
)z2 ,
h22

w1′ s′1 = Im exp(j arg(h11 ))x1 ,

w2′ s′2 = Im exp(j arg(h21 ))x2 ,
h21  
)ξ2 ,
n′1 = Im(ξ1 ), n′2 = Im exp(j arg
h22

(3a)
(3b)
(3c)
(3d)
(3e)
(3f)
(3g)
(3h)

where Re(·) and Im(·) are the real and imaginary parts of the
complex number, respectively. We also assume that s1 ∈ AM1 ,
s′1 ∈ AM1′ , sent by S1 , and s2 ∈ AM2 , s′2 ∈ AM2′ , transmitted
by S2 , are the information-bearing symbols, which are drawn
from standard PAM constellation with equal probability, in
M/2
which AM , {±(2k − 1)}k=1 is a M -ary PAM constellation
set. Moreover, the scaling factors w1 , w2 , w1′ , and w2′ are
real positive scalars that determine the minimum Euclidean
distance of the corresponding PAM constellation set.

From (3), we obtain exp j arg(h
11 ) x1 = (w1 s1 +

w1′ s′1 j) ∈ Q1 and exp j arg(h21 ) x2 = (w2 s2 + w2′ s′2 j) ∈

4

(4a)

if M1 = M1′ and M2 = M2′ , then the two sub-channels are
identical. It is worth noting that, similar design for BC or MAC
can be included as a special case of our proposed design for
the considered ZC.

(4b)

A. Problem Formulation

Q2 , where
Q1 , {±w1 (2k − 1) ± w1′ (2ℓ − 1)j :
k = 1, . . . , M1 /2, ℓ = 1, . . . , M1′ /2},
Q2 , {±w2 (2k − 1) ±

w2′ (2ℓ

− 1)j :

k = 1, . . . , M2 /2, ℓ = 1, . . . , M2′ /2},

are M1 M1′ - and M2 M2′ -ary QAM constellations, respectively.
If w1 = w1′ and w2 = w2′ , we call Q1 and Q2 symmetric QAM
constellations. Otherwise, we are using unsymmetric signaling [45], [48], [55]. In addition, n1 , n2 , n′1 , n′2 ∼ N (0, σ 2 )
are i.i.d. real additive white Gaussian since the complex noise
terms are assumed to be CSCG. Then, the in-phase and
quadrature sub-channels of (1), as illustrated in Fig. 1(b), can
be reformulated by
y1 = |h11 |w1 s1 + |h21 |w2 s2 + n1 ,
y2 = |h22 |w2 s2 + n2 ,

(5a)

y1′
y2′

(5b)

=
=

|h11 |w1′ s′1
|h22 |w2′ s′2

+
+

|h21 |w2′ s′2
n′2 .

+

n′1 ,

The transmitted signals over both subchannels should still be
subject to average power constraints, i.e., E[w12 |s1 |2 ] ≤ p1 ,
E[w22 |s2 |2 ] ≤ p2 , E[w1′2 |s′1 |2 ] ≤ p′1 , E[w2′2 |s′2 |2 ] ≤ p′2 such that
p1 +p′1 = P1 and p2 +p′2 = P2 . The following power allocation
among the in-phase and quadrature components is normally
performed to balance the minimum Euclidean distance of the
two PAM constellations [51, Ch. 6.1.4], i.e.,
(M 2 − 1)P2
(M12 − 1)P1
, p2 = 2 2
,
2
′2
M1 + M1 − 2
M2 + M2′2 − 2
(M2′2 − 1)P2
(M ′2 − 1)P1
′
,
p
=
.
p′1 = 2 1
M1 + M1′2 − 2 2
M22 + M2′2 − 2
p1 =

(6a)
(6b)

It can be observed that, if square-QAM constellations are used
at both transmitters with M1 = M1′ and M2 = M2′ , we have
p1 = p′1 = P1 /2 and p2 = p′2 = P2 /2.
An important problem for the considered ZC is that for
any given QAM constellation sizes of both messages, how
to optimize the values of scaling coefficients w1 , w2 , w1′
and w2′ to minimize the average error probability at both
receivers, subject to the individual average power constraint
at both transmitters. By leveraging the decomposable property
of the complex Gaussian ZC and the symmetry of the two
subchannels, we can simply focus on the design for one of
the two real Gaussian ZCs with PAM constellation sets, which
will be elaborated in the next section2 .
III. T HE C ONSTELLATION D ESIGN
G AUSSIAN ZC

FOR THE

R EAL

In this section, we consider the constellation design problem, i.e., finding the optimal values of w1 and w2 for the
in-phase real Gaussian ZC characterized by (5a). The optimal
solution to the quadrature component can be obtained in a
similar fashion and hence omitted for brevity. In particular,
2 It should be pointed out that this design is a practical but not necessarily
optimal approach, which has been widely adopted in practice [36], [45], [48],
[52]–[57].

As the first effort towards the design of NOMA with finitealphabet inputs in ZC, in this paper we concentrate on the case
that M1 = M2 = M and M1′ = M2′ = M ′ . As a result, we
M/2
2
2
have s1 , s2 ∈ AM = {±(2k − 1)}k=1 . As E[wq
1 |s1 | ] ≤ p1

1
and E[w22 |s2 |2 ] ≤ p2 , we thus have 0 < w1 ≤ M3p
2 −1 and
q
2
0 < w2 ≤ M3p
2 −1 .
In our scheme, the transmitted signal from S1 and S2
are superimposed together at D1 , which is inherently a nonorthogonal transmission. In line with [36], we use a joint
decoding3 at the receiver D1 since the error performance
is dominated by the minimum Euclidean of the resulting
sum-constellation at D1 . We assume that each receiver uses
a coherent maximum-likelihood (ML) detector to estimate
the transmitted signals in a symbol-by-symbol fashion4 . For
receivers D1 and D2 , the estimated signals can be expressed
as

(ŝ1 , ŝ2 ) = arg min

(s1 ,s2 )

y1 − (|h11 |w1 s1 + |h21 |w2 s2 ) ,

ŝ′2 = arg min y2 − |h22 |w2 s2 .
s2

By applying the nearest neighbor approximation
method [51, Ch. 6.1.4] at high SNRs for the ML receiver, the
average error rate is dominated by the minimum Euclidean
distance of the received constellation points owing to the
exponential decaying of the Gaussian distribution. To balance
the error performance of both receivers, in this paper, we
aim to devise the optimal value of w1 and w2 by applying
the max-min fairness criterion on the minimum Euclidean
distance of the received signal constellation points among
both receivers.
The Euclidean distance between two received signals
y1 (s1 , s2 ) and y1 (s̃1 , s̃2 ) at D1 and that between y2 (s2 )
and y2 (s̃2 ) at D2 for the transmitted signal vectors (s1 , s2 )
and (s̃1 , s̃2 ) at S1 and S2 in the noise-free case are given,
respectively, by
|y1 (s1 , s2 ) − y1 (s̃1 , s̃2 )| = |h11 |w1 (s1 − s̃1 )
− |h21 |w2 (s̃2 − s2 ) ,
|y2 (s2 ) − y2 (s̃2 )| = |h22 |w2 |s2 − s̃2 |.

(7a)
(7b)

3 We note that, for ZIC, the joint decoder used by D may not necessarily
1
be the most efficient one. Instead, we should use a treat-interference-as-noise
(TIN) receiver when the channel gain of the cross link is very low and use
a successive-interference-cancellation (SIC) receiver when the channel gain
of cross link is very strong compared with the direct link [36], [43], [44].
In general, a joint decoder can be used when the cross link is moderately
strong [36], [43], which will result in a similar design as our case. However,
how to extend our design to ZIC with finite-alphabet input is still an open
problem and has been left as a future work.
4 Since we are doing a symbol-by-symbol detection, the decoding complexity is O(M 2 ) with M being the PAM constellation size of S1 and S2 ,
respectively. Although we can use the message passing algorithm (MPA) [58]
to further decrease the decoding complexity, however, our method is feasible.

5

Note that s1 , s̃1 , s2 and s̃2 are all odd numbers, and we thus
can let s1 −s̃1 = 2n and s̃2 −s2 = 2m, in which m, n ∈ ZM−1
with ZM−1 , {0, ±1, . . . , ±M − 1} denoting the set containing all the possible differences. Similarly, we also define
Z2M−1 , {(a, b) : a, b ∈ ZM−1 }, NM−1 , {0, 1, · · · , M − 1}
and N2M−1 , {(a, b) : a, b ∈ NM−1 }. From the above
definition, (s1 , s2 ) 6= (s̃1 , s̃2 ) is equivalent to (m, n) 6= (0, 0).
Here, by (m, n) 6= (0, 0), we mean that m 6= 0 or n 6= 0. To
proceed, we define
1
d1 (m, n) = |y1 (s1 , s2 ) − y1 (s̃1 , s̃2 )|
2
= |h11 |w1 n − |h21 |w2 m , for (m, n) ∈ Z2M−1 \ {(0, 0)},
1
d2 (m) = |y2 (s2 ) − y2 (s̃2 )|
2
= |h22 |w2 |m|, for m ∈ ZM−1 \ {0}.
We are now ready to formally formulate the following maxmin optimization problem:
Problem 1 (Optimal Design of NOMA in real scalar ZC
with PAM constellation): Find the optimal value of (w1∗ , w2∗ )
subject to the individual average power constraint such that the
minimum value of the minimum Euclidean distances of the
received signal constellation points over both received signals
is maximized, i.e.,
n
(w1∗ , w2∗ ) = arg max min
min
d1 (m, n),
2
(w1 ,w2 )

(m,n)∈ZM −1 \{(0,0)}

|

r

{z
T1

min

m∈ZM −1 \{0}

|

{z
T2

}

(9a)

r

3p2
. (9b)
M2 − 1

Note that the inner optimization problem of finding the
minimum Euclidean distances is a discrete one, while the outer
optimization problem on (w1 , w2 ) is a continuous problem.
In other words, Problem 1 is a mixed continuous-discrete
optimization problem and it is in general hard to solve. To
the best of our knowledge, only numerical solutions to such
kind of problems are available in the literature, see e.g., [30]–
[32] and references therein.
To optimally and systematically solve this problem, we
now develop a novel framework based on the Farey sequence
(also known as Farey series) [49], which can divide the
entire feasible region of (w1 , w2 ) into a finite number of
mutually exclusive sub-intervals. Then for each sub-interval,
the formulated optimization problem can be solved optimally
with a closed-form solution, and subsequently the overall
maximum value of Problem 1 can be attained by taking the
maximum value of the objective function among all the subintervals.
For the inner optimization problem of T2 given in (9a), it
can be observed that
s.t. 0 < w1 ≤

min

m∈ZM −1 \{0}

3p1
and 0 < w2 ≤
M2 − 1

d2 (m)
}

o

d2 (m) =

min

m∈ZM −1 \{0}

|h22 |w2 |m|

= |h22 |w2 , with m = 1.

However, for T1 , we have
min

(m,n)∈Z2M −1 \{(0,0)}

=

min

d1 (m, n)

(m,n)∈Z2M −1 \{(0,0)}

|h11 |w1 n − |h21 |w2 m . (10)

We should point out that the closed-form solution to the
optimal (m, n) is not trivial, since the solution depends on
the values of |h11 | and |h22 |, which can span the whole
positive real axis. Moreover, the value of w1 and w2 can not
be determined beforehand. Actually, the problem in (10) is
essentially equivalent to finding a real rational number with
finite order to approximate a real irrational number as closely
as possible. This naturally leads us to resorting to the Farey
sequence, which particularly plays a critical role in solving
such kind of problems [49]. In the subsequent section, we
will introduce the definition and some important properties of
Farey sequence.
B. Farey Sequence
The Farey sequence characterizes the relationship between
two positive integers and the formal definition is given as
follows:
Definition 1 (Farey sequence): [49] The Farey sequence FK
of order K is the ascending sequence of irreducible fractions
between 0 and 1 whose denominators are less than or equal
to K.

|FK |
By the definition, FK = abkk k=1 is a sequence of fractions
bk
ak such that 0 ≤ bk ≤ ak ≤ K and hak , bk i = 1 arranged in
an increasing order, where ha, bi denotes the largest
P common
divider of non-negative integers a, b. |FK | = 1 + K
m=1 ϕ(m)
is the cardinality of FK with ϕ(·) being the Euler’s totient
function [49]. Some examples of Farey sequences are given
as follows:
Example 1: F5 is the ordered sequence
0 1 1 1 2 1 3 2 3 4 1
.
, , , , , , , , , ,
1 5 4 3 5 2 5 3 4 5 1

It can be observed that each Farey sequence begins with
number 0 (fraction 01 ) and ends with 1 (fraction 11 ). The series
of breakpoints after 11 is the reciprocal version of the Farey
sequence. We call the Farey number sequence together with
its reciprocal version as the extended Farey sequence which is
formally defined as follows:
Definition 2 (Extended Farey sequence): The extended Farey
sequence SK of order K is the sequence of ascending irreducible fractions, where the maximum value of the numerator
and denominator do not exceed K.

|SK |
From the definition, we have SK = abkk k=1 with
P
hak , bk i = 1 and |SK | = 1 + 2 K
m=1 ϕ(m). We have the
following example:
Example 2: S5 is the sequence
0 1 1 1 2 1 3 2 3 4 1 5 4 3 5 2 5 3 4 5 1
.
, , , , , , , , , , , , , , , , , , , ,
1 5 4 3 5 2 5 3 4 5 1 4 3 2 3 1 2 1 1 1 0
It can be observed that the extended Farey sequence starts with
number 0 (fraction 10 ) and end with ∞ (fraction 10 ).

6

The positive real axis can be divided by the extended
Farey
PK sequence SK into a finite number (i.e., |SK | − 1 =
2 m=1 ϕ(m)) of intervals. In this paper, we call the fractions
consisting of adjacent terms in the extended Farey sequence
as a Farey pair, and the interval between the Farey pair is
referred to as a Farey interval. We then have the Farey interval
set formally defined as follows:
Definition 3 (Farey interval set): A Farey interval set SK of
order K is the set containing all the Farey intervals generated
by the Farey pair of the extended Farey sequence SK .


 |SK |
By the definition, we have SK = abkk , abk+1
,
where
k=1
k+1
PK
|SK | = |SK | − 1 = 2 m=1ϕ(m). Note that, with a slight
b
abuse of notations, abkk , ak+1
denotes the interval between
k+1

rather than a vector, and this will be
end nodes abkk and abk+1
k+1
clear from the context.
Example 3: The Farey interval set S5 is the set given by
n 0 1   1 1   1 1 
 3 4   4 5   5 1 o
,
,...,
,
,
.
,
, ,
,
,
,
,
1 5
5 4
4 3
1 1
1 1
1 0
The Farey interval set
can be further divided into two
′

b′ b
L
subsets UK
= { ak′ , ak+1
∈ SK : b′k + b′k+1 ≥ L}
′
k
k+1
′′

′′
b
b
L
∈ SK : b′′k + b′′k+1 < L} for
and VK
= { ak′′ , ak+1
′′
k

n
∈ SK }, where SK
Proposition 1: Let F2K = {(m, n) : m
is the extended Farey number sequence of order K, and then
we have

min

(m,n)∈Z2K \{(0,0)}

C. The Minimum Euclidean Distance of the Received Signal
Constellation Points
We are now ready to solve the problem in (10) to find
the constellation point pairs (m, n) that have the minimum
Euclidean distance. To that end, we first have the following
preliminary propositions.

min

(m,n)∈F2K

d1 (m, n).


The proof is given in Appendix-A.

n1 n2
Proposition 2: Consider the Farey interval m
∈ SK ,
,
1 m2
|h21 |w2
n1
n2
n1 n2
with K ≥ 1, and m
<
.
Then,
for
∈
(
m2
|h11 |w1
m1 , m2 )
1
and d1 (m, n) = |h11 |w1 n − |h21 |w2 m , we have
n1 +n2
21 |w2
1) If |h
|h11 |w1 = m1 +m2 , then d1 (m1 , n1 ) = d1 (m2 , n2 );

n1 +n2
n1
21 |w2
2) If |h
<
∈
|h11 |w1
m1 , m1 +m2 , then d1 (m1 , n1 )
d1 (m2 , n2 );

n2
n1 +n2
21 |w2
3) If |h
>
∈
|h11 |w1
m1 +m2 , m2 , then d1 (m1 , n1 )
d1 (m2 , n2 ).

The proof can be found in Appendix-B.
n1 n2 n3 n4
Proposition 3: Consider m
, m2 , m3 , m4 ∈ SK , with K ≥
n2
n3 1
n4
n2 n3
n1
where m
,
form one
2, such that m1 < m2 < m3 < m
4
2 m3
Farey pair,
n2
n2 +n3
21 |w2
1) If |h
|h11 |w1 ∈ ( m2 , m2 +m3 ), then
min

k+1

L
L
L
L
L = 1, 2, . . . , 2K such that SK = UK
∪VK
and UK
∩VK
= ∅.
1
1
2K
In particular, UK = SK and VK = ∅ while UK = ∅ and
2K
VK
= SK .
Example 4: For the Farey interval set S5 , we have
n 1 3   3 2   2 3   3 4   4 1 
U54 =
, ,
,
,
,
,
,
,
,
,
2 5
5 3
3 4
4 5
5 1
1 5 5 4 4 3 3 5 5 2
,
,
,
,
,
,
,
,
,
,
1 4
4 3
3 2
2 3
3 1
 2 5   5 3   3 4   4 5   5 1 o
,
,
,
,
,
,
,
,
,
,
1 2
2 1
1 1
1 1
1 0
and
n 0 1   1 1   1 1   1 2   2 1 o
V54 =
,
, ,
,
,
.
,
,
,
,
1 5
5 4
4 3
3 5
5 2
We now review some elementary properties of Farey sequences [49] which are also true for extended Farey sequences.
n2
n1
and m
are two adjacent terms (Farey
Property 1: If m
1
2
n2
n1
< m
, then,
pairs) of SK with K ≥ 1, such that m
1
2
1) m1 n2 − m2 n1 = 1.

m1 +m2
n1 n2
m2 m1
2
2) mn11 +n

+m2 ∈ m1 , m2 , n1 +n2 ∈ n2 , n1 .
n2
n3
n1
,
and
are
three
consecutive
terms
Property 2: If m
m2
m3
1
n1
n2
n3
n2
of SK with K ≥ 1 such that m
< m
< m
, then, m
=
1
2
3
2
n1 +n3
.

m1 +m3
n1 n2 n3 n4
Property 3: Given K ≥ 2, we assume m
,
,
,
∈
1 m2 m3 m4
n2
n3
n4
n2
n3
n1
form
one
SK and m1 < m2 < m3 < m4 . If m2 and m
3
n2
n3
n2 +n4
3
Farey pair, then mn11 +n

+m3 ≤ m2 and m3 ≤ m2 +m4 .

d1 (m, n) =

(m,n)∈F2M −1

d1 (m, n) = d1 (m2 , n2 )
= |h21 |w2 m2 − |h11 |w1 n2 .

2) If

|h21 |w2
|h11 |w1

+n3
∈ ( mn22 +m
, n3 ), then
3 m3

min

(m,n)∈F2M −1

d1 (m, n) = d1 (m3 , n3 )
= |h11 |w1 n3 − |h21 |w2 m3 .


The proof is provided in Appendix-C.
Proposition 4: For the Farey interval set SM−1 =
 bk bk+1  |SM −1 |
ak+1
21 |
, if |h
ak , ak+1
|h22 | ≥ bk + bk+1 , then bk+1 +
k=1
+ak+1
≤ abkk +b
≤ abkk − bk|h|h2221| | .
k+1
The proof is given in Appendix-D.
|h22 |
bk+1 |h21 |



21 |w2
D. Optimal Solution to Problem 1 for |h
|h11 |w1 in Certain Farey
Interval
|h21 |w2
In this section, we solve Problem 1 by restricting |h
11 |w1
into a certain Farey interval where a closed-form solution
is attainable. We consider the Farey interval set SM−1
 b bk+1  |SM −1 |
k
where |SM−1 | =
given by SM−1 =
ak , ak+1
k=1
PM−1
21 |w2
2 m=1 ϕ(m). Now we consider the case |h
∈
|h11 |w1

b
bk
k+1
ak , ak+1 , k = 1, 2, . . . , |SM−1 | and we aim to find the
optimal (w1∗ (k), w2∗ (k)) such that
b b

n
k
k+1
g
= max min
min
d1 (m, n),
,
(w1 ,w2 )
ak ak+1
(m,n)∈F2M −1
o
min d2 (m)
(11a)
m∈Z\{0}
r
|h21 |w2
bk+1
3p1
bk
<
≤
, 0 < w1 ≤
s.t.
ak
|h11 |w1
ak+1
M2 − 1
r
3p2
and 0 < w2 ≤
. (11b)
M2 − 1

7

By applying the propositions in last subsections, we manage
to attain the following lemma in terms of the optimal solution
to problem (11).
Lemma 1: The optimal solution to (11) is given as follows.
|h21 |
1) If |h
≤ bk + bk+1 , the following statements are true:
22 |
√

p (a +ak+1 )
|h11 |
bk bk+1
√ 2 k
≥
=
a) If |h
p1 (bk +bk+1 ) , then g ak , ak+1
21 | q
|h21 |
3p2
∗
∗
and (w1 (k), w2 (k))
=
2
q
bk +bk+1 M −1q

(ak +ak+1 )|h21 |
3p2
3p2
(bk +bk+1 )|h11 |
M 2 −1 ,
M 2 −1 .
√

p (a +ak+1 )
|h11 |
bk bk+1
< √p21 (bkk +bk+1
=
b) If |h
) , then g ak , ak+1
21 | q
|h11 |
3p1
and (w1∗ (k), w2∗ (k))
=
ak +ak+1
M 2 −1
q
q

(bk +bk+1 )|h11 |
3p1
3p1
.
M 2 −1 , (ak +ak+1 )|h21 |
M 2 −1
2) If

|h21 |
|h22 |

a) If

b)

> bk + bk+1 , we have the following results:
|h11 |
|h21 |

|h22 | 
,
bk+1 |h
q21 |
3p2
|h22 | M 2 −1


√
p
√ 2
p1

≥
then
and
q

bk+1 

ak+1
bk+1

g abkk , ak+1
(w1∗ (k), w2∗ (k))

q


3p2
2
M −1 .

ak+1 |h21 |+|h22 |
3p2
bk+1 |h11 |
M 2 −1 ,
√
p
ak+1
|h11 |
√ 2
If
<
|h21 |
p1 bk+1


|h22 |
then
g abkk , abk+1
bk+1 |h21 | ,
k+1
q
bk+1 |h11 ||h22 |
3p1
∗
∗
ak+1 |h21 |+|h22 |
M 2 −1 and (w1 (k), w2 (k))
q

q
bk+1 |h11 |
3p1
3p1
.
M 2 −1 , ak+1 |h21 |+|h22 |
M 2 −1

+
=
=

+
=
=


The proof of Lemma 1 can be found in Appendix-E.
Remark 1: We have the following insights from the above
lemma,
1) It can be observed from Lemma 1 that at least one transmitter should employ the maximum allowable power,
since otherwise we could scale up both transmitted
power without violating the power constraint such that
the minimum Euclidean distance is enlarged.
2) We can see from Lemma 1 that the optimal value of
the objective function together with (w1∗ (k), w2∗ (k)) substantially depend on the relative strength of the channel
coefficients. Inspired by [36], we divide the considered
real Gaussian ZC into three scenarios:
21 |
a) Gaussian ZC with weak cross link, i.e., |h
|h22 | ∈
(0, 1];
|h21 |
b) Gaussian ZC with strong cross link, i.e., |h
∈
22 |
(1, 2M );
c) Gaussian ZC with very strong cross link, i.e.,
|h21 |
|h22 | ∈ [2M, ∞).
Then, a compact closed-form expression for Problem 1
can be established for three complimentary scenarios,
which constitutes main contents of the subsequent subsection.

E. The Optimal NOMA Design with PAM Constellation for
the Gaussian ZC
Now we are ready to give the closed-form optimal solution
of (w1∗ , w2∗ ) to Problem 1 which maximizes the minimum

Euclidean distance over all the Farey intervals for the aforementioned three scenarios.
1) Scenario 1: ZC with Weak Cross Link: For this case, we
21 |
have |h
|h22 | ∈ (0, 1]. Consider the Farey interval set SM−1 =
 b bk+1  |SM −1 |
+ak+2
k
. By Property 1, we have abk+1
<
ak , ak+1
k=1
k+1 +bk+2
ak+1
ak +ak+1
<
,
and
therefore
the
positive
real
axis
can
be
bk+1
bk +bk+1
divided into the |SM−1 | + 1 intervals in an increasing order:
(
 a
 a
|SM −1 | + a|SM −1 |+1
|SM −1 | + a|SM −1 |+1
0,
,
,
b|SM −1 | + b|SM −1 |+1
b|SM −1 | + b|SM −1 |+1
)
a + a

a|SM −1 |−1 + a|SM −1 | 
1
2
,··· ,
,∞ .
b|SM −1 |−1 + b|SM −1 |
b1 + b2
In particular, we have

a1 +a2
b1 +b2

= M and

1
M.

a|SM −1 | +a|SM −1 |+1
b|SM −1 | +b|SM −1 |+1

=

Theorem 1 (Gaussian ZC with weak cross link): Suppose
|h21 |
that |h
∈ (0, 1]. Then, the optimal power scaling factors to
22 |
Problem 1 are explicitly
determined as follows:
√
p
|h11 |
∗
√2 , then, we have (w1
≤
, w2∗ ) =
1) If |h
M qp1
q 21 |

M|h11 |
3p1
3p1
M 2 −1 , |h21 |
M 2 −1 ;
√
M p2
|h11 |
√
≥
we have (w1∗ , w2∗ )
|h21 | q
p1 , then,
q

M|h21 |
3p2
3p2
;
,
2
|h11 |
M 2 −1
 √ M −1

√
p2 (aℓ1 +1 +aℓ1 +2 )
p2 (aℓ1 +aℓ1 +1 )
|h11 |
√
√
,
Let |h21 | ∈
p1 (bℓ1 +1 +bℓ1 +2 )
p1 (bℓ1 +bℓ1 +1 )

2) If

3)

=

for

some ℓ1 = 1, · · · , |SM−1 | − 1. If we let ℓ̃a =
arg mink {(a1 + a2 ), · · · , (aℓ1 + aℓ1 +1 )} and ℓ̃b =
arg mink {(bℓ1 +1 + bℓ1 +2 ), · · · , (b|SM −1 | + b|SM −1 |+1 )},
then we have
q
q

(bℓ̃a +bℓ̃a +1 )|h11 |
3p1
3p1


,
2 −1 (a
2 −1 ,

)|h
|
M
M
+a
21

ℓ̃a
ℓ̃a +1 √


)
p2 (a +a
|h11 |


if |h
≥ √p1 (bℓ̃a +b ℓ̃a +1) ;
21 |
ℓ̃b +1
ℓ̃b
∗
∗
(w1 , w2 ) =  (aℓ̃ +aℓ̃ +1 )|h21 | q 3p q 3p 
2
2
b
b


 (bℓ̃b +bℓ̃b +1 )|h11 | M 2 −1 , M 2 −1 ,

√


p2 (aℓ̃a +aℓ̃a +1 )
|h11 |


√
if |h
).
| <
p (b +b
21

1

ℓ̃b

ℓ̃b +1


The proof of Theorem 1 is provided in Appendix-F.
2) Scenario 2: ZC with Strong Cross Link: In this case,
|h21 |
|h21 |
|h22 | ∈ (1, 2M ). We suppose that L − 1 < |h22 | ≤ L.
Then, the optimal solution to Problem 1 can be obtained
by considering the following
two cases: the Farey interval
′

b′ b
L
set UM−1
= { ak′ , ak+1
∈ SM−1 : b′k + b′k+1 ≥ L} and
′
k
′′
k+1
b′′ b
L
VM−1
= { ak′′ , ak+1
∈ SM−1 : b′′k + b′′k+1 < L}. The whole
′′
k
k+1
discussions on them can be summarized as the following
theorem.
Theorem 2 (ZC with Strong Cross Link): The optimal
solution to Problem 1 in this case can be obtained by finding
the maximum value of the objective functions in the following
two cases which are explicitly attained as follows:
′

b′ b
L
∈ UM−1
. Then, the following three
1) Let ak′ , ak+1
′
k
k+1
statements are true.
√
′
′
a) If
then

|h11 |
|h21 |

we

≤

p2 (a|U L
+a|U L
)
|
|+1
M −1
M −1
√
),
+b′ L
p1 (b′ L
|
|U
M −1

have

(wu∗ 1 , wu∗ 2 )

|+1
|U
M −1

=

8

q

q

b)

c)



(b′ℓ̃ +b′ℓ̃ +1 )|h11 |
3p1
3p1
c
c
,
2
)|h21 |
M −1 (a′ +a′
M 2 −1 , where ℓ̃c =
ℓ̃c +1
ℓ̃c
arg mink {(a′1 + a′2 ), . . . , (a′|U L | + a′|U L |+1 )}.
M −1
M −1
√
M p2
11 |
∗
∗
√
≥
,
then
we
have
(w
If |h
u1 , wu2 ) =
p1 q
 |h21 | q

M|h21 |
3p2
3p2
|h11 |
M 2 −1 ,
M 2 −1 .
|h11 |
∈
Suppose
that
21 | 
 √p (a′ +a′ ) √p (a′ +a|h
′
2
2
ℓ2 +1 )
ℓ2 +1
ℓ2 +2
ℓ2
√
√
for
p1 (b′ℓ +1 +b′ℓ +2 ) , p1 (b′ℓ +b′ℓ +1 )
2
2
2
2
′
L
some ℓ2
=
1, · · · , |UM−1 |. If we let

ℓ̃d = arg mink {(a1 + a2 ), . . . , (aℓ2 + aℓ2 +1 )} and
ℓ̃e = arg mink {(bℓ2 +1 + bℓ2 +2 ), . . . , (b|UM
L
+
−1 |
L
b|UM
)},
then
we
have
−1 |+1
q
(b′ℓ̃ +b′ℓ̃ +1 )|h11 | q 3p 
3p1
1
d
d


,
2

M −1 (a′ +a′
)|h21 |
M 2 −1 ,

ℓ̃d
ℓ̃d +1

√
′
′


 if |h11 | ≥ √p2 (aℓ̃d +aℓ̃d +1 ) ;

′
′
)
|h21 |
p1 (b +b
ℓ̃e +1
ℓ̃e
q
(wu∗ 1 , wu∗ 2 )= (a′ +a′ )|h21 | q


3p
3p2
ℓ̃
+1
ℓ̃
2
e
e

,
′
′

2
)|h11 |
(b +b
M −1
M 2 −1

ℓ̃e +1 √
ℓ̃e



p (a′ +a′
)

 if |h11 | < √ 2 ℓ̃′ d ′ℓ̃d +1 .
)
|h21 |
p1 (b +b
ℓ̃e

b′′
k−1
a′′
k−1

b′′
k
a′′
k

ℓ̃e +1



L
,
∈ VM−1
. Then, the following two
2) Let
statements are true.√
p
|h11 |
a) If |h
≤ √p21 , then we have (wv∗1 , wv∗2 ) =
q
q 22 |
′′
b|V
|+1 |h11 ||h22 |
1
1
( M3p2 −1
, a′′ M −1 |h21 |+|h22 | M3p2 −1
).
|VM −1 |+1
√
′′
′′

a
a
|h | p1
22 |
√
b) If |h11
∈ b′′ℓ3 +1 + b′′ |h22|h|21 | , b′′ℓ3 + b′′|h|h
21 | p2
21 |
ℓ3 +1
ℓ3 +1
ℓ3
ℓ3
for some ℓ3 = 1, · · · , |VM−1 |, then we have
 q
q

b′′
ℓ3 |h11 |
3p1
3p1


,
2 −1 a′′ |h |+|h |
2 −1 ,

M
M
21
22

ℓ

√3
′′



 if |h11 | ≥ √p2 a′′ℓ3 + ′′|h22 | ;
|h21 |
p1 bℓ
bℓ |h21 |
∗
∗
3q
(wv1 , wv2 ) =
q3

a′′
ℓ3 +1 |h21 |+|h22 |
3p2

2

, M3p2 −1
′′
2

bℓ +1 |h11 |
M2 −1

2
3

√

p2 aℓ3 ′′
|h11 |
|h22 | 

 if |h
√
| <
p
b′′ + b′′ |h | .
21

1

ℓ3

ℓ3

21


The proof of Theorem 2 is very similar to that of Theorem 1
and the following Theorem 3 and thus, is omitted due to space
limitation.
3) Scenario 3: ZC with Very Strong Cross Link: In this
21 |
case, |h
|h22 | ∈ [2M, ∞). Likewise, we consider the Farey
 b bk+1  |SM −1 |
21 |
k
. Note that |h
interval set SM−1 =
ak , ak+1
|h22 | ≥
k=1
2M > bk + bk+1 for k = 1, . . . , |SM−1 |. Then, by Property 1
ak+1
ak +ak+1
|h22 |
and Proposition 4, we have bk+1
+ bk+1
|h21 | ≤ bk +bk+1 <
< abkk + bk|h|h2221| | . As a result, the positive real axis can be
divided into the following intervals in increasing order:
ak
bk

n

0,

a|SM −1 |+1

|h22 |
+
,
b|SM −1 |+1
b|SM −1 |+1 |h21 |
|h22 | a1
|h22 | o
a2
+
+
,
,
...,
b2
b2 |h21 | b1
b1 |h21 |

where ab11 + b1|h|h2221| | = ∞.
Theorem 3 (Gaussian ZC with very strong cross link): Let
|h21 |
|h22 | ∈ [2M, ∞). Then, the optimal solution to Problem 1 is
given below:

√
p

11 |
1) If |h
≤ √p21 , then we have (w1∗ , w2∗ ) =
q
q |h22 |
|h11 |
3p1
1
( M3p
2 −1 , |h |
M 2 −1 ).
22
√

a +1
a
|h | p1
22 |
22 |
√
, bℓℓ4 + bℓ|h|h
∈ bℓℓ4+1
+ bℓ |h
for some
2) If |h11
|h
|
|
21 | p2
21
21
+1
4
4
4
4
ℓ4 = 1, . . . , |SM−1 |, then we have
q
 q

bℓ4 |h11 |
3p1
3p1

,
2

M −1 a√
M 2 −1 ,
ℓ4 |h21 |+|h22 |




 if |h11 | ≥ √p2 aℓ4 + |h22 | ;
|h21 |
p1 q
bℓ4
bℓ4 |h
|
∗
∗
21
q
(w1 , w2 ) =

a
|h21 |+|h22 |
3p2
3p2

 ℓ4 +1
,
2

bℓ4 +1 |h11
|
M −1
M 2 −1 ,

√

 if |h11 | < √p2 aℓ4 + |h22 | .
|h21 |
p1 bℓ
bℓ |h21 |
4

4



The proof of Theorem 3 is provided in Appendix-G.
IV. S IMULATION R ESULTS AND D ISCUSSIONS
In this section, computer simulations are carried out to
demonstrate the effectiveness of our proposed NOMA design under different channel configurations. More precisely,
we compare our proposed NOMA design with CR based
NOMA [30], time-division multiple access (TDMA) and
frequency-division multiple access (FDMA) approaches. Without loss of generality, we set P1 = P2 = 1. For simplicity,
in the simulations, we assume that the same square-QAM
constellation is adopted by both users, i.e., M = M ′ and
according to (6), we have p1 = p′1 = P1 /2 = 1/2 and
p2 = p′2 = P2 /2 = 1/2.
A. The Resulting Optimal Sum Constellation at Receiver D1
For Several Deterministic Channels
Here, we consider several deterministic channels corresponding to the three scenarios of weak cross link, strong
cross link and very strong cross link. We assume that 16-QAM
constellations are used by both users with M = M ′ = 4. We
discuss the results of these cases one by one as follows:
1) Weak Cross Link: In this case, we assume that h11 = 1,
h21 = 1/2 and h22 = 1. Based on the derived expressions
for the optimal solution provided in the previous section,
we can readily obtain that (w1∗ , w2∗ ) = (0.4472, 0.2236).
The corresponding received constellation at D1 is plotted in
Fig. 2(a). It can be observed that the sum-constellation at D1 is
a regular 256-QAM generated by the superposition of two 16QAM. Hereafter, we call the signal constellation with smaller
minimum Euclidean distance as the satellite constellation. By
observing that w1∗ h11 = 0.4472 and w2∗ h21 = 0.1118 (i.e.,
w2∗ h21 = 41 w1∗ h11 ), we can deduce that the constellation used
by S2 forms the satellite constellation of the sum constellation
at D1 .
2) Strong Cross Link: We investigate two channel realizations for this scenario. For the first realization, we let
h11 = 1, h21 = 3/2 and h22 = 1. We then have (w1∗ , w2∗ ) =
(0.1677, 0.4472) and the resulting sum-constellation at D1 is
also regular, as illustrated in Fig. 2(b). Since w1∗ h11 = 0.1677
and w2∗ h21 = 0.6708 (i.e., w2∗ h21 = 4w1∗ h11 ) in this case,
the constellation used by S1 forms the satellite constellation
at D1 . For the second realization, we set h11 = 1, h21 = 3
and h22 = 1, leading to (w1∗ , w2∗ ) = (0.3354, 0.4472). The

9

2.5
1.5
2

1

1.5

1
0.5
Quadrature

Quadrature

0.5

0

0

−0.5
−0.5
−1

−1

−1.5

−2
−1.5
−2.5
−1.5

−1

−0.5

0
In−Phase

0.5

1

1.5

−2.5

−2

−1.5

−1

−0.5

(a)

0
0.5
In−Phase

1

1.5

2

2.5

(b)

5
6

4

3

4

2
2
Quadrature

Quadrature

1

0

0

−1
−2
−2
−4
−3

−4

−6

−5
−5

−4

−3

−2

−1

0
In−Phase

1

2

3

4

5

−6

−4

−2

0
In−Phase

2

4

6

(c)
(d)
Fig. 2. Optimal resulting constellations at receiver D1 for various cases: (a) Weak cross link with h11 = 1, h21 = 1/2 and h22 = 1. (b) Strong cross link
with h11 = 1, h21 = 3/2 and h22 = 1. (c) Strong cross link with h11 = 1, h21 = 3 and h22 = 1 and (d) Very strong cross link with h11 = 1, h21 = 5
and h22 = 1/2.

resulting constellation plotted in Fig. 2(c) is also uniform as
in the previous two scenarios. We have w1∗ h11 = 0.3354
and w2∗ h21 = 1.3416 (i.e., w2∗ h21 = 4w1∗ h11 ). Thus, the
constellation used by the transmitter S1 forms the satellite
constellation at D1 .
3) Very Strong Cross Link: In this case, we suppose that
h11 = 1, h21 = 5 and h22 = 1/2, generating (w1∗ , w2∗ ) =
(0.2236, 0.4472). The obtained constellation at D1 is shown in
Fig. 2(d). In this case, we have w1∗ h11 = 0.2236 and w2∗ h21 =
2.236 (i.e., w2∗ h21 = 10w1∗ h11 ). However, it can be observed
that w2∗ h22 = 0.2236, i.e., w1∗ h11 = w2∗ h22 .
B. Average Error Performance Comparison in Rayleigh Fading Channels
We now compare the average BER of the proposed NOMA
scheme with that of three existing methods, including TDMA,
FDMA, and NOMA with CR [30] methods, over Rayleigh
2
2
fading channels with h11 ∼ CN (0, δ11
), h21 ∼ CN (0, δ21
)

2
and h22 ∼ CN (0, δ22
). Recall that we use error performance
(i.e., BER) as the design criterion for the NOMA in ZCs
with finite-alphabet inputs using fixed transmission rate (i.e.,
fixed constellation size). However, we are unable to compare
the error performance of the considered system using finitealphabet inputs with that of Gaussian inputs. This is because
for Gaussian input, it is intractable to evaluate the BER for
uncoded system since its input signal is continuous. Moreover,
the BER for coded system with Gaussian input is hard to
simulate due to the huge storage capacity requirement for the
large codebook and the high computational complexity [59,
Ch. 9].
For TDMA, we assume that both users transmit alternatively
by using half of the total time slots and thus no interference
occurs at the destination side. More importantly, the individual
instantaneous power constraints on both users S1 and S2
remain unchanged. On the other hand, for FDMA, each
user uses only half of the available bandwidth. Due to the

10

0

0

10

10

NOMA, proposed
FDMA
TDMA
NOMA, CR

NOMA, proposed
FDMA
TDMA
NOMA, CR

−1

−1

10

BER

BER

10

−2

10

−2

10

−3

10

−3

10
0

10

20

30

40

50

60

0

10

20

30

40

50

60

70

SNR, dB

70

SNR, dB

(a)
(a)
0

10

NOMA, proposed
FDMA
TDMA
NOMA, CR

0

10

NOMA, proposed
FDMA
TDMA
NOMA, CR
−1

10
−1

BER

BER

10

−2

10

−2

10

−3

10

0

10

20

30

40
SNR, dB

50

60

70

80

−3

10

0

10

20

30

40
SNR, dB

50

60

70

80

(b)
Fig. 3. Comparison of the proposed NOMA, CR based NOMA, TDMA and
2
2
2
FDMA with (δ11 , δ21 , δ22 ) = (1, 1, 1): (a) 16-QAM is used for proposed
NOMA while 16-PSK is used for CR based NOMA. (b) 64-QAM is used for
proposed NOMA while 64-PSK is used for CR based NOMA.

orthogonality between different frequency band, there is also
no interference occurring at the destination side. Note that,
in FDMA, the bandwidth occupied by each user is halved
and the noise arises at the receiver is assumed to be white
Gaussian. Therefore, the variance of the noise is also halved. In
addition, for the CR based NOMA as proposed in [30], we let
each user transmit at the maximum allowable power by using
j2πℓ+jπ
N −1
−1
constellations {exp( j2πk
)}N
ℓ=0 for
N )}k=0 and {exp(
N
user S1 and S2 , respectively.
In Fig. 3, we consider that the variances of all channels are
2
2
2
the same, i.e., (δ11
, δ21
, δ22
) = (1, 1, 1), and the average BER
over both receivers of all the methods are plotted against the
SNR ρ = 2σ1 2 . In Fig. 3(a), 16-QAM is used for the proposed
NOMA scheme while 16-PSK is employed by the NOMA with
CR. Since only half of the total time slots or total bandwidths
are available for each transmitter, to maintain the same data
rate for each user in each block compared with NOMA
methods, we should increase the constellation size by using
256-QAM in both TDMA and FDMA. It can be observed that,
our proposed method has significant BER gain over TDMA
and FDMA methods, which confirms the effectiveness of the
NOMA scheme. From the simulation results, we also find
that FDMA has a smaller BER compared with TDMA. This
is because the variance of the effective noise is smaller than

(b)
Fig. 4. Comparison of the proposed NOMA, CR based NOMA, TDMA and
2 , δ 2 , δ 2 ) = (1, 1/4, 1): (a) 16-QAM is used for proposed
FDMA with (δ11
21 22
NOMA while 16-PSK is used for CR based NOMA. (b) 64-QAM is used for
proposed NOMA while 64-PSK is used for CR based NOMA.

that of TDMA. As the rotation based method uses the PSK
constellation, which is not spectrally efficient, it has the worst
BER performance. Then, in Fig. 3(b), the average BER of
all the cases are plotted against the SNR where 64-QAM is
used by each user for the proposed NOMA and 64-PSK is
used by CR based NOMA while 4096-QAM constellations
are used by TDMA and FDMA methods. We also simulate
another two cases with unequal channel variances. Specifically,
the average BER of both receivers for all the considered
methods is plotted in Fig. 4, wherein the variances of three
2
2
2
channels are set as (δ11
, δ21
, δ22
) = (1, 1/4, 1). The case with
2
2
2
the variances of three channels (δ11
, δ21
, δ22
) = (1, 4, 1) is
provided in Fig. 5. Similar observations can be seen as the
previous case with equal channel variances, which further
verifies the effectiveness of our proposed NOMA design.
V. C ONCLUSIONS
In this paper, we developed a novel and practical design
framework for the non-orthogonal multiple access (NOMA)
in a classical two-transmitter two-receiver Z-channel with
widely-used quadrature amplitude modulation (QAM) and
max-min user fairness. Specifically, we formulated a max-min
optimization problem to jointly optimize the scaling factors at
both transmitters to maximize the smaller minimum Euclidean
distance among the two resulting signal constellations at both

11

0

10

NOMA, proposed
FDMA
TDMA
NOMA, CR

−1

BER

10

−2

10

−3

10

0

10

20

30
SNR, dB

40

50

60

(a)

By a similar argument on S2K,3 and S2K,4 , it follows that
min(m,n)∈S2K,1 d1 (m, n) = min(m,n)∈Z2K,1 \{(0,0)} d1 (m, n).
Then, we further divide S2K,1 into F2K and S2K,1 \ F2K . In
what follows, we will show that min (m,n)∈F2K d1 (m, n) =
min (m,n)∈S2K,1 d1 (m, n). This can be proved by contradiction. Suppose that min (m,n)∈S21 \F2K d1 (m, n) =
min (m,n)∈S2K,1 d1 (m, n), where the minimum is achieved
by (m∗ , n∗ ) ∈ S2K,1 \ F2K such that hm∗ , n∗ i = ℓ > 1 by the
∗
∗
definition of Farey sequences. Then, we can find ( mℓ , nℓ ) ∈
∗
∗
F2K such that d1 ( mℓ , nℓ ) = 1ℓ d1 (m∗ , n∗ ) < d1 (m∗ , n∗ ),
which contradicts the assumption. This completes the proof. 
B. Proof of Proposition 2

0

10

Recall that d1 (m, n) = |h11 |w1 n − |h21 |w2 m . Therefore,

n1 n2
21 |w2
for |h
|h11 |w1 ∈ m1 , m2 , we have: d1 (m1 , n1 ) − d1 (m2 , n2 ) =

n1 +n2
21 |w2
(m1 + m2 )|h11 |w1 |h
|h11 |w1 − m1 +m2 . The results presented
in the proposition can be readily obtained, and we complete
the proof.


NOMA, proposed
FDMA
TDMA
NOMA, CR

−1

BER

10

−2

10

C. Proof of of Proposition 3
−3

10

0

10

20

30

40

50

60

70

SNR, dB

(b)
Fig. 5. Comparison of the proposed NOMA, CR based NOMA, TDMA and
2
2
2
FDMA with (δ11 , δ21 , δ22 ) = (1, 4, 1): (a) 16-QAM is used for proposed
NOMA while 16-PSK is used for CR based NOMA. (b) 64-QAM is used for
proposed NOMA while 64-PSK is used for CR based NOMA.

receivers, subject to the individual average power constraint
at each transmitter. The formulated mixed continuous-discrete
problem was successfully resolved in compact closed-form by
strategically applying the Farey sequences and their unique
properties. Simulation results verified the correctness of our
analytical derivations and showed that the proposed NOMA
design significantly outperforms the existing orthogonal multiple access and NOMA schemes, especially at high signal-tonoise ratio. Furthermore, the performance gap of the proposed
scheme over its existing counterparts can be further enlarged
when the size of constellations used at the transmitter side
becomes larger.
A PPENDIX
A. Proof of Proposition 1
We first divide the feasible region Z2K \ {(0, 0)} into four
subsets given by S2K,1 = {(m, n) : m, n ∈ NK , (m, n) 6=
(0, 0)}, S2K,2 = {(−m, n) : m, n ∈ NK , (m, n) 6=
(0, 0)}, S2K,3 = {(m, −n) : m, n ∈ NK , (m, n) 6=
(0, 0)}, S2K,4 = {(−m, −n) : m, n ∈ NK , (m, n) 6=
(0, 0)}. Then clearly we have Z2K \ {(0, 0)} = ∪4k=1 S2K,k .
For (−m, n) ∈ S2K,2 with m, n ≥ 0, we can always
find (m, n) ∈ S2K,1 such that d1 (m, n) = |h11 |w1 n −
|h21 |w2 m ≤ |h11 |w1 n − |h21 |w2 (−m) = d1 (−m, n).
Hence, min(m,n)∈S2K,1 d1 (m, n) ≤ min(m,n)∈S2K,2 d1 (m, n).

|h21 |w2
n2
∈ (m
, n2 +n3 ). By
We first consider the case |h
11 |w1
2 m2 +m3
using Proposition 2, we have d1 (m2 , n2 ) < d1 (m3 , n3 ).
|h21 |w2
Then, with the help of Property 1 and 3, we have |h
∈
11 |w1
n2 +n3
n2 n3
n2
n2 +n4
n2
( m2 , m2 +m3 ) ⊂ ( m2 , m3 ) ⊆ ( m2 , m2 +m4 ) and using Property 2 again, we have d1 (m2 , n2 ) ≤ d1 (m4 , n4 ). On the other
n2
n2 +n3
21 |w2
hand, by Property 1 and 3, we have |h
|h11 |w1 ∈ ( m2 , m2 +m3 ) ⊂
n3
n2 n3
3
,
) ⊆ ( mn11 +n
(m
+m3 , m3 ) and using Property 2 again, we
2 m3
have d1 (m3 , n3 ) ≤ d1 (m1 , n1 ).
n4
n1
and m
are randomly picked entry in SK , or
As m
1
4
equivalently (m1 , n1 ) and (m4 , n4 ) are randomly picked in
21 |w2
F2M−1 \ {(m2 , n2 ), (m3 , n3 )}, this proves that for |h
|h11 |w1 ∈
n2
+n3
( m2 , mn22 +m
), min(m,n)∈F2M −1 d1 (m, n) = d1 (m2 , n2 ) =
3
|h21 |w2 m2 − |h11 |w1 n2 . The other case can also be proved
by a similar argument and hence omitted. We completes the
proof.


D. Proof of Proposition 4
ak +ak+1
bk +bk+1 −
|h21 |
|h22 | ≥ bk +

First of all, we calculate the following difference

|h21 |−|h22 |(bk +bk+1 )
(bk +bk+1 )bk+1 |h21 | . Hence, if
ak+1
ak +ak+1
|h22 |
bk+1 , then bk+1 + bk+1
|h21 | ≤ bk +bk+1 . Similarly, we have
ak +ak+1
|h21 |−|h22 |(bk +bk+1 )
|h22 |
ak
bk − bk |h21 | − bk +bk+1 =
bk (bk +bk+1 )|h21 | . As a result,
ak +ak+1
|h21 |
if |h22 | ≥ bk + bk+1 , then bk +bk+1 ≤ abkk − bk|h|h2221| | . This
ak+1
bk+1

|h22 |
− bk+1
|h21 | =

completes the proof.



E. Proof of Lemma 1

b
=
According to proposition 3 and notice that abkk , ak+1
k+1

bk +bk+1 bk+1 
bk bk +bk+1
ak , ak +ak+1 ∪ ak +ak+1 , ak+1 , problem in (11) can be
further divided into the following two sub-problems and the
overall solution is the maximum value of the two problems:

12

Problem 2 (Sub-problem 1): We aim to solve the following
optimization problem:
b b

k
k+1
g1
= max min {|h11 |w1 bk+1
,
w1 ,w2
ak ak+1
− |h21 |w2 ak+1 , |h22 |w2 } (12a)
|h21 |w2
bk+1
bk + bk+1
≤
≤
,
(12b)
s.t.
ak + ak+1
|h11 |w1
ak+1
r
r
3p1
3p2
0 < w1 ≤
, 0 < w2 ≤
.
(12c)
M2 − 1
M2 − 1

Problem 3 (Sub-problem 2): The optimization problem is
stated as follows:
b b

k
k+1
= max min{|h21 |w2 ak
g2
,
w1 ,w2
ak ak+1
− |h11 |w1 bk , |h22 |w2 }
(13a)
|h21 |w2
bk + bk+1
bk
≤
<
,
(13b)
s.t.
ak
|h11 |w1
ak + ak+1
r
r
3p1
3p2
0 < w1 ≤
,
0
<
w
≤
.
(13c)
2
M2 − 1
M2 − 1

We first consider Sub-problem 1 and Problem 2, which can
be divided into the following two case:
1) Case 1: Receiver 1 has smaller Euclidean distance.
b b

k
k+1
g11
= max bk+1 |h11 |w1 − ak+1 |h21 |w2
,
w1 ,w2
ak ak+1
s.t. bk+1 |h11 |w1 − ak+1 |h21 |w2 ≤ |h22 |w2 ,
(14a)
|h21 |w2
bk+1
bk + bk+1
≤
≤
,
(14b)
ak + ak+1
|h |w
ak+1
r
r 11 1
3p1
3p2
, 0 < w2 ≤
. (14c)
0 < w1 ≤
M2 − 1
M2 − 1
2) Case 2: Receiver 2 has smaller Euclidean distance.
b b

k
k+1
g12
= max |h22 |w2
,
w1 ,w2
ak ak+1
s.t. |h22 |w2 < bk+1 |h11 |w1 − ak+1 |h21 |w2 ,
(15a)
|h21 |w2
bk+1
bk + bk+1
≤
<
,
(15b)
ak + ak+1
|h |w
ak+1
r 11 1
r
3p1
3p2
0 < w1 ≤
, 0 < w2 ≤
. (15c)
2
M −1
M2 − 1
For Case-1 of Sub-problem 1, constraint (14a) is equiv|h11 |
alent to ak+1bk+1
|h21 |+|h22 | w1 ≤ w2 and constraint (14b) means
bk+1 |h11 |
(bk +bk+1 )|h11 |
(ak +ak+1 )|h21 | w1 ≤ w2 ≤ ak+1 |h21 | w1 . Also we notice that if
(bk +bk+1 )|h11 |
bk+1 |h11 |
|h21 |
|h22 | ≥ bk + bk+1 then (ak +ak+1 )|h21 | w1 ≤ ak+1 |h21 |+|h22 | w1 .

Hence, the optimization problem can be further divided into
two cases:
21 |
1) If |h
|h22 | ≤ bk + bk+1 , then
b b

k
k+1
g11
= max bk+1 |h11 |w1 − ak+1 |h21 |w2
,
w1 ,w2
ak ak+1
bk+1 |h11 |
(bk + bk+1 )|h11 |
w1 ≤ w2 ≤
,
(16a)
s.t.
(ak + ak+1 )|h21 |
ak+1 |h21 |
r
r
3p1
3p2
0 < w1 ≤
, 0 < w2 ≤
. (16b)
2
M −1
M2 − 1

(bk +bk+1 )|h11 |
(ak +ak+1 )|h21 | w1 , then the objective
11 |w1
. In this case, (16b) is
function is a|h
k +ak+1
q
(a +ak+1 )|h21 |
3p2
equivalent to w1 ≤ (bkk +bk+1
)|h11 |
M 2 −1 and
q

3p1
bk bk+1
w1 ≤
=
M 2 −1 and hence, we have g11 ak , ak+1
q
q
o
n
(ak +ak+1 )|h21 |
|h11 |
3p2
3p1
11 |
,
min ak|h
2
2
+ak+1 (bk +bk+1 )|h11 |
M −1 ak +ak+1
M −1 .

We let w2 =

As a consequence, we have

g11

b

k

,

bk+1 
ak+1
q

ak

|h21 |
3p2


bk +bk+1
M 2 −1 , where q

q



(ak +ak+1 )|h21 |

3p2
3p2

(w
,
w
)
=
,
1
2
2

(bk +bk+1 )|h11 |
M −1
M 2 −1 ,

√


p (a +ak+1 )
11 |

√ 2 k
if |h
|h21 | ≥
p1 (bk +bk+1 ) ;
q
=
|h
|
3p
11
1



ak +ak+1
M 2 −1 , where

q
q



(bk +bk+1 )|h11 |
3p1
3p1

(w1 , w2 ) =
,
2

M −1 (ak +ak+1 )|h21 |
M 2 −1 ,

√


p2 (ak +ak+1 )
11 |
√
if |h
|h21 | <
p1 (bk +bk+1 ) ,

2) If

|h21 |
|h22 |

> bk + bk+1 , then

bk+1 
= max bk+1 |h11 |w1 − ak+1 |h21 |w2
w1 ,w2
ak ak+1
bk+1 |h11 |
bk+1 |h11 |
w1 ≤ w2 ≤
w1 , (18a)
s.t.
ak+1 |h21 | + |h22 |
ak+1 |h21 |
r
r
3p1
3p2
0 < w1 ≤
, 0 < w2 ≤
. (18b)
2
M −1
M2 − 1

g11

b

k

,

|h11 |
|h11 |
< abk+1
w1
We first notice that ak+1bk+1
|h21 |+|h22 | w1
k+1 |h21 |
and hence the problem is always feasible. By letting
|h11 |
w2 = ak+1bk+1
|h21 |+|h22 | w1 , the objective function can be
bk+1 |h11 ||h22 |
ak+1 |h21 |+|h22 | w1 .

In this case, the constraints
q
3p1
in (18b) are equivalent to w1 ≤
M 2 −1 and w1 ≤
q
ak+1 |h21 |+|h22 |
3p2
bk+1 |h11 |
M 2 −1 . Therefore, the solution is
written by

g11

b

k

,

bk+1 
ak+1

ak
 b |h ||h | q
3p1
11
22
k+1

 ak+1 |h21 |+|h22 | M 2 −1 , where

q
q



bk+1 |h11 |
3p1
3p1
(w , w ) =

,
2
1
2

a
|h
|+|h
|
M 2 −1 ,
M1 −1
21
22
k+1

√


p
ak+1
|h22 | 
|h11 |

≤ √p21 bk+1
+ bk+1
if |h
|
|h21 | ;
21
q
=
3p
2

|h |
, where


 22 M22 −1
q
q



3p2

2
(w1 , w2 ) = ak+1b |h21|h|+|h| 22 | M3p2 −1
,
2 −1 ,

M
11
k+1

√


p
|h22 | 
|h11 |
k+1
> √p21 abk+1
+ bk+1
if |h
|h21 | .
21 |

For Case-2 of Sub-problem 1, constraint (15a) is equiva|h11 |
lent to w2 < ak+1bk+1
|h21 |+|h22 | w1 and constraint (15b) implies
(bk +bk+1 )|h11 |
(ak +ak+1 )|h21 | w1

≤ w2 ≤

bk+1 |h11 |
ak+1 |h21 | w1 .

By noticing that

13

bk+1 |h11 |
ak+1 |h21 |+|h22 |

b

bk+1 |h11 |
ak+1 |h21 | ,

<

the problem is equivalent to



bk+1
= max |h22 |w2
w1 ,w2
ak ak+1
bk+1 |h11 |
(bk + bk+1 )|h11 |
w1 ≤ w2 ≤
w1 ,
s.t.
(ak + ak+1 )|h21 |
ak+1 |h21 | + |h22 |
(20a)
r
r
3p1
3p2
0 < w1 ≤
, 0 < w2 ≤
.
(20b)
M2 − 1
M2 − 1

g12

k

,

21 |
Constraint (20a) is feasible if |h
|h22 | ≥ bk + bk+1 . In this case,
the solution is
b b

k
k+1
g12
,
ak ak+1
q

bk+1 |h11 ||h22 |
3p1


ak+1 |h21 |+|h22q
|
M 2 −1 , where

q



bk+1 |h11 |

3p1
3p1

(w
,
w
)
=
2 −1 , a
1
2

M
|h
|+|h
|
M 2 −1 ,

21
22
k+1
√


p2 ak+1
|h22 | 
11 |

√
if |h
|h21 | ≤
p1 bk+1 + bk+1 |h21 | ;
q
=
2

, where
|h22 | M3p

2

2 −1

q
q



|h21 |+|h22 |
3p2
3p2

,
,
(w1 , w2 ) = ak+1bk+1

2

|h
|
M2 −1
M22 −1
11

√

p2 ak+1

|h11 |
|h22 | 
if |h21 | > √p1 bk+1 + bk+1 |h21 | .

The solution of Sub-problem 2 can be attained in a similar
fashion as Sub-problem 1, and hence is omitted.
Then, for

|h21 |w2
bk bk+1
∈
,
,
the
solution
to
the subinterval division |h
ak ak+1
11 |w1
both Sub-problems can be summarized as follows:
+bk+1
|h21 |w2
≤ |h
:
≤ abk+1
1) Scenario 1, abkk +a
11 |w1
k+1
k+1
|h21 |
|h22 |

≤ bk + bk+1

k bk+1
,
g11
ak ak+1
q

|h21 |
3p2


bk +bk+1
M 2 −1 , where

q
q



(ak +ak+1 )|h21 |

3p2
3p2

(w
,
w
)
=
2 −1 ,
1
2

(b
+b
)|h
|
M
M 2 −1 ,
11
k

√ k+1


p (a +ak+1 )
|h11 |

≥ √p21 (bkk +bk+1
if |h
);
21 |
q
=
|h11 |
3p1



ak +ak+1
M 2 −1 , where

q
q



(bk +bk+1 )|h11 |
3p1
3p1


(w1 , w2 ) =
,
2

M −1 (ak +ak+1 )|h21 |
M 2 −1 ,

√


p (a +ak+1 )
|h11 |
< √p21 (bkk +bk+1
if |h
),
21 |

a) If

b

|h21 |
|h22 |

> bk + bk+1

k bk+1
,
g11
ak ak+1
q

bk+1 |h11 ||h22 |
3p1


ak+1 |h21 |+|h22 |
M 2 −1 , where (w1 , w2 ) =

q
q



bk+1 |h11 |

3p1
3p1

2 −1 , a

M
|h
|+|h
|
M 2 −1 ,
21
22

k+1 √


p
|h11 |
|h22 | 
k+1

if |h
≤ √p21 abk+1
+ bk+1
|h21 | ;
21 |
q
=
3p2


, where (w1 , w2 ) =
|h |

 22 M 2 −1 q
q



3p2

2
 ak+1b |h21|h|+|h| 22 | M3p2 −1
,

M 2 −1 ,
11
k+1

√

p

|h22 | 
|h11 |
k+1
> √p21 abk+1
+ bk+1
if |h
|h21 | .
21 |

b) If

b

+bk+1
|h21 |w2
2) Scenario 2, abkk ≤ |h
:
≤ abkk +a
11 |w1
k+1
a) Case 1: Receiver 1 has smaller Euclidean distance.

i) If

|h21 |
|h22 |

g21

b

≤ bk + bk+1 , the solution is

k

,

bk+1 
ak+1
q

ak

|h21 |
3p2

, where (w1 , w2 ) =

 bk +bk+1 M 2 −1q
q



3p2
3p2

k +ak+1 )|h21 |
 (a
,
2

(bk +bk+1 )|h11 |
M −1
M 2 −1 ,

√


p (a +ak+1 )
|h11 |

≥ √p21 (bkk +bk+1
if |h
|
);
21
q
=
|h11 |
3p1


(w1 , w2 ) =

aq
M 2 −1 , whereq
k +ak+1




(bk +bk+1 )|h11 |
3p1
3p1


2 −1 , (a +a
2 −1 ,

M
)|h
|
M
21
k
k+1

√


p (a +ak+1 )
|h11 |
√ 2 k
<
if |h
|
p1 (bk +bk+1 ) .
21

ii) If

|h21 |
|h22 |

b

≥ bk + bk+1 ,

bk+1 
ak ak+1
q

3p2

|h
|
22

M 2 −1 , where (w1 , w2 ) =

q
q




3p2
3p2
 |h21 |ak −|h22 |
,

2
2 −1 ,
bk |h11 |
M
−1
M

√



p
|h11 |

≥ √p21 abkk − |h|h2122|b|k ;
if |h
21 |q
=
bk |h11 ||h22 |
3p1


(w1 , w2 ) =

|h
M 2 −1 , where
21 |ak −|h22 |

q
q



bk |h11 |
3p1
3p1


2 −1 , |h |a −|h |
2 −1 ,

M
M
21
22
k

√


p

|h11 |
< √p21 abkk − |h|h2122|b|k .
if |h
21 |

g21

k

,

b) Case 2: Receiver 2 has smaller Euclidean distance.
21 |
i) If |h
|h22 | ≥ bk + bk+1 , the solution is
g22

b

k

,

bk+1 
ak+1

ak
q

(bk +bk+1 )|h11 ||h22 |
3p1

 (ak +ak+1 )|h21 |
M 2 −1 , where (w1 , w2 )

q
q



(bk +bk+1 )|h11 |
3p1
3p1
=

,
2

M −1 (ak +ak+1 )|h21 |
M 2 −1 ,

√


 if |h11 | ≤ √p2 (ak +ak+1 ) ;
|h21 |
p1 (bk +bk+1 )
q
=
3p2


,
where
(w1 , w2 ) =
|h
|
22

M 2 −1 q

q



(ak +ak+1 )|h21 |
3p2
3p2


,
2

(bk +bk+1 )|h11 |
M −1
M 2 −1 ,

√


p (a +ak+1 )
|h11 |
if |h
> √p21 (bkk +bk+1
).
21 |

bk bk+1
ak , ak+1



bk bk+1
ak , ak+1 .
|h22 |
ak
21 |
By Proposition 4, for |h
k+1 , we have bk − |h21 |bk ≥
|h22 | ≥ bk +b
√

p
ak +ak+1
|h11 |
|h22 |
2
ak
√
bk +bk+1 . Also, for |h21q
| ≤
p1 bk − bk |h21 | , we have
q
bk |h11 |
|h21 |
3p2
3p1
M 2 −1 ≥ |h21 |ak −|h22 |
M 2 −1 . In addition, for |h22 | ≥ bk +
(bk +bk+1 )
bk+1 , we have (ak +a
≥ |h21 |abkk−|h22 | . Hence, we have
k+1 )|h21 |


b
b
21 |
≥ g21 abkk , ak+1
. We notice that for |h
g22 abkk , ak+1
|h22 | ≤
k+1
k+1


|h21 |
bk + bk+1 , g11 abkk , abk+1
= g21 abkk , abk+1
. Then for |h
≥
k+1
k+1

22 |
bk bk+1
bk bk+1
bk + bk+1 , we combine g11 ak , ak+1 and g22 ak , ak+1 . By
√
p
ak+1
21 |
√ 2
Proposition 4, for |h
≥
b
+
b
,
we
have
k
k+1
|h22 |
p1 bk+1 +
√

p (a +ak+1 )
|h22 |
|h21 |
√ 2 k
k+1 , we atbk+1 |h21 | ≤
p1 (bkq
+bk+1 ) . Also, for |h22 | ≥ bk +b
q
bk+1 |h11 ||h22 |
(bk +bk+1 )|h11 ||h22 |
3p1
3p1
tain ak+1 |h21 |+|h22 | M 2 −1 ≥ (ak +ak+1 )|h21 |
. In
2
qM −1
√
p (a +ak+1 )
3p2
11 |
√ 2 k
addition, for |h
|h21 | ≤
p1 (bk +bk+1 ) , we have |h22 |
M 2 −1 ≥

Now, we aim to combine g21

and g22

14

q

(bk +bk+1 )|h11 ||h22 |
3p1
(ak +ak+1 )|h21 |
M 2 −1 . In conclusion,


b
b
≥ g22 abkk , ak+1
.
bk+1 , g11 abkk , ak+1
k+1
k+1

for

|h21 |
|h22 |

≥ bk +

and
 (a

else
ℓ̃b +aℓ̃b +1 )|h21 |

(bℓ̃ +bℓ̃

With the above discussion, we have the result in Lemma 1
and we complete the proof.


b

b +1

)|h11 |

we
q

3p2
M 2 −1 ,

This completes the proof.

have
q

3p2
M 2 −1

∗
∗
(w1 , w2 )
.

=


G. Proof of Theorem 3
F. Proof of Theorem 1

1) If

For the weak cross link, we have
k = 1, . . . , |SM−1 |.
√
p
√2 ,
p1

|h21 |
|h22 |

≤ 1 ≤ bk + bk+1 for
√
p (a +ak+1 )
√ 2 k
p1 (bk +bk+1 ) ,

bk bk+1
=
ak , ak+1

|h11 |

√
p |h |
√ 2 22
p1 |h21 |
√
p2 a
√ ( k+1
p1 bk+1

√
p
√ 2 |h22 |,
p1
√
p a|SM −1 |+1
√ 2(
p1 b|SM −1 |+1
|h22 |
bk+1 |h21 | ) for

≤
=

we
+

have
|h22 |

|h11 |
|h21 |

b|SM −1 |+1 |h21 | )

≤
≤

1, . . . , |SM−1 |.

b
=
According to Lemma 1, we have g abkk , ak+1
k+1
q
and then by Lemma 1, we attain g
b
|h
||h
|
q
3p1
11
22
k+1
|h11 |
3p1
ak+1 |h21 |+|h22 |
M 2 −1 , for k = 1, . . . , |SM−1 |. We
for k
=
1, . . . , |SM−1 |.
ak +ak+1
M 2 −1
also have
Also, note that |SM−1 | = arg mink {(a1 +
n b |h ||h | r 3p
a2 ), . .n. , (a|SMq−1 |
+
a|SM −1 |+1 )}, q theno
2 11
22
1
max
,...,
2−1
|h11 |
|h11 |
3p1
3p1
a
|h
|
+
|h
|
M
2
21
22
=
max a1 +a2 M 2 −1 , · · · , a|S
2
+a|SM −1 |+1
M −1
r
r
M −1 |
q
b|SM −1 |+1 |h11 ||h22 |
3p1 o
3p1
3p1
∗
∗
=
|h
|
,
|h11 |
,
with
(w1 , w2 )
=
11
2
a|SM −1 |+1 |h21 | + |h22 | M 2 − 1
M2 − 1
q

q M −1
M|h11 |
3p1
3p1
q
q
M 2 −1 , |h21 |
M 2 −1 .
|h11 |
3p1
1
√
√
, |h
with (w1∗ , w2∗ ) = ( M3p2 −1
M p2
(ak +ak+1 ) p2
|h11 |
|h11 |
M 2 −1 ).
22 |
2) If |h21 | ≥ √p1 , we have |h21 | ≥ (bk +bk+1 )√p1
√
a
a
p
|h22 |
|h22 | 
ℓ4 +1
ℓ4
11 |

√ 2
for
2) If |h
|h21 | ∈
p1 bℓ4 +1 + bℓ4 +1 |h21 | , bℓ4 + bℓ4 |h21 |
=
and then, by using Lemma 1, we have g abkk , abk+1
k+1
q
ℓ
=
1,
.
.
.
,
|S
|.
4
M−1
|h21 |
3p2

for k = 1, . . . , |SM−1 |. Also, note that
bk +bk+1
M22 −1
=
Then, according to Lemma 1, we have g abkk , abk+1
k+1
q
1 = arg mink {(b1 + b2 ), . . . , (b|SM −1 | + b|SM −1 |+1 )}.
bk+1 |h11 ||h22 |
3p1
= 1, . . . , ℓ4 − 1 and
ak+1 |h21 |+|h22 |
M 2 −1 , for k
Then,
q

b
2
n |h | r 3p
for k = ℓ4 , . . . , |SM−1 |.
= |h22 | M3p2 −1
g abkk , ak+1
k+1
2
21
,
·
·
·
,
max

q
2
b1 + b2 M − 1
3p1
 bk+1 |h11 ||h22 |
r
r

o
b
b
k
ak+1 |h21 |+|h22 |
M 2 −1 for k = 1, . . . , ℓ4 − 1,
k+1
|h21 |
3p2
3p2
q
g
,
=
= |h21 |
,
3p2
|h22 |
ak ak+1
b|SM −1 | + b|SM −1 |+1 M 2 − 1
M2 − 1
M 2 −1 for k = ℓ4 , . . . , |SM−1 |.
q
q


Then, we have,
3p2
3p2
21 |
.
,
with (w1∗ , w2∗ ) = M|h
2
2
|h11 |
M −1
M −1 
√
n b |h ||h | r 3p
√
p2 (aℓ1 +1 +aℓ1 +2 )
p2 (aℓ1 +aℓ1 +1 )
|h11 |
2 11
22
1
√
√
∈
,
3) If |h
,
for
ℓ
=
max
,...,
1
p1 (bℓ1 +1 +bℓ1 +2 )
p1 (bℓ1 +bℓ1 +1 )
21 |
a2 |h21 | + |h22 | M 2 − 1
r
r
1, . . . , |SM−1 | − 1, then, with the help of Lemma 1, we
bℓ4 |h11 ||h22 |
3p1
3p2 o
have
,
|h
|
22
2
aℓ4 |h21 | + |h22 | M − 1
M2 − 1
b b

k
k+1
r
r
n b |h ||h |
g
,
3p1
3p2 o
22
ℓ4 11
ak ak+1
= max
.
,
|h
|
22

q
aℓ4 |h21 | + |h22 | M 2 − 1
M2 − 1
3p1
 |h11 |
for
k
=
1,
.
.
.
,
ℓ
,
1
2
√
ak +ak+1

aℓ 4
p
q M −1
|h11 |
22 |
=
√ 2
As a result, if |h
≥
+ bℓ|h|h
,
3p2
 |h21 |
|
p
b
|
21
1
21
ℓ
4q
4
q
bk +bk+1
M 2 −1 for k = ℓ1 + 1, . . . , |SM−1 |.
bℓ4 |h11 ||h22 |
1
2
then aℓ |h21 |+|h22 | M3p
≥
|h22 | M3p2 −1
and
2 −1
4
q
q

Note that, ℓ̃a = arg mink {(a1 + a2 ), . . . , (aℓ1 + aℓ1 +1 )}
bℓ4 |h11 |
3p1
3p1
∗
∗
Else, we
M 2 −1 , aℓ4 |h21 |+|h22 | q
M 2 −1 . q
and ℓ̃b = arg mink {(bℓ1 +1 + bℓ1 +2 ), . . . , (b|SM −1 | + (w1 , w2 ) =

a
|h
|+|h
|
21
22
ℓ4 +1
3p2
3p2
b|SM −1 |+1 )}, hence we have
can attain (w1∗ , w2∗ ) =
bℓ4 +1 |h11 |
M 2 −1 ,
M 2 −1 .
r
n |h | r 3p
This completes the proof.

3p1
|h11 |
11
1
1) If

|h11 |
|h21 |

≤

M

we have

|h11 |
|h21 |

≤

+

k

=

max

,··· ,
,
a1 + a2 M 2 − 1
aℓ1 + aℓ1 +1 M 2 − 1
r
r
|h21 |
|h21 |
3p2
3p2 o
,...,
2
bℓ1 +1 + bℓ1 +2 M − 1
b|SM −1 | + b|SM −1 |+1 M 2 − 1
r
r
n
|h11 |
3p1
3p2 o
|h21 |
= max
.
,
aℓ̃a + aℓ̃a +1 M 2 − 1 bℓ̃ + bℓ̃ +1 M 2 − 1
b

Therefore, if
q

|h11 |
|h21 |

|h11 |
aℓ̃a +aℓ̃a +1

3p1
M 2 −1

(w1∗ , w2∗ )

=

≥
q

≥

b

√
p2 (aℓ̃a +aℓ̃a +1 )
√
p1 (bℓ̃ +bℓ̃ +1 ) ,

qb

we have

b

|h21 |
3p2
bℓ̃ +bℓ̃ +1
M 2 −1 and hence,
b
b
q

(bℓ̃a +bℓ̃a +1 )|h11 |
3p1
3p1
,
2
M −1 (aℓ̃a +aℓ̃a +1 )|h21 |
M−1

R EFERENCES

[1] D. Tse and P. Viswanath, Fundamentals of Wireless Communications.
Cambridge, U.K.: Cambridge Univ. Press, 2005.
[2] L. Dai, B. Wang, Y. Yuan, S. Han, C. l. I, and Z. Wang, “Non-orthogonal
multiple access for 5G: solutions, challenges, opportunities, and future
research trends,” IEEE Commun. Mag., vol. 53, pp. 74–81, Sept. 2015.
[3] W. Shin, M. Vaezi, B. Lee, D. J. Love, J. Lee, and H. V. Poor, “Nonorthogonal multiple access in multi-cell networks: Theory, performance,
and practical challenges,” CoRR, vol. abs/1611.01607, 2016.
[4] Z. Ding, Z. Yang, P. Fan, and H. V. Poor, “On the performance of
non-orthogonal multiple access in 5G systems with randomly deployed
users,” IEEE Signal Process. Lett., vol. 21, pp. 1501–1505, Dec. 2014.

15

[5] 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, pp. 185–191, Feb. 2017.
[6] Z. Ding, P. Fan, and H. V. Poor, “Impact of user pairing on 5G
nonorthogonal multiple-access downlink transmissions,” IEEE Trans.
Veh. Technol., vol. 65, pp. 6010–6023, Aug. 2016.
[7] Y. Saito, A. Benjebbour, Y. Kishiyama, and T. Nakamura, “System-level
performance evaluation of downlink non-orthogonal multiple access
(NOMA),” in Proc. IEEE Personal Indoor and Mobile Radio Commun.
(PIMRC’13), pp. 611–615, Sept. 2013.
[8] T. Cover, “Broadcast channels,” IEEE Trans. Inf. Theory, vol. 18, pp. 2–
14, Jan. 1972.
[9] G. Caire and S. Shamai, “On the achievable throughput of a multiantenna
Gaussian broadcast channel,” IEEE Trans. Inf. Theory, vol. 49, pp. 1691–
1706, July 2003.
[10] A. Wyner, “Recent results in the shannon theory,” IEEE Trans. Inf.
Theory, vol. 20, pp. 2–10, Jan. 1974.
[11] R. S. Cheng and S. Verdu, “Gaussian multiaccess channels with ISI:
capacity region and multiuser water-filling,” IEEE Trans. Inf. Theory,
vol. 39, pp. 773–785, May 1993.
[12] A. Benjebbovu, A. Li, Y. Saito, Y. Kishiyama, A. Harada, and T. Nakamura, “System-level performance of downlink NOMA for future LTE
enhancements,” in Proc. Globecom Workshops (GC Wkshps’13), pp. 66–
70, Dec. 2013.
[13] Y. Saito, Y. Kishiyama, A. Benjebbour, T. Nakamura, A. Li, and
K. Higuchi, “Non-orthogonal multiple access (NOMA) for cellular future radio access,” in Proc. IEEE 77th Vehi. Tech. Conf. (VTC Spring’13),
pp. 1–5, June 2013.
[14] Q. C. Li, H. Niu, A. T. Papathanassiou, and G. Wu, “5G network
capacity: Key elements and technologies,” IEEE Veh. Technol. Mag.,
vol. 9, pp. 71–78, Mar. 2014.
[15] P. Wang, J. Xiao, and L. P, “Comparison of orthogonal and nonorthogonal approaches to future wireless cellular systems,” IEEE Veh.
Technol. Mag., vol. 1, pp. 4–11, Sept. 2006.
[16] Z. Wei, J. Yuan, D. W. K. Ng, M. Elkashlan, and Z. Ding, “A
survey of downlink non-orthogonal multiple access for 5G wireless
communication networks,” CoRR, vol. abs/1609.01856, 2016.
[17] J. G. Andrews, S. Buzzi, W. Choi, S. V. Hanly, A. Lozano, A. C. K.
Soong, and J. C. Zhang, “What will 5G be?,” IEEE J. Sel. Areas
Commun., vol. 32, pp. 1065–1082, June 2014.
[18] 3rd Generation Partnership Project (3GPP), “Study on downlink multiuser superposition transmission (MUST) for LTE (release 13),” 3GPP
TR 36.859 V13.0.0, Dec. 2015.
[19] H. Lee, S. Kim, and J. H. Lim, “Multiuser superposition transmission
(MUST) for LTE-A systems,” in Proc. IEEE Int. Conf. Commun.
(ICC’16), pp. 1–6, May 2016.
[20] Q. Sun, S. Han, C. L. I, and Z. Pan, “On the ergodic capacity of MIMO
NOMA systems,” IEEE Wireless Commun. Lett., vol. 4, pp. 405–408,
Aug. 2015.
[21] Y. Liu, G. Pan, H. Zhang, and M. Song, “On the capacity comparison between MIMO-NOMA and MIMO-OMA,” IEEE Access, vol. 4,
pp. 2123–2129, May 2016.
[22] Z. Ding, M. Peng, and H. V. Poor, “Cooperative non-orthogonal multiple
access in 5G systems,” IEEE Commun. Lett., vol. 19, pp. 1462–1465,
Aug. 2015.
[23] Z. Ding, H. Dai, and H. V. Poor, “Relay selection for cooperative
NOMA,” IEEE Wireless Commun. Lett., vol. 5, pp. 416–419, Aug. 2016.
[24] Y. Liu, Z. Ding, M. Elkashlan, and H. V. Poor, “Cooperative nonorthogonal multiple access with simultaneous wireless information and
power transfer,” IEEE J. Sel. Areas Commun., vol. 34, pp. 938–953, Apr.
2016.
[25] Y. Liu, Z. Ding, M. Elkashlan, and J. Yuan, “Nonorthogonal multiple
access in large-scale underlay cognitive radio networks,” IEEE Trans.
Veh. Technol., vol. 65, pp. 10152–10157, Dec. 2016.
[26] A. Lozano, A. M. Tulino, and S. Verdu, “Optimum power allocation
for parallel Gaussian channels with arbitrary input distributions,” IEEE
Trans. Inf. Theory, vol. 52, pp. 3033–3051, July 2006.
[27] T. Kasami and S. Lin, “Coding for a multiple-access channel,” IEEE
Trans. Inf. Theory, vol. 22, pp. 129–137, Mar. 1976.
[28] T. Kasami and S. Lin, “Bounds on the achievable rates of block coding
for a memoryless multiple-access channel,” IEEE Trans. Inf. Theory,
vol. 24, pp. 187–197, Mar. 1978.
[29] R. Ahlswede and V. Balakirsky, “Construction of uniquely decodable
codes for the two-user binary adder channel,” IEEE Trans. Inf. Theory,
vol. 45, pp. 326–330, Jan. 1999.

[30] J. Harshan and B. Rajan, “On two-user Gaussian multiple access
channels with finite input constellations,” IEEE Trans. Inf. Theory,
vol. 57, pp. 1299–1327, Mar. 2011.
[31] J. Harshan and B. S. Rajan, “A novel power allocation scheme for
two-user GMAC with finite input constellations,” IEEE Trans. Wireless
Commun., vol. 12, pp. 818–827, Feb. 2013.
[32] Y. Wu, C. K. Wen, C. Xiao, X. Gao, and R. Schober, “Linear precoding
for the MIMO multiple access channel with finite alphabet inputs and
statistical CSI,” IEEE Trans. Wireless Commun., vol. 14, pp. 983–997,
Feb 2015.
[33] J. Choi, “On the power allocation for a practical multiuser superposition
scheme in noma systems,” IEEE Commun. Lett., vol. 20, pp. 438–441,
March 2016.
[34] S. L. Shieh and Y. C. Huang, “A simple scheme for realizing the
promised gains of downlink nonorthogonal multiple access,” IEEE
Trans. Commun., vol. 64, pp. 1624–1635, April 2016.
[35] S. Vishwanath, N. Jindal, and A. Goldsmith, “The “Z” channel,” in Proc.
IEEE Global Commu. Conf. (GLOBECOM’03), vol. 3, pp. 1726–1730
vol.3, Dec. 2003.
[36] H. F. Chong, M. Motani, and H. K. Garg, “Capacity theorems for the
“Z” channel,” IEEE Trans. Inf. Theory, vol. 53, pp. 1348–1365, Apr.
2007.
[37] N. Liu and S. Ulukus, “On the capacity region of the Gaussian Zchannel,” in Proc. IEEE Global Commu. Conf. (GLOBECOM’04), vol. 1,
pp. 415–419 Vol.1, Nov. 2004.
[38] V. R. Cadambe, S. A. Jafar, and S. Vishwanath, “The capacity region
of a class of deterministic Z channels,” in Proc. IEEE Int. Symp. Inform. Theory (ISIT’09), pp. 2634–2638, June 2009.
[39] S. Salehkalaibar and M. R. Aref, “On the capacity region of the degraded
Z channel,” in Proc. IEEE Inform. Theory Workshop (ITW’10), pp. 1–5,
Aug. 2010.
[40] J. Lee, K. Kim, S. C. Cho, and S. H. Lee, “Linear beamforming with
one-way cooperation for MISO interference channels,” in Proc. IEEE
77th Vehi. Tech. Conf. (VTC Spring’13), pp. 1–5, June 2013.
[41] R. Duan, Y. Liang, A. Khisti, and S. S. Shitz, “State-dependent Gaussian
Z-channel with mismatched side-information and interference,” in Proc.
IEEE Inform. Theory Workshop (ITW’13), pp. 1–5, Sept. 2013.
[42] N. Liu and W. Kang, “The capacity region of a class of Z channels
with degraded message sets,” IEEE Transactions on Information Theory,
vol. 60, pp. 6144–6153, Oct. 2014.
[43] T. S. Han and K. Kobayashi, “A new achievable rate region for the
interference channel,” IEEE Trans. Inf. Theory, vol. 27, pp. 49–60, Jan.
1981.
[44] M. Costa, “On the Gaussian interference channel,” IEEE Trans. Inf.
Theory, vol. 31, pp. 607–615, Sept. 1985.
[45] E. Kurniawan and S. Sun, “Improper Gaussian signaling scheme for
the Z-interference channel,” IEEE Trans. Wireless Commun., vol. 14,
pp. 3912–3923, July 2015.
[46] R. Kolte, A. Özgür, and H. Permuter, “The capacity region of a class
of deterministic state-dependent Z-interference channels,” in Proc. IEEE
Int. Symp. Inform. Theory (ISIT’14), pp. 656–660, June 2014.
[47] M. Vaezi and H. V. Poor, “Simplified Han-Kobayashi region for onesided and mixed Gaussian interference channels,” in Proc. IEEE Int.
Conf. Commun. (ICC’16), pp. 1–6, May 2016.
[48] S. Lagen, A. Agustin, and J. Vidal, “On the superiority of improper
Gaussian signaling in wireless interference MIMO scenarios,” IEEE
Trans. Commun., vol. 64, pp. 3350–3368, Aug. 2016.
[49] G. Hardy and E. Wright, An Introduction to the Theory of Numbers.
Oxford Univ. Press, 4 ed., 1975.
[50] K. Ramchandran, A. Ortega, K. Uz, and M. Vetterli, “Multiresolution
broadcast for digital HDTV using joint source/channel coding,” IEEE J.
Sel. Areas Commun., vol. 11, pp. 6–23, Jan. 1993.
[51] A. Goldsmith, Wireless Communications. Cambridge University Press,
2005.
[52] V. R. Cadambe and S. A. Jafar, “Interference alignment and degrees of
freedom of the K-user interference channel,” IEEE Trans. Inf. Theory,
vol. 54, pp. 3425–3441, Aug. 2008.
[53] R. H. Etkin, D. Tse, and H. Wang, “Gaussian interference channel
capacity to within one bit,” IEEE Trans. Inf. Theory, vol. 54, pp. 5534–
5562, Dec. 2008.
[54] G. Bresler, A. Parekh, and D. Tse, “The approximate capacity of the
many-to-one and one-to-many Gaussian interference channels,” IEEE
Trans. Inf. Theory, vol. 56, pp. 4566–4592, Sept. 2010.
[55] V. R. Cadambe, S. A. Jafar, and C. Wang, “Interference alignment with
asymmetric complex signaling: Settling the Høst-Madsen-Nosratinia
conjecture,” IEEE Trans. Inf. Theory, vol. 56, pp. 4552–4565, Sept 2010.

16

[56] V. R. Cadambe, S. A. Jafar, and S. Shamai, “Interference alignment on
the deterministic channel and application to fully connected Gaussian
interference networks,” IEEE Trans. Inf. Theory, vol. 55, pp. 269–274,
Jan. 2009.
[57] R. Prasad, S. Bhashyam, and A. Chockalingam, “On the Gaussian manyto-one X channel,” IEEE Trans. Inf. Theory, vol. 62, pp. 244–259, Jan.
2016.
[58] R. Hoshyar, F. P. Wathan, and R. Tafazolli, “Novel low-density signature
for synchronous CDMA systems over AWGN channel,” IEEE Trans.
Signal Process., vol. 56, pp. 1616–1626, Apr. 2008.
[59] T. M. Cover and J. A. Thomas, Elements of Information Theory. WileyInterscience, 2nd ed., 2006.

