Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Semi-analytical approach for performance evaluation of direct-detection OFDM optical communication systems

Open Access Open Access

Abstract

A semi-analytical method to evaluate the bit error ratio (BER) in direct-detection (DD) optical fibre transmission systems employing orthogonal frequency division multiplexing (OFDM) and optically preamplified receivers is proposed. The method considers a Gaussian approach for the signal at the equalizer output and allows evaluating accurately the BER of each OFDM subcarrier for the receiver structure considered and for practical optical and electrical filters shapes, being a powerful tool to perform the optimization of these systems. The results obtained by the proposed method have shown excellent agreement with Monte Carlo estimates for DD-OFDM ultra-wideband radio signals and for two different DD-OFDM signals proposed for long-haul systems.

©2009 Optical Society of America

1. Introduction

Transmission of orthogonal frequency-division multiplexing (OFDM) signals in optical communication systems has received an enhanced interest in the past few years due to its robustness to intersymbol interference (ISI) [1]. Transmission of OFDM signals over optical fibre has been experimentally proposed for different systems such as radio-over-fibre for access networks [2, 3], metro networks [4] and long-haul systems [5]. In addition to experimental work, some important aspects of the transmission of OFDM signals over optical fibre have been also analyzed through numerical simulation [6, 7].

The bit error ratio (BER) is the figure of merit commonly used to assess the transmission performance of OFDM signals over optical fibre. Three approaches have been used to calculate the BER of direct-detection (DD) OFDM systems: (i) Monte Carlo (MC) simulation, where the signal and the noise introduced in the system are generated and propagated together along the transmission system and the BER is obtained by direct error counting [6, 8, 9]; (ii) semi-analytical approaches, where the signal is generated and propagated along the transmission system to assess the channel induced distortions on the signal waveform, whereas the noise and the BER are characterized by analytical means using the information of the noise statistical properties [10]; and (iii) analytical methods, where the signal and the noise at the system output are characterized analytically [11]. Though accurate BER estimation is achieved byMCsimulation, low BER values are remarkably difficult to simulate as they require unacceptable computation time. Instead, the assessment of the system performance by using analytical methods provide fast results independently of the BER range. However, this is achieved at the expense of the performance estimation accuracy as closed-form expressions are usually obtained by neglecting a few sources of signal distortion. Semi-analytical approaches provide a compromise between the other approaches. The semi-analytical method reported in [10] allows quickly evaluating the BER in each subcarrier from the signal-to-noise ratio. However, the bandwidth limitation and the distortion induced by the OFDM receiver on the noise and on the OFDM signal are neglected and it is not clear the way how the distortions induced on each OFDM subcarrier by the optical link are considered as the formulation used to evaluate the BER considers absence of ISI and intercarrier interference (ICI).

In this work, a semi-analytical Gaussian approach (SAGA) to evaluate the BER through numerical simulation in OFDM optical transmission systems for optically pre-amplified/DD receivers is proposed. Although the SAGA has been initially developed to evaluate the transmission performance of DD-OFDM ultra-wideband (UWB) radio signals over optical fiber, it is shown that it applies also to other DD-OFDM systems. The SAGA presents several novel contributions relative to the previous semi-analytical approach reported in [10]: (i) it allows obtaining the variance of the I and Q components of each subcarrier at OFDM receiver output; (ii) the dependence of the signal-noise and noise-noise beat terms contribution to the variance of each subcarrier on the optical and electrical filters shape and on the FFT block are considered; (iii) it takes into account eventual phase differences of the local oscillator in the I and Q branches of the electrical receiver; (iv) it considers the influence of the equalizer transfer function on the signal-noise and noise-noise beat terms; and (v) the statistical distribution of the ISI between consecutive OFDM symbols and of the ICI are correctly taken into account.

2. Theory

In this section, the optically pre-amplified DD receiver of the OFDM optical transmission system is presented. Afterwards, closed-form expressions for the mean and the variance of each subcarrier at system output are derived and the SAGA used to evaluate the system performance is presented.

2.1. System modeling

 figure: Fig. 1.

Fig. 1. Block diagram of the optically pre-amplified DD-OFDM receiver.

Download Full Size | PDF

Fig. 1 shows a scheme of the optically pre-amplified DD-OFDM receiver. The receiver includes an optical amplifier (OA), an optical filter (OF) and the PIN photodetector. The OF is characterized by its impulse response, ho(t), and transfer function, Ho(ν), and the PIN photodetector by a quadratic detection characteristic with responsivity Rλ. The electrical part of the receiver consists in the I and Q demodulator branches, the FFT block, the equalizer and the symbol demapping block. The phase of the local oscillator for the I and Q branches, identified in Fig. 1 as ϕ 1 and ϕ 2, respectively, can represent either the inaccuracy of the local oscillator or different propagation delays between the I and Q branches of the electrical receiver. In the FFT block, N represents the total number of subcarriers transmitted in each OFDM symbol.

2.2. Gaussian distribution characterization

To characterize analytically the signal as a Gaussian distribution, closed-form expressions for the mean and variance of the signal at the equalizer output are derived in this subsection.

Considering that the signal field at OF input is polarized and the noise field is completely unpolarized, the Jones vectors of the OFDM signal field, r s(t), and the amplified spontaneous emission (ASE) noise field, rn(t), at OF input can be written as r s(t)=r s,‖(t)u and rn(t)=r n,‖(t)u +r n,⊥(t)u, respectively, where r s,‖(t), r n,‖(t) and r n,⊥(t) are the OFDM signal and the ASE noise field components in the parallel (‖) and perpendicular (⊥) polarizations directions defined by u and u , respectively. The derivation of the SAGA assumes that the received signal can be demodulated by the DD-OFDM receiver shown in Fig. 1, independently of the OFDM transmitter employed. OFDM-UWB radio-over-fiber and single-sideband (SSB) OFDM signals used in access [2] and in long-haul DD systems [14, 15], respectively, are examples of different OFDM transmitters that can use the receiver of Fig. 1 to demodulate the received signal. The ASE noise field due to the OA is assumed zero mean additive white Gaussian noise with power spectral density (PSD) in each polarization direction (parallel and perpendicular to the signal) given by S ASE,(‖,⊥)= 0 nsp(g-1), where h is the Planck’s constant, nsp is the spontaneous emission noise factor, ν 0 is the optical carrier frequency and g is the OA gain. Considering that the OF operation is independent of the polarization direction, the field at the OF output can be expressed as a function of the signal and noise field components as e(t)=s(t)+n(t), where s(t)=r s(t) * ho(t), n(t)=rn(t)*ho(t). The low-pass equivalent (LPE) of the field at OF output, el(t), is defined as e(t)=√2ℜ{el(t)exp[j2πν 0 t]}=√2ℜ{[sl(t)+nl(t)]exp[j2πν 0 t]}, where sl(t) and nl(t) are the LPE of s(t) and n(t), respectively, and ℜ{z} is the real part of z. The LPE of the OF impulse response, ho,l(t), and transfer function, Ho,l(f), are defined as ho(t)=2ℜ{ho,l(t)exp[j2πν 0 t]} and by Ho(ν)=Ho,l(ν-ν 0)+H*o,l(-ν-ν 0), respectively. Considering the representation of the LPE of each noise polarization in their I and Q components, the LPE of the field at the OF output is given by:

e1(t)=[sl,(t)+nl,I,(t)+jnl,Q,(t)]u+[nl,I,(t)+jnl,Q,(t)]u

where s l,‖(t), n l,‖(t) and n l,⊥(t) are the ‖ and ⊥ components of sl(t) and nl(t). Hereinafter, the subscript l is dropped for the sake of notation simplicity, as the analytical derivations are performed with the LPE of the signals. As the OFDM signal field is completely polarized on the parallel direction, the subscript ‖ in the signal description is also dropped.

