Bit-Interleaved Coded Modulation for Hybrid RF/FSO Systems by Xiaohui He B.ENG., McMaster University, 2005 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 December 2008 c Xiaohui He, 2008 Abstract In this thesis, we propose a novel architecture for hybrid radio frequency (RF)/free–space optics (FSO) wireless systems. Hybrid RF/FSO systems are attractive since the RF and FSO sub–systems are affected differently by weather and fading phenomena. We give a thorough introduction to the RF and FSO technology, respectively. The state of the art of hybrid RF/FSO systems is reviewed. We show that a hybrid system robust to different weather conditions is obtained by joint bit–interleaved coded modulation (BICM) of the bit streams transmitted over the RF and FSO sub–channels. An asymptotic performance analysis reveals that a properly designed convolutional code can exploit the diversity offered by the independent sub–channels. Furthermore, we develop code design and power assignment criteria and provide an efficient code search procedure. The cut–off rate of the proposed hybrid system is also derived and compared to that of hybrid systems with perfect channel state information at the transmitter. Simulation results show that hybrid RF/FSO systems with BICM outperform previously proposed hybrid systems employing a simple repetition code and selection diversity. ii Contents Abstract .............................................. ii Contents ............................................. iii List of Figures ......................................... v List of Abbreviations and Symbols ........................ vii Acknowledgments ...................................... ix 1 Introduction ........................................ 1 1.1 Background and Motivation . . . . . . . . . . . . . . . . . . . . 1 1.2 Contributions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 1.3 Thesis Organization . . . . . . . . . . . . . . . . . . . . . . . . . 4 2 Introduction to FSO and 60-GHz RF ................... 5 2.1 FSO Fundamentals . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.1.1 Brief Description of Optical Communication Systems . . 6 2.1.2 Factors Affecting FSO . . . . . . . . . . . . . . . . . . . 7 2.1.3 FSO Fading Model . . . . . . . . . . . . . . . . . . . . . 9 2.2 60-GHz Millimeter-Wave Fundamentals . . . . . . . . . . . . . . 11 2.2.1 Attraction of 60–GHz MMW . . . . . . . . . . . . . . . . 13 2.2.2 60-GHz Band over 2.4-GHz and 5 GHz . . . . . . . . . . 14 2.2.3 Atmospheric Absorption . . . . . . . . . . . . . . . . . . 15 2.2.4 60-GHz MMW Fading Model . . . . . . . . . . . . . . . 16 iii 2.3 State-of-the-art of Hybrid RF/FSO Systems . . . . . . . . . . . 17 3 Hybrid FSO/RF System Model ........................ 21 3.1 System Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 3.2 Transmitter Structure . . . . . . . . . . . . . . . . . . . . . . . 23 3.3 Receiver Structure . . . . . . . . . . . . . . . . . . . . . . . . . 26 4 Performance Analysis and System Optimization ........... 27 4.1 Asymptotic Pairwise Error Probability . . . . . . . . . . . . . . 27 4.2 Asymptotic Codeword Error Probability . . . . . . . . . . . . . 30 4.3 Transmit Power Assignment . . . . . . . . . . . . . . . . . . . . 30 4.4 Bound on Rc . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32 5 Code and Assignment Pattern Search ................... 33 5.1 Search Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . 33 5.2 Computation Method . . . . . . . . . . . . . . . . . . . . . . . . 34 5.3 CC and AP Search Results . . . . . . . . . . . . . . . . . . . . . 35 6 Cut-off Rate Analysis ................................ 37 7 Simulation Results .................................. 40 7.1 Outage Probability . . . . . . . . . . . . . . . . . . . . . . . . . 40 7.2 Frame Error Rate (FER) . . . . . . . . . . . . . . . . . . . . . 41 8 Conclusions and Future Work ......................... 51 8.1 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 8.2 Recommendations for Future Work . . . . . . . . . . . . . . . . 52 Bibliography .......................................... 53 iv List of Figures 2.1 Schematic diagram of point-to-point FSO communication system[1]. 2.2 Sample log-normal and Gamma-Gamma fading intensity pdfs . 12 2.3 Sample Rice fading intensity pdfs . . . . . . . . . . . . . . . . . 17 3.1 SNR of FSO (γ1 ) and RF (γ2) systems vs. link distance L for transmit powers of P1 = P2 = 10 mW and different weather conditions. The values for a1 and arain have been taken from Table 2.1. All other link parameters have been chosen as in Table 3.1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25 3.2 Transmitter and receiver structure of proposed hybrid RF/FSO system. S/P: Serial–to–parallel conversion. P/S: Parallel–to– serial conversion. VD: Viterbi decoder. . . . . . . . . . . . . . . 26 7.1 Outage probability for target rate of Rtarget = 1 Gbit/s for hybrid RF/FSO system with joint BICM (proposed), independent BICM with CSI at the transmitter, and independent BICM without CSI at the transmitter, respectively. 16–QAM in RF link and L = 2.5 km. . . . . . . . . . . . . . . . . . . . . . . . . 42 7.2 Outage probability for target rate of Rtarget = 500 Mbit/s for hybrid RF/FSO system with joint BICM (proposed), independent BICM with CSI at the transmitter, and independent BICM without CSI at the transmitter, respectively. 16–QAM in RF link and L = 2.5 km. . . . . . . . . . . . . . . . . . . . . . . . . 43 v 7 7.3 FER of hybrid RF/FSO system with joint BICM (16—QAM, Rc = 1/2, q = 7, N = 6, equal power for RF and FSO sub– channels, cf. Table 5.1) and RC–SD (uncoded RF and uncoded FSO sub–channels), respectively, and uncoded FSO. L = 2.5 km and data rate of 1 Gbit/s. . . . . . . . . . . . . . . . . . . . 45 7.4 FER of hybrid RF/FSO system with joint BICM (16—QAM, Rc = 1/2, q = 7, N = 6, cf. Table 5.1) and RC–SD (uncoded RF and uncoded FSO sub–channels), respectively, and uncoded FSO. P2,max = 10 dBm, L = 2.5 km, and data rate of 1 Gbit/s. . 46 7.5 FER of hybrid RF/FSO system with joint BICM (Rc = 1/4, q = 3, N = 12, equal power for RF and FSO sub–channels, cf. Table 5.1) and RC–SD (16–QAM, Rc = 1/2, q = 3, in RF and FSO sub–channels), respectively, and coded RF (Rc = 1/2, q = 3) and coded FSO (Rc = 1/2, q = 3). L = 2.5 km and data rate of 500 Mbit/s. . . . . . . . . . . . . . . . . . . . . . . . . . 48 7.6 FER of hybrid RF/FSO system with joint BICM (4–QAM, Rc = 1/3, q = 4, N = 9, equal power for RF and FSO sub–channels, cf. Table 5.1) and RC–SD (uncoded RF sub–channel with 4– QAM, Rc = 1/2, q = 4 in FSO sub–channel), respectively, and coded FSO (Rc = 1/2, q = 4). L = 2.5 km and data rate of 1 Gbit/s. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 7.7 FER of hybrid RF/FSO system with joint BICM with different combination of constraint length, assignment pattern length and power limitation (coded RF subchannel with 16–QAM Rc = 1/2 and coded FSO subchannel Rc = 1/2 ). L = 2.5 km and data rate of 1 Gbit/s . . . . . . . . . . . . . . . . . . . . . . . . . . . 50 vi List of Abbreviations and Symbols Acronyms AP Assignment pattern AWGN Additive white Gaussian noise BICM Bit-interleaved coded modulation CC Convolutional code CEP Codeword error probability CSI Channel state information EDFA Erbium doped fiber amplifier EIRP Equivalent isotropic radiated power FSO Free-space optics GTF Generated transfer function IM/DD Intensity modulation and direct detection LD Laser diode LDPC Low density parity check LED Light-emitting diode MIMO multiple input multiple output MMW Millimeter wave OOK On-off keying pdf Probability density function PEP Pairwise error probability vii QAM Quadrature amplitude modulation RC-SD Repetition code and selection diversity RF Radio frequency S.I. Scintillation index SNR Signal-to-noise ratio VD Viterbi decoder Operators and Notations |·| Absolute value of a complex number ·⊗· Convolution δ(·) Dirac delta function E {·} Expectation det (·) Determinant tr (·) Trace [·]∗ Complex conjugate [·]T Matrix or vector transposition [·]H Matrix or vector Hermitian transposition [x]+ max(x, 0) 0m All–zero column vector of length m Im Identity matrix with dimension m × m In (·) The modified Bessel function of the first kind and n-th order viii Acknowledgments Firstly and foremost, I would like to take this chance to give my warm and grateful thanks to my supervisor, Dr. Robert Schober, for providing me with the opportunity to work in the field of free space optics and wireless communication, and nurturing my knowledge with amazing generosity and patience. I always benefited from his overall guidance, and support of my future career. I would also like to extend my thanks to the people I have met at UBC Data Communication Lab during my research: Yangwen Liang, Harry Chen, Ali Nezampour, Ehsan Bayaki, and others, for their help and encouragement. Sometimes they have shared the knowledge with me, some other times they have just provided moral support to me but, one way or the other, they have always been most helpful. Finally, I would like to thank my family. They gracefully endured these years of voluntary exile in Vancouver and gave me their unlimited support. Xiaohui He The University of British Columbia Vancouver, Canada September 2008 ix Chapter 1 Introduction The following section provides an overview of the background information and motivation for this work. We also review the related work that has been proposed by other researchers in this field. The contributions of this work are briefly summarized in the second section of this chapter, and the concluding section outlines the organization of the thesis. 1.1 Background and Motivation Free–space optics (FSO) and millimeter wave (MMW) radio frequency (RF) systems have received considerable attention for high rate (e.g. 100 Mbit/s– 10 Gbit/s) wireless communication over distances of up to a few kilometers [2]–[5]. Besides the high data rates, FSO and 60–GHz RF systems offer similar advantages such as flexibility of deployment, license–free operation, and inherent security due to high link attenuation. However, both technologies are affected quite differently by atmospheric and weather effects. FSO links suffer from extremely high attenuation in the presence of fog but are less affected by rain. In contrast, fog has practically no effect on 60–GHz systems but rain significantly increases link attenuation. Similarly, while atmospheric turbulence caused by variations in the refractive index is the main cause of small–scale fading in FSO links [6, 7], RF links are impaired by Rayleigh or Ricean fading 1 due to multipath propagation [4, 8, 9]. The complementary nature of FSO and RF channels has led to various proposals for hybrid RF/FSO systems. The combination of high–rate FSO with low–rate RF systems operating in the 2.4 GHz frequency band was proposed in [10] and experimentally validated in [11, 12]. However, ultimately the achievable data rate of this type of hybrid system is limited by the small RF bandwidths available in the 2.4 GHz band. This problem can be overcome if the RF system operates in the MMW frequency band which supports data rates similar to those of FSO [13]–[15]. If channel state information (CSI) is available at the transmitter, a simple CSI–dependent switching between the FSO and the RF system can be applied [10]. An efficient coding scheme based on low density parity check codes for hybrid RF/FSO channels with CSI at the transmitter has been proposed in [16]. To avoid the overhead necessary for CSI acquisition and to be robust against instantaneous interruptions caused by e.g. bird fly–through in the FSO system, transmitting the same data over both links has been proposed in [13]. In this scheme, the receiver detects the data streams received over the FSO and the RF links independently, and selects that data stream which is deemed more reliable on a frame–by–frame basis. This approach may be viewed as a simple repetition code across both links. However, the complimentary nature of both links is still not fully exploited. More details about the scheme in [13] can be found in Section 2.3. 1.2 Contributions In this thesis, we introduce a robust coding scheme for hybrid RF/FSO channels for the case when CSI is not available at the transmitter. From a communication theory point of view, the hybrid RF/FSO channel can be modeled as a heterogeneous block fading channel with two independent sub–channels having different fading statistics and favoring different modulation schemes. Coding for homogeneous block fading channels having identical statistical prop- 2 erties has been extensively studied in the literature [17, 18]. Bit–interleaved coded modulation (BICM) [19] was shown to achieve close–to–optimum performance in this type of channel [20] and efficient code designs have been reported [21, 22]. As will be shown in this thesis, BICM is also well suited for the heterogeneous block fading channel typical for hybrid RF/FSO systems. In particular, its convenient separation of coding and modulation facilitates the application of different modulation schemes over the RF and FSO sub– channels. The main contributions of the present research work are as follows: • For the case when CSI is not available, a joint BICM scheme is designed for hybrid FSO/RF channels. Coded sequences are transmitted simultaneously over both channels according to a bit assignment pattern and combined via standard Viterbi decoding at the receivers. For FSO channel a Gamma-Gamma fading model is adopted, while Ricean fading is assumed for RF channel fading. • Analytical expressions for the asymptotic pairwise error probability (PEP) and the asymptotic codeword error probability (CEP) are presented. In particular, union bounds for asymptotic PEP and CEP are derived. Bounds on code rates for both channels are determined. This is an extension of the work in [19]. In [19], PEP and CEP are derived when all parallel channels are identical. In our work, we extend the theory to nonidentical, heterogeneous parallel channels. • Algorithms and computational methods for searching the optimum combination of code and assignment pattern for different code constraint lengths and pattern lengths are presented. Results are also presented. • A cut-off rate analysis for the proposed joint BICM for FSO/RF channels is presented. A simulation comparison between the proposed design and independent BICM with CSI, shows that the cut-off rate of joint BICM is only slightly lower than that of independent BICM with CSI while the system complexity is greatly reduced because of the removal of CSI 3 estimation. • The simulation results for the FER show that for a given bit rate, the proposed joint BICM scheme outperforms the scheme proposed in [13], namely repetition coding with selection diversity, as well as the scheme when only FSO channel is in use. We compare the performance under weather conditions such as clear air, haze and rain. The results of our work are summarized in the following papers: • B. He and R. Schober. Bit–Interleaved Coded Modulation for Hybrid RF/FSO Systems. Accepted subject to revisions for publication in the IEEE Transactions on Wireless Communications, July 2008. • B. He and R. Schober. Bit-Interleaved Coded Modulation for Hybrid RF/FSO Systems. Submitted to IEEE International Conference on Communication (ICC) 2009,Dresden, Germany, July 2008. 1.3 Thesis Organization To explain the above findings in detail, the thesis is organized as follows. In Chapter 2, an introduction to both FSO and 60-GHz MMW channels is given. Fading models for both types of channels are provided. In Chapter 3, the system model including transmitter and receiver structure is given. We describe the algorithm for joint code and assignment pattern search in Chapter 5. Search results are also presented in this section. The analysis of the cut-off rate for joint BICM hybrid RF/FSO systems is given in Chapter 6. Simulation results are provided in Chapter 7. The results also include comparisons between our joint BICM system and previously proposed systems. Finally, conclusions and discussion of future work are provided in Chapter 8. 4 Chapter 2 Introduction to FSO and 60-GHz RF 2.1 FSO Fundamentals Free space optics (FSO) can be defined as unguided or wireless communications in the optical band emission, both coherently and incoherently. Historically, conveying information optically through the atmospheric channel was one of the most primitive forms of communication known to mankind. Classical examples include Euclid’s Optica and Hero’s Catopricia which discussed how the ancient Greeks used the sun’s reflection from metal disks to signal over vast distances. In the 1790s, Claude Chappie invented the Optical Telegraph which uses a series of semaphores-mounted towers to relay messages from one to another. In 1880, Alexander Graham Bell used a photophone to transmit telephone signals via an intensity-modulated optical beam 200 m through air to a distant receiver. The resurrection of modern optical communication can be associated with the invention of lasers in the late 1950s. The properties of laser radiation, such as monochrome, coherence, and intenseness, enable the possibility of transmitting high data rate at around 2.5 Gbps [1, 10]. Early development of unguided laser communication was directed towards 5 satellite and remote military applications. One example is the Low Earth Orbit (LEO) or the Geosynchronous (GEO) satellite-based communication system, with an optical link operating at a typical range of 45,000 km. Recently, although fiber-optic cabling is still the preferred media for long haul, high-bandwidth transport, its drawbacks such as high installation cost and the unremovable nature have opened a room for the development of commercial FSO research and deployment. Currently FSO is regarded as one of the most potent last mile technologies, and already deployed as point-to-point Giga bit connection for high-rise building. Other applications of FSO also include metropolitan area network extensions, enterprise/local area network connectivity, fiber backup, back-haul for wireless cellular networks, redundant link and disaster recovery [2, 3]. 2.1.1 Brief Description of Optical Communication Systems Figure 2.1 illustrates a conventional point-to-point FSO communication system [1]. The transmitter comprises a photoemitter, e.g., a laser diode (LD) or a light-emitting diode (LED), as light source, an erbium doped fiber amplifier (EDFA) which often has a wavelength around 980 nm and sometimes around 1550 nm,, and a telescope assembly. Input data is modulated by Intensity Modulation (IM), and carried by the emitted light from the photoemitter. The intensity, and hence the instantaneous power is proportional to the modulating current at the photoemitter. The power of the emitted light may be low and is usually amplified by an EDFA. The telescope can be designed by using either lenses or a parabolic mirror, and is used to narrow the light beam and to project it toward the receiver. At the receiver side, the transmitted light is picked up using a lens or a mirror and then focused on a direct detection (DD) photodetector, e.g., a semiconductor photodiode, a PIN (P-type, intrinsic, N-type) diode, or an avalanche 6 photodiode (APD). An electrical current is generated by photoabsorption at a rate which is proportional to the instantaneous optical power incident on the active detector surface. The intrinsically discrete nature of the current gives rise to shot noise observed in low light level detection; additional shot noise results from external sources of radiation, called background radiation, as well as from spontaneously generated charge particles, called dark current [6]. The signal is hence determined by the photocurrent. Data Data EDFA APD Receiver Telescope 1550nm Laser Path Control Figure 2.1: Schematic diagram of point-to-point FSO communication system[1]. 2.1.2 Factors Affecting FSO Several factors are associated with FSO link impairments and cause transmission power attenuation and degradation of system reliability. These factors include scintillation, obstruction, misalignment, weather effects, and window attenuation. For our research, we assume that a Line-of-Sight (LOS) link is always available. This is a practical assumption for most FSO applications. Also, we always assume perfect alignment. Hence in the following, we would focus on scintillation and weather effect. 7 Atmospheric Turbulence In hot, dry climates, turbulence effects may degrade transmission performance. As the ground heats up in the sun, the air heats up as well. Some air cells or air pockets heat up more than others. This causes changes in the index of refraction, which in turn changes the path of light propagation. This fluctuation of the index of refraction results in atmospheric turbulence, also known as scintillation, and is a major reason of signal fading. In the literature, the normalized variance of intensity, more commonly known as scintillation index (S.I.), is used to denote the statistical characteristic of the light intensity [1]: S.I. = E I2 −1 (E {I})2 (2.1) Here I denotes the light intensity. In general, a low S.I. corresponds to less scattering, which . Many statistical models for the intensity fluctuation in FSO channels have been proposed over the last two decades [7]. For weak turbulence, the lognormal distribution has been proposed as a suitable model for the intensity fluctuation. However, this model is not valid for strong turbulence. In recent research on FSO channel modeling, the Gamma-Gamma distribution was found more suitable for weak, medium and strong turbulence under various conditions [23, 7, 24]. More details will be given in Sub-section 2.1.3. Notice that the parameters which determine the fading models can be directly expressed in terms of the S.I.. It can be confidently expected that the error performance will degrade with increased S.I., as decision error events are typically correlated with low instantaneous fading intensity values. Weather Effects The FSO link quality is highly dependent on the weather conditions. FSO links suffer additional path loss under conditions such as rain, fog, and haze. The path loss due to low-visibility weather conditions is considered the major drawback in terms of reliability of FSO links and a significant barrier for 8 deployment. The path loss is a result of scattering and hence depends on the particle size and visibility in the atmosphere. For example, the impact of rain is significantly less than that of fog because the radius of rain drops (200-2000 µm) is significantly larger than the wavelength of typical FSO light source. The attenuation of the laser power can be described by Beers-Lambert Law [25] g1 [dB] a1 L1 , (2.2) where a1 denotes the weather–dependent attenuation coefficient (in dB/km) which varies depending on the weather, while L is the link distance. The attenuation coefficient a1 can be further expressed as a function of visibility V [km], the size distribution of the scattering particles q, and wavelength λ [km][26]: a1 = 3.91 λ ( )−q , V 550nm The values of q depends on V : 1.6 for high visibility (V >50 km) 1.3 for average visibility (6 km< V <50 km) q= 0.16V + 0.34 for haze visibility (1 km< V <6 km) V − 0.5 for mist visibility (0.5 km< V <1 km) 0 for fog visibility (V <0.5 km) (2.3) (2.4) For our case, we choose λ = 1550 nm as the FSO wavelength. According to 2.4, with V provided by [26], we can calculate the values of a1 for clear air, haze, light fog, moderate fog, heavy fog, moderate rain and heavy rain are (in [dB/km]) 0.43, 4.2, 20, 42.2, 125, 5.8, and 9.2, respectively, see also Table 2.1 at the end of this chapter, which we summarize the parameters we used in for our later simulation. 2.1.3 FSO Fading Model As mentioned in Section 2.1.2, we present two specific fading intensity models. Without loss of generality, all types of power density function (pdf) have been 9 normalized such that µI = E {I} = 1. Note that the distributions are only valid for I ≥ 0. Log-normal fading is the most widely accepted fading model for weak turbulence coditions. Its pdf is given by [27],[7] pI (I) = 1 (ln(I) − mx )2 exp − 2σx2 2πσx2 I (2.5) Notice that, given a specific scintillation index, the distribution parameters are set according to mx = −0.5ln(S.I. + 1) and σx2 = ln(S.I. + 1), in order to satisfy E {I} = 1. The Gamma-Gamma fading distribution, a more recently proposed FSO fading intensity model, is considered to be more suitable for fading channels with moderate to strong turbulence, as well as weak turbulence. The fading is considered to be two stochastic processes, in which small-scale atmospheric fluctuations are modulated by large-scale fluctuations; both types of atmospheric fluctuations are modeled as being gamma distributed [23],[7]. The pdf of the intensity is given by pI (I) = 2(αβ)(α+β)/2 (α+β)/2−1 I Kα−β 2 Γ(α)Γ(β) αβI , I ≥ 0, (2.6) where Γ {·} is the gamma function, Kv (·) is the vth-order modified Bessel function of the second kind, and parameters α > 0 and β > 0 can be adjusted to achieve a good agreement between pI (I) and measurement data for a wide range of turbulence conditions (weak to strong). Alternatively, assuming spherical wave propagation, α and β can be directly linked to physical parameters via [7, 23] α = β = exp 0.49χ2 (1 + 0.18d2 + 0.56χ12/5 )7/6 exp 0.51χ2 (1 + 0.69χ12/5 )−5/6 (1 + 0.9d2 + 0.62d2χ12/5 )5/6 −1 −1 , (2.7) −1 −1 , (2.8) where χ2 d 0.5Cn2 κ7/6 L11/6 , (2.9) (κD2 /4L)1/2 , (2.10) 10 and κ 2π/λ1 . (2.11) Here, λ1 and D are the wavelength and the diameter of the receiver’s aperture, respectively. Cn2 stands for the altitude–dependent index of refraction structure parameter. Here we adapt the Hufnagel-Valley model [6] for Cn2 profile model Cn2 (h) = 0.00594(v/27)2(10−5 h)10 exp(h/1000) +2.7 × 10−6 exp(−h/1500) + A · exp(−h/1000) (2.12) where h is the altitude in meters (m), v is the rms windspeed in meters per second (m/sec) and A is a nominal value of Cn2 (0) at the ground in m−2/3 . For FSO links near the ground, Cn2 can be taken approximately 1.7 × 10−14 m−2/3 during daytime and 8.4 × 10−15 m−2/3 at night. Generally, Cn2 varies from 10−13 m−2/3 for strong turbulence to 10−17 m−2/3 for weak turbulence with 10−15 m−2/3 as a typical average value [28]. The scintillation index can be calculated from α and β according to S.I. = α−1 + β −1 + (α · β)−1 (2.13) For our research we adopt the Gamma-Gamma fading model for the rest of this thesis. Figure 2.2 depicts sample pdf’s of the log-normal fading model when S.I= 0.4 and the Gamma-Gamma fading model when S.I.= 0.88 . 2.2 60-GHz Millimeter-Wave Fundamentals In response to the increasing demand for wireless bandwidth far in excess of the currently available bands at 2.4-2.5 and 5.2-5.8 GHz, in 2001, the Federal Communications Commission (FCC) allocated 7 GHz of bandwidth in the 5764 GHz band for unlicensed use. This new allocation triggered great interest in developing affordable 60-GHz radios. 11 1 Lognormal, S.I.=0.4 Gamma−Gamma, S.I.=0.88 0.9 0.8 0.7 0.6 pI (I) 0.5 0.4 0.3 0.2 0.1 0 0 0.5 1 1.5 2 2.5 3 3.5 I Figure 2.2: Sample log-normal and Gamma-Gamma fading intensity pdfs 12 4 2.2.1 Attraction of 60–GHz MMW There are a few major benefits associated with the 60-GHz band. Easy to deploy First of all, the majority of the countries have allocated 7 GHz of continuous unlicensed spectrum at varying points between 57–66 GHz. For example, besides the USA, the spectrum allocations for Australia, Japan, and Europe are 59.4 GHz-62.9 GHz, 59 GHz-66 GHz, 57-66 GHz, respectively. More importantly, the 60-GHz band is unlicensed. Commercially, this means companies can deploy universal products globally and consumers can operate equipment in this band without purchasing a license. Achievable high data rates The 60-GHz band provides 5 GHz to 7 GHz available spectrum worldwide, and each of the channels shares 2500 MHz bandwidth, this is beyond the capacities of IEEE 801.11n (520 MHz) and UWB (40 MHz). Also, due to the oxygen absorption characteristics at 60 GHz, various regulators across Asia, Europe, and the Americas allow high Equivalent Isotropic Radiated Power (EIRP) for transmissions in this band. According to Shannon’s Law [29] C = B · log2 (1 + SNR) (b/s) (2.14) where C, B and SNR denote the channel capacity, the bandwidth of the channel and the signal to noise ratio. The large bandwidth in the 60-GHz band, in addition to the high transmit power, enables data rates of up to 25,000 Mbps, which is much greater than 802.11n (1,100 Mbps) and UWB (80 Mbps) [30]. Interference-free and highly secure Oxygen absorption has a peak around 60 GHz. This introduces benefits for 60 GHz and is one of the main reasons why this band is chosen due to its resis- 13 tance to interference. This is because that a particular radio link in this band is quickly reduced to a level that will not interfere with other 60-GHz links, giving the possibility of operating several 60-GHz devices at same location. In addition, radios operating at this frequency require high power to combat oxygen absorption, and hence the signals are typically very directional. Directivity is a measure of how well an antenna focuses its energy in an intended direction. A narrow beam sizes gives immunity to interference and interception. 2.2.2 60-GHz Band over 2.4-GHz and 5 GHz In fact, 60-GHz is not the only license-free band. There are also other bands reserved internationally for the use of RF electromagnetic fields for industrial, scientific and medical (ISM) purposes other than communication [31]. The local wireless access based IEEE 802.11 standard operates in the 2.4-GHz and 5-GHz bands bandwidth of nearly 100 MHz and 150 MHz respectively [32]. In Europe, the high performance radio local area network (HIPERLAN) standard operates in the 5-GHz band [33]. The maximum net bit rates of these bands can be up to 54 Mbps for 802.11g [32] and HIPERLAN/2 [33]. However, there exists a multitude of multimedia applications calling for data rates up to a range of hundreds of megabits per second, which cannot be accommodated by current’s local wireless access standards [34]. Examples include wireless TV (150-270 Mbps), wireless video conferencing (10-100 Mbps), and wireless IEEE 1394 (up to 400 Mbps) [34]. In principle, there are two ways to improve the channel capacity: increasing the spectral efficiency and/or using more bandwidth. spectral efficiency is defined as the amount of information that can be transmitted over a given bandwidth in a specific communication system, and is measure of how efficiently a limited frequency spectrum is utilized [34]. Methods to increase spectral efficiency include using higher order modulation methods are multiple transmitters/receivers with space-time coding. However, these improvements have tradeoff in different ways. For example, using higher order modulation comes with higher transmit power, while using multiple an- 14 tennas increases the complexity in terms of processing power and number of transceivers. If spectral efficiency stays about moderate level, increasing the bandwidth at a higher frequency is the only solution to higher data rates. This makes 60-GHz band the perfect solution for its up to 7 GHz bandwidth. 2.2.3 Atmospheric Absorption In terms of large-scale fading and path loss, although there are losses due to phenomena such as water vapor, foliage, and vegetation [9], for 60-GHz channel, the major loss is due to oxygen absorption and rain attenuation. Oxygen absorption As mentioned before, oxygen absorption is a unique characteristic of the 60GHz band. This gives the 60-GHz band some benefits, along with disadvantages such as path loss. The peak of oxygen absorption is very localized in frequency. The related attenuation can take high value near 60 GHz, and decreases by an order of magnitude when the frequency is near 66 GHz. For simplification, according to [35, 36], the path loss per distance can be modeled by an expression with respect to frequency f for the band 60 GHz–66 GHz with maximum error of 5%: aoxy [dB/km](f[GHz] ) = 15.10 − 0.104 × (f − 60)3.26 11.35 + (f − 60)2.25 − 5.33 × (f − 63)1.27 60 ≤ f ≤ 63 (2.15) 63 ≤ f ≤ 66 In our calculation, we adopt 15.1 dB/km as the loss coefficient [9]. Rain attenuation Clouds of ice crystals and snow do not cause appreciable attenuation [9]. However, when an electromagnetic wave propagates though a rain cell, it encounters a great amount of water droplets with different radii. The cumulative effect of electromagnetic wave penetration causes a considerable attenuation. Similar 15 to oxygen absorption, the per kilometer attenuation could also be modeled with respect to the rain fall rate R and frequency f [36] arain [dB/km](f[GHz] ) = a × Rb a = 4.09 · 10−2 f 0.669 b = 2.63f −0.272 for 54 ≤ f ≤ 180 (2.16) for 25 ≤ f ≤ 164 In our model, we adopt 12.5 mm/h and 25 mm/h as the fall rate for moderate and heavy rain respectively [10]. The corresponding values of arain are 5.6 dB/km and 10.2 dB/km. 2.2.4 60-GHz MMW Fading Model Small-scale fading in multipath channels occurs due to the incoherent superposition of a great number of multipath components, each having a different phase variation over time or frequency, or both. Among with all fading models used for 60 GHz in literature, the most popular models for possibly time-variant multipath channels are Rayleigh or Rice fading model [9]. The difference between Rayleigh distribution and Rice distribution is that Rayleigh models fading channels without line-of-sight (LOS) path, while Rice fading describes the case when there is a dominant incoming wave which is the LOS component. One can view the Rayleigh distribution as one extreme case of the Rice distribution when the dominant path is insignificant. For simplicity, we assume the a LOS always exists since the FSO link requires a LOS. Hence, we choose the Rice distribution as our fading model and use it throughout this thesis. If we denote 2σr2 the mean power of the sum of the non-specular (non-LOS) component and s2 denotes the LOS component, then the pdf can be expressed as: pR (r) = r r 2 + s2 r·s exp(− )I0 ( 2 ), 2 2 σr 2σr σr r ∈ [0, ∞) (2.17) where In (·) denotes the modified Bessel function of the first kind and n-th order. 16 0.5 K=1 K=2 K=3 0.45 0.4 0.35 0.3 pR(r) 0.25 0.2 0.15 0.1 0.05 0 0 1 2 3 4 5 6 7 8 r Figure 2.3: Sample Rice fading intensity pdfs Let’s define the Rice factor K = s2 /2σr2 , the ratio of the power of the dominant multipath component to the average fading power of the remaining, non-specular multipath [9]. The normalized pdf of (2.17) can be re-stated as pR (r) = 2(K + 1)r · exp(−K − (K + 1)r 2 )I0 (2r K(K + 1)) (2.18) Notices that, when K −→ 0, the pdf becomes Rayleigh distribution. Figure 2.3 depicts the pdf of Rice distributions with different K values. 2.3 State-of-the-art of Hybrid RF/FSO Systems In order to guarantee the connectivity of a data network, commonly link availability is required to meet some requirement for different deployment pro- 17 poses. For enterprise applications, link availability is required to be greater than 99%, called enterprise-class availability. For carrier or internet service providers (ISP), the availability needs to be greater than 99.999%. This is called carrier-class availability. In comparison with other last mile technologies such as fiber optic, the main barrier for FSO deployment is that its availability is variable and difficult to predict due to weather. For example, fiber optic cable attenuates the transmission power at a constant, predictable rate. Current multi-mode fiber attenuates at 2 to 3 dB/km, and single-mode fiber attenuates at 0.2 to 0.5 dB/km. On the other hand, atmospheric attenuation can vary from 0.2 dB/km in clear weather to 250 dB/km in very heavy fog [10]. The large attenuation values in certain weather conditions can reduce the up-time or link availability of FSO systems. The link availability can be used for link budget calculation. From Section 2.2 it is known that the large-scale attenuation depends on both the distance of the link and the weather conditions. Hence, for the same transmitted power, the trade off is between the link distance and the link availability. The maximum link distance for specific types of link availability requirements can be determined for system transmitted power and weather conditions. Reference [10] gives a detailed examination of the distance limitation of FSO systems for both carrier and enterprise applications. It is concluded from [10] that under some weather conditions such as heavy fog, the maximum link distance is only at the range of a few hundred meters, regardless of the value of the transmitted power or the size of the receive optics. To overcome this major drawback, the most cost-effective solution is to incorporate an RF back-up. The only weather condition which negatively affects both links is heavy rain and heavy fog simultaneously. However, this does not happen since rain droplets would absorb the suspended fog water droplets [14]. A great amount of research has been done on hybrid FSO/RF links [11, 12, 14, 15, 37, 38]. One is to use RF link as a complementary channel [11]. In this 18 scheme, data is sent exclusively through the FSO link until a drop of the link quality is observed. In this case, data will switch to the RF link completely. This requires the feedback of the link availability to the transmitter. Another scheme uses RF link as a redundant link [13]. Identical information is sent into both links. At the receiver, the data at the link with higher signal quality will be picked up. This does not require link condition knowledge before hand, with the trade-off of redundant data sent in both links. All these schemes proposed or used do not give consideration on utilizing the diversity of the two channels to maximize the capacity. Recently, a scheme proposed by Vangala et al. [37] utilizes hybrid channel codes to achieve a high efficiency and exploits diversity of FSO/RF channels. A hybrid channel code is the designed that a sequence first goes through a non-uniform encoder followed by a rate-adaptive code convertor. Non-uniform codes are used for data transmission over a set of parallel independent sub-channels [39]. In [37], a codeword of length n is generated by a single low-density parity check (LDPC) encoder. The codeword is then transmitted over the set of channels such that n(j) bits in any codeword are transmitted over the jth channel. Additionally, in order to maintain an acceptable quality of service when channels are time-variant, Vangala et al. proposes to using punctured coded so that code rate is adaptive to channel conditions. However, Vangala et al. only provides the concept of his design, no actual code searching method and results are presented. Also, the proposed system is not robust to any weather conditions since feedback of the weather condition is required. In this thesis, we could propose a joint BICM scheme which splits each codeword accordingly into two bit streams and transmits them via both FSO/RF channels. A diversity is reached. 19 Table 2.1: Weather–dependent parameters of FSO and RF sub–systems [9, 40, 25, 26]. Weather conditions a1 [dB/km] arain [dB/km] Cn2 Clear air 0.43 0 5 · 10−14 Haze 4.2 0 1.7 · 10−14 Light fog 20 0 3 · 10−15 Moderate fog 42.2 0 2 · 10−15 Heavy fog 125 0 1 · 10−15 Moderate rain (12.5 mm/h) 5.8 5.6 5 · 10−15 Heavy rain (25 mm/h) 9.2 10.2 4 · 10−15 20 Chapter 3 Hybrid FSO/RF System Model In Chapter 2, we have described the channel fading models for FSO channel and RF channel, respectively. In this chapter, we derive the signal-to-noise ratio (SNR) for both channels separately based on the channel fading models detailed in Chapter 2, as well as the link budget calculation. The link budget calculation is dependent on weather condition and link distance. Here, we adopt the weather attenuation models and values of attenuation coefficient provided in Chapter 2. We later describe the transmitter and receiver structures of your system. 3.1 System Model In this thesis, we assume that both the FSO and the RF sub–systems have one transmit and one receive aperture, respectively. Furthermore, we assume that both the FSO and RF channels are constant during one codeword. Because of the very high data rates considered, this is a reasonable assumption even for mobile hybrid RF/FSO systems [12, 38]. FSO Channel: We consider an FSO system with intensity modulation and direct detection (IM/DD) [6]. The photocurrent generated by the photodetector is integrated over one symbol interval and the constant bias due to background radiation is removed. Under these conditions, for high–energy, 21 shot noise (caused by background radiation) limited FSO systems the sufficient statistic in the k1 th symbol interval can be modeled as [27]–[40] r1 [k1 ] = P1 g1 h1 ρx1 [k1 ] + n1 [k1 ], where P1 , g1 , h1 , ρ, x1 [k1 ] ∈ X1 {0, 2} (M1 (3.1) |X1 | = 2), and n1 [k1 ] are the transmitted optical power, the average gain of the FSO link, the FSO fading gain with E{h1 } = 1, the responsivity of the photodetector, the transmitted OOK symbol, and real–valued additive white Gaussian noise (AWGN) with variance σ12 E{|n1 [k1 ]|2 }, respectively. The Beers–Lambert Law dictates g1 . For the probability density function (pdf) f1 (h1 ) of h1 we adopt the popular Gamma–Gamma fading model 2.6. Based on (3.1) we define the electrical SNR of the FSO link as γ1 P12 g12ρ2 (E{P1g1 h1 ρx1 [k1 ]})2 = . σ12 σ12 (3.2) RF Channel: The received signal of the RF channel in the k2 th symbol interval is modeled as r2 [k2 ] = √ P2 g2 h2 x2 [k2 ] + n2 [k2 ], (3.3) where P2 , g2 , h2 , x2 [k2 ] ∈ X2 , and n2 [k2 ] denote the RF transmit power, the average power gain of the RF link, the RF fading gain with E{h22} = 1, the modulated RF symbol with E{|x2 [k2 ]|2 } = 1 taken from M2 –ary symbol alphabet X2 (e.g. quadrature amplitude modulation (QAM)), and complex AWGN with variance σ22 E{|n2 [k2 ]|2 }, respectively. Assuming a carrier frequency of 60 GHz, the effective gain can be modeled as [9] g2 [dB] Gt + Gr − 20 log10 4πL λ2 − aoxy L − arain L, (3.4) where Gt and Gr denote the transmit and receive antenna gains, respectively, λ2 is the wavelength of the RF system, and aoxy and arain are the attenuations caused by oxygen absorbtion and rain (both in dB/km), respectively. The RF noise variance is given by σ22 [dB] = BN0 + NF , where B, N0 , and NF denote 22 the RF bandwidth, the noise power spectral density (in dBm/MHz), and the noise figure of the receiver, respectively [9]. The fading gain h2 can be modeled as Ricean distributed [4, 8, 9], where the Ricean factor K depends on various factors such as link distance, antenna height, and the environment, and may also change with time in mobile applications [12, 38, 41]. The SNR of the RF link is defined as γ2 √ E{| P2 g2 h2 x2 [k1 ]|2 } P2 g 2 = 2 . 2 σ2 σ2 (3.5) Since both a1 and arain are weather dependent cf. Chapter 2, we show the SNRs γ1 and γ2 in Fig. 3.1 for different weather conditions, where we have adopted P1 = P2 = 10 mW and the values given in Tables 3.1 and 2.1 for all relevant parameters of the FSO and RF sub–systems. The values in Tables 3.1 and 2.1 were taken from [8, 9, 23, 40, 25, 26]. The values for a1 and arain are summary of those stated in Section 2.1.2 and Section 2.2.3 3.2 Transmitter Structure The transmitter structure of the proposed hybrid RF/FSO system is shown in the top part of Fig. 3.2. A sequence of Kc information bits b [b[1] b[2] . . . b[Kc ] is encoded with a convolutional encoder (CC) of rate Rc to produce the coded sequence c [c[1] c[2] . . . c[Nc ]] of length Nc Kc /Rc . Subsequently, c is divided into two sub–sequences, c1 and c2 , according to assignment pattern (AP) a {0, 1}N , which contains n1 “1”s and n2 = N − n1 “0”s. For this purpose c is grouped into blocks of N bits and within each block bits in positions where a has a “1” are assigned to c1 , while bits in positions where a has a “0” are assigned to c2 . Thus, c1 and c2 are sequences of lengths n1 Nc /N and n2 Nc /N, respectively. The bits in c1 and c2 are interleaved with interleavers Π1 and Π2 , respectively, and the interleaved bit sequences are mapped to symbols x1 [k1 ] and x2 [k2 ] using Gray labeling, respectively.1 The symbol rates of 1 We note that since binary modulation is used over the FSO link, interleaver Π1 is not needed for the considered channel model. However, we have included Π1 since, in practice, 23 Table 3.1: Parameters of FSO and RF sub–systems [8, 9, 23, 40]. FSO System Parameter Symbol Value Wavelength λ1 1550 nm Symbol rate 1/T1 1 Gbit/s Responsivity ρ 0.5 V1 Noise variance σ12 10−14 A2 RF System Parameter Symbol Value Carrier frequency fc 60 GHz Bandwidth B 250 MHz Symbol rate 1/T2 250 Msym/s Transmit antenna gain Gt 44 dBi Receive antenna gain Gr 44 dBi Attenuation (Oxygen) aoxy 15.1 dB/km Ricean factor K 6 dB Noise power spectral density N0 -114 dBm/MHz Receiver noise figure NF 5 dB 24 60 Clear air Haze Moderate rain Heavy rain Light fog Moderate fog Heavy fog RF 50 40 SNR [dB] 30 20 FSO 10 0 −10 0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 L [m] Figure 3.1: SNR of FSO (γ1 ) and RF (γ2 ) systems vs. link distance L for transmit powers of P1 = P2 = 10 mW and different weather conditions. The values for a1 and arain have been taken from Table 2.1. All other link parameters have been chosen as in Table 3.1. both sub–systems are given by 1 Tν nν log2 (Mν )Rc NTb symbols , s ν ∈ {1, 2}, (3.6) where T1 and T2 denote the symbol durations of the FSO and the RF sub– systems, respectively, and 1/Tb is the information bit rate. Therefore, assuming FSO and RF sub–channels with coded bit rates of 1/T1 and log2 (M2 )/T2 , respectively, n1 and n2 have to fullfill n1 T2 = . n2 log2 (M2 )T1 (3.7) burst errors may be caused by e.g. bird fly–through in the FSO link. In this case, the interleaver would be beneficial. 25 c1 b c CC MOD Π2 MOD Π−1 1 DEMOD Π−1 2 DEMOD x2[k2] λ1 (κ1, cˆ[κ1 ]) P/S (a) r2[k2 ] x1[k1] S/P (a) c2 r1[k1 ] Π1 λ(κ, cˆ[κ]) ˆ b VD λ2 (κ2, cˆ[κ2 ]) Figure 3.2: Transmitter and receiver structure of proposed hybrid RF/FSO system. S/P: Serial–to–parallel conversion. P/S: Parallel–to–serial conversion. VD: Viterbi decoder. We note that both the CC and the AP a have to be optimized for maximum performance, cf. Sections 4 and 5. 3.3 Receiver Structure To simplify our notation, we define the normalized received signals rν [kν ] rν [kν ]/σν . As shown in the bottom part of Fig. 3.2, the received symbol streams of both links are demodulated to obtain the bit metrics [19] λν (κν , cˆ[κν ]) min i {ξν |rν [kν ] − xν ∈Xν (ˆ c[κν ]) √ γν hν xν |2 }, ν ∈ {1, 2}, (3.8) where Xνi (c) denotes the subset of all symbols xν ∈ Xν whose label has the value c ∈ {0, 1} in position i, 1 ≤ i ≤ log2 Mν , ξ1 1/2 and ξ2 1 account for the fact that n1 [k1 ] and n2 [k2 ] are real and complex AWGN, respectively, and κν is determined by kν and i. The bit metrics are deinterleaved with Π−1 1 and Π−1 2 , respectively, and parallel–to–serial converted based on AP a to produce a single sequence of bit metrics λ(κ, cˆ[κ]), 1 ≤ κ ≤ Nc , corresponding to the coded sequence c. Finally, a standard Viterbi decoder (VD) is used to generate ˆ for the transmitted bit sequence based on these bit metrics. an estimate b 26 Chapter 4 Performance Analysis and System Optimization 4.1 Asymptotic Pairwise Error Probability For calculation of the union bound on the CEP, the PEP is needed. Consider an error event of length d, where d1 erroneous bits correspond to symbols x1 [k1,i ] ∈ X1 , 1 ≤ i ≤ d1 , transmitted over the FSO link and d2 d − d1 bits correspond to symbols x2 [k2,i ] ∈ X2 , 1 ≤ i ≤ d2 , transmitted over the RF link. The corresponding metric difference can be expressed as [19] 1 ∆ = 2 d1 i=1 d2 + i=1 |r1 [k1,i ] − |r2 [k2,i ] − √ γ1 h1 x¯1 [k1,i ]|2 − |r1 [k1,i ] − √ γ2 h2 x¯2 [k2,i ]|2 − |r2 [k2,i ] − √ γ1 h1 x1 [k1,i ]|2 √ γ2 h2 x2 [k2,i ]|2 ,(4.1) where xν [kν,i ] = x¯ν [kν,i ] ∈ Xν , ν ∈ {1, 2}. Based on (4.1) the conditional PEP can be obtained as ¯ |h1 , h2 ) Pe (x → x Pr{∆ < 0} = Q ¯ )/2 , γ1 h21 d1 + γ2 h22 d22 (x, x (4.2) ¯ contain the involved xν [kν,i ] and x¯ν [kν,i ], respectively, where vectors x and x and d2 ¯) d22 (x, x i=1 |x2 [k2,i ] − x¯2 [k2,i ]|2 27 (4.3) While the exact average PEP can be expressed as ¯) Pe (x → x ¯ |h1 , h2 )} Eh1 ,h2 {Pe (x → x (4.4) , and cannot be calculated in closed form, we derive the average PEP for asymptotic case (γ1 1 and γ2 1). First by using the modified version of the Gaussian Q–function 1 Q(x) = π π/2 e−x 2 /(2 sin2 θ) dθ (4.5) 0 and (4.2), the average PEP can be expressed as π/2 1 ¯) = Pe (x → x π I1 (θ) I2 (θ) dθ, (4.6) 0 where Eh1 {exp(−γ1 h21 d21 /(2 sin2 θ))} (4.7) ¯ )/(4 sin2 θ))} Eh2 {exp(−γ2 h22 d22 (x, x (4.8) I1 (θ) and I2 (θ) To arrive at a closed–form solution, we concentrate on the high SNR regime (γ1 1 and γ2 1). In this case, it has been shown in [42, 24] that the pdf fν (hν ) may be replaced by the first term of its Taylor series expansion fν (hν ) = mν hlνν + o(hlνν ), where m1 hν ≥ 0, (αβ)β Γ(α − β)/[Γ(α)Γ(β)], l1 ν ∈ {1, 2}, β − 1, α > β, α − β = ZZ for 2(1 + K)e−K , l2 Gamma–Gamma fading [24] and m2 (4.9) 1 for Ricean fading [42]. We note that the conditions α > β and α − β = ZZ are valid in practice, cf. (2.7), (2.8), [7]. Exploiting (4.9) and the integral ∞ xm−1 exp(−sx2 ) dx = s−m/2 Γ(m/2) (4.10) 0 for m > 0, s > 0 [43, Eq. (3.4782)], For I1 (θ) we obtain ∞ I1 (θ) = 0 exp − γ1 d1 h2 2 sin2 θ 28 m1 hβ−1 dh (4.11) h , sin θ Let u = 4.11 becomes ∞ I1 (θ) = m1 exp − 0 = 2β/2 m1 Γ(β/2) β/2 β/2 γ 1 d1 γ 1 d1 u 2 2 uβ−1 sinβ−1 θdu sinβ θ (4.12) For I2 (θ) we obtain ∞ I2 (θ) = 0 h , sin θ Let u = exp − γ2 d2 h2 4 sin2 θ m2 hdh (4.13) 4.13 becomes ∞ I2 (θ) = m2 sin θ exp{ 0 ¯ )u2 γ2 d22 (x, x }du 4 4m2 = sin2 θ. ¯) γ2 d22 (x, x (4.14) Combining (4.6) with (4.12) and (4.14) and exploiting π/2 sinm−1 x dx = Γ(1/2)Γ(m/2)/[2Γ((m + 1)/2)] (4.15) 0 [43, Eq. (3.6211)] leads to . ¯) = Pe (x → x C β/2 γ1 γ2 · 1 β/2 ¯) d1 d22 (x, x ˙ ≤ C β/2 d2min γ1 γ2 · 1 β/2 d1 d2 , (4.16) where C 2β/2+2 (1 + K)e−K (αβ)β Γ(α − β)Γ(β/2)Γ(β/2 + 3/2) √ , [ πΓ(α)Γ(β)Γ(β/2 + 2)] (4.17) dmin denotes the minimum Euclidean distance of X2 , and d1 > 0 and d2 > 0 were assumed for derivation of (4.16). We also exploited the inequality ¯ ) ≥ d2min d2 in (4.16). Eq. (4.16) shows that codes with d1 > 0 and d2 > 0 d22 (x, x guarantee that the asymptotic PEP decreases with increasing SNRs γ1 and γ2 . In this case, the RF and the FSO sub–systems both contribute to the reliability of the transmission and diversity is achieved. Thus, in this thesis, we will only consider codes which achieve d1 > 0 and d2 > 0. We note that Fig. 3.1 shows 29 that in practice it might not be possible to simultaneously achieve γ1 and γ2 1 1, which was assumed for derivation of (4.16), under all weather conditions and for all link distances. Nevertheless, codes guaranteeing d1 > 0 and d2 > 0 are robust against low SNR in one sub–channel and are able to exploit the high SNR in the other sub–channel. 4.2 Asymptotic Codeword Error Probability From (4.16) we can obtain an upper bound on the asymptotic CEP as [22] ˙ Pw ≤ F CF β/2 d2minγ1 γ2 w(d1, d2 ) β/2 d1 d2 d1 d2 (4.18) (4.19) where w(d1, d2 ) denotes the number of error events with Hamming distances d1 and d2 for the FSO and RF sub–channels, respectively. For the considered problem w(d1 , d2 ) depends on both the CC and the AP a. w(d1 , d2 ) is given by the coefficients of the generalized transfer function (GTF) [21, 22] w(d1 , d2)D1d1 D2d2 . T (D1 , D2 ) = d1 (4.20) d2 We will explain in Section 5 how T (D1 , D2 ) can be efficiently computed for the problem at hand. It is well known that the union bound in (4.18) is not necessarily tight. However, despite this deficiency, it was found in [22] that for homogeneous block fading channels codes optimized according to the union bound still outperform codes optimized according to other criteria [17, 18]. Therefore, in Chapter 5, we will introduce a CC and AP search procedure which is based on (4.18). 4.3 Transmit Power Assignment The total transmit power is given by Pt = P1 + P2 and should be divided in an optimal way between P1 and P2 . Furthermore, in practice, the maximum 30 FSO transmit power P1,max and the maximum RF transmit power P2,max are both limited [3, 44]. From (3.2), (3.5), and (4.18) we observe that the upper bound on Pw is proportional to 1/(P1β P2 ). Since our goal is the minimization of the upper bound on Pw , we obtain the following optimization problem: P1−β P2−1 min P1 ,P2 (4.21) s.t. 0 ≤ P1 ≤ P1,max (4.22) 0 ≤ P2 ≤ P2,max (4.23) P1 + P2 ≤ Pt (4.24) The optimum values for P1 and P2 can be found with the standard Karush– Kuhn–Tucker approach [45] for constraint optimization problems. Assuming there is a feasible solution, i.e., Pt ≤ P1,max + P2,max , the optimum power assignment is given by P1 = β P, 1+β t P2 = 1 P, 1+β t if Pt ≤ min P1 = P1,max , P2 = Pt − P1,max , if 1+β P1,max β 1+β P1,max , β (1 + β)P2,max < Pt ≤ (1 + β)P2,max P1 = Pt − P2,max , P2 = P2,max , if (1 + β)P2,max < Pt ≤ (4.25) 1+β P1,max β For the case Pt ≤ min{(1 + 1/β)P1,max , (1 + β)P2,max , P1,max + P2,max } it is interesting to compare the optimum power assignment according to (4.25) with the pragmatic assignment P1 = P2 = Pt /2. From (4.18) and (4.25) we observe that the asymptotic power gain achievable with the optimum power assignment is β β 1+β GP (β) = 2 . (4.26) 1+β Not surprisingly, for β = 1 (4.26) yields GP (1) = 1 since in this case equal power assignment is optimum, cf. (4.25). For β = 2, β = 5, and β = 10 we obtain from (4.26) power gains of 0.24 dB, 1.15 dB, and 1.69 dB, respectively. Furthermore, it can be shown that GP (β) ≤ 2, ∀β > 0, i.e., the maximum power gain is limited to 3 dB, where the maximum gain is achieved for β → 0 and β → ∞. It is interesting to note that provided d1 > 0 and d2 > 0 the optimum power assignment is independent of the CC and the AP but only depends on the statistics of the FSO channel. 31 4.4 Bound on Rc For the CC and AP search it is of interest to find a bound on Rc that guarantees that codes with d1 > 0 and d2 > 0 can be found. In this context, we note that the overall code rate including the effect of a is Rν NRc /nν for sub–channel ν ∈ {1, 2}. It is clear that dν > 0 can only be achieved if Rν ≤ 1 [29]. Thus, we obtain the necessary condition Rc ≤ Rc,max min ν∈{1, 2} nν N = min ξ∈{n1 /n2 , n2 /n1 } 1 1+ξ . (4.27) A simple code with Rc = Rc,max that achieves dν > 0, ν ∈ {1, 2}, is the repetition code with an arbitrary AP having n1 ”1”s and n2 ”0”s . Thus, (4.27) is sufficient and necessary for the existence of a full diversity code (i.e., dν > 0, ν ∈ {1, 2}). Note that by exploiting (3.7) Rc,max may equivalently be formulated in terms of T1 , T2 , and log2 (M2 ). If both the FSO and the RF sub–channels support the same coded bit rate, n1 /n2 = 1 follows from (3.7), and we obtain from (4.27) Rc ≤ 1/2. If the FSO and RF sub–channels support different coded bit rates, i.e., n1 /n2 = 1, the maximum code rate Rc,max is lower than 1/2. 32 Chapter 5 Code and Assignment Pattern Search In this chapter, we explain the proposed CC and AP search procedure. Furthermore, we outline how the GTF is computed and present some CC and AP search results. 5.1 Search Algorithm Since we adopt the upper bound on the CEP in (4.18) as optimality criterion, the optimum CC and AP have to minimize F in (4.19). To limit the complexity of the search procedure and following the reasoning in [22] we limit d1 and d2 in the sums in (4.19) to d1 + d2 ≤ 2df − 1, where df is the free distance of the code. We assume that parameters T1 , T2 , and M2 of the hybrid RF/FSO system are given. For the proposed CC and AP search we first have to select the desired constraint length q and rate Rc = kc /nc ≤ Rc,max (kc and nc are integers) of the CC and the length N of the AP. Subsequently, we generate all possible CCs with the prescribed q and Rc and all possible APs meeting (4.27). For all possible combinations of the generated CCs and APs we calculate the corresponding F , and retain that combination of code and AP as the optimum one which has the minimum F . 33 Now, some remarks on how to choose q and N are in order. The decoding complexity is only affected by q but not by N. Larger q will generally yield a better performance at the expense of a higher decoding complexity. The complexity of the proposed exhaustive search procedure increases with both q and N. For larger q there are more generator polynomials to be examined and for larger N there are more possible APs for each combination of n1 and n2 . On the other hand, since larger N provide more flexibility in choosing the AP, larger N generally yield a superior performance as compared to smaller N. 5.2 Computation Method In the proposed search procedure, F has to be calculated for each considered combination of CC and AP. Therefore, the GTF T (D1 , D2 ) has to be computed for each combination of CC and AP. For simplicity, we assume ideal interleavers and neglect any edge effects caused by finite length codewords. As mentioned before, the coefficients w(d1, d2 ) of the GTF denote the number of error events with Hamming distances d1 and d2 in the FSO and RF sub–channels, respectively. Without loss of generality, we may assume that an all–zero sequence was transmitted and the decoded codeword constitutes the corresponding error event. Thus, w(d1, d2 ) is the number of codewords c1 and c2 with weights d1 and d2 , respectively. Based on this observation, we can efficiently compute T (D1 , D2 ) with the approach proposed in [22]. For this purpose, we consider the trellis diagram of a CC with rate Rc = kc /nc and 2q states Sj , 1 ≤ j ≤ 2q , with S1 being the all–zero state. For the tth trellis segment we define the elements Aij of 2q ×2q matrix At (D1 , D2 ) as follows [22, 46]. If a trellis branch connects state Si with state Sj , Aij = D1w1 D2w2 where w1 and w2 denote the number of “1”s assigned to c1 and c2 in the considered trellis segment according to the AP, respectively. If states Si and Sj are not connected by a trellis branch, Aij = 0. Considering a codeword of length T T0 · lcm(nc , N), where T0 is an integer and lcm(nc , N) is the lowest 34 common multiple of nc and N, the GTF can be obtained as [46, 22] T (D1 , D2 ) = A11 (D1 , D2 ) − 1, (5.1) where A11 (D1 , D2 ) is the first element of the first column of matrix lcm(nc ,N ) T A(D1 , D2 ) t=1 At (D1 , D2 ) = t=1 T0 At (D1 , D2 ) . (5.2) In (5.2), we have exploited the fact that the CC trellis including the effect of the AP is periodic with period lcm(nc , N). Based on (5.1) and (5.2) the weights w(d1, d2 ), d1 +d2 ≤ 2df −1, can be efficiently computed using standard software such as MATLABTM . 5.3 CC and AP Search Results A few example CCs and APs found with the proposed search method are provided in Table 5.1 for RF channels with 16–QAM (Rc = 1/2 and 1/4) and 4–QAM (Rc = 1/3), respectively. These codes were optimized for K = 6 dB and α = 2.5 and β = 2, which are valid for haze and L = 3.5 km. We note that the optimum CCs and APs are not very sensitive to the exact values of K, α, and β. Thus, we employed the codes provided in Table 5.1 for all simulation results presented in Section 7. We also provide the free distance df and the resulting value for F , cf. (4.19), in Table 5.1. For each code in Table 5.1 two different AP lengths N are considered. As expected, the larger N generally leads to a smaller F , and thus better performance, due to the larger flexibility in assigning the coded bits to the sub–channels. However, for the considered examples the potential gains by using the larger N are comparatively small. 35 Table 5.1: Convolutional codes and assignment patterns optimized for α = 2.5, β = 2, and K = 6 dB. The generator polynomials g and the assignment patterns are both in octal representation. Rc = 1/2 and 16–QAM q g N a df F 3 (5, 7) 4 (21) 5 58.0 3 (5, 7) 6 (46) 5 55.1 4 (15, 17) 4 (30) 6 105.0 4 (15, 17) 6 (51) 6 93.2 5 (43, 35) 4 (22) 7 170.0 5 (43, 35) 6 (62) 7 162.8 6 (53, 75) 4 (50) 8 368.3 6 (53, 75) 6 (52) 8 346.5 7 (133, 171) 4 (41) 10 3274.8 7 (133, 171) 6 (25) 10 3047.2 Rc = 1/4 and 16–QAM q g N a df F 3 (5, 7, 7, 7) 8 (306) 10 31.0 3 (5, 7, 7, 7) 12 (3306) 10 29.8 8 (230) 13 41.9 4 (13, 15, 15, 17) 12 (6314) 13 39.8 4 (13, 15, 15, 17) Rc = 1/3 and 4–QAM q g N a df F 3 (5, 7, 7) 3 (3) 8 66.6 3 (5, 7, 7) 6 (53) 8 59.8 4 (13, 15, 17) 6 (55) 10 92.3 4 (13, 15, 17) 9 (665) 10 87.8 5 (25, 33, 37) 15 (33535) 12 214.2 6 (47, 53, 75) 15 (76266) 13 242.0 36 Chapter 6 Cut-off Rate Analysis Since the cut–off rate constitutes the ultimate performance limit for convolutionally coded systems [29], in this section, we evaluate the cut–off rate of the proposed hybrid RF/FSO system. To simplify the derivation, we assume random APs assigning on average n1 and n2 bits to the FSO and RF channels, respectively. For simplicity of notation, we drop time indices k1 and k2 in this section. Taking also into account the bit interleaver, we have 1 + log2 (M2 ) statistically independent parallel channels [19] for which the cut– off rate R0 (h1 , h2 ) for given channel gains h1 and h2 can be obtained as [47] R0 (h1 , h2 ) = (R1,max + R2,max )[1 − log2 (1 + B(h1 , h2 ))] [bit/s], where R1,max 1/T1 , R2,max (6.1) log2 (M2 )/T2 , and B(h1 , h2 ) n1 n2 B1 (h1 ) + B2 (h2 ) N N (6.2) with B1 (h1 ) Er1 ,c p1 (r1 |x1 = c¯, h1 )/p1 (r1 |x1 = c, h1 ) = exp −h21 r1 /2 , (6.3) log2 (M2 ) B2 (h2 ) i=1 37 B2,i (h2 ) , log2 (M2 ) (6.4) , and B2,i (h2 ) Er2 ,c p2 (r2 |x2 , h2 ) [19]. p2 (r2 |x2 , h2 ) / x2 ∈X i (c) x2 ∈X i (¯ c) 2 2 (6.5) Here, c ∈ {0, 1}, c¯ denotes the complement of c, and p1 (r1 |x1 , h1 ) and p2 (r2 |x2 , h2 ) are the conditional pdfs of received signals r1 and r2 , respectively. We note that the Bhattacharyya factors B1 (h1 ) and B2,i (h2 ), 1 ≤ i ≤ log2 (M2 ), can be efficiently computed via Monte–Carlo simulation. For comparison, the cut–off rates of the FSO and RF links are given by R0,1 (h1 ) = R1,max [1 − log2 (1 + B1 (h1 ))] [bit/s]. (6.6) R0,2 (h2 ) = R2,max [1 − log2 (1 + B2 (h2 ))] [bit/s]. (6.7) and respectively. Thus, a system which knows the sub–channel SNRs at the transmitter and performs independent BICM over both channels can achieve a cut– off rate of R0,SNR (h1 , h2 ) = R0,1 (h1 ) + R0,2 (h2 ). (6.8) Note that the sub–channel SNRs are weather dependent and do not have to be known if the proposed hybrid RF/FSO structure is adopted. This advantage over a system which knows the sub–channel SNRs comes at the expense of some potential performance degradation. In particular, assuming identical power assignment in both cases the inequality R0 (h1 , h2 ) ≤ R0,SNR (h1 , h2 ) can be established by exploiting the convexity of the logarithm, log2 (1 + B(h1, h2 )) ≥ n1 log2 (1 + B1 (h1 ))/N + n2 log2 (1 + B2 (h2 ))/N, and (3.6) and (6.1)–(6.8). On the other hand, the rate of a system which does not know the SNR at the transmitter but performs independent BICM over both sub–channels is limited by R0,ind (h1 , h2 ) = 2 min{R0,1 (h1 ), R0,2 (h2 )}, (6.9) since such a system has to assign the same rates to both sub–channels, and thus, will be limited by the weaker sub–channel. 38 Since we consider a block fading channel, it is of interest to define the outage probability Pout = Pr{R0 (h1 , h2 ) < Rtarget }, (6.10) i.e., the probability that the considered system cannot support a rate of Rtarget . The outage probabilities for systems with independent BICM over both sub– channels with known and unknown sub–channel SNR at the transmitter are defined analogously to (6.10). 39 Chapter 7 Simulation Results In this section, we study the performance of the proposed hybrid RF/FSO system with joint BICM over both links under various weather conditions. As performance measures we consider the outage probability as defined in Section 6 and the frame error rate (FER). The adopted values for the various system and channel parameters are given in Tables 3.1 and 2.1. 7.1 Outage Probability In this subsection, we adopt 16–QAM for the RF link, i.e., the coded bit rate of both sub–channels is R1,max = R2,max = 1 Gbit/s, cf. Table 3.1. Figs. 7.1 and 7.2 show the outage probability as a function of the transmit power Pt for the proposed hybrid RF/FSO system with joint BICM over both sub–channels for target rates of Rtarget = 1 Gbit/s and Rtarget = 500 Mbit/s, respectively. Performance under clear air, haze, and moderate rain are considered for a link distance of L = 2.5 km. For comparison, the outage probabilities for hybrid RF/FSO systems performing independent BICM over both sub–channels are also shown in the presence and absence of CSI in form of the sub–channel SNRs at the transmitter, respectively. The corresponding cut–off rates are given in (6.8) and (6.9), respectively. For all considered schemes equal power assignment to the RF and FSO sub–channels is assumed. Since with de- 40 creasing visibility the FSO link attenuation increases while the severity of the fading decreases, cf. Table 2.1, for all considered schemes clear air yields the best performance for small Pt but haze results in a better performance for larger Pt due to the steeper slope of the outage curves. For all considered weather conditions independent BICM without CSI performs relatively poorly as performance is limited by the weaker sub–channel. Independent BICM with known CSI yields the best performance since the rates of both sub–channels can be perfectly adapted to the channel conditions. However, the proposed joint BICM scheme, which does not require CSI at the transmitter, achieves an excellent performance, especially for Rtarget = 500 Mbit/s, where the loss compared to independent BICM with CSI is negligible in the considered range of outage probabilities. The improved robustness of joint BICM for lower target rates is due to the stronger channel codes that can be used in this case, which provide better protection against low SNRs in one sub–channel. 7.2 Frame Error Rate (FER) Throughout this section we assume that a frame contains 200 information bits. Besides the proposed joint BICM scheme we also show results for repetition coding over both channels with selection diversity at the receiver (RC–SD) as proposed in [13]. In RC–SD the same data is sent over both links and the receiver selects that frame which is deemed more reliable. For simplicity, we assume that this selection is perfect. For the proposed hybrid RF/FSO scheme we adopt the CCs and APs provided in Table 5.1 for all considered weather conditions. Optimizing the CCs and APs for the actual weather conditions would result in small additional performance gains at the expense of a considerable increase in signaling overhead (estimation of α, β, and K) and complexity (on–line optimization of CCs and APs). If the maximum RF and FSO transmit powers are not constrained (Figs. 7.3, 7.5, 7.6), we assign equal powers to both sub–channels as for the considered channel models the additional gains 41 0 Pr{R0(h1, h2) < 1 Gbit/s} 10 −1 10 −2 10 Clear Air Haze Rain Joint BICM (proposed) Independent BICM (with CSI) Independent BICM (without CSI) −3 10 0 5 10 15 20 Pt [dBm] 25 30 35 40 Figure 7.1: Outage probability for target rate of Rtarget = 1 Gbit/s for hybrid RF/FSO system with joint BICM (proposed), independent BICM with CSI at the transmitter, and independent BICM without CSI at the transmitter, respectively. 16–QAM in RF link and L = 2.5 km. 42 0 Pr{R0(h1, h2) < 500 Mbit/s} 10 −1 10 −2 10 Clear Air Haze Rain Joint BICM (proposed) Independent BICM (with CSI) Independent BICM (without CSI) −3 10 0 5 10 15 20 25 Pt [dBm] Figure 7.2: Outage probability for target rate of Rtarget = 500 Mbit/s for hybrid RF/FSO system with joint BICM (proposed), independent BICM with CSI at the transmitter, and independent BICM without CSI at the transmitter, respectively. 16– QAM in RF link and L = 2.5 km. 43 possible with the optimum power assignment in (4.25) are negligible. Fig. 7.3 depicts the FER vs. Pt for hybrid RF/FSO systems with joint BICM (16–QAM, Rc = 1/2, q = 7, N = 6, cf. Table 5.1) and RC–SD (uncoded RF and uncoded FSO link), respectively. Clear air fog and moderate rain are considered for a link distance of L = 2.5 km. For comparison, we also show results for an uncoded FSO system, which yields a better performance than the uncoded RF system (not shown here). Note that all considered schemes have a data rate of 1 Gbit/s. Fig. 7.3 shows that joint BICM (solid lines) and RC–SD (dashed lines) have the same diversity orders and both schemes achieve a diversity gain compared to an uncoded FSO system (dotted lines) by exploiting the RF link. However, joint BICM substantially outperforms RC–SD. For example, at FER = 10−4 joint BICM yields gains of 4.3 dB,4.0 dB and 4.0 dB compared to RC–SD for clear air fog, and moderate rain, respectively. Furthermore, while joint BICM yields a performance similar to or better than that of uncoded FSO for all considered scenarios, RC–SD suffers from a 3 dB performance loss compared to uncoded FSO for moderate rain. For moderate rain the SNR in the FSO channel is substantially better than the SNR in the RF channel. Therefore, the simple RC–SD scheme chooses almost always the packets transmitted over the FSO link and the power assigned to the RF link is wasted. For Fig. 7.4 we adopted the same parameters as for Fig. 7.3. However, now the maximum RF power P2,max was limited to 10 dBm, which is the legal limit in Japan and Australia, cf. [5], and for the proposed joint BICM scheme the RF and FSO transmit powers were assigned using the assignment scheme outlined in Section 4.3. For RC–SD the powers were assigned in a similar fashion. Note that the maximum FSO output power is not important here as it lies outside the considered Pt range, cf. [10]. The effect of RF power limitation is clearly visible in Fig. 7.4 as the FER curves for the two considered hybrid RF/FSO schemes are asymptotically parallel to the FER curves of the uncoded FSO link. Nevertheless, at FER = 10−3 the proposed joint BICM scheme yields a gain of 9.5 dB and 2.7 dB compared to 44 0 10 −1 10 −2 FER 10 −3 10 Clear Air Haze Rain Joint BICM (proposed) RC−SD FSO −4 10 0 5 10 15 20 25 30 35 40 Pt [dBm] Figure 7.3: FER of hybrid RF/FSO system with joint BICM (16—QAM, Rc = 1/2, q = 7, N = 6, equal power for RF and FSO sub–channels, cf. Table 5.1) and RC– SD (uncoded RF and uncoded FSO sub–channels), respectively, and uncoded FSO. L = 2.5 km and data rate of 1 Gbit/s. 45 0 10 −1 10 −2 FER 10 −3 10 Clear Air Haze Rain Joint BICM (proposed) RC−SD FSO −4 10 0 5 10 15 20 25 30 35 40 Pt [dBm] Figure 7.4: FER of hybrid RF/FSO system with joint BICM (16—QAM, Rc = 1/2, q = 7, N = 6, cf. Table 5.1) and RC–SD (uncoded RF and uncoded FSO sub–channels), respectively, and uncoded FSO. P2,max = 10 dBm, L = 2.5 km, and data rate of 1 Gbit/s. the FSO link for clear air and moderate rain, respectively. It is also interesting to note that for moderate rain both joint BICM and RC–SD benefit from the power limitation in the RF link as more power is assigned to the FSO link which has a higher SNR. Therefore, the respective FER curves in Fig. ?? show an improved performance compared to those in Fig. 7.3. In Fig. 7.5, we show the FER of hybrid RF/FSO systems with joint BICM (16–QAM, Rc = 1/4, q = 3, N = 12, cf. Table 5.1) and RC–SD and the FER of pure FSO and RF systems. For RC–SD and the RF system 16–QAM is adopted for the RF link. For RC–SD we employed independent BICM for the RF and FSO sub– channels using the code with Rc = 1/2 and q = 3 from [29, Tab. 8.2-1]. The same code was also used for the pure FSO and RF systems. The data rate 46 of all four considered systems is 500 Mbit/s. In Fig. 7.5, performance under light fog and heavy rain is considered for L = 2.5 km. In light fog, the FSO system fails completely because of the high link attenuation. The RF system outperforms both hybrid systems by about 3 dB as the power assigned to the FSO link is essentially wasted in this case. Note that, in the considered Pt range, the FER curves for moderate and heavy fog are practically identical to those for light fog. In contrast, in heavy rain, the RF system fails and the FSO system has the best performance. The hybrid RF/FSO systems still achieve an acceptable performance and perform about 3 dB worse than the FSO system. Fig. 7.5 shows that because of its ability to provide protection against low SNR in one sub–channel by exploiting the high SNR in the other sub–channel, the proposed joint BICM approach offers a robust performance (without requiring CSI at the transmitter) even under extreme weather conditions. Fig. 7.6 shows the FER of hybrid RF/FSO systems with joint BICM and RC–SD with 4–QAM in the RF link, respectively, and the FER of an FSO system. Haze and moderate rain scenarios are considered for L = 2.5 km. If full diversity is desired, the maximum code rate for joint BICM is Rc = 1/3, cf. Section 4.4. For joint BICM we adopt the code with Rc = 1/3, q = 4, and N = 9 from Table 5.1. For RC–SD the maximum free distance code with Rc = 1/2 and q = 4 from [29, Tab. 8.2-1] was used for the FSO link, while the RF link remained uncoded. The same code was also used for the FSO system. The data rate of all considered schemes is 500 Mbit/s. Fig. 7.6 confirms the superior performance of the proposed joint BICM approach as compared to RC–SD and FSO also for the case where the RF link supports a lower data rate than the FSO link. Figure 7.7 shows the FER of hybrid FSO/RF with joint BICM (RF: 16– QAM, FSO: 16–QAM, Rc = 1/2) for code with different constraint length K and different assignment pattern N, according to the Table 5.1, for a distance of L = 2.5 Km, and a signal rate of 1 Gbit/s, under haze condition. Here we compare K from 3 to 7, with N between 4 and 6. Also, we compare the 47 0 10 −1 10 −2 FER 10 Fog −3 10 Moderate Rain Heavy Rain Joint BICM (proposed) RC−SD RF FSO −4 10 0 5 10 15 20 25 30 35 40 Pt [dBm] Figure 7.5: FER of hybrid RF/FSO system with joint BICM (Rc = 1/4, q = 3, N = 12, equal power for RF and FSO sub–channels, cf. Table 5.1) and RC–SD (16– QAM, Rc = 1/2, q = 3, in RF and FSO sub–channels), respectively, and coded RF (Rc = 1/2, q = 3) and coded FSO (Rc = 1/2, q = 3). L = 2.5 km and data rate of 500 Mbit/s. 48 0 10 −1 10 −2 FER 10 −3 10 Clear Air Haze Rain Joint BICM (proposed) RC−SD FSO −4 10 0 5 10 15 20 25 30 35 40 Pt [dBm] Figure 7.6: FER of hybrid RF/FSO system with joint BICM (4–QAM, Rc = 1/3, q = 4, N = 9, equal power for RF and FSO sub–channels, cf. Table 5.1) and RC– SD (uncoded RF sub–channel with 4–QAM, Rc = 1/2, q = 4 in FSO sub–channel), respectively, and coded FSO (Rc = 1/2, q = 4). L = 2.5 km and data rate of 1 Gbit/s. 49 0 10 −1 10 −2 FER 10 K=3,N=4,Power Unlimited K=3, N=6, Power Unlimited −3 10 K=4, N=4, Power Unlimited K=4,N=6,Power Unlimited K=5,N=4 Power Unlimited K=5,N=6 Power Unlimited K=7,N=6,Optimal Assignment K=7,N=6,Power unlimited −4 10 10 15 20 25 Pt [dBm] Figure 7.7: FER of hybrid RF/FSO system with joint BICM with different combination of constraint length, assignment pattern length and power limitation (coded RF subchannel with 16–QAM Rc = 1/2 and coded FSO subchannel Rc = 1/2 ). L = 2.5 km and data rate of 1 Gbit/s cases between 2 power assignment schemes Optimal power ssignment vs. equal power assignment) for K = 7, N = 6. From the figure we observe that the asymptotic FERs for different scenarios are different. The length of assignment pattern has less effect than constraint length in terms of performance improvement. It is also interesting to see that two power assignment scenarios have similar performances asymptotically. 50 Chapter 8 Conclusions and Future Work This chapter concludes the thesis with some general comments on the joint BICM hybrid system proposed in this work, followed by a discussion of possible future work for further investigation of the general our design. 8.1 Conclusions In this thesis, a novel architecture for hybrid RF/FSO systems has been proposed. Joint BICM is applied to the bit streams transmitted over the RF link and the FSO link resulting in a robust hybrid system. The provided asymptotic performance analysis has shown that properly designed CCs and APs can exploit the diversity offered by the independent RF and FSO sub–channels. An efficient search procedure for corresponding CCs and APs has been provided. In this context, an upper bound for the code rate has been derived which guarantees the existence of codes achieving full diversity. Furthermore, the optimum power assignment if the maximum RF and/or FSO transmit power is limited has been discussed, and it has been shown that the RF link can be beneficial even if the maximum RF transmit power is limited. FER and cut–off rate results for typical system and channel parameters have confirmed the excellent performance of the proposed hybrid RF/FSO scheme for various weather conditions. 51 8.2 Recommendations for Future Work We believe that the research work we initiated here on hybrid FSO/RF systems only scratches as the tip of the iceberg and many important questions remain to be answered. We list some recommendations for future work as follows: • In the system model, we assume that the fading models are GammaGamma distribution for the FSO channel and Rice fading for the RF channel. The union bound of CEP is derived based on these assumptions. However, it will be interesting to also take into account other fading models for both channels (for example, log-normal model for FSO for weak turbulence case). New union bounds for the CEP and new code/assignment pattern search algorithms should be derived for these channel models. • Furthermore, we only discuss the case when a single transmitter and a single receiver is presented for each channel. We know multiple input multiple output (MIMO) systems can be applied to both systems. It would be interesting to investigate how our proposed design would change if MIMO was taken in to consideration . • In the thesis, we use convolutional codes for channel coding. It is shown that turbo coding [48] has a performance approaching the Shannon limit. It is also suggested in [19] that BICM can be designed using binary parallel turbo coding. The high decoding complexity of turbo coding can be alleviated by using iterative decoding. It would be interesting to utilize turbo coding with BICM in our system and compare the CEP, the capacity, and the FER performance with the results with convolutional coding. 52 Bibliography [1] H. Willeband and B. Ghuman. Free Space Optics: Enabling Optical Connectivity in Today’s Network. Sams publishing, Indianapolis, 2002. [2] D. Kedar and S. Arnon. Urban Optical Wireless Communication Networks: The Main Challenges and Possible Solutions. IEEE Communications Magazine, 42(5):2–7, 2004. [3] Free Space Optics. http://www.free-space-optics.org. [4] F. Giannetti, M. Luise, and R. Reggiannini. Mobile and Personal Communications in the 60 GHz Band: A Survey. Wireless Personal Communications, 10:207–243, July 1999. [5] N. Guo, R. Qiu, S. Mo, and K. Takahashi. 60-GHz Millimeter-Wave Radio: Principle, Technology, and New Results. EURASIP J. Wireless Commun. and Networking, pages 1–7, January 2007. [6] L. Andrews, R. Phillips, and C. Hopen. Laser Beam Scintillation with Applications. EURASIP J. Wireless Commun. and Networking, 2001. [7] M. Al-Habash, L. Andrews, and R. Phillips. Mathematical Model for the Irradiance Probability Density Function of a Laser Beam Propagating Through Turbulent Media. Optical Engineering, 40(8):1554–1562, August 2001. 53 [8] P. Vasconcelos and L. Correia. Fading Characterization of the Mobile Radio Channel at the Millimetre Waveband. In Proc. IEEE Vehicular Technol. Conf. (VTC), pages 999–1003, May 1997. [9] J. Sch¨other. Fading Characterization of the Mobile Radio Channel at the Millimetre Waveband. Ist-2001-32686, Broadway, 2001. [10] I. Kim and E. Korevaar. Availability of Free Space Optics (FSO) and Hybrid FSO/RF Systems. Optical Wireless Communications IV, August 2001. [11] A. Akbulut, G. Ilk, and F. Ari. Design, Availability and Reliability Analysis on an Experimental Outdoor FSO/RF Communication System. In Proc. Intern. Transparent Optical Networks Conf. (ICTON), pages 403– 406, July 2005. [12] J. Derenick, C. Thorne, and J. Spletzer. On the Deployment of a Hybrid Free-space Optic/Radio Frequency (FSO/RF) Mobile Ad-hoc Network. In Proc. IEEE Intern. Conf. Commun. (ICC), pages 3990–3996, August 2005. [13] S. Bloom and W. Hartley. The Last-Mile Solution: Hybrid FSO Radio, May 2002. White Paper, AirFiber,. Inc., 802-0008-000 M-A1. [14] E. Leitgeb and M. Gebhart. High Availability of Hybrid Wireless Networks. In Proc. SPIE, volume 5465, pages 238–249, 2004. [15] H. Wu, B. Hamzeh, and M. Kavehrad. Achieving Carrier Class Availability of FSO Link via a Complementary RF Link. In Proc. Asilomar Conf. Signals, Systems and Computers, volume 2, pages 1483–1487, November 2004. [16] S. Vangala and H. Pishro-Nik. A Highly Reliable FSO/RF Communication System Using Efficient Codes. In Proc. IEEE Global Telecommun. Conf. (Globecom), pages 2232–2236, Washington, DC, November 2007. 54 [17] E. Malkam¨aki and H. Leib. Coded Diversity on Block-Fading Channels. IEEE Trans. Inform. Theory, 45:771–781, March 1999. [18] R. Knopp and P. Humblet. On Coding for Block Fading Channels. IEEE Trans. Inform. Theory, 46:189–205, January 2000. [19] G. Caire, G. Taricco, and E. Biglieri. Bit-Interleaved Coded Modulation. IEEE Trans. Inform. Theory, 44:927–946, May 1998. [20] A. Guillen i Fabregas and G. Caire. Coded Modulation in the BlockFading Channel: Coding Theorems and Code Construction. IEEE Trans. Inform. Theory, 52:91–114, January 2006. [21] Y. Leung, S. Wilson, and J. Ketchum. Multifrequency Trellis Coding with Low Delay for Fading Channels. IEEE Trans. Commun., 41:1450–1459, October 1993. [22] M. Chiani, A. Conti, and V. Tralli. Further Results on Convolutional Code Search for Block-Fading Channels. IEEE Trans. Inform. Theory, 50:1312–1318, June 2004. [23] M. Uysal, J. Li, and M. Yu. Error Rate Performance Analysis of Coded Free-Space Optical Links over Gamma-Gamma Atmospheric Turbulence Channels. IEEE Trans. Wireless Commun., 5(6):1229–1233, June 2006. [24] E. Bayaki, R. Schober, and R. Mallik. Performance Analysis of FreeSpace Optical Systems in Gamma-Gamma Fading. In Accepted for Presentation at the IEEE Global Telecommun. Conf. (Globecom), [Online] http://www.ece.ubc.ca/∼ehsanb/GC08.pdf, 2008. [25] S. Muhammad, P. K¨ohldorfer, and E. Leitgeb. Channel Modeling for Terrestrial Free Space Optical Links. In Proc. IEEE Intern. Conf. Transparent Optical Networks (ICTON), pages 407–410, July 2005. 55 [26] I. Kim, B. McArthur, and E. Korevaar. Comparison of Laser Beam Propagation at 785 nm and 1550 nm in Fog and Haze for Optical Wireless Communications. In Proc. SPIE, volume 4214, pages 26–37, 2001. [27] X. Zhu and J. M. Kahn. Free-Space Optical Communication Through Atmospheric Turbulence Channels. IEEE Trans. Commun., 50(8):1293– 1300, August 2002. [28] J. Goodman. Statistical optics. John Wiley and Son, 1985. [29] J. Proakis. Digital Communications. McGraw–Hill, New York, 2000. [30] Lianne Caetano and Sheung Li. Benefits of 60 GHz: Right Frequency, Right Time. SiBEAM Inc., November 2005. [31] ISM band. http://en.wikipedia.org/wiki/ISMband. [32] IEEE 802.11. http://en.wikipedia.org/wiki/802.11. [33] HIPERLAN. http://en.wikipedia.org/wiki/HIPERLAN. [34] P. Smulders. Exploiting The 60 GHz Band For Local Wireless Multimedia Access: Prospects and Future Directions . IEEE Communication magazine, 40(1):140–147, January 2002. [35] K. Gosseand K. Taffin and N. Whinnett. Propagation and Deployment of Mobile Broadband Systems in The 60 GHz Band. Technical report, Motorola CRM internal, March 1999. v.1.0. [36] P. Lusina. Ka-Broadband Satellite Communication using Cyclostationary Parabolic Beamforming. Master’s thesis, Queen’s University, Kingston, Ontario, Canada, August 1997. [37] S. Vangala and H. Pishro-Nik. Optimal Hybrid RF-wireless Optical Communication for Maximum Efficiency and Reliability. In Information sciences and systems, 2007. CISS ’07. 41st Annual Conference on , pages 684–689, March 2007. 56 [38] S. Milner and C. Davis. Hybrid Free Space Optical/RF Networks for Tactical Operations. In Proc. IEEE Military Commun. Conf. (MILCOM), pages 409–415, Monterey, CA, October 2004. [39] H. Pishro-Nik and H. F. Fekri. Results on Non-uniform Error Correction Using Low-Density Parity Check Codes. In Proc. IEEE global telecommun. conf. (Globecom) , pages 2041–2045, December 2003. [40] A. A. Farid and S. Hranilovic. Outage Capacity Optimization for FreeSpace Optical Links With Pointing Errors. IEEE J. Lightwave Techn., pages 1702–1710, July 2007. [41] W. Sch¨afer. Channel Modelling of Short-Range Radio Links at 60 GHz for Mobile Intervehicle Communication. In Proc. IEEE Veh. Techn. Conf. (VTC), pages 314–319, St. Louis, USA, May 1991. [42] Z. Wang and G. Giannakis. A Simple and General Parameterization Quantifying Performance in Fading Channels. IEEE Trans. Commun., 51:1389–1398, August 2003. [43] I. Gradshteyn and I. Ryzhik. Table of Integrals, Series, and Products. Academic Press, New York, 2000. [44] S. Yong and C.-C. Chong. An Overview of Multigigabit Wireless through Millimeter Wave Technology: Potentials and Technical Challenges. EURASIP J. Wireless Commun. and Networking, pages 1–10, January 2007. [45] T. Moon and W. Stirling. Mathematical Methods and Algorithms for Signal Processing. Prentice Hall, New York, 2000. [46] J. Wolf and A. Viterbi. On the Weight Distribution of Linear Block Codes Formed From Convolutional Codes. IEEE Trans. Commun., 44:1049– 1051, September 1996. 57 [47] I. Sason and I. Goldenberg. Coding for Parallel Channels: Gallager Bounds and Application to Turbo-Like Codes. IEEE Trans. Inform. Theory, 53:2394–2428, July 2007. [48] C. Berrou and A. Glavieux. Near Optimum Error Correcting Coding and Decoding: Turbo-Codes. IEEE Trans. Commun., 44(10):1261–1271, October 1996. 58
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Bit-interleaved coded modulation for hybrid rf/fso...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Bit-interleaved coded modulation for hybrid rf/fso systems He, Xiaohui 2008
pdf
Page Metadata
Item Metadata
Title | Bit-interleaved coded modulation for hybrid rf/fso systems |
Creator |
He, Xiaohui |
Publisher | University of British Columbia |
Date Issued | 2008 |
Description | In this thesis, we propose a novel architecture for hybrid radio frequency (RF)/free–space optics (FSO) wireless systems. Hybrid RF/FSO systems are attractive since the RF and FSO sub–systems are affected differently by weather and fading phenomena. We give a thorough introduction to the RF and FSO technology, respectively. The state of the art of hybrid RF/FSO systems is reviewed. We show that a hybrid system robust to different weather conditions is obtained by joint bit–interleaved coded modulation (BICM) of the bit streams transmitted over the RF and FSO sub–channels. An asymptotic performance analysis reveals that a properly designed convolutional code can exploit the diversity offered by the independent sub–channels. Furthermore, we develop code design and power assignment criteria and provide an efficient code search procedure. The cut–off rate of the proposed hybrid system is also derived and compared to that of hybrid systems with perfect channel state information at the transmitter. Simulation results show that hybrid RF/FSO systems with BICM outperform previously proposed hybrid systems employing a simple repetition code and selection diversity. |
Extent | 405005 bytes |
Subject |
FSO RF BICM Physical layer Channel coding Cut-off rate Pairwise Error Probability Codeword Error Probability Ricean fading Gamma-Gamma fading |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2008-12-11 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivatives 4.0 International |
DOI | 10.14288/1.0064860 |
URI | http://hdl.handle.net/2429/2868 |
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 |
Graduation Date | 2009-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/4.0/ |
Aggregated Source Repository | DSpace |
Download
- Media
- 24-ubc_2009_spring_he_xiaohui.pdf [ 395.51kB ]
- Metadata
- JSON: 24-1.0064860.json
- JSON-LD: 24-1.0064860-ld.json
- RDF/XML (Pretty): 24-1.0064860-rdf.xml
- RDF/JSON: 24-1.0064860-rdf.json
- Turtle: 24-1.0064860-turtle.txt
- N-Triples: 24-1.0064860-rdf-ntriples.txt
- Original Record: 24-1.0064860-source.json
- Full Text
- 24-1.0064860-fulltext.txt
- Citation
- 24-1.0064860.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-0064860/manifest