STUDY OF BUFFER-AIDED PHYSICAL-LAYER NETWORK CODING WITH ML AND
MMSE LINEAR CODE DESIGNS FOR COOPERATIVE DS-CDMA SYSTEMS
Jiaqi Gu z, Rodrigo C. de Lamare ♣, z and Mario Huemer ♠

arXiv:1707.06333v1 [] 20 Jul 2017

z Communications Research Group, Department of Electronics, University of York, United Kingdom
♣ CETUC, PUC-Rio, Rio de Janeiro, Brazil
♠ Johannes Kepler University Linz, Austria
Emails: jg849@york.ac.uk, rodrigo.delamare@york.ac.uk, Mario.Huemer@jku.at
ABSTRACT
In this work, we propose buffer-aided physical-layer network coding
(PLNC) techniques for cooperative direct-sequence code-division
multiple access systems. In the proposed buffer-aided PLNC
schemes, a relay pair selection algorithm is employed to obtain
the relay pair and the packets in the buffer entries with the best
performance and the associated link combinations used for data
transmission. We also devise PLNC schemes based on the design
of linear network coding matrices using maximum likelihood and
minimum mean-square error criteria to generate the network coded
symbols that are sent to the destination. Simulation results show that
the proposed techniques significantly outperform prior art.
Index Terms— DS–CDMA networks, cooperative communications, physical-layer network coding, buffer-aided schemes.
1. INTRODUCTION
Interference mitigation is one of the key problems in wireless communications and many approaches have been devised for this task
in the last decades or so. Unlike traditional strategies that treat interference as a nuisance to be avoided, physical-layer network coding (PLNC) techniques take advantage of the superposition of radio
signals and embrace the interference to improve throughput performance [1]. PLNC techniques have generated a number of fertile theoretical and application-oriented studies, and are foreseen to be successfully implemented in future wireless applications [2, 3, 4, 5, 6].
In wireless and cognitive radio networks, the physical superposition
of the signals can be seen as a benefit rather than an interference
and exploited for coding at the physical-layer level [7]. The communication can be protected against attacks from malicious nodes
[8], eavesdropping entities [9], and impairments such as noise and
information losses [10] thanks to the property of the network acting
as an encoder. In peer-to-peer networks, the well-known problem of
the missing block at the end of the download can be avoided by the
distribution of a number of encoded versions of the source data [11].
PLNC has significant advantages in wireless multi-hop networks. Multiple relay nodes are employed in the network to transmit
data from sources to the destination [12]. It allows a node to exploit
as far as possible all signals that are received simultaneously, rather
than treating them as interference [12]. Additionally, instead of
decoding each incoming data stream separately, a node detects and
forwards the function of the incoming data streams [13]. There are
several different network coding techniques, namely, the XOR mapping schemes and linear network coding designs [1, 14, 2, 15, 3, 6].
In cooperative relaying systems, strategies that employ relays
have been investigated in [16, 17, 18, 19, 20, 21, 22]. In order to
This work is funded by the ESII consortium in the UK under Task 26 for
low-cost wireless ad hoc and sensor networks.

further increase the quality and reliability of cooperative schemes,
the concept of buffer has been introduced and used to equip relay
nodes in cooperative relaying scenarios [23, 24, 25, 26, 27, 28]. In
this context, relay selection algorithms must be employed to obtain
the combination of links in the buffer entries that optimize some criterion and result in performance improvements.
In this work, we propose buffer-aided PLNC techniques for
cooperative direct-sequence code-division multiple access (DSCDMA) systems. In the proposed buffer-aided PLNC schemes, a
relay pair selection algorithm is employed to obtain the relay pair
and the packets in the buffer entries with the highest SINR and
the associated link combinations used for data transmission. We
also devise PLNC schemes based on the design of linear network
coding matrices using maximum likelihood (ML) and minimum
mean-square error (MMSE) criteria to generate the network coded
symbol (NCS) that is then sent to the destination. Simulation results
show that the proposed techniques significantly outperform prior art.
This paper is structured as follows: Section 1 introduces the cooperative DS-CDMA system model. Section 2 presents the proposed
ML and MMSE linear network code designs, whereas Section 3 describes the proposed buffer-aided PLNC transmission scheme. Section 4 presents the simulations and Section 5 gives the conclusions.
2. COOPERATIVE DS-CDMA SYSTEM MODEL
Source
Relay
k1
l1

