Top Banner
1 Synchronous Physical-Layer Network Coding: A Feasibility Study Yang Huang, Student Member, IEEE, Shiqiang Wang, Student Member, IEEE, Qingyang Song, Member, IEEE, Lei Guo, Member, IEEE, and Abbas Jamalipour, Fellow, IEEE Abstract—Recently, physical-layer network coding (PNC) at- tracts much attention due to its ability to improve throughput in relay-aided communications. However, the implementation of PNC is still a work in progress, and synchronization is a signifi- cant and difficult issue. This paper investigates the feasibility of synchronous PNC with M-ary quadrature amplitude modulation (M-QAM). We first propose a synchronization scheme for PNC. Then, we analyze the synchronization errors and overhead of po- tential synchronization techniques, which includes phase-locked loop (PLL) and maximum likelihood estimation (MLE) based synchronization schemes. Their effects on the average symbol error rate and the goodput are subsequently discussed. Based on the analysis, we perform numerical evaluations and reveal that synchronous PNC can outperform conventional network coding (CNC) even when taking synchronization errors and overhead into account. The theoretical throughput gain of PNC over CNC can be approached when using the MLE based synchronization method with optimized training sequence length. The results in this paper provide some insights and benchmarks for the implementation of synchronous PNC. Index Terms—Communication; denoise-and-forward (DNF); physical-layer network coding (PNC); synchronization; two-way relay networks. I. I NTRODUCTION Relay-aided communications are widely adopted when di- rect communications among end nodes cannot be performed. Physical-layer network coding (PNC) [2] is considered as a promising technology to improve the throughput performance of relay networks. It employs the natural network coding ability introduced by the superposition of electromagnetic waves. Between the two methods of PNC, i.e. amplify-and- forward [3] and denoise-and-forward (DNF), the DNF method shows more performance advantages because it avoids noise amplification [4]. Hence, DNF has attracted much interest in Some preliminary ideas of this paper have been presented in IEEE GLOBECOM 2012 [1]. This work was supported in part by the National Natural Science Foundation of China (61172051), the Fok Ying Tung Education Foundation (121065), the Fundamental Research Funds for the Central Universities (N110204001, N110804003, N120804002, N120404001), the Program for New Century Excellent Talents in University (NCET-12-0102), and the Specialized Research Fund for the Doctoral Program of Higher Education (20110042110023, 20110042120035, 20120042120049). Y. Huang, Q. Song (corresponding author) and L. Guo are with School of Information Science and Engineering, Northeastern University, Shenyang 110819, P. R. China. Email: [email protected], {songqingyang, guolei}@ise.neu.edu.cn. S. Wang is with Department of Electrical and Electronic Engineer- ing, Imperial College London, SW7 2AZ, United Kingdom. Email: [email protected]. A. Jamalipour is with School of Electrical and Information Engineering, University of Sydney, NSW, 2006, Australia. Email: [email protected]. recent research, and this paper considers the DNF scheme. We use DNF and PNC interchangeably in subsequent discussions. Research regarding PNC has been carried out focusing on two aspects: phase-asynchronous PNC [5]–[7] and phase- synchronous PNC [8]–[12]. The basic idea of asynchronous PNC is to map the superposed signal with arbitrary phase differences to encoded symbols. However, these schemes require knowledge of the instantaneous phases of the signals superposing at the relay, and imperfect channel information may also degrade the performance of asynchronous PNC [13]. Hence, tracking phase variations 1 during data packet transmission is necessary (although synchronization is not needed), which can be difficult especially for the superposed signal. The complexity of obtaining symbol mapping under various phase differences can also be high, in particular with high-level modulations [5]. Compared with asynchronous PNC, synchronous PNC allows more efficient constellation design [8] and can make use of capacity-approaching channel codes [11]. The capacity region of the Gaussian two-way relay channel can also be reached with synchronous PNC [12]. Further, [14] shows that compared with other schemes, phase synchronization can maximize the minimum distance between adjacent points in the constellation for superposed signals, provided that the signals are with the same modulation and amplitude at the relay. Therefore, we focus on the phase- synchronous PNC in this paper. Synchronization is a significant issue for synchronous PNC, which, however, has not been adequately studied. Existing works [8]–[12] generally assume that superposing signals arrive in-phase at the relay. However, these works have not addressed how to achieve such a synchronization, to the best of our knowledge. Although [15] investigated the impact of imperfect synchronization for binary phase-shift keying (BPSK) modulated PNC, it did not explicitly introduce a synchronization scheme and also did not investigate the interaction between synchronization and data transmission. In the literature, some phase synchronization schemes for distributed beamforming have been studied [16]–[19]. Al- though both PNC and distributed beamforming make use of signal superposition, the goal of PNC is to increase network throughput, while distributed beamforming is for increasing the signal strength at the receiver. Meanwhile, the end nodes cannot communicate with each other when using PNC (oth- erwise relaying is unnecessary); while in beamforming, the 1 Note that frequency errors accumulate over time and may cause the phase difference between the two superposed signals change continuously. © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.
10

Synchronous Physical-Layer Network Coding: A Feasibility Studysw4410/papers/PNCSyncJournal.pdf · 2013. 6. 16. · 1 Synchronous Physical-Layer Network Coding: A Feasibility Study

Jan 26, 2021

Download

Documents

