Outage performance of MIMO FSO links over strong ... - OSA Publishing

Outage performance of MIMO FSO links over strong turbulence and misalignment fading channels Antonio Garc´ıa-Zambrana,1,∗ Carmen Castillo-Vázquez,2 and Beatriz Castillo-Vázquez1 1 Department 2 Department

of Communications Engineering, University of Málaga, E-29071 Málaga, Spain of Statistics and Operations Research, University of Málaga, E-29071 Málaga, Spain ∗ [email protected]

Abstract: Atmospheric turbulence produces fluctuations in the irradiance of the transmitted optical beam, which is known as atmospheric scintillation, severely degrading the performance over free-space optical (FSO) links. Additionally, since FSO systems are usually installed on high buildings, building sway causes vibrations in the transmitted beam, leading to an unsuitable alignment between transmitter and receiver and, hence, a greater deterioration in performance. In this paper, the outage probability as a performance measure for multiple-input/multiple-output (MIMO) FSO communication systems with intensity modulation and direct detection (IM/DD) over strong atmospheric turbulence channels with pointing errors is analyzed. Novel closed-form expressions for the outage probability as well as their corresponding asymptotic expressions are presented when the irradiance of the transmitted optical beam is susceptible to either strong turbulence conditions, following a negative exponential distribution, and pointing error effects, following a misalignment fading model where the effect of beam width, detector size and jitter variance is considered. Obtained results show that the diversity order is independent of the pointing error when the equivalent beam radius at the receiver is at least twice the value of the pointing error displacement standard deviation at the receiver. Simulation results are further demonstrated to confirm the analytical results. Additionally, since proper FSO transmission requires transmitters with accurate control of their beamwidth, asymptotic expressions here obtained for different diversity techniques are used to find the optimum beamwidth that minimizes the outage performance. © 2011 Optical Society of America OCIS codes: (010.1330) Atmospheric turbulence; (060.2605) Free-space optical communication; (060.4510) Optical communications.

References and links 1. J. M. Kahn and J. R. Barry, “Wireless infrared communications,” Proc. IEEE 85, 265–298 (1997). 2. L. B. Stotts, L. C. Andrews, P. C. Cherry, J. J. Foshee, P. J. Kolodzy, W. K. McIntire, M. Northcott, R. L. Phillips, H. A. Pike, B. Stadler, and D. W. Young, “Hybrid optical RF airborne communications,” Proc. IEEE 97(6), 1109–1127 (2009). 3. W. Lim, C. Yun, and K. Kim, “BER performance analysis of radio over free-space optical systems considering laser phase noise under gamma-gamma turbulence channels,” Opt. Express 17(6), 4479–4484 (2009). 4. L. Andrews, R. Phillips, and C. Hopen, Laser Beam Scintillation with Applications (SPIE Press, 2001).

#146253 - $15.00 USD

(C) 2011 OSA

Received 19 Apr 2011; revised 6 Jun 2011; accepted 10 Jun 2011; published 28 Jun 2011

4 July 2011 / Vol. 19, No. 14 / OPTICS EXPRESS 13480