k2

Destination
l2

k3

.
.
.
kK

D

.
.
.
lL

Direct link
Source-relay
Relay-destination

Fig. 1. Uplink of a cooperative DS-CDMA system.
We consider the uplink of a two-phase cooperative synchronous
DS-CDMA system with K users, L relays equipped with finite-size
buffers capable of storing J packets and N chips per symbol that experiences channels with flat fading, as shown in Fig. 1. The system
is equipped with a cooperative protocol at each relay and the transmit data are organized in packets of P symbols. The received signals
are filtered by a matched filter and sampled at chip rate to obtain sufficient statistics and organized into N × 1 vectors ysr , ysd and yrd ,

which represent the signals received from the sources to the relays,
the sources to the destination and the relay to the destination, respectively. Every set of m users and m relays are randomly assigned into
a group, forming pairs of users and relays. The signal received from
the remaining users to the target relays are seen as interference.
In the first phase, the signals received at the destination and the
l-th relay can be described as
K
X
aksd sk hsd,k bk + nsd ,
(1)
ysd =
k=1

ysrl =

K
X

aksrl sk hsrl ,k bk

+ nsrl ,

(2)

k=1

where bk correspond to the transmitted symbols. The quantities aksd
and aksrl are the k-th user’s amplitude from the source to the destination and from the source to relay l, respectively. The vector
sk = [sk (1), sk (2), ...sk (N )]T is the signature sequence for user
k. The quantities hsd,k and hsrl ,k are the complex channel coefficients from user k to the destination and from user k to the l-th relay,
respectively. The vectors nsd and nsrl are the noise vectors that contain samples of zero mean complex Gaussian noise with variance σ 2 .
After the detected symbols of interest are generated, a PLNC
scheme is then employed at each relay. When the bit-wise XOR
operation (the modulo-2 sum in the binary field) is considered, we
assume that the mapping from user symbol to a modulated symbol
is denoted as MB . This mathematical mapping is then expressed as
bk = MB = 1 − 2ck .
(3)
We then perform the above mapping from b̂rl d,k , which represents
the decoded symbol for user k at the output of relay l using a relay
protocol, to its corresponding user symbol as described by
(4)
b̂rl d,k → ckl .
Similarly, when linear network coding is adopted at the l-th relay,
the corresponding mapping LB is given by
K
X
gkl b̂rl d,k .
(5)
bl = L B =
k=1

In the following, we briefly review existing PLNC schemes that employ linear network codes at each relay. When the signals are sent to
the relays, we obtain the detected
 symbol as given by
b̂rΩ(l) d,Υ(k) = Q (wsΥ(k) rΩ(l) )H ysrl ,

(6)

where Υ is a 1 × m user set and Ω is the corresponding 1 × m
relay set since we combine every m users and m relays into a sub
transmission group. b̂rΩ(l) d,Υ(k) is the detected result for user Υ(k)
at relay Ω(l), wsΥ(k) rΩ(l) refers to the detector for user Υ(k) at relay
Ω(l). After that, the detected symbol is mapped through (3) to obtain
Υ(k)
cΩ(l) = (1 − b̂rΩ(l) d,Υ(k) )/2.
(7)
If linear network coding techniques [29, 30, 31] are adopted at the
relays, every relay group will be allocated a unique m × m linear
network coding matrix G whose structure is given by


g11
g12 . . . g1m
 g21
g22 . . . g2m 


(8)
G=
.
..


.
gm1 gm2 . . . gmm
Thus, the following linear combination can be obtained as
m
X
(9)
gkl b̂rΩ(l) d,Υ(k) ,
bΩ(l) =
k=1

where b̂rΩ(l) d,Υ(k) is the detected symbol for user Υ(k) at relay
Ω(l). After that, in order to forward the encoded information to

the destination, the following mapping operation is required
bΩ(l) = 1 − 2cΩ(l) .
(10)
Consequently, in the second phase the encoded symbols are stored
and prepared to be sent to the destination. The NCS are then transmitted from relay set Ω to the destination as described by
m
X
hrΩ(l) d bΩ(l) + nrd ,
(11)
yrΩ d =
l=1