The derivation of the mean of the I and Q electrical components of the k-th subcarrier of the γ-th OFDM symbol at equalizer output from the field at OF output is performed in Appendix A. The mean of those components can be written, respectively, as:

me,I(γ)[k]={E[se(γ)[k]]}me,Q(γ)[k]={E[se(γ)[k]]}

where ℑ{z} is the imaginary part of z and E[s (γ) e[k]] is the mean of the signal at equalizer output given by Eq. 17 in Appendix A.

The variance of the k-th subcarrier of the γ-th OFDM symbol for the I and Q components at the FFT block output is also derived in Appendix A and it can be written as a function of the signal-noise beat terms variance and the noise-noise beat terms variance of each polarization:

(σFFT,(I,Q)(γ)[k])2=y=14(σsASE,(I,Q),y(γ)[k])2+(σASEASE,,(I,Q),y(γ)[k])2+(σASEASE,,(I,Q),y(γ)[k])2

where the first and the second subscripts in (I;Q) should be used for either the I or the Q component of each subcarrier. Each one of the four signal-noise beat terms are the variances of the beats between the I and Q components of the signal and the noise fields while the four noise-noise beat terms are the variances of the beats between the I and Q components of the noise field in each polarization direction. As shown in Appendix B, the variance of each signal-noise beat term of the k-th subcarrier of the γ-th OFDM symbol at FFT output is given by:

(σsASE,(I,Q),y(γ)[k])2=Rλ24+Sy(f1){exp(j2πfRFt0(γ))2H1,m,(I,Q)(γ)(f1+fRF,k)+
+exp(j2πfRFt0(γ))2H2,m,(I,Q)(γ)(f1fFR,k)}{exp(j2πfRFt0(γ))2H2,p,(I,Q)(γ)(f1fRF,k)+
+exp(j2πfRFt0(γ))2H1,p,(I,Q)(γ)(f1+fRF,k)}*df1

where the function H (γ) i,m,(I,Q)(f,k) represents the effects of the different phases of the local oscillator in the I and Q branches of the OFDM receiver and the bandwidth limitation imposed by the electrical low-pass filter (LPF) and the FFT block on the signal-noise beat terms. H (γ) i,m,(I,Q)(f,k) is defined by Eq. 33 in Appendix B. t (γ) 0 is the first sampling time instant of the γ-th OFDM symbol (see Appendix A) and Sy(f) is the PSD (or cross-PSD, depending on which signal-noise beat term is being evaluated) of the I and the Q noise components at PIN input.

Considering the same PSD for the I and Q noise components of each polarization direction and the same PSD for the noise in each polarization direction, we can drop the noise-noise beat terms dependence on the polarization in Eq. 3. As shown in Appendix C, the variance of each noise-noise beat term of each subcarrier at FFT block output is given by:

(σASEASE,(I,Q),y(γ)[k])2=18+[Sy(f1)*Sy(f1)]Heq,(I,Q)(γ)(f1,k)2df1

where H (γ) eq,(I,Q)(f 1,k) represents the equivalent transfer function of the OFDM receiver and describes the influence of the phase of the local oscillator and the bandwidth limitation imposed by the electrical LPF and the FFT block on the noise-noise beat term.

The variance of the I and Q components of each subcarrier at equalizer output can be expressed in terms of the variance of the signal at FFT output (as shown in Appendix D) as:

(σe,(I,Q)(γ)[k])2=(Ae(γ)[k])2(σFFT,(I,Q)(γ)[k]))2

where A (γ) e [k] is the amplitude response of the equalizer. Eq. 6 shows that the variance of the I and Q components of the signal at equalizer output is independent of the equalizer phase.

2.3. Performance evaluation

In order to take into account correctly the statistical distribution of the ISI and the ICI, the BER can be evaluated using an exhaustive Gaussian approach (GA) [12] for the I and Q components of each OFDM subcarrier at the equalizer output:

BER(I,Q)[k]=1Ns[γa(I,Q)(γ)=0=1NsQ(F(I,Q)[k]me,(I,Q),0(γ)[k]σe,(I,Q)(γ)[k])+γa(I,Q)(γ)=1=1NsQ(me,(I,Q),1(γ)[k]F(I,Q)[k])σe,(I,Q)(γ)[k])]

where m (γ) e,(I,Q),0[k] and m (γ) e,(I,Q),1[k] are the mean and σ (γ) e,(I,Q)[k] is the standard deviation of the k-th subcarrier of the γ-th OFDM symbol at the equalizer output when the bit transmitted in the (I,Q) component, a (γ) (I,Q), is 0 or 1, respectively, and F (I,Q)[k] is the decision threshold level. Although Eq. 7 considers binary phase shift keying (BPSK) or quadrature phase shift keying (QPSK), it can be easily generalized to other more efficient mappings. Notice that k∈Ω, where Ω is the subset of the subcarriers indexes corresponding only to the information subcarriers, and 1≤γNs, γ∈ℕ, where Ns is the number of simulated OFDM symbols.

If BPSK symbol mapping is used (with symbols j and -j), the means and the standard deviations considered in Eq. 7 are relative to the Q component only. Instead, if QPSK is used and assuming that the I and Q signal components at the equalizer output are uncorrelated and that Gray mapping is used, the BER of each OFDM subcarrier is given by BER[k]=12[1(1BERI[k])(1BERQ[k])], where BERI [k] and BERQ[k] are the BER of the I and Q components of each subcarrier at the equalizer output given by Eq. 7. The overall BER is evaluated averaging the BER of each subcarrier over all information subcarriers, Ni:

BER=k=1(kΩ)N×BER[k]Ni

3. Results and analysis for OFDM-UWB radio signals

In this section, the accuracy of the SAGA is assessed and discussed considering the transmission of OFDM-UWB radio signals [13]. The OFDM-UWB baseband signal is composed of 128 subcarriers: 100 are used as information subcarriers, 12 are used as pilot subcarriers, 10 are used as guard subcarriers, and 6 are used as null subcarriers. The OFDM-UWB signal bandwidth is 528 MHz and the time interval of each OFDM-UWB symbol is 312.5 ns. The I and Q components of the generated OFDM-UWB signal are filtered by 6th order Bessel LPFs with a -3 dB bandwidth of 400 MHz to reduce the aliasing components. The electro-optic conversion is modeled by a Mach-Zehnder modulator (MZM) biased at quadrature point and the analysis performed considers OFs tuned to the optical carrier. Fig. 2(a) shows the spectra of the double sideband OFDM-UWB signal with fRF=3.4 GHz applied to MZM arms and at MZM output.

Null phases for the I and Q branches of the OFDM receiver are considered (ϕ 1=ϕ 2=0). The results are obtained in back-to-back, with R λ=1 A/W. In almost all the results presented, BPSK symbol mapping is used. In other cases, an explicit reference to the symbol mapping used is presented. The transfer function of the equalizer is estimated by using the information provided by the channel distortion induced on the OFDM-UWB pilots subcarriers (spread along the OFDM signal band). A linear regression method is implemented to estimate the equalizer transfer function along the band of the transmitted subcarriers.

3.1. Probability density function of the signal at equalizer output

In order to assess the accuracy of the proposed GA, the probability density function (PDF) of the I and Q components of each subcarrier at the equalizer output is evaluated by MC simulation and compared with the one provided by the GA. 32 OFDM-UWB symbols are considered in numerical results and fRF=3.4 GHz is used to generate the OFDM-UWB radio signal. Two 6th order Bessel LPFs with a -3 dB bandwidth of 400 MHz are used at the OFDM-UWB electrical receiver. The OF is modeled by a 2nd order super-Gaussian filter with -3 dB bandwidth of 30 GHz. The estimation of the PDF for each subcarrier is performed over 3.7×106 noise runs.

