Performance Evaluation ofMPPM-CodedWirelessOpticalMIMO System with Combined Effects over Correlated Fading Channel

-e performance of wireless optical MIMO system with multiple pulse position modulation (MPPM) over correlated fading channel is investigated. -e combined effects of atmospheric attenuation, atmospheric turbulence, and pointing error are taken into consideration.-e bit error rate (BER) and the ergodic channel capacity are analyzed by utilizing the Poisson counting model and the exponential correlation model. Moreover, their approximate expressions are derived. -e simulation results demonstrate that the pointing error is the most prominent influence factor over weak correlated channel.-e performance degradation caused by a high channel correlation coefficient is more than that of pointing error in strong correlated channel. -erefore, the use of pointing, acquisition, and tracking (PAT) system and reasonable arrangement of the number and spacing of antennas at the transceiver are the keys to improve system performance.


Introduction
Wireless optical communication (WOC) has the advantages of both optical fiber links and wireless communications which has attracted a lot of attention in recent years [1,2]. One of the major challenges in WOC is channel fading caused by atmospheric attenuation and atmospheric turbulence [3]. Another major challenge is pointing error which cannot be neglected among the fading due to dynamic wind loads, thermal expansion, and internal vibrations [4]. To eliminate the fading caused by atmospheric channel and pointing error, a large fraction of research has focused on multiple-input multiple-output (MIMO) technology to improve the system performance [5]. e bit error rate (BER) performance and ergodic channel capacity of the WOC MIMO system with atmospheric attenuation and atmospheric turbulence have been studied in references [6,7]. Furthermore, the influence of misalignment characterized by the pointing error on outage probability, channel capacity, and average BER has been studied in references [8][9][10].
However, the aforementioned investigations have neglected the spatial correlation of MIMO subchannels. In practical WOC MIMO applications, the true benefits provided by multiple lasers may be limited by spatial correlation [11,12]. e spatial correlation is quantified as a function of beam separation, turbulence strength, propagation distance, and receiver aperture for the multiple-input single-output (MISO) system in [13]. Moreover, the relationship of the correlation coefficient and receiving aperture size for the MIMO system has been analyzed in reference [14]. ereafter, the BER and the outage probability of correlated gamma-gamma fading channel have been studied in reference [15]. e results indicate that the system performance is very sensitive to the presence of spatial correlation. erefore, spatial correlation needs to be considered to analyze WOC MIMO performance.
On the contrary, background radiation is one of the major noise sources in WOC MIMO links, and the background radiation that impacts on BER and outage capacity of typical system has been considered in reference [16]. In addition, most of the researches utilize on-off keying (OOK) and pulse position modulation (PPM) as their modulation format. As is known, OOK has lower transmission efficiency and is difficult to satisfy the requirement of high-speed transmission [17,18]. PPM has lower bandwidth efficiency and higher synchronization requirement [19]. However, MPPM gets much attention due to its advantages of higher optical power utilization and background radiation tolerance than OOK and PPM in high-speed intensity modulation/direct detection (IM/DD) communications [20]. erefore, it is necessary to quantify the performance of MPPM-coded WOC MIMO systems. e symbol error rate (SER) or the average BER performance of the MIMO system with the MPPM has been analyzed over turbulence fading channel in references [21][22][23]. But the spatial correlation and the background radiation are not considered in references [20][21][22][23].
Based upon the abovementioned analysis, it could be concluded that the factors that can make impacts on the WOC MIMO applications include atmospheric attenuation, atmospheric turbulence, pointing error, spatial correlation, background radiation, and modulation format. However, to the best of our knowledge, the performance of the MPPMcoded WOC MIMO system under combined effects over correlated fading channel has not been considered in previous studies. Concerning this issue, we focus on the BER performance and the ergodic channel capacity of the MPPMcoded WOC MIMO system over correlated log-normal channel when taking the background radiation and the combined effects into consideration. Also, their approximate expressions with low complexity are derived to achieve the performance bound of the WOC MIMO system in practical communication link. e rest of this paper is organized as follows: in Section 2, theoretical analysis is proposed. In Section 3, BER and ergodic channel capacity are analyzed, and their approximate expressions with combined effects over correlated fading channel are derived. Numerical simulation results are given in Section 4, where the Monte Carlo method is utilized to investigate the performance of our proposal. Finally, a conclusion is given in Section 5.