where hrΩ(l) d = aΥ
rΩ(l) d sΥ hrΩ(l) d denotes the N × 1 channel vector from relay Ω(l) to the destination, aΥ
rΩ(l) d is the amplitude for
combined source (user set Υ) from the Ω(l)-th relay to the destination, hrΩ(l) d is the complex channel fading coefficient from the
Ω(l)-th relay to the destination, sΥ is the N × 1 spreading sequence
assigned to the NCS bΩ(l) and nrd is the N × 1 zero mean complex
Gaussian noise with variance σ 2 .
At the destination, the detected symbol for users from set Υ after
linear network coding at each
by
 of the relays is obtained

b̂Ω(l) = Q G−1 (wrΩ d )H yrΩ d .

(12)

The decoding process for linear network coding is expressed by
b̂Ω(1) = gΥ(1)Ω(1) b̂Υ(1) + gΥ(2)Ω(1) b̂Υ(2) + ... + gΥ(m)Ω(1) b̂Υ(m)
..
.
b̂Ω(m) = gΥ(1)Ω(m) b̂Υ(1) + gΥ(2)Ω(m) b̂Υ(2) + ... + gΥ(m)Ω(m) b̂Υ(m)
(13)
The m unknown detected symbols b̂Υ(k) , k = 1, 2, ...m for each
user in user set Υ can be obtained by solving the above m equations. However, in this scenario, since we have direct transmission,
there is another decoding alternative. Specifically, we can obtain the
detected symbol for the user b̂Υ(k) as given by:
Υ(1)
Υ(m) 
b̂Ω(l) − gΥ(1)Ω(l) b̂sd + . . . + gΥ(m)Ω(l) b̂sd
b̂Υ(k) =
. (14)
gΥ(k)Ω(l)
With direct transmissions, we can effectively reduce the bit error rate. In particular, there is a possibility that the detection for
b̂Ω(l) , l ∈ Ω is incorrect, which could directly lead to a problem
when solving the system of equations in (13), causing incorrect
decisions. However, this is not a problem when applying (14) as we
can select another detected symbol b̂Ω(s) , s 6= l as the minuend in
(14) instead of the possibly incorrect detected symbol b̂Ω(l) .

3. PROPOSED COOPERATIVE LINEAR NETWORK
CODING SCHEMES
In this section, we present novel linear network coding designs for
G. The idea is inspired by the work in [32] but applied in a completely different context. In the cooperative PLNC transmission scenario, the simplest way to construct G is to generate it randomly
[33], however, this approach often does not result in the optimum
performance as it does not satisfy an optimality criterion. Therefore,
in order to obtain performance advantages, we devise methods based
on the ML and MMSE criteria.

3.1. Maximum Likelihood Design
When a linear network coding scheme is adopted at the relay, we
introduce the m × m linear network code matrix G, which is used to
generate the NCS. An ML design of G that requires the evaluation

2

of all possible (2m ) binary matrices is given by

 H


wrΩ(1) d yrΩ(1) d
bΩ(1)



 .
..
−1

GML = arg min
−G 
 ..


.
2
1≤j≤2m
H
bΩ(m)
wrΩ(m) d yrΩ(m) d

2

,

(15)
where wrl d is the linear detector used at relay l, yrl d is the signal
transmitted from relay l. In summary, the principle of this approach
is to obtain an m×m code matrix GML to match G (or GML ). Thus,
the minimum distance between the transmitted symbols and detected
symbols can be obtained. Note that this approach corresponds to a
2
combinatorial problem where we must test all 2m possible candiML
dates to obtain G associated with the best performance. The computational cost is one of the disadvantages of this approach but the
design of GML can be carried out online. At the destination side, we
employ the same matrix GML that is applied at the relays in order to
obtain better detection performance.
3.2. MMSE design
The ML approach in (15) does not consider the influence brought
by noise. This fact motivates us to seek another efficient decoding method that can consider both interference and noise. In order
to further exploit the minimum distance and improve the transmission performance, we introduce a code matrix G̃ at the destination
to match the generated binary matrix G and perform the linear network coding operation. The proposed MMSE linear network coding
matrix GMMSE is obtained
 H
 as follows:
 

wrΩ(1) d yrΩ(1) d
bΩ(1)

  .
 2

..
 ,
GMMSE = arg min E 
 ..
 − G̃ 


 
