# Modulation Transfer Function of a Gaussian Beam Based on the Generalized Modified Atmospheric Spectrum

Authors: Chao Gao; Xiaofeng Li
Journal: International Journal of Optics
 (2016)

---

## Abstract

This paper investigates the modulation transfer function of a Gaussian beam propagating through a horizontal path in weak-fluctuation non-Kolmogorov turbulence. Mathematical expressions are obtained based on the generalized modified atmospheric spectrum, which includes the spectral power law value of non-Kolmogorov turbulence, the finite inner and outer scales of turbulence, and other optical parameters of the Gaussian beam. The numerical results indicate that the atmospheric turbulence would produce less negative effects on the wireless optical communication system with an increase in the inner scale of turbulence. Additionally, the increased outer scale of turbulence makes a Gaussian beam influenced more seriously by the atmospheric turbulence.


---

## Body

## 1. Introduction

Optical wireless communication technology has drawn much attention for its significant technological challenges and prospective applications. It uses beams of laser propagating in the atmosphere to wirelessly transmit data at high speed. However, the atmosphere is full of numerous turbulence eddies, which has great degrading impacts on the performance of the communication system. The degrading effects of atmospheric turbulence on the communication system can be characterized statistically by the modulation transfer function (MTF) [1]. In the past few decades, various power spectrum models of refractive index have been proposed to analyze the MTF for different situations. Generally speaking, these turbulence power spectrum models can be classified into two typical categories: Kolmogorov and non-Kolmogorov models. The former have a fixed power law value of 11/3, while the latter allow the power law value to vary in the range from three to four. Most non-Kolmogorov models can be generalized from their corresponding Kolmogorov models, and thus the Kolmogorov models can be regarded as specific cases of the non-Kolmogorov models [2]. Among these models, the generalized modified atmospheric spectrum not only considers the variable spectral power law value between the ranges from 3 to 4, but also takes the finite inner and outer scales of turbulence into account [3]. Besides, the generalized modified atmospheric spectrum features the small rise at a high wavenumber, which is clearly seen in temperature data recorded by sensors. These properties make the generalized modified atmospheric spectrum suitable and unique in the investigation of the MTF for plane and spherical waves [4].In this study, the generalized modified atmospheric spectrum is used to investigate the MTF of a Gaussian beam in non-Kolmogorov turbulence along a horizontal path. The Gaussian beam, whose transverse electric field and intensity are normally distributed, is a typical kind of electromagnetic wave [5]. The rest of the paper is organized as follows. Section 2 introduces the generalized modified atmospheric spectrum and the MTF of a Gaussian beam. Section 3 presents a detailed expression reduction. The influences of the inner and outer scales of turbulence on the MTF of a Gaussian beam are analyzed in Section 4, followed by conclusions in Section 5.


## 2. Theoretical Models

### 2.1. Generalized Modified Atmospheric Spectrum

The generalized modified atmospheric spectrum takes the form [3](1)
Φ
n
κ
=
A
α
C
n
2
κ
-
α
f
κ
,
α
,where κ
∈
0
,
+
∞ is the angular wavenumber of the turbulence scale, α
∈
3,4 is the spectral power law value, C
n
2 is the generalized atmospheric structure parameter, l
0
≥
0 is the inner scale of turbulence, and L
0
≥
l
0 is the outer scale of turbulence. A
α in (1) is a function related to α:(2)
A
α
=
Γ
α
-
1
4
π
2
sin
⁡
π
α
-
3
2
,where Γ
x is the gamma function.For the convenience of mathematical analysis, let(3)
C
α
=
3
-
α
3
Γ
3
-
α
2
+
a
4
-
α
3
Γ
4
-
α
2
-
b
3
+
β
-
α
3
Γ
3
+
β
-
α
2
,where the constant coefficients a, b, and β in (3) are usually set as(4)
a
=
1.802
,
b
=
0.254
,
β
=
7
6
.It must be pointed out that the values of these coefficients are based on the experiments for the classic Kolmogorov turbulence but are widely used for theoretical analyses of non-Kolmogorov turbulence [1, 4]. Nevertheless, f
κ
,
α in (1) takes the form (5)
f
κ
,
α
=
exp
⁡
-
κ
2
κ
l
2
×
1
-
exp
⁡
-
κ
2
κ
0
2
×
1
+
a
κ
κ
l
-
b
κ
κ
l
β
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
κ
p
i
exp
⁡
-
d
j
2
κ
2
,where(6)
κ
0
=
4
π
L
0
,
κ
l
=
π
A
α
C
α
1
/
α
-
5
l
0
.And the coefficients are c
1
=
1, c
2
=
a
/
κ
l, c
3
=
-
b
/
κ
l
β, p
1
=
0, p
2
=
1, p
3
=
β, d
1
=
1
/
κ
l
2, and d
2
=
1
/
κ
l
2
+
1
/
κ
0
2.