Fig. 2(b) and 2(c) show the PDFs of the Q component of the signal at the equalizer output obtained by MC and by the GA, for modulation indexes of 19% and 64%, respectively, and for an optical signal-to-noise ratio (OSNR), defined in a reference bandwidth of 0.1 nm, of 7 dB. The modulation index is defined as m=VRMS/Vb, where VRMS is the root-mean-square (rms) voltage of the OFDM-UWB signal applied to the MZM arms and Vb is the MZM bias point.

 figure: Fig. 2.

Fig. 2. a) Spectra of the OFDM-UWB radio signal applied to MZM arms and at the MZM output, with ν0 the optical carrier frequency. PDF of the Q component of two subcarriers of the signal at the equalizer output for a) a modulation index of 19% and b) a modulation index of 64%. PDF obtained by using the GA (lines) and MC simulation (marks).

Download Full Size | PDF

Fig. 2(b) and 2(c) show the PDFs of two subcarriers of the first OFDM-UWB received symbol that carry a negative and a positive BPSK symbol. The PDF of the OFDM-UWB signal is well approximated by the proposed GA. Particularly, a very good accuracy can be observed in the tail regions for both modulation indexes, confirming that the statistics of the signal at equalizer output is adequately modeled by a Gaussian distribution.

3.2. Mean and variance of the signal at equalizer output

In this section, the mean and the variance of each subcarrier at the equalizer output is evaluated by MC simulation and compared with the mean and the variance given by the GA. The OFDM-UWB signal is generated with the same features as the ones described in the previous section and 2×104 runs are used to estimate the mean and the variance of each subcarrier. The results are evaluated for the first OFDM-UWB symbol and for the OFDM-UWB signal located into the UWB sub-band 1 (fRF=3.4 GHz) or the UWB sub-band 14 (fRF=10.3 GHz).

 figure: Fig. 3.

Fig. 3. a) Mean of the Q component of the first OFDM symbol at equalizer output as a function of the transmitted information subcarrier. Results obtained by using the GA (×) and MC simulation (circles). b) Variance contributions to the total variance of the Q component of the first OFDM symbol at the equalizer output as a function of the transmitted information subcarrier. Results obtained by using the GA (lines) and MC simulation (marks).

Download Full Size | PDF

Fig. 3(a) shows the mean of the Q component of the first OFDM symbol at the equalizer output as a function of the information subcarriers index. The subcarriers with indexes from 1 to 50 correspond to subcarriers transmitted in the frequencies below fRF while subcarriers with indexes from 51 to 100 correspond to subcarriers transmitted in frequencies above fRF (indexes 1 and 100 correspond to the information subcarriers transmitted on the edges of the OFDM spectrum). An excellent agreement between the results obtained by MC and by the GA can be observed, confirming that the mean at the equalizer output is well described by Eq. 2. Fig. 3(b) depicts the variance of the first OFDM symbol and the different beat terms contributions at the equalizer output obtained by MC simulation and by the GA as a function of the information subcarriers. In the GA, the noise-noise beat term (in both polarizations) at equalizer output is obtained by using Eq. 3, 5 and 6. The optical carrier-noise beat term is obtained by switching-off the voltage signal applied to the MZM arms and by subtracting the noise-noise beat contributions from the total variance. The comparison between the MC and the GA results shows that the beat terms contributions to the total variance of each subcarrier at equalizer output are very well estimated by the closed-form expressions given by the GA. Fig. 3(b) shows also that the main contribution to the total variance of each subcarrier comes from the optical carrier-noise beat term. Fig. 3(b) does not show the OFDM signal-noise beat as its contribution to the total variance is negligible.

Fig. 3(b) shows that the variance at the equalizer output is higher for the subcarriers with lower and higher indexes (subcarriers located on the edges of the OFDM-UWB spectrum). The filters shape reduces the amplitude of the OFDM subcarriers located at the edges of the spectrum requiring an higher amplitude response of the equalizer transfer function on those frequencies. Therefore, though lower noise PSD is observed for the subcarriers located on the edges of the OFDM-UWB spectrum at the equalizer input due to the non-rectangular frequency response of the optical and, mainly, of the electrical filters, the noise PSD at the equalizer output is higher for those subcarriers leading to lower signal-to-noise ratio. This is in agreement with the experimental results shown in [5]. The results of Fig. 3(b) are different from the ones presented in [10] (the PSD of the noise decreases with the increase of the subcarrier frequency) as, in that case, the noise PSD is affected only by the optical filter and the influence of the equalizer transfer function on the noise PSD has not been considered.

 figure: Fig. 4.

Fig. 4. Variance of the Q component of the first OFDM symbol at the equalizer output as a function of the transmitted information subcarrier for a) different types of LPF and modulation indexes values and b) different types of optical filters. Variance obtained by the GA (lines) and the MC simulation (marks).

Download Full Size | PDF

Fig. 4(a) shows the variance of the Q component of the first OFDM symbol at the equalizer output as a function of the transmitted information subcarrier considering the OF modeled by a rectangular filter with a -3dB bandwidth of 30 GHz and OSNR=7 dB. Fig. 4(a) shows that the influence of the electrical LPF of the OFDM-UWB receiver on the variance of each subcarrier at the equalizer output given by Eq. 6 is correctly taken into account. In addition, Fig. 4(a) shows also that the variance changes considerably from subcarrier to subcarrier when different types of LPF are used. This change is due to the combined effect of two main features: i) different LPF shapes lead to different gains of each OFDM subcarrier. These gains are compensated by the equalizer transfer function affecting the noise variance at the equalizer output, as indicated by Eq. 6. ii) As shown in Appendix B, the noise power at the equalizer output is imposed by an equivalent filter given by the magnitude of the LPF transfer function and by the magnitude of the equivalent FFT transfer function, T (γ) (1,2)(f,k).

Fig. 4(b) shows results similar to Fig. 4(a), but considering different types of OFs and 6th order Bessel LPFs with a -3 dB bandwidth of 400 MHz at the electrical receiver. These results show that the impact of the OF shape on the variance of each subcarrier at the equalizer output is also correctly described by Eq. 6. The different variance levels obtained for the different OF types is due the influence of the OF shape on the noise power and on the subcarriers amplitude.

3.3. BER

 figure: Fig. 5.

Fig. 5. BER at the equalizer output as a function of a) the information subcarriers and b) the modulation index. Results obtained by using the GA (lines) and MC simulation (marks).

Download Full Size | PDF

Fig. 5(a) depicts the BER evaluated by MC and by the GA as a function of the information subcarriers for different optical and electrical filters. In MC simulation, the BER is computed when 100 errors occur in the subcarrier with better performance. The BER changes considerably from subcarrier to subcarrier mainly due to the dependence of the variance of each subcarrier at equalizer output on the amplitude of the frequency response of the optical and electrical filters. As shown in [5], the subcarriers more affected by the amplitude reduction induced by the channel response present higher BER levels. Fig. 5(b) shows the BER obtained by MC and by Eq. 8, as a function of the modulation index for BPSK and QPSK symbol mappings. 2nd order super-Gaussian OF with a -3 dB bandwidth of 30 GHz and 6th order Bessel LPFs with a -3 dB bandwidth of 400 MHz are considered. The results are presented for two different OSNR values. Fig. 5(b) indicates that the SAGA provides excellent BER estimates for all modulation indexes, i. e., in case of performance dominantly impaired by noise or by interference, or simultaneously by both effects.

4. Results and analysis of OFDM signals proposed for long-haul networks

In this section, the accuracy of the SAGA is studied for two different OFDM transmitters (both using DD in the receiver) recently proposed for long-haul networks: the RF tone assisted gapped OFDM system (referred henceforth as A-OFDM transmitter) [14] and ii) the OFDM system proposed in [15] (referred henceforth as B-OFDM transmitter).