.
G̃
bΩ(m)
wH
rΩ(m) d yrΩ(m) d
(16)

Υ(k)
bΩ(1)

Υ(k)
bΩ(2)

Υ(k)
bΩ(m) , k

where ideally we have
=
= ... =
This design problem can be recast as the following:
h
2i
GMMSE = arg min E a − G̃b
,

∈ [1, m].
(17)

G̃

where the quantities in the argument are a = [bΩ(1) . . . bΩ(m) ]T and
T
H
b = [wH
rΩ(1) d yrΩ(1) d . . . wrΩ(m) d yrΩ(m) d ] . By taking the gradient
of the cost function with respect to G̃ and equating the terms to zero,
we obtain the MMSE linear network code matrix given by
(18)
GMMSE = Pab R−1
b ,
where the statistical quantities in the MMSE code matrix are the
cross-correlation matrix Pab = E[abH ] and the auto-correlation
matrix Rb P
= E[bbH ]. The elements of Pab are given by
2
H
[Pab ]k,j = ( m
i=1 gik gik σi )hrΩ(j)d wrΩ(j)d , for k, j = 1, . . . , m,
whereas the main diagonal entries of Rb are described by [Rb ]j,j =
2
wrHΩ(j)d (hrΩ(j)d hH
rΩ(j)d + σ I)wrΩ(j)d , for j = 1, . . . , m.
4. PROPOSED BUFFER-AIDED PLNC SCHEME
In this section, we consider groups of m = 2 users and m = 2
relays and develop a buffer-aided PLNC scheme, where each relay
is equipped with a buffer so the received data can be stored and wait
until the link pair associated with the best performance is selected.
Consequently, encoded data are stored at the buffer entries and then
forwarded to the destination when the appropriate time interval
comes. Moreover, the destination is also equipped with a buffer

so that the detected symbols from the direct transmissions can be
stored. After that, the detected symbols are obtained through PLNC
decoding and mapping operations at the appropriate time instants.
Transmission between a selected user-relay pair (m = 2)
At the relays, the received signal is processed by linear network
coding and forwarded to the destination, which employs a decoding
matrix and a linear receiver. Specifically, the proposed algorithm
starts with a selection procedure using all possible link combinations associated with all relay pairs from both source-relay and
relay-destination phases. Since any linear network coding technique
can be adopted, every (m = 2) relays are combined into a group
and paired with (m = 2) users. Note that arbitrary numbers of users
and relays can also be considered but for the sake of simplicity,
we only consider m = 2 relays and users. In this case, the signal
transmitted from other users are seen as the interference component.
All possible links are considered and their SINR are then calculated:
K P
m
P
wH
sk rΩ(l) ρsk rΩ(l) wsk rΩ(l)
k=1 l=1
SINRsrΩ = K L
,
m
P
P P H
σ 2 wH
wsk rj ρsk rj wsk rj +
sk rΩ(l) wsk rΩ(l)
k=1 j=1
j6∈Ω

l=1

(19)
SINRrΩ d =

K P
m
P

(wkrΩ(l) d )H ρkrΩ(l) d wkrΩ(l) d
k=1 l=1
K P
L
m
P
P
σ 2 (wkrΩ(l) d )H wkrΩ(l) d
(wkrj d )H ρkrj d wkrj d +
k=1 j=1
l=1
j6∈Ω

(20)
where ρsk rΩ(l) = hH
sk rΩ(l) hsk rΩ(l) is the correlation coefficient of
the desired user k between the source and relay Ω(l), ρkrΩ(l) d =
(hkrΩ(l) d )H hkrΩ(l) d is the correlation coefficient for user k from relay
Ω(l) to the destination. hsk rΩ(l) = ask rΩ(l) sk hsk rΩ(l) is the channel vector from user k to relay Ω(l). In (19), SINRsrΩ denotes the
SINR for the combined paths from all users to relay set Ω, wsk rΩ(l)
is the detector used at the relay Ω(l). When the RAKE receiver is
adopted at the corresponding relay, wsk rΩ(l) is expressed as
wsk rΩ(l) = hsk rΩ(l) .
(21)
Similarly, if the linear MMSE receiver [34, 35, 36, 37, 38] is employed at the relays, wsk rΩ(l) is equal to
X
−1
K
H
2
wsk rΩ(l) =
hsk rΩ(l) hsk rΩ(l) + σ I
hsk rΩ(l) ,
(22)
k=1

