Performance Analysis of Electrical MMSE Linear Equalizers in Optically Amplified OOK Systems
Performance Analysis of Electrical MMSE Linear Equalizers in Optically Amplified OOK Systems
Journal of the Optical Society of Korea. 2011. Sep, 15(3): 232-236
Copyright ©2011, Optical Society of Korea
  • Received : July 21, 2011
  • Accepted : August 02, 2011
  • Published : September 25, 2011
Export by style
Cited by
About the Authors
Jang-Woo, Park
Won-Zoo, Chung

We analyze the linear equalizers used in optically amplified on-off-keyed (OOK) systems to combat chromatic dispersion (CD) and polarization mode dispersion (PMD), and we derive the mathematical minimum mean squared error (MMSE) performance of these equalizers. Currently, the MMSE linear equalizer for optical OOK systems is obtained by simulations using adaptive approaches such as least mean squared (LMS) or constant modulus algorithm (CMA), but no theoretical studies on the optimal solutions for these equalizers have been performed. We model the optical OOK systems as square-law nonlinear channels and compute the MMSE equalizer coefficients directly from the estimated optical channel, signal power, and optical noise variance. The accuracy of the calculated MMSE equalizer coefficients and MMSE performance has been verified by simulations using adaptive algorithms.
As data traffic demands increase rapidly, optical communication systems are widely used, providing high data rates with a large number of channels per fiber. However, optical fiber communication systems suffer from various linear and nonlinear transmission impairments [1] . Chromatic dispersion(CD) and polarization mode dispersion (PMD) are two important factors degrading the performance of high speed optical fiber transmission systems [2] . Intersymbol interference(ISI) caused by CD, PMD, and other impairments increases bit error rate (BER). As the transmission rate increases, ISI mitigation becomes indispensable and many schemes for CD and PMD compensation have been proposed. These compensation techniques can be separated into two approaches;optical domain equalization [3] and electrical domain equalization [2 , 4 , 5] . Optical equalization schemes can achieve perfect compensation using inverse system response, while the performance of electrical domain equalization is limited due to the nonlinear channel effect caused by the photodetector,which is a square-law detector. Nonetheless, electrical domain equalization schemes are widely used thanks to the multiple advantages, including compactness, flexibility, and low cost, that are brought by high-speed integrated-circuits technology.
Until now, most work on electronic domain equalizers has been focused on experimental studies [2 , 5 , 6] . The equalizer coefficients are determined by adaptive algorithms such as least mean square (LMS) algorithms and recursive least square (RLS) algorithms [7] . The obtained minimum mean square error (MMSE) solutions for optical on-off-keyed(OOK) systems have not been confirmed by theoretical studies.In this paper, we present theoretical closed form expressions for the MMSE linear equalizer coefficients and their MMSE performance. Through simulation, we verified that the calculated equalizer coefficients coincide with estimated equalizer coefficients using adaptive algorithm. Furthermore, once the MMSE equalizer coefficients are obtained, it is possible to estimate MSE performance of the MMSE linear equalizers in square-law nonlinear channels.
The rest of the paper is organized as follows. The next section describes the system model and computation of MMSE linear equalizer in square-law nonlinear channels.Simulation results using optical communication system are discussed in Section III. Section IV concludes the paper.
Lager Image
Fiber communication system using NRZ-OOK.
Consider a typical non-return-to-zero on-off-keyed (NRZ-OOK)optical communication model in Figure 1 .
The bit sequence { bk } with a bit rate T is pulse shaped by a continuous time pulse p(t) (for example, a raised cosine pulse with roll-off factor 1),
Lager Image
and transmitted over an optical fiber after laser modulation.We assume the fiber is a lossless linear channel with CD and first order PMD distortion. The CD and PMD are generally non-linear distortions, but often modeled as a linear distortion by a first order approximation in the optical field domain.The effect of dispersion can be modeled as a linear filtering process given by
Lager Image
where So(t) , ro(t) , and ho(t) are the transmitted, received signals, and the fiber impulse response in optical field domain,respectively, and ⨂ denotes convolution.
For example, the first-order PMD including CD distortion can be characterized as the following frequency response [2]
Lager Image
where L denotes the fiber length, λ the wavelength, D the dispersion parameter at λ, and c speed of light, γ the power splitting ratio, telling how the power of the input light of the fiber is divided onto the two input principle state of polarization, and τ denotes the differential group delay between the two polarizations.
The received optical signals having undergone these dispersions are demodulated to electronic signal by a squarelaw detector. In contrast to the signal detection process in wireless communication systems, this process is nonlinear.Hence, the linear convolution relation in optical field domain is not preserved in the baseband electronic domain.The received signal in electronic domain is
Lager Image
Electronic domain equalization model of opticalOOK systems.
Lager Image
where wo(t) denotes noise in the optical domain due to fiber amplifier and we(t) denotes the noise in the electronic domain. In the discrete time domain, by denoting rk:=re(kT) with ideal sampling timing phase [8] ,
Lager Image
where co(t) denotes the combined impulse response with the pulse shaping filter p(t) and PMD/CD ho(t) . We assume that the amplified-spontaneous-emission (ASE) noise from the erbium-doped fiber amplifiers dominates and the electronic domain noise can be neglected, i.e we(t) =0. Then, the whole process can be modeled as a linear channel followed by a square law detector as illustrated in Figure 2 .
We assume that the transmitted signal xk consists of {0, σ x }to model OOK (σ x =1 for the conventional OOK), and discrete time channel model c = [ c0 ,… cNc-1]T represent the optical domain channel, where [·] T denotes transpose operator.Furthermore, wk denotes the zero-mean optical noise with variance σ2n and yk denotes the output of square law detector. It is assumed that the transmitted signal xk and noise wk are statistically independent. In the next section we will derive MMSE equalizer applied to yk under this channel mode.
In contrast to the case in conventional wireless communication systems, the source signal xk in our channel model has nonzero mean and the following properties which make the derivation of MMSE equalizer difficult:
Lager Image
where E[·] denotes the expectation.
Denoting xk = [ xk ,… xk-(Nc-1) ] T a vector consisted of the source signals, the optical received signal before square law detector is expressed as
Lager Image
And the input to the linear equalizer is given by
Lager Image
Let f = [ f0, …fNt-1 ] T be the linear equalizer of length Nf .The MMSE equalizer coefficients of f is obtained by minimizing MMSE cost function
Lager Image
where yk = [ yk, …, yk-(Nk-1) ] T and Δ is decision delay. It is well known that the MMSE solution of f is given by the following [7]
Lager Image
The goal of this paper is to express E[ yk ykT ] and E[ yΔ yk ]in terms of the channel c = [ c0,…, cNc-1 ] T , the signal magnitude σx , and the noise variance σ2x . In order to obtain E[ yk ykT ] -1 , we need to calculate E[ yk ykT ] as shown below.
Lager Image
We know that E[ ykyk ]= E[ yk-δyk-δ ] = E[ yk-δyk ]. Because E[ ykykT ] is symmetric Toeplitz matrix [7] , we have only to compute one row or column of E[ ykykT ]. Note that E[ ykykT ] -1 is also a symmetric matrix. The elements of E[ ykykT ], E[ ykym ], are very complicated to compute since the transmitted signal xk has nonzero mean and the equalizer input yk is the magnitude square of received signal rk . Hence, it is impossible to provide a simple linear solution in contrast to the linear channel case. E[ ykym ] is composed of the following four terms:
Lager Image
Let’s define the first term of E[yk ym],
Lager Image
as a function α ( k, m, c, σ2x, σ2n ) :
Lager Image
When k = m
Lager Image
When ' k ?m '=δ (≠ 0)
Lager Image
The second term is a product of two expectations
Lager Image
where 1 =[1,…,1] T is the vector of length N c consisting of 1’s. Note that
Lager Image
The third term, the noise term, can be expressed as the square of noise variance.
Lager Image
The last term becomes the second term, when k = m .
Lager Image
Now, we have obtained the element of E[ ykykT ],
Lager Image
Consequently, we have obtained E[ ykykT ] in terms of the channel c = [ c0,…,cnc-1 ] T , the signal magnitude σx and the noise variance σ2x as desired and consequently E[ ykykT ] -1 can be computed by matrix inversion. Finally, we can obtain the MMSE equalizer tap f by computing E[ yΔ yk ]as shown below.
Lager Image
Once equalizer tap coefficients are obtained, the MSE performance of the linear equalizer can be calculated by substituting equalizer tap f into the MMSE cost function.
Lager Image
Once optical channel, optical noise variance, and the signal magnitude are given one can determine the MMSE equalizer and its MMSE performance. Numerical program languages such as Matlab are quite helpful to utilize the above results.
Lager Image
Comparison of equalizer taps obtained by calculationmethod and adaptive algorithm (L=100 km SNR=20 dB).
We consider mainly CD distortion to demonstrate the analysis on the electronic MMSE linear equalizers for optical OOK. In the frequency domain, the optical channel is given as
Lager Image
where the fiber CD parameter D is set to 17ps/nm/km, L is the fiber length, c is the speed of light, and λ is the signal wavelength, which is assumed to be 1550nm. By taking the inverse Fourier transform, the time domain transfer function can be calculated and expressed as [9]
Lager Image
In order to produce a discrete time channel impulse response, continuous time transfer function is sampled at a rate of 1/ T where T is the symbol period.
Fig. 3 ~ 5 compare the MMSE equalizer coefficients of Nf = 21 computed from analytical method and computed by adaptive algorithm for various degree of CD controlled by fiber length. The x- axis denotes the equalizer coefficient index k =1,…,21 and the y- axis draws the corresponding coefficient value f k . The agreement of two equalizer coefficients indicates the validity of the previous analysis that the computed linear equalizer is MMSE equalizer in squarelaw nonlinear channels.
Fig. 6 shows that the MSE performances of two equalizers also agree. This result suggests that the MMSE performance of a linear equalizer for CD can be analytically computed instead of adaptive simulation.
Lager Image
Comparison of equalizer taps obtained by calculationmethod and adaptive algorithm (L=150 km SNR=20 dB).
Lager Image
Comparison of equalizer taps obtained by calculationmethod and adaptive algorithm (L=200 km SNR=20 dB).
Lager Image
MSE performance of MMSE equalizer and adaptiveequalizer (SNR=20 dB).
We have performed analysis on the MMSE linear equalizer in square-law nonlinear channels found in optical OOK system and presented the MSE performance of equalizer. It has been verified that the analytically driven equalizer is indeed the MMSE solution by comparing with simulated equalizers by adaptive algorithms. As a result,we are able to theoretically bound MMSE performance of the electronic domain linear equalizer applied to optical OOK systems.
Agrawal G. P 2001 Nonlinear Fiber Optics Academic Press San Diego USA
Wang J , Kahn J. M 2004 Performance of electrical equalizers in optically amplified OOK and DPSK systems IEEE Photon. Technol. Lett. 16 1397 - 1399    DOI : 10.1109/LPT.2004.826222
Gnauck A. H , Cimini L. J , Stone J , Stulz L. W 1990 Optical equalization of fiber chromatic dispersion in a 5-Gb/s transmission system IEEE Photon. Technol. Lett. 8 585 - 587
2005 Performance of single-mode fiber links using electronic feed-forward and decision feedback equalizers IEEE Photon. Technol. Lett. 17 2206 - 2208    DOI : 10.1109/LPT.2005.856326
Weiss A. J 2003 On the performance of electrical equalization in optical fiber transmission systems IEEE Photon.Technol. Lett. 15 1225 - 1227    DOI : 10.1109/LPT.2003.816120
Ibragimov E 2006 Limits of optical dispersion compensation using linear electrical equalizers IEEE Photon. Technol.Lett. 18 1427 - 1429    DOI : 10.1109/LPT.2006.877225
Proakis J. G 2001 Digital Communication McGraw-Hill New York USA
Kim K. S , Lee J , Chung W , Kim S. C 2008 An electronic domain chromatic dispersion monitoring scheme insensitive to OSNR using kurtosis J. Opt. Soc. Korea 12 249 - 254    DOI : 10.3807/JOSK.2008.12.4.249
Khafaji M , Gustat H , Ellinger F , Scheytt C 2010 General time domain representation of chromatic dispersion in singlemode fibers IEEE Photon. Technol. Lett. 22 314 - 316    DOI : 10.1109/LPT.2009.2038355