Correlated Channel.
e channel coefficient under the combined effects of atmospheric attenuation, atmospheric turbulence, and pointing errors can be written as where h l � exp(− σL) is the laser power attenuation with transmission distance L, σ is the attenuation coefficient, h l can be considered as a constant scaling factor during a long period of time [24], h a � exp(χ) is the turbulence fading coefficient, and χ is a Gaussian random variable with mean μ χ and variance σ 2 χ . e variance σ 2 χ � − 2μ χ can be rewritten as σ 2 χ � 1.23C 2 n k 7/6 L 11/6 [25]. C 2 n is the refractive-index structure parameter, k � 2π/λ is the wave number, and λ is the wavelength. h p ≈ A 0 exp(− 2R 2 /ω 2 ) is the loss caused by zero-boresight misalignment [8], A 0 � [erf(v)] 2 is the fraction of the collected power at R � 0, andR is the radial displacement at the receiver plane.
, ω L is the beam waist at distance L, and a is the receiving aperture radius.
It should be noted that in practical WOC MIMO systems, there is a certain spatial correlation among subchannels. erefore, the channel coefficient matrix of the correlated channel can be defined as where H � [h nm ] NM is the intensity fading coefficient matrix when the subchannels are independent, R r is the N × N receiver correlation matrix, and R t is the M × Mtransmitter correlation matrix. In order to make a simple description of spatial correlation, the exponential correlation model [26] can be used to define R r and R t as R r � [r |c− b| r ] c,b�1,...,N and R t � [r |c− b| t ] c,b�1,...,M , where r r (0 ≤ r r ≤ 1) and r t (0 ≤ r t ≤ 1) are the correlation coefficients of receiver and transmitter, respectively.
When the transmitter and receiver are correlated (r r ≠ 0, r t ≠ 0), the channel coefficient matrix can then be expressed as

MIMO Channel. A WOC MIMO system is composed of
M lasers diodes (LDs) and N photodetectors (PDs), as shown in Figure 1. Assume the total transmitting power is E s , then the average power of each laser is E s /M. When using the Q − ary MPPM modulation format, a symbol period is divided into Q slots. Each symbol contains ζ "on" slots (ζ ≤ Q) and , where A m denotes light intensity in an "on" slot and 0 denotes an "off" slot. Assume the mean of counted photons in "on" slots is λ on|g nm and in "off" slots is λ off , respectively. When the Poisson statistics model and the equal gain combination (EGC) method are utilized at the receiver end, then there are where η is the photoelectric conversion efficiency, f is the center frequency of optical carriers, ε is the Planck constant, E b is the noise power at the PD, and g nm is the channel coefficient from the m − th LD to the n − th PD. It can be seen from equation (4) that the simplification of M m�1 N n�1 g nm is important for the effective calculation of photons at the receiver. When the subchannels are correlated, the assuming that variable D is equal to the M m�1 N n�1 g nm in equation (4) and D should be written as where l 0 � exp(− σL), δ � 2R 2 /ω 2 and its probability density , and σ s is the standard deviation of the jitter offset at the receiver. e jitter variance σ 2 s is generally used to define the intensity of pointing error. K nm � χ nm − 2(R 2 nm − R 2 )/ω 2 follows normal distribution; therefore, its mean and variance can be written as , respectively, where ϕ n (ϕ n ∈ ψ, ψ is the coordinate set of all receiving apertures) denotes the central coordinate vector of the n − th received laser beam.
On the contrary, when K nm multiplies R t and R r , the result still follows the normal distribution. erefore, variable D can be expressed as where V k follows normal distribution, and its mean and respectively. e sum of log-normal variables can be approximated to another log-normal variable [27], that is, where U � W − δ, W is a normal distributed variable and their PDF can be written as As a result, λ on|g nm can be simplified as e probability distribution of U can be expressed as where Substituting f δ (δ) and equation (10) into equation (9), the result is erefore, the PDF of variable D can then be expressed as It should be noted that the cumulative distribution of intensity attenuation coefficients is no longer a log-normal distribution under the combined effects. It follows a distribution composed of coefficients, complementary error function, and exponential function. International Journal of Antennas and Propagation observations, where z nq denotes the observed photoelectrons of n − thPD in time slot q. When employing the EGC method, there is z q � N n�1 z nq . If one of the Q ζ binary patterns in the x k is launched, then the value detected by ML can be written as

Performance Analysis
where Q (k) on and Q (k) off denote the k − th"on" and "off" slot, respectively. z q !, exp(− λ on|g nm ), and exp(− λ off ) are constant to x k . erefore, equation (13) can be simplified by taking the logarithm as Since M m�1 N n�1 g nm is simplified to D, equation (14) can be written as e results and discussion may be presented separately, or in one combined section, and may optionally be divided into headed subsections.

Bit Error Rate
Performance. Assume the signal passes through a slow fading channel, which means the channel fading coefficient is constant in a symbol period. Taking the background radiation into consideration, the possible errors fall into two categories: the definite error and the indefinite error.

Definite Error Probability.
e definite error probability can be defined as that the number of photoelectrons detected in an "off" slot is higher than that in an "on" slot. Let S on,l , l � 1, . . . , ζ and S off,l , l � 1, . . . , Q − ζ denote the sum of photoelectrons count in "on" and "off" slots, respectively. e definite error probability can be expressed as Since the received photoelectrons in "on" and "off" slots follow Poisson distribution, which can be described as us, equation (16) can be rewritten as