hsk rΩ(l) = ask rΩ(l) sk hsk rΩ(l) is the effective signature vector from
user k to the relay Ω(l). In (20), SINRrΩ d represents the SINR for
the combined paths from relay set Ω to the destination. The receive
filter wkrΩ(l) d is employed by the detector used at the destination.
When the RAKE receiver is adopted at the destination, wkrΩ(l) d is
expressed as
wkrΩ(l) d = hkrΩ(l) d .
(23)
In an analog way, if the linear MMSE receiver is employed at the
relays, wkrΩ(l) d is equal to
−1
X
K
k
k
H
2
k
hrΩ(l) d (hrΩ(l) d ) + σ I
wrΩ(l) d =
hkrΩ(l) d .
(24)
k=1

Both RAKE and MMSE receivers are considered here due to their
reasonably low complexity and it should be mentioned that other
detectors [39, 40, 41, 42, 43] including the ML detector can also

,

5. SIMULATIONS

be used. Alternatively, transmit processing techniques can also be
employed [44, 45, 46, 47].
After the computation of receive filters and SINR values, we sort
all these results according to a decreasing power level and choose the
relay pair (m = 2) with the highest SINR as given by
(25)
SINRi,j = arg max {SINRsrΩ , SINRrΩ d },
Ω∈[1,2,...,L]

where SINRi,j denotes the highest SINR associated with the relay i
and relay j. After the relay pair (m = 2) with the highest SINR is
selected, the signal is transmitted through the corresponding channels. In this case, two different situations of the buffer mode need to
be considered as follows.
Transmission mode:
If the link combinations associated with the selected relay set belongs to the relay-destination phase, the buffers are turned to the
transmission mode. A buffer space check needs to be conducted first
to ensure the corresponding buffer entries are not empty, namely:
Φbuffer
6= ∅, i ∈ [1, 2, ..., L],
(26)
i
and
6= ∅, j ∈ [1, 2, ..., L],
Φbuffer
j

In this section, a simulation study of the proposed buffer-aided
PLNC scheme and linear network coding techniques is carried out.
The DS-CDMA network uses randomly generated spreading codes
of length N = 16. The corresponding channel coefficients are modeled as complex Gaussian random variables. We assume perfectly
known channels at the relays and receivers and remark that results
with channel estimation have the same performance hierarchy. We
consider K = 6 users, L = 6 relays, equal power allocation and
packets with 1000 BPSK symbols in the transmissions. For PLNC
techniques, we first compare the XOR against linear network coding
schemes using different designs of the matrix G with and without
buffers (J = 4) in Fig. 2. The results show that the use of buffers
can provide a significant gain in performance and that linear network
coding techniques outperform XOR-based approaches.
0

10

(27)
−1

where Φbuffer
and Φbuffer
denote the buffers equipped at relay i and
i
j
relay j. If the corresponding buffer entries are not empty, we transmit the NCS according to (11). On the other hand, if the buffer
condition does not satisfy the transmission requirements, namely,
the buffer entries are empty, the selected relay pair cannot help to
forward the NCS.
In this case, we drop the current relay pair i and j and choose
another relay pair with the second highest SINR as given by
SINRpre
(28)
i,j = SINRi,j
SINRu,v ∈ max{SINRsrΩ , SINRrΩ d } \ SINRpre
i,j ,

(29)

BER

10

−2

10

PLNC with XOR
Buffer−aided PLNC with XOR
Random

PLNC with G

random

Buffer−aided PLNC with G
−3

10

5

10

15

20

25

SNR (dB)

Fig. 2. Performance of linear network coding techniques with and
without buffers using BPSK modulation.

where SINRu,v denotes the second highest SINR associated with
the updated relay pair u and v. {SINRsrΩ , SINRrΩ d } \ SINRpre
i,j
represents a complementary set, where we drop the SINRpre
i,j from
the set of SINR links {SINRsrΩ , SINRrΩ d }. Consequently, the
above process repeats until the buffer condition achieves the transmission requirement.

0

10

−1

10

BER

Reception mode:

0

−2

10

random

If the link combination associated with the highest relay pair SINR
belongs to the source-relay phase, the buffers are switched to reception mode. Similarly, the buffer is checked to ensure that the
corresponding buffer entries are not full. We then have
Φbuffer
6= U, i ∈ [1, 2, ..., L],
(30)
i
and

