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

- Received : July 21, 2011
- Accepted : August 02, 2011
- Published : September 25, 2011

Download

PDF

e-PUB

PubReader

PPT

Export by style

Share

Article

Metrics

Cited by

TagCloud

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.
Fiber communication system using NRZ-OOK.
b_{k}
} 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),
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
where
S_{o}(t)
,
r_{o}(t)
, and
h_{o}(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]
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
Electronic domain equalization model of opticalOOK systems.
where
w_{o}(t)
denotes noise in the optical domain due to fiber amplifier and
w_{e}(t)
denotes the noise in the electronic domain. In the discrete time domain, by denoting
r_{k}:=r_{e}(kT)
with ideal sampling timing phase
[8]
,
where
c_{o}(t)
denotes the combined impulse response with the pulse shaping filter
p(t)
and PMD/CD
h_{o}(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
w_{e}(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
x_{k}
consists of {0, σ
_{x}
}to model OOK (σ
_{x}
=1 for the conventional OOK), and discrete time channel model c = [
c_{0}
,…
c_{Nc-1}]^{T}
represent the optical domain channel, where [·]
^{T}
denotes transpose operator.Furthermore,
w_{k}
denotes the zero-mean optical noise with variance
σ^{2}_{n}
and
y_{k}
denotes the output of square law detector. It is assumed that the transmitted signal
x_{k}
and noise
w_{k}
are statistically independent. In the next section we will derive MMSE equalizer applied to
y_{k}
under this channel mode.
x_{k}
in our channel model has nonzero mean and the following properties which make the derivation of MMSE equalizer difficult:
where E[·] denotes the expectation.
Denoting xk = [
x_{k}
,…
x_{k-(Nc-1)}
]
^{T}
a vector consisted of the source signals, the optical received signal before square law detector is expressed as
And the input to the linear equalizer is given by
Let
f
= [
f_{0}, …f_{Nt-1}
]
^{T}
be the linear equalizer of length
N_{f}
.The MMSE equalizer coefficients of
f
is obtained by minimizing MMSE cost function
where
y_{k}
= [
y_{k}, …, y_{k-(Nk-1)}
]
^{T}
and Δ is decision delay. It is well known that the MMSE solution of
f
is given by the following
[7]
The goal of this paper is to express E[
y_{k} y_{k}^{T}
] and E[
y_{Δ} y_{k}
]in terms of the channel
c
= [
c_{0},…, c_{Nc-1}
]
^{T}
, the signal magnitude
σ_{x}
, and the noise variance
σ^{2}_{x}
. In order to obtain E[
y_{k} y_{k}^{T}
]
^{-1}
, we need to calculate E[
y_{k} y_{k}^{T}
] as shown below.
We know that E[
y_{k}y_{k}
]= E[
y_{k-δ}y_{k-δ}
] = E[
y_{k-δ}y_{k}
]. Because E[
y_{k}y_{k}^{T}
] is symmetric Toeplitz matrix
[7]
, we have only to compute one row or column of E[
y_{k}y_{k}^{T}
]. Note that E[
y_{k}y_{k}^{T}
]
^{-1}
is also a symmetric matrix. The elements of E[
y_{k}y_{k}^{T}
], E[
y_{k}y_{m}
], are very complicated to compute since the transmitted signal
x_{k}
has nonzero mean and the equalizer input
y_{k}
is the magnitude square of received signal
r_{k}
. Hence, it is impossible to provide a simple linear solution in contrast to the linear channel case. E[
y_{k}y_{m}
] is composed of the following four terms:
Let’s define the first term of E[yk ym],
as a function
α
(
k, m, c, σ^{2}_{x}, σ^{2}_{n}
) :
When
k = m
When '
k ?m
'=δ (≠ 0)
The second term is a product of two expectations
where
1
=[1,…,1]
^{T}
is the vector of length N
_{c}
consisting of 1’s. Note that
The third term, the noise term, can be expressed as the square of noise variance.
The last term becomes the second term, when
k = m
.
Now, we have obtained the element of E[
y_{k}y_{k}^{T}
],
Consequently, we have obtained E[
y_{k}y_{k}^{T}
] in terms of the channel
c
= [
c_{0},…,c_{nc-1}
]
^{T}
, the signal magnitude
σ_{x}
and the noise variance
σ^{2}_{x}
as desired and consequently E[
y_{k}y_{k}^{T}
]
^{-1}
can be computed by matrix inversion. Finally, we can obtain the MMSE equalizer tap f by computing E[
y_{Δ} y_{k}
]as shown below.
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.
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.
Comparison of equalizer taps obtained by calculationmethod and adaptive algorithm (L=100 km SNR=20 dB).
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]
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
N_{f}
= 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.
Comparison of equalizer taps obtained by calculationmethod and adaptive algorithm (L=150 km SNR=20 dB).
Comparison of equalizer taps obtained by calculationmethod and adaptive algorithm (L=200 km SNR=20 dB).
MSE performance of MMSE equalizer and adaptiveequalizer (SNR=20 dB).

Optical communications
;
Square-law detector
;
Electrical equalizers
;
Nonlinear channels
;
(060.2360) Fiber optics links and subsystems
;
(060.4510) Optical communications
;
(060.0060) Fiber optics and optical communications
;
(060.2330) Fiber optics communications

I. INTRODUCTION

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

II. SYSTEM MODEL

Consider a typical non-return-to-zero on-off-keyed (NRZ-OOK)optical communication model in
Figure 1
.
The bit sequence {
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

III. MMSE EQUALIZER

In contrast to the case in conventional wireless communication systems, the source signal
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

IV. SIMULATION RESULT

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

Lager Image

Lager Image

Lager Image

Lager Image

V. CONCLUSION

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

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

Proakis J. G
2001
Digital Communication
McGraw-Hill
New York USA

Citing 'Performance Analysis of Electrical MMSE Linear Equalizers in Optically Amplified OOK Systems
'

@article{ E1OSAB_2011_v15n3_232}
,title={Performance Analysis of Electrical MMSE Linear Equalizers in Optically Amplified OOK Systems}
,volume={3}
, url={http://dx.doi.org/10.3807/JOSK.2011.15.3.232}, DOI={10.3807/JOSK.2011.15.3.232}
, number= {3}
, journal={Journal of the Optical Society of Korea}
, publisher={Optical Society of Korea}
, author={Park, Jang-Woo
and
Chung, Won-Zoo}
, year={2011}
, month={Sep}