dariahiddleston
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
  • 1

    Synchronous Physical-Layer Network Coding:A Feasibility Study

    Yang Huang, Student Member, IEEE, Shiqiang Wang, Student Member, IEEE, Qingyang Song, Member, IEEE,Lei Guo, Member, IEEE, and Abbas Jamalipour, Fellow, IEEE

    Abstract—Recently, physical-layer network coding (PNC) at-tracts much attention due to its ability to improve throughputin relay-aided communications. However, the implementation ofPNC is still a work in progress, and synchronization is a signifi-cant and difficult issue. This paper investigates the feasibility ofsynchronous PNC with M -ary quadrature amplitude modulation(M -QAM). We first propose a synchronization scheme for PNC.Then, we analyze the synchronization errors and overhead of po-tential synchronization techniques, which includes phase-lockedloop (PLL) and maximum likelihood estimation (MLE) basedsynchronization schemes. Their effects on the average symbolerror rate and the goodput are subsequently discussed. Based onthe analysis, we perform numerical evaluations and reveal thatsynchronous PNC can outperform conventional network coding(CNC) even when taking synchronization errors and overheadinto account. The theoretical throughput gain of PNC over CNCcan be approached when using the MLE based synchronizationmethod with optimized training sequence length. The resultsin this paper provide some insights and benchmarks for theimplementation of synchronous PNC.

    Index Terms—Communication; denoise-and-forward (DNF);physical-layer network coding (PNC); synchronization; two-wayrelay networks.

    I. INTRODUCTION

    Relay-aided communications are widely adopted when di-rect communications among end nodes cannot be performed.Physical-layer network coding (PNC) [2] is considered as apromising technology to improve the throughput performanceof relay networks. It employs the natural network codingability introduced by the superposition of electromagneticwaves. Between the two methods of PNC, i.e. amplify-and-forward [3] and denoise-and-forward (DNF), the DNF methodshows more performance advantages because it avoids noiseamplification [4]. Hence, DNF has attracted much interest in

    Some preliminary ideas of this paper have been presented in IEEEGLOBECOM 2012 [1].

    This work was supported in part by the National Natural Science Foundationof China (61172051), the Fok Ying Tung Education Foundation (121065),the Fundamental Research Funds for the Central Universities (N110204001,N110804003, N120804002, N120404001), the Program for New CenturyExcellent Talents in University (NCET-12-0102), and the Specialized ResearchFund for the Doctoral Program of Higher Education (20110042110023,20110042120035, 20120042120049).

    Y. Huang, Q. Song (corresponding author) and L. Guo are withSchool of Information Science and Engineering, Northeastern University,Shenyang 110819, P. R. China. Email: [email protected], {songqingyang,guolei}@ise.neu.edu.cn.

    S. Wang is with Department of Electrical and Electronic Engineer-ing, Imperial College London, SW7 2AZ, United Kingdom. Email:[email protected].

    A. Jamalipour is with School of Electrical and Information Engineering,University of Sydney, NSW, 2006, Australia. Email: [email protected].

    recent research, and this paper considers the DNF scheme. Weuse DNF and PNC interchangeably in subsequent discussions.

    Research regarding PNC has been carried out focusingon two aspects: phase-asynchronous PNC [5]–[7] and phase-synchronous PNC [8]–[12]. The basic idea of asynchronousPNC is to map the superposed signal with arbitrary phasedifferences to encoded symbols. However, these schemesrequire knowledge of the instantaneous phases of the signalssuperposing at the relay, and imperfect channel informationmay also degrade the performance of asynchronous PNC[13]. Hence, tracking phase variations1 during data packettransmission is necessary (although synchronization is notneeded), which can be difficult especially for the superposedsignal. The complexity of obtaining symbol mapping undervarious phase differences can also be high, in particularwith high-level modulations [5]. Compared with asynchronousPNC, synchronous PNC allows more efficient constellationdesign [8] and can make use of capacity-approaching channelcodes [11]. The capacity region of the Gaussian two-wayrelay channel can also be reached with synchronous PNC[12]. Further, [14] shows that compared with other schemes,phase synchronization can maximize the minimum distancebetween adjacent points in the constellation for superposedsignals, provided that the signals are with the same modulationand amplitude at the relay. Therefore, we focus on the phase-synchronous PNC in this paper.

    Synchronization is a significant issue for synchronous PNC,which, however, has not been adequately studied. Existingworks [8]–[12] generally assume that superposing signalsarrive in-phase at the relay. However, these works havenot addressed how to achieve such a synchronization, tothe best of our knowledge. Although [15] investigated theimpact of imperfect synchronization for binary phase-shiftkeying (BPSK) modulated PNC, it did not explicitly introducea synchronization scheme and also did not investigate theinteraction between synchronization and data transmission.In the literature, some phase synchronization schemes fordistributed beamforming have been studied [16]–[19]. Al-though both PNC and distributed beamforming make use ofsignal superposition, the goal of PNC is to increase networkthroughput, while distributed beamforming is for increasingthe signal strength at the receiver. Meanwhile, the end nodescannot communicate with each other when using PNC (oth-erwise relaying is unnecessary); while in beamforming, the

    1Note that frequency errors accumulate over time and may cause the phasedifference between the two superposed signals change continuously.

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.

  • 2

    MA

    BC BC

    MA

    A R B

    Fig. 1. PNC over a two-way relay network.

    end nodes (sensors) may communicate with each other. Thedifference between these two techniques can make synchro-nization schemes for beamforming infeasible for PNC. Thelimited feedback-based synchronization for beamforming suchas [16] may cause large synchronization overhead due to theiterative process, which violates the intention of PNC, since theoverhead can reduce the goodput (i.e. effective throughput).Open-loop schemes as in [17]–[19] are also inapplicable forPNC, because they require the end nodes to communicatewith each other. Moreover, synchronization schemes for PNCdo not need to consider large-sized networks, because onlytwo nodes (rather than multiple nodes as in beamforming)are generally involved in the PNC process [20]. To considerthe requirements of PNC, in this paper, we propose a phasesynchronization scheme for PNC.

    Based on the proposed synchronization scheme, we discusssynchronization errors arising during the phase synchroniza-tion process and their impacts on the symbol error rate(SER) and network goodput in this paper. In terms of SERanalysis for PNC, [21] derived the SER for PNC with perfectsynchronization and unequal power of the superposing signals.Assuming the knowledge of channel gains, the SER for PNCwith decoding methods that do not require phase synchroniza-tion are discussed in [22] and [23], which respectively focus onminimum distance estimation and maximum a posteriori baseddecoding methods. The above existing works did not considerphase variations that may result from synchronization errors.In our preliminary work [24], we focused on SER of PNCwith deterministic phase deviations. In this paper, we focuson random phase deviations due to random synchronizationerrors. We derive analytical expressions of the average SERfor PNC with M -ary quadrature amplitude modulation (M -QAM), and subsequently study the impact of synchronizationerrors and overhead to the network goodput.

    We consider a two-way relay network as shown in Fig. 1.The main contribution of this paper is outlined as follows:

    1) We propose a phase synchronization scheme for PNC,which takes into account the characteristics and require-ments of PNC as aforementioned. The synchronizationerrors of the proposed synchronization scheme are thenanalyzed by considering potential frequency and phaseestimation techniques, namely, analog phase-locked loop(PLL), which is a conventional approach, and maximumlikelihood estimation (MLE), which is a more sophisti-cated but accurate approach.

    2) We derive analytical expressions and their approximatesolutions of the average SER for M -QAM modulatedPNC under the presence of synchronization errors. Ran-dom synchronization errors which accumulate and varyover time are considered. The analytical results are thenverified via simulations.

    3) We consider the joint operation of synchronization and

    data transmission, and study the goodput of the two-way relay network. The feasibility of phase-synchronousPNC is shown by numerical results.

    In summary, we present a phase synchronization scheme forPNC and study the interactions between the synchronizationoverhead, accuracy, SER, and network goodput, under estima-tion methods with PLL and MLE. Such a study enables usto understand whether phase-synchronous PNC is feasible orbeneficial when incorporating with the synchronization pro-cedure that uses common estimation methods. The analyticalresults also allow us to optimize the length of the trainingsequence that is used for synchronization (as will be discussedin Section V-B). Meanwhile, the framework that we use foranalysis can be applied when other estimation methods and/ornoise sources are considered.

    The remainder of this paper is organized as follows. SectionII illustrates the system model of this paper. Section III intro-duces the phase-level synchronization scheme and analyzeserrors with different estimation methods. In Section IV, theaverage SER under the impact of synchronization errors isdiscussed. The goodput of synchronous PNC is analyzed inSection V. Conclusions are drawn in Section VI.

    II. SYSTEM MODEL

    We consider a typical bidirectional relay network with flatfading channels, and the relay node R performs DNF relaying,as shown in Fig. 2. The DNF process includes multipleaccess (MA) phase and broadcast (BC) phase. Without lossof generality, we focus on square M -QAM modulated PNCin this paper, and end nodes A and B simultaneously transmitsquare M -QAM modulated data to the relay in the MA phase.The case of some common non-square M -QAM modulations(such as 32-QAM) can be treated similarly as square M -QAM,as discussed in [10]. The signal YR received by R is given by

    YR = SA + SB + Zn,R , (1)

    where SA and SB denote M -QAM signals from A and Brespectively, and Zn,R is the additive white Gaussian noise(AWGN) at R. In this paper, we consider the case where theaverage powers of SA and SB are equal.

    The minimum distance estimation is employed at the relayR to map the superposed signal YR to a network-codedsymbol. In this paper, PNC is performed with phase-levelsynchronization to maximize Euclidean distances, i.e. eachconstellation point (ideally) appears in the center of thecorresponding decision region.

    Because a M -QAM signal can be viewed as a complex√M -ary pulse amplitude modulation (

    √M -PAM) signal,

    its in-phase component IR(mΣ) and quadrature componentQR(nΣ) can be extracted from the superposed constellationpoint SmΣ,nΣ , i.e. SmΣ,nΣ = IR(mΣ) + jQR(nΣ). Thescalar values of these components are given by IR(mΣ) =2(mΣ −

    √M)d0 and QR(nΣ) = 2(nΣ −

    √M)d0, where

    mΣ, nΣ ∈ {1, 2, · · · , 2√M − 1} denote indices of constel-

    lation points for the superposed signal, and d0 representsthe Euclidean distance between two adjacent points in the

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.

  • 3

    A R B

    A R B

    A R B

    A R B

    A R B

    A R B

    Slot 1

    Slot 2

    Slot 3

    Slot 4MA

    BC

    Sync. Time

    (Tsync)

    Trans. Time

    (Ttrans)

    (Tsync and Ttransalter periodically.)

    Fig. 2. Network topology and timing diagram. Synchronization (sync.) andtransmission (trans.) alternate over the multiple access (MA) phase, and theyare performed periodically.

    constellation diagram for√M -PAM. The minimum distance

    estimation for ŜmΣ,nΣ is given as

    (m̂Σ, n̂Σ) = arg minmΣ,nΣ

    |YR − (IR(mΣ) + jQR(nΣ))|, (2)

    ŜmΣ,nΣ = IR(m̂Σ) + jQR(n̂Σ), (3)

    where | · | stands for the modulus (absolute value).The estimated ŜmΣ,nΣ will be mapped into a network-coded

    symbol with the approach proposed in [10].

    III. PHASE SYNCHRONIZATIONThis section firstly introduces a round-trip estimation based

    carrier synchronization scheme for PNC. Afterwards, we ana-lyze phase synchronization errors when performing phase andfrequency estimation using PLL and MLE, respectively.

    A. Synchronization Process

    As depicted in Fig. 2, the synchronization phase (whoselength is denoted as the synchronization time Tsync) is dividedinto four timeslots.

    In timeslot 1, the relay R broadcasts a beacon b0(t) =a0 cos(ωct + ϕ0), where a0 represents the amplitude of thissinusoidal signal, ωc denotes the reference angular frequency,and ϕ0 is the initial phase at t = 0. The received beaconbR,A(t) at end node A (because the case for node B is similar,we only focus on node A in the subsequent discussions) isgiven by

    bR,A(t) = aR,A cos(ωct+ ϕR,A) + Zn,A , (4)

    where aR,A and ϕR,A respectively denote the amplitude andphase of the received signal, and Zn,A denotes the AWGN atnode A. Upon receiving bR,A(t), node A estimates the value ofωc as ω̂c. Then, node A adjusts its local oscillator to generate asinusoidal signal with frequency ω̂c. The same estimation andrecovering process is performed at node B. By this means, weachieve frequency synchronization between nodes A and B;and the remaining timeslots are for phase synchronization.

    In timeslot 2, the recovered beacon at A is bounced backto the relay R. The signal that is received by node R is givenby

    bA,R(t) = aA,R cos(ω̂ct+ ϕA,R) + Zn,R , (5)

    where aA,R, ϕA,R, and Zn,R respectively denote the ampli-tude, phase, and AWGN at node R. The relay R estimates thephase ϕA,R of the received signal, and the estimation resultis denoted by ϕ̂A,R. The process is similar for node B intimeslot 3.

    In timeslot 4, the relay R transmits the difference betweenthe estimated phase ϕ̂A,R and a reference phase ϕref back tothe end node A for compensation. The reference phase ϕref canbe set to an arbitrary value (for instance ϕ0), because we onlyrequire that the signals arrive in-phase at R. The operation fornode B is same as the above. After compensation, the signalsfrom nodes A and B arrive in-phase (both aligned to ϕref) atthe relay R.

    In the transmission phase (whose length is denoted as thetransmitting time Ttrans) that follows, the recovered signal isused as the carrier signal. Unfortunately, the frequency estima-tion error causes the phase error increase with time. Therefore,as shown in Fig. 2, synchronization needs to be performedperiodically over the MA phase. The synchronization periodalso needs to be within the duration that channel state remainsalmost unchanged.

    B. Synchronization Errors

    Estimation errors occur during synchronization, becausereceived beacons are interfered with AWGNs as in (4) and (5).Thus, ωc = ∆ωc+ω̂c and ϕA,R = ∆ϕA,R+ϕ̂A,R, where ∆ωcand ∆ϕA,R represent corresponding error terms. The error∆ωc occurs at the end node, and the error ∆ϕA,R occurs atthe relay, as discussed in Section III-A. The frequency error∆ωc also results in a linearly increasing phase error duringdata transmission, which makes the phases of the two signalsmisalign at the relay and hence increases the average SER.

    The errors vary with different estimation methods. In thesubsequent discussion, we focus on error analysis for esti-mation with PLL and MLE, respectively. Note that, althoughfrequency and phase estimation are respectively (not concur-rently) performed at the end nodes and the relay, we analyzeboth frequency and phase errors in the subsequent discussion.The reason is that PLL and MLE can estimate both frequencyand phase. Meanwhile, in a general network, each node mayhave both roles of end node and relay [25]. The specific roledepends on the traffic pattern of the network. In such cases, theestimation module can be reused for estimating the frequencyand phase. When necessary, we use subscripts “PLL” and“MLE” to represent variables in the corresponding cases.

    C. Synchronization Error with PLL Based Estimation

    In this subsection, we consider the scenario that a PLL isadopted in the nodes to track the frequency and phase. Wederive analytical expressions of the variances of estimationerrors through the transfer function of a linearized PLL model.As depicted in Fig. 3, the PLL model consists of a phasedetector (PD), a loop filter, and a voltage-controlled oscillator(VCO). The phase of the input (in the S-domain) is denoted byϕin(s) and the phase of the VCO output is denoted by ϕout(s);Kd and K0 respectively denote the phase-detector gain and theVCO gain; HLF(s) is the transfer function of the loop filter. In

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.

  • 4

    d LF 0

    n, rcv

    out

    in

    VCO

    n, PLL

    Fig. 3. Linearized PLL model.

    timeslot 1, the PLL works in the closed-loop mode, to track thephase and frequency of the reference carrier sent by the relayR. In the remaining timeslots of the synchronization processand also during data transmission, the oscillating frequencyωVCO(s) of the VCO is captured by a sample and hold circuit,and the PLL operates in the open-loop mode without furthertracking the input signal. The output of the VCO is then usedto modulate the data symbols for transmission. Note that thephase difference for compensation can be obtained with anadditional phase detector with ϕout and ϕref as the input; phasecompensation (as discussed in Section III-A) can be performedon the baseband, i.e. by rotating the signal constellation.

    Considering an input such as (4), as discussed in [26], theadditive noise term Zn,rcv at the receiver can be equivalent toZ ′n,rcv as shown in Fig. 3. The power spectral density (PSD)of Z ′n,rcv is 2N0/a

    2rcv = N0Ts/Es, where N0 denotes the

    PSD of Zn,A, arcv stands for the received signal amplitude atthe receiver, Ts denotes the symbol duration, and Es denotesthe energy per symbol. Meanwhile, Z ′n,rcv is a narrow bandnoise signal with bandwidth ωB , because the received signalis processed by a bandpass filter at the receiver. For an idealreceiver that maximizes the bandwidth efficiency, we haveωB = 2π/(2Ts) for one dimensional signal. The additionalnoise from the components inside the PLL is denoted byZn,PLL, which can be conservatively regarded as AWGN withPSD Np [17].

    The value of ωVCO that is captured by the sample andhold component corresponds to the estimated carrier frequencyω̂c,PLL. Hence, to investigate the error ∆ωc,PLL of frequencyestimation, we need to study the noise component at ωVCO. Wenote that the noise components Z ′n,rcv and Zn,PLL can also beregarded as the input of the PLL, as shown in Fig. 3. Therefore,the transfer function for noise signal can be evaluated by

    H(s) =ωVCO(s)

    ϕin(s)=

    sKdK0HLF(s)

    s+KdK0HLF(s). (6)

    Considering the respective PSD and bandwidth of Z ′n,rcvand Zn,PLL, we can obtain the variance of the frequency error∆ωc,PLL:

    σ2ωc,PLL =2N0a2rcv

    · 12π

    ∫ ωB0

    |H(jω)|2 dω

    +Np ·1

    ∫ ∞0

    |H(jω)|2 dω , (7)

    where H(jω) denotes the system frequency response. BecauseZ ′n,rcv and Zn,PLL are Gaussian noises, ∆ωc,PLL conforms toa zero-mean Gaussian distribution given by N (0, σ2c,PLL).

    In the case of a second-order PLL with lag filter (which isfrequently used in a wireless repeater [26], for instance), H(s)can be rewritten as

    H(s) =ω2ns

    s2 + 2ξωns+ ω2n, (8)

    where ωn and ξ respectively denote the natural frequency anddamping ratio. Then, the integral terms in (7) can be evaluated2

    as follows:∫ ωB0

    |H(jω)|2 dω = ω3n

    4ξ(f1 + f2 − f3) , (9)

    where

    f1 =arctan

    (ωB + ωn

    √1− ξ2

    ξωn

    ),

    f2 =arctan

    (ωB − ωn

    √1− ξ2

    ξωn

    ),

    f3 =ξ

    2√1− ξ2

    ln

    (ω2B + 2ωBωn

    √1− ξ2 + ω2n

    ω2B − 2ωBωn√

    1− ξ2 + ω2n

    );

    and ∫ ∞0

    |H(jω)|2 dω = πω3n

    4ξ. (10)

    The phase error can be derived in a similar method byevaluating the transfer function between ϕout(s) and ϕin(s).For the second-order PLL with lag filter, this transfer functionis

    H ′(s) =ϕout(s)

    ϕin(s)=

    KdK0HLF(s)

    s+KdK0HLF(s)=

    ω2ns2 + 2ξωns+ ω2n

    .

    (11)The variance σ2ϕPLL of the phase error can be evaluated in thesame way as (7), with∫ ωB

    0

    |H ′(jω)|2 dω = ωn4ξ

    (f1 + f2 + f3) , (12)

    and ∫ ∞0

    |H ′(jω)|2 dω = πωn4ξ

    . (13)

    The natural frequency ωn is related to the necessary trainingtime Ttrain, which is the duration that the PLL spends on adjust-ing frequencies, also known as the settling time of PLL. Fora second-order PLL with lag filter, we have ωn ≈ 4/(ξTtrain)[28]. Because estimation needs to be performed in timeslots1, 2, and 3, we have Tsync = 3Ttrain +Tctrl, where Tctrl denotesthe duration of control data transmission in timeslot 4.

    D. Synchronization Error with MLE

    In this subsection, we consider the case where nodesestimate the frequency and phase with the MLE method.Although more sophisticated maximum a posteriori (MAP)based algorithms such as in [29] have been proposed, thispart analyzes estimation errors based on the MLE algorithmproposed in [30] which is believed to be more feasible andrelaxes the need of huge computational complexity [31], due

    2We employ Maple [27] to evaluate some sophisticated integrals.

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.

  • 5

    to practical considerations. Different from [30], we considerarbitrary symbol duration (Ts) in our discussion, to betterrelate the analysis to actual data transmission.

    For the received beacon brcv(t), when the symbol timingis accurate [32], [33], putting the signal into a pair of or-thogonal matched filters and sampling the resulting signal ata time interval of Ts yields a complex signal b̃rcv(kTs) (k =1, 2, 3, ...), where

    ∣∣b̃rcv(kTs)∣∣2 and arg (b̃rcv(kTs)) respectivelycorrespond to the energy and average phase of brcv(t) over Ts.Then, likelihood function can be written as

    L(ωc, ϕ)

    =

    (1

    πN0

    )Ntrainexp

    (−

    Ntrain−1∑k=0

    ∣∣b̃rcv(kTs)−√Esej(ωckTs+ϕ)∣∣2N0

    ),

    (14)

    where N0 is the variance of AWGN after traversing thematched filter and Ntrain denotes the length of the trainingsequence. Similarly with [30], by solving

    ∂ lnL(ωc, ϕ)∂ωc

    = 0 and∂ lnL(ωc, ϕ)

    ∂ϕ= 0 , (15)

    we obtain the maximum-likelihood estimators for ωc and ϕ as

    ω̂c,MLE =

    Ntrain−1∑k=0

    kUVNtrain−1∑k=0

    U −Ntrain−1∑k=0

    UVNtrain−1∑k=0

    kU

    TsNtrain−1∑k=0

    k2UNtrain−1∑k=0

    U − Ts(

    Ntrain−1∑k=0

    U

    )2(16)

    and

    ϕ̂MLE=

    Ntrain−1∑k=0

    kUVNtrain−1∑k=0

    kU −Ntrain−1∑k=0

    k2UNtrain−1∑k=0

    UV(Ntrain−1∑k=0

    kU

    )2−

    Ntrain−1∑k=0

    k2UNtrain−1∑k=0

    U

    ,

    (17)where U = |brcv(kTs)| and V = arg

    (brcv(kTs)

    ).

    The variances of estimation errors are bounded by theCramér-Rao lower bounds by

    σ2ωc,MLE ≥6N0

    EsNtrain(N2train − 1)T 2s(18)

    and

    σ2ϕMLE ≥N0(2Ntrain − 1)

    EsNtrain(Ntrain + 1). (19)

    The lower bounds in (18) and (19) can be attained when theSNR is relatively high, as discussed in [30]. Hence, we usethese values as to approximate the variances when using MLEin subsequent discussions.

    For MLE, we have Ttrain = NtrainTs and Ttrans = NtransTs,where Ntrans denotes the number of transmitted symbolsover the transmitting time. Similar to the case of PLL,Tsync = 3Ttrain + Tctrl, ∆ωc,MLE and ∆ϕMLE conform to zero-mean Gaussian distributions respectively given by ∆ωc,MLE ∼N (0, σ2ωc,MLE) and ∆ϕMLE ∼ N (0, σ

    2ϕMLE

    ). The impacts ofthese errors will be analyzed in subsequent sections.

    IV. SYMBOL ERROR RATE WITH ESTIMATION ERRORS

    This section analyzes the SER at the relay under the impactof estimation errors studied in the previous section. We firststudy the SER for M -QAM and quadrature phase shift keying(QPSK) with arbitrary deterministic phase deviations. Then,analytical expression of the average SER over a period of timewith random phase deviations is derived. Because a receiverusually performs channel estimation through preambles [34],we assume that the receiver only tracks the phase fromknowledge of the preamble at the beginning of each dataframe. The receiver is unaware of subsequent phase variationscaused by frequency offsets (i.e. ∆ωc) in data carrying signalsover the transmitting time [35].

    A. SER with Deterministic Phase Deviations

    To ensure unique decodability for PNC with M -QAM,points in any

    √M by

    √M square in the constellation for

    superposed signals have to be mapped into different symbols[10]. When M is large enough, it is of low probability thatthe noise can let the superposed signal step over severaldecision regions and reach a region that should be mappedto a coded symbol that is identical with the correct symbol.Accordingly, we neglect the correct probability of this case inour discussion.

    When power control and synchronization are performed, theminimum distance estimation in the 2-dimensional space canbe separately performed in the in-phase channel (I-channel)and the quadrature channel (Q-channel). Assume that thetransmitted symbols are equiprobable, the error probabilitiescalculated in both I-channel and Q-channel are equal. Fordifferent intervals of decision regions, the error probabilitiesin the I-channel can be approximated by [24]:

    Ps

    ∣∣∣mA,mB

    Q

    (d0 + µ0 − µ

    σ0

    ), if mA,mB = 1

    Q

    (d0 + µ− µ0

    σ0

    ), if mA,mB =

    √M

    Q

    (d0 + µ− µ0

    σ0

    )+Q

    (d0 + µ0 − µ

    σ0

    ), else

    (20)

    where mA, nA,mB , nB ∈ {1, 2, · · · ,√M} respectively rep-

    resent indices of the M -QAM constellation points in the I-channel and Q-channel from nodes A and B; σ0 =

    √N0/2

    denotes the standard deviation of AWGN in the I-channel; µ0denotes the original constellation point without phase deviationin the I-channel and it is given by µ0 = 2(mA + mB −1 −

    √M)d0; and µ denotes the constellation point when

    suffering phase deviation in the I-channel, which is given byµ = (2mA−1−

    √M)d0 cosψA+(2mB−1−

    √M)d0 cosψB−

    (2nA − 1 −√M)d0 sinψA − (2nB − 1 −

    √M)d0 sinψB .

    Variables ψA and ψB represent instantaneous phase deviations(with respect to strict synchronization when the deviations arezero) of SA and SB . According to [36], d0 can be obtainedby

    d0 =

    (3Eb log2

    √M

    M − 1

    )1/2, (21)

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.

  • 6

    where Eb represents the average energy per bit of the receivedsignal at the relay R. For equiprobable symbols, any combina-tion of (mA, nA,mB , nB) shares the same probability 1/M2.Hence the error probability in the I-channel is

    Ps

    ∣∣∣I-channel

    =1

    M2

    √M∑

    n′A,n′B=1

    (Ps

    ∣∣∣mA,mB=1

    + Ps

    ∣∣∣mA,mB=

    √M

    +

    √M∑

    mA,mB=1

    Ps

    ∣∣∣mA+mB ̸=2,2

    √M

    ). (22)

    Then, the SER for M -QAM modulated PNC with determin-istic phase deviations can be evaluated by

    Ps = 1−(1− Ps

    ∣∣∣I−channel

    )2. (23)

    When using QPSK, the approximated results (which neglectconstellation points that are mapped to identical symbols) canbecome inaccurate, because there is only one other decisionregion between those regions that are to be mapped to the samesymbol. Therefore, we evaluate the exact SER for QPSK. Thein-phase component of the superposed constellation is givenby IR(mΣ) ∈ {−2d0, 0, 2d0}, and the mapping rule is that{−2d0, 2d0} is mapped to bit “0” (or, correspondingly, “1”)and {0} is mapped to bit “1” (or, correspondingly, “0”). Thus,cases of mA,mB = 1 and mA,mB =

    √M in (20) can be

    combined as

    P ′s

    ∣∣∣mA,mB=1 or

    √M

    = Q

    (µ− d0σ0

    )−Q

    (µ+ d0σ0

    ). (24)

    Let (24) be the substitutes for cases of mA,mB = 1 andmA,mB =

    √M in (20), the exact SER for QPSK modulated

    PNC with phase deviation can be calculated with (23).

    B. Average SER with Random Phase over A Segment of TimeThe phase deviation accumulates with time due to the

    presence of frequency estimation error. Because the transmit-ting time is usually much longer than the duration of thetraining sequence, the phase deviation can accumulate to avalue which is much larger than the initial phase estimationerror. Therefore, we mainly focus on phase deviation causedby frequency error in this subsection.

    Remark that in the following analysis, we only focus onψA due to the similarity between ψA and ψB . As depicted inFig. 4, the instantaneous phase deviation ψA(t) is given byψA(t) = tψA,max/Ttrans, where ψA,max denotes the maximumphase deviation at the end of each data transmission. The phasedeviation process is a cyclostationary process with Tsync +Ttrans as the period. Due to the relationship given by ψA,max =∆ωcTtrans, both ψA,max and ψA(t) follow zero-mean Gaussiandistributions. The variance of ψA,max is denoted by σ2A,max, andσ2A,max = T

    2transσ

    2ωc . It follows that the instantaneous variance

    of ψA(t) is

    σ2A(t) =t2

    T 2transσ2A,max. (25)

    Then, the expectation of the SER at time instant t is

    Ps(t)=

    ∫ +∞−∞

    ∫ +∞−∞

    Ps(ψA, ψB)p(ψA, ψB , t) dψAdψB , (26)

    sync trans

    A,max

    A

    Fig. 4. Phase deviation at end node A. Phase deviation increases linearlydue to the frequency estimation error that is generated during synchronization.Random frequency errors cause different deviations in different transmissionphases. Similar phenomenon can be observed at end node B.

    where Ps is calculated with (23) under different values of ψAand ψB , p(ψA, ψB, t) stands for the joint probability densityfunction of ψA(t) and ψB(t). Because ψA(t) and ψB(t) areindependently distributed, and ψA(t) ∼ N (0, σ2A(t)), ψB(t) ∼N (0, σ2B(t)), we have

    p(ψA, ψB , t) =1

    2πσA(t)σB(t)e− ψ

    2A

    2σ2A

    (t)− ψ

    2B

    2σ2B

    (t) (27)

    The average SER for over the whole transmitting timeduring the MA phase is then given by

    Ps,MA =1

    Ttrans

    ∫ Ttrans0

    Ps(t) dt. (28)

    C. Approximate Analytical Solution

    Due to the absence of explicit expressions for (26) and (28)and the complexity when calculating numerical integrations,in this subsection, we derive an approximate solution to (26)and (28).

    Assume that the instantaneous phase deviations are small,i.e. ψ(t) ≈ 0, we have sin(ψ(t)) ≈ ψ(t) and cos(ψ(t)) ≈ 1.Substituting these approximations into (20), and recalling thatψA(t) and ψB(t) are Gaussian random variables, ∆µ(t) =µ(t) − µ0(t) can be regarded as a Gaussian random variablewith zero mean and variance σ2µ(nA, nB , t) = (2nA − 1 −√M)2d20σ

    2A(t) + (2nB − 1−

    √M)2d20σ

    2B(t). Further, by ig-

    noring the square terms in (23), we achieve Ps ≈ 2Ps∣∣I−channel,

    and the integral3 in (26) can be performed on each termcorresponding to one Q-Function in (20). Considering that

    Q(x) ≈ Qapprox(x) ,{

    12e

    − x22 if x ≥ 01− 12e

    − x22 if x < 0, (29)

    3Note that the integral can be written as a one dimensional integral now,because we consider a single Gaussian variable ∆µ here.

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.

  • 7

    the integrated value for each term in (20) is4

    F (t) =1√

    2πσµ(t)

    ∫ +∞−∞

    Qapprox

    (d0 ±∆µ

    σ0

    )e− (∆µ)

    2

    2σ2µ(t) d(∆µ)

    =

    σ0erf

    (d0σ0

    σµ(t)√

    2σ20+2σ2µ(t)

    )e− d

    20

    2σ20+2σ2µ(t)

    2√σ20 + σ

    2µ(t)

    +Q

    (d0σµ(t)

    ).

    (30)

    Summing up the result in (30) for all the indicesmA, nA,mB , nB ∈ {1, 2, ...,

    √M} as in (22) and multiplying

    by two yields the approximate result for (26).To obtain an approximate result for (28), we perform an

    asymptotic analysis. When σ0 → 0, the first term in (30)vanishes to zero. Again, using Q(x) ≈ 12e

    − x22 for x ≥ 0and σµ(t) = tσµ,max/Ttrans, where σµ,max denotes the standarddeviation of ∆µ at the end of each data transmission, we have

    F (t)∣∣∣σ0→0

    ≈ 12e− d

    20T

    2trans

    2t2σ2µ,max . (31)

    Taking its logarithm yields

    ln (2F (t))∣∣∣σ0→0

    ≈ − d20T

    2trans

    2t2σ2µ,max. (32)

    It follows that

    ln (2F (t))

    ln (2F (Ttrans))

    ∣∣∣∣σ0→0

    ≈ T2trans

    t2, (33)

    and

    F (t)|σ0→0≈(2F (Ttrans))

    T2transt2

    2

    ∣∣∣∣∣∣σ0→0

    . (34)

    Relaxing the constraint of σ0 → 0, the average value of oneterm in (20) can be approximated by

    F ≈ 12Ttrans

    ∫ Ttrans0

    (2F (Ttrans))T2transt2 dt

    = F (Ttrans)−√

    −π · Fln Q(√

    −2Fln), (35)

    where Fln = ln(2F (Ttrans)). The approximate result for (28)can then be evaluated by summing up the result in (35) for allthe indices mA, nA,mB , nB ∈ {1, 2, ...,

    √M} as in (22) and

    multiplying by two.

    D. Numerical Results

    We perform Monte Carlo simulations to verify the analyticalresults. Figs. 5 and 6 show the comparison among simulationresults, analytical results evaluated by (28) using numerical in-tegration, and approximate analytical results derived in SectionIV-C. We consider the case where σ2A,max = σ

    2B,max = σ

    2max.

    The results indicate agreements between analytical results,approximate analytical results, and simulation results. It can beobserved that with both 16-QAM and QPSK, the average SERcurves do not always fall as SNR increases, but level off andconverge to stable values at some values of σmax. The reason is

    4For simplicity, we omit the variables nA and nB .

    0 5 10 15 20 25 3010

    −6

    10−5

    10−4

    10−3

    10−2

    10−1

    100

    Eb / N

    0 (dB)

    Avera

    ge S

    ER

    Analytical Results

    Approx. Ana. Results

    Simulated Results

    0 5 10 15 20 25 3010

    −6

    10−5

    10−4

    10−3

    10−2

    10−1

    100

    Eb / N

    0 (dB)

    Avera

    ge S

    ER

    σmax

    = 2°

    σmax

    = 3°

    σmax

    = 5°

    σmax

    = 4°

    σmax

    = 8°

    Fig. 5. Average SER for 16-QAM modulated PNC.

    0 5 10 15 20 25 3010

    −6

    10−5

    10−4

    10−3

    10−2

    10−1

    100

    Eb / N

    0 (dB)

    Avera

    ge S

    ER

    Analytical Results

    Approx. Ana. Results

    Simulated Results

    0 5 10 15 20 25 3010

    −6

    10−5

    10−4

    10−3

    10−2

    10−1

    100

    Eb / N

    0 (dB)

    Avera

    ge S

    ER

    σmax

    = 7°

    σmax

    = 5° σmax = 10°

    σmax

    = 15°

    σmax

    = 20°

    Fig. 6. Average SER for QPSK modulated PNC.

    that in high SNR regions, the symbol error is mainly caused byphase deviations, therefore the SER does not decrease muchwith increasing SNR as long as the value of σmax remainsunchanged. The approximate analytical results are not veryaccurate when σmax is large, as shown in Fig. 6, because theassumption sinψ ≈ ψ only holds for small phase deviations.

    V. GOODPUT PERFORMANCE ANALYSIS

    This section investigates the goodput (i.e. the amount ofsuccessfully transmitted information) performance for PNCunder the joint impact from synchronization overhead andincreased SER caused by phase deviations.

    A. Analytical Evaluation

    Recall that in a two-way relay network with bidirectionalflows, the ideal throughput for conventional network coding(CNC) and PNC are respectively 2/3 log2M and log2M [2].Considering synchronization overhead and packet loss, the

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.

  • 8

    goodput for PNC is

    GPNC =1

    2(1− ρ)(1− Ps,MA)Npk(1− Ps,BC,A)Npk log2M

    +1

    2(1− ρ)(1− Ps,MA)Npk(1− Ps,BC,B)Npk log2M, (36)

    where ρ denotes the synchronization overhead in percentages.When the transmission rate remains unchanged, the transmit-ting time in the MA phase equals that in the BC phase. Then,we have

    ρ = Tsync/(Tsync + 2Ttrans). (37)

    The value of (1−Ps,MA)Npk is the packet success rate at R overthe MA phase, where Npk denotes the packet length and Ps,MAis the average SER during MA phase which can be evaluatedby (28). Likewise, (1 − Ps,BC,A)Npk and (1 − Ps,BC,B)Npk arerespectively the packet success rates at nodes A and B in theBC phase, where Ps,BC,A and Ps,BC,B respectively representSERs for common M -QAM at nodes A and B. The SER forM -QAM is given by [37]:

    PM−QAM = 4

    (√M − 1√M

    )Q

    (√3Es

    N0(M − 1)

    )

    −4

    (√M − 1√M

    )2Q2

    (√3Es

    N0(M − 1)

    ).(38)

    Likewise, the goodput for CNC is

    GCNC =1

    2· 2 log2M

    3(1− Ps,B,R)Npk(1− Ps,R,A)Npk

    +1

    2· 2 log2M

    3(1− Ps,A,R)Npk(1− Ps,R,B)Npk , (39)

    where Ps,B,R, Ps,A,R, Ps,R,A, and Ps,R,B respectively denotethe SERs for corresponding uplinks (B → R and A→ R) anddownlinks (R → A and R → B), and these probabilities canalso be evaluated by (38).

    B. Impact of Training Sequence Time-Length

    The training sequence time-length Ttrain has a trade-offeffect on the goodput when using PNC. Recall that Tsync =3Ttrain + Tctrl as discussed in Section III-B, a larger value ofTtrain yields longer synchronization time, which may increasethe overhead. However, a larger Ttrain also results in a moreprecise phase and frequency estimation, which could increasethe packet success rate and, subsequently, the goodput.

    Therefore, an appropriate value of Ttrain should be selectedto maximize the goodput. This can be formulated as thefollowing optimization problem:

    maxTtrain

    GPNC

    s.t. 0 ≤ Ttrain ≤1

    3(Tc − Tctrl), (40)

    where Tc denotes the period of the synchronization cycle. Wesolve (40) using numerical evaluation methods in MATLAB.Fig. 7 shows the optimal Ttrain under different values of Eb/N0(i.e. SNR per bit), where Eb denotes the energy per bit,when using the MLE method and the approximate solutionas discussed in Section IV-C for evaluation.

    6 10 20 30 40 500

    2.5

    5

    7.5

    10

    12.5

    15

    17.5

    Eb / N

    0 (dB)

    Optim

    ized T

    train (

    ms)

    QPSK

    16−QAM

    Fig. 7. Optimized Ttrain values under different SNRs when using the MLEmethod.

    C. Numerical Results

    The goodput performance of synchronous PNC is evaluatednumerically in this subsection. We consider PNC with bothPLL and MLE based synchronization methods (notated asPLL-PNC and MLE-PNC in the following discussions), andalso compare with the goodput of CNC.

    In our simulations, we set Tc = 64 ms, which correspondsto the channel coherence time (i.e. the time that the channelalmost remains unchanged) of fixed nodes with 2.4 GHzradio transceivers in fast varying environments [35]. Thetransmitting time Ttrans is then Ttrans = Tc−Tsync. The symbolduration Ts is set to 1µs, and the packet length Npk is set to1024 bytes. The duration of control data Tctrl is set to 0.3 ms,which is enough for transmitting several hundred bits. For thePLL, the values of ξ and Np are respectively set to 0.707 and7.0× 10−11 Hz−1 [17], [26].

    Regarding the value of Ttrain, we consider both fixed andoptimal value settings. For the fixed value setting, we setTtrain = 5 ms and evaluate the performance of PLL-PNC andMLE-PNC, respectively. We select Ttrain = 5 ms because itis close to the optimal Ttrain corresponding to the minimumEb/N0 requirement for QPSK and 16-QAM, as shown in Fig.7. For the optimal value setting, we set Ttrain to the optimalvalues as in Fig. 7 and only evaluate the performance ofMLE-PNC. We do not evaluate the performance of PLL-PNCwith optimal Ttrain, because the settling time of the PLL is adesigned hardware parameter which is difficult to adjust basedon Eb/N0. However, when using the MLE based estimationmethod, it is possible to adapt the training sequence length toEb/N0.

    The goodputs when using different techniques are shownin Fig. 8. It can be observed that, when Ttrain = 5 ms, thegoodputs of MLE-PNC and PLL-PNC both converge to stablevalues that correspond to a goodput gain of approximately 1.30over CNC, for both 16-QAM and QPSK modulations. Such aconvergence is because, at high SNR values, the packet loss isvery low so that the goodput does not vary much with the SNR.The difference between the observed goodput gain and themaximal throughput gain (which is 1.5) is due to the overhead.With our simulation settings, according to (37), the overheadρ = 15.3/(15.3+2×48.7) = 0.136. The goodput gain with thegiven overhead can be evaluated by 1.5(1− ρ) = 1.30, whichmatches with the numerical results. At medium SNR values,

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.

  • 9

    Fig. 8. Comparison between the goodput of synchronous PNC and that ofCNC.

    6 10 20 30 40 500

    0.5

    1

    1.5

    2

    2.5

    3

    3.5

    4

    Eb / N

    0 (dB)

    Go

    od

    pu

    t (b

    ps/H

    z)

    Optim. QPSK MLE

    Optim. 16−QAM MLE

    Approx. Optim. QPSK MLE

    Approx. Optim. 16−QAM MLE

    Fig. 9. Actual and approximated goodput of MLE-PNC with optimized Ttrain.

    we can see that MLE-PNC outperforms PLL-PNC, becauseMLE provides higher estimation accuracy and the goodput isaffected by the packet loss in this SNR region.

    For both 16-QAM and QPSK, MLE-PNC with optimizedTtrain outperforms the other schemes, and its goodput keepsincreasing with Eb/N0. This is because the value of Ttrain isoptimized based on the SNR. At high SNRs, Ttrain can be con-siderably small, yielding a very small overhead. The goodputat some higher SNR values is plotted in Fig. 9. We can observethat, when Eb/N0 = 50 dB, the goodput gain is approximately1.48, which is very close to the maximal throughput gain.Also, the goodputs evaluated with the analytical approximatesolutions as discussed in IV-C matches with their actual values.

    VI. CONCLUSIONS

    In this paper, we have analyzed the feasibility of PNCwith phase-level synchronization. We have proposed a syn-chronization scheme for PNC. Subsequently, we have re-

    vealed analytical relationships among the goodput, averageSER, synchronization overhead, and estimation errors, whenusing either PLL or MLE based synchronization techniques.Numerical results show that the goodput of a two-way re-lay network can benefit from synchronous PNC, and MLEbased synchronization schemes can attain more goodput gainthan PLL based schemes. Our study also reveals that highergoodput can be obtained by adjusting the training sequencelength according to the SNR. The goodput evaluated in thispaper is based on symbols without channel coding. We wouldforesee that the goodput performance of synchronous PNCcould be further improved when channel coding is performed.Although the error analysis in this paper focuses on phaseand frequency estimation errors, it can be easily generalizedto incoporate some other error terms, using the same analyticalframework. The results in this paper provide some insights andbenchmarks for the implementation of synchronous PNC. Inthe future, we will focus on the impact of estimation errors onasynchronous PNC schemes, because asynchronous PNC alsorequires phase and frequency tracking (although adjustment isnot needed), which introduces estimation errors similarly assynchronous PNC.

    REFERENCES

    [1] Y. Huang, Q. Song, S. Wang, and A. Jamalipour, “Phase-level synchro-nization for physical-layer network coding,” in Proc. IEEE GLOBECOM2012, Dec. 2012, pp. 4423–4428.

    [2] S. Zhang, S. C. Liew, and P. P. Lam, “Hot topic: Physical-layer networkcoding,” in Proc. ACM MobiCom 2006, Sep. 2006, pp. 358–365.

    [3] P. Popovski and H. Yomo, “The anti-packets can increase the achievablethroughput of a wireless multi-hop network,” in Proc. IEEE ICC 2006,Jun. 2006, pp. 3885–3890.

    [4] K. Lee and L. Hanzo, “Resource-efficient wireless relaying protocols,”IEEE Wireless Commun. Mag., vol. 17, no. 2, pp. 66–72, Apr. 2010.

    [5] T. Koike-Akino, P. Popovski, and V. Tarokh, “Optimized constellationsfor two-way wireless relaying with physical network coding,” IEEE J.Sel. Areas Commun., vol. 27, no. 5, pp. 773–787, Jun. 2009.

    [6] L. Lu, S. C. Liew, and S. Zhang, “Optimal decoding algorithm forasynchronous physical-layer network coding,” in Proc. IEEE ICC 2011,Jun. 2011, pp. 1–6.

    [7] T. Yang and I. Collings, “Asymptotically optimal error-rate performanceof linear physical-layer network coding in rayleigh fading two-way relaychannels,” IEEE Commun. Lett., vol. 16, no. 7, pp. 1068 –1071, Jul.2012.

    [8] M. Noori and M. Ardakani, “On symbol mapping for binary physical-layer network coding with PSK modulation,” IEEE Trans. WirelessCommun., vol. 11, no. 1, pp. 21–26, Jan. 2012.

    [9] H. J. Yang, Y. Choi, and J. Chun, “Modified high-order PAMs for binarycoded physical-layer network coding,” IEEE Commun. Lett., vol. 14,no. 8, pp. 689 –691, Aug. 2010.

    [10] S. Wang, Q. Song, L. Guo, and A. Jamalipour, “Constellation mappingfor physical-layer network coding with M-QAM modulation,” in Proc.IEEE GLOBECOM 2012, Dec. 2012, pp. 4429–4434.

    [11] M. P. Wilson, K. Narayanan, H. D. Pfister, and A. Sprintson, “Jointphysical layer coding and network coding for bidirectional relaying,”IEEE Trans. Inf. Theory, vol. 56, no. 11, pp. 5641–5654, Nov. 2010.

    [12] W. Nam, S.-Y. Chung, and Y. H. Lee, “Capacity of the gaussian two-way relay channel to within 1/2 bit,” IEEE Trans. Inf. Theory, vol. 56,no. 11, pp. 5488–5494, Nov. 2010.

    [13] K. Yasami, A. Razi, and A. Abedi, “Analysis of channel estimation errorin physical layer network coding,” IEEE Commun. Lett., vol. 15, no. 10,pp. 1029–1031, Oct. 2011.

    [14] T. Koike-Akino, P. Popovski, and V. Tarokh, “Adaptive modulation andnetwork coding with optimized precoding in two-way relaying,” in Proc.IEEE GLOBECOM 2009, Nov. 2009, pp. 1–6.

    [15] S. Zhang, S.-C. Liew, and P. P. Lam, “On the synchronization ofphysical-layer network coding,” in Proc. IEEE Information TheoryWorkshop 2006, Oct. 2006, pp. 404–408.

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.

  • 10

    [16] S. Sigg, R. M. E. Masri, and M. Beigl, “Feedback-based closed-loopcarrier synchronization: A sharp asymptotic bound, an asymptoticallyoptimal approach, simulations, and experiments,” IEEE Trans. MobileComput., vol. 10, no. 11, pp. 1605–1617, Nov. 2011.

    [17] R. Mudumbai, G. Barriac, and U. Madhow, “On the feasibility ofdistributed beamforming in wireless networks,” IEEE Trans. WirelessCommun., vol. 6, no. 5, pp. 1754–1763, May 2007.

    [18] D. R. Brown and H. V. Poor, “Time-slotted round-trip carrier synchro-nization for distributed beamforming,” IEEE Trans. Signal Process.,vol. 56, no. 11, pp. 5630–5643, Nov. 2008.

    [19] D. R. Brown, B. Zhang, B. Svirchuk, and M. Ni, “An experimentalstudy of acoustic distributed beamforming using round-trip carriersynchronization,” in Proc. IEEE ARRAY 2010, Oct. 2010, pp. 316–323.

    [20] S. Wang, Q. Song, X. Wang, and A. Jamalipour, “Rate and poweradaptation for analog network coding,” IEEE Trans. Veh. Technol.,vol. 60, no. 5, pp. 2302–2313, Jun. 2011.

    [21] K. Lu, S. Fu, Y. Qian, and H.-W. Chen, “SER performance analysis forphysical layer network coding over AWGN channels,” in Proc. IEEEGLOBECOM 2009, Dec. 2009, pp. 1–6.

    [22] M. Park, I. Choi, and I. Lee, “Exact BER analysis of physical layernetwork coding for two-way relay channels,” in Proc. IEEE VTC 2011,May 2011, pp. 1–5.

    [23] M. Ju and I.-M. Kim, “Error performance analysis of BPSK modulationin physical-layer network-coded bidirectional relay networks,” IEEETrans. Commun., vol. 58, no. 10, pp. 2770–2775, Oct. 2010.

    [24] Y. Huang, Q. Song, S. Wang, and A. Jamalipour, “Symbol error rateanalysis for M-QAM modulated physical-layer network coding withphase errors,” in Proc. IEEE PIMRC 2012, Sep. 2012, pp. 2003–2008.

    [25] S. Wang, Q. Song, X. Wang, and A. Jamalipour, “Distributed MACprotocol supporting physical-layer network coding,” IEEE Trans. MobileComput., vol. 12, no. 5, pp. 1023–1036, May 2013.

    [26] F. M. Gardner, Phaselock Techniques, 3rd ed. Hoboken: John Wiley& Sons, 2005.

    [27] Maple 15. [Online]. Available: http://www.maplesoft.com[28] F. Golnaraghi and B. C. Kuo, Automatic Control System. New York:

    John Wiley & Sons, 2003.[29] H. Fu and P. Y. Kam, “Weighted phase averager for frequency estimation

    of a noisy single sinusoid: Application of the observation phase noisemodel,” in Proc. IEEE PIMRC 2009, Sep. 2009, pp. 1923–1927.

    [30] ——, “MAP/ML estimation of the frequency and phase of a singlesinusoid in noise,” IEEE Trans. Signal Process., vol. 55, no. 3, pp. 834–845, Mar. 2007.

    [31] D. Rife and R. Boorstyn, “Single tone parameter estimation fromdiscrete-time observations,” IEEE Trans. Inf. Theory, vol. 20, no. 5, pp.591–598, Sep. 1974.

    [32] G. Tavares, L. Tavares, and A. Petrolino, “On the true cramer-rao lowerbound for data-aided carrier-phase-independent frequency offset andsymbol timing estimation,” IEEE Trans. Commun., vol. 58, no. 2, pp.442–447, Feb. 2010.

    [33] S. Chang and B. Kelley, “An efficient time synchronization scheme forbroadband two-way relaying networks based on physical-layer networkcoding,” IEEE Commun. Lett., vol. 16, no. 9, pp. 1416 –1419, Sep. 2012.

    [34] W. U. Bajwa, J. Haupt, A. M. Sayeed, and R. Nowak, “Compressedchannel sensing: A new approach to estimating sparse multipath chan-nels,” Proc. IEEE, vol. 98, no. 6, pp. 1058–1076, Jun. 2010.

    [35] E. Aryafar, M. Khojastepour, K. Sundaresan, S. Rangarajan, and E. W.Knightly, “ADAM: An adaptive beamforming system for multicastingin wireless LANs,” in Proc. IEEE INFOCOM 2012, Mar. 2012, pp.1467–1475.

    [36] J. G. Proakis, Digital Communications, 4th ed. Boston: McGraw-Hill,2001.

    [37] M. K. Simon and M.-S. Alouini, Digital Communication over FadingChannels, 2nd ed. New York: Wiley, 2005.

    Yang Huang received the Bachelor’s degree fromNortheastern University, China, in 2011. Currently,he is pursuing the Master’s degree in Electronicsand Communication Engineering, at NortheasternUniversity, China. His general research interests liein communication systems, cooperative communica-tions, network coding, and radio resource manage-ment. He is a student member of the IEEE.

    Shiqiang Wang received the BEng and MEng de-grees from Northeastern University, China, in 2009and 2011, respectively. He is currently working to-ward the PhD degree in the Department of Electricaland Electronic Engineering, Imperial College Lon-don, United Kingdom. His research interests includenetwork coding, protocol design, optimization, andprototyping for wireless networks. He has a dozenscholarly publications in international journals andconferences. He served on the program committee ofIEEE VTC 2012-Fall, 2013-Spring, and 2013-Fall.

    Qingyang Song received the PhD degree in telecom-munications engineering from the University of Syd-ney, Australia. She is an associate professor at North-eastern University, China. She has authored morethan 30 papers in major journals and internationalconferences. These papers have been cited more than500 times in scientific literature. Her current researchinterests are in radio resource management, networkcoding, cognitive radio networks, and cooperativecommunications. She is a member of the IEEE.

    Lei Guo received the Ph.D. degree in communica-tion and information systems from School of Com-munication and Information Engineering, Universityof Electronic Science and Technology of China,Chengdu, China, in 2006. He is currently a professorin College of Information Science and Engineer-ing, Northeastern University, Shenyang, China. Hisresearch interests include optical networks, accessnetworks, network optimization and wireless com-munications. He has published over 200 technicalpapers in the above areas on international journals

    and conferences, such as IEEE Trans. Commun., IEEE/OSA J. LightwaveTechnol., IEEE Commun. Lett., IEEE GLOBECOM, IEEE ICC, etc. Dr. Guois a member of IEEE and OSA, and is also a senior member of China Instituteof Communications. He is now serving as an editor for three internationaljournals.

    Abbas Jamalipour (S’86-M’91-SM’00-F’07) re-ceived the Ph.D. degree from Nagoya University,Nagoya, Japan. He is the Chair Professor of Ubiqui-tous Mobile Networking with the School of Elec-trical and Information Engineering, University ofSydney, Sydney, NSW, Australia. He is a Fellowof the Institute of Electrical, Information, and Com-munication Engineers (IEICE) and the Institutionof Engineers Australia, an IEEE Distinguished Lec-turer, and a Technical Editor of several scholarlyjournals. He has been a Chair of several international

    conferences, including the IEEE International Conference on Communicationsand the IEEE Global Communications Conference, General Chair of the 2010IEEE Wireless Communications and Networking Conference, as well as beingthe technical program chair of IEEE PIMRC2012 and IEEE ICC2014. He isthe Vice President - Conferences and a member of Board of Governors ofthe IEEE Communications Society (ComSoc). He is the recipient of severalprestigious awards, including the 2010 IEEE ComSoc Harold Sobol Awardfor Exemplary Service to Meetings and Conferences, the 2006 IEEE ComSocDistinguished Contribution to Satellite Communications Award, and the 2006IEEE ComSoc Best Tutorial Paper Award.

    © 2013 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.