### 2.2. MTF of a Gaussian Beam

The MTF is relative to the wave structure function (WSF). Based on the Rytov approximation, the WSF of Gaussian beam takes the simple form [1](7)
D
ρ
=
8
L
k
2
π
2
∫
0
1
d
ξ
∫
0
+
∞
d
κ
×
κ
Φ
n
κ
exp
⁡
-
L
Λ
κ
2
ξ
2
k
×
I
0
Λ
ρ
κ
ξ
-
J
0
ρ
κ
1
-
Θ
¯
ξ
,where ρ is the scalar separation between two observation points and L is the propagation optical path length. k in (7) is the angular wavenumber of Gaussian beam wave(8)
k
=
2
π
λ
,where λ is the wavelength of Gaussian beam. Both Λ and Θ
¯ in (7) are optical parameters of the Gaussian beam at the receiver(9)
Θ
=
Θ
0
Θ
0
2
+
Λ
0
2
,
Θ
¯
=
1
-
Θ
,
Λ
=
Λ
0
Θ
0
2
+
Λ
0
2
,where Θ
0 is the curvature parameter of Gaussian beam at transmitter and Λ
0 is the Fresnel ratio of Gaussian beam at transmitter(10)
Θ
0
=
1
-
L
R
0
,
Λ
0
=
2
L
k
W
0
.In (10), R
0 is the phase front radius of Gaussian beam at transmitter, and W
0 is the radius of Gaussian beam at transmitter. I
0
x in (7) is the modified Bessel function of the first kind with zero order, and J
0
x in (7) is the Bessel function of the first kind with zero order [6](11)
I
0
x
=
∑
n
=
0
+
∞
1
n
!
2
x
2
2
n
,
J
0
x
=
∑
n
=
0
+
∞
-
1
n
n
!
2
x
2
2
n
.The atmospheric turbulence MTF takes the form [1](12)
M
T
F
μ
=
exp
⁡
-
1
2
D
μ
d
,where μ is the normalized spatial frequency and d is the receiver aperture diameter. It is clear that the value range of the MTF is the interval from 0 to 1.


## 2.1. Generalized Modified Atmospheric Spectrum