QPSK symbol mapping and 10 Gbit/s OFDM signals are generated in both transmitters. A-OFDM and B-OFDM transmitters use 128 and 512 data subcarriers, respectively, with 2-times oversampling used to obtain brick-wall spectra. The electrical OFDM signal is converted to the optical domain by a linear IQ modulator and back-to-back operation with R λ=1 A/W are considered in the analysis. It should be stressed that SSB-OFDM spectra are obtained at the IQ modulator output in both implementations. The receiver used for both systems is the one presented in Fig. 1. Rectangular optical and electrical filters with a bandwidth of 15 GHz and 12 GHz, respectively, are considered in the receiver of the system employing the A-OFDM transmitter, while the receiver of the system employing the B-OFDM transmitter uses rectangular optical and electrical filters with a bandwidth of 25 GHz and 3 GHz. In both cases, the OF is centered at ν 0=193.1 THz. OFDM signals with 5 GHz bandwidth (centered around 7.5 GHz) at the PIN output are obtained in both configurations with a 5 GHz gap to accommodate for the distortion induced by the subcarriers-subcarriers beat term [14, 15].

For the system employing the B-OFDM transmitter, the down conversion process at the electrical receiver is performed using fRF=7.5 GHz. For the A-OFDM transmitter, the down conversion process is performed using two different approaches: i) using fRF=5 GHz and ii) using an FFT implemented over twice the number of points used in the IFFT of the B-OFDM transmitter. The latter approach avoids the I and Q down conversion process presented in Fig. 1. In the proposed GA, this is equivalent to set fRF=0. The equalizer transfer function is estimated using training symbols, and in both systems, the BER given by the GA and DEC is evaluated as a function of the carrier-to-signal power ratio (CSPR) considering ϕ 1=ϕ 2=0.

 figure: Fig. 6.

Fig. 6. a) Spectra of the OFDM signal at the IQ modulator output (I and III) and at the PIN output (II and IV) considering the A-OFDM and the B-OFDM transmitters. b) BER at the equalizer output as a function of CSPR. Results obtained by using the GA (lines) and MC simulation (marks), with the approach i) of OFDM-A (diamond marks and dashed lines) and the approach ii) of OFDM-A (square marks and dashed-dotted lines).

Download Full Size | PDF

Fig. 6(a) depicts the spectra of the OFDM signals at the IQ modulator output (I and III) and at the PIN output (II and IV) for the A-OFDM and B-OFDM transmitters. The spectra at the IQ modulator output confirm that both transmitters generate SSB-OFDM signals while the spectra at the PIN output show the unwanted subcarriers-subcarriers components falling close to DC and the useful OFDM signal located around 7.5 GHz [14, 15]. Fig. 6(b) shows the BER obtained using Eq. 8 as a function of the CSPR for two OSNR levels and for systems employing the A-OFDM and B-OFDM transmitters. The comparison between the results given by the GA and DEC clearly shows that the system performance is very well estimated by the GA, independently of the OFDM transmitter considered. Furthermore, Fig. 6(b) shows that for low OSNR levels the optimum CSPR is around 0 dB, as presented in [16].

5. Conclusion

A SAGA for each subcarrier of each OFDM symbol at the equalizer output, that allows evaluating the BER of each OFDM subcarrier in optically pre-amplified and DD receivers has been proposed. The approach is valid for arbitrary OFDM transmitters as long as the OFDM signal can be demodulated by the considered receiver structure. Closed-form analytical expressions for the mean and variance of the signal-noise and noise-noise beat terms of each subcarrier at equalizer output have been derived. These expressions take into account the distortions induced by the frequency response of the optical and electrical filters, the different phases of the local oscillator in the I and Q branches of the OFDM receiver and the equalizer transfer function.

The results obtained by the proposed method have shown excellent agreement with MC estimates for OFDM-UWB radio signals used in access networks and for two SSB-OFDM signals proposed for long-haul systems.

Acknowledgments

The work of Tiago Alves was supported by Fundação para a Ciência e a Tecnologia from Portugal under contract SFRH/BD/29871/2006 and the projects TURBO-PTDC/EEA-TEL/104358/2008 and SAMURAI-PTDC/EEA-TEL/108165/2008. This work was also supported in part by the UCELLS-FP7-IST-1-216785 project.

A. Mean and variance of the signal at the equalizer output

Considering a quadratic characteristic for the PIN, the current at the PIN output is given by:

sPIN(t)=Rλel(t)·el*(t)=Rλ{s(t)2+2sI(t)nI,(t)+2sQ(t)nQ,(t)+nI,2(t)+nQ,2(t)+
+nI,2(t)+nQ,2(t)}

with sI(t)=ℜ{s(t)}, sQ(t)=ℑ{s(t)} and |s(t)|2=s 2 I(t)+s 2 Q(t). Considering the scheme of the OFDM receiver depicted in Fig. 1, the current at LPF output in the I and Q branches of the electrical OFDM receiver can be expressed as:

sLPF,I(t)=[sPIN(t)2cos(2πfRFt+ϕ1)]*hr(t)=Rλ2+[s(τ)2×+2sI(τ)nI,(τ)+
+2sQ(τ)nQ,(τ)+nI,2(τ)+nQ,2(τ)+nI,2(τ)+nQ,2(τ)]cos(2πfRFτ+ϕ1)hr(t-τ)dτ
sLPF,Q(t)=[sPIN(t)2sin(2πfRFt+ϕ2)]*hr(t)=Rλ2+[s(τ)2+2sI(τ)nI,(τ)+
+2sQ(τ)nQ,(τ)+nI,2(τ)+nQ,2(τ)+nI,2(τ)+nQ,2(τ)]sin(2πfRFτ+ϕ2)hr(t-τ)dτ

where hr(t) is the impulse response of the LPF of the electrical receiver and the 1/2 factor is due to the current division by the two branches of the electrical receiver.

Prior to the FFT block, the signal is sampled and the guard time is removed. Thus, the continuous time, t, is replaced by t (γ) 0+(i-1)Tc, where t (γ) 0=t 0+(γ-1)T is the first sampling time instant of the γ-th OFDM symbol, T is the OFDM symbol time, Tc=(T-Tg)/N is the chip time and Tg is the guard time. The index i stands for the i-th chip and thus, 1≤iN, i∈ℕ. The signal at FFT block input can be expressed as:

sLPF(t0(γ)+(i1)Tc)=+sPIN(τ)2[cos(2πfRFτ+ϕ1)jsin(2πfRFτ+ϕ2)]×
×hr(t0(γ)+(i1)Tcτ)

Hereinafter, we define s (γ) LPF[i]=sLPF(t (γ) 0+(i-1)Tc) for the sake of notation simplicity. The signal at FFT block output corresponding to the k-th subcarrier of the γ-th symbol of the received signal can then be written as:

sFFT(γ)[k]=sFFT,I(γ)[k]+jsFFT,Q(γ)[k]=i=1NsLPF(γ)[i]exp[jWN(k,i)]

where WN(k,i)=2π(k-1)(i-1)/N. From Eq. 12, the I and Q components of the signal at the FFT output are given by:

sFFT,(I,Q)(γ)[k]=12i=1N+sPIN(τ)f(I,Q)(τ,k,i)hr(t0(γ)+(i1)Tcτ)

where

fI(τ,k,i)=cos(2πfRFτ+ϕ1)cos[WN(k,i)]sin(2πfRFτ+ϕ2)sin[WN(k,i)]
fQ(τ,k,i)=cos(2πfRFτ+ϕ1)sin[WN(k,i)]sin(2πfRFτ+ϕ2)cos[WN(k,i)]

From Eq. 12 and by evaluating the mean of the current at PIN output using Eq. 9, the mean of the signal at FFT block output is given by:

E[sFFT(γ)[k]]=Rλ2i=1N{{s(t0(γ)+(i1)Tc)2[cos(2πfRF(t0(γ)+(i1)Tc)+ϕ1)
jsin(2πfRF(t0(γ)+(i1)Tc)+ϕ2)]}*hr(t0(γ)+(i1)Tc)+
+Pn[{Hr(fRF)exp[j2πfRF[t0(γ)+(i1)Tc]]exp(1)}{Hr(fRF)×
×exp[j2πfRF[t0(γ)+(i1)Tc]]exp(2)}]}exp[jWN(k,i)]

where Pn is the optical noise power at the OF output. Eq. 16 shows that if the LPF bandwidth is narrow when compared with the RF carrier frequency the contribution of the noise-noise beat term to the mean of the signal at FFT block output can be neglected.

Considering the equalizer transfer function of the k-th subcarrier of the γ-th OFDM symbol, H (γ) e[k]=A (γ) e[k]exp(-jΦ(γ) e [k]), the mean of the signal at the equalizer output can be expressed as:

E[se(γ)[k]]=E[sFFT(γ)[k]]Ae(γ)[k]exp(jΦe(γ)[k])

Hence, the mean of the I and Q components of each subcarrier at equalizer output are given by:

me,I(γ)[k]={E[se(γ)[k]]}me,Q(γ)[k]={E[se(γ)[k]]}

The variance of the k-th subcarrier of the γ-th OFDM symbol for the I and Q components at the FFT block output can be evaluated as:

(σFFT,(I,Q)(γ)[k])2=E[(sFFT,(I,Q)(γ)[k])2]{E[sFFT,(I,Q)(γ)[k]]}2

From Eq. 16, the mean of the I and Q components of the k-th subcarrier of the γ-th OFDM symbol at FFT block output can be expressed as:

E[sFFT,(I,Q)(γ)[k]]=12i=1N+E[sPIN(τ)]f(I,Q)(τ,k,i)hr(t0(γ)+(i1)Tcτ)

By using Eq. 9, the mean of the current at PIN output is given by:

E=[sPIN(τ)]=Rλ{s(τ)2+E[nI,2(τ)]+E[nQ,2(τ)]+E[nI,2(τ)]+E[nQ,2(τ)]}

Considering the statistical independence of the noise over both polarization directions, i. e., E{n x,‖ n y,⊥]=E[n x,‖] E[n y,⊥] where x and y stands for the I and Q components, the odd order moments of Gaussian processes with zero mean are null and that the fourth order moments of Gaussian processes with zero mean are given by E[X 1 X 2 X 3 X 4]=E[X 1 X 2] E[X 3 X 4]+E[X 1 X 3] E[X 2 X 4]+E[X 1 X 4] E[X 2 X 3], the variance of the I and Q components of each subcarrier at FFT output can be written as:

(σFFT,(I,Q)(γ)[k])2=Rλ2i=1Nn=1N++{+[sI(τ1)sI(τ2)+SnI,nI,(f1)+sI(τ1)×
×sQ(τ2)SnI,nQ,(f1)+sQ(τ1)sI(τ2)SnQ,nI,(f1)+sQ(τ1)sQ(τ2)SnQ,nQ,(f1)]×
×exp[j2πf1(τ1τ2)]df1+12++[SnI,nI,(f2)SnI,nI,(f3)+SnI,nQ,(f2)SnI,nQ,(f3)+
+SnQ,nI,(f2)SnQ,nI,(f3)+SnQ,nQ,(f2)SnQ,nQ,(f3)+SnI,nI,(f2)SnI,nI,(f3)+
+SnI,nQ,(f2)SnI,nQ,(f3)+SnQ,nI,(f2)SnQ,nI,(f3)+SnQ,nQ,(f2)SnQ,nQ,(f3)]×
×exp[j2π(f2+f3)(τ1τ2)]df2df3}g(I,Q)(τ1,τ2,k,i,n)dτ1dτ2

where SnI,ynI,y(f), SnQ,ynQ,y(f), SnI,ynQ,y(f) and SnQ,ynI,y(f) are the PSD of nI,y(t) and nQ,y(t), and the cross PSD between nI,y(t) and nQ,y(t) and between nQ,y(t) and nI,y(t), respectively, at the OF output. y represents either the ‖ or the ⊥ polarization direction. The derivation of Eq. 22 used the Wiener-Kintchine theorem. For tuned OF, the PSD of nI(t) or nQ(t) can be evaluated from the PSD at the OF input as Sn(I,Q)n(I,Q),y(f)=0.5×SASE,yHo,l(f)2, while the cross PSD between both noise field components at the OF output is null. If detuned OFs are used, the relation between the noise PSD at the OF output and at the OF input is described in [17]. The functions g (I,Q)(τ 1,τ 2,k, i,n) are given by:

g(I,Q)(τ1,τ2,k,i,n)=f(I,Q)(τ1,k,i)f(I,Q)(τ2,k,n)hr(t0(γ)+(i1)Tcτ1)×
×hr(t0(γ)+(n1)Tcτ2)

Eq. 22 can still be written as shown in Eq. 3.

B. Signal-noise beat contribution

From the variance of the I and Q components of the k-th subcarrier of the γ-th OFDM symbol at FFT block output given by Eq. 22 of Appendix A, the contribution of the I and Q signal-noise beat to the variance are given by:

(σsASE,(I,Q),y(γ)[k])2=Rλ2i=1Nn=1N+++sm(τ1)sp(τ2)Sy(f1)exp[j2πf1(τ1τ2)]df1×
×g(I,Q)(τ1,τ2,k,i,n)dτ1dτ2

where

