ADOC ch04 pp135-154
ADOC ch04 pp135-154
4 Transmission Performance
Fundamentals
Signals and noise are important parameters to determine the quality of any engineering system,
especially in the telecommunications field. The modulated signals under different modulation for-
mats will be treated in subsequent chapters. This chapter presents an overview of the noise process
and mechanism in terms of equivalent noise power and contribution to the optoelectronic process at
the optical receiver. Further, the integration of the noise and the statistical process for determination
of the receiver sensitivity under transmission of different modulation formats are described together
with modeling techniques in the MATLAB and Simulink platform. Statistical analyses are also
given for the evaluation of the transmission performance, especially when there are no uniform
probability distributions.
4.1 INTRODUCTION
Optical detection and the noise interference in such processes are critical to the performance of
optical communications systems. A fundamental understanding of the noise and sensitivity of the
receiver, the minimum optical power available at the photodetector, requires evaluation of the trans-
mission performance of different modulation formats. Optical receivers have evolved from binary
digital direct detection to coherent mixing of signals and local lightwave oscillator and detection
with optical preamplification. This chapter presents an overview of the noise process and mecha-
nism in terms of the equivalent noise power and contribution.
A schematic diagram of a single-channel DPSK system is illustrated in Figure4.1.
As the bit or symbol rate of the optical transmission systems is increased, the demand for mod-
eling is intense, especially for a modeling platform that can simulate the real physical photonic
subsystems. Inorder to enhance the effective transmission capacity with minimum renovation of
the photonic and electronic subsystems, there are two possible solutions: (1) increasing the base
bit rate or (2) employing a multilevel modulation technique such as the M-ary amplitude and/or
phase shift keying [1]. Thelatter solution would be preferred because for the same as the bit rate,
the multilevel offers a much higher transmission capacity (in bits/sec.) but with the same baud rate
without increasing the complexity of the transmitters and receivers as well as the dispersion toler-
ance of the transmission fibers which can be single-mode fiber types ITU G652 or G655. However
due to the nonlinear threshold of the fiber the highest level is limited and so an increase in the level
of the multilevel modulation system will lower the opening of the eye diagram of sublevels. Thus
this may reduce bit error rate.
One possible technique that could offer insights into possible solutions that facilitate the simula-
tion of ultra-high-capacity and ultra-high-bit-rate transmission systems is to develop a comprehensive
modeling platform. Further, the modeling platform should take advantage of any user-friendly soft-
ware platform that is popular, easy to use, and easy to extend. This platform would offer the research
community of optical communication engineering a basis for extension and enhance the linkages
between research groups.
135
136 Advanced Digital Optical Communications
Thus, one of the principal objectives of this chapter is to present the development of MATLAB
and Simulink 7.0 [2] for computer experiments in optical transmission systems under advanced
modulation formats, especially the amplitude and/or phase shift keying modulation of binary or
multilevel. Indeed, duobinary in association with an appropriate low-pass filter can offer multiple
rates much higher than the basic rate. The advantage of duobinary is that the detection is much
simpler than that of DPSK, DQPSK, or M-ary PSK, in which the receiver is a direct detection type
receiver. In the photonic domain, the combination of phase and amplitude is taken into account for
coding modulation schemes that involve the tri- or higher-order levels. These features can be imple-
mented with ease on the Simulink platform.
We base our system bit rate on a 40 Gbps per channel. The modulation formats of binary RZ,
NRZ, and CSRZ-ASK-DPSK are demonstrated with the modeling of the transmission over stan-
dard single-mode optical fiber (SSMF), as well as the photonic decoding of this optical signaling.
Thetransmission performance of some optical systems is also given in terms of the bit error rate
(BER), the -factor (usually known as Q factor), and the receiver sensitivity.
In this chapter, we present, to the best of our knowledge, the first MATLAB and Simulink
platform simulation test bed for modeling of the transmission of the amplitude and differential
phase modulation incorporating RZ, NRZ, and carrier-suppressed RZ formats for advanced high-
capacity and long-haul optical fiber transmission systems. A novel, modified fiber propagation algo-
rithm has been developed and optimized to minimize the simulation processing time and enhance
its accuracy. The performance of optical transmission systems can be automatically and accurately
evaluated with various methods. Analysis and methodologies for future development of the simula-
tion test bed are also presented.
In recent years, the capacity and reach of optical transmission systems have dramatically
increased owing to the accelerating growth of data usage demand (Internet, peer-to-peer
network,).Therefore, the necessity of upgrading of current dense wavelength division multiplex-
ing (DWDM) 10 Gbps systems to DWDM 40Gbps or even higher bit rates becomes crucial to tele-
communications service providers. Accordingly, several modulation formats have been proposed
and investigated as the alternatives to the current onoff keying (OOK) intensity modulation, which
is severely degraded at high bit rates owing to dispersion and nonlinear effects of the transmis-
sion fibers. Apart from the requirements of robustness to the transmission impairments, the cost-
effectiveness of system upgrades is also significant. Among the candidates, differential binary and/
or quadrature phase shift keying (DPSK/DQPSK) have recently attracted much attention owing to
the following advantages: (1) a 3 dB improvement on receiver sensitivity (if the balanced receiving
technique is used) [1,2]; (2) high tolerance to fiber nonlinearities, especially to intrachannel nonlin-
ear effects, cross-phase modulation (XPM), and four-wave mixing (FWM) [2,3]; (3) superior spec-
tral efficiency (DQPSK), and hence high tolerance to optical filtering [3]; and finally, (4)advantages
in all-optical networks incorporating optical add-drop multiplexers or optical cross-connects [4].
Recently, several experimental demonstrations of DPSK/DQPSK long-haul transmission DWDM
Optical Receivers and Transmission Performance 137
systems for 10 and 40 Gbps have been reported [5,6]. Therefore, a simulation test bed is necessary for
detailed design, investigation, and verification of the benefits and shortcomings of these advanced
modulation formats on fiber-optic transmission systems. In this chapter, we present, to the best of
our knowledge, the first MATLAB-and Simulink-based simulation package for photonic DWDM
systems. The simulator is still in the first phase of the development, with the focus on single-channel
systems, and it is being continuously improved and updated.
This chapter is organized as follows: (i) in Section4.2, receiver noise sources are outlined; noise
analyses and performance of basic binary digital optical receivers are then given in Section4.3;
(ii)Section4.3 presents the architecture and operational principles of the simulator; (iii) Section4.4
presents simulation results of the up-to-date simulator, and the corroboration of the developed
models in comparison with experimental results are presented in Section 4.5 for both ASK and
DPSK modulation formats. Finally, concluding remarks are presented.
2
iNshot = 2qPin Be (4.1)
2
iND 2qI D Be (4.2)
2
iEA ( I eq )2 Be (4.3)
d (iI2 )
sI 2qI in A2 /Hz (4.4)
df
where q is the electronic charge. The current iI represents the noise current generated due to the
biasing current I.
d is2
sQ 2q is2 (4.5)
df
This is the signal-dependent noise, and it is a unique feature of optical communications. When
the avalanche photo-detector (APD) is used, the noise spectral density is given by
d is2
sQ 2q is2 Gn2 (4.6)
df
It is noted here that the dark currents are generated in the photodetector is a natural phenom-
enon in semiconductor material operating at room temperature. This dark current noise must be
included in the total equivalent noise current as referred to the input. These currents are generated
even in the absence of an optical signal. These dark currents can be eliminated by cooling the
photodetector to at least below the temperature of liquid nitrogen (77K).
sR =
d iR2( ) = 2k T B
(4.7)
df R
where k B is the Boltzmanns constant, T is the absolute temperature (in K), R is the resistance in
ohms, and iR denotes the noise current due to R.
where Gop = amplifier gain; nsp = spontaneous emission factor; m = number of polarization
modes (1 or 2); Bo = mean noise in bandwidth; optical signal-to-noise ratio (OSNR) at the output
ofEDFA.
Optical Receivers and Transmission Performance 139
OSNR in
FN = 2nsp for Gop 1 (4.9)
OSNR out
Modern optical amplifiers are optimized, and nsp reaches unity; so the NF can reach 3 dB.
2 Be
isig-ASE = 2(qGop )2 (2nsp ) Pin (4.10)
h
where nsp,1 is the noise as seen from the input of the amplifier, and nsp,1 is the equivalent noise factor
at the output of the amplifier. Then, the accumulated ASE noise over the transmission spans is given
140 Advanced Digital Optical Communications
by the sum of all the noise sources of the amplifiers. At the output of the final Nth stage, without
taking into account the other noise sources, the equivalent ASE accumulated noise is given as
=
nsp,N ni =1
sp,i
(4.12)
e GN
nsp,N
nsp,i = + nsp,i
L
under the assumption that the gain of optical amplifiers has equalized all the losses of transmission
and dispersion-compensating fibers (DCFs).
When an optical preamplified receiver and two optical filters are placed before and after the
optical preamplifier, the effective ASE noise is given by
nsp,NGN 2 (4.13)
nspe = + nsp
L 1
where 1 and 2 are the time constants of the two optical filters.
The NF of the equivalent accumulated noise is then given by
This is also the amount of degradation of the OSNR due to the accumulated ASE noise at the
front of the photodetector.
p(t ) = a h (t jT )
j
j p b (4.15)
The average output voltage at the output of the electronic preamplifier is thus given by
Tb / 2
1
vo = Gn j
aj
Tb hp (t jTb ) dt RI A (4.16)
Tb / 2
Optical Receivers and Transmission Performance 141
where Tb is the bit period, and hp(tjTb) is the impulse response of the system evaluated at each time
interval. RI is the input resistance of the overall amplifier of the system including both the front end
and the linear channel amplifier. It is assumed that the overall amplifier has a flat gain response
Aover the bandwidth of the system.
bo 0
aj t=0 = (4.17)
b1
with b 0 being the energy when a transmitted 0 is received, and b1 is the energy when a transmitted 1
is received. The sum over a number of periods is necessary to take into account the contribution of
adjacent optical pulses. We now have to distinguish between two cases when a 0 or a 1 is transmit-
ted and received.
bo
vo o
= voo = GRI A 0 (4.18)
Tb
2
with the total equivalent noise voltage at the output, vNTo is
2
vNTo 2
vNA 2
iNeq RI2 A2 (4.19)
where iNeq is the total equivalent current at the input of the electronic preamplifier. Appendix 2
(Section 4.12) and Annex 4 gives a method for estimating this noise current for any preamplifier
whose equivalent Y-parameters are known.
b1
vo 1
= v01 = Gn RI A (4.20)
Tb
with the total noise equivalent mean voltage at the output given by
2 2 2
vNT1 = voSN + vNA (4.21)
where v2osN is the signal-dependent shot noise. v2NA is the amplifier noise at the output and is given by
2 2
vNA iNeq R12 A2 B (4.22)
The signal-dependent noise is, in fact, the quantum shot noise and is given by
B
2
voSN =
o
2q is 1
Gn2 RI2 A2 df (4.23)
142 Advanced Digital Optical Communications
where B is the 3 dB bandwidth of the overall amplifier, <is>1 is the average photocurrent received
when a 1 was transmitted. This current can be estimated as follows
Tb / 2
is 1
= Tb
1
b Tb / 2
hp (t jTb ) dt (4.24)
or
b1
is 1
=
Tb
hp (t ) dt (4.25)
with a normalization that is
Equation 4.21 becomes
hp (t ) dt = 1.
b1
is 1
= (4.26)
Tb
1 ( v0 v00 )2
p [ v0 0 ] = 2 1/2 exp 2 (4.27)
(2v NTO ) 2vNT0
1 ( v0 v01 ) 2
p [ v0 1 ] = 1/2 exp 2 (4.28)
(2 vNT1 ) 2v NT1
where p(1) and p(0) are the probabilities of receiving a 1 and a 0, respectively, and p(1/0) andp(0/1)
are the probabilities of deciding 1 when a 0 is transmitted and vice versa. In an OOK bit stream,
1and 0 are likely to occur equally, that is, p(1) = p(0) = 0.5; then Equation 4.25 becomes
1
BER = p ( 0/1) + p (1/0 ) (4.30)
2
Optical Receivers and Transmission Performance 143
Thus, as an equal probability of transmitting a 0 and a 1 is assumed, for a decision voltage level
of d, the total probability of error PE is the sum of the errors of deciding 0 or 1. Integrating the pdf
over the overlapping regions of the two pdfs gives
d
1 1
p [ v0 0 ] dv0 +
PE =
2
d
2 p[ v
0 1 ] dv0 (4.31)
Substituting for the probability distribution using Equations 4.26 and 4.27 leads to
d
x2 x2
1 1
BER =
2 e
d v 00
2 dx +
2 e
v 01 d
2 dx (4.32)
vNT 0 v NT1
The functions in Equation 4.28 have the standard form of the complementary error function Q()
defined as
1 2
Q ( a) =
2
e x /2
dx (4.33)
1 d v d
BER = Q + Q 01 (4.34)
2 vNT0 vNT1
d v d
= 01 o (4.35)
vNT0 vNT1
The Marcum Q() function is a standard function, and this curve is shown in Figure4.2. Note
that for a BER = 10 9, the value of is about 6, which is the normal standard for communications
at bit rates of 134 Mb/s40 Gbps.
Thus, by eliminating the decision-level variable d from Equation 4.31, we obtain
v01 vNT0
= (4.37)
v NT1
144 Advanced Digital Optical Communications
100
102
104
Probability error
106
108
1010
1012
0 1 2 3 4 5 6 7
or alternatively
1/2
b1 b1 2
Gn R1 A = 2q Gn BR12 A2 + v 2NA + vNA (4.39)
Tb Tb
However, the amplifier noise voltage vNA is given by vNA = iNeqRI AB1/2. Thus, by substituting
this noise voltage, eliminating RI A, and solving for the energy required for the 1 transmitted and
received at the photodetector b1, we obtain
2iNeqTb
b1 = q 2G + (4.40)
qG
where we have used the approximation <Gn> = G and <Gn2> = G 2+x with x being the factor depen-
dent on the ionization ratio in an APD. For a pin photodetector, G = 1. The optical receiver sensitiv-
ity in dBm, denoted as RS in dBm, can thus be obtained as
pav
RS 10Log10 in dBm (4.41)
p0
where Pav = b1/Tb and Po = 1.0 mW. This is the optical receiver sensitivity defined as the mini-
mum optical power required for the receiver to operate reliably with a BER below a specific value.
InEquation 4.36, there are two terms clearly specifying the dependence of the signal dependence
Optical Receivers and Transmission Performance 145
and the amplifier noise contribution at a certain data rate and for a certain modulation format.
Recall that the term b1 represents the optical energy of the 1 required for the optical receiver to
detect with a certain BER. The term q/R in the right-hand side of Equation 4.36 is equivalent to the
optical power required to generate one electron, or the number of photons required to generate one
electron. A typical measure for the optical receiver is the number of photons required for it to oper-
ate with a specified BER.
Note that the common term of the energy b1 is q/R, which is equivalent to about one photon
energy; we can rewrite Equation 4.36 by using q/R = h/
h 2 2iNeq Tb
b1 = G+ (4.42)
qG
Thus, we can observe from equation 4.42 when we have an ideal electronic amplifier, that is, a
noiseless amplifier iNeq = 0, then it requires only a number of photon energy of 2 G for error-free
detection. The second term in the bracket is thus the number of photons required to overcome the
amplifier noise. Note that the determination of the BER using the eye diagram when the pdf is not
Gaussian is inaccurate. However, a much more accurate statistical analysis is developed and pre-
sented later in Section4.9.
For an average output power P0 of the EDFA, the OSNR can be determined by
P0
OSNR = (4.44)
2nsp hfT (G 1) B0
where nsp is the number of spontaneous emissions of the EDFA, h is Planks constant, and G is the
linear power gain coefficient of the EDFA. Then, if the signal ASE beat noise dominates the noises
of the detection process, an NF of the EDFA can be used to determine the ratio of the OSNR at the
input and output as this is much more practical. Thus, the NF can be approximated by Fn ~ 2nsp.
Thus, the OSNR can be rewritten as
P0
OSNR = (4.45)
Fn hfT (G 1) B0
the optical gain is much greater than unity, the OSNR can be written for an operating wavelength
of 1550 nm as
If an eye diagram is used, then the eye opening penalty (EOP) can be determined for the same
BER. The EOP can be written as
EO
EOP = !10Log10 .....dB (4.47)
EO n
after thephotodetection and at the input of the electronic amplifier. Note that coherent systems are
identified with the mixing of the optical signals and local oscillator whose polarization directions
are aligned with each other.
Electronic
preamplifier
S1,0(t)
(a)
(b)
Electronic
preamplifier
S1,0(t)
OA
(c)
(d)
FIGURE4.3 Schematic diagram of (a) direct detection, (b) coherent detection, (c) direct detection with OA,
and (d) coherent detection with OA. OA=optical amplifier, FC = fiber coupler, PD = photodetector, DSP =
digital signal processor, ADC = analog-to-digital converter. (Continued)
148 Advanced Digital Optical Communications
Modulated lightwave
channels
Decision
Local
oscillator Elect. receiver
(e)
Lightwave channels
Sampler
DSP
ADC
Local
Electronic
oscillator
preamp.
(f )
Decision
Lightwave
channels
PD DSP
FC Electronic
preamp.
Fiber one-bit delay
(g)
FIGURE4.3 (Continued) Schematic diagram of (e) fiber version of coherent receiver, (f) coherent receiver
using balance detection with two photodetectors connected back to back and a local oscillator (laser),
(g) self-homodyne reception balanced receiver with one-bit delay photonic phase comparator. OA=optical
amplifier, FC = fiber coupler, PD = photodetector, DSP = digital signal processor, ADC = analog-to-digital
converter.
For the ASK system, the signals for the 1 and 0 can be expressed in terms of the amount of photon
energy over the entire bit period, contained within a bit period T as
np
for 1
s1,0 (t ) = T 0 t T (4.48)
0 for 0
where np is the number of photons, and the energy of the lightwave is normalized with a single
photon energy at the operating wavelength. Thus, the total optical energy is
Es = nph.T (4.49)
Thus, the electronic current generated after the photodetector, with R being the responsivity of
the photodetection, is
q
is = np h.T = np h.T (4.50)
h
Optical Receivers and Transmission Performance 149
Thus, we could say that np is the number of photons per bit required for the detection of a1 if
there is no noise contributed by the electronic amplifier or detection.
With the probability of 1 and 0 being equal (50%), the probability of error of the detection is
1 np
Pe = e (4.51)
2
Thus, for a BER of 1e9 the argument np = (32)2 or np = 18 with an allowance of a factor of
thesingle-sided estimation so that np = 20 for the full detection error. This is the super quantum
limit. We also assume a unity responsibility of the photodetection.
2 q 1 1 (4.52)
iN(LO) = 2qPLO = 2q 2 PLO
h T h T
and the SNR is given by
is2 4 2 ps (t ) PLO /T
SNR = 2
= = 2np (4.53)
iN(LO) 2qPLO / T
Pe erfc(2np ) (4.54)
where d is the signal power separation between the average level of the 1 and 0 for binary sys-
tems and the equidistance between the constellation points of the modulation scheme as shown in
Figure4.4. Let E1 and E0 be the field amplitudes of the signals 1 and 0; then, the Euclidean distance
d is given by
d 2 = E12 + E02 2E1E0
(4.56)
T
2
with =
T s (t ) s (t ) dt
0
1 0
is the correlation coefficient between the two logic levels or, alternatively, the Euclidean angle
between the two vector signals as represented on the scattering plane of the constellation.
150 Advanced Digital Optical Communications
d
0
FIGURE4.4 Signal constellation and energy level and the geometrical distance between 1 and 0.
d
0 1 and E1
FIGURE 4.5 Signal constellation, energy level, and the geometrical distance between 1 and 0 of ASK
system.
2np
cos IF t for 1
s1,0 (t ) = T 0t T (4.57)
0 for 0
where the amplitude of the lightwave-modulated signal is expressed in terms of the energy of the
photons over the time interval, and thus the square root of this amount is the amplitude of the field
Optical Receivers and Transmission Performance 151
of the lightwave. Naturally, the characteristic impedance of medium is set at unity. Thedistance is
then d = np.
Thus, the BER is given by
np
BER = erfc (4.58)
2N0
2np
cos 1t for 1
T
s1,0 (t ) = 0t T (4.59)
2np
cos 0 t for 0
T
The FSK modulation scheme with two distinct frequencies f1 and f 2 can be represented with a
constant envelope and variation in the carrier frequency, or continuous phase between the two states
as shown in Figure4.6.
The modulation index can be defined as
1 0
m= T (4.60)
2
0 1 and E1
FIGURE4.6 Signal constellation and energy level evolution of the signal envelope and the continuous phase
between the 1 and 0 of FSK system.
152 Advanced Digital Optical Communications
If the two frequencies are large enough and the second and higher harmonics are outside the
detection region, the signal correlation coefficient is given by
T T
2 2
=
T
0
s1 ( t ) s0 ( t ) dt =
T cos ( t ) cos ( t ) dt
1 0
(4.61)
sin(2 m)
2 m
1 n sin(2 m)
BER = erfc p 1 (4.62)
2 2 2 m
np sin(2m)
BER = 10 9 1 =3 2
2 2m
(4.63)
36
then, np =
sin(2m)
1
2m
For MSK, the modulation index is 0.25, leading to the required number of photon energy per
bit of 6070, much higher than that of 0.8 at which only 30 photons energy is required per bit. Itis,
however, shown that the MSK can be optimum for the transmission over a dispersive medium owing
to the optimum bandwidth of the modulation scheme and hence minimum dispersive effects on
thephase of the carrier.
The FSK can be implemented with continuous phase frequency shift keying (CPFSK); that is,the
phase of the carrier is continuously chirped. Assume that the phase is linearly chirped such that
thephase variation for the 1 and 0 are given by
2np
cos 1 (t ) for 1
T
s1,0 (t ) =
2np
cos 0 (t ) for 0
T
where
m
IF t + t for 1 (4.64)
T T
1,0 (t ) = 0t
m 2m
IF t t for 0
T
IF t + for 1
2 T
1,0 (t ) = tT
2m
IF t for 0
2
Optical Receivers and Transmission Performance 153
1 1 1
BER = 10 9 = erfc np 1 np 1
=3 2
2 4 m 4 m
(4.65)
18
then, np =
1
1
4m
The variation in the photon energy with the modulation index for this linear CPFSK is shown in
Figure4.7.
The FSK can be modified with the control of the relative phase of the two carriers between the
two bits; thus, the data bits can be written as
2np
cos(1t + ) for 1
T
s1,0 (t ) = 0t T (4.66)
2np
cos 0t for 0
T
m for p m 2 p + 1 p = 0,1,2....
1,0 (t ) = (4.67)
m for 2 p + 1 m 2( p + 1)
50
45
40
No. of photon energy
35
30
25
20
0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
Modulation index (m)
FIGURE4.7 Fundamental limits of the coherent phase and frequency shift keying detection for (a) mid-gray
o CPFSK, (b) light gray + FSK (1 and 2) with optimized correlation coefficient with phase control, and
(c) black * linear DMPSK.
154 Advanced Digital Optical Communications
T T
2 2
=
T
0
s1 (t )s0 (t )dt =
T cos( t + ) cos( t)dt
0
1 0
(4.68)
sin(m)
m
1 n sin ( m )
BER = erfc p 1 +
2 2 m
(4.69)
36
np = for a....BER = 10 9
sin( m)
1+
m
The required number of photon is plotted against the modulation index and shown in Figure4.7b
black*.
Electrical Data
filter recovery
Tb
one-bit delay