Abstract

Recently, a new and generalized statistical model, called ℳ or Málaga distribution, was proposed to model the irradiance fluctuations of an unbounded optical wavefront (plane and spherical waves) propagating through a turbulent medium under all irradiance fluctuation conditions in homogeneous, isotropic turbulence. Málaga distribution was demonstrated to have the advantage of unifying most of the proposed statistical models derived until now in the bibliography in a closed-form expression providing, in addition, an excellent agreement with published plane wave and spherical wave simulation data over a wide range of turbulence conditions (weak to strong). Now, such a model is completed by including the adverse effect of pointing error losses due to misalignment. In this respect, the well-known effects of aperture size, beam width and jitter variance are taken into account. Accordingly, after presenting the analytical expressions for the combined distribution of scintillation and pointing errors, we derive its centered moments of the overall probability distribution. Finally, we obtain the analytical expressions for the average bit error rate performance for the ℳ distribution affected by pointing errors. Numerical results show the impact of misalignment on link performance.

© 2012 OSA

1. Introduction

Atmospheric optical communication (AOC) is receiving considerable attention recently for use in high data rate wireless links. However, even in clear sky conditions, wireless optical links experience fluctuations in both the intensity and the phase of an optical wave propagating through this medium [1] due to time varying inhomogeneities in the refractive index of the atmosphere. The reliability of an optical system operating in such an environment can be deduced from a mathematical model for the probability density function (pdf) of the randomly fading irradiance signal (scintillation). For that reason, one of the goals in studying optical wave propagation through turbulence is the identification of a tractable pdf of the irradiance under all intensity fluctuation regimes.

Over the years, many irradiance pdf models have been proposed with different degrees of success. Perhaps the most successful models are the lognormal [2] and the gamma-gamma [3] ones. The scope of the lognormal model is restricted under weak irradiance fluctuations. Conversely, the gamma-gamma pdf was suggested by Andrews et al. as a reasonable alternative to Beckmann’s pdf [4] because of its much more tractable mathematical model. Nevertheless, both lognormal and gamma-gamma models are particular cases of the new and recently proposed ℳ (Málaga) distribution model [5]. ℳ distribution was validated by comparing its pdf with published simulation data, showing that it unifies in an analytical expression most of the irradiance statistical models proposed in literature by the scientific community for more than four decades. So, for instance: Rice-Nakagami [1], gamma [3], shadowed-Rician [6], K [3, 7], homodyned-K [7], exponential [3] or Gamma-Rician [5] pdfs are considered as particular cases of this Málaga distribution.

In addition to the scintillation effect, misalignment between transmitter and receiver due to building sway causes vibrations of the transmitted beam and, thus, pointing errors that limit the performance of atmospheric free-space optical links. The impact of pointing error (jitter) was widely investigated for intersatellite space-based free-space optical links [8, 9] although without providing closed-form expressions. In a further step, Kiasaleh [10] deduced approximated expressions for the pdf of optical signal intensity in an optical communication channel impaired by motion-induced beam jitter and turbulence, assuming that the beam scintillation is governed by a lognormal pdf and a K distribution. Recently, a new procedure to derive the statistical distribution of pointing errors was proposed in [11]. Results derived in this last work were employed by Sandalidis et al. [12] to study the average bit error rate (ABER) expression for an intensity-modulation/direct detection (IM/DD) atmospheric optical system with on-off keying (OOK), assuming pointing errors, but restricted to strong turbulence.

Now, through this paper we widely generalize the results obtained in [12] for all turbulence regimes, thus completing our works reported in [5, 13] with ℳ distributed turbulence, by including the impact of misalignment fading on the performance of AOC systems. In this sense, we present a statistical AOC channel model considering the joint effects of ℳ distributed turbulence and pointing errors. As detailed in [5], ℳ distribution is valid under any turbulence regime, unifying in a closed-form expression the link performance of most of the irradiance statistical models proposed in literature. In addition, we derive the centered moments of the overall probability distribution and, finally, we provide closed-form expressions for the ABER performance of an IM/DD system over the combined effect of Málaga turbulence and misalignment fading, showing the impairment in performance of pointing errors in atmospheric optical links.

2. System model

In this paper, we consider IM/DD channels using on-off keying (OOK) modulation, which is widely employed in practical systems. Laser beams propagate along a horizontal path through a turbulence channel governed by an ℳ distribution in the presence of pointing errors. An additive white Gaussian noise (AWGN) is assumed that includes any shot noise caused by ambient light that may be much stronger than the desired signal as well as front-end thermal noise in the electronics following the photodetector. The receiver integrates the photocurrent signal which is related to the incident optical power by the detector responsivity, R, for each bit period. As explained above, the received optical power, PR(t), suffers from a fluctuation in signal intensity due to atmospheric turbulence and misalignment, as well as additive noise, and can be modeled, from [11], as:

PR(t)=h(t)RPT(t)+N(t),
with PT (t) being the transmitted optical power, with R being the responsivity, whereas N(t) is the AWGN with variance σN2, and h(t) is the normalized channel fading coefficient representing the effect of the intensity fluctuations on the transmitted signal due to the combined actions of scintillation and misalignment, in the way explained in [11]. This channel state is considered to be a product of two random components with a third and deterministic factor. The random components that take part in h are geometric spread and pointing errors, denoted by hp and atmospheric turbulence, represented by ha, following the same notation employed in [11]. The deterministic component that must be included in h arises due to atmospheric path loss, hl. This atmospheric attenuation, hl, is described by the exponential Beers-Lambert Law and can be expressed in terms of the visibility [11]. However, although weather-induced attenuation can also degrade the performance of AOC systems in the way shown in [14], due to its deterministic nature hl acts as a scaling factor as shown in [15] and for this reason, it is not considered in this manuscript. Thus, through this paper, h(t) = ha(t) ·hp(t), and each of these two factors are explained in detail through Sections 3 and 4, respectively.

Finally, and as in [11], it is considered that the pointing error is due to building sway, with sub-Hz bandwidth (correlation times of a few seconds), so that the scintillation process must have a smaller correlation time (10–100 ms) as compared with building sways. Hence, the two above impairments are stated to be independent.

3. Generalized atmospheric distribution model

In this section, we detail the analytical development for the factor ha introduced in Section 2. The starting point is the new statistical model for the irradiance fluctuations of an unbounded optical wavefront propagating through a turbulent medium under all irradiance fluctuation conditions in homogeneous, isotropic turbulence that was presented in [5]. For this reason, a brief description of such a model is achieved through this section in order to improve the overall understanding of this paper. The proposed model, labeled Málaga or simply ℳ, is based on a physical model of the scattering processes that is seen as an extension of the previous work developed by Churnside and Clifford [4]. The main advantage of the Málaga distribution associated to the proposed physical model is that it does reduce to a simple closed-form analytical formula valid from every turbulence regime.

Thus, assume an electromagnetic wave propagating through a turbulent atmosphere with a random refractive index. As the wave passes through this medium, part of the energy is scattered and the form of the irradiance probability distribution is determined by the type of scattering involved. The basis of the ℳ distribution is the conception of an alternative physical model in the generation of small-scale fluctuations. The propagation scheme is illustrated in Fig. 1. As detailed in [5] the observed field at the receiver is supposed to consist of three terms: the first one is the line-of-sight (LOS) contribution, UL, the second one is the component which is quasi-forward scattered by the eddies on the propagation axis, USC, and coupled to the LOS contribution; whereas the third term, USG, is due to energy which is scattered to the receiver by off-axis eddies, this latter contribution being statistically independent from the previous two other terms. The inclusion of the coupled to the LOS scattering component, USC, is the main novelty of the model and it can be justified by the high directivity and the narrow beamwidths of laser beams in atmospheric optical communications.

 