PLNC with G

ML

PLNC with G
−3

MMSE

10

PLNC with G

random

Buffer−aided PLNC with G

ML

Buffer−aided PLNC with G

MMSE

Buffer−aided PLNC with G
−4

10

0

5

10

15

20

25

SNR (dB)

(31)

Fig. 3. Performance of linear network coding techniques with and
without buffers using BPSK modulation.

where U represents a full buffer set. If the buffers are not full, then
the sources send the data to the selected relay pair i and j according
to (2). Otherwise, the algorithm reselects a new relay pair as in (28)
and (29). The re-selection process stops when the buffer entries are
not full. The size J of buffers plays an important role in the system
performance. When we increase the buffer size, better channels can
be selected as a relatively larger candidate pool is generated. Therefore, extra degrees of freedom in the system are also available.

The BER performance of buffer-aided (J = 4) PLNC schemes
with the proposed and existing linear network coding techniques is
shown in Fig. 3. The results show that the introduction of GMMSE
provides the best performance, followed by the introduction of GML
and Grandom . In particular, the adoption of GMMSE results in a gain
of up to 1.5dB in SNR over GML and up to 3 dB in SNR as compared to Grandom for the same BER performance. The results also

Φbuffer
6= U, j ∈ [1, 2, ..., L],
j

demonstrate that with the use of buffers the overall system performance significantly improves, achieving a gain of up to 5 dB in terms
of SNR for the same BER performance as compared to the system
without buffers even though the diversity order is the same.
6. CONCLUSIONS
We have presented a buffer-aided PLNC scheme and devised novel
ML and MMSE linear network coding algorithms for cooperative
DS-CDMA systems with different relay pair selection techniques.
Simulation results show that the performance of the proposed
scheme and algorithms can offer significant gains as compared
to previously reported techniques.
7. REFERENCES
[1] S. Zhang, S. C. Liew, and P. K. Lam, “Hot topic: Physical-layer network coding,” in Proc. Annual Int. Conf. on Mobile Computing and
Networking(MobiCom)., LA, USA, Sep. 2006.
[2] M. Sanna and E. Izquierd, “A survey of linear network coding and network error correction code constructions and algorithms,” International
Journal of Digital Multimedia Broadcasting, vol. 2011, 2011.
[3] J. He and S. C. Liew, “Building blocks of physical-layer network coding,” IEEE Transactions on Wireless Communications, vol. 14, no. 5,
pp. 2711–2728, May 2015.
[4] L. Yang, T. Yang, Y. Xie, J. Yuan, and J. An, “Linear physical-layer network coding and information combining for the k-user fading multipleaccess relay network,” IEEE Transactions on Wireless Communications, vol. 15, no. 8, pp. 5637–5650, Aug 2016.
[5] M. Nokleby and B. Aazhang, “Cooperative compute-and-forward,”
IEEE Transactions on Wireless Communications, vol. 15, no. 1, pp.
14–27, Jan 2016.
[6] H. Zhang, L. Zheng, and L. Cai, “Design and analysis of heterogeneous
physical layer network coding,” IEEE Transactions on Wireless Communications, vol. 15, no. 4, pp. 2484–2497, April 2016.
[7] P. A. Chou and Y. Wu, “Network coding for the internet and wireless
networks,” in Tech. Rep.,Microsoft Research, 2007.
[8] S. Jaggi, M. Langberg, S. Katti, D. K. T. Ho, M. Medard, and M. Effros,
“Resilient network coding in the presence of byzantine adversaries,”
IEEE Trans. Inf. Theory, vol. 54, no. 6, pp. 2596–2603, June 2008.
[9] N. Cai and T. Chan, “Theory of secure network coding,” Proceedings
of the IEEE, vol. 99, no. 3, pp. 421–437, 2011.
[10] R. W. Yeung and N. Cai, “Network error correction. I. basic concepts
and upper bounds,” Communications in Information & Systems, vol. 6,
no. 1, pp. 19–35, 2006.
[11] K. Jain, L. Lovasz, and P. A. Chou, “Building scalable and robust peerto-peer overlay networks for broadcasting using network coding,” Distributed Computing, vol. 19, no. 4, pp. 301–3115, 2007.
[12] A. Burr and D. Fang, “Linear physical layer network coding for multihop wireless networks,” in 22nd European Signal Processing Conference (EUSIPCO), September 2014, pp. 1153–1157.
[13] ——, “Linear physical-layer network coding for 5G radio access networks,” in 5G for Ubiquitous Connectivity (5GU), 2014 1st International Conference on, November 2014, pp. 116–121.
[14] R. Ahlswede, N. Cai, S. Y. R. Li, and R. W. Yeung, “Network information flow,” IEEE Trans. Inf. Theory., vol. 46, no. 4, p. 12041216,
2000.