5. X. Zhu and J. M. Kahn, “Free-space optical communication through atmospheric turbulence channels,” IEEE Trans. Commun. 50(8), 1293–1300 (2002). 6. E. J. Lee and V. W. S. Chan, “Part 1: optical communication over the clear turbulent atmospheric channel using diversity,” IEEE J. Sel. Areas Commun. 22(9), 1896–1906 (2004). 7. I. B. Djordjevic, S. Denic, J. Anguita, B. Vasic, and M. Neifeld, “LDPC-coded MIMO optical communication over the atmospheric turbulence channel,” J. Lightwave Technol. 26(5), 478–487 (2008). 8. M. Simon and V. Vilnrotter, “Alamouti-type space-time coding for free-space optical communication with direct detection,” IEEE Trans. Wireless Commun. 4(1), 35–39 (2005). 9. C. Abou-Rjeily and W. Fawaz, “Space-time codes for MIMO ultra-wideband communications and MIMO freespace optical communications with PPM,” IEEE J. Sel. Areas Commun. 26(6), 938–947 (2008). 10. T. A. Tsiftsis, H. G. Sandalidis, G. K. Karagiannidis, and M. Uysal, “Optical wireless links with spatial diversity over strong atmospheric turbulence channels,” IEEE Trans. Wireless Commun. 8(2), 951–957 (2009). 11. E. Bayaki and R. Schober, “On space-time coding for free-space optical systems,” IEEE Trans. Commun. 58(1), 58–62 (2010). 12. A. Garcia-Zambrana, C. Castillo-Vazquez, B. Castillo-Vazquez, and A. Hiniesta-Gomez, “Selection transmit diversity for FSO links over strong atmospheric turbulence channels,” IEEE Photon. Technol. Lett. 21(14), 1017– 1019 (2009). 13. A. Garc´ıa-Zambrana, C. Castillo-Vázquez, and B. Castillo-Vázquez, “Space-time trellis coding with transmit laser selection for FSO links over strong atmospheric turbulence channels,” Opt. Express 18(6), 5356–5366 (2010). 14. S. Arnon, “Optimization of urban optical wireless communication systems,” IEEE Trans. Wireless Commun. 2(4), 626–629 (2003). 15. S. Arnon, “Effects of atmospheric turbulence and building sway on optical wireless-communication systems,” Opt. Lett. 28(2), 129–131 (2003). 16. A. A. Farid and S. Hranilovic, “Outage capacity optimization for free-space optical links with pointing errors,” J. Lightwave Technol. 25(7), 1702–1710 (2007). 17. H. G. Sandalidis, T. A. Tsiftsis, G. K. Karagiannidis, and M. Uysal, “BER performance of FSO links over strong atmospheric turbulence channels with pointing errors,” IEEE Commun. Lett. 12(1), 44–46 (2008). 18. H. G. Sandalidis, “Coded free-space optical links over strong turbulence and misalignment fading channels,” IEEE Trans. Commun. 59(3), 669–674 (2011). 19. D. K. Borah and D. G. Voelz, “Pointing error effects on free-space optical communication links in the presence of atmospheric turbulence,” J. Lightwave Technol. 27(18), 3965–3973 (2009). 20. W. Gappmair, S. Hranilovic, and E. Leitgeb, “Performance of PPM on terrestrial FSO links with turbulence and pointing errors,” IEEE Commun. Lett. 14(5), 468–470 (2010). 21. A. A. Farid and S. Hranilovic, “Diversity gains for MIMO wireless optical intensity channels with atmospheric fading and misalignment,” in Proc. IEEE GLOBECOM Workshops (GC Wkshps, 2010), pp. 1015–1019. 22. M. A. Al-Habash, L. C. Andrews, and R. L. Phillips, “Mathematical model for the irradiance probability density function of a laser beam propagating through turbulent media,” Opt. Eng. 40, 1554–1562 (2001). 23. A. Garc´ıa-Zambrana, C. Castillo-Vázquez, and B. Castillo-Vázquez, “Rate-adaptive FSO links over atmospheric turbulence channels by jointly using repetition coding and silence periods,” Opt. Express 18(24), 25422–25440 (2010). 24. B. Castillo-Vazquez, A. Garcia-Zambrana, and C. Castillo-Vazquez, “Closed-Form BER expression for FSO links with transmit laser selection over exponential atmospheric turbulence channels,” Electron. Lett. 45(23), 1185–1187 (2009). 25. M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th ed. (Dover, 1970). 26. M. K. Simon and M.-S. Alouini, Digital Communications over Fading Channels, 2nd ed. (Wiley-IEEE Press, 2005). 27. I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products, 7th ed. (Academic Press Inc., 2007). 28. Z. Wang and G. B. Giannakis, “A simple and general parameterization quantifying performance in fading channels,” IEEE Trans. Commun. 51(8), 1389–1398 (2003). 29. S. M. Navidpour, M. Uysal, and M. Kavehrad, “BER performance of free-space optical transmission with spatial diversity,” IEEE Trans. Wireless Commun. 6(8), 2813–2819 (2007). 30. H. A. David and H. N. Nagaraja, Order Statistics, 3rd ed. (John Wiley and Sons Inc., 2003). 31. H. G. Sandalidis, “Optimization models for misalignment fading mitigation in optical wireless links,” IEEE Commun. Lett. 12(5), 395–397 (2008). 32. Wolfram Research, Inc., Mathematica, version 8.0.1. ed. (Wolfram Research, Inc., Champaign, Illinois, 2011).

#146253 - $15.00 USD

(C) 2011 OSA



1.

Introduction

Optical wireless communications using intensity modulation and direct detection (IM/DD) can provide high-speed links for a variety of applications [1], providing an unregulated spectral segment and high security. Recently, the use of atmospheric free-space optical (FSO) transmission is being specially interesting to solve the “last mile” problem, as well as a supplement to radio-frequency (RF) links [2, 3]. However, atmospheric turbulence produces fluctuations in the irradiance of the transmitted optical beam, which is known as atmospheric scintillation, severely degrading the link performance [4, 5]. Additionally, since FSO systems are usually installed on high buildings, building sway causes vibrations in the transmitted beam, leading to an unsuitable alignment between transmitter and receiver and, hence, a greater deterioration in performance. Error control coding as well as diversity techniques can be used over FSO links to mitigate turbulence-induced fading [6–11]. In [12, 13], selection transmit diversity is proposed for FSO links over strong turbulence channels, where the transmit diversity technique based on the selection of the optical path with a greater value of irradiance has shown to be able to extract full diversity as well as providing better performance compared to general FSO spacetime codes (STCs) designs, such as conventional orthogonal space-time block codes (OSTBCs) and repetition codes (RCs). The combined effect of atmospheric and misalignment fading is analyzed in the case of single-input/single-output (SISO) FSO channels [14, 15]. In [16], the effects of atmospheric turbulence and misalignment considering aperture average effect were considered to study the outage capacity for SISO links. In [17, 18] the error rate performance for uncoded and coded FSO links over strong turbulence and misalignment fading channels is studied. The capacity calculation and the analysis of error rate performance for FSO links over log-normal and gamma-gamma turbulence and misalignment fading channels is presented in [19]. In [20] the error rate performance of PPM on terrestrial FSO Links with gamma-gamma turbulence and pointing errors is analyzed. However, to the best of our knowledge, the performance of MIMO FSO channels under the effects of atmospheric and misalignment fading has only been studied in [21], assuming repetition coding on the transmitter side and equal gain combining on the receiver side. In this work [21], the study of the outage probability and diversity gain has been considered for multiple-input/multiple-output (MIMO) FSO communication systems impaired by log-normal atmospheric turbulence and misalignment fading, showing that the diversity gain is independent of the number of transmitters and receivers, being conditioned only by misalignment parameters. In this paper, the outage probability as a performance measure for MIMO FSO communication systems using IM/DD over strong atmospheric turbulence channels with pointing errors is analyzed, assuming different configurations as repetition coding and transmit laser selection on the transmitter side and equal gain combining and selection combining on the receiver side. Novel closed-form expressions for the outage probability as well as their corresponding asymptotic expressions are presented when the irradiance of the transmitted optical beam is susceptible to either strong turbulence conditions, following a negative exponential distribution, and pointing error effects, following a misalignment fading model, as in [16–18, 21], where the effect of beam width, detector size and jitter variance is considered. In this strong turbulence FSO scenario, obtained results show that the diversity order is independent of the pointing error when the equivalent beam radius at the receiver is at least twice the value of the pointing error displacement standard deviation at the receiver, showing the same slope of the outage performance versus average signal-to-noise ratio (SNR) as in a similar FSO scenario where misalignment fading is not considered. However, different coding gain, i.e. the horizontal shift in the outage performance in the limit of large SNR, is achieved as a consequence of the severity of the pointing errors effects and MIMO configuration assumed. Additionally, a significant improvement in performance is demonstrated when MIMO FSO links based on transmit laser