Fig. 1 Proposed propagation geometry for a laser beam in a Málaga model to form the small-scale fluctuations [5].

Download Full Size | PPT Slide | PDF

Accordingly, the total observed field when atmospheric turbulence is the only adverse effect considered in h can be written as:

U=(UL+USC+USG)exp(χ+jS)
where UL=GΩexp(jϕA), USC=Gρ2b0exp(jϕB) and USG=(1ρ)US, being USC and USG statistically independent stationary random processes. Of course, UL and USG are also independent random processes. In Eq. (2), G is a real variable following a gamma distribution with E[G]=1. It represents the slow fluctuation of the LOS component. Following the notation of [6], the parameter Ω =E[|UL|2] represents the average power of the LOS term whereas the average power of the total scatter components is denoted by 2b0=E[|USC|2+|USG|2]. ϕA and ϕB are the deterministic phases of the LOS and the coupled-to-LOS scatter terms, respectively. On another note, 0≤ρ≤1 is the factor expressing the amount of scattering power coupled to the LOS component. Finally, US is a circular Gaussian complex random variable, and χ and S are real random variables representing the log-amplitude and phase perturbation of the field induced by the atmospheric turbulence, respectively. A plausible justification for the coupled-to-LOS scattering component, USC, is given in [16]. There, it is said that if the turbulent medium is so thin that multiple scattering can be ignored, the multipath delays of the scattered radiation collected by a diffraction-limited receiver will usually be small relative to the signal bandwidth. Then the scattered field will combine coherently with the unscattered field and there will be no-“interfering” signal component of the field, in a similar way as USC combines with UL in our proposed model. Of course, when the turbulent medium becomes so thick, then the unscattered component of the field can be neglected.

Now, and from Eq. (2), the observed irradiance of the proposed propagation model (ha, as we have adopted the notation given in [11]) can be expressed as:

ha=|UL+USC+USG|2exp(2χ)=YX,{Y|UL+USC+USG|2(small-scalefluctuations)Xexp(2χ)(large-scalefluctuations);
whereas the lowpass-equivalent complex envelope can be rewritten as:
R(t)=(UL+USC+USG)=G[Ωexp(jϕA)+ρ2b0exp(jϕB)]+(1ρ)US,
so that we have the identical shadowed Rice single model employed in [6], composed by the sum of a Rayleigh random phasor (the independent scatter component, US) and a Nakagami distribution, G. Then, we can apply the same procedure exposed in [6] consisting in calculating the expectation of the Rayleigh component with respect to the Nakagami distribution and then deriving the pdf of the instantaneous power. Hence, the pdf of Y is given by:
fY(y)=1γ[γβγβ+Ω]βexp[yγ]1F1(β;1;1γΩ(γβ+Ω)y)
where β ≜ (E[G])2/Var[G] is the amount of fading parameter with Var[·] being the variance operator, where | · | denotes the absolute value. For simplicity, we have denoted γ=E[|USG|2]=2b0(1ρ) whereas Ω=Ω+ρ2b0+22b0Ωρcos(ϕAϕB)represents the average power from the coherent contributions, as detailed in [5]. Finally, 1F1 (a; c; x) is the Kummer confluent hypergeometric function of the first kind.

Otherwise, the large-scale fluctuations, X ≜ exp (2χ), is widely accepted to be a lognormal amplitude [4] but, as in [3, 6], this distribution is approximated by a gamma one, this latter with a more favorable analytical structure. For further details about the fact that the gamma distribution closely approximates the lognormal distribution, we suggest to consult references [17, 18]. Then:

fX(x)=ααΓ(α)xα1exp(αx),
where α is a positive parameter related to the effective number of large-scale cells of the scattering process, as in [3].

Thus to obtain the statistical characterization of the intensity, ha, the mixture of fX (x) and fY (y) must be accomplished. After the corresponding and non-obvious algebraic manipulation, extensively detailed in [5], the pdf of the irradiance is represented by:

fha(ha)=Ak=1βakhaα+k21Kαk(2αβhaγβ+Ω)
where
{A2αα2γ1+α2Γ(α)(γβγβ+Ω)β+α2;ak(β1k1)(γβ+Ω)1k2(k1)!(Ωγ)k1(αβ)k2.
with β being a natural number. In Eq. (7), Kν (·) is the modified Bessel function of the second kind and order ν.

A generalized expression of Eq. (7) was also given in [5] with β being a real number:

fha(ha)=A(G)k=1ak(G)haα+k21Kαk(2αhaγ),
where
{A(G)2αα2γ1+α2Γ(α)(γβγβ+Ω)β;ak(G)(β)k1(αγ)k2[(k1)!]2γk1(Ω+γβ)k1.

In Eq. (10), (β)k represents the Pochhammer symbol.

4. Misalignment fading model

As explained in Section 2, the channel fading coefficient, h(t) is considered through this paper to be a product of two random components: the scintillation, ha, detailed in Section 3, and geometric spread and pointing errors, denoted by hp. This last factor is now analytically detailed through this section. Thus, in atmospheric optical links, pointing accuracy is an important issue in determining link performance. For this reason, Farid and Hranilovic [11] derived its statistical nature, fhp (hp), assuming a Gaussian spatial intensity profile of beam waist, wz, at distance z from the transmitter:

Ibeam(ρ;z)=2πwz2exp(2ρ2wz2),
where ρis the radial vector from the beam center. In addition, a circular detection aperture of radius a was considered and a Gaussian beam profile at the receiver, Ibeam. Then, the fraction of the collected power at a receiver of radius a in the transverse plane of the incident wave when a pointing error, r, is present can be approximated as the Gaussian form:
hp(r;z)A0exp(2r2wzeq2),
where r = ||r|| is the radial distance. In particular, due to the symmetry of the beam shape and the detector area, hp depends only on this magnitude, r. On another note, A0 = [erf(v)]2 is the fraction of the collected power at r = 0, where v=(πa)/(2wz), with erf(·) denoting the error function, whereas the square of the equivalent beam width is given by:
wzeq2=wz2πerf(v)2vexp(v2).

By considering independent identical Gaussian distributions for the elevation and the horizontal displacement (sway), the radial displacement, r, at the receiver follows a Rayleigh distribution.

fr(r)=rσs2exp(r22σs2),r>0;
where σs2 is the jitter variance at the receiver. Then, from Eqs. (12) and (14), the pdf of hp is given by:
fhp(hp)=g2A0g2hpg21,0hpA0;
where g = wzeq/(2σs) is the ratio between the equivalent beam radius at the receiver, wzeq, and the pointing error displacement standard deviation at the receiver, σs. See [11] for more details.

5. Combined channel statistical model

In this section, we derive a complete statistical model of a stochastic atmospheric optical channel taking into account both turbulence-induced scintillation and misalignment-induced fading. Therefore, the unconditional pdf, fh(h), for the channel state, h, is obtained by calculating the mixture of the two distributions presented above in Eqs. (7) or (9) and (15):

fh(h)=fh|ha(h|ha)fha(ha)dha,
where fh|ha (h|ha) is the conditional probability given a turbulence state, ha, and it is expressed as:
fh|ha(h|ha)=1hafhp(hha)=g2A0g2ha(hha)g21,0hA0ha;
and where fha (ha) was written in Eqs. (7) and (9) depending on whether parameter β is a natural number or a real one, respectively. For the first case, we substitute Eq. (7) into the integral expression given by Eq. (16) to finally obtain fh(h):
fh(h)=g2AA0g2hg21k=1βakh/A0haα+k21g2Kαk(2αβhaγβ+Ω)dha
where the parameters A and ak were already defined in Eq. (8). From [19, Eq. (07.34.03.0605.01)], the modified Bessel function of the second kind, Kv(·) can be expressed as a special case of the Meijer G function, given by the following relationship:
G0,22,0(x|a,b)=2x(a+b)/2Kab(2x),
where a definition of the Meijer G function can be found at [19, Eq. (07.34.02.0001.01)]. Taking a = −b = (αk)/2 in Eq. (19) and substituting such an equation into Eq. (18), then:
fh(h)=g2A2A0g2hg21k=1βakh|A0haα+k21g2G0,22,0(αβhaγβ+Ω|αk2,αk2)dha.

Finally, to solve the last integrate, we can employ Eq. (07.34.21.0085.01) in [19] and, thus, the unconditional pdf, fh(h), for the channel state, h, can be written in a closed-form expression:

fh(h)=g2A2h1k=1βak(αβγβ+Ω)α+k2G1,33,0(αβγβ+ΩhA0|g2+1g2,α,k),
just after having simplified the G-Meijer function using the relationship
xrGp,qm,n(x|apcq)=Gp,qm,n(x|ap+rcq+r)
gathered in [19, Eq. (07.34.17.0011.01)].

For the generic case of Eq. (9) with β being a real number, the obtention of fh(h) is not so evident because it is firstly mandatory to verify if the infinite summation of Eq. (9) and the integrate operator that appears when building fh(h) can be interchanged. To prove that, by normalization, it is possible to express Eq. (9) as a discrete mixture [20] in the form of f(ha)=k=1ωkfk(ha) where ωk0k,k=1ωk=1 and with fk(·) being a Gamma-Gamma pdf. Since the cumulative distribution function (CDF) can be written as F(ha)=k=1ωkFk(ha), where Fk(ha) is the CDF of each fk(·), then we can build a new CDF in the form: ϕn(ha)=k=1nωkFk(ha). By applying the ratio criterion [21] it is easy to test that Fk(ha) converges absolutely so ϕn(ha) is said to converge weakly to F(ha). This condition and the fact that Eq. (17) is bounded and continuous allow us to apply the Helly-Bray theorem [22]. By proceeding in an analogous manner as in [13], we can conclude that all the conditions needed to apply the Helly-Bray theroem [22] are satisfied. Thus, operating in a similar way to obtain Eq. (21), then fh(h) is written as an infinite series of Meijer-G functions as

fh(h)=g2A(G)2h1k=1ak(G)(αγ)α+k2G1,33,0(αγhA0|g2+1g2,α,k).

5.1. Moments of the combined ℳ probability distribution

In this subsection, the kth moment of the channel state, h, following a combined ℳ probability distribution, fh(h), as given in Eq. (23) is derived. As we indicated in Section 2, the channel state is considered to be a product of hp and ha, where hp represents the random effect of pointing errors whereas ha denotes the turbulence-induced scintillation suffered by the observed irradiance. Furthermore, ha can be decomposed into two other components, as indicated in Eq. (3): ha = YX. The first component, denoted by Y and following a shadowed-Rician distribution [6], is associated to small-scale fluctuations of the irradiance, i.e., the small-scale contributions to scintillation associated with turbulent cells smaller than either the first Fresnel zone or the transverse spatial coherence radius, whichever is smallest; whereas the second component of ha, denoted by X and following a gamma distribution, is related to to large-scale fluctuations of the irradiance that are generated by turbulent cells larger than that of either the Fresnel zone or the so-called “scattering disk”, whichever is largest.

The kth centered moments of ha, denoted by mk(ha), following an ℳ distribution was derived in [5]. Based on assumptions of statistical independence for the underlying random processes, X and Y, then:

mk(ha)E[hak]=E[Xk]E[Yk]=mk(X)mk(Y),
whereas the operator E[·] stands for ensemble average. From [23, Eq. (2.23)], the moment of a Nakagami-m pdf is given by:
mk(X)=Γ(α+k)Γ(α)αk,

In addition, the moment of a shadowed-Rician distribution with the restriction of β being a natural number was derived in [5]. Then:

mk(Y)=1γ(γβγβ+Ω)βr=0β1(β1r)1r!(Ωγ(γβ+Ω))rΓ(k+r+1)(βγβ+Ω)k+r+1.

Hence, when performing the product of Eq. (25) by Eq. (26), we certainly obtain:

mk(ha)=Γ(α+k)Γ(α)αkγ(γβγβ+Ω)βr=0β1(β1r)1r!(Ωγ(γβ+Ω))rΓ(k+r+1)(βγβ+Ω)k+r+1,
as indicated in [5].

And for the generic case of Eq. (9) with the parameter β being a real number, the moment of a shadowed-Rician distribution was given in [23, Eq. (2.69)]:

mk(Y)=(γβγβ+Ω)βγkΓ(k+1)2F1(k+1,β;1;Ωγβ+Ω),
and, thus, Eq. (27) is now generalized as:
mk(G)(ha)=Γ(α+k)Γ(α)αk(γβγβ+Ω)βγkΓ(k+1)2F1(k+1,β;1;Ωγβ+Ω).

On the other hand, as discussed in [11] and indicated in Section 2, it is considered that pointing errors and turbulence-induced scintillation are independent processes. As in the case of ha, we can obtain the kth centered moments of the channel state, h, as

mk(h)E[hk]=E[hak]E[hpk]=mk(ha)mk(hp).

To calculate mk (hp) we must solve:

mk(hp)=E[hpk]=hpkfhp(hp)dhp.

Inserting Eq. (15) into Eq. (31), then

mk(hp)=g2A0g20A0hpg2+k1dhp=g2g2+kA0k.

Hence, for the case of β being a natural number, we can obtain the kth centered moments of h by performing the product of Eqs. (27) and (32):

mk(h)=Γ(α+k)Γ(α)αkγ(γβγβ+Ω)βg2g2+kA0k××r=0β1(β1r)1r!(Ωγ(γβ+Ω))rΓ(k+r+1)(βγβ+Ω)k+r+1.

Finally, for the generic case of β being a real number, the kth centered moments of h is obtained by the product of Eqs. (29) and (32):

mk(h)=g2A0kg2+kΓ(α+k)Γ(α)αk(γβγβ+Ω)βγkΓ(k+1)2F1(k+1,β;1;Ωγβ+Ω).

6. Average BER

The study of the ABER of the ℳ probability distribution in the absence of pointing errors was derived in [13]. Now, in this section, we complete such a previous work by detailing the procedure to derive the analytical expressions for the ABER of the ℳ probability distribution in the presence of misalignment fading. In this sense, an IM/DD link using OOK signalling technique is considered again. For such links, each bit symbol is transmitted by pulsing the light source either on or off during each bit time. Since in most practical systems, the receiver signal-to-noise ratio (SNR) is limited by shot noise caused by ambient light that may be much stronger than the desired signal and/or by thermal noise in the electronics following the photodetector, then noise can usually be modeled to high accuracy as AWGN that is statistically independent of the desired signal [2]. Because of AWGN, errors may be made at the receiver in determining the actual symbols transmitted. The overall probability of error, Pb(e) is expressed as the following weighted sum:

Pb(e)=p0p(e|0)+p1p(e|1),
where p0 is the transmission probability of a binary “0” and p1 is the transmission probability of a binary “1”, as indicated in [24]. On another note, p(e|0) and p(e|1) denote the conditional bit error probabilities when the transmitted bit is “0” and “1”, respectively. We assume that either symbol (binary “0” or binary “1”) are equally likely to be transmitted, i.e., p0 = p1 = 1/2. Considering that p(e|0) = p(e|1), then it is straightforward to show that the conditional BER is expressed, from [13, 24], as
Pb(e|h)=p(e|0,h)=p(e|1,h)=Q(PRhσN)
where P is the average of optical power, with R being the responsivity and σN2 being the AWGN variance while Q(·) is the Gaussian Q function which is related to the complementary error function, erfc(·), by erfc(X)=2Q(2x). The ABER, Pb(e), can be obtained by averaging Pb(e|h) over fh(h):
Pb(e)=0Pb(e|h)fh(h)dh,
as shown in [12, 13].

Thus, for the particular case of β being a natural number, the pdf of the channel state, fh(h), is taken from Eq. (21). Now we can express the Gaussian Q function appeared in Eq. (36) in terms of the complementary error function, erfc(x). Next, we can write both the erfc(·) and the modified Bessel integrands as Meijer G functions by employing [19, Eqs. (07.34.03.0619.01), (07.34.03.0605.01) ], respectively, as in [12, 13]. Such relationships are written here for added convenience: in fact, we have already written in Eq. (19) the relationship between a modified Bessel function and a Meijer G function. In the same way, the erfc(·) and the Meijer G function are related by:

G1,22,0(z|aa1,a12)=πza1erfc(z).

Finally, using [19, Eq. (07.34.21.0013.01)], the ABER is derived in a closed-form expression as:

Pb(e)=2αg2ABα232ππk=1β2kBk2akG7,42,6(8R2P2A02σN2B2|1g22,2g22,1α2,2α2,1k2,2k2,10,12,g22,1g22)
with B=γβ+Ωαβ.

For the generic case of the irradiance pdf represented by Eq. (9) with β being a real number, we verify again that the infinite summation and the integrate operator can be interchanged. Operating in a similar way as to obtain Eq. (39), the ABER is written as:

Pb(e)=2αg2A(G)32ππ(γα)α2k=12k(γα)k2ak(G)××G7,42,6(8R2P2A02σN2γ2α2|1g22,2g22,1α2,2α2,1k2,2k2,10,12,g22,1g22).

7. Numerical results and discussions

Some numerical results are shown in Figs. 23, where the IM/DD AOC system is evaluated for its error performance capabilities against average SNR, and in terms of the Rytov variance, σR2, as in [12, 13]. Effectively, as indicated in [24], atmospheric turbulence is well characterized by such a Rytov variance when using the optical wave models of an infinite plane wave or a spherical wave (point source). The Rytov variance is defined, for a plane wave, by [24, 11]:

σR2=1.23Cn2k7/6L11/6.
where Cn2(m2/3) is the index of refraction structure parameter of the atmosphere, k = 2π/λ is the optical wave number, λ (m) is the wavelength, and L (m) is the propagation path length between transmitter and receiver. To obtain typical values of σR2, we have employed some experimental data of Cn2 provided by University of Waseda, in Japan, on October 15th, 2009, obtained from atmospheric parameters. Thus, near midday, Cn2 parameter oscillated from 2.8 to 11 · 10−14 m−2/3, and the temperature oscillated from 22 °C to 26.4 °C, and the received power was recorded from −22.52 dBm to −20.48 dBm. At night (near 1 a.m.), Cn2 parameter registered its minimum value (∼ 7 ·10−15 m−2/3), for a temperature of 15.1 °C and a received power of −22.91 dBm. At sunrise (6.45 a.m.), Cn2=1.21014m2/3, for a temperature of 15.9 °C and with a received power of −23.47 dBm. The wavelength employed in all cases was 785 nm and the receiver aperture had a diameter of 100 mm. The propagation path length is L = 1 km at a height of 25 m. The optical power transmitted is 11.5 dBm with a responsivity of 0.8 A/W. With such a set of parameters, if we consider Cn2 to 7.2 ·10−15 m−2/3 (night), then σR20.32; secondly, for the case of sunrise when Cn2=1.2m2/3, then σR20.52; finally, for a Cn2=2.81014m2/3 (near midday), then σR21.2. This will be the values of Rytov variance employed through this section. As a final remark, all cases presented were obtained by employing the same normalized average optical power (Ω + 2b0 = 1).

 

Fig. 2 Average optical BER against average SNR for different values of ρ and wz/a. In all curves, α = 10, β = 5 and the transmitted power is normalized, i.e., Ω + 2b0 = 1. The case of ρ = 1 corresponds to the Gamma-Gamma distribution.

Download Full Size | PPT Slide | PDF

 

Fig. 3 Average optical BER against average SNR for different values of α, β, ρ assuming σR2=0.36 and wz/a = 20. In all curves, the transmitted power is normalized, i.e., Ω + 2b0 = 1.

Download Full Size | PPT Slide | PDF

Figure 2 shows the ABER for various values of the normalized beamwidth: wz/a=10, 20, and 25. Parameters of ℳ turbulence model were fixed to: α = 10 and β = 5 in all cases. Moreover, for each wz/a value, three different behaviors were obtained: ρ= 1, which corresponds to a Gamma-Gamma distribution; ρ = 0.75, and ρ = 0.25. For such values, the associated Rytov variance is σR2=0.32,0.52 and 1.2, respectively, coinciding with the values selected above. As a reference, the behavior of the system in the absence of misalignment fading is also provided (dotted curves) for the three intensities of turbulence mentioned before [13]. Hence, it is straightforward to check the impairment in performance when pointing errors are considered: for instance, the penalties in SNR with respect to the ideal case where no pointing errors are considered are 17.03, 23.02 and 24.95 optical dB, for wz/a = 10, 20 and 25, respectively, for a bit error rate of 10−6.

In Fig. 3, different ABERs are displayed for wz/a = 20 and a same intensity of turbulence ( σR2=0.36, Cn2=0.831014m2/3, measured on October 15th at 23:10), a typical value in terrestrial horizontal links. It turns out that the behavior of each curve is different depending on ρ. Thus, when ρ = 1 the overall scattering power travels through USC, so it is entirely coupled to the LOS component, minimizing its adverse effect, as displayed in Fig. 3. When ρ is getting lower then the scattering power is distributed both into USC and into USG. The presence of a higher power from the USG component has a more harmful effect in the performance of an AOC system because it is not correlated to the LOS contribution, as corroborated in Fig. 3. Obviously, obtained results are consistent with the ones achieved in [13], being the results shown in this paper a 23.03 optical dB displaced version with reference to the ones displayed in Fig. 3 of [13]. For this reason, same conclusions derived there are applicable when pointing errors are considered.

On the other hand, all BER curves tend to a same diversity order, where the diversity order is defined in [25].

Finally, we must recall that Eq. (7) is a particularization of Eq. (9) valid when β is restricted to a natural number. Nevertheless, due to the high degree of freedom inherent to the ℳ distribution, demonstrated in the set of results shown in [5], this particularization allows us to reproduce every turbulent scenario, avoiding the necessity of employing the infinite summation included in Eqs. (9) or (40).

8. Concluding remarks

In summary, in this paper, two analytical closed-form representations for the ABER performance of an AOC system operating over a generalized turbulence in the presence of pointing errors are derived. Due to the high degree of freedom inherent to the ℳ distribution [5], the particularization given in Eq. (39) let us reproduce every turbulent scenario, avoiding the infinite summation included in Eq. (40).

As discussed in Section 7, the effect of pointing errors on the final performance of an AOC system can be seen as a penalty in SNR with respect to the ideal case where no pointing errors are considered, but with the special feature that such a penalty in SNR is seen as a constant value for every different turbulence regime when the ratio wz/a is fixed, and for a given value of bit error rate.

Finally, as can be deduced, the consideration of both factors, turbulence and misalignment, affecting the AOC systems performance, is an essential AOC link design criterion. Through this paper, we have derived a generalized analytical expression to represent the ABER of an AOC system including the presence of pointing errors, so that a designer can obtain in a straightforward manner the ABER related to a particular location affected for both turbulence-induced scintillation and pointing errors due to misalignment between the transmitter and receiver.

Acknowledgment

This work was supported by the Spanish Ministerio de Ciencia e Innovación, Project TEC2008-06598. The authors thank Prof. Mitsuji Matsumoto from University of Waseda, in Japan, for his kind support and for providing some experimental measurements of the Cn2 parameter needed to calculate realistic values of Rytov variance.

References and links

1. L. C. Andrews and R. L. Phillips, Laser Beam Propagation Through Random Media (SPIE, 1998).

2. X. Zhu and J. M. Kahn, “Free-space optical communication through atmospheric turbulence channels,” IEEE Trans. Commun. 50, 1293–1300 (2002). [CrossRef]  

3. 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). [CrossRef]  

4. J. H. Churnside and S. F. Clifford, “Log-normal Rician probability-density function of optical scintillations in the turbulent atmosphere,” J. Opt. Soc. Am. A 4, 1923–1930 (1987). [CrossRef]  

5. A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, and A. Puerta-Notario, “A unifying statistical model for atmospheric optical scintillation,” in: Numerical Simulations of Physical and Engineering Processes (Intech, 2011). [CrossRef]  

6. A. Abdi, W. C. Lau, M. S. Alouini, and M. A. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” IEEE Trans. Wireless Commun. 2, 519–528 (2003). [CrossRef]  

7. E. Jakerman, “On the statistics of K-distributed noise,” J. Phys. A 13, 31–48 (1980). [CrossRef]  

8. S. Arnon and N. S. Kopeika, “Laser satellite communication network-vibration effect and possible solutions,” Proc. IEEE 85, 1646–1661 (1997). [CrossRef]  

9. S. Arnon, “Optimization of urban optical wireless communication systems,” IEEE Trans. Wireless Commun. 2, 626–629 (2003). [CrossRef]  

10. K. Kiasaleh, “On the probability density function of signal intensity in free-space optical communications systems impaired by pointing jitter and turbulence,” Opt. Eng. 33, 3748–3757 (1994). [CrossRef]  

11. A. A. Farid and S. Hranilovic, “Outage capacity optimization for free-space optical links with pointing errors,” J. Lightwave Technol. 25, 1702–1710 (2007). [CrossRef]  

12. 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, 44–46 (2008). [CrossRef]  

13. A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, M. Castillo-Vázquez, and A. Puerta-Notario, “General analytical expressions for the bit error rate of atmospheric optical communication systems,” Opt. Lett. 36, 4095–4097 (2011). [CrossRef]   [PubMed]  

14. M. Al Naboulsi and H. Sizun, “Fog attenuation prediction for optical and infrared waves,” Opt. Eng. 43, 319–329 (2004). [CrossRef]  

15. ITU-R Report F.2106-1 “Fixed service applications using free-space optical links,” Nov. 2010.

16. R. S. Kennedy, “Communication through optical scattering channels: an introduction,” Proc. IEEE 58, 1651–1665 (1970). [CrossRef]  

17. T. Rappaport, Wireless Communications: Principles and Practice, 2nd ed. (Prentice Hall, 2001.).

18. J. R. Clark and S. Karp, “Approximations for lognormally fading optical signals,” Proc. IEEE 58, 1964–1965 (1970). [CrossRef]  

19. Wolfram, http://functions.wolfram.com/

20. C. A. Charalambides, Combinatorial Methods in Discrete Distributions (John Wiley & Sons, 2005) [CrossRef]  

21. L. C. Andrews, Special Functions of Mathematics for Engineers, 2nd ed. (SPIE, 1998).

22. P. Billingsley, Convergence of Probability Measures, 2nd ed. (John Wiley & Sons, 2005)

23. M. K. Simon and M. S. Alouini, Digital Communications over Fading Channels, 2nd ed. (Wiley-Interscience, 2005).

24. L. C. Andrews, R. L. Phillips, and C. Y. Hopen, Laser Beam Scintillation with Applications (SPIE, 2001). [CrossRef]  

25. L. Zheng and D. N. C. Tse, “Diversity and multiplexing: A fundamental tradeoff in multiple antenna channels,” IEEE Trans. Inf. Theory 49, 1073–1096 (2003). [CrossRef]  

References

  • View by:
  • |
  • |
  • |

  1. L. C. Andrews and R. L. Phillips, Laser Beam Propagation Through Random Media (SPIE, 1998).
  2. X. Zhu and J. M. Kahn, “Free-space optical communication through atmospheric turbulence channels,” IEEE Trans. Commun. 50, 1293–1300 (2002).
    [CrossRef]
  3. 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).
    [CrossRef]
  4. J. H. Churnside and S. F. Clifford, “Log-normal Rician probability-density function of optical scintillations in the turbulent atmosphere,” J. Opt. Soc. Am. A 4, 1923–1930 (1987).
    [CrossRef]
  5. A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, and A. Puerta-Notario, “A unifying statistical model for atmospheric optical scintillation,” in: Numerical Simulations of Physical and Engineering Processes (Intech, 2011).
    [CrossRef]
  6. A. Abdi, W. C. Lau, M. S. Alouini, and M. A. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” IEEE Trans. Wireless Commun. 2, 519–528 (2003).
    [CrossRef]
  7. E. Jakerman, “On the statistics of K-distributed noise,” J. Phys. A 13, 31–48 (1980).
    [CrossRef]
  8. S. Arnon and N. S. Kopeika, “Laser satellite communication network-vibration effect and possible solutions,” Proc. IEEE 85, 1646–1661 (1997).
    [CrossRef]
  9. S. Arnon, “Optimization of urban optical wireless communication systems,” IEEE Trans. Wireless Commun. 2, 626–629 (2003).
    [CrossRef]
  10. K. Kiasaleh, “On the probability density function of signal intensity in free-space optical communications systems impaired by pointing jitter and turbulence,” Opt. Eng. 33, 3748–3757 (1994).
    [CrossRef]
  11. A. A. Farid and S. Hranilovic, “Outage capacity optimization for free-space optical links with pointing errors,” J. Lightwave Technol. 25, 1702–1710 (2007).
    [CrossRef]
  12. 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, 44–46 (2008).
    [CrossRef]
  13. A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, M. Castillo-Vázquez, and A. Puerta-Notario, “General analytical expressions for the bit error rate of atmospheric optical communication systems,” Opt. Lett. 36, 4095–4097 (2011).
    [CrossRef] [PubMed]
  14. M. Al Naboulsi and H. Sizun, “Fog attenuation prediction for optical and infrared waves,” Opt. Eng. 43, 319–329 (2004).
    [CrossRef]
  15. ITU-R Report F.2106-1 “Fixed service applications using free-space optical links,” Nov.2010.
  16. R. S. Kennedy, “Communication through optical scattering channels: an introduction,” Proc. IEEE 58, 1651–1665 (1970).
    [CrossRef]
  17. T. Rappaport, Wireless Communications: Principles and Practice, 2nd ed. (Prentice Hall, 2001.).
  18. J. R. Clark and S. Karp, “Approximations for lognormally fading optical signals,” Proc. IEEE 58, 1964–1965 (1970).
    [CrossRef]
  19. Wolfram, http://functions.wolfram.com/
  20. C. A. Charalambides, Combinatorial Methods in Discrete Distributions (John Wiley & Sons, 2005)
    [CrossRef]
  21. L. C. Andrews, Special Functions of Mathematics for Engineers, 2nd ed. (SPIE, 1998).
  22. P. Billingsley, Convergence of Probability Measures, 2nd ed. (John Wiley & Sons, 2005)
  23. M. K. Simon and M. S. Alouini, Digital Communications over Fading Channels, 2nd ed. (Wiley-Interscience, 2005).
  24. L. C. Andrews, R. L. Phillips, and C. Y. Hopen, Laser Beam Scintillation with Applications (SPIE, 2001).
    [CrossRef]
  25. L. Zheng and D. N. C. Tse, “Diversity and multiplexing: A fundamental tradeoff in multiple antenna channels,” IEEE Trans. Inf. Theory 49, 1073–1096 (2003).
    [CrossRef]

2011 (1)

2008 (1)

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, 44–46 (2008).
[CrossRef]

2007 (1)

2004 (1)

M. Al Naboulsi and H. Sizun, “Fog attenuation prediction for optical and infrared waves,” Opt. Eng. 43, 319–329 (2004).
[CrossRef]

2003 (3)

L. Zheng and D. N. C. Tse, “Diversity and multiplexing: A fundamental tradeoff in multiple antenna channels,” IEEE Trans. Inf. Theory 49, 1073–1096 (2003).
[CrossRef]

S. Arnon, “Optimization of urban optical wireless communication systems,” IEEE Trans. Wireless Commun. 2, 626–629 (2003).
[CrossRef]

A. Abdi, W. C. Lau, M. S. Alouini, and M. A. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” IEEE Trans. Wireless Commun. 2, 519–528 (2003).
[CrossRef]

2002 (1)

X. Zhu and J. M. Kahn, “Free-space optical communication through atmospheric turbulence channels,” IEEE Trans. Commun. 50, 1293–1300 (2002).
[CrossRef]

2001 (1)

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).
[CrossRef]