[17] P. Clarke and R. C. de Lamare, “Transmit diversity and relay selection
algorithms for multi-relay cooperative MIMO systems,” IEEE Trans.
Veh. Technol, vol. 61, no. 3, pp. 1084–1098, Mar 2012.
[18] S. Talwar, Y. Jing, and S. Shahbazpanahi, “Joint relay selection and
power allocation for two-way relay networks,” IEEE Signal Process.
Lett, vol. 18, no. 2, pp. 91–94, Feb 2011.
[19] R. C. De Lamare, “Joint iterative power allocation and linear interference suppression algorithms for cooperative ds-cdma networks,” IET
communications, vol. 6, no. 13, pp. 1930–1942, 2012.
[20] J. Gu and R. de Lamare, “Joint sic and multi-relay selection algorithms
for cooperative DS-CDMA systems,” in Proceedings of the 22nd European Signal Processing Conference (EUSIPCO), Sept 2014, pp. 556–
560.
[21] T. Peng, R. de Lamare, and A. Schmeink, “Adaptive distributed spacetime coding based on adjustable code matrices for cooperative MIMO
relaying systems,” IEEE Trans. on Commun., vol. 61, no. 7, pp. 2692–
2703, July 2013.
[22] J. Gu and R. C. de Lamare, “Joint interference cancellation and relay selection algorithms based on greedy techniques for cooperative
ds-cdma systems,” EURASIP Jounal on Wireless Communications and
Networking, 2016.
[23] N. Zlatanov, A. Ikhlef, T. Islam, and R. Schober, “Buffer-aided cooperative communications: opportunities and challenges,” IEEE Communications Magazine, vol. 52, no. 4, pp. 146–153, May 2014.
[24] N. Zlatanov, R. Schober, and P. Popovski, “Throughput and diversity
gain of buffer-aided relaying,” IEEE Globecom, pp. 1–6, Dec 2011.
[25] I. Krikidis, T. Charalambous, and J. Thompson, “Buffer-aided relay
selection for cooperative diversity systems without delay constraints,”
IEEE Trans. Wireless Commun., vol. 11, no. 5, pp. 1592–1967, May
2012.
[26] A. Ikhlef, D. S. Michalopoulos, and R. Schober, “Max-max relay selection for relays with buffers,” IEEE Trans. Wireless Commun., vol. 11,
no. 5, pp. 1124–1135, Jan 2012.
[27] T. Peng and R. C. de Lamare, “Adaptive buffer-aided distributed spacetime coding for cooperative wireless networks,” IEEE Transactions on
Communications, vol. 64, no. 5, pp. 1888–1900, May 2016.
[28] J. Gu, R. C. de Lamare, and M. Huemer, “Buffer-aided physical-layer
network coding with optimal linear code designs for cooperative networks,” IEEE Transactions on Communications, 2017.
[29] Y. Cai, R. C. De Lamare, and R. Fa, “Switched interleaving techniques
with limited feedback for interference mitigation in ds-cdma systems,”
IEEE Transactions on Communications, vol. 59, no. 7, pp. 1946–1956,
2011.
[30] A. G. Uchôa, C. Healy, R. C. de Lamare, and R. D. Souza, “Design
of ldpc codes based on progressive edge growth techniques for block
fading channels,” IEEE Communications Letters, vol. 15, no. 11, pp.
1221–1223, 2011.
[31] C. T. Healy and R. C. de Lamare, “Design of ldpc codes based on multipath emd strategies for progressive edge growth,” IEEE Transactions
on Communications, vol. 64, no. 8, pp. 3208–3219, 2016.
[32] M. Huemer, C. Hofbauer, and J. B. Huber, “Non-systematic complex
number rs coded ofdm by unique word prefix,” IEEE Trans. Signal Processing, vol. 60, no. 1, pp. 285–299, Jan. 2012.
[33] T. Ho, M. Medard, R. Koetter, D. R. Karger, M. Effros, J. Shi, and
B. Leong, “A random linear network coding approach to multicast,”
IEEE Trans. Inf. Theory, vol. 52, no. 10, pp. 4413–4430, 2006.