#146253 - $15.00 USD

(C) 2011 OSA



selection with equal gain combining are adopted. Simulation results are further demonstrated to confirm the analytical results. Here, not only rectangular pulses are considered but also on-off keying (OOK) formats with any pulse shape, corroborating the advantage of using pulses with high peak-to-average optical power ratio (PAOPR). In contrast to [21], wherein it is concluded that the diversity gain for MIMO FSO communication systems impaired by log-normal atmospheric turbulence and misalignment fading is conditioned only by misalignment parameters, it is shown in this paper that the outage diversity of MIMO FSO links over strong turbulence and misalignment fading channels is independent of pointing errors when transmitters are adequately designed, providing an accurate control of their beamwidth in order to guarantee an equivalent beam radius at the receiver of at least twice the value of the pointing error displacement standard deviation. 2.

System and channel model

We adopt a multiple-input/multiple-output array based on L laser sources, assumed to be intensity-modulated only and all pointed towards a distant array of M photodetectors, assumed to be ideal noncoherent (direct-detection) receivers. The transmit and receive apertures are physically situated so that all transmitters are simultaneously observed by each receiver. The use of infrared technologies based on IM/DD links is considered, where the instantaneous current ylm (t) in the mth receiving photodetector corresponding to the information signal transmitted from the lth laser can be written as ylm (t) = η ilm (t)x(t) + zm (t)

(1)

where η is the detector responsivity, assumed hereinafter to be the unity, X x(t) represents the optical power supplied by the lth source and Ilm ilm (t) the equivalent real-valued fading gain (irradiance) through the optical channel between the lth transmit and the mth receive aperture. Additionaly, the fading experienced between source-detector pairs Ilm is assumed to be statistically independent. Zm zm (t) is assumed to include any front-end receiver thermal noise as well as shot noise caused by ambient light much stronger than the desired signal at the mth detector. In this case, the noise can usually be modeled to high accuracy as AWGN with zero mean and variance σm2 = N0 /2M, i.e. Zm ∼ N(0, N0 /2M), independent of the on/off state of the received bit. The factor M is consequence of the fact that the sum of the M receive aperture areas is the same as the aperture area of a system with no receive diversity, wherein the noise is usually modelled as AWGN with zero mean and variance σ 2 = N0 /2. This allows the systems to be compared fairly [6]. Since the transmitted signal is an intensity, X must satisfy ∀t x(t) ≥ 0. Due to eye and skin safety regulations, the average optical power is limited and, hence, the average amplitude of X is limited. The received electrical signal Ylm ylm (t), however, can assume negative amplitude values. We use Ylm , X, Ilm and Zm to denote random variables and ylm (t), x(t), ilm (t) and zm (t) their corresponding realizations. The irradiance is susceptible to either strong atmospheric turbulence conditions and pointing error effects. In this case, it is considered to be a product of two independent random variables, (a) (p) (a) (p) i.e. Ilm = Ilm Ilm , representing Ilm and Ilm the attenuation due to atmospheric turbulence and the attenuation due to geometric spread and pointing errors, respectively, between lth transmitter and mth receiver. Although the effects of turbulence and pointing are not strictly independent, for smaller jitter values they can be approximated as independent [19]. Considering a limiting case of strong turbulence conditions [4, 22, 23], the turbulence-induced fading is modelled as a multiplicative random process which follows the negative exponential distribution, whose probability density function (PDF) is given by f

#146253 - $15.00 USD

(C) 2011 OSA

(a)

Ilm

(i) = e−i , i ≥ 0.

(2)