Indefinite Error
Probability. e indefinite error probability means that one or more of the "off" slots has the same photoelectrons count as the lowest-valued "on" slots.
As a result, the receiver is likely to make a wrong decision. e indefinite error probability can be expressed as As a result, the conditional probability of symbol error with background radiation and channel fading can be written as International Journal of Antennas and Propagation Equation (20) shows the conditional probability of symbol error depends on the sum of channel fading coefficients M m�1 N n�1 g nm . erefore, the theoretical probability of symbol error is Because the sum of channel coefficients g nm can be simplified to a variable D, the approximate probability of symbol error with background radiation can be written as Considering that P b � [Q/2(Q − 1)]P s , the approximate BER can be expressed as Equation (23) reveals that the approximate BER upper bound related to the parameters of Q, ζ, M, N, E s , E b , c, and D. It should be noted that the parameters l 0 , r t , r r , σ 2 s , and σ 2 χ nm are factors which can influence D.

Ergodic Channel Capacity.
According to equation (4), the channel capacity of the wireless optical MIMO system can be defined by where x and y denote the launched and received signal. e instantaneous channel capacity can be written as Considering the channel coefficient is a random variable. erefore, the ergodic channel capacity is usually utilized to analyze the system capacity. As a result, the theoretical MPPM-coded ergodic channel capacity of the WOC MIMO system can be expressed as 6 International Journal of Antennas and Propagation It is complex to perform the M × N − fold integration. As mentioned in Section 2.2, the sum of M m�1 N n�1 g nm can be approximated as variable D. us, equation (26) can be simplified as e approximation is much simpler than its former version since it only requires a single integration. Substituting equation (12) into equation (27), the approximate ergodic channel capacity under combined effects can be further expressed as Equation (28) is a simplification of equation (26). It is evident from equation (28) that the ergodic channel capacity is related to parameters as Q, ζ, ε, f, E s , M, and D, where l 0 , σ 2 χ nm , σ 2 s and r r , r t are factors which can influence D.

