International Journal of Optics

Volume 2016 (2016), Article ID 2613816, 8 pages

http://dx.doi.org/10.1155/2016/2613816

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

School of Astronautics and Aeronautics, University of Electronic Science and Technology of China, 2006 Xiyuan Ave, West Hi-Tech Zone, Chengdu 611731, China

Received 26 May 2016; Accepted 3 August 2016

Academic Editor: Sulaiman Wadi Harun

Copyright © 2016 Chao Gao and Xiaofeng Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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.

#### 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]where is the angular wavenumber of the turbulence scale, is the spectral power law value, is the generalized atmospheric structure parameter, is the inner scale of turbulence, and is the outer scale of turbulence. in (1) is a function related to :where is the gamma function.

For the convenience of mathematical analysis, letwhere the constant coefficients , , and in (3) are usually set asIt 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, in (1) takes the form whereAnd the coefficients are , , , , , , , and .

##### 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]where is the scalar separation between two observation points and is the propagation optical path length. in (7) is the angular wavenumber of Gaussian beam wavewhere is the wavelength of Gaussian beam. Both and in (7) are optical parameters of the Gaussian beam at the receiverwhere is the curvature parameter of Gaussian beam at transmitter and is the Fresnel ratio of Gaussian beam at transmitterIn (10), is the phase front radius of Gaussian beam at transmitter, and is the radius of Gaussian beam at transmitter. in (7) is the modified Bessel function of the first kind with zero order, and in (7) is the Bessel function of the first kind with zero order [6]The atmospheric turbulence MTF takes the form [1]where is the normalized spatial frequency and 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 thatFor mathematical convenience, letThus, (13) can be presented by

##### 3.1. Reduction of

Substituting (11) into (14), is rewritten asIn most situations, . 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 expressionConsider the iterated integral in (17). According to (5), there isBased on the equation for and [7],we can getWithout loss of generality, the integrand in (17) takes the formwhere and . Based on the equation for and [7],we can getwhere is the Gaussian hypergeometric function [6]. Thus, can be computed by (17) and (23) with .

##### 3.2. Reduction of

Following similar procedures as presented in Section 3.1, in (14) is rewritten asExpanding (24) by the binomial theorem, it follows thatThus, could be computed by (25) and (23) with .

#### 4. Numerical Simulations

The following simulations are conducted by the Gaussian beam with these settings: m, m, rad/m, , m, , and m. Of course, other values can also be chosen.

Figure 1 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 m and m, respectively. As shown in Figure 1(a), the atmospheric turbulence apparently produces more effects on the propagation of the convergent Gaussian beam () 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 () in Figure 1(b) and the divergent Gaussian beam () in Figure 1(c).