This PDF has also been adopted in different works [8, 9, 24] to describe turbulence-induced fading, leading to an easier mathematical treatment. Regarding to the impact of pointing errors, we use the general model of misalignment fading given in [16] by Farid and Hranilovic, wherein the effect of beam width, detector size and jitter variance is considered. Assuming a Gaussian spatial intensity profile of beam waist radius, ωz , on the receiver plane at distance z from the (p) transmitter and a circular receive aperture of radius r, the PDF of Ilm is given by f

(p) Ilm

(i) =

ϕ2 ϕ2 A0

iϕ

2 −1

, 0 ≤ i ≤ A0

(3)

where ϕ = ωzeq /2σs is the ratio between the equivalent beam radius at the receiver and√ the pointing error displacement √ standard deviation (jitter) at the receiver, √ ωz2eq = ωz2 π erf(v)/2v exp(−v2 ), v = π r/ 2ωz , A0 = [erf(v)]2 and erf(·) is the error function [25, eqn. (7.1.1)]. Here, independent identical Gaussian distributions for the elevation and the horizontal displacement (sway) are considered, being σs2 the jitter variance at the receiver. Using the previous PDFs for turbulence and misalignment fading, a closed-form expression of the combined PDF of Ilm was derived in [18] as fIlm (i) =

ϕ2 ϕ2 A0

iϕ

2 −1

Γ(1 − ϕ 2 ,

i ), i ≥ 0 A0

(4)

where Γ(·, ·) is the upper incomplete Gamma function [25, eqn. (6.5.3)]. Here, it must be commented that the pointing error model in Eq. (3) is not general enough to be directly applied to FSO systems with more than one photodetector (i.e. M > 1) since it depends on locations of the transmitter and receiver, requiring to take into account the relative position for each receive aperture regarding initial pointing in undisturbed position. This represents a fixed pointing error called boresight. As previously indicated, we adopt a MIMO array wherein the transmit and receive apertures are physically situated so that all transmitters are simultaneously observed by each receiver. For the sake of simplicity, it is here assumed that the receivers are placed very closely together and beam radius is large enough so that the relative displacement for each receiver regarding the initial pointing can be considered negligible, and hence the pointing error model in Eq. (3) can be approximately applied. This assumption can be reinforced by the fact that strong turbulence channels are more robust to boresight [19]. Additionally, as shown in the appendix, we consider OOK formats with any pulse shape and reduced duty cycle, allowing the increase of the PAOPR parameter [12, 23]. 3.

Outage performance analysis

In this section, the outage probability as a performance measure is evaluated for different MIMO FSO configurations over strong atmospheric turbulence channels with pointing errors. Together with the average error rate, outage probability, Pout , is another standard performance criterion of communications systems operating over fading channels [26]. It is defined as the probability that the instantaneous combined SNR, γT , falls below a certain specified threshold, γth , which represents a protection value of the SNR above which the quality of the channel is satisfactory, i.e., Pout = Prob(γT ≤ γth ) =

γth

0

fγT (γT )d γT = FγT (γth )

(5)

where fγT (γT ) and FγT (γT ) are the PDF and the cumulative distribution function (CDF) of γT , respectively. Next, the outage performance of a SISO FSO system is firstly evaluated in order to be considered as a benchmark in the analysis of the remaining MIMO configurations, #146253 - $15.00 USD

(C) 2011 OSA



corresponding to a multiple-input/single-output (MISO) system with a L × 1 array, a singleinput/multiple-output (SIMO) system with a 1 × M array and, finally, a MIMO FSO system with a L × M array. 3.1.

SISO FSO system

In this subsection, the outage probability for the FSO system model previously presented with L = M = 1 is evaluated. According to Eq. (1) and the OOK signaling described in the appendix, a constellation of two equiprobable points in a one-dimensional space with an Euclidean distance of d, the statistical channel model corresponding to the SISO configuration can be written as (6) Y = XI + Z, X ∈ {0, d}, Z ∼ N(0, N0 /2). The received electrical SNR can be defined, as in [16], as

γTSISO (i) =

2 T ξ 1 d 2 2 4Popt b 2 i = i = γξ i2 2 N0 /2 N0

(7)