[15] S. Y. R. Li, Q. T. Sun, and Z. Shao, “Linear network coding: Theory
and algorithms,” Proceedings of the IEEE, vol. 99, no. 3, pp. 372–387,
March 2011.

[34] R. de Lamare and R. Sampaio-Neto, “Adaptive reduced-rank processing based on joint and iterative interpolation, decimation, and filtering,”
IEEE Trans. Signal Processing, vol. 57, no. 7, pp. 2503–2514, July
2009.

[16] Y. Jing and H. Jafarkhani, “Single and multiple relay selection schemes
and their achievable diversity orders,” IEEE Trans. Wireless Commun.,
vol. 8, no. 3, pp. 1084–1098, Mar 2009.

[35] R. C. De Lamare and R. Sampaio-Neto, “Reduced-rank adaptive filtering based on joint iterative optimization of adaptive filters,” IEEE
Signal Processing Letters, vol. 14, no. 12, pp. 980–983, 2007.

[36] ——, “Reduced-rank space–time adaptive interference suppression
with joint iterative least squares algorithms for spread-spectrum systems,” IEEE transactions on vehicular technology, vol. 59, no. 3, pp.
1217–1228, 2010.
[37] A. Uchoa, C. Healy, and R. de Lamare, “Iterative detection and decoding algorithms for mimo systems in block-fading channels using ldpc
codes,” IEEE Transactions on Vehicular Technology, 2016.
[38] H. Ruan and R. de Lamare, “Robust adaptive beamforming based on
low-rank and cross-correlation techniques,” IEEE Transactions on Signal Processing, vol. 64, no. 15, pp. 3919–3932, 2016.
[39] R. de Lamare and R. Sampaio-Neto, “Minimum mean-squared error
iterative successive parallel arbitrated decision feedback detectors for
ds-cdma systems,” IEEE Trans. on Commun., vol. 56, no. 5, pp. 778–
789, May 2008.
[40] R. C. De Lamare, R. Sampaio-Neto, and A. Hjorungnes, “Joint iterative
interference cancellation and parameter estimation for cdma systems,”
IEEE Communications Letters, vol. 11, no. 12, 2007.
[41] P. Li, R. C. De Lamare, and R. Fa, “Multiple feedback successive interference cancellation detection for multiuser mimo systems,” IEEE
Transactions on Wireless Communications, vol. 10, no. 8, pp. 2434–
2439, 2011.
[42] R. de Lamare, “Adaptive and iterative multi-branch mmse decision
feedback detection algorithms for multi-antenna systems,” IEEE Transactions on Wireless Communications, vol. 12, no. 10, pp. 5294–5308,
2013.
[43] Y. Cai, R. C. de Lamare, B. Champagne, B. Qin, and M. Zhao, “Adaptive reduced-rank receive processing based on minimum symbol-errorrate criterion for large-scale multiple-antenna systems,” IEEE Transactions on Communications, vol. 63, no. 11, pp. 4185–4201, 2015.
[44] R. C. de Lamare, “Massive mimo systems: Signal processing challenges and future trends,” Radio Science Bulletin, 2013.
[45] W. Zhang, H. Ren, C. Pan, M. Chen, R. C. de Lamare, B. Du, and
J. Dai, “Large-scale antenna systems with ul/dl hardware mismatch:
Achievable rates analysis and calibration,” IEEE Transactions on Communications, vol. 63, no. 4, pp. 1216–1229, 2015.
[46] K. Zu, R. de Lamare, and M. Haardt, “Generalized design of lowcomplexity block diagonalization type precoding algorithms for multiuser mimo systems,” IEEE Transactions on Communications, no. 99,
pp. 1–11, 2013.
[47] ——, “Multi-branch tomlinson-harashima precoding design for mumimo systems: Theory and algorithms,” IEEE Transactions on Communications, vol. 62, no. 3, pp. 939–951, 2014.