1997 (1)

S. Arnon and N. S. Kopeika, “Laser satellite communication network-vibration effect and possible solutions,” Proc. IEEE 85, 1646–1661 (1997).
[CrossRef]

1994 (1)

K. Kiasaleh, “On the probability density function of signal intensity in free-space optical communications systems impaired by pointing jitter and turbulence,” Opt. Eng. 33, 3748–3757 (1994).
[CrossRef]

1987 (1)

1980 (1)

E. Jakerman, “On the statistics of K-distributed noise,” J. Phys. A 13, 31–48 (1980).
[CrossRef]

1970 (2)

R. S. Kennedy, “Communication through optical scattering channels: an introduction,” Proc. IEEE 58, 1651–1665 (1970).
[CrossRef]

J. R. Clark and S. Karp, “Approximations for lognormally fading optical signals,” Proc. IEEE 58, 1964–1965 (1970).
[CrossRef]

Abdi, A.

A. Abdi, W. C. Lau, M. S. Alouini, and M. A. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” IEEE Trans. Wireless Commun. 2, 519–528 (2003).
[CrossRef]

Al Naboulsi, M.

M. Al Naboulsi and H. Sizun, “Fog attenuation prediction for optical and infrared waves,” Opt. Eng. 43, 319–329 (2004).
[CrossRef]

