Multiple-Diﬀerential Encoding for Multi-Hop Amplify-and-Forward Relaying in IR-UWB Systems by Maziyar Hamdi B.Sc., Sharif University of Technology, Tehran, Iran, 2007 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF APPLIED SCIENCE in The Faculty of Graduate Studies (Electrical and Computer Engineering) THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver) February 2010 c Maziyar Hamdi 2010 Abstract In this thesis, we propose to improve the performance and coverage of impulse-radio ultra-wideband (IR-UWB) systems by means of cooperative multi-hop relaying. With regard to a simple practical realization, we focus on a non-coherent system setup in conjunction with amplify-and-forward (A&F) relaying. In particular, we propose to employ a multiple-diﬀerential encoding scheme at the source node and single diﬀerential decoding at each relay and at the destination node, respectively, so as to eﬃciently limit intersymbol-interference (ISI) eﬀects at the destination node. A thorough performance analysis of the proposed scheme is provided, along with a closed-form optimization of the transmit power allocation to the source node and the relays. Simulation results illustrate the excellent performance of the proposed scheme, which is also compared to alternative coherent and noncoherent multi-hop IR-UWB schemes based on A&F relaying and decodeand-forward (D&F) relaying. ii Table of Contents Abstract . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii Table of Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . iii List of Tables vi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . List of Figures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . vii List of Abbreviations Acknowledgments Dedication . . . . . . . . . . . . . . . . . . . . . . . . . ix . . . . . . . . . . . . . . . . . . . . . . . . . . . xi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xii 1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1.1 Ultra-wideband Communication 1.2 UWB Receivers 1 . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . 3 1.3 Cooperative Diversity in UWB Systems . . . . . . . . . . . . 6 1.4 Thesis Contributions . . . . . . . . . . . . . . . . . . . . . . 7 1.5 Thesis Organization . . . . . . . . . . . . . . . . . . . . . . . 9 1.6 Related Publications . . . . . . . . . . . . . . . . . . . . . . 10 iii Table of Contents 2 System Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 2.1 IR-UWB Transmission Format . . . . . . . . . . . . . . . . . 11 2.2 Proposed Structure of Source Node . . . . . . . . . . . . . . 14 2.3 Channel Model . . . . . . . . . . . . . . . . . . . . . . . . . . 16 2.4 Proposed A&F Relay Structure 18 2.5 Proposed Structure of the Destination Node . . . . . . . . . . . . . . . . . . . . . . . . . 21 3 Dual-hop Case . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 3.1 Eﬀective SNR at the Destination Node for the Dual-Hop Case 24 3.2 Optimized Transmit Power Allocation for the Dual-Hop Case 4 Multi-hop Case 4.1 . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 Eﬀective SNR at the Destination Node for the Multi-Hop Case . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4.2 28 30 Optimized Transmit Power Allocation for the Multi-Hop Case 33 5 Performance Analysis of Dual-hop Transmission with Rake Receiver . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 5.1 Eﬀective SNR at the Destination Node . . . . . . . . . . . . 38 5.2 Optimized Transmit Power Allocation . . . . . . . . . . . . . 44 6 Numerical Performance Results . . . . . . . . . . . . . . . . . 47 7 Conclusions and Future Work . . . . . . . . . . . . . . . . . 59 7.1 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . 59 7.2 Future Work . . . . . . . . . . . . . . . . . . . . . . . . . . . 60 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62 iv Table of Contents Appendices A Calculation of the Variance of the Noise at Each Relay for Multi-hop Transmission . . . . . . . . . . . . . . . . . . . . . . 68 A.1 Calculation of σz2 . . . . . . . . . . . . . . . . . . . . . . . 68 A.2 Calculation of σz2i+1 . . . . . . . . . . . . . . . . . . . . . . . 69 i+1 v List of Tables 6.1 Transmit power allocation factors αi∗ for the three-hop case (m = 3); the ﬁrst column indicates the relative positions of the relays (normalized w.r.t. the source-destination link length). . 6.2 57 Near-optimum transmit power allocation factors αi∗ for the four-hop case (m = 4); the ﬁrst column indicates the relative positions of the relays (normalized w.r.t. the sourcedestination link length). . . . . . . . . . . . . . . . . . . . . . 58 vi List of Figures 1.1 FCC spectral mask for indoor communication systems . . . . 4 2.1 Eﬀects of multi-hop relaying in straightforward A&F relaying 13 2.2 Block-Diagram . . . . . . . . . . . . . . . . . . . . . . . . . . 14 2.3 Diﬀerent channel realizations for CM1 model. . . . . . . . . . 17 2.4 Power decay proﬁle versus delay for CM1 channel model. . . 18 4.1 (a) Optimizing α1∗ using a recursive approach. (b) Optimization of α2∗ for known α1∗ . . . . . . . . . . . . . . . . . . . . . . 35 6.1 Pulse Shape (Time Domain). . . . . . . . . . . . . . . . . . . 48 6.2 Pulse Shape (Frequency Domain). . . . . . . . . . . . . . . . 49 6.3 SNR at the destination versus α1 for three diﬀerent values for the source-relay distances and the optimal value of α1 for each distance. . . . . . . . . . . . . . . . . . . . . . . . . . . . 6.4 50 BER at the destination versus Eg /N0 for double diﬀerential transmission and Rake receivers, respectively. A&F and D&F protocols. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52 vii List of Figures 6.5 BER at the destination versus Eg /N0 for double diﬀerential encoding and Rake receivers with both A&F and D&F protocols in coded case (hard input Viterbi decoding). . . . . . . 6.6 53 BER at the destination versus Eg /N0 for double diﬀerential encoding and Rake receivers with both A&F and D&F protocols in coded case (soft input Viterbi decoding). . . . . . . 54 6.7 BER vs. ρ. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55 6.8 Eﬀective SNR at the destination node for the three-hop and the four-hop case. . . . . . . . . . . . . . . . . . . . . . . . . . 56 viii List of Abbreviations AWGN Additive white Gaussian noise A&F Amplify-and-forward BER Bit error rate CDMA Code division multiple access CIR Channel impulse response D&F Decode-and-forward DF Diﬀerential DS-SS Direct sequence spread spectrum EGC Equal gain combining FCC Federal Communications Commission FEC Forward error correction GPS Global positioning system IPI Interpulse interference IR Impulse radio ISI Intersymbol-interference MAI Multiple access interference MB-OFDM Multiband orthogonal frequency-division multiplexing MC Multi-carrier MRC Maximum ratio combining ix List of Abbreviations PSD Power spectral density RF Radio frequency S-D Source-Destination SNR Signal-to-noise ratio TH Time-hopping TR Transmitted reference UWB Ultra wideband WPAN Wireless personal area network x Acknowledgments I owe my deepest gratitude to my supervisor, Dr. Robert Schober, for the constant guidance and support that he has provided. This thesis would not have been possible without his helpful suggestions, advice and constant encouragement. I would also like to thank Dr. Jan Mietzner for his support of and helpful comments on my research. I would like to extend my thanks as well to all of my friends and colleagues in the communication group for creating a friendly and productive environment. Furthermore, a heartfelt and deep thanks go to my family, and especially my brother Mazda, who have always encouraged me and given their unconditional support throughout my studies at the University of British Columbia. Without their support I could not have completed this thesis. Maziyar Hamdi The University of British Columbia February 2010 xi Dedication To my parents... xii Chapter 1 Introduction In this chapter, ﬁrst a brief overview of ultra-wideband (UWB) technology is provided in Section 1.1. After that, in Section 1.2, diﬀerent types of UWB receivers are discussed and some approaches in the literature for performance improvement of these receivers are summarized. Section 1.3 gives a literature survey about cooperative diversity with emphasis on its applications in UWB systems. Then, in Sections 1.4 and 1.5, respectively, we summarize the contributions and outline the organization of thesis. Finally, publications related to this thesis are listed in Section 1.6. 1.1 Ultra-wideband Communication UWB is a wireless spectral underlay technology for transmitting signals with a bandwidth larger than 500 MHz or a fractional bandwidth of more than 20%. Here, the fractional bandwidth is deﬁned by B/fc , where B and fc denote the −10 dB bandwidth and the centre frequency [1–4], respectively. Furthermore, B = fH −fL , and fc = fH +fL , 2 where fH and fL are the higher and lower frequencies for −10 dB gain [3], respectively. The above deﬁnition of UWB signals implies that in UWB communication systems with a centre frequency, fc , greater than 2.5 GHz the signal bandwidth should be larger 1 1.1. Ultra-wideband Communication than 500 MHz. For lower centre frequencies, the bandwidth has to satisfy B/fc > 0.02 [5]. Owing to their large bandwidth, UWB transmission techniques are, for example, envisioned for short-range high-speed indoor communications. Due to the favorable characteristics of UWB communications such as low-cost, simple implementation, and high data rates, UWB-based systems are being used in many applications such as wireless personal area networks (WPAN), sensor networks, and vehicular radar systems. Currently, there are three different standards or standard proposals for UWB radio [6], which are quite diﬀerent in nature, namely (i) the ECMA-368 standard [7], which is based on a multiband orthogonal frequency-division multiplexing (MB-OFDM) modulation scheme, (ii) the direct-sequence (DS) UWB standard proposal [8], which is based on a DS spread-spectrum (DS-SS) scheme, and (iii) the IEEE 802.15.4a standard [9], which is based on the (rather traditional) impulseradio (IR) technique. In multicarrier (MC) based UWB systems, the concepts in direct-sequence spread spectrum (DS-SS) or code-division multipleaccess (CDMA) are expanded with MC techniques. Conventional DS-SS can also be used to produce a UWB signal if a suﬃciently high chip rate issued. MC schemes have a higher receiver complexity because of multi-stage mixing or Fourier transform [6]. Impulse radio is another form to implement time-modulated UWB. In IR-based UWB (IR-UWB) systems, the transmitted signal consists of a train of pulses of very short duration [10]. By this means, an ultrawide bandwidth is generated, without requiring any additional carrier modulation. One of the most attractive characteristics of IR-UWB systems is their 2 1.2. UWB Receivers low complexity because of their baseband nature . An UWB transmitter can send data by modulating short pulses and radio frequency (RF) mixing is not required [11]. So the implementation of IR-UWB systems is low-cost compared with conventional communications systems. As a result of these ultra short pulses, IR-UWB is widely used for high speed and high data rate communications. Hundreds of mega bit per second (Mbps) are achievable with IR-UWB communication systems. Due to its simple practical realization and its robustness to multipath fading and intersymbol-interference (ISI) eﬀects, IR-UWB has attracted considerable attention. Moreover, by deploying pseudo-random time-hopping (TH) sequences, IR-UWB can also support multiuser communications with a minimum of multiple-access interference (MAI). 1.2 UWB Receivers The US Federal Communications Commission (FCC) has introduced tight limitations (below -41.3 dBm/Hz) on the transmitted power spectral density (PSD) of UWB systems to restrict interference to incumbent wireless services [5]. The FCC spectral mask, Fig. 1.1, allows UWB systems to operate in frequency bands that have been allocated to licensed narrowband systems. For example, in the 960 MHz-1.61 GHz band, where the global positioning system (GPS) operates, UWB system can send data with a power level lower than −76.3 dBm. Because of these limitations, it is essential to capture at the receiver most of the signal energy provided by the large number of resolvable multipath 3 1.2. UWB Receivers UWB EIRP Emission Level [dBm/MHz] −40 −45 −50 −55 −60 −65 Indoor limit Part 15 limit −70 −75 9 10 f [Hz] 10 10 Figure 1.1: FCC spectral mask for indoor communication systems components. A favorable property of IR-UWB systems is that they allow for an eﬃcient energy combining at the receiver [12] – either by means of coherent Rake combining [13] or non-coherent energy detection schemes [14, 15]. In the case of Rake combining, the received signal is correlated with an appropriate template signal, which is designed for coherent combination of the energies of the individual multipath signal components. An ideal all-Rake (A-Rake) receiver captures the energy of all resolvable multipath components and thus oﬀers optimum performance in terms of bit error rate (BER). However, it requires accurate channel estimation at the receiver and 4 1.2. UWB Receivers precise timing synchronization, which might be challenging in practice. In particular, a large number of “Rake ﬁngers” is typically required in order to capture the entire signal energy [16]. As opposed to this, partial Rake (P-Rake) and selective Rake (S-Rake) combiners collect only part of the signal energy, thus providing suboptimal performance at the beneﬁt of a reduced receiver complexity. Still, accurate channel knowledge and precise timing synchronization are indispensable. In S-Rake receivers, the Ls most signiﬁcant multi-path components are collected while the P-Rake receiver combines the Lp ﬁrst multi-path components. As opposed to this, non-coherent energy detection schemes relieve the receiver from any channel estimation task and are thus easier to realize in practice. They capture the energy of the multipath components by means of an autocorrelation of the received signal, followed by an integrate-and-dump operation. Among these techniques, diﬀerential (DF) and transmitted- reference (TR) schemes are the most popular options [14]. In the TR scheme, pulses are transmitted in pairs, where the ﬁrst pulse serves as reference pulse and the second pulse, which follows the reference pulse with a ﬁxed delay, is modulated by the transmitted information bit sequence. At the receiver, the reference pulse serves as a noisy template for the autocorrelation process, so that the transmitted information bit modulated onto the second pulse can be recovered by means of the subsequent integrate-and-dump operation. In the DF scheme, however, the information bit sequence is diﬀerentially encoded at the transmitter, and the autocorrelation process at the receiver (in conjunction with the integrate-and-dump operation) performs the corresponding diﬀerential decoding step. A notable advantage of the DF scheme 5 1.3. Cooperative Diversity in UWB Systems is that it oﬀers a higher data rate compared to a TR scheme [14], since the additional template pulse can be discarded. In the literature, several methods have been proposed to improve the performance of DF schemes, such as reference ﬁltering [17], weighted correlation [18], and multiple symbol detection [19]. 1.3 Cooperative Diversity in UWB Systems Cooperative diversity is another way to overcome the limited coverage of (coherent or non-coherent) UWB systems [20, 21]. The idea of cooperative diversity was ﬁrst proposed in [22] and [23]. Through cooperative diversity, wireless network terminals can create redundancy by relaying message for each other. The receiver at the destination can use these redundant signals to overcome multi-path fading and shadowing. Cooperative diversity improves network performance and improves the achievable rate by use of diﬀerent cooperative diversity protocols. Two popular schemes for cooperative relaying are amplify-and-forward (A&F) and decode-and-forward (D&F) relaying. A&F relays simply amplify and re-transmit the received signal, whereas D&F relays ﬁrst decode and then re-encode the received signal, before re-transmission is performed [24]. Correspondingly, D&F relaying is more complex than A&F relaying, especially if a forward-error-correction (FEC) code is employed. The diversity gain of A&F and D&F relay was investigated in [25]. Coherent and non-coherent receivers for these schemes were described in [26] and [27], respectively. End-to-end performance, BER, and outage probability in Rayleigh fading channels for a simple two-hop sys- 6 1.4. Thesis Contributions tem were investigated in [28]. The performance analysis has been extended to Nakagami fading in [29]. 1.4 Thesis Contributions In this thesis, we investigate the use of relays to overcome the limited coverage of (coherent or non-coherent) UWB systems [20, 21, 30–32]. Instead of transmitting signals directly from a source node, S, to a destination node, D, the signal is (in the simplest case) received by an intermediate relay, R1 , which forwards the received signal from the source node to the destination node. In a more complicated multi-hop scenario consisting of m subsequent links, the source signal is forwarded between (m−1) relays R1 , ..., Rm−1 , before it ﬁnally reaches the destination node. With relaying, substantial path-loss gains can be provided due to shorter link lengths. Available relays can either be dedicated relays, which do not disseminate any data of their own, or temporarily inactive IR-UWB devices acting as relays, so as to assist the current source-destination link. Previous work on relaying for UWB systems has focused on the ECMA368 UWB standard [20, 21], space time-code design for coherent IR-UWB systems [30], coherent A&F and D&F relaying requiring a Rake combiner at each relay [31, 32], and dual-hop D&F relaying for IR-UWB with singlediﬀerential encoding at the source [32]. For the sake of a simple practical realization, we will focus on a combination of the DF IR-UWB scheme with A&F relaying. A straightforward combination of the two techniques with diﬀerential encoding at the source 7 1.4. Thesis Contributions node, simple A&F relaying at the intermediate relay(s), and diﬀerential decoding at the destination node has the major drawback that the length of the eﬀective overall channel impulse response (CIR) increases with each hop. Compared to direct transmission (i.e., without any relaying), signiﬁcantly larger guard intervals between the transmitted pulses are therefore required, in order to achieve a similar level of ISI. Otherwise, the increased amount of ISI will severely limit the overall performance and can even compromise the achieved path-loss gains. Note that an increased guard interval will signiﬁcantly lower the eﬀective data rate compared to direct transmission. As an alternative, we propose a multiple-diﬀerential encoding setup with m-fold diﬀerential encoding at the source node in conjunction with single diﬀerential decoding at each A&F relay and at the destination node. By this means, the same level of ISI is achieved as in the case of direct transmission, without requiring an extended guard interval. To the best of our knowledge, such a use of multiple-diﬀerential encoding is novel and quite diﬀerent from narrowband systems, where, for example, double-diﬀerential encoding is employed to mitigate carrier frequency oﬀsets [33]. The main contributions of this thesis are as follows: • We propose the use of double diﬀerential encoding at the transmitter for dual-hop relaying to limit ISI without any loss in terms of data rate. • For the dual-hop case we derive a simple closed-form expression for the power allocation factors for the transmitter and the relay, which only depends on the length of source-relay and relay-destination links. 8 1.5. Thesis Organization • We generalize the idea of using diﬀerential encoding in order to mitigate ISI to multi-hop UWB systems and propose to use multiple differential encoding at the transmitter. • We derive a recursive formula for the signal-to-noise ratio (SNR) at the destination for the multi-hop case. • We propose a sub-optimal, recursive approach for optimizing the power allocation factors in multi-hop UWB systems. 1.5 Thesis Organization The rest of the thesis is organized as follows. In Chapter 2, we describe the system setup. The structures of transmitter, channel, relay, and receiver for multi-hop transmission are discussed in this chapter. In Chapter 3, we consider a dual-hop system and analyze its performance in terms of SNR and derive a closed-form expression for the optimum power allocation which maximizes the SNR at the destination. A recursive formula for the SNR of multi-hop IR-UWB systems and a sub-optimal approach for optimizing the corresponding power allocation factors at the source and the relays are provided in Chapter 4. In Chapter 5, the structure of the Rake receiver is reviewed and the optimum power allocation for dual-hop transmission with Rake combining is derived. Simulation results, which illustrate the excellent performance of our scheme, are presented in Chapter 6. Finally, Chapter 7 concludes the thesis and presents possible directions for future work. 9 1.6. Related Publications 1.6 Related Publications • M. Hamdi, J. Mietzner, and R. Schober: Double Diﬀerential Encoding for Dual-Hop Amplify-and-Forward Relaying in IR-UWB Systems. Accepted for presentation at IEEE Vehicular Technology Conference (VTC). Taipei, Taiwan, March 2010 • M. Hamdi, J. Mietzner, and R. Schober: Multiple Diﬀerential Encoding for Multi-Hop Amplify-and-Forward Relaying in IR-UWB Systems. To be submitted to IEEE Transactions on Wireless Communications. Manuscript is 95 % complete. 10 Chapter 2 System Model Throughout this thesis, we consider a serial multi-hop DF IR-UWB system consisting of a source node S, (m−1) A&F relays R1 , ..., Rm−1 , and a destination node D. Here, m > 1 denotes the total number of hops. The link between source S and relay R1 is in the sequel denoted as Link 1, the link between relay Ri−1 and relay Ri is denoted as Link i (2 ≤ i ≤ m−1), and the link between the last relay Rm−1 and the destination node D is denoted as Link m. For the ease of exposition, we focus on the single-user case here. Before we provide a detailed description of the system setup, including the underlying channel model as well as the structures of the source node, the relay node(s), and the destination node, we brieﬂy recapitulate the IRUWB transmission format and highlight the advantages of our proposed multiple-diﬀerential A&F relaying scheme. 2.1 IR-UWB Transmission Format In IR-UWB systems, the transmitted signal consists of a train of ultra-short pulses (on the order of nanoseconds), which are modulated by the transmitted information symbols [10]. Within the scope of this thesis, focus will be on pulse-amplitude modulation (PAM). As illustrated in Fig. 2.1 (a), Nf ≥ 1 11 2.1. IR-UWB Transmission Format frames are used to convey a single information symbol (in the example, we have Nf = 2). By this means, the eﬀective symbol energy at the receiver can be increased, however, at the expense of a decreased throughput. After transmission, the IR-UWB signal is convolved with the UWB CIR, which can have a length on the order of hundreds of nanoseconds [4], depending on the radio environment under consideration. Correspondingly, the received pulses are considerably spread out in time, cf. Fig. 2.1 (b). In the DF setup, the energy of the underlying multipath signal components is collected at the receiver by means of the integrate-and-dump operation (for details see Section 2.4). If a straightforward combination of the IR-UWB transmission format with A&F relaying is employed, the eﬀective CIR seen at the destination node results from a convolution of the UWB CIRs of all intermediate links. Correspondingly, the time spread of the received pulses will grow with each hop, leading to signiﬁcant interpulse interference (IPI) and ISI, as the received pulses associated with one information symbol (i.e., with one set of Nf frames) interfere with those associated with the subsequent symbol. In order to keep the level of IPI and ISI comparable to the single-hop case, the only option is to increase the guard intervals between subsequent pulses and symbols, which will (further) decrease the eﬀective throughput. In the following, we present the details of our proposed multiple-diﬀerential A&F relaying scheme. A major advantage of our scheme is that the level of IPI and ISI is not increased compared to the single-hop case. In particular, since diﬀerential decoding (in conjunction with the integrate-and-dump operation) is performed at each relay, the time spread of the received pulses 12 2.1. IR-UWB Transmission Format Figure 2.1: Eﬀects of relaying in straightforward dual-hop A&F IR-UWB: (a) transmitted signal, (b) received signal after one hop, and (c) received signal after two consecutive hops without intermediate diﬀerential decoding. depends only on a single UWB CIR, so that the same guard interval size can be employed as in the case of a single hop. 13 2.2. Proposed Structure of Source Node 2.2 Proposed Structure of Source Node Fig. 2.2 shows the proposed structures of the source node S, the ith A&F relay Ri , and the destination node D. Throughout this thesis, the employed relays R1 , ..., Rm−1 are assumed to have identical structures. We start with a description of the source node structure. Figure 2.2: Block diagram of the (a) transmitter, (b) relay, and (c) destination. At the source node, the transmitted information bits, b[κ] ∈ {0, 1}, are 14 2.2. Proposed Structure of Source Node (possibly) ﬁrst encoded by a forward-error-correction (FEC) encoder (see Fig. 2.2 (a)), which yields coded bits c[k]. These are then mapped onto binary antipodal symbols q1 [k] ∈ {±1}, according to the mapping rule 0 → −1, 1 → +1. Next, the information symbols q1 [k] are m-times diﬀerentially encoded according to qr [k] = qr−1 [k]qr [k − 1], 2 ≤ r ≤ m + 1, (2.1) where qr [k] ∈ {±1} denotes the kth intermediate symbol after the (r −1)th diﬀerential encoder (2 ≤ r ≤ m+1). The symbols qm+1 [k] are then modulated onto a train of (real-valued) short pulses wtx (t) with duration Tp , according to [3] ∞ s1 (t) = Nf −1 qm+1 [k] wtx (t − jTf − kTs ). α1 Eg (2.2) k=−∞ j=0 Here, s1 (t) denotes the transmitted signal, α1 the transmit power allocation factor for the source node, Eg the energy per pulse ( +∞ 2 −∞ wtx (t) dt := 1), Nf the number of frames used for conveying a single information symbol, Tf the frame duration, which is chosen much larger than the pulse duration (i.e., Tf Tp ), and Ts := Nf Tf the symbol duration. The eﬀective bit rate is given by Rb := 1/Ts . The length of the guard interval between subsequent pulses and subsequent symbols, Tg := Tf − Tp , is typically chosen longer than the (expected) length of the UWB CIR, so as to circumvent IPI and ISI eﬀects. 15 2.3. Channel Model 2.3 Channel Model In the following, Link i denotes the link between node (i) and node (i + 1). For all links under consideration, we employ the IEEE 802.15.3a channel models [34] for UWB personal area networks. Consequently, the passband version of the CIR of the rth link consists of LC clusters of LR rays and is modeled as LC LR hr (t) = (r) (r) λ(r) μ,ν δ(t − Tν − τμ,ν ), (2.3) ν=0 μ=0 (r) where λμ,ν models the random multipath gain coeﬃcient of the μth ray (r) of the νth cluster, Tν (r) the delay of the νth cluster, τμ,ν the delay of the μth ray of the νth cluster, and δ(·) denotes a Dirac impulse. The multipath gain coeﬃcients are normalized such that LC ν=0 2(r) LR μ=0 λμ,ν = 1. In [34], four diﬀerent parameter sets are speciﬁed for the various parameters in (2.3). The resulting channel models, CM1-CM4, represent diﬀerent usage scenarios and entail diﬀerent amounts of IPI and ISI. Fig. 2.3 shows one hundred realizations of the CIR of an UWB link for model CM1. Fig. 2.4 depicts the average power decay proﬁle versus delay time. As shown in these ﬁgures, the maximum excess delay is about 60 ns so with a guard interval of about 70 ns, eﬀects of ISI can be mitigated. The channel gain is aﬀected by log-normal fading and path loss, which is usually modeled in the UWB literature according to [35, 36] G(d) = G0 + 10 · p log(d0 /d) + χ, (2.4) where G(d) denotes the channel gain in dB, d the link length in meter, G0 16 2.3. Channel Model Impulse response realizations 1 0.8 0.6 0.4 Amplitude 0.2 0 −0.2 −0.4 −0.6 −0.8 −1 0 20 40 60 80 100 120 t[ns] Figure 2.3: Diﬀerent channel realizations for CM1 model. the channel gain resulting for some reference distance d0 (e.g., d0 = 1m), p the path-loss exponent, and χ a log-normal fading term. In the following, the source-destination (S-D) link will serve as a reference link for the path loss, in order to allow for a fair comparison between the proposed multihop relaying setup and the case of direct transmission. Assuming omnidirectional antennas at all nodes, the relative channel gain associated with the ith link is thus modeled by a factor [21] Ai := χi dS−D di p , (2.5) 17 2.4. Proposed A&F Relay Structure Average Power Decay Profile 0 −10 Average Power [dB] −20 −30 −40 −50 −60 0 20 40 60 80 100 Delay [ns] Figure 2.4: Power decay proﬁle versus delay for CM1 channel model. where χi models the log-normal fading associated with the ith link, and dS−D and di denote the lengths of the source-destination link and the ith link, respectively. The path-loss exponent p is typically between 1.7 ≤ p ≤ 3.5 [34]. Throughout this thesis, we assume that the lognormal shadowing terms χi and χi associated with two diﬀerent links i = i are uncorrelated. 2.4 Proposed A&F Relay Structure At the receiver front-end of each relay, the received signal is ﬁrst passed through a bandpass ﬁlter hBP (t) with one-sided bandwidth W , so as to 18 2.4. Proposed A&F Relay Structure eliminate out-of-band noise (see Fig. 2.2 (b)). The ﬁltered received signal at the ﬁrst relay is given by ri (t) = Ai hi (t) ∗ hBP (t) ∗ si (t) + ni (t) ∞ = Nf −1 ξi [k] wrx,i (t − jTf − kTs ) + ni (t), (2.6) Ai αi Eg k=−∞ j=0 where Ai is the relative channel gain associated with Link i, hi (t) the corresponding CIR, si (t) the transmitted signal of the source node (i = 1) or the (i−1)th relay (1 < i < m − 1), ni (t) the ﬁltered additive white Gaussian noise (AWGN) with zero mean and single-sided noise PSD N0 2 , αi the trans- mit power allocation factor for the source node (i = 1) or the (i−1)th relay (1 < i ≤ m − 1), wrx,i (t) = wtx (t) ∗ hi (t) ∗ hBP (t) the received pulse associated with Link i, and ‘∗’ denotes linear convolution. Moreover, ξi [k] := ⎧ ⎪ ⎨ qm+1 [k] for i = 1 ⎪ ⎩ x ˆi−1 [k] for 1 < i < m − 1 , (2.7) where qm+1 [k] denotes the kth transmitted symbol of the source node and x ˆi−1 [k] the soft estimate of the kth intermediate symbol qm−i+2 [k] formed at the (i−1)th relay (1 < i < m − 1), see details below. For the numerical results presented in Chapter 6, the bandwidth W of the bandpass hBP (t) was optimized numerically for maximization of the SNR. After the bandpass ﬁlter, (single) diﬀerential decoding of the ﬁltered receive signal ri (t) is performed. To this end, signal ri (t) is ﬁrst delayed by a symbol duration Ts and is then multiplied by itself. The resulting 19 2.4. Proposed A&F Relay Structure signal ri (t)ri (t − Ts ) is passed through an integrator with basic integration duration Ti .1 Throughout this thesis, we assume perfect synchronization of the relays and the destination node with respect to the employed frame structure [14]. In particular, the overall integration time is composed of Nf separated sub-region intervals (one per received frame), as shown in Fig. 2.1. By this means, the negative impact of the additive noise can be eﬃciently limited (e.g., in comparison to a single integration interval of length Nf Tf ). The integrator yields the discrete-time output sample Nf −1 x ˆi [k] = kTs +jTf +Ti kTs +jTf j=0 ri (t)ri (t − Ts )dt, (2.8) which can be interpreted as a (soft) estimate of the intermediate symbol qm−i+1 [k] = qm−i+2 [k] · qm−i+2 [k − 1] (2.9) (due to the diﬀerential decoding step, see Chapter 3 for further details). Similar to (2.2), the estimated symbols x ˆi [k] are ﬁnally modulated onto a signal ∞ si+1 (t) = Nf −1 x ˆi [k] wtx (t − jTf − kTs ), αi+1 Eg (2.10) k=−∞ j=0 which is then re-transmitted to the next relay (1 ≤ i ≤ m − 2) or the destination node (i = m − 1). Here, αi+1 denotes the power allocation 1 Similar to the ﬁlter bandwidth W , the integration duration Ti has to be optimized such that (on average) most of the signal energy is captured, while the collected noise energy is kept to a minimum. 20 2.5. Proposed Structure of the Destination Node factor for the ith relay node. In the case of D&F relaying, a hard decision on the symbols xˆi [k] would be performed prior to re-transmission, followed by a FEC decoding step. Note that D&F relaying is thus more complex than A&F relaying. 2.5 Proposed Structure of the Destination Node The receiver structure of the destination node is identical to that of the relays. In particular, we assume an identical bandpass ﬁlter hBP (t) for simplicity. Similar to (2.6), the ﬁltered received signal is given by rm (t) = Am (hm (t) ∗ hBP (t)) ∗ sm (t) + nm (t) ∞ = Nf −1 x ˆm−1 [k] wrx,m (t − jTf − kTs ) + nm (t), Am αm Eg k=−∞ j=0 (2.11) where wrx,m (t) := wtx (t) ∗ hm (t) ∗ hBP (t) and hm (t) denotes the CIR of the ﬁnal link (Link m). In order to recover the symbols q1 [k] transmitted by the source node, the destination node performs another diﬀerential decoding step followed by an integrate-and-dump operation. This yields estimated symbols Nf x ˆm [k] = j=0 kTs +jTf +Ti kTs +jTf rm (t)rm (t − Ts )dt (2.12) (cf. Fig. 2.2 (c)), which are used for decoding of the transmitted information bits. At the destination node, the Viterbi algorithm is used for decoding the output of the integration-and-dump unit. We consider both hard and soft 21 2.5. Proposed Structure of the Destination Node input Viterbi decoding. In hard decoding, xˆm [k] is passed through a slicer, in order to arrive at the hard estimates qˆ1 [k] of the transmitted symbols q1 [k]. Subsequently, Viterbi decoding is performed based on qˆ1 [k] (‘FEC−1 ’), in order to obtain estimates ˆb[κ] for the transmitted information bits b[κ]. In contest, in the soft input Viterbi decoding, the soft estimates, xˆm [k], are used to obtain estimates of the transmitted bits. The squared Euclidean distance is used as a metric for soft input Viterbi decoding. In Chapter 4, it will be shown that the input-output behavior of the system can be described as x ˆm [k] = β˜m q1 [k] + zm [k], where β˜m and zm [k] are a gain factor and a noise sample, respectively (c.f. Section 4.1). Thus, we use the following branch metric for the soft input Viterbi decoder [37] ι[k] = |ˆ xm [k] − β˜m q˜1 [k]|2 , (2.13) where q˜1 [k] ∈ {±1} denotes a trial symbol. 22 Chapter 3 Dual-hop Case In the following, we analyze the performance of dual-hop IR-UWB with double diﬀerential encoding at the source. In this case, we have only one relay i.e., m = 2. We show that the output of this system can be written in terms of the input as x ˆ3 [k] = β˜3 q1 [k] + z3 [k], (3.1) where β˜3 and z3 [k] represent a gain factor and an eﬀective noise sample, respectively. Note that w3 [k] is in general non-Gaussian, so ﬁnding a closedform expression for the BER seems diﬃcult. Therefore, we focus on the eﬀective SNR at the destination node in Section 3.1. In particular, a closedform expression for the SNR at the destination node is derived. In Section 3.2, keeping the total power ﬁxed, an optimum solution for power allocation factors α1 and α2 is devised to maximize the SNR at the destination. Based on the results of this section, we propose a method for optimization of the power allocation factors of more complex multi-hop systems in Chapter 4. 23 3.1. Eﬀective SNR at the Destination Node for the Dual-Hop Case 3.1 Eﬀective SNR at the Destination Node for the Dual-Hop Case We start by substituting the ﬁltered received signal r1 (t) at the relay, cf. (2.6) (i = 1), into (2.8) and get =q2 [k] x ˆ1 [k] = β1 q3 [k]q3 [k − 1] +z1 [k], Ti 0 where β1 := Nf A1 α1 ε1 , ε1 := Eg (3.2) 2 (t)dt, and z [k] := z [k] + z wrx,1 1 1,1 1,2 + z1,3 [k] with z1,1 [k], z1,2 , and z1,3 [k] being zero-mean noise terms deﬁned as follows: Nf −1 z1,1 [k] := q3 [k] A1 α1 Eg j=0 kTs +jTf +Ti kTs +jTf n1 (t − Ts )wrx,1 (t − jTf − kTs ) dt, (3.3) Nf −1 z1,2 [k] := q3 [k − 1] A1 α1 Eg j=0 kTs +jTf +Ti kTs +jTf n1 (t)wrx,1 (t − jTf − kTs ) dt, (3.4) Nf −1 z1,3 [k] := j=0 kTs +jTf +Ti kTs +jTf n1 (t)n1 (t − Ts ) dt. (3.5) As described earlier, n1 (t) is obtained by ﬁltering an AWGN process with single-sided noise PSD N0 2 with a bandpass ﬁlter with one-sided bandwidth W . The autocorrelation function (ACF) φ1 (τ ) of n1 (t) is thus given by [37] φ1 (τ ) = E{n1 (t)n1 (t − τ )} = N0 sin(πW τ ) · cos(2πf0 τ ), 2 πW τ (3.6) 24 3.1. Eﬀective SNR at the Destination Node for the Dual-Hop Case where f0 denotes the center frequency of the bandpass ﬁlter, and E{·} denotes statistical expectation. Assuming that the bandwidth W is chosen suﬃciently large, such that the frequency response of the received pulse wrx,1 (t) falls completely inside the PSD Φ1 (f ) of n1 (t), and the PSD Φ1 (f ) is suﬃciently ﬂat in the area of interest, the ACF φ1 (τ ) can be replaced by N0 2 δ(τ ). Thus, the variances of the noise terms z1,1 [k], z1,2 [k], and z1,3 [k] can be approximated as σz2 1,1 2 = E{z1,1 [k]} = Nf A1 α1 Eg ≈ Nf A1 α1 Eg σz2 1,2 σz2 1,3 Ti N0 2 0 Ti 0 Ti 0 wrx,1 (t)wrx,1 (τ )φ1 (t − τ )dtdτ 2 wrx,1 (t)dt = β1 N0 , 2 (3.7) 2 = E{z1,2 [k]} = σz2 , (3.8) 1,1 2 = E{z1,3 [k]} = Nf Ti = Nf 0 Ti 0 Ti 0 φ21 (t − τ )dtdτ = Nf Ti 0 Ti Ti −t φ21 (u)dtdu, W Nf Ti N02 N02 · 2 · W dt = · 4 2 (3.9) In (3.9), we have exploited the fact that the integrand φ21 (u) is Dirac-like, i.e., the integral vanishes outside [−t, Ti − t], and we have employed Parseval’s theorem to calculate the integral. Altogether, the variance of the noise term z1 [k] = z1,1 [k] + z1,2 + z1,3 [k] in (3.2) is given by σ12 := β1 N0 + W Nf Ti N02 . 2 Note that z1 [k] is not Gaussian distributed, since z1,3 [k] is not Gaussian. Along the same lines, the integrator output of the destination node can be analyzed. Substituting (2.11) into (2.12) (i = 1, m = 2) yields x ˆ2 [k] = β2 x ˆ1 [k]ˆ x1 [k − 1] + z2 [k], (3.10) 25 3.1. Eﬀective SNR at the Destination Node for the Dual-Hop Case where β2 := Nf A2 α2 ε2 , ε2 := Eg Ti 0 2 (t)dt, and the noise term z [k] = wrx,2 2 z2,1 [k] + z2,2 [k] + z2,3 [k] is similarly deﬁned as in (3.2) and (3.3)–(3.5), with ˆ1 [k], A1 by A2 , α1 by α2 , n1 (t) by n2 (t), and wrx,1 (t) q3 [k] being replaced by x by wrx,2 (t). Plugging x ˆ1 [k] from (3.2), into (3.10) we get =q1 [k] x ˆ2 [k] = β2 β12 q2 [k]q2 [k − 1] +z2 [k] + z2,4 [k] + z2,5 [k] + z2,6 [k], (3.11) where z2,4 [k], z2,5 [k], and z2,6 [k] are deﬁned as z2,4 [k] = β2 β1 q2 [k]z1 [k − 1], (3.12) z2,5 [k] = β2 β1 q2 [k − 1]z1 [k], (3.13) z2,6 [k] = β2 z1 [k]z1 [k − 1]. (3.14) The variances of z2,1 [k], z2,2 [k], and z2,3 [k] can be calculated similarly to (3.7)–(3.9). One obtains σz2 2,1 ≈ E{ˆ x21 [k]}β2 σz2 2,2 N0 N0 = (β12 + σ12 )β2 , 2 2 = σz2 , σz2 2,1 2,3 = σz2 . (3.15) 1,3 The variances of z2,4 [k], z2,5 [k], and z2,6 [k] can be calculated from (3.12)– (3.14): σz2 2,4 = σz2 2,5 = (β2 β1 )2 σ12 and σz2 2,6 ≈ β22 σ14 . (3.16) For computing the variance of z2,6 [k], we have assumed that z1 [k] and 26 3.1. Eﬀective SNR at the Destination Node for the Dual-Hop Case z1 [k − 1] are statistically independent. Based on (3.11)–(3.16), we are now ready to obtain the input-output behavior of our system according to (3.1), where β˜2 := β2 β12 and z2 [k] := z2 [k] + z2,4 [k]+ z2,5 [k] + z2,6 [k]. In order to calculate the variance of z2 [k], denoted as σ22 in the sequel, we note that z2,4 [k] and z2,5 [k] are not mutually independent, i.e., their correlation cannot be ignored. In particular, one ﬁnds that E{z2,4 [k]z2,5 [k]} = (β2 β1 )2 E{q2 [k]q2 [k − 1]z1 [k]z1 [k − 1]} = (β2 β1 )2 E{q2 [k]q2 [k − 1]z1,1 [k]z1,2 [k − 1]}. (3.17) A careful look at (3.3)–(3.3) reveals that q3 [k −2]z1,2 [k −1] = q3 [k]z1,1 [k]. Thus, we have E{q2 [k]q2 [k − 1]z1,1 [k]z1,2 [k − 1]} = E{q3 [k]q3 [k − 2]z1,1 [k]z1,2 [k − 1]} 2 [k]} = β1 = E{z1,1 N0 2 (3.18) and get E{z2,4 [k]z2,5 [k]} = β22 β13 N0 . 2 (3.19) Based on (3.1) and (3.11)–(3.19), the eﬀective SNR at the destination node can be calculated as SNR = 2 β˜2 = 2 σ22 σz 2,4 ≈ + σz2 2,5 + σz2 2,6 (β2 β12 )2 + 2E{z2,4 [k]z2,5 [k]} + σz2 (β2 β12 )2 β22 2β12 σ12 + β13 N0 + σ14 + β2 N0 (β12 + σ12 ) + 2 W Nf Ti N02 2 ·(3.20) 27 3.2. Optimized Transmit Power Allocation for the Dual-Hop Case 3.2 Optimized Transmit Power Allocation for the Dual-Hop Case Based on (3.20), we can now optimize the transmit power allocation factors α1 and α2 for the source node and the A&F relay, respectively. We aim to maximize the eﬀective SNR at the destination node, under the constraint of keeping the total transmit power ﬁxed. It can be shown that the total transmit power constraint can be expressed as ! α1 + α2 (β12 + σ12 ) = 1. (3.21) Concerning (3.20), we can make a high-SNR approximation, according to SNR ≈ β2 (2σ12 β2 β12 , + β1 N0 ) + N0 where we have used that, for high SNR values, 2β12 σ12 +β13 N0 and W Nf Ti N02 2 (3.22) σ14 , β12 σ12 , becomes negligible. Assuming again high SNR, we can now formulate the Lagrange problem Λ(α1 , α2 , γ) = β2 β12 − γ(α1 + α2 β12 − 1), (3β1 β2 + 1)N0 (3.23) where γ denotes the Lagrange multiplier. In order to arrive at (3.23), we have used the high-SNR approximations β12 σ12 and σ12 ≈ β1 N0 . Based on (3.23), the optimal transmit power allocation results as α1,opt = 1 1+ A1 3A2 , α2,opt = 1 − α1,opt , β12 + σ12 (3.24) 28 3.2. Optimized Transmit Power Allocation for the Dual-Hop Case which only depends on the link gains A1 and A2 . The numerical results in Chapter 6 will reveal that (3.24) is, in fact, very close to the optimal solution (also for moderate SNR values). 29 Chapter 4 Multi-hop Case In the following, we turn to the case of multiple hops (m > 2). In particular, we derive a closed-form expression for the eﬀective SNR at the destination node and use this expression for optimization of the power allocation factors for the source node and the A&F relays. 4.1 Eﬀective SNR at the Destination Node for the Multi-Hop Case Similar to (3.1), the input-output behavior of the system can be described as x ˆm [k] = β˜m q1 [k] + zm [k], (4.1) where β˜m and zm [k] represent a gain factor and an eﬀective noise sample, respectively. We note again that zm [k] is in general non-Gaussian. In order to obtain an expression for the gain factor β˜m , the noise term zm [k], and 2 (and ﬁnally the eﬀective SNR β ˜2 /σ 2 the corresponding noise variance σm m m at the destination node), the following recursion steps are required (similar to the derivation for the dual-hop case in Section 3.1): • Assume that we have already obtained the input-output relation for 30 4.1. Eﬀective SNR at the Destination Node for the Multi-Hop Case the soft estimate x ˆi [k] formed by the ith relay, according to x ˆi [k] = βi qm−i+1 [k] + zi [k] (4.2) x ˆi [k] = β˜i qm−i+1 [k] + zi [k] (4.3) (i = 1) or (1 < i < m), cf. (3.1) and (3.2) for the dual-hop case. In particular, assume that we have calculated the gain factor βi or β˜i and the variance σi2 of the noise term zi [k]. In the sequel, we focus on the case i > 1 for simplicity. For i = 1, parameter β˜i needs to be replaced by βi in the following equations. • For the integrator output of the (i + 1)th relay (i < m − 1) or the destination node (i = m−1), we obtain the relation x ˆi+1 [k] = βi+1 x ˆi [k]ˆ xi [k − 1] + zi+1 [k], where βi+1 := Nf Ai+1 αi+1 εi+1 , εi+1 := Eg Ti 0 (4.4) 2 wrx,i+1 (t)dt, and zi+1 [k] represents a zero-mean noise term with components deﬁned similar to ˆi [k], A1 by Ai+1 , α1 by αi+1 , (3.3–3.5), with q3 [k] being replaced by x n1 (t) by ni+1 (t), and wrx,1 (t) by wrx,i+1 (t). The variance of the noise term zi+1 [k] is calculated2 as σz2 i+1 ≈ (β˜i2 + σi2 )βi+1 N0 + W Ti N02 . 2 • Plugging the expression for x ˆi [k] from (4.3) into (4.4), we obtain a new 2 Derivation is described in Appendix A.1 31 4.1. Eﬀective SNR at the Destination Node for the Multi-Hop Case representation for the soft estimate x ˆi+1 [k], according to x ˆi+1 [k] = β˜i+1 qm−i [k] + zi+1 [k]. (4.5) The gain factor β˜i+1 is obtained via the recursion β˜i+1 = βi+1 β˜i2 . (4.6) Moreover, making similar approximations as in Section 3.1 and taking the correlations between the involved noise terms into account, 2 the variance σi+1 of the noise term zi+1 [k] can be calculated3 via the recursion 2 2 σi+1 2β˜i2 (σi2 + ηi ) + σi4 + βi+1 N0 (β˜i2 + σi2 ) + ≈ βi+1 W Nf Ti N02 , 2 (4.7) where 2 2 2 N0 ηi ≈ βi2 β˜i−1 σi−1 . + 2ηi−1 + βi−1 β˜i−1 2 (4.8) Here, we have used the approximation that for a small to moderate number of intermediate A&F relays the correlation between the noise terms zi+1 [k] and zi+1 [k], cf. (4.4) and (4.5), can be neglected. 3 Calculation is explained in Appendix A.2 with more details 32 4.2. Optimized Transmit Power Allocation for the Multi-Hop Case At the end of the recursion (i = m−1), we thus obtain the following expression for the eﬀective SNR at the destination node: SNR = ≈ 2 β˜m 2 σm 2 )2 (βm β˜m−1 4 2 2 2β ˜2 ˜2 (σ 2 βm m−1 m−1 + ηm−1 ) + σm−1 + βm N0 (βm−1 + σm−1 ) + W Nf Ti N02 2 · (4.9) Based on the above results, in the following, we will propose an algorithm for obtaining a suboptimal, yet eﬃcient solution for the transmit power allocation factors αi for the source node (i = 1) and the relay nodes (2 ≤ i ≤ m), with the goal to maximize the eﬀective SNR at the destination node. 4.2 Optimized Transmit Power Allocation for the Multi-Hop Case A closed-form solution for the optimal transmit power allocation factors αi (1 ≤ i ≤ m) does not seem feasible in the multi-hop case (m > 2) because of the involved expression for the SNR (4.9). Thus, in the following, we present a heuristic algorithm for ﬁnding near-optimal values for the transmit power allocation factors αi . The proposed algorithm is based on the optimization presented in Section 3.2 for the dual-hop case. Essentially, we decompose the m–hop transmission into m − 1 two-hop transmissions. In particular, we ﬁrst assume that there are only three nodes, namely the source node, the last relay (i = m), and the destination node and ﬁnd the corresponding 33 4.2. Optimized Transmit Power Allocation for the Multi-Hop Case optimal transmit power factor for the source-relay link. In the next step, the transmit power allocated to the source node in the previous step is regarded as the overall available transmit power, and the optimum power allocation between the source node and the (m−2)th relay (Node m − 1) is calculated assuming that the (m−1)th relay (Node m) is the destination. We apply this method repeatedly until we ﬁnd the overall power allocation for the source node and subsequently calculate the powers allocated to each relay in a similar manner. In the following, let Ai,j denote the gain factor associated with the link from node i to node j, cf. (2.5) and Fig. 4.1, where i = 1 represents the source node and j = m+1 the destination node. Assuming we have only the source (Node 1), the last relay (Node m), and the destination (Node m + 1), 1 , is given by according to (3.24), the power factor for the source, ζm 1 1 = ζm 1+ A1,m 3Am,m+1 . (4.10) The fraction of power allocated to the (m − 1)th relay (Node m) is given 1 . However, we cannot calculate the associated power allocation by 1 − ζm factor αm since the received power of the relay is not known yet. Instead, 1 to be the available power and take in the next step, we consider Nf Eg ζm into account that the received signal at the (m − 1)th relay (Node m) actually originates from the (m − 1)th relay. To ﬁnd the corresponding power allocation, we consider only the source, relay i = m − 1, and relay i = m, with the latter playing the role of the destination. Therefore, the fraction 1 ζ1 of power allocated to the source is now given by ζm m−1 . This procedure is 34 4.2. Optimized Transmit Power Allocation for the Multi-Hop Case Figure 4.1: (a) Optimizing α1∗ using a recursive approach. (b) Optimization of α2∗ for known α1∗ . repeated until we arrive at ﬁrst relay (i = 2) , which leads to the following ﬁnal power allocation for the source: m α1∗ ζi1 = (4.11) i=2 At each step we will have the following expression for the fraction of 35 4.2. Optimized Transmit Power Allocation for the Multi-Hop Case power allocated to each link 1 ζi1 = 1 3 1+ · A1,i Ai,i+1 · (4.12) Knowing the power allocated to the source, the available power for allocation to the m − 1 relays is given by (1 − α1∗ )Nf Eg . The power allocation factor for relay 1 is constrained by ¯2 = 1 − α1∗ , α2 Erx,1 + E (4.13) ¯2 is where Erx,1 := β12 + σ12 denotes the received power of ﬁrst relay and E the remaining power to be allocated to the other m − 2 relays later on. The power allocation factor α2 of the ﬁrst relay (i = 2) is now obtained using the same recursive approach as for the source before, i.e., the ﬁrst relay plays the role of the source. Subsequently, the same steps are repeated for relays i = 3 to i = m. For the jth relay, 1 ≤ j ≤ m − 1, this leads to power allocation ∗ αj+1 = κj · Erx,j m ζij+1 , (4.14) i=j+2 with ζij = κj = Erx,j 1 . Aj,i 3Ai,i+1 κj−1 − αj∗ Erx,j−1 (4.15) 1+ = (Aj αi Erx,j−1 Nf Ej )2 + σj2 (4.16) (4.17) 36 4.2. Optimized Transmit Power Allocation for the Multi-Hop Case where the received power at ith relay, Erx,i , and normalization factor κi , are calculated recursively for i ≥ 2. For initialization, we have κ0 := 1 and Erx,0 = 1. It is easy to check that the suboptimal power allocation in (4.14) fulﬁlls the overall power constraint Nf Eg α1∗ + m ∗ i=2 αi Nf Eg Erx,i−1 = Nf Eg . Although the proposed heuristic power allocation scheme for the multihop case is suboptimal, our results in Chapter 6 will conﬁrm its near-optimal performance. Furthermore, the proposed power allocation algorithm allows for a semi-distributed implementation, where the source has to have only access to estimates of its own path loss to all relays and the path-losses between all neighboring relays. In particular, for the proposed semi-distributed implementation, the source computes the factors fj := j m i=j ζi , 1 ≤ j ≤ m, feeds back fj+1 to jth relay, computes α1∗ , and feeds back κ1 = 1−α1∗ to ﬁrst relay. The jth relay, 1 ≤ j ≤ m − 1, estimates its own received power Erx,j , computes αj∗ based on (4.14), computes κj based on (4.16), and feeds back κj to the next (j + 1)th relay. The last relay only has to estimate Erx,m−1 ∗ . and compute αm 37 Chapter 5 Performance Analysis of Dual-hop Transmission with Rake Receiver In Chapter 6, the performance of the proposed scheme is compared with a Rake receiver for dual-hop relaying. To have a fair comparison between Rake receivers and diﬀerential receivers, the power allocation factors for Rake reception should be optimized for maximization of the SNR at the destination. So here, we ﬁrst describe the model of the Rake receiver and then calculate the SNR of the destination for Rake reception at the relay and the destination, respectively. Then, based on the expression for the SNR, the power allocation factors are optimized for dual-hop transmission. 5.1 Eﬀective SNR at the Destination Node The performance of a pulse based Rake receiver has been studied in [13]. Here, we want to give a brief overview of the analytical basis for Rake receiver systems. In Rake combiners, the received signal is correlated with an 38 5.1. Eﬀective SNR at the Destination Node appropriate template signal which is designed to capture energy of the signal eﬃciently. The output of the correlator is passed through an integrationand-dump block. The output of the integration-and-dump unit is either transmitted to the next node (in A&F relay) or used for decision making (in D&F relay or destination node). The transmitted signal at the source can be described as Nf −1 ∞ s1 (t) = b[k] wtx (t − jTf − kTs ), α1 Eg (5.1) k=−∞ j=0 where b[k] ∈ {±1} denotes the information bit, α1 is the power allocation factor for the transmitter, and Eg represents the energy per pulse. wtx (t) is the transmited pulse shape with duration Tp such that Tp 0 2 (t)dt = 1, as wtx before. As described earlier, the multipath model of channel model is deﬁned based on [34] as (2.2) LC LR (r) (r) λ(r) μ,ν δ(t − Tν − τμ,ν ). hr (t) = (5.2) ν=0 μ=0 To simplify our analysis, for performance evaluation, it is assumed that the multipath arrival times are multiples of the minimum resolvable path interval, Tp . Thus, we can rewrite the channel model as L (r) λl δ(t − (l − 1)Tp ). hr (t) = (5.3) l=1 If there is no arrival at an speciﬁc time ((k−1)Tp ), the corresponding channel 39 5.1. Eﬀective SNR at the Destination Node (r) coeﬃcient gain will be zero (λk = 0). The received signal at the relay is the convolution of the transmitted signal with the CIR and described by Nf −1 ∞ r1 (t) = (1) b[k]wrx (t − jTf − kTs ) + nR (t), α1 A1 Eg (5.4) k=−∞ j=0 where nR (t) is zero-mean white Gaussian noise with variance Ω2n and A1 is (1) the channel’s relative gain which is given in (2.5) and wrx (t) is deﬁned as [13] L (1) (t) wrx (1) λl wtx (t − (l − 1)Tp ). = (5.5) l=1 Assuming a Rake receiver at the relay, a template signal is used for capturing and combining the energy of the received signal. The template signal is deﬁned as Nf −1 ϑ(t − jTf − kTs ), stemp,k (t) = (5.6) j=0 where L υl wtx (t − (l − 1)Tp ), ϑ(t) = (5.7) l=1 and υ = [υ1 , . . . , υl ] are the weights for Rake combining. An S-Rake scheme is used in our analysis where the Ls most signiﬁcant multi-path components are captured, causing only these components to be non-zero. There are diﬀerent schemes for combining the signal components available in the literature, e.g. Equal Gain Combining (EGC) and Maximum Ratio Combining 40 5.1. Eﬀective SNR at the Destination Node (MRC) [38]. In EGC all the combining weights are equal, and they have the same weight in contributing to the ﬁnal energy of the signal. However, in MRC we choose the Rake combining weights equal to the corresponding (i) multipath coeﬃcients i.e., for the ith realization of the channel, υl = λl . At the Rake receiver, the correlation between the received signal and the template signal is calculated. The output of the Rake receiver can be obtained as Nf −1 y1 [k] = j=0 kTs +jTf +Ti r(t)stemp,k (t)dt. kTs +jTf (5.8) Replacing (5.4) in (5.8), y1 [k] can be rewritten as follows Nf −1 y1 [k] = j=0 kTs +jTf +Ti kTs +jTf ⎛ ⎝ α1 A1 Eg ×⎝ Nf −1 k=−∞ n=0 ⎛ Nf −1 L ∞ ⎞ (1) b[k]wrx (t − nTf − kTs ) + nR (t)⎠ ⎞ υl wtx (t − (l − 1)Tp − iTf )⎠ dt i=0 l=1 Nf −1 = α1 A1 Eg ⎛ ×⎝ j=0 ∞ kTs kTs +jTf +Ti kTs +jTf Nf −1 L L υl wtx (t − (l − 1)Tp − iTf ) i=(k−1)Ts l=1 ⎞ ⎠ λ(1) r wtx (t − (r − 1)Tp − jTf − kTs ) + nR (t) dt. k=−∞ j=0 r=1 (5.9) Knowing that the duration of wtx (t) is Tp , (5.9) can be simpliﬁed to L y1 [k] = b[k] α1 A1 E1 (1) λl υl + a ˆ + nR,1 [k], (5.10) l=1 41 5.1. Eﬀective SNR at the Destination Node Ti 2 2 0 wtx (t)dt) . The noise Nf −1 kTs +jTf +Ti nR (t)stemp,k (t)dt. j=0 kTs +jTf where a ˆ is an IFI term4 and E1 = Eg (Nf nR,1 [k], is given by nR,1 [k] := term, Here we assume that LTp < Tf so there is no IFI (the assumption was made for the diﬀerential receiver). Therefore (5.10) can be rewritten as L y1 [k] = b[k] α1 A1 E1 (1) λl υl + nR,1 [k], (5.11) l=1 where nR,1 [k] is AWGN with zero mean. The variance of nR,1 [k] can be calculated as ⎛ Nf −1 Ω21 = E{n2R,1 [k]} = E ⎛ =E ⎝ ⎝ j=0 Nf −1 j=0 kTs +jTf +Ti kTs +jTf kTs +jTf +Ti (k+1)Nf −1 L kTs +jTf j=kNf ⎞2 nR (t)stemp,k (t)dt⎠ ⎞2 nR (t)υl wtx (t − jTf − (l − 1)Tp )dt⎠ l=1 L = Nf Ω2n υl2 . (5.12) l=1 Relay retransmits y1 [k] toward the destination. The transmitted signal from the relay can be expressed in a similar way as (5.1) as ∞ s2 (t) = Nf −1 y1 [k] wtx (t − jTf − kTs ). α2 Eg (5.13) k=−∞ j=0 The received signal at the second Rake receiver is the convolution of s2 (t) with the CIR of the second link plus additive white Gaussian noise, nR (t), with the same statistical properties as nR (t). The output of second 4 Performance analysis for the general case, including the case LTp > Ts , can be found in [39]. 42 5.1. Eﬀective SNR at the Destination Node Rake combiner can be derived along the same lines as (5.8)-(5.11). Finally, at the destination, we have L (2) y2 [k] = y1 [k] α2 A2 E2 λl υl + nR,2 [k], l=1 L L (1) λl υl · = b[k] α2 A2 E2 α1 A1 E1 l=1 (2) λl υl l=1 L (2) + nR,1 [k] α1 A1 E1 λl υl + nR,2 [k], (5.14) l=1 where E2 = Eg (Nf Ti 0 2 (t)dt)2 and n wtx R,2 [k] is deﬁned as follows Nf −1 nR,2 [k] := j=0 kTs +jTf +Ti kTs +jTf nR (t)stemp,k (t)dt. (5.15) The variance of nR,2 [k] can be calculated with the same approach used for obtaining (5.12), which yeilds Ω22 = E{n2R,2 [k]} = Nf Ω2n L 2 l=1 υl . The total noise at destination, nR,3 [k] is the sum of two noise terms, L nR,3 [k] = nR,1 [k] α1 A1 E1 (2) λl υl + nR,2 [k]. (5.16) l=1 The variance of the noise at the destination can be calculated as follows Ω23 = E{n23 [k]} = α2 A2 E2 Ω22 + Ω21 . (5.17) Splitting (5.14) into a signal part and a noise part, the SNR at the destina- 43 5.2. Optimized Transmit Power Allocation tion can be derived as SNR = 5.2 α2 A2 E2 · α1 A1 E1 ( N0 2 L 2 l=1 υl 1+ (2) (1) L L 2 2 l=1 λl υl ) ( l=1 λl υl ) · (2) 2 α2 A2 E2 ( L λ υ ) l l=1 l (5.18) Optimized Transmit Power Allocation Based on (5.18), we can optimize the power allocation factors for the sourcerelay and relay-destination links. Here, the eﬀective SNR at the destination node is optimized, under the constraint of keeping the total transmit power ﬁxed. It can be shown that the total transmit power constraint can be expressed as ⎛ α1 + α2 ⎝α1 A1 E1 L 2 (1) λl υl l=1 N0 + 2 L ⎞ υl2 ⎠ = 1. (5.19) l=1 To simplify our analysis, three variables are introduced as follows L K1 := A1 E1 2 (1) λl υl , l=1 L K2 := A2 E2 2 (2) λl υl , l=1 K3 := N0 2 L υl2 . (5.20) l=1 Using these variables, the power constraint can be re-written as α1 + α2 (α1 K1 + K3 ) = 1, (5.21) 44 5.2. Optimized Transmit Power Allocation and the SNR simpliﬁes to SNR = α1 α2 K1 K2 · K3 (1 + α2 K2 ) (5.22) Based on (5.21) and (5.22), a Lagrange problem can be formulated. Using Lagrange multipliers, power allocation factors, α1 and α2 are optimized for maximization of the SNR at destination. Based on this approach, the result for the optimal power factor, α1,opt , is as follows α1,opt = (K1 + K3 )(K2 + K3 ) − (K2 + K3 ) . K1 − K2 (5.23) For high SNR, K3 is negligible compared to K1 and K2 because K3 is a noise term and K1 and K2 are signal terms. So, assuming high SNR, (5.23) yields α1,opt ≈ K1 K2 − 1 K1 K2 − 1 1 = 1+ K1 K2 , (5.24) and α2 can be calculated from (5.19) as α2,opt ≈ 1 − α1,opt . α1,opt K1 + K3 If the CIRs of the ﬁrst and the second links are similar, (5.25) K1 K2 ≈ A1 A2 So the optimal power allocation factor for the source is α1,opt ≈ follows. 1 1+ A1 A2 , which only depends on the link gains A1 and A2 (similar to the DF case cf. (3.24)). Since the comparison between Rake and DF receivers is made for the 45 5.2. Optimized Transmit Power Allocation dual-hop case in Chapter 6, we only consider the power allocation for this case here. However, extending our analysis to the multi-hop case is possible by following the same steps as in Chapter 4. 46 Chapter 6 Numerical Performance Results In the following, we present numerical performance results, which illustrate the excellent performance of our proposed scheme and corroborate our analysis in Chapters 3 and 4. We start with considering the dual-hop case without channel coding. Afterwards, we present results for the multi-hop case and consider the beneﬁts of an outer FEC scheme. In the following, the information symbols q1 [k] ∈ {±1} are transmitted in blocks of 1000 symbols. The CIR is assumed to remain static for the duration of an entire block (70 μs). As an example, we focus on channel model CM1 in the sequel and assume a path-loss exponent of p = 3 unless stated otherwise. One frame is used for transmitting a single information symbol (Nf = 1), and the frame length is chosen such that the guard interval between subsequent pulses is larger than the root-mean-square (rms) delay spread of the channel (Tf = 70 ns), so as to circumvent ISI eﬀects. For the transmitted pulse wtx (t), we employ the widely used second derivative of a 47 Chapter 6. Numerical Performance Results 0.3 0.25 0.2 0.15 Amplitude 0.1 0.05 0 −0.05 −0.1 −0.15 0 1 2 3 t[s] 4 5 6 7 −10 x 10 Figure 6.1: Transmitted pulse, wtx (t), in time domain. Gaussian pulse, i.e., 2 wtx (t) = [1 − 4π(t − vp )/vm ] exp[−2π((t − vp )/vm )2 ], where vp = 0.35 ns and vm = 0.2877 ns (Tp = 0.7 ns). Fig 6.1 shows the employed transmitted pulse in time domain. The integration time Ti has been optimized such that on average the maximum eﬀective SNR at the integrator output is obtained (Ti = 5.25 ns). Moreover, the bandwidth W of the bandpass ﬁlter is optimized such that the maximum received SNR is obtained (W = 5 GHz). Fig. 6.2 depicts the frequency response of the transmitted pulse. 48 Chapter 6. Numerical Performance Results 20 10 0 −10 Amplitude [dB] −20 −30 −40 −50 −60 −70 −80 0 2 4 6 8 f [Hz] 10 9 x 10 Figure 6.2: Transmitted pulse in frequency domain, 20 log(|wtx (f )|). Fig. 6.3 shows the eﬀective SNR at the destination node for the doublediﬀerential A&F-relaying scheme for Eg /N0 = 9 dB, considering three different positions of the relay (ρ := d1 /dS−D = {0.2, 0.4, 0.8}), where d1 and dS−D denote the lengths of the source-relay and the source-destination link, respectively. Analytical results based on (3.20) are represented by lines, and corresponding simulation results are represented by markers. Moreover, the cross signs indicate the optimum value for the transmit power allocation factor α1 for the source node, which was found based on (3.24). As can be seen, the analytical and the simulation results ﬁt well, and considering the fact that our (approximate) formula (3.24) was derived for high SNR values 49 Chapter 6. Numerical Performance Results 9 ρ = 0.8 (Simulation) ρ = 0.8 (Analysis) ρ = 0.4 (Simulation) ρ = 0.4 (Analysis) ρ = 0.2 (Simulation) ρ = 0.2 (Analysis) α1,opt 8 7 SNR [dB] 6 5 4 3 2 1 0 0 0.2 0.4 α1 0.6 0.8 1 Figure 6.3: SNR at the destination versus α1 for three diﬀerent values for the source-relay distances and the optimal value of α1 for each distance. while the SNR considered in Fig. 6.3 is rather moderate, the optimal power allocation factor α1,opt obtained with (3.24) oﬀers a remarkable accuracy. Fig. 6.4 illustrates the performance of the proposed non-coherent A&F relaying scheme with double-diﬀerential encoding at the source node, obtained by means of Monte-Carlo simulations over a large number of independent CIR realizations for ρ = 0.2. The proposed scheme is compared with (i) direct transmission from the source to the destination node, (ii) non-coherent A&F relaying with single diﬀerential encoding at the source node, (iii) coherent A&F and D&F relaying with S-Rake combining at the relay and at the destination node (using L = 3, 5 Rake-ﬁngers), and (iv) 50 Chapter 6. Numerical Performance Results non-coherent D&F relaying scheme with double-diﬀerential encoding at the source node. The transmit power allocation factors α1 and α2 were optimized based on (3.24) for the proposed A&F relaying scheme with doublediﬀerential encoding and based on (5.24) for the A&F relaying scheme with Rake reception. For the other relaying schemes, power allocation factors are optimized based on simulations. As can be seen, all considered schemes oﬀer substantial performance improvements over direct transmission. For suﬃciently high SNR values, our proposed scheme outperforms the coherent A&F and D&F relaying schemes, unless a relatively large number of Rake ﬁngers is employed (e.g., L ≥ 5). The double-diﬀerential D&F relaying scheme oﬀers a small performance advantage of about 0.5 dB compared to the A&F version, at the expense of an increased relay complexity. A&F relaying with straightforward single-diﬀerential encoding at the source, ampliﬁcation at the relay, and diﬀerential decoding at the destination suﬀers from a signiﬁcant loss in performance compared to the proposed scheme. We have doubled the integration time Ti at the receiver for A&F relaying with single-diﬀerential encoding since the eﬀective CIR seen by the receiver is the convolution of the CIRs of the source-relay and the relay-destination channels. We note that the performance of A&F relaying with single-diﬀerential encoding could be improved by increasing the frame duration Tf at the expense of a loss in data rate. Now, we consider the case when an outer FEC scheme is used. The FEC block is placed before the diﬀerential encoder and includes a convolutional encoder and a bit interleaver. The data is encoded with a 64-state, rate 1/2 binary convolutional code based on octal generators [133, 171] [37]. Fig. 6.5 51 Chapter 6. Numerical Performance Results 0 10 −1 10 −2 BER 10 −3 A&F, Double Differential D&F, Double Differential A&F, Single Differential A&F, S−Rake, L=3 D&F, S−Rake, L=3 A&F, S−Rake, L=5 D&F, S−Rake, L=5 Direct Transmission 10 −4 10 −5 10 −6 10 4 5 6 7 8 9 Eg /N0[dB] 10 11 12 Figure 6.4: BER at the destination versus Eg /N0 for double diﬀerential transmission and Rake receivers, respectively. A&F and D&F protocols. shows the performance of the proposed non-coherent dual-hop A&F relaying scheme with double-diﬀerential encoding at the source node if an outer FEC scheme and hard input Viterbi decoding are used. As can be seen, using FEC improves the performance of all schemes signiﬁcantly, especially that of the proposed scheme such that we can achieve BER in order of 10−5 for less than 8 dB. Similar to Fig. 6.4, we observe also in Fig. 6.5 that the performance of the Rake receiver is better than that of double-diﬀerentially transmission with A&F relaying only for lower SNRs. Fig. 6.6 shows the performance of the proposed scheme if a soft input 52 Chapter 6. Numerical Performance Results 0 10 −1 10 −2 BER 10 −3 10 A&F Double Differantial D&F Double Differential A&F Single Differential A&F S−Rake, L=3 D&F S−Rake, L=3 Direct Transmission −4 10 −5 10 0 1 2 3 4 Eg /No[dB] 5 6 7 8 Figure 6.5: BER at destination versus Eg /N0 for double diﬀerential encoding and Rake receivers with both A&F and D&F protocols in coded case (hard input Viterbi decoding). Viterbi decoder is used at the destination. As can be seen, using soft decision decoding, the performance of the proposed non-coherent dual-hop A&F relaying scheme with double-diﬀerential encoding at the source node improves by more than 1 dB at the BER of 10−5 . Also the performance of other schemes improves signiﬁcantly compared to hard decoding. Similar to Fig. 6.5, the slope of the curves in Fig. 6.6 suggest that the proposed scheme out performs coherent transmission with Rake receivers at high enough SNR. We note that Viterby decoding with squared Euclidean distance branch metric 53 Chapter 6. Numerical Performance Results 0 10 −1 10 BER −2 10 −3 10 A&F Double Differential D&F Double Differential A&F Single Differential A&F S−Rake, L=3 D&F S−Rake, L=3 Direct Transmission −4 10 0 1 2 3 4 5 Eg /N0[dB] 6 7 8 Figure 6.6: BER at destination versus Eg /N0 for double diﬀerential encoding and Rake receivers with both A&F and D&F protocols in coded case (soft input Viterbi decoding). is not necessary optimal for the proposed double-diﬀerential encoding since that the noise at the destination is non-Gaussian. As can be seen in Fig. 6.6, for the case of double diﬀerential encoding, the A&F scheme has a signiﬁcantly better performance than the D&F scheme in the presence of soft input Viterbi decoding. This occurs because of the unreliable decision at the relay in D&F relaying which make the squared Euclidean distance metric staying suboptimal. Fig. 6.7 depicts the BER versus ρ for a dual-hop IR-UWB system with 54 Chapter 6. Numerical Performance Results 0.012 0.01 BER 0.008 0.006 0.004 0.002 0 0.1 0.2 0.3 0.4 ρ 0.5 0.6 0.7 0.8 0.9 Figure 6.7: BER vs. ρ. double diﬀerential encoding at the transmitter for a ﬁxed Eg /N0 (Eg /N0 = 9.7 dB). The power allocation factors are optimized based on (3.24). As can be seen, the optimum position for the A&F relay is the mid-point between the source and the destination (similar to the conventional single diﬀerential case [27]). Another result which is shown in Fig. 6.7, is that the performance of dual-hop system is robust to the position of the relay around the optimal point if the optimal power allocation factors are used. Next, we turn to the case of multi-hop relaying. Fig. 6.8 shows the effective SNR at the destination node versus Eg N0 for three-hop (m = 3) and four-hop (m = 4) transmission. In particular, the simulated eﬀective SNR 55 Chapter 6. Numerical Performance Results Eﬀective SNR at the Destination [dB] 30 25 20 15 3−hop case 1 , Simulation 3−hop case 2 , Simulation 4−hop case 1 , Simulation 4−hop case 2 , Simulation Analysis 10 5 0 5 10 15 20 25 30 Eg /N0[dB] Figure 6.8: Eﬀective SNR at the destination node for the three-hop and the four-hop case. values at the destination node (markers) and the corresponding analytical values are compared for two diﬀerent arrangements of relays. The arrangement of the relays in case 1 for m = 3 and m = 4 are 0.5 − 0.4 − 0.1 and 0.25 − 0.25 − 0.25 − 0.25 (normalized with respect to the length of S-D link), respectively. These arrangements in case 2 are 0.6 − 0.3 − 0.1 and 0.1 − 0.3 − 0.2 − 0.4 for three-hop and four-hop relaying, respectively. As Fig. 6.8 shows, the SNR calculated based on (4.9) and the simulated eﬀective SNR are in excellent agreement for both three-hop and four-hop transmission. 56 Chapter 6. Numerical Performance Results Table 6.1 and Table 6.2 show the (near-) optimum transmit power allocation factors αi obtained by means of the algorithm proposed in Section 4.2 for three-hop (m = 3) and four-hop (m = 4) relaying, respectively, for diﬀerent relay positions. Eg /N0 is 12.8 dB for both tables, and we chose p = 3 and p = 4 for the three-hop case and the four-hop case, respectively. The resulting values are compared with those obtained from an exhaustive search. As can be observed, the results obtained with (4.14) are close to optimal. Table 6.1: Transmit power allocation factors αi∗ for the three-hop case (m = 3); the ﬁrst column indicates the relative positions of the relays (normalized w.r.t. the source-destination link length). Relative position of relays 0.1 − 0.1 − 0.8 0.1 − 0.2 − 0.7 0.1 − 0.6 − 0.3 0.33 − 0.33 − 0.33 0.4 − 0.4 − 0.2 Type of optimization Proposed formula Exhaustive search Equal power Proposed formula Exhaustive search Equal power Proposed formula Exhaustive search Equal power Proposed formula Exhaustive search Equal power Proposed formula Exhaustive search Equal power α1∗ 0.11 0.1 0.33 0.12 0.11 0.33 0.09 0.08 0.33 0.52 0.5 0.33 0.59 0.57 0.33 α2∗ εrx1 0.06 0.05 0.33 0.18 0.15 0.33 0.75 0.8 0.33 0.30 0.31 0.33 0.34 0.37 0.33 SNR 14.95 14.99 10.36 15.52 15.59 12.53 13.73 13.69 10.54 14.20 14.24 13.63 12.90 12.93 11.40 57 Chapter 6. Numerical Performance Results Table 6.2: Near-optimum transmit power allocation factors αi∗ for the fourhop case (m = 4); the ﬁrst column indicates the relative positions of the relays (normalized w.r.t. the source-destination link length). Relative position of relays 0.25 − 0.25 − 0.25 − 0.25 0.1 − 0.3 − 0.2 − 0.4 Type of optimization Proposed formula Exhaustive search Equal power Proposed formula Exhaustive search Equal power α1∗ 0.52 0.52 0.25 0.11 0.10 0.25 α2∗ εrx1 0.27 0.26 0.25 0.52 0.52 0.25 α3∗ εrx2 0.14 0.12 0.25 0.11 0.11 0.25 58 SNR 15.38 15.40 13.95 17.48 17.48 16.00 Chapter 7 Conclusions and Future Work 7.1 Conclusions In this thesis, we proposed a double-diﬀerential encoding scheme for two-hop A&F relaying in IR-UWB systems. In conventional A&F relaying systems, the length of the eﬀective overall CIR increases in each hop. As a result, the amount of ISI impairing the received signal increases in each hop leading to a poor performance. By using double-diﬀerential encoding at the transmitter, the amount of ISI can be limited to the same level as for direct transmission, without any need to increase guard intervals between transmitted pulses. This means that with the proposed scheme, the overall system beneﬁts from path-loss gains without any loss in data rate. After a thorough performance analysis of the proposed scheme, a closed-form solution for the optimum power allocation to the source node and the relay was derived. Simulation results showed that the proposed A&F-relaying scheme oﬀers an excellent performance and can even compete with coherent A&F and D&F-relaying schemes that are based on S-Rake combining. Subsequently, we extended 59 7.2. Future Work the presented framework to the case of multi-hop relaying IR-UWB systems. A multiple-diﬀerential encoding at the source node and a single-diﬀerential decoding at each relay and at the destination node was proposed for the multi-hop case. Similar to the dual-hop case, by this means, the ISI is limited in each hop without requiring an extended guard interval. Finally, a recursive formula for the SNR at the destination node and a sub-optimal approach for ﬁnding the power allocation factors was presented for the multihop case. Simulation results showed the excellent performance of the nearoptimum power allocation factors in the more challenging case of multi-hop relaying as well as in the dual-hop case. 7.2 Future Work There are several possible extensions of the work presented in this thesis and these are listed below. • In this thesis, we have assumed that only one relay participates in each hop. Considering more than one relay e.g., two relays in each hop is a topic for future work. • A second area for further development is assuming multiple antennas at the transmitter and receivers of all nodes and employing space-time codes to improve the reliability of data transmission. • The eﬀective noise at the destination is in general non-Gaussian and thus, ﬁnding a closed-form expression for BER seems diﬃcult. Therefore, we have focused on the eﬀective SNR at the destination node as 60 7.2. Future Work performance measure in this thesis. Finding the distribution of the noise at the destination is an interesting area for future research. 61 Bibliography [1] D. Porcino and W. Hirt. Ultra-wideband radio technology: Potential and challenges ahead. IEEE Commun. Mag., 41(7):66–74, July 2003. [2] S. Roy, J. R. Foerster, V. S. Somayazulu, and D. G. Leeper. Ultrawideband radio design: The promise of high-speed, short-range wireless connectivity. Proc. IEEE, 2(92):295–311, Feb. 2004. [3] Liuqing Yang and G.B. Giannakis. Ultra-wideband communications: An idea whose time has come. IEEE Signal Processing Mag., 21(6):26– 54, Nov. 2004. [4] R. A. Scholtz, D. M. Pozar, and W. Namgoong. Ultra-wideband radio. EURASIP J. Applied Signal Processing, 2005(3):252–272, 2005. [5] Federal Communications Commission (FCC). Revision of part 15 of the commission’s rules regarding ultra-wideband transmission systems. First Report and Order, ET Docket 98–153, FCC 02–48; Adopted: February 2002; Released: March 2003. [6] Maria-Gabriella Di Benedetto, Thomas Kaiser, Andreas F. Molisch, Ian Oppermann, Christian Politano, and Domenico Porcino. UWB Com- 62 Bibliography munication Systems – A Comprehensive Overview. Hindawi Publishing Corporation, 2006. [7] ECMA. Standard ECMA-368: High rate ultra wideband PHY and MAC standard. First Edition, Dec. 2005. [8] R. Fisher, R. Kohno, M. McLaughlin, and M. Welbourn. DS-UWB physical layer submission to IEEE 802.15 Task Group 3a. Document Number P802.15-04/0137r4, IEEE P802.15, Jan. 2005. [9] IEEE P802.15–TG4a. Part 15.4: Wireless medium access control (MAC) and physical layer (PHY) speciﬁcations for low-rate wireless personal area networks (LR-WPANs): Amendment to add alternate PHY. Jan. 2007. [10] M. Win and R. Scholtz. Impulse radio: how it works. IEEE Commun. Lett., 2(1):36–38, Jan. 1998. [11] Time Domain Corporation. Comments of time domain corporation, docket 98-154. in the matter of revision of part 15 of the FCC’s rules regarding ultra-wideband transmission systems. 1998. [12] M.Z. Win and R.A. Scholtz. On the energy capture of ultrawide bandwidth signals in dense multipath environments. IEEE Commun. Lett., 2(9):245–247, Sep 1998. [13] A. Rajeswaran, V.S. Somayazulu, and J.R. Foerster. Rake performance for a pulse based UWB system in a realistic UWB indoor channel. In 63 Bibliography Proc. IEEE Int. Conf. Commun. (ICC), volume 4, pages 2879–2883, May 2003. [14] Ning He and C. Tepedelenlioglu. Performance analysis of non-coherent UWB receivers at diﬀerent synchronization levels. IEEE Trans. Wireless Commun., 5(6):1266–1273, June 2006. [15] A.A. D’Amico and L. Taponecco. A diﬀerential receiver for UWB systems. IEEE Trans. Wireless Commun., 5(7):1601–1605, July 2006. [16] D. Cassioli, M.Z. Win, F. Vatalaro, and A.F. Molisch. Performance of low-complexity RAKE reception in a realistic UWB channel. In Proc. IEEE Int. Conf. Commun. (ICC), volume 2, pages 763–767, May 2002. [17] Shiwei Zhao, Huaping Liu, and Zhi Tian. A decision-feedback autocorrelation receiver for pulsed ultra-wideband systems. In Proc. IEEE Radio and Wireless Conf., pages 251–254, Sept. 2004. [18] J. Romme and G. Durisi. Transmit reference impulse radio systems using weighted correlation. In Proc. Int. Workshop on Ultra Wideband Systems, pages 141–145, May 2004. [19] Nan Guo and R.C. Qiu. Improved autocorrelation demodulation receivers based on multiple-symbol detection for UWB communications. IEEE Trans. on Wireless Commun., 5(8):2026–2031, Aug. 2006. [20] W.P. Siriwongpairat, Weifeng Su, Zhu Han, and K.J.R. Liu. Employing cooperative diversity for performance enhancement in UWB communi- 64 Bibliography cation systems. In Proc. IEEE Wireless Commun. Networking Conf. (WCNC), volume 4, pages 1854–1859, April 2006. [21] J. Kan, J. Mietzner, C. Snow, and R. Schober. Enhancement of the ECMA-368 UWB system by means of compatible relaying techniques. In Proc. IEEE Int. Conf. on Ultra-Wideband (ICUWB), volume 2, pages 63–67, Sept. 2008. [22] A. Sendonaris, E. Erkip, and B. Aazhang. User cooperation diversityPart I. System description. IEEE Trans. on Commun., 51(11):1927– 1938, Nov. 2003. [23] A. Sendonaris, E. Erkip, and B. Aazhang. User cooperation diversityPart II. Implementation aspects and performance analysis. IEEE Trans. on Commun., 51(11):1939–1948, Nov. 2003. [24] J. Nicholas Laneman. Cooperative Diversity in Wireless Networks: Algorithms and Architectures. PhD thesis, Massachusetts Institute of Technology, September 2002. [25] J.N. Laneman, D. N. C. Tse, and G.W. Wornell. Cooperative diversity in wireless networks: Eﬃcient protocols and outage behavior. IEEE Trans. on Inf. Theory, 50(12):3062–3080, Dec. 2004. [26] J.N. Laneman and G.W. Wornell. Energy-eﬃcient antenna sharing and relaying for wireless networks. In IEEE Wireless Communications and Networking Conference, (WCNC), volume 1, pages 7–12 vol.1, 2000. 65 Bibliography [27] Q. Zhao and H. Li. Diﬀerential modulation for cooperative wireless systems. IEEE Trans. on Signal Processing, 55(7):2276–2283, May 2007. [28] A. Sendonaris, E. Erkip, and B. Aazhang. End-to-end performance of transmission systems with relays over Rayleigh-fading channels. IEEE Trans. on Commun., 6(2):1126–1131, Nov. 2003. [29] M.O. Hasna and M.-S. Alouini. Harmonic mean and end-to-end performance of transmission systems with relays. IEEE Trans. on Commun., 52(1):130–135, Jan. 2004. [30] C. Abou-Rjeily and N. Daniele J.-C. Belﬁore. On the amplify-andforward cooperative diversity with time-hopping ultra-wideband communications. IEEE Trans. Commun., 56(4):630–641, April 2008. [31] S. Zhu, K. Leung, and A. Constantinides. Distributed cooperative data relaying for diversity in impulse-based UWB ad-hoc networks. IEEE Trans. Wireless Commun., 8(8):4037–4047, August 2009. [32] K. Maichalernnukul, T. Kaiser, and F. Zheng. On the performance of coherent and noncoherent UWB detection systems using a relay with multiple antennas. IEEE Trans. Wireless Commun., 8(7):3407–3414, July 2009. [33] M.K. Simon and D. Divsalar. On the implementation and performance of single and double diﬀerential detection schemes. IEEE Trans. Commun., 40(2):278–291, Feb 1992. [34] A.F. Molisch, J.R. Foerster, and M. Pendergrass. Channel models 66 Bibliography for ultrawideband personal area networks. IEEE Wireless Commun., 10(6):14–21, Dec. 2003. [35] S. Ghassemzadeh and V. Tarokh. A statistical model for in-home UWB channels. In Proc. IEEE Int. Conf. on Ultra Wideband Systems and Technology (ICUWB), pages 59–64, 2002. [36] J Kunisch and J Pamp. Measurements results and modeling aspects for the UWB radio channels. In Proc. IEEE Int. Conf. on Ultra Wideband Systems and Technology (ICUWB), pages 19–23, 2002. [37] John Proakis. Digital Communications. New York: McGraw Hill, 4 edition, 2001. [38] Marvin K. Simon and Mohamed-Slim Alouini. Digital Communication Over Fading Channels. Hoboken, New Jersey: John Wiley & Sons, Inc., 2 edition, 2005. [39] S. Gezici, H. Kobayashi, H. V. Poor, and A. F. Molisch. Performance evaluation of impulse radio UWB systems with pulse-based polarity randomization. IEEE Transactions on Signal Processing, 53:2537–2549, 2005. 67 Appendix A Calculation of the Variance of the Noise at Each Relay for Multi-hop Transmission A.1 Calculation of σz2 i+1 From (2.6), (2.8), and (4.4) we can divide zi+1 into its components as follows zi+1 [k] = zi+1,1 [k] + zi+1,2 [k] + zi+1,3 [k], (A.1) where zi+1,1 [k], zi+1,2 [k], and zi+1,1 [k] are deﬁned by zi+1,1 [k] = x ˆi+1 [k] Ai+1 αi+1 Eg Nf −1 × j=0 kTs +Ti kTs (A.2) ni+1 (t − Ts )wrx,i+1 (t − jTf − kTs )dt, 68 Appendix A.2. Calculation of σz2i+1 zi+1,2 [k] = x ˆi+1 [k − 1] Ai+1 αi+1 Eg Nf −1 kTs +Ti × kTs +Ti zi+1,3 [k] = ni+1 (t)wrx,i+1 (t − jTf − (k − 1)Ts )dt, kTs j=0 kTs (A.3) ni+1 (t)ni+1 (t − Ts )dt. (A.4) Replacing x ˆi+1 from (4.1) in (A.3), through the same lines, the variance of zi+1 [k] can be calculated. The variance of the noise terms zi+1,1 [k], zi+1,2 [k], and zi+1,3 [k] can be approximated as σz2 i+1,1 N0 N0 2 = (β˜i + σi2 )βi+1 , 2 2 W Ti N02 . = σz2 , σz2 = σz2 = 1,3 i+1,1 i+1,3 2 ≈ E{ˆ x2i+1 [k]}βi+1 σz2 i+1,2 Thus, the variance of zi+1 , σz2 i+1 W Ti N02 . 2 A.2 , is calculated as σz2 i+1 (A.5) ≈ (β˜i2 +σi2 )βi+1 N0 + Calculation of σz2i+1 Plugging x ˆν from (4.1) into (4.4), yields qm−i [k] 2 x ˆi+1 [k] = βi+1 (β˜i qm−i+1 [k]qm−i+1 [k − 1] +β˜i (qm−i+1 [k]zi [k − 1] +qm−i+1 [k − 1]zi [k]) + zi [k]zi [k − 1]) + zi+1 [k]. (A.6) 69 Appendix A.2. Calculation of σz2i+1 Comparing (A.6) with (4.5), we obtain zi+1 [k] = βi+1 (β˜i · (qm−i+1 [k]zi [k − 1] + qm−i+1 [k − 1]zi [k]) +zi [k]zi [k − 1]) + zi+1 [k]. (A.7) The ﬁrst two terms in (A.7) are not mutually independent, and we need to calculate their correlation in order to ﬁnd the variance of zi+1 [k]. From (A.7) we observe that qm−i+1 [k]zi [k] = βi (β˜i−1 (qm−i+2 [k − 1]zi−1 [k − 1] + qm−i+2 [k]zi−1 [k]) + qm−i+1 [k]zi−1 [k]zi−1 [k − 1]) + qm−i+1 [k]zi [k], (A.8) and qm−i+1 [k − 1]zi [k − 1] =βi (β˜i−1 (qm−i+2 [k − 2]zi−1 [k − 2] + qm−i+2 [k − 1]zi−1 [k − 1]) + qm−i+1 zi−1 [k − 1]zi−1 [k − 2]) + qm−i+1 [k − 1]zi [k − 1]. (A.9) Using the same approach as in Section 3.1, ηi := E{qm−i+1 [k]zi [k]qm−i+1 [k− 1]wi [k − 1]} can be calculated as follows 70 Appendix A.2. Calculation of σz2i+1 ηi = E{qm−i+1 [k]zi [k]qm−i+1 [k − 1]wi [k − 1]} 2 2 zi−1 [k − 1]} ≈ E{βi2 β˜i−1 2 + 2E{βi2 β˜i−1 qm−i+1 [k − 1]zi [k − 1]qm−i+1 [k − 2]zi [k − 2]} N0 2 βi · + β˜i−1 2 (A.10) The last term is obtained from E{qm−i+1 [k]zi [k]qm−i+1 [k −1]zi [k −1]}. Similar to the Section 3.1, a careful look at (A.2) shows qm−i+2 [k]zi,1 [k] and qm−i+2 [k −2]zi,2 [k −1] are not mutually independent. Thus similar to (3.18), we obtain E{qm−i+1 [k]zi [k]qm−i+1 [k−1]zi [k−1]} = E{qm−i+2 [k]zi,1 [k]qm−i+2 [k− 2 β N0 . 2]zi,2 [k − 1]} = β˜i−1 i 2 Finally ηi can be calculated recursively as 2 2 2 2 ηi ≈ βi2 β˜i−1 σi−1 + 2βi2 β˜i−1 ηi−1 + β˜i−1 βi N0 · 2 (A.11) Here, we have used the approximation that for a small to moderate number of intermediate A&F relays the correlation between the noise terms zi+1 [k] and zi+1 [k], cf. (4.4) and (4.5), can be neglected. According to (A.5), (A.7), and (A.11), the variance of the noise at the ith relay can also be obtained recursively as W Ti N02 2 2 σi+1 2β˜i2 (σi2 + ηi ) + σi4 + βi+1 N0 (β˜i2 + σi2 ) + . ≈ βi+1 2 (A.12) 71
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Multiple-differential encoding for multi-hop amplify-and-forward...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Multiple-differential encoding for multi-hop amplify-and-forward relaying in IR-UWB systems Hamdi, Maziyar 2010
pdf
Page Metadata
Item Metadata
Title | Multiple-differential encoding for multi-hop amplify-and-forward relaying in IR-UWB systems |
Creator |
Hamdi, Maziyar |
Publisher | University of British Columbia |
Date Issued | 2010 |
Description | In this thesis, we propose to improve the performance and coverage of impulse-radio ultra-wideband (IR-UWB) systems by means of cooperative multi-hop relaying. With regard to a simple practical realization, we focus on a non-coherent system setup in conjunction with amplify-and-forward (A&F) relaying. In particular, we propose to employ a multiple-differential encoding scheme at the source node and single differential decoding at each relay and at the destination node, respectively, so as to efficiently limit intersymbol-interference (ISI) effects at the destination node. A thorough performance analysis of the proposed scheme is provided, along with a closed-form optimization of the transmit power allocation to the source node and the relays. Simulation results illustrate the excellent performance of the proposed scheme, which is also compared to alternative coherent and non-coherent multi-hop IR-UWB schemes based on A&F relaying and decode-and-forward (D&F) relaying. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2010-03-01 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial 3.0 Unported |
DOI | 10.14288/1.0064926 |
URI | http://hdl.handle.net/2429/21267 |
Degree |
Master of Applied Science - MASc |
Program |
Electrical and Computer Engineering |
Affiliation |
Applied Science, Faculty of Electrical and Computer Engineering, Department of |
Degree Grantor | University of British Columbia |
GraduationDate | 2010-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Rights URI | http://creativecommons.org/licenses/by-nc/3.0/ |
AggregatedSourceRepository | DSpace |
Download
- Media
- 24-ubc_spring_2010_hamdi_maziyar.pdf [ 471.83kB ]
- Metadata
- JSON: 24-1.0064926.json
- JSON-LD: 24-1.0064926-ld.json
- RDF/XML (Pretty): 24-1.0064926-rdf.xml
- RDF/JSON: 24-1.0064926-rdf.json
- Turtle: 24-1.0064926-turtle.txt
- N-Triples: 24-1.0064926-rdf-ntriples.txt
- Original Record: 24-1.0064926-source.json
- Full Text
- 24-1.0064926-fulltext.txt
- Citation
- 24-1.0064926.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.24.1-0064926/manifest