Numerical Investigation
In this section, the BER and the ergodic channel capacity of the WOC MIMO system with combined effects over correlated fading channel have been numerically investigated based on the foregoing analysis. Especially, the BER curves are obtained by Monte Carlo simulations, and the curves of ergodic channel capacity are based on the equations. Moreover, all the derived expressions are verified via simulations. Numerical simulations are carried out by MATLAB with parameters as η � 0.5, Q � 5, ζ � 2, λ � 1550 nm, l 0 � 0.9, L � 1000 m, E b � − 170 dBJ, σ 2 s � 0.1, C 2 n � 2.23 × 10 − 14 , and σ 2 χ nm � 0.3. e turbulence fading and misalignment fading are characterized by parameters σ 2 χ nm and σ 2 s . Table 1 demonstrates the 5-ary MPPM symbol constellation with 2 active time slots. It is evident from the table that only 2 time slots need to be activated when sending 3 bit information during a symbol period. Figure 2 shows the BER curves under combined effects, where CE denotes the condition of combined effects in independent identically distributed (iid) channel, cor denotes the correlated channel, appro denotes the approximate value of the proposed method, and exact denotes the theoretical value. It is observed that the approximate BER performances are slightly worse than the corresponding theoretical value. Note that the computation complexity is nearly 75% lower than that of the theoretical one. erefore, such an approximation is appropriate for practical application. Figure 3 shows the approximate BER performance of a 2 × 2 WOC MIMO system, where Per denotes the ideal condition, AT denotes the condition only with atmospheric turbulence, and AT&AA denotes the condition with both the atmospheric turbulence and atmospheric attenuation. e result shows that the CE seriously impacts on the system performance. When BER � 10 − 4 , the required energy increases by about 24 dBJ compared with the Per case. It should be noted that the pointing error leads to the most significant impact on the BER performance among all the influencing factors. For example, when BER � 10 − 4 , the required energy for CE is 22.5 dBJ more than that of AT&AA. erefore, the pointing error is the most prominent influence factor over independent channel. Figure 4 demonstrates the approximate BER curves versus E s under different correlation coefficients in a 2 × 2 WOC MIMO system. It can be observed that the higher the International Journal of Antennas and Propagation correlation coefficient, the worse the BER performance. When BER � 10 − 4 and the correlated coefficients are 0.3, 0.5, and 0.8, the required energy increases by about 10 dBJ, 15 dBJ, and 28 dBJ, respectively, compared with independent channel. Particularly, the BER goes to its maximum when the correlation coefficient is 1, which can be regarded as the upper bound of the system BER. Moreover, it is evident that the strong channel correlation can seriously impact on the system BER performance. For example, when BER � 10 − 4 , the power loss caused by the combined effects is 24 dBJ, while the power loss caused by strong correlation (r ≥ 0.8) is 28 dBJ. erefore, high channel correlation coefficientcaused performance degradation is more than that of combined effects in strong correlated channel. According to Figures 3 and 4, the BER performance of the WOC MIMO system is seriously affected by pointing error over independent and weak correlated channel, while the spatial correlation is the most prominent influence factor over strong correlated channel. Figure 5 shows the approximate BER performance under combined effects over correlated channel with different turbulence (σ 2 χ nm ) and pointing error jitter variance (σ 2 s ). It is evident from the plot that the increased σ 2 χ nm and σ 2 s lead to worse BER, where σ 2 s plays a more prominent role in the process. For example, at BER � 10 − 4 , the required energy increases by about 0.1 dBJ and 10 dBJ, when σ 2 χ nm and σ 2 s increase from 0.1 to 0.3, respectively. erefore, the increase of σ 2 s has a greater impact on the BER performance of the WOC MIMO system than the increase of σ 2 χ nm . Figure 6 shows the BER performance of PPM and MPPM formats in a 2 × 2 MIMO system. It is evident from the plot that the BER performance of (Q � 5, ζ � 2) MPPM is better than that of 4-PPM and 8-PPM formats over correlated channel. e reason for this is that, for activated pulse combination of MPPM, the redundancy resulting from the constellation causes stronger anti-interference ability. In addition, the spectrum efficiency of (Q � 5, ζ � 2) MPPM is 0.6 bits/s/Hz, which is 0.1 bits/s/Hz and 0.225 bits/s/Hz higher than that of 4-PPM and 8-PPM, respectively. Figure 7 shows the ergodic channel capacity of PPM and MPPM formats in a 2 × 2 MIMO system. It is clear that the ergodic channel capacity of (Q � 5, ζ � 2) MPPM is 3.4 bits/ channel, which is 1.4 bits/channel and 0.4 bits/channel higher than that of 4-PPM and 8-PPM, respectively. Comparing the BER, ergodic channel capacity, and spectrum efficiency of PPM and MPPM, it can be obtained that MPPM has advantages in the wireless optical MIMO communication system. Figure 8 shows the theoretical and approximate ergodic channel capacity with different number of "on" slots. Both theoretical (equation (26)) and approximate (equation (28)) expressions are depicted in the figure. It is observed that the performance of approximation is slightly worse than the corresponding theoretical one. Note that the computation complexity of approximation is much lower. erefore, the slight performance loss is acceptable for the approximation.    International Journal of Antennas and Propagation Figure 9 shows the approximate ergodic channel capacity of 2 × 2 and 4 × 4 WOC MIMO system over independent channel. It is evident from the plot that the CE will seriously impact on the ergodic channel capacity. For example, when E s � − 190 dBJ, the ergodic channel capacity reduced by 2.64 bits/channel compared with the Per case in the 4 × 4 system. It should be noted that the influence of pointing error is more significant than atmospheric attenuation and atmospheric turbulence. e reduction caused by pointing error is about 2.27 bits/channel. Figure 10 shows     − 135 dBJ, respectively. Furthermore, when the ergodic channel capacity reaches its saturation, the power loss caused by the combined effects is 20 dBJ, while the power loss caused by the strong correlation is 22.5 dBJ, which is greater than that of combined effects. erefore, it can be considered that strong correlation has more effect on ergodic channel capacity than combined effects. Figure 11 shows the approximate ergodic channel capacity with different atmospheric turbulence and pointing error jitter variance over correlated channel. e ergodic channel capacity decreases with the increasing of σ 2 χ nm and σ 2 s , where σ 2 s plays a more important role in the process. For example, at E s � − 172.5 dBJ, the ergodic channel capacity decreases by about 0.034 bits/channel and 0.467 bits/channel, when σ 2 χ nm and σ 2 s increase from 0.1 to 0.3, respectively. As a result, the increase of σ 2 s has a greater impact on the ergodic channel capacity of the WOC MIMO system than the increase of σ 2 χ nm .

Conclusions
e MPPM-coded WOC MIMO system performance over correlated fading channel under combined effects of atmospheric attenuation, atmospheric turbulence, and pointing error is discussed. e approximate expressions of BER and ergodic channel capacity are proposed to reduce the computational complexity. e research results indicate that (1) the pointing error is the prominent influence factor in independent and weak correlated channel and (2) the performance degradation caused by high channel correlation coefficient is more than that of pointing error in strong correlated channel. erefore, in practical applications, the number and spacing of lasers should be reasonably designed to minimize the spatial correlation of subchannels. Additionally, the system performance can be further improved by   Data Availability e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper. International Journal of Antennas and Propagation 11