Al-Habash, M. A.

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).
[CrossRef]

Alouini, M. S.

A. Abdi, W. C. Lau, M. S. Alouini, and M. A. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” IEEE Trans. Wireless Commun. 2, 519–528 (2003).
[CrossRef]

M. K. Simon and M. S. Alouini, Digital Communications over Fading Channels, 2nd ed. (Wiley-Interscience, 2005).

Andrews, L. C.

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).
[CrossRef]

L. C. Andrews and R. L. Phillips, Laser Beam Propagation Through Random Media (SPIE, 1998).

L. C. Andrews, R. L. Phillips, and C. Y. Hopen, Laser Beam Scintillation with Applications (SPIE, 2001).
[CrossRef]

L. C. Andrews, Special Functions of Mathematics for Engineers, 2nd ed. (SPIE, 1998).

Arnon, S.

S. Arnon, “Optimization of urban optical wireless communication systems,” IEEE Trans. Wireless Commun. 2, 626–629 (2003).
[CrossRef]

S. Arnon and N. S. Kopeika, “Laser satellite communication network-vibration effect and possible solutions,” Proc. IEEE 85, 1646–1661 (1997).
[CrossRef]

Billingsley, P.

P. Billingsley, Convergence of Probability Measures, 2nd ed. (John Wiley & Sons, 2005)