The generalized modified atmospheric spectrum takes the form [3](1)
Φ
n
κ
=
A
α
C
n
2
κ
-
α
f
κ
,
α
,where κ
∈
0
,
+
∞ is the angular wavenumber of the turbulence scale, α
∈
3,4 is the spectral power law value, C
n
2 is the generalized atmospheric structure parameter, l
0
≥
0 is the inner scale of turbulence, and L
0
≥
l
0 is the outer scale of turbulence. A
α in (1) is a function related to α:(2)
A
α
=
Γ
α
-
1
4
π
2
sin
⁡
π
α
-
3
2
,where Γ
x is the gamma function.For the convenience of mathematical analysis, let(3)
C
α
=
3
-
α
3
Γ
3
-
α
2
+
a
4
-
α
3
Γ
4
-
α
2
-
b
3
+
β
-
α
3
Γ
3
+
β
-
α
2
,where the constant coefficients a, b, and β in (3) are usually set as(4)
a
=
1.802
,
b
=
0.254
,
β
=
7
6
.It must be pointed out that the values of these coefficients are based on the experiments for the classic Kolmogorov turbulence but are widely used for theoretical analyses of non-Kolmogorov turbulence [1, 4]. Nevertheless, f
κ
,
α in (1) takes the form (5)
f
κ
,
α
=
exp
⁡
-
κ
2
κ
l
2
×
1
-
exp
⁡
-
κ
2
κ
0
2
×
1
+
a
κ
κ
l
-
b
κ
κ
l
β
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
κ
p
i
exp
⁡
-
d
j
2
κ
2
,where(6)
κ
0
=
4
π
L
0
,
κ
l
=
π
A
α
C
α
1
/
α
-
5
l
0
.And the coefficients are c
1
=
1, c
2
=
a
/
κ
l, c
3
=
-
b
/
κ
l
β, p
1
=
0, p
2
=
1, p
3
=
β, d
1
=
1
/
κ
l
2, and d
2
=
1
/
κ
l
2
+
1
/
κ
0
2.


## 2.2. MTF of a Gaussian Beam

The MTF is relative to the wave structure function (WSF). Based on the Rytov approximation, the WSF of Gaussian beam takes the simple form [1](7)
D
ρ
=
8
L
k
2
π
2
∫
0
1
d
ξ
∫
0
+
∞
d
κ
×
κ
Φ
n
κ
exp
⁡
-
L
Λ
κ
2
ξ
2
k
×
I
0
Λ
ρ
κ
ξ
-
J
0
ρ
κ
1
-
Θ
¯
ξ
,where ρ is the scalar separation between two observation points and L is the propagation optical path length. k in (7) is the angular wavenumber of Gaussian beam wave(8)
k
=
2
π
λ
,where λ is the wavelength of Gaussian beam. Both Λ and Θ
¯ in (7) are optical parameters of the Gaussian beam at the receiver(9)
Θ
=
Θ
0
Θ
0
2
+
Λ
0
2
,
Θ
¯
=
1
-
Θ
,
Λ
=
Λ
0
Θ
0
2
+
Λ
0
2
,where Θ
0 is the curvature parameter of Gaussian beam at transmitter and Λ
0 is the Fresnel ratio of Gaussian beam at transmitter(10)
Θ
0
=
1
-
L
R
0
,
Λ
0
=
2
L
k
W
0
.In (10), R
0 is the phase front radius of Gaussian beam at transmitter, and W
0 is the radius of Gaussian beam at transmitter. I
0
x in (7) is the modified Bessel function of the first kind with zero order, and J
0
x in (7) is the Bessel function of the first kind with zero order [6](11)
I
0
x
=
∑
n
=
0
+
∞
1
n
!
2
x
2
2
n
,
J
0
x
=
∑
n
=
0
+
∞
-
1
n
n
!
2
x
2
2
n
.The atmospheric turbulence MTF takes the form [1](12)
M
T
F
μ
=
exp
⁡
-
1
2
D
μ
d
,where μ is the normalized spatial frequency and d is the receiver aperture diameter. It is clear that the value range of the MTF is the interval from 0 to 1.


## 3. Expression Reduction

The calculation equation (7) will spend too much time because of its improper iterated integral. As an alternative, the closed-form expression of (7) can replace the improper iterated integral with special functions, which has corresponding packages in frequently used software. This section mainly discusses the reduction of (7).Substituting (1) into (7), it follows that(13)
D
ρ
=
8
A
α
C
n
2
L
k
2
π
2
∫
0
1
d
ξ
∫
0
+
∞
d
κ
×
κ
1
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
×
I
0
Λ
ρ
κ
ξ
-
J
0
ρ
κ
1
-
Θ
¯
ξ
.For mathematical convenience, let(14)
D
I
=
∫
0
1
d
ξ
∫
0
+
∞
d
κ
×
κ
1
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
×
I
0
Λ
ρ
κ
ξ
-
1
,
D
J
=
∫
0
1
d
ξ
∫
0
+
∞
d
κ
×
κ
1
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
×
J
0
ρ
κ
1
-
Θ
¯
ξ
-
1
.Thus, (13) can be presented by(15)
D
ρ
=
8
A
α
C
n
2
L
k
2
π
2
×
D
I
-
D
J
.