sm(τ1)={sI(τ1)ify=1ory=2sQ(τ1)ify=3ory=4sp(τ2)={sI(τ2)ify=1ory=3sQ(τ2)ify=2ory=4
Sy(f1)={SnI,nI,(f1)ify=1SnI,nQ,(f1)ify=2SnQ,nI,(f1)ify=3SnQ,nQ,(f1)ify=4

Using Eq. 23, (σsASE,(I,Q),y(γ)[k])2 can be written as:

(σsASE,(I,Q),y(γ)[k])2=Rλ2i=1Nn=1N++Sy(f1)Sm(f2)F(I,Q)(f3,k,i)Hr(f4)Sp(f5)×
×F(I,Q)(f6,k,n)Hr(f7)exp[j2πf4(t0(γ)+(i1)Tc)]exp[j2πf7(t0(γ)+(n1)Tc)]×
×exp[j2π(f1+f2+f3f4)τ1]exp[j2π(f1+f5+f6f7)τ2]df1df7dτ1dτ2

where S (m,p)(f)=𝓕{s (m,p)(τ)}, F (I,Q)(f,k,i)=𝓕{f (I,Q)(τ,k,i)}and Hr(f)=𝓕 {hr(τ)g. From Eq. 14 and 15, F (I,Q)(f,k,i) are given by:

FI(f,k,i)=𝓕{f1(f,k,i)}=12[A(k,i)exp(1)+jB(k,i)exp(2)]δ(ffRF)+
+12[A(k,i)exp(1)jB(k,i)exp(2)]δ(f+fRF)
FQ(f,k,i)=𝓕{fQ(f,k,i)}12[B(k,i)exp(1)+jA(k,i)exp(2)]δ(ffRF)+
+12[B(k,i)exp(1)jA(k,i)exp(2)]δ(f+fRF)

where A(k,i)=cos[WN(k,i)] and B(k,i)=sin [WN(k,i)]. Thus, after performing some algebra, (σsASE,(I,Q),y(γ)[k])2 can be written as:

(σsASE,(I,Q),y(γ)[k])2Rλ24+Sy(f1){Rl,m,(I,Q)(γ)(f1fRF,k)+R2,m,(I,Q)(γ)(f1+fRF,k)}×
×{R2,p,(I,Q)(γ)(f1+fRF,k)+R1,p,(I,Q)(γ)(f1fRF,k)}*df1

where (z)* represents the complex conjugate of z. The functions R (γ) l,m,(I,Q)(-f 1-fRF,k) and R (γ) 2,m,(I,Q)(-f 1+fRF,k) are given by:

R1,m,(I,Q)(γ)(f1fRF,k)=exp(j2πfRFt0(γ))2H1,m,(I,Q)(γ)(f1+fRF,k)
R2,m,(I,Q)(γ)(f1+fRF,k)=exp(j2πfRFt0(γ))2H2,m,(I,Q)(γ)(f1fRF,k)

where the function H (γ) i,m,(I,Q)(f,k) is given by:

Hi,m,(I,Q)(γ)(f,k)=𝓕{νm(γ)(t)gi,(I,Q)(γ)(t,k)}gi,(I,Q)(γ)(t,k)=𝓕1{Gi,(I,Q)(γ)(f,k)}
νm(γ)=𝓕1{Vm(γ)(f)}=𝓕1{exp(j2πft0(γ))sm(f)}

The same procedure can be performed to obtain R (γ) 1,p,(I,Q)(-f 1-fRF,k) and R (γ) 2,p,(I,Q)(-f 1+fRF,k), respectively, by replacing m by p. The transfer functions G (γ) i,(I,Q)(f,k) can be represented as the elements of a column vector, G (γ) (I,Q), given by:

G(I,Q)(γ)=[G1,(I,Q)(γ)(f,k)G2,(I,Q)(γ)(f,k)]=Hr(f)P(I,Q)T(γ)

where the vector T (γ)=[T (γ) 1(f,k) T (γ) 2 (f,k)]T represents the effects of the FFT operation and is given by:

T1(γ)(f,k)=exp{j2πN12[fTc+(k1)N]}[sin{2π[fTc+(k1)N]N2}sin{π[fTc+(k1)N]}]
T2(γ)(f,k)=exp{j2πN12[fTc(k1)N]}[sin{2π[fTc(k1)N]N2}sin{π[fTc(k1)N]}]

PI and P Q are the I and Q matrixes that describe the mismatch of the phases of the local oscillator:

PI=[exp(1)+exp(2)exp(1)exp(2)exp(1)exp(2)exp(1)+exp(2)]
PQ=j[exp(1)+exp(2)exp(1)+exp(2)exp(1)exp(2)exp(1)exp(2)]

The PI and P Q matrixes allow assessing the orthogonality between the I and Q branches of the OFDM receiver and, consequently, the orthogonality between the OFDM subcarriers at FFT output. PI and P Q become diagonal matrixes if there is not phase difference between the I and Q branches (ϕ 1=ϕ 2).

The transfer functions G (γ) i,(I,Q)(f,k) describe the effects of the different phases of the local oscillator in the I and Q branches, of the electrical filtering and of the FFT operation of the OFDM receiver on the signal-noise beat terms. It should be stressed that Eq. 35, 36 and 37 indicate that the noise at the equalizer output is limited by the LPF and also by the FFT block due to the magnitude response of terms T (γ) 1(f,k) and T (γ) 2(f,k).

Substitution of Eq. 31 and 32 in Eq. 30 allow writing the signal-noise beat terms variance of the k-th subcarrier of the γ-th OFDM symbol at FFT block output as shown in Eq. 4.

C. Noise-noise beat contribution

From Eq. 3 and 22, the I and Q noise-noise beat contribution to the variance of the k-th subcarrier of the γ-th OFDM symbol at FFT output are given by:

(σASEASE,(I,Q),y(γ)[k])2=Rλ22i=1Nn=1N++++Sy(f2)Sy(f3)×
×exp[j2π(f2+f3)(τ1τ2)]df2df3g(I,Q)(τ1,τ2,k,i,n)dτ1dτ2

The simplification of Eq. 40 follows a similar procedure to the one used in Appendix B to simplify the signal-noise beat contribution. Each noise-noise beat term can be written as:

(σASEASE,(I,Q),y(γ)[k])=18+[Sy(f1)*Sy(f1)]Heq,(I,Q)(γ)(f1,k)2df1

where H (γ) eq,(I,Q)(f 1,k)=H (γ) eq,1,(I,Q)(f 1+fRF,k)+H (γ) eq,2,(I,Q)(f 1-fRF,k) and with:

[Heq,1,(I,Q)(γ)(f,k)Heq,2,(I,Q)(γ)(f,k)]=Hr(f)2exp[j2πft0(γ)]P(I,Q)T(γ)

P (I,Q) and T (γ) have been already defined in Appendix B.

D. Variance at the equalizer output

The signal at equalizer output can be written as s (γ) e[k]=s (γ) FFT [k]H (γ) e[k]. Thus, the I and Q components of the each subcarrier at the equalizer output are given by:

se,I(γ)[k]=Ae(γ)[k]cos(Φe(γ)[k])sFFT,I(γ)[k]+Ae(γ)[k]sin(Φe(γ)[k])sFFT,Q(γ)[k]
se,Q(γ)[k]=Ae(γ)[k]cos(Φe(γ)[k])sFFT,Q(γ)[k]Ae(γ)[k]sin(Φe(γ)[k])sFFT,I(γ)[k]

The variance of the I and Q components of each subcarrier can then be evaluated following a similar procedure to the one of Eq. 19. Assuming that the frequency dependency of the ASE noise PSD along the OFDM signal bandwidth is small and using the relation between the elements of the PI and P Q matrixes defined in Appendix B, it can be shown that the covariance between the I and Q components of each subcarrier at FFT output is null and thus, E[s (γ) FFT,(I,Q)[k]s (γ) FFT,(Q,I)[k]]=E[s (γ) FFT,(I,Q)[k]]E[s (γ) FFT,(Q,I)[k]]. Hence, the variance of the I and Q components of the k-th subcarrier of the γ-th OFDM symbol at the equalizer output can be written (α (γ) e,(I,Q)[k])2=(A (γ) e[k])2 (σ (γ) FFT,(I,Q)[k])2 where we have used that (σ (γ) FFT,I[k])2=(σ (γ) FFT,Q[k])2.

References and links

1. J. Armstrong, “OFDM for optical communications”, J. Lightwave Technol. 27, pp. 189–204 (2009). [CrossRef]  

2. R. Llorente, T. Alves, M. Morant, M. Beltran, J. Perez, A. Cartaxo, and J. Marti, “Ultra-wideband radio signals distribution in FTTH networks”, Photon. Technol. Lett. 20, pp. 945–947 (2008). [CrossRef]  

3. C. Lin, Y. Lin, J. Chen, S. Dai, P. Shih, P. Peng, and S. Chi, “Optical direct-detection OFDM signal generation for radio-over-fibre link using frequency doubling scheme with carrier supression”, Opt. Express 16, pp. 6056–6063 (2008). [CrossRef]  

4. J. Tang, P. Lane, and K. Shore, “30 Gbit/s transmission over 40 km directly modulated DFB laser-based SMF links without optical amplification and dispersion compensation for VSR and metro applications”, in Optical Fibre Communication Conference, OSA Technical Digest Series (CD) (Optical Society of America, 2006), paper JThB8. [CrossRef]  

5. B. Schmidt, A. Lowery, and J. Armstrong, “Experimental demonstration of 20 Gbit/s direct-detection optical OFDM and 12 Gbit/s with a colorless transmitter”, in Optical Fibre Communication Conference, OSA Technical Digest Series (CD) (Optical Society of America, 2007), paper PDP18.

6. M. Schuster, S. Randel, A. Bunge, S. Lee, F. Breyer, B. Spinnler, and K. Petermann, “Spectrally efficient compatible single-sideband modulation for OFDM transmission with direct detection”, Photon. Technol. Lett. 20, pp. 670–672 (2008). [CrossRef]  

7. T. Alves and A. Cartaxo, “Performance degradation due to OFDM-UWB radio signal transmission along dispersive single-mode fiber”, Photon. Technol. Lett. 21, pp. 158–160 (2009). [CrossRef]  

8. F. Buchali and R. Dischler, “Optimized sensitivity direct detection O-OFDM with multi level subcarrier modulation”, in Optical Fibre Communication Conference, OSA Technical Digest Series (CD) (Optical Society of America, 2008), paper OMU5.

9. A. Lowery and J. Armstrong, “Orthogonal-frequency-division multiplexing for dispersion compensation of long-haul optical systems”, Opt. Express 14, pp. 2079–2084 (2006). [CrossRef]   [PubMed]  

10. W. Peng, K. Feng, S. Chi, and A. Willner, “Bit error rate calculation for a single sideband OFDM signal with direct detection optically pre-amplified receivers”, in Conference on Lasers and Electro-Optics, OSA Technical Digest Series (CD) (Optical Society of America, 2008), paper CWN3.

11. A. Lowery, “Amplified-spontaneous noise limit of optical OFDM lightwave systems”, Opt. Express 16, pp. 860–865 (2008). [CrossRef]   [PubMed]  

12. J. Rebola and A. Cartaxo, “Gaussian approach for performance evaluation of optically preamplified receivers with arbitrary optical and electrical filters”, IEE Proc.-Optoelectronics 148, pp. 135–142 (2001). [CrossRef]  

13. High Rate UltraWideband PHY and MAC Standard, 2nd ed. Geneve, Switzerland: ECMA Int. (2007).

14. W. Peng, X. Wu, V. Arbab, K. Feng, B. Shamee, L. Christen, J. Yang, A. Willner, and S. Chi, “Theoretical and experimental investigations of direct-detected RF-tone-assisted optical OFDM systems”, J. Lightwave Technol. 27, pp. 1332–1339, May 15 (2009).

15. B. Schmidt, A. Lowery, and L. Du, “Low sample rate transmitter for direct-detection optical OFDM”, in Optical Fibre Communication Conference, OSA Technical Digest Series (CD) (Optical Society of America, 2009), paper OWM4.

16. S. Jansen, I. Morita, and H. Tanaka, “Carrier-to-signal power ratio in fiber-optic SSB-OFDM transmission systems”, in IEICE General Conference, Nagoya, Japan, (Institute of Electronics, Information and Communication Engineers, 2007), paper B-10-24.

17. A. Carlson, P. Crilly, and J. Rutledge, Communication Systems: An Introduction to Signals and Noise in Electrical Communication, (McGraw-Hill, New York, 2002).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1.
Fig. 1. Block diagram of the optically pre-amplified DD-OFDM receiver.
Fig. 2.
Fig. 2. a) Spectra of the OFDM-UWB radio signal applied to MZM arms and at the MZM output, with ν0 the optical carrier frequency. PDF of the Q component of two subcarriers of the signal at the equalizer output for a) a modulation index of 19% and b) a modulation index of 64%. PDF obtained by using the GA (lines) and MC simulation (marks).
Fig. 3.
Fig. 3. a) Mean of the Q component of the first OFDM symbol at equalizer output as a function of the transmitted information subcarrier. Results obtained by using the GA (×) and MC simulation (circles). b) Variance contributions to the total variance of the Q component of the first OFDM symbol at the equalizer output as a function of the transmitted information subcarrier. Results obtained by using the GA (lines) and MC simulation (marks).
Fig. 4.
Fig. 4. Variance of the Q component of the first OFDM symbol at the equalizer output as a function of the transmitted information subcarrier for a) different types of LPF and modulation indexes values and b) different types of optical filters. Variance obtained by the GA (lines) and the MC simulation (marks).
Fig. 5.
Fig. 5. BER at the equalizer output as a function of a) the information subcarriers and b) the modulation index. Results obtained by using the GA (lines) and MC simulation (marks).
Fig. 6.
Fig. 6. a) Spectra of the OFDM signal at the IQ modulator output (I and III) and at the PIN output (II and IV) considering the A-OFDM and the B-OFDM transmitters. b) BER at the equalizer output as a function of CSPR. Results obtained by using the GA (lines) and MC simulation (marks), with the approach i) of OFDM-A (diamond marks and dashed lines) and the approach ii) of OFDM-A (square marks and dashed-dotted lines).