Castillo-Vázquez, M.

Charalambides, C. A.

C. A. Charalambides, Combinatorial Methods in Discrete Distributions (John Wiley & Sons, 2005)
[CrossRef]

Churnside, J. H.

Clark, J. R.

J. R. Clark and S. Karp, “Approximations for lognormally fading optical signals,” Proc. IEEE 58, 1964–1965 (1970).
[CrossRef]

Clifford, S. F.

Farid, A. A.

Garrido-Balsells, J. M.

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, M. Castillo-Vázquez, and A. Puerta-Notario, “General analytical expressions for the bit error rate of atmospheric optical communication systems,” Opt. Lett. 36, 4095–4097 (2011).
[CrossRef] [PubMed]

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, and A. Puerta-Notario, “A unifying statistical model for atmospheric optical scintillation,” in: Numerical Simulations of Physical and Engineering Processes (Intech, 2011).
[CrossRef]

Hopen, C. Y.

L. C. Andrews, R. L. Phillips, and C. Y. Hopen, Laser Beam Scintillation with Applications (SPIE, 2001).
[CrossRef]

Hranilovic, S.

Jakerman, E.

E. Jakerman, “On the statistics of K-distributed noise,” J. Phys. A 13, 31–48 (1980).
[CrossRef]

Jurado-Navas, A.

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, M. Castillo-Vázquez, and A. Puerta-Notario, “General analytical expressions for the bit error rate of atmospheric optical communication systems,” Opt. Lett. 36, 4095–4097 (2011).
[CrossRef] [PubMed]

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, and A. Puerta-Notario, “A unifying statistical model for atmospheric optical scintillation,” in: Numerical Simulations of Physical and Engineering Processes (Intech, 2011).
[CrossRef]

Kahn, J. M.

X. Zhu and J. M. Kahn, “Free-space optical communication through atmospheric turbulence channels,” IEEE Trans. Commun. 50, 1293–1300 (2002).
[CrossRef]

Karagiannidis, G. K.

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, 44–46 (2008).
[CrossRef]

Karp, S.

J. R. Clark and S. Karp, “Approximations for lognormally fading optical signals,” Proc. IEEE 58, 1964–1965 (1970).
[CrossRef]

Kaveh, M. A.

A. Abdi, W. C. Lau, M. S. Alouini, and M. A. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” IEEE Trans. Wireless Commun. 2, 519–528 (2003).
[CrossRef]

Kennedy, R. S.