### 3.1. Reduction ofD
I

Substituting (11) into (14), D
I is rewritten as(16)
D
I
=
∫
0
1
d
ξ
∫
0
+
∞
d
κ
×
∑
n
=
1
+
∞
1
n
!
2
Λ
ρ
κ
ξ
2
2
n
×
κ
1
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
.In most situations, ρ
≪
1. This is because MTF will quickly converge to zero when ρ approaches one; that is, MTF is significantly larger than zero when ρ approaches zero. Thus, (16) could be approximated by the simpler expression(17)
D
I
≈
Λ
2
ρ
2
4
∫
0
1
ξ
2
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
.Consider the iterated integral in (17). According to (5), there is(18)
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
κ
3
-
α
+
p
i
×
exp
⁡
-
d
j
2
+
L
Λ
ξ
2
k
κ
2
.Based on the equation for u
>
-
1 and v
>
0 [7],(19)
∫
0
+
∞
x
u
exp
⁡
-
v
x
2
d
x
=
1
2
v
-
u
+
1
/
2
Γ
u
+
1
2
,we can get(20)
∫
0
+
∞
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
d
κ
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
2
d
j
2
+
L
Λ
ξ
2
k
-
4
-
α
+
p
i
/
2
×
Γ
4
-
α
+
p
i
2
.Without loss of generality, the integrand in (17) takes the form(21)
ξ
n
∫
0
+
∞
κ
p
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
d
κ
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
2
Γ
1
+
p
+
p
i
2
×
L
Λ
k
-
1
+
p
+
p
i
/
2
ξ
n
k
d
j
2
L
Λ
+
ξ
2
-
1
+
p
+
p
i
/
2
,where n
=
2 and p
=
3
-
α. Based on the equation for u
>
0 and w
>
0 [7],(22)
∫
0
1
x
u
-
1
w
2
+
x
2
v
d
x
=
1
u
w
2
v
F
1
2
-
v
,
u
2
;
u
+
2
2
;
-
1
w
2
,we can get(23)
∫
0
1
ξ
n
d
ξ
∫
0
+
∞
d
κ
×
κ
p
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
2
d
j
1
+
p
+
p
i
n
+
1
×
Γ
1
+
p
+
p
i
2
×
F
1
2
1
+
p
+
p
i
2
,
n
+
1
2
;
n
+
3
2
;
-
L
Λ
k
d
j
2
,where F
1
2
a
,
b
;
c
;
z is the Gaussian hypergeometric function [6]. Thus, D
I can be computed by (17) and (23) with n
=
2.


### 3.2. Reduction ofD
J

Following similar procedures as presented in Section3.1, D
J in (14) is rewritten as(24)
D
J
≈
-
ρ
2
4
∫
0
1
1
-
Θ
¯
ξ
2
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
.Expanding (24) by the binomial theorem, it follows that(25)
D
J
=
-
Θ
¯
2
ρ
2
4
∫
0
1
ξ
2
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
+
Θ
¯
ρ
2
2
∫
0
1
ξ
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
-
ρ
2
4
∫
0
1
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
.Thus, D
J could be computed by (25) and (23) with n
=
0,1
,
2.


## 3.1. Reduction ofD
I