Equations (67)

Equations on this page are rendered with MathJax. Learn more.

e1(t)=[sl,(t)+nl,I,(t)+jnl,Q,(t)]u+[nl,I,(t)+jnl,Q,(t)] u
me,I(γ)[k]={E[se(γ)[k]]}me,Q(γ)[k]={E[se(γ)[k]]}
(σFFT,(I,Q)(γ)[k])2=y=14(σsASE,(I,Q),y(γ)[k])2+(σASEASE,,(I,Q),y(γ)[k])2 + (σASEASE,,(I,Q),y(γ)[k])2
(σsASE,(I,Q),y(γ)[k])2 = Rλ24 +Sy(f1) { exp(j2πfRFt0(γ))2 H1,m,(I,Q)(γ) (f1+fRF,k) +
+exp(j2πfRFt0(γ))2H2,m,(I,Q)(γ)(f1fFR,k)}{exp(j2πfRFt0(γ))2H2,p,(I,Q)(γ)(f1fRF,k)+
+exp(j2πfRFt0(γ))2H1,p,(I,Q)(γ)(f1+fRF,k)}*df1
(σASEASE,(I,Q),y(γ)[k])2=18 +[Sy(f1)*Sy(f1)]Heq,(I,Q)(γ)(f1,k)2df1
(σe,(I,Q)(γ)[k])2=(Ae(γ)[k])2 (σFFT,(I,Q)(γ)[k]))2
BER(I,Q)[k]=1Ns[γa(I,Q)(γ)=0=1NsQ(F(I,Q)[k]me,(I,Q),0(γ)[k]σe,(I,Q)(γ)[k])+γa(I,Q)(γ)=1=1NsQ(me,(I,Q),1(γ)[k]F(I,Q)[k])σe,(I,Q)(γ)[k])]
BER= k=1(kΩ)N × BER[k]Ni
sPIN(t)=Rλel(t)·el*(t)=Rλ { s(t)2 + 2sI (t) nI, (t)+2sQ (t) nQ, (t) + nI,2 (t) + nQ,2 (t)+
+ nI,2 (t) + nQ,2 (t) }
sLPF,I(t)=[sPIN(t)2cos(2πfRFt+ϕ1)] * hr (t) = Rλ2 +[s(τ)2×+2sI(τ)nI,(τ) +
+ 2sQ (τ) nQ, (τ) + nI,2 (τ) + nQ,2 (τ)+nI,2(τ) + nQ,2 (τ) ] cos(2πfRFτ+ϕ1)hr(t-τ)dτ
sLPF,Q (t) = [sPIN(t)2sin(2πfRFt+ϕ2)] * hr (t)=Rλ2 +[s(τ)2+2sI(τ)nI,(τ)+
+ 2sQ (τ)nQ,(τ)+nI,2(τ)+nQ,2(τ)+nI,2 (τ) + nQ,2 (τ)]sin(2πfRFτ+ϕ2)hr(t-τ)dτ
sLPF (t0(γ)+(i1)Tc)=+sPIN(τ)2[cos(2πfRFτ+ϕ1)jsin(2πfRFτ+ϕ2)]×
×hr(t0(γ)+(i1)Tcτ)
sFFT(γ) [k] = sFFT,I(γ) [k]+jsFFT,Q(γ)[k]=i=1NsLPF(γ)[i]exp[jWN(k,i)]
sFFT,(I,Q)(γ) [k]=12i=1N+sPIN(τ)f(I,Q)(τ,k,i)hr(t0(γ)+(i1)Tcτ)
fI (τ,k,i)=cos(2πfRFτ+ϕ1)cos[WN(k,i)]sin(2πfRFτ+ϕ2)sin[WN(k,i)]
fQ (τ,k,i)=cos(2πfRFτ+ϕ1)sin[WN(k,i)]sin(2πfRFτ+ϕ2)cos[WN(k,i)]
E [sFFT(γ)[k]] = Rλ2 i=1N{{s(t0(γ)+(i1)Tc)2[cos(2πfRF(t0(γ)+(i1)Tc)+ϕ1)
j sin (2πfRF(t0(γ)+(i1)Tc)+ϕ2)]}*hr(t0(γ)+(i1)Tc)+
+ Pn [{Hr(fRF)exp[j2πfRF[t0(γ)+(i1)Tc]]exp(1)}{Hr(fRF)×
×exp[j2πfRF[t0(γ)+(i1)Tc]]exp(2)}]}exp[jWN(k,i)]
E [se(γ)[k]] = E [sFFT(γ)[k]] Ae(γ) [k]exp(jΦe(γ)[k])
me,I(γ) [k]={E[se(γ)[k]]}me,Q(γ)[k]={E[se(γ)[k]]}
(σFFT,(I,Q)(γ)[k])2 = E [(sFFT,(I,Q)(γ)[k])2] {E[sFFT,(I,Q)(γ)[k]]}2
E [sFFT,(I,Q)(γ)[k]] = 12 i=1N+E[sPIN(τ)]f(I,Q)(τ,k,i)hr(t0(γ)+(i1)Tcτ)
E = [sPIN(τ)]=Rλ{s(τ)2+E[nI,2(τ)]+E[nQ,2(τ)]+E[nI,2(τ)]+E[nQ,2(τ)]}
(σFFT,(I,Q)(γ)[k])2 = Rλ2 i=1Nn=1N++{+[sI(τ1)sI(τ2)+SnI,nI,(f1)+sI(τ1)×
× sQ (τ2)SnI,nQ,(f1)+sQ(τ1)sI(τ2)SnQ,nI,(f1)+sQ(τ1)sQ(τ2)SnQ,nQ,(f1)]×
×exp[j2πf1(τ1τ2)]df1+12++[SnI,nI,(f2)SnI,nI,(f3)+SnI,nQ,(f2)SnI,nQ,(f3)+
+ SnQ,nI, (f2)SnQ,nI,(f3)+SnQ,nQ,(f2)SnQ,nQ,(f3)+SnI,nI,(f2)SnI,nI,(f3)+
+ SnI,nQ, (f2)SnI,nQ,(f3)+SnQ,nI,(f2)SnQ,nI,(f3)+SnQ,nQ,(f2)SnQ,nQ,(f3)]×
×exp[j2π(f2+f3)(τ1τ2)]df2df3}g(I,Q)(τ1,τ2,k,i,n)dτ1dτ2
g(I,Q) (τ1,τ2,k,i,n)=f(I,Q)(τ1,k,i)f(I,Q)(τ2,k,n)hr(t0(γ)+(i1)Tcτ1)×
×hr(t0(γ)+(n1)Tcτ2)
(σsASE,(I,Q),y(γ)[k])2 = Rλ2 i=1Nn=1N+++sm(τ1)sp(τ2)Sy(f1)exp[j2πf1(τ1τ2)]df1×
×g(I,Q)(τ1,τ2,k,i,n)dτ1dτ2
sm (τ1)={sI(τ1)ify=1ory=2sQ(τ1)ify=3ory=4 sp(τ2)={sI(τ2)ify=1ory=3sQ(τ2)ify=2ory=4
Sy (f1)={SnI,nI,(f1)ify=1SnI,nQ,(f1)ify=2SnQ,nI,(f1)ify=3SnQ,nQ,(f1)ify=4
(σsASE,(I,Q),y(γ)[k])2 =Rλ2i=1Nn=1N++Sy(f1)Sm(f2)F(I,Q)(f3,k,i)Hr(f4)Sp(f5)×
×F(I,Q)(f6,k,n)Hr(f7)exp[j2πf4(t0(γ)+(i1)Tc)]exp[j2πf7(t0(γ)+(n1)Tc)]×
×exp[j2π(f1+f2+f3f4)τ1]exp[j2π(f1+f5+f6f7)τ2]df1df7dτ1dτ2
FI (f,k,i)=𝓕{f1(f,k,i)}=12[A(k,i)exp(1)+jB(k,i)exp(2)]δ(ffRF)+
+ 12 [A(k,i)exp(1)jB(k,i)exp(2)]δ(f+fRF)
FQ (f,k,i)=𝓕{fQ(f,k,i)}12[B(k,i)exp(1)+jA(k,i)exp(2)]δ(ffRF)+
+ 12 [B(k,i)exp(1)jA(k,i)exp(2)]δ(f+fRF)
(σsASE,(I,Q),y(γ)[k])2 Rλ24 +Sy(f1){Rl,m,(I,Q)(γ)(f1fRF,k)+R2,m,(I,Q)(γ)(f1+fRF,k)}×
×{R2,p,(I,Q)(γ)(f1+fRF,k)+R1,p,(I,Q)(γ)(f1fRF,k)}*df1
R1,m,(I,Q)(γ) (f1fRF,k)=exp(j2πfRFt0(γ))2 H1,m,(I,Q)(γ) (f1+fRF,k)
R2,m,(I,Q)(γ) (f1+fRF,k)=exp(j2πfRFt0(γ))2 H2,m,(I,Q)(γ) (f1fRF,k)
Hi,m,(I,Q)(γ) (f,k)=𝓕 {νm(γ)(t)gi,(I,Q)(γ)(t,k)}gi,(I,Q)(γ)(t,k)=𝓕1{Gi,(I,Q)(γ)(f,k)}
νm(γ) = 𝓕1 {Vm(γ)(f)}=𝓕1{exp(j2πft0(γ))sm(f)}
G(I,Q)(γ) = [G1,(I,Q)(γ)(f,k)G2,(I,Q)(γ)(f,k)] = Hr (f)P(I,Q)T(γ)
T1(γ) (f,k)=exp{j2πN12[fTc+(k1)N]}[sin{2π[fTc+(k1)N]N2}sin{π[fTc+(k1)N]}]
T2(γ) (f,k)=exp{j2πN12[fTc(k1)N]}[sin{2π[fTc(k1)N]N2}sin{π[fTc(k1)N]}]
PI = [exp(1)+exp(2)exp(1)exp(2)exp(1)exp(2)exp(1)+exp(2)]
PQ = j [exp(1)+exp(2)exp(1)+exp(2)exp(1)exp(2)exp(1)exp(2)]
(σASEASE,(I,Q),y(γ)[k])2 = Rλ22 i=1Nn=1N++++Sy(f2)Sy(f3)×
×exp[j2π(f2+f3)(τ1τ2)]df2df3g(I,Q)(τ1,τ2,k,i,n)dτ1dτ2
(σASEASE,(I,Q),y(γ)[k])=18+[Sy(f1)*Sy(f1)]Heq,(I,Q)(γ)(f1,k)2df1
[Heq,1,(I,Q)(γ)(f,k)Heq,2,(I,Q)(γ)(f,k)]=Hr(f)2exp[j2πft0(γ)]P(I,Q)T(γ)
se,I(γ) [k]=Ae(γ)[k]cos(Φe(γ)[k])sFFT,I(γ)[k]+Ae(γ)[k]sin(Φe(γ)[k])sFFT,Q(γ)[k]
se,Q(γ) [k]=Ae(γ)[k]cos(Φe(γ)[k])sFFT,Q(γ)[k]Ae(γ)[k]sin(Φe(γ)[k])sFFT,I(γ)[k]
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.