R. S. Kennedy, “Communication through optical scattering channels: an introduction,” Proc. IEEE 58, 1651–1665 (1970).
[CrossRef]

Kiasaleh, K.

K. Kiasaleh, “On the probability density function of signal intensity in free-space optical communications systems impaired by pointing jitter and turbulence,” Opt. Eng. 33, 3748–3757 (1994).
[CrossRef]

Kopeika, N. S.

S. Arnon and N. S. Kopeika, “Laser satellite communication network-vibration effect and possible solutions,” Proc. IEEE 85, 1646–1661 (1997).
[CrossRef]

Lau, W. C.

A. Abdi, W. C. Lau, M. S. Alouini, and M. A. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” IEEE Trans. Wireless Commun. 2, 519–528 (2003).
[CrossRef]

Paris, J. F.

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, M. Castillo-Vázquez, and A. Puerta-Notario, “General analytical expressions for the bit error rate of atmospheric optical communication systems,” Opt. Lett. 36, 4095–4097 (2011).
[CrossRef] [PubMed]

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, and A. Puerta-Notario, “A unifying statistical model for atmospheric optical scintillation,” in: Numerical Simulations of Physical and Engineering Processes (Intech, 2011).
[CrossRef]

Phillips, R. L.

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).
[CrossRef]

L. C. Andrews and R. L. Phillips, Laser Beam Propagation Through Random Media (SPIE, 1998).

L. C. Andrews, R. L. Phillips, and C. Y. Hopen, Laser Beam Scintillation with Applications (SPIE, 2001).
[CrossRef]

Puerta-Notario, A.

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, M. Castillo-Vázquez, and A. Puerta-Notario, “General analytical expressions for the bit error rate of atmospheric optical communication systems,” Opt. Lett. 36, 4095–4097 (2011).
[CrossRef] [PubMed]

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, and A. Puerta-Notario, “A unifying statistical model for atmospheric optical scintillation,” in: Numerical Simulations of Physical and Engineering Processes (Intech, 2011).
[CrossRef]

Rappaport, T.

T. Rappaport, Wireless Communications: Principles and Practice, 2nd ed. (Prentice Hall, 2001.).

Sandalidis, H. G.

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, 44–46 (2008).
[CrossRef]

Simon, M. K.

M. K. Simon and M. S. Alouini, Digital Communications over Fading Channels, 2nd ed. (Wiley-Interscience, 2005).

Sizun, H.

M. Al Naboulsi and H. Sizun, “Fog attenuation prediction for optical and infrared waves,” Opt. Eng. 43, 319–329 (2004).
[CrossRef]

Tse, D. N. C.

L. Zheng and D. N. C. Tse, “Diversity and multiplexing: A fundamental tradeoff in multiple antenna channels,” IEEE Trans. Inf. Theory 49, 1073–1096 (2003).
[CrossRef]

Tsiftsis, T. A.

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, 44–46 (2008).
[CrossRef]

Uysal, M.

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, 44–46 (2008).
[CrossRef]

Zheng, L.

L. Zheng and D. N. C. Tse, “Diversity and multiplexing: A fundamental tradeoff in multiple antenna channels,” IEEE Trans. Inf. Theory 49, 1073–1096 (2003).
[CrossRef]

Zhu, X.

X. Zhu and J. M. Kahn, “Free-space optical communication through atmospheric turbulence channels,” IEEE Trans. Commun. 50, 1293–1300 (2002).
[CrossRef]

IEEE Commun. Lett. (1)

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, 44–46 (2008).
[CrossRef]

IEEE Trans. Commun. (1)

X. Zhu and J. M. Kahn, “Free-space optical communication through atmospheric turbulence channels,” IEEE Trans. Commun. 50, 1293–1300 (2002).
[CrossRef]

IEEE Trans. Inf. Theory (1)

L. Zheng and D. N. C. Tse, “Diversity and multiplexing: A fundamental tradeoff in multiple antenna channels,” IEEE Trans. Inf. Theory 49, 1073–1096 (2003).
[CrossRef]

IEEE Trans. Wireless Commun. (2)

A. Abdi, W. C. Lau, M. S. Alouini, and M. A. Kaveh, “A new simple model for land mobile satellite channels: first- and second-order statistics,” IEEE Trans. Wireless Commun. 2, 519–528 (2003).
[CrossRef]

S. Arnon, “Optimization of urban optical wireless communication systems,” IEEE Trans. Wireless Commun. 2, 626–629 (2003).
[CrossRef]

J. Lightwave Technol. (1)

J. Opt. Soc. Am. A (1)

J. Phys. A (1)

E. Jakerman, “On the statistics of K-distributed noise,” J. Phys. A 13, 31–48 (1980).
[CrossRef]

Opt. Eng. (3)

K. Kiasaleh, “On the probability density function of signal intensity in free-space optical communications systems impaired by pointing jitter and turbulence,” Opt. Eng. 33, 3748–3757 (1994).
[CrossRef]

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).
[CrossRef]

M. Al Naboulsi and H. Sizun, “Fog attenuation prediction for optical and infrared waves,” Opt. Eng. 43, 319–329 (2004).
[CrossRef]

Opt. Lett. (1)

Proc. IEEE (3)

R. S. Kennedy, “Communication through optical scattering channels: an introduction,” Proc. IEEE 58, 1651–1665 (1970).
[CrossRef]

S. Arnon and N. S. Kopeika, “Laser satellite communication network-vibration effect and possible solutions,” Proc. IEEE 85, 1646–1661 (1997).
[CrossRef]

J. R. Clark and S. Karp, “Approximations for lognormally fading optical signals,” Proc. IEEE 58, 1964–1965 (1970).
[CrossRef]

Other (10)

Wolfram, http://functions.wolfram.com/

C. A. Charalambides, Combinatorial Methods in Discrete Distributions (John Wiley & Sons, 2005)
[CrossRef]

L. C. Andrews, Special Functions of Mathematics for Engineers, 2nd ed. (SPIE, 1998).

P. Billingsley, Convergence of Probability Measures, 2nd ed. (John Wiley & Sons, 2005)

M. K. Simon and M. S. Alouini, Digital Communications over Fading Channels, 2nd ed. (Wiley-Interscience, 2005).

L. C. Andrews, R. L. Phillips, and C. Y. Hopen, Laser Beam Scintillation with Applications (SPIE, 2001).
[CrossRef]

A. Jurado-Navas, J. M. Garrido-Balsells, J. F. Paris, and A. Puerta-Notario, “A unifying statistical model for atmospheric optical scintillation,” in: Numerical Simulations of Physical and Engineering Processes (Intech, 2011).
[CrossRef]

T. Rappaport, Wireless Communications: Principles and Practice, 2nd ed. (Prentice Hall, 2001.).

ITU-R Report F.2106-1 “Fixed service applications using free-space optical links,” Nov.2010.

L. C. Andrews and R. L. Phillips, Laser Beam Propagation Through Random Media (SPIE, 1998).

Cited By

OSA participates in CrossRef's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (3)

Fig. 1
Fig. 1

Proposed propagation geometry for a laser beam in a Málaga model to form the small-scale fluctuations [5].

Fig. 2
Fig. 2

Average optical BER against average SNR for different values of ρ and wz/a. In all curves, α = 10, β = 5 and the transmitted power is normalized, i.e., Ω + 2b0 = 1. The case of ρ = 1 corresponds to the Gamma-Gamma distribution.

Fig. 3
Fig. 3