Substituting (11) into (14), D
I is rewritten as(16)
D
I
=
∫
0
1
d
ξ
∫
0
+
∞
d
κ
×
∑
n
=
1
+
∞
1
n
!
2
Λ
ρ
κ
ξ
2
2
n
×
κ
1
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
.In most situations, ρ
≪
1. This is because MTF will quickly converge to zero when ρ approaches one; that is, MTF is significantly larger than zero when ρ approaches zero. Thus, (16) could be approximated by the simpler expression(17)
D
I
≈
Λ
2
ρ
2
4
∫
0
1
ξ
2
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
.Consider the iterated integral in (17). According to (5), there is(18)
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
κ
3
-
α
+
p
i
×
exp
⁡
-
d
j
2
+
L
Λ
ξ
2
k
κ
2
.Based on the equation for u
>
-
1 and v
>
0 [7],(19)
∫
0
+
∞
x
u
exp
⁡
-
v
x
2
d
x
=
1
2
v
-
u
+
1
/
2
Γ
u
+
1
2
,we can get(20)
∫
0
+
∞
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
d
κ
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
2
d
j
2
+
L
Λ
ξ
2
k
-
4
-
α
+
p
i
/
2
×
Γ
4
-
α
+
p
i
2
.Without loss of generality, the integrand in (17) takes the form(21)
ξ
n
∫
0
+
∞
κ
p
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
d
κ
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
2
Γ
1
+
p
+
p
i
2
×
L
Λ
k
-
1
+
p
+
p
i
/
2
ξ
n
k
d
j
2
L
Λ
+
ξ
2
-
1
+
p
+
p
i
/
2
,where n
=
2 and p
=
3
-
α. Based on the equation for u
>
0 and w
>
0 [7],(22)
∫
0
1
x
u
-
1
w
2
+
x
2
v
d
x
=
1
u
w
2
v
F
1
2
-
v
,
u
2
;
u
+
2
2
;
-
1
w
2
,we can get(23)
∫
0
1
ξ
n
d
ξ
∫
0
+
∞
d
κ
×
κ
p
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
=
∑
i
=
1
3
∑
j
=
1
2
-
1
j
-
1
c
i
2
d
j
1
+
p
+
p
i
n
+
1
×
Γ
1
+
p
+
p
i
2
×
F
1
2
1
+
p
+
p
i
2
,
n
+
1
2
;
n
+
3
2
;
-
L
Λ
k
d
j
2
,where F
1
2
a
,
b
;
c
;
z is the Gaussian hypergeometric function [6]. Thus, D
I can be computed by (17) and (23) with n
=
2.


## 3.2. Reduction ofD
J

Following similar procedures as presented in Section3.1, D
J in (14) is rewritten as(24)
D
J
≈
-
ρ
2
4
∫
0
1
1
-
Θ
¯
ξ
2
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
.Expanding (24) by the binomial theorem, it follows that(25)
D
J
=
-
Θ
¯
2
ρ
2
4
∫
0
1
ξ
2
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
+
Θ
¯
ρ
2
2
∫
0
1
ξ
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
-
ρ
2
4
∫
0
1
d
ξ
∫
0
+
∞
d
κ
×
κ
3
-
α
f
κ
,
α
exp
⁡
-
L
Λ
κ
2
ξ
2
k
.Thus, D
J could be computed by (25) and (23) with n
=
0,1
,
2.


## 4. Numerical Simulations