where γ represents the received electrical SNR in absence of turbulence when the classical rectangular pulse shape is adopted for OOK formats, Tb parameter is the bit period, Popt is the average optical power transmitted and ξ represents the square of the increment in Euclidean distance due to the use of a pulse shape of high PAOPR, as explained in a greater detail in the appendix. Using Eq. (5), the outage probability for a SISO FSO system can be written as γth 1 γth 2 2 (8) ) = FI Pout = Prob(γξ i ≤ γth ) = Prob(i ≤ ξγ ξ γ where the CDF of I, FI (i), can be easily derived from Eq. (4) as

i FI (i) = 1 − A0

ϕ 2

i , i≥0 ϕ 2 Γ −ϕ 2 , A0

(9)

by applying the differential relation in [25, eqn. (6.5.26)] and the fact that the series expansion corresponding to the upper incomplete Gamma function can be simplified by az + O(z2 )) [27, eqn. (8.354.2)], where Γ(·) is the well-known Γ(a, z) ∝ Γ(a) − (za /a)(1 − 1+a Gamma function. The results corresponding to this FSO scenario are illustrated in Fig. 1, where rectangular pulse shapes with ξ = 1 are used, assuming different values of normalized beamwidth, ωz /r, and normalized jitter, σs /r (i.e. ωz /r = 5 with σs /r = {1, 3, 4, 5} and ωz /r = 10 with σs /r = {1, 4, 7, 11}). Outage simulation results are furthermore included as a reference, demonstrating an excellent agreement with the analytical result in Eq. (8) for different misalignment fading conditions. Additionally, we can take advantage of this series expansion in order to quantify the outage probability at high SNR, showing that the asymptotic performance of this metric as a function of the average SNR is characterized by two parameters: the diversity and coding gains. Therefore, for large enough SNR (γ → ∞), the asymptotic outage performance can be expressed after some algebraic manipulations as . Pout =

#146253 - $15.00 USD

(C) 2011 OSA

ϕ4

γth 2 2 (1 − ϕ 2 ) A0 ξ γ

1/2

⎞ϕ 2 /2 ⎛

−2ϕ 2 Γ 1 − ϕ2 γ th ⎠ +⎝ . γ A20 ξ

(10)



0

10

−1

σs /r=10

10

−2

Outage probability

10

σs /r=5

−3

10

σs /r=4 −4

10

σs /r=1

σs /r=7

−5

10

no pointing errors

−6

10

ωz /r=5

−7

σs /r=3

ωz /r=10 Simulation Asymptotic analysis

10

−8

10

−9

10

0

50

100

150

200

Normalized average SNR, γ/γth (dB) Fig. 1. Probability of outage versus normalized average SNR in FSO IM/DD links over the exponential atmospheric turbulence channel with pointing errors, assuming different values of normalized beamwidth, ωz /r, and normalized jitter, σs /r. Here, an even greater simplification can be deduced, leading to different asympotic expressions depending on the value of ϕ , as follows . Pout =

ϕ4

γth 2 2 2 (1 − ϕ ) A0 ξ γ

1/2 ,

⎛

⎞ϕ 2 /2 2 2 2/ϕ γth ⎠ . ⎝Γ 1−ϕ Pout = , γ A20 ξ

ϕ >1

ϕ 1. Once this condition is satisfied, taking into account the coding gain in Eq. (11a), the impact of the pointing error effects translates into a coding gain disadvantage, D[dB], relative to strong atmospheric turbulence without misalignment fading given by ϕ2 . (12) D[dB] 20 log10 A0 (ϕ 2 − 1) According to this expression, it can be observed in Fig. 1 that coding gain disadvantages of 23.7, 34.4 and 42.7 decibels are achieved for values of (ωz /r, σs /r) = (5, 1), (ωz /r, σs /r) = (10, 1) and (ωz /r, σs /r) = (10, 4), respectively. 3.2.

MISO FSO system with a L × 1 array

In this subsection, the outage probability for the FSO system model previously presented where the information signal is transmitted via L apertures and received by only one aperture (i.e. M = 1) is evaluated. Here, two different transmit diversity strategies are considered: repetition coding (RC) and transmit laser selection (TLS). 3.2.1.

Transmit diversity using repetition coding

It is assumed that the same symbol is transmitted from the L transmitters at a given bit period [6, 10, 29]. According to Eq. (1), the statistical channel model corresponding to this MISO configuration can be written as Y =X

1 L ∑ Il + Z, L l=1

X ∈ {0, d},

Z ∼ N(0, N0 /2)

(13)

where Il represents the equivalent irradiance through the optical channel between the lth transmit aperture and the photodetector. Here, the division by L is considered so as to maintain the average optical power in the air at a constant level of Popt , being transmitted by each laser an average optical power of Popt /L. In this way, the total transmit power is the same as in a FSO system with no transmit diversity. The resulting received electrical SNR can be defined as

γTMISO−RC

1 = γξ 2 L

L

∑ il

l=1

2 = γξ

1 2 i L2 T

(14)

and, hence, the outage probability for a MISO FSO system with repetition coding can be written as 1 γth (15) Pout = FIT L ξ γ Next, in order to address the cumbersome statistical problem related to the CDF of the sum of L variates, the PDF in Eq. (4) is approximated by a single polynomial term as fIlm (i) ≈ aib , based on the fact that the asymptotic behavior of the system performance is dominated by the behavior of the PDF near the origin, i.e. fIlm (i) at i → 0 determines high SNR performance [28]. Using the series expansion corresponding to the upper incomplete Gamma function [27, eqn. (8.354.2)], different asympotic expressions for Eq. (4), depending on the value of ϕ , can be

#146253 - $15.00 USD

(C) 2011 OSA



written as

ϕ2 , ϕ >1 A0 (ϕ 2 − 1)

2 2 2Γ 1−ϕ fIlm (i) ≈ ϕ iϕ −1 , ϕ < 1 ϕ2 A0

fIlm (i) ≈

(16a) (16b)

Since the variates are independent, knowing that the resulting PDF of their sum is the convolution of their corresponding PDFs, obtained via single-sided Laplace and its inverse transforms, approximate expressions for the CDF, FIT (i), of the combined variate can be easily derived as L ϕ2 1 iL , ϕ > 1 (17a) FIT (i) ≈ Γ(L + 1) A0 (ϕ 2 − 1) L

L

Γ 1 − ϕ2 Γ 1 + ϕ2 2 iL ϕ , ϕ < 1 (17b) FIT (i) ≈ ϕ2 Γ (1 + Lϕ 2 ) A 0

and, then, the corresponding approximate expressions for the outage probability for a MISO FSO system with repetition coding can be obtained from Eq. (15) as follows L/2 ϕ4 L2 γth . , ϕ >1 (18a) Pout = Γ(L + 1)2/L A20 (1 − ϕ 2 )2 ξ γ ⎛ ⎞Lϕ 2 /2 2

2

2 2 2/ϕ Γ 1 − ϕ 2 2/ϕ γth ⎠ . ⎝L Γ 1+ϕ Pout = , ϕ 1

⎛

⎞Lϕ 2 /2 2 2 2/ϕ ϕ Γ 1 − γ . th ⎠ Pout = ⎝ , γ A20 ξ

ϕ 1 (22a) AT LS [dB] 10 log10 Γ(L + 1)2/L ⎛ 2/ϕ 2 ⎞

L2 Γ 1 + ϕ 2 ⎠, ϕ < 1 (22b) AT LS [dB] 10 log10 ⎝ 2 Γ (1 + Lϕ 2 )2/Lϕ As observed in Fig. 2, it can be deduced from expressions in Eq. (22) that the greater the number of transmit apertures, the better performance in terms of coding gain advantage. Additionally, this superiority is even more significant when ϕ < 1 and the value of ϕ is decreasing. 3.3.

SIMO FSO system with a 1 × M array

In this subsection, the outage probability for the FSO system model previously presented is evaluated when the information signal is transmitted via only one aperture (i.e. L = 1) and received by M apertures. Here, two different receive diversity strategies are considered: equal gain combining (EGC) and selection combining (SC). 3.3.1.

Receive diversity using equal gain combining

Following this combining technique, the receiver adds the information corresponding to M photodetectors and, hence, according to Eq. (1), the statistical channel model for this SIMO configuration can be written as Y =X

1 M ∑ Im + ZEGC , M m=1

X ∈ {0, d},

ZEGC ∼ N(0, N0 /2)

(23)

where Im represents the equivalent irradiance through the optical channel between the transmit aperture and the mth photodetector. Here, the division by M is considered to ensure that the area of the receive aperture in SISO links has the same size as in the sum of M receive aperture areas of SIMO links [6]. It can be noted that the resulting expression is equivalent to the expression in Eq. (13), obtained for MISO FSO links assuming RC at the transmitter side [10]. Therefore, we can interchange L by M in Eq. (15) and Eq. (18) to obtain the outage probability and its corresponding asympotic expressions for SIMO FSO links using equal gain combining.

#146253 - $15.00 USD

(C) 2011 OSA



3.3.2.

Receive diversity using selection combining

This type of combining is based on the selection of the receive aperture corresponding to the optical path with a greater value of fading gain (irradiance) [6,10]. In this way, in a similar way as previously presented for MISO FSO links with transmit laser selection, our SIMO system model can be considered as an equivalent SISO system model where the channel irradiance corresponding to the SC scheme, Imax , can be written as Imax = maxm=1,2,···M Im and, hence, the statistical channel model can be written as 1 Y = X Imax + ZSC , X ∈ {0, d}, ZSC ∼ N(0, N0 /2M) (24) M where M represents the same role as previously explained when EGC is adopted. It can be observed that the noise variance is given by N0 /2M since only one of the M receive apertures is used. Taking into account Eq. (7), the resulting received electrical SNR can be defined as 1 2 I (25) M max and, hence, the outage probability for a SIMO FSO system with selection combining can be written as M M γth M γth = FI . (26) Pout = FImax ξ γ ξ γ

γTSIMO−SC = γξ

In the same way, its corresponding asympotic expressions, depending on the value of ϕ , can be easily derived from Eq. (11) as M/2 ϕ4 γth . , ϕ >1 (27a) Pout = M (1 − ϕ 2 )2 A20 ξ γ ⎛ ⎞Mϕ 2 /2

2 2 2/ϕ γth ⎠ . ⎝ Γ 1−ϕ Pout = M , ϕ 1 and the superiority of both SC scheme or EGC scheme when ϕ < 1, depending on the values of M and ϕ . This improvement in performance of the EGC scheme when ϕ > 1 translates into a coding gain advantage, AEGC [dB], relative to the SC scheme as follows Γ(M + 1)2/M , ϕ >1 (28) AEGC [dB] 10 log10 M corroborating the fact that the greater the number of apertures, the higher coding gain advantage. Nonetheless, it can be noted that there is no difference between both diversity techniques for a value of M = 2. In a similar way, the superiority of SC scheme when ϕ < 1 translates into a coding gain advantage, ASC [dB], relative to the EGC scheme as follows ⎛ 2/ϕ 2 ⎞

MΓ 1 + ϕ 2 ⎠, ϕ < 1 ASC [dB] 10 log10 ⎝ (29) 2 Γ (1 + M ϕ 2 )2/M ϕ #146253 - $15.00 USD

(C) 2011 OSA



Nonetheless, it can be concluded from this expression that the superiority of the SC scheme is only valid for low values of ϕ and when a small number of receive apertures is considered, being only justified the adoption of this receive diversity technique in the context of very severe pointing errors. As depicted in Fig. 3 for a normalized beamwidth of ωz /r = 10, better performance for the EGC scheme is again achieved as the number of receive apertures increases and the value of the normalized jitter, σs /r, is lower and lower. So, the value of ϕ becomes higher and higher, and hence it is closer to the unity.

Coding gain advantage, ASC [dB]

3

ωz /r=10 2

1

0 −1

σs /r = 6 −2

σs /r = 7 σs /r = 8

−3

σs /r = 9

−4 2

4

6

8

10

12

14

16

number of receive apertures, M Fig. 3. Coding gain advantage for the SC scheme relative to the EGC scheme, as-

suming a normalized beamwidth of ωz /r = 10 and values of normalized jitter of 6, 7, 8 and 9, corresponding to values of ϕ of 0.83, 0.71, 0.62 and 0.55, respectively.

3.4.

MIMO FSO system with a L × M array

According to results in previous subsections, the outage probability for MIMO FSO IM/DD links is here evaluated when transmit laser selection and equal gain combining are the diversity strategies adopted in the transmit and receive side, respectively. The statistical channel model for this MIMO configuration can be written as Y =X

1 M ∑ Im + ZEGC , M m=1

X ∈ {0, d},

ZEGC ∼ N(0, N0 /2)

(30)

where Im represents the equivalent irradiance through the optical channels between the L transmit apertures and the mth photodetector, based on the selection of the optical path with a greater value of irradiance, i.e. Im = maxl=1,2,···L Ilm . The resulting received electrical SNR can be defined as 2 M 1 1 γTMIMO = γξ 2 ∑ im = γξ 2 i2T (31) M m=1 M

#146253 - $15.00 USD

(C) 2011 OSA



and, hence, the outage probability is given by Pout = FIT

M

1 γth ξ γ

(32)

where IT represents the sum of M variates following the Lth-order statistics distribution for L observations from the distribution in Eq. (4). As in previous subsections, in order to address the cumbersome statistical problem related to the CDF of the sum of these M variates and making use of Eq. (16), we use an approximated expression for the PDF of Im , which is given by L ϕ2 iL−1 , ϕ > 1 (33a) fIm (i) ≈ L (ϕ 2 − 1) A0 L 2 2 2 Γ(1 − ϕ ) iLϕ −1 , ϕ < 1 (33b) fIm (i) ≈ Lϕ ϕ2 A0 according to order statistics [30]. Since the variates are independent, knowing that the resulting PDF of their sum is the convolution of their corresponding PDFs, obtained via single-sided Laplace and its inverse transforms, approximate expressions for the CDF, FIT (i), of the combined variate can be easily derived as LM ϕ2 (LΓ(L))M FIT (i) ≈ iLM , ϕ > 1 (34a) Γ(LM + 1) (ϕ 2 − 1) A0 LM

M

Γ 1 − ϕ2 Γ Lϕ 2 + 1 2 iLMϕ , ϕ < 1 (34b) FIT (i) ≈ 2 2 ϕ Γ (LM ϕ + 1) A 0

and, then, the corresponding approximate expressions for the outage probability for a MIMO FSO system with transmit laser selection and equal gain combining can be obtained from Eq. (32) as follows

LM 2 M 2 Γ(L + 1)2/L ϕ4 γth , ϕ >1 Γ(LM + 1)2/LM (ϕ 2 − 1)2 A20 ξ γ ⎛ ⎞LM ϕ 2 2 2ϕ 2

2 2/Lϕ 2

2 2 Γ 1−ϕ M Γ Lϕ + 1 γth ⎠ . Pout = ⎝ , 2 2 2/LM ϕ γ A0 ξ Γ (LM ϕ 2 + 1)

. Pout =

(35a)

ϕ 1, being ϕ a main parameter to consider in order to optimize the outage performance. Once this condition is satisfied, comparing different diversity techniques and their corresponding asymptotic expressions, it can be noted that better performance is achieved when increasing the number of transmit apertures, i.e. L, instead of the number of receive apertures, i.e. M, in order to guarantee a same diversity order, as shown in Fig. 4 when assuming a normalized beamwidth of ωz /r = 5 with a normalized jitter of σs /r = 1 and a diversity order of 8 with values of {L, M} = {4, 2} and {L, M} = {2, 4}. Taking into account expressions in Eq. (35a) and Eq. (21a), this can be quantified as a coding gain disadvantage, DMIMO [dB], for the MIMO FSO system relative to the MISO FSO system using laser selection as follows MΓ(Od /M + 1)M/Od . (36) DMIMO [dB] 20 log10 Γ(Od + 1)1/Od where Od and M denote outage diversity and the number of receive apertures, respectively. In agreement with this expression, it can be seen in Fig. 4 that a difference of 2.13 decibels is achieved for a diversity order of 8 and values of M = 2 and M = 4 when a normalized beamwidth of ωz /r = 5 with a normalized jitter of σs /r = 1 is assumed. Lastly, since proper FSO transmission requires transmitters with accurate control of their beamwidth, the optimization procedure is finished by finding the optimum beamwidth, ωz /r, that gives the minimum #146253 - $15.00 USD

(C) 2011 OSA



Optimum normalized beamwidth, ωz /r

outage [16,17,31]. From Eq. (35a), it can be observed that this is equivalent to minimize the expression in Eq. (12) corresponding to the SISO configuration, being independent of the MIMO configuration adopted. In this way, the optimum beamwidth can be achieved using numerical optimization methods for different values of normalized jitter, σs /r [32]. Numerical results for the optimum beamwidth are illustrated in Fig. 5 for a range of values in the normalized jitter from σs /r = 1 to σs /r = 10 in discrete steps of 0.5 when the stochastic function minimizer simulated annealing is used. From this figure, it can be deduced that the outage optimization

30

Numerical evaluation Approximated expression

25 20 15 10 5 0 1

2

3

4

5

6

7

8

9

10

Normalized jitter, σs /r Fig. 5. Optimum normalized beamwidth versus normalized jitter, σs /r in FSO IM/DD links over the exponential atmospheric turbulence channel with pointing errors.

provides numerical results following a linear performance. This leads to easily obtain a firstdegree polynomial as follows

ωz /roptimum ≈ 2.85(σs /r − 1) + 2.6

(37)

considering the values of optimum beamwidth corresponding to normalized jitters of σs /r = 1 and σs /r = 10, respectively. As can be seen in this figure, it is clearly depicted that the approximative analytical expression remains very accurate to numerical results. The use of this expression is shown in Fig. 4, where results assuming the optimum beamwidth corresponding to values of jitter of σs /r = 1 and σs /r = 7 are also included for L = 4 and M = 2. As expected, a greater improvement in performance is achieved when σs /r = 7 compared to the improvement in performance corresponding to the normalized jitter of σs /r = 1 since not only coding gain but also diversity gain is obtained, making outage diversity performance independent of the impact of misalignment fading. 5.

Conclusions

In this paper, the outage probability as a performance measure for MIMO FSO communication systems using IM/DD over strong atmospheric turbulence channels with pointing errors is analyzed, assuming different configurations as repetition coding and transmit laser selection on the

#146253 - $15.00 USD

(C) 2011 OSA



transmitter side and equal gain combining and selection combining on the receiver side. Novel closed-form expressions for the outage probability as well as their corresponding asymptotic expressions are presented when the irradiance of the transmitted optical beam is susceptible to either strong turbulence conditions, following a negative exponential distribution, and pointing error effects, following a misalignment fading model where the effect of beam width, detector size and jitter variance is considered. In this strong turbulence FSO scenario, obtained results show that the diversity order is independent of the pointing error effects when the equivalent beam radius at the receiver is at least twice the value of the pointing error displacement standard deviation at the receiver, showing the same slope of the outage performance versus average SNR as in a similar FSO scenario where misalignment fading is not considered. However, different coding gain, i.e. the horizontal shift in the outage performance in the limit of large SNR, is achieved as a consequence of the severity of the pointing errors effects and MIMO configuration assumed. Additionally, a significant improvement in performance is demonstrated when MIMO FSO links based on transmit laser selection with equal gain combining are adopted. Simulation results are further demonstrated to confirm the analytical results. The asympotic expressions here obtained are especially useful for performance analysis of MIMO FSO systems, providing an easier understanding of performance limiting factors when operating under different misalignment fading conditions and various diversity techniques are adopted. From the relevant results here obtained, investigating in this FSO scenario the impact on the diversity order of different receive configurations, e.g., rectangular and circular arrays of receive apertures, wherein the displacement for each photodetector regarding the initial pointing is considered as well as the possible correlation among their corresponding shifting, is an interesting topic for future research. Appendix We consider OOK formats with any pulse shape and reduced duty cycle, allowing the increase Eg of the PAOPR parameter [12, 23]. A new basis function φ (t) is defined as φ (t) = g(t) where g(t) represents any normalized pulse shape satisfying the non-negativity constraint, with ∞ 2 g (t)dt is the electrical energy. In 0 ≤ g(t) ≤ 1 in the bit period and 0 otherwise, and Eg = −∞ this way, an expression for the optical intensity can be written as x(t) =

∞

∑

ak

k=−∞

2Tb Popt g (t − kTb ) G( f = 0)

(38)

where G( f = 0) represents the Fourier transform of g(t) evaluated at frequency f = 0, i.e. the area of the employed pulse shape, and Tb parameter is the bit period. The random variable (RV) ak follows a Bernoulli distribution with parameter 1/2, taking the values of 0 for the bit “0” (off pulse) and 1 for the bit “1” (on pulse). From this expression, it is easy to deduce that the points average optical power transmitted is Popt , defining a constellation of two equiprobable (x0 = 0 and x1 = d) in a one-dimensional space with an Euclidean distance of d = 2Popt Tb ξ where ξ = Tb Eg /G2 ( f = 0) represents the square of the increment in Euclidean distance due to the use of a pulse shape of high PAOPR, alternative to the classical rectangular pulse. Acknowledgments The authors would like to thank the anonymous reviewers for their useful comments that helped to improve the presentation of the paper. The authors are grateful for financial support from the Junta de Andaluc´ıa (research group “Communications Engineering (TIC-0102)”).

#146253 - $15.00 USD

(C) 2011 OSA