Average optical BER against average SNR for different values of α, β, ρ assuming σ R 2 = 0.36 and wz/a = 20. In all curves, the transmitted power is normalized, i.e., Ω + 2b0 = 1.

Equations (41)

Equations on this page are rendered with MathJax. Learn more.

P R ( t ) = h ( t ) R P T ( t ) + N ( t ) ,
U = ( U L + U S C + U S G ) exp ( χ + j S )
h a = | U L + U S C + U S G | 2 exp ( 2 χ ) = Y X , { Y | U L + U S C + U S G | 2 ( small - scale fluctuations ) X exp ( 2 χ ) ( large - scale fluctuations ) ;
R ( t ) = ( U L + U S C + U S G ) = G [ Ω exp ( j ϕ A ) + ρ 2 b 0 exp ( j ϕ B ) ] + ( 1 ρ ) U S ,
f Y ( y ) = 1 γ [ γ β γ β + Ω ] β exp [ y γ ] 1 F 1 ( β ; 1 ; 1 γ Ω ( γ β + Ω ) y )
f X ( x ) = α α Γ ( α ) x α 1 exp ( α x ) ,
f h a ( h a ) = A k = 1 β a k h a α + k 2 1 K α k ( 2 α β h a γ β + Ω )
{ A 2 α α 2 γ 1 + α 2 Γ ( α ) ( γ β γ β + Ω ) β + α 2 ; a k ( β 1 k 1 ) ( γ β + Ω ) 1 k 2 ( k 1 ) ! ( Ω γ ) k 1 ( α β ) k 2 .
f h a ( h a ) = A ( G ) k = 1 a k ( G ) h a α + k 2 1 K α k ( 2 α h a γ ) ,
{ A ( G ) 2 α α 2 γ 1 + α 2 Γ ( α ) ( γ β γ β + Ω ) β ; a k ( G ) ( β ) k 1 ( α γ ) k 2 [ ( k 1 ) ! ] 2 γ k 1 ( Ω + γ β ) k 1 .
I beam ( ρ ; z ) = 2 π w z 2 exp ( 2 ρ 2 w z 2 ) ,
h p ( r ; z ) A 0 exp ( 2 r 2 w z eq 2 ) ,
w z eq 2 = w z 2 π erf ( v ) 2 v exp ( v 2 ) .
f r ( r ) = r σ s 2 exp ( r 2 2 σ s 2 ) , r > 0 ;
f h p ( h p ) = g 2 A 0 g 2 h p g 2 1 , 0 h p A 0 ;
f h ( h ) = f h | h a ( h | h a ) f h a ( h a ) d h a ,
f h | h a ( h | h a ) = 1 h a f h p ( h h a ) = g 2 A 0 g 2 h a ( h h a ) g 2 1 , 0 h A 0 h a ;
f h ( h ) = g 2 A A 0 g 2 h g 2 1 k = 1 β a k h / A 0 h a α + k 2 1 g 2 K α k ( 2 α β h a γ β + Ω ) d h a
G 0 , 2 2 , 0 ( x | a , b ) = 2 x ( a + b ) / 2 K a b ( 2 x ) ,
f h ( h ) = g 2 A 2 A 0 g 2 h g 2 1 k = 1 β a k h | A 0 h a α + k 2 1 g 2 G 0 , 2 2 , 0 ( α β h a γ β + Ω | α k 2 , α k 2 ) d h a .
f h ( h ) = g 2 A 2 h 1 k = 1 β a k ( α β γ β + Ω ) α + k 2 G 1 , 3 3 , 0 ( α β γ β + Ω h A 0 | g 2 + 1 g 2 , α , k ) ,
x r G p , q m , n ( x | a p c q ) = G p , q m , n ( x | a p + r c q + r )
f h ( h ) = g 2 A ( G ) 2 h 1 k = 1 a k ( G ) ( α γ ) α + k 2 G 1 , 3 3 , 0 ( α γ h A 0 | g 2 + 1 g 2 , α , k ) .
m k ( h a ) E [ h a k ] = E [ X k ] E [ Y k ] = m k ( X ) m k ( Y ) ,
m k ( X ) = Γ ( α + k ) Γ ( α ) α k ,
m k ( Y ) = 1 γ ( γ β γ β + Ω ) β r = 0 β 1 ( β 1 r ) 1 r ! ( Ω γ ( γ β + Ω ) ) r Γ ( k + r + 1 ) ( β γ β + Ω ) k + r + 1 .
m k ( h a ) = Γ ( α + k ) Γ ( α ) α k γ ( γ β γ β + Ω ) β r = 0 β 1 ( β 1 r ) 1 r ! ( Ω γ ( γ β + Ω ) ) r Γ ( k + r + 1 ) ( β γ β + Ω ) k + r + 1 ,
m k ( Y ) = ( γ β γ β + Ω ) β γ k Γ ( k + 1 ) 2 F 1 ( k + 1 , β ; 1 ; Ω γ β + Ω ) ,
m k ( G ) ( h a ) = Γ ( α + k ) Γ ( α ) α k ( γ β γ β + Ω ) β γ k Γ ( k + 1 ) 2 F 1 ( k + 1 , β ; 1 ; Ω γ β + Ω ) .
m k ( h ) E [ h k ] = E [ h a k ] E [ h p k ] = m k ( h a ) m k ( h p ) .
m k ( h p ) = E [ h p k ] = h p k f h p ( h p ) d h p .
m k ( h p ) = g 2 A 0 g 2 0 A 0 h p g 2 + k 1 d h p = g 2 g 2 + k A 0 k .
m k ( h ) = Γ ( α + k ) Γ ( α ) α k γ ( γ β γ β + Ω ) β g 2 g 2 + k A 0 k × × r = 0 β 1 ( β 1 r ) 1 r ! ( Ω γ ( γ β + Ω ) ) r Γ ( k + r + 1 ) ( β γ β + Ω ) k + r + 1 .
m k ( h ) = g 2 A 0 k g 2 + k Γ ( α + k ) Γ ( α ) α k ( γ β γ β + Ω ) β γ k Γ ( k + 1 ) 2 F 1 ( k + 1 , β ; 1 ; Ω γ β + Ω ) .
P b ( e ) = p 0 p ( e | 0 ) + p 1 p ( e | 1 ) ,
P b ( e | h ) = p ( e | 0 , h ) = p ( e | 1 , h ) = Q ( P R h σ N )
P b ( e ) = 0 P b ( e | h ) f h ( h ) d h ,
G 1 , 2 2 , 0 ( z | a a 1 , a 1 2 ) = π z a 1 erfc ( z ) .
P b ( e ) = 2 α g 2 A B α 2 32 π π k = 1 β 2 k B k 2 a k G 7 , 4 2 , 6 ( 8 R 2 P 2 A 0 2 σ N 2 B 2 | 1 g 2 2 , 2 g 2 2 , 1 α 2 , 2 α 2 , 1 k 2 , 2 k 2 , 1 0 , 1 2 , g 2 2 , 1 g 2 2 )
P b ( e ) = 2 α g 2 A ( G ) 32 π π ( γ α ) α 2 k = 1 2 k ( γ α ) k 2 a k ( G ) × × G 7 , 4 2 , 6 ( 8 R 2 P 2 A 0 2 σ N 2 γ 2 α 2 | 1 g 2 2 , 2 g 2 2 , 1 α 2 , 2 α 2 , 1 k 2 , 2 k 2 , 1 0 , 1 2 , g 2 2 , 1 g 2 2 ) .
σ R 2 = 1.23 C n 2 k 7 / 6 L 11 / 6 .

Metrics