The following simulations are conducted by the Gaussian beam with these settings:λ
=
1.55
×
1
0
-
6 m, L
=
1000 m, k
≈
4.0537
×
1
0
6 rad/m, C
n
2
=
1.7
×
1
0
-
14
m
3
-
α, W
0
=
0.1 m, Λ
0
≈
0.0493, and d
=
0.1 m. Of course, other values can also be chosen.Figure1 depicts the effects of spectral power law value on MTF for different types of Gaussian beams. In this calculation, the inner and outer scales of turbulence are set as l
0
=
0.01 m and L
0
=
5 m, respectively. As shown in Figure 1(a), the atmospheric turbulence apparently produces more effects on the propagation of the convergent Gaussian beam (Θ
0
=
0.5) with an increase in the normalized spatial frequency μ, which acts in accordance with common sense. Besides, from Figure 1(a), it can be found that the non-Kolmogorov atmospheric turbulence would bring more effects on the wireless optical communication system when the spectral power law value α decreases. The same trends are obtained for the collimated Gaussian beam (Θ
0
=
1) in Figure 1(b) and the divergent Gaussian beam (Θ
0
=
2) in Figure 1(c).Figure 1
Effects of spectral power law value on MTF for different types of Gaussian beams.
(a)
ConvergentΘ
0
=
0.5
(b)
CollimatedΘ
0
=
1
(c)
DivergentΘ
0
=
2To analyze the effects of the turbulence inner scale on MTF, the spectral power law value and the outer scale of turbulence are fixed to constant values asα
=
3.6 and L
0
=
5 m. Several inner scales of turbulence are used, and calculation results are depicted in Figure 2 for different types of Gaussian beams. It can be seen that, with an increase in the inner scale of turbulence, the value of MTF also increases. This can be physically explained by the change of inertial subrange of turbulence. When the inner scale of turbulence increases, the frequency’s upper bound of inertial subrange would move to a lower position, and thus the atmospheric turbulence would bring less effects on the propagation of the Gaussian beam.Figure 2
Effects of inner scale on MTF for different types of Gaussian beams.
(a)
ConvergentΘ
0
=
0.5
(b)
CollimatedΘ
0
=
1
(c)
DivergentΘ
0
=
2The influences of outer scale of turbulence on MTF are depicted in Figure3 for different types of Gaussian beams. For the real atmospheric turbulence, the outer scale of turbulence is usually in the order of meters. Hence, it is set to 1 m, 5 m, and 25 m, respectively. The spectral power law value and the inner scale of turbulence are set to α
=
3.6 and l
0
=
0.01 m as example. It can be seen that, with an increase in the outer scale of turbulence, the value of MTF decreases and thus the quality of the Gaussian beam is degraded severely by the atmospheric turbulence. This is because WSF is mostly influenced by the large-scale turbulence eddies, which are relevant to the low-frequency part of the atmospheric turbulence spectrum. A larger turbulence outer scale would lead to a larger range of inertial subrange.Figure 3
Effects of outer scale on MTF for different types of Gaussian beams.
(a)
ConvergentΘ
0
=
0.5
(b)
CollimatedΘ
0
=
1
(c)
DivergentΘ
0
=
2For further discussions and analyses, the inner and outer scales of turbulence are assigned to constant valuesl
0
=
0.01 m and L
0
=
5 m, respectively. The spectral power law value still uses the default value α
=
3.6. Some typical values of wavelength in the near infrared region, λ
=
850 nm, λ
=
1060 nm, and λ
=
1550 nm, are investigated in this simulation. Figure 4 depicts MTF for different Gaussian beams as a function of μ with different λ. It is obvious that the value of MTF increases with an increase in λ for certain type of Gaussian beam if other optical parameters are fixed. This phenomenon may be caused by the fact that the larger the beam wavelength, the more pronounced the diffraction. Thus, a laser beam with larger wavelength can be less affected by turbulence eddies.Figure 4
Effects of wavelength on MTF for different types of Gaussian beams.
(a)
ConvergentΘ
0
=
0.5
(b)
CollimatedΘ
0
=
1
(c)
DivergentΘ
0
=
2


## 5. Conclusions

In this paper, a theoretical expression of the MTF is derived for a Gaussian beam propagating through the non-Kolmogorov atmospheric turbulence along a horizontal path. This expression contains a variable spectral power law value, finite inner and outer scales of turbulence, and other important optical parameters of a Gaussian beam. Numerical simulations indicate that the atmospheric turbulence would produce less degrading effects on the wireless optical communication system with an increase in the spectral power law value. The decreased inner scale of turbulence makes a Gaussian beam influenced more seriously by the atmospheric turbulence. With an increase in the outer scale of turbulence, the quality of a Gaussian beam is degraded more severely by the atmospheric turbulence. A laser beam with larger wavelength can be less affected by turbulence eddies.



---

*Source: 2613816-2016-09-26.xml*