Abstract

We present a symbol-by-symbol coherent optical receiver, which employs a novel, complex-weighted, decision-aided, maximum-likelihood (CW-DA-ML) carrier phase and frequency offset estimator. The CW-DA-ML carrier estimator uses a CW transversal filter to generate a carrier reference phasor, and the filter weights are automatically adapted on-line by linear regression on the observed signals. A complete modulo-R reduced frequency offset estimation (FOE) range of ±R/2 is achieved, independent of modulation format, where R is the symbol rate. Carrier phase and frequency tracking is achieved rapidly. The acquisition speed of frequency offset in quaternary phase-shift keying (4-PSK) signals is more than 5 times faster than that of differential FOE. A constant penalty of approximately 1 dB at bit-error rate of 10−4 is demonstrated for all frequency offsets in 4-PSK signals with laser-linewidth-symbol-duration product of 8×10−5.

© 2012 OSA

1. Introduction

Traditionally, a major challenge in homodyne coherent receivers is to phase- and frequency-lock the local oscillator (LO) laser to the received optical signal. Recently, digital signal processing (DSP) algorithms have been widely used to perform post-detection carrier synchronization, thereby giving birth to the intradyne coherent receiver featuring an intermediate frequency (IF) which allows for a free-running LO laser [1]. A prevalent DSP based block Mth power carrier estimation (CE) was proposed in [2], but is meant for M-ary phase-shift keying (MPSK) formats only and requires nonlinear operations which increase system latency. A computationally linear, decision-aided, maximum-likelihood (DA-ML) CE with comparable performance to block Mth power CE was derived in [3] and proposed for use in coherent optical communications in [4]. In an intradyne receiver, the frequency offset, ∆f, between LO and transmitter laser can be as large as ±5 GHz [5]. However, block Mth power and DA-ML CE were designed to estimate the carrier phase without considering any frequency offset. Block Mth power and DA-ML CE have limited frequency offset tolerance of 10−3 times the symbol rate at laser-linewidth-symbol-duration product of 6.35×10−5 [2, 6]. Hence, it is imperative for a digital CE to incorporate frequency offset estimation (FOE) capability.

A differential FOE, using the sum of differenced signal phases over an observation interval, was suggested for the block Mth power CE in [7]. It utilizes the nonlinear operations of arctan(·) and Mth power to remove data modulation. Consequently, it is only applicable to MPSK formats and the FOE range is limited to ±R/2M, where R is the symbol rate. With higher order modulation formats, the FOE range shrinks while the complexity of raising to the Mth power increases. This differential FOE has an R/M frequency estimate ambiguity, requiring differential encoding (DE) which leads to a DE-induced performance penalty [8].

Extending upon the DA-ML CE of [3, 4], we propose here a modulation-format-independent symbol-by-symbol receiver employing a novel, complex-weighted (CW), DA-ML CE to estimate the unknown carrier phase and frequency offset. A complete modulo-R reduced frequency offset range of ±R/2 can be estimated for arbitrary modulation formats. CW-DA-ML CE achieves near-ideal FOE with a low training overhead cost and no a priori knowledge of the system statistics. Carrier phase and frequency are acquired rapidly. The FOE for 4-PSK signals is more than 5 times faster than that of the differential FOE of [7].

The paper is organized as follows. In Section 2, the received signal in a phase and polarization diversity intradyne coherent receiver is modeled. Section 3 provides the derivation of the symbol-by-symbol receiver structure consisting of CW-DA-ML CE and data detection (DD). Section 4 examines the CE performance via simulations, and includes a discussion on its complexity and parallel implementation. Section 5 concludes the paper. Throughout this paper, Ε[·]is the expectation, |·| is the modulus operator, boldface characters represent matrices, and δis the Kronecker delta. Superscript , T, and H denote conjugate, transpose, and complex conjugate transpose, respectively. All phase quantities are treated modulo-2π to account for their circular nature and frequency quantities are treated modulo-R, which is justified later in Subsection 4.3.

2. Phase and polarization diversity intradyne coherent receiver

The opto-electronic front end followed by DSP modules in a phase and polarization diversity intradyne coherent receiver is shown in Fig. 1 . The entire received optical field, Er(t), within the receiver bandwidth is mapped into x- and y-polarized, in-phase, I, and quadrature-phase, Q, photocurrents by mixing Er(t) with an LO using two 90° optical hybrids. The LO laser is free running, in contrast to a homodyne coherent receiver where the LO laser needs to be locked to the received optical signal. Each photocurrent branch i(t) is sampled at symbol rate R which amounts to sampling the complex signal at 2R, therefore satisfying the Nyquist theorem for fully preserving the optical signal and noise statistics.

 

Fig. 1 Phase and polarization diversity intradyne coherent receiver. LO: local oscillator, PBS: polarization beam splitter

Download Full Size | PPT Slide | PDF

We assume ideal signal conditioning by the first four DSP blocks, namely, analog-to-digital conversion (ADC), clock recovery and retiming, chromatic dispersion compensation, and polarization demultiplexing and polarization-mode dispersion compensation. The received signal samples rx(k) and ry(k) over the kth symbol interval [kT,(k+1)T) (T = symbol duration) are clock-synchronized with one complex sample per symbol. They are assumed to be free of intersymbol interference (ISI), polarization crosstalk, and nonlinear distortions. The r(k) in one of the polarizations can be shown to be [9]

r(k)=iI(k)+jiQ(k)=χPrPLOej(2πΔfTk+ϕ(k)+θ(k))+χPLOn(k)ej(2πΔfTkθLO(k))+ish
where χ,Pr,PLO, and ϕ(k) denote responsivity of the photodiode, average received signal power, LO power, and data phase in the kth interval, respectively. Here, θ(k) is the combined laser phase noise of the transmitter and LO lasers, whereas θLO(k) is the LO phase noise. Most long-haul optical fiber systems are limited by amplified spontaneous emission (ASE) noise from inline optical amplifiers, which is represented by n'(k) in Eq. (1). The ish is the LO shot noise, assuming PLOPr. Thermal noise is ignored in Eq. (1). On the right-hand side of Eq. (1), the first term is the signal-LO beat term and the second term is the LO-ASE beat noise. Using a pair of balanced photodiodes to detect each photocurrent branch, compared to a single photodiode in [1], eliminates the signal-ASE and ASE-ASE beat noise terms.

A canonical model of the received signal in Eq. (1) can be written as

r(k)=m(k)ej(2πΔfTk+θ(k))+n(k),k=0,1,2,
where m(k)=χPrPLOejϕ(k) is the kth data symbol. The LO-ASE beat noise and LO shot noise are modeled as additive white Gaussian noise (AWGN), and are combined to form n(k)=χPLOn(k)ej(2πΔfTkθLO(k))+ish. The set {n(k)} is a sequence of circularly symmetric, AWGN with mean zero and covariance function Ε[n(k)n*(ki)]=N0δ(i). The carrier phase θ(k) has a Gaussian random-walk model of θ(k)=η(k)+θ(k1), where {η(k)}is a set of independent, identically distributed, Gaussian random variables with mean zero and variance σp2=2πΔνT [9]. Here, ∆ν is the combined 3-dB linewidth of the transmitter and LO lasers. The carrier phase θ(k) has a temporal correlation characterized by Ε[θ(i)θ(j)]=σp2min[i,j]. Signal-to-noise ratio (SNR) per bit is defined as γb=Ε[|m(k)|2]/N0log2M. In an intradyne receiver, Er(t) is downconverted to an IF of ∆f as seen in Eq. (2). Optical frequency bands around fLO+Δf and fLOΔf map to the same IF, where fLO is the LO frequency. In order to avoid crosstalk in a dense wavelength-division-multiplexing system and to avoid excess optical noise from unwanted image bands, an optical filter of a single channel bandwidth is required before the intradyne receiver.

3. Carrier estimation and data detection

In this paper, we consider the detection of an uncoded data symbol sequence {m(k)} transmitted over an AWGN channel with unknown carrier phase and frequency offset, as modeled in Eq. (2). Because the data is uncoded, the symbols of {m(k)} are independent, and each symbol can assume with equal probability any point Si, i = 0,…,M-1, in the signal constellation. The elements of {θ(k)} are correlated, whereas the assumption of an AWGN channel with no ISI makes the elements of {n(k)} independent. The elements of {r(k)} are rendered independent when conditioned on given values of uncoded data symbol sequence {m(k)}, carrier phase process {θ(k)}, and frequency offset, ∆f. Hence, each data symbol m(k) in the sequence {r(k)} will be detected individually, i.e., symbol-by-symbol with minimum symbol error probability. We assume mutual statistical independence among {m(k)}, {θ(k)}, and ∆f, which leads to the separation of the CE problem from the DD problem [10]. Therefore, at high SNR, the optimum symbol-by-symbol receiver structure consists of a CE followed by a coherent DD, as illustrated in Fig. 1 [10].

We develop our new CE algorithm in Subsection 3.1 for carrier phase and frequency offset estimation. The estimated phase and frequency are treated as if they were the true values of the carrier and are used in (partially) coherent DD in Subsection 3.2. Separate symbol-by-symbol receiver structures, consisting of CE and DD, can be applied independently on each polarization channel as shown in Fig. 1. All equations and quantities expressed hereafter are thus meant for one polarization channel.

3.1 Complex-weighted DA-ML carrier estimation

In DA-ML CE [3, 4], the carrier phase is estimated using a complex reference phasor (RP), where the argument of the RP is the ML phase estimate of the carrier. It is assumed that the carrier phase process {θ(k)} is fluctuating slowly compared to the symbol rate such that it can be approximated to be piecewise constant over intervals longer than LT, where L is the estimator filter length. The RP V(k+1) without considering any frequency offset is computed in [3, 4] using the immediate past L received signals as

V(k+1)=C(k)l=1Lr(kl+1)m^*(kl+1)
where m^(k) is the decision on kth symbol from DD and C1(k)=l=1L|m^(kl+1)|2 is a factor to normalize the RP's magnitude in the event of a non-constant-energy signal constellation.

Considering the presence of an unknown carrier frequency offset, we propose here a new RP V'(k+1) that is formed by

V(k+1)=C(k)l=1Lwl(k)r(kl+1)m^*(kl+1)
where each wl(k) is a complex weight. Momentarily ignoring the carrier phase θ(k) and AWGN n(k), the filter input terms with a carrier frequency offset appear as {r(kl+1)m^*(kl+1)=|m(kl+1)|2ejΔω(kl+1),1lL} where we have assumed m^(k)=m(k), i.e., all receiver decisions are correct, and Δω=2πΔfT is the angular frequency offset. The complex weight wl(k) is designed to rotate r(kl+1)m^*(kl+1) by Δωl such that all filter input terms are aligned with identical angular frequency offset of Δω(k+1). This leads to the RP V'(k+1) being angular-frequency-offset aligned to the next received signal r(k+1)=m(k+1)ejΔω(k+1), which enables an enhanced signal demodulation. The filter in Eq. (4) operates without the knowledge of SNR, carrier phase, and frequency offset statistics.

Ideally, one should have wl(k)=ejΔωl, but ∆ω is unknown in practice. Thus, we propose to choose the complex weights automatically and adaptively at each time k based on the observations {r(l),0lk} to minimize the cost function J(k) given by [11]

J(k)=l=1k|e(l)|2=l=1k|r(l)m^(l)C(l1)wT(k)y(l1)|2
where w(k)=[w1(k)wL(k)]T is the L-by-1 filter-weight vector and y(k)=[r(k)m^*(k)r(kL+1)m^*(kL+1)]T is the L-by-1 filter-input vector at time k. The error e(l) is the difference between the desired response r(l)/m^(l) and the RP output of Eq. (4) at time l1 using the latest set of filter coefficients w(k). Minimization of the cost function J(k) forces RP V(l) to track the normalized term r(l)/m^(l), and thus forces wl(k) to track ejΔωl. Adaptation of filter weights using a least-squares criterion, as opposed to a mean-square-error criterion J˜(k)=Ε[|e(k)|2], requires no statistical information about the AWGN, carrier phase, or frequency offset.

Since J(k) is a real valued function of w(k) and w*(k), we solve the equation: J(k)/w*(k)=0, for w(k) using the complex vector differentiation identity: [awHb]/w*=b, where a is a scalar and b is a vector. The optimum w^(k) is obtained as the solution of a least-squares normal equation

w^(k)=Φ1(k)z(k),k1
where Φ(k)=l=1kC2(l1)y*(l1)yT(l1) is the L-by-L time-average autocorrelation matrix and z(k)=l=1kC(l1)y*(l1)r(l)/m^(l) is the L-by-1 time-average cross-correlation vector. The Φ1(k) matrix and the least-squares w^(k) at each time k can be obtained recursively using the matrix inversion lemma [12]. Therefore, CW-DA-ML CE avoids any matrix inversion. The optimum w^(k) can respond to changing channel conditions, as it depends on the observations {r(k)}.

3.2 Symbol-by-symbol data detection

Treating V(k) as the true phasor ej(Δωk+θ(k)), we multiply r(k) by V*(k) to compensate for the carrier phase and frequency offset. The derotated r(k) is plugged into the ML minimum-distance DD for an AWGN channel given by [12]

m^(k)=argmin0iM1|r(k)V*(k)Si|.
Assuming unconstrained symbol energy, the DD is simplified by expanding the equation, dropping terms that do not depend on the trial data symbol Si, and rearranging them to yield
m^(k)=argmin0iM1(|r(k)V*(k)|2+|Si|2Re[2r(k)V*(k)Si*])=argmin0iM1(|Si|2Re[2r(k)V*(k)Si*])=argmax0iM1Re[r(k)V*(k)Si*12|Si|2].
Finally, our symbol-by-symbol receiver performs (partially) coherent DD by forming the decision statistic as in Eq. (8).

Our symbol-by-symbol receiver algorithm employing the novel CW-DA-ML CE with recursive updating of Φ1(k) matrix is outlined in Table 1 . The Φ1(k) matrix is Hermitian symmetric, thus only the upper triangle of Φ1(k) needs to be computed whereas the lower triangle is obtained through diagonal reflection as signified by the Τri{·} operator in Table 1. In operating the filter of Eq. (4), we set V(0)=1 and w^(0)=[100]T to give a maximum gain of one on the first filter input r(0)m^*(0). An initial preamble of N known symbols is used to enable w^(k) to settle to a steady state and for V(k) to acquire tracking of the phasor ej(Δωk+θ(k)). The filter operates in decision-directed mode subsequently.

Tables Icon

Table 1. Symbol-by-symbol receiver employing CW-DA-ML CE

4. Numerical results and discussion

In all subsequent Monte-Carlo simulations, a single-polarization received signal modeled after Eq. (2) is used with a fixed bit rate of 50 Gbit/s for all modulation formats. A filter length of L = 5 is maintained for all CE algorithms tested. The modulo-2π reduced ∆ω is assumed to have a probability density function given by p(∆ω) = 1/2π for ∆ω[-π, π), where ∆ω is random but time invariant. The models of the carrier phase,θ(k), and angular frequency offset, ∆ω, are not known to the receiver.

4.1 Operation of CW-DA-ML CE

In this subsection, we analyze the statistical performance of CW-DA-ML CE by simulating carrier recovery for a 4-PSK signal with a preamble of N = 50 and a laser linewidth of 2 MHz.

Figure 2 plots the acquisition of the argument of the complex weights, {arg(w^l(k)),1lL}, from Eq. (6) at an γb of 10 dB and a frequency offset of 3 GHz. The arg(w^l(k)) converges quickly to the actual value of lω. Hence, wl(k) in Eq. (4) effectively removes the angular frequency offset difference between consecutive filter input terms {r(kl+1)m^*(kl+1),1lL} and helps to produce a more accurate RP compared to that produced by Eq. (3).

 

Fig. 2 Adaptation of arg(w^l(k)) from Eq. (6) in one sample simulation run.

Download Full Size | PPT Slide | PDF

The mean-square error (MSE) learning curve of CW-DA-ML CE algorithm given by J'(k)=Ε[|r(k)/m^(k)V(k)|2] is empirically evaluated. The ensemble-average squared error curves averaged over 104 runs at frequency offsets of 0, 0.03, and 3 GHz are plotted for SNR values of 7, 10, and 13 dB in Fig. 3 . Ideal decision feedback is assumed, i.e., m^(k)=m(k) for all k. Using Eq. (2), the MSE J'(k) can be expanded into

J(k)=Ε[|n(k)m(k)|2+|ej(Δωk+θ(k))V(k)|2+(ej(Δωk+θ(k))V(k))(n(k)m(k))*+(ej(Δωk+θ(k))V(k))*(n(k)m(k))].
The V(k) is independent of n(k) and m(k), as V(k) depends only on {n(l),m(l),l<k}. Moreover, ∆f, θ(k), n(k), and m(k) are mutually independent . Thus, the third and fourth terms on the right hand side of Eq. (9) amount to zero because Ε[n(k)/m(k)]=0. Equation (9) then reduces to J(k)=Jmin(k)+Jex(k), where Jmin(k)=Ε[|n(k)/m(k)|2] is the minimum achievable MSE and Jex(k)=Ε[|ej(Δωk+θ(k))V(k)|2] is the excess MSE due to the MSE of the RP. It can be shown that Jmin=1/γblog2Mfor 4-PSK signals. Figure 3 shows the values of Jmin, and also Jex which was obtained by averaging from time k = 100 to k = 200 the average of |ej(Δωk+θ(k))V(k)|2 over 104 runs at each time point k. A Jex value less than 2.4×10−2 demonstrates the good tracking of ej(Δωk+θ(k)) by V(k), for SNR values above 7 dB. The excess MSE Jex decreases with increasing SNR but is indifferent to varying frequency offset, attesting that FOE by CW-DA-ML CE algorithm is unbiased to the frequency offset present.

 

Fig. 3 Ensemble-average squared error curves with three different ∆f and SNR values.

Download Full Size | PPT Slide | PDF

CW-DA-ML CE algorithm converges quickly yielding rapid carrier phase and frequency tracking. We note that a short preamble of approximately twice the filter length, N ≈2L, is sufficient to aid the complex weight w^l(k) in acquiring tracking of ejΔωl, thus keeping the training overhead cost low. The MSE of CW-DA-ML CE for 4-PSK signals at SNR values above 7 dB settles to a steady state, implying acquisition of the frequency offset, within 100 observed samples. Whereas the differential FOE for 4-PSK signals in [7] used a minimum observation sample size of 500 for acquisition of the frequency offset. Hence, CW-DA-ML FOE is more than 5 times faster than that of the differential FOE of [7] for 4-PSK signals at SNR values above 7 dB.

4.2 Bit-error rate performance with CW-DA-ML CE

Bit-error rate (BER) curves of block Mth power [2], DA-ML [3, 4], and CW-DA-ML CE are plotted in Fig. 4 for 4- and 8-PSK, as a function of SNR per bit, γb, in the presence of different frequency offsets. Filter length and laser linewidth was set to 5 and 0, respectively. A preamble of N = 10 known symbols is used for DA-ML and CW-DA-ML CE. DE according to [8] is employed to resolve phase ambiguity in block Mth power CE, and to counter runaways due to decision errors in DA-ML and CW-DA-ML CE.

 

Fig. 4 BER curves of block Mth power, DA-ML, and CW-DA-ML CE with different frequency offsets for (a) 4-PSK, and (b) 8-PSK.

Download Full Size | PPT Slide | PDF

It is evident that block Mth power and DA-ML CE are intolerant to frequency offsets. For instance, in block Mth power CE of 4-PSK signals, as the frequency offset increases from 0 to 100 and 200 MHz, the γb penalty rises from 1.13 to 1.33 and 1.80 dB, respectively. Likewise for DA-ML CE, the γb penalty rises from 0.89 to 1.30 and 2.14 dB. SNR per bit penalties are referenced to the perfect coherent (ASE noise limited) receiver at BER = 10−4. Moreover, block Mth power and DA-ML CE become increasingly sensitive to frequency offset with increasing modulation order. For a fixed frequency offset of 100 MHz, as modulation format increases from 4- to 8-PSK, block Mth power CE incurs an increasing γb penalty of 1.33 and 3.18 dB, respectively. Likewise, DA-ML CE incurs an increasing γb penalty of 1.30 and 3.63 dB. Frequency offset intolerance of block Mth power and DA-ML CE originates from the violation of their assumption that {r(l),kL+1lk} have identical angular frequency offsets.

On the contrary, the γb penalty for CW-DA-ML CE remains unchanged at 0.89 and 1.13 dB regardless of the frequency offset present in 4- and 8-PSK, respectively, at BER = 10−4. Regardless of ∆f, BER of CW-DA-ML CE over all tested SNR values closely replicate that of DA-ML CE at zero frequency offset. Hence, CW-DA-ML CE approaches ideal FOE, without requiring any statistical knowledge of AWGN, carrier phase, or frequency offset.

Figure 4 also shows CW-DA-ML CE at ∆f = 100 MHz with ideal decision feedback, i.e., m^(k)=m(k) for all k. In the range of SNR shown, the performance difference between actual and ideal decision feedback is not noticeable for 4-PSK, and is slightly more pronounced at the lower SNR region for 8-PSK. The effect of decision errors on the receiver performance becomes negligible for BER's below 10−3.

4.3 FOE range of CW-DA-ML CE

SNR per bit, γb, penalty of block Mth power CE [2], DA-ML CE [3, 4], differential FOE [7], and CW-DA-ML CE are plotted in Fig. 5 for 4- and 8-PSK, as frequency offset is varied between ±R/2, in the presence of different laser linewidths, ∆ν. Any integer-multiple-of-R frequency offset, ∆f, translates into the same modulo-2π reduced ∆ω, thus having the same rotational effect on r(k). Hence, it is sufficient to consider only the primary modulo-R reduced frequency offset range ±R/2. The differential FOE utilizes a subsequent block Mth power CE to correct any residual carrier phase, DE to overcome its R/M frequency estimate ambiguity, and a sample size of 500 to acquire the frequency offset, following that of [7].

 

Fig. 5 FOE range of block Mth power CE, DA-ML CE, differential FOE, and CW-DA-ML CE for (a) 4-PSK, and (b) 8-PSK. Insets show enlarged sections of the frequency offset axis.

Download Full Size | PPT Slide | PDF

CE techniques, such as block Mth power and differential FOE, which erase data modulation by raising the received signal to the Mth power, and extract the phase using the arctan(·) operator are inherently disadvantaged in three ways. First, they are restricted to MPSK signals. Second, their FOE range is limited, modulation-format dependent, and decreases with increasing modulation order. For a fixed bit rate of Rb, block Mth power CE has a maximum allowable frequency offset of ±Rb/2LMlog2(M) in the absence of carrier phase variations [7], and the differential FOE has an FOE range of ±Rb/2Mlog2(M) [13]. Third, they require phase unwrapping in estimating the carrier phase and frequency offset, which is difficult to perform accurately at low SNR and is prone to the highly nonlinear cycle slipping phenomenon [14].

On the other hand, our CW-DA-ML CE achieves a complete modulo-R reduced, modulation-format-independent, FOE range of ±R/2 as witnessed in Fig. 5. This is attributed to the use of the RP, V(k), having an unambiguous phase tracking range of [0, 2π). The use of the RP also eliminates the need for phase unwrapping in carrier phase and frequency offset estimation. The circular-2π nature of the RP results in frequency-estimate duplicity of R, but this duplicity does not call for pilot assistance or DE to ensure correct DD. From Fig. 5(a), it is notable that a penalty of only approximately 1-dB is incurred at BER = 10−4 for a 50 Gbit/s 4-PSK signal having a laser linewidth of 2 MHz, regardless of the frequency offset present.

4.4 Complexity and parallel implementation of CW-DA-ML CE

The computational load, defined as the number of operations required to estimate ej(Δωk+θ(k)) per symbol, of the CW-DA-ML CE algorithm presented in Table 1 is O(L2) real multiplications and additions. However, by implementing Eq. (4) in a lattice structure, which is not presented here for brevity, the complexity reduces to simply O(L) real multiplications and additions [12]. Differential FOE (inclusive of subsequent block Mth power CE) requires O(log2M) real multiplications and additions, and (2L + 1)/L accesses to a read-only-memory to map arctan(·) operations. The 1/L factor arises from sharing of the common arctan(·) operation among L symbols during the block Mth power processing. The higher complexity of CW-DA-ML CE compared to differential FOE is traded off with CW-DA-ML CE’s wider and faster frequency estimation, and its applicability to arbitrary modulation formats. The update of w(k) need be performed only until the frequency offset is acquired, say, until the 100th iteration for 4-PSK signals according to Fig. 3. Thereafter, the same weight vector w(k) can be used in Eq. (4) at each time point k, i.e., skip steps 4 to 8 in Table 1. Then, the complexity of CW-DA-ML CE becomes simply 4L + 6 and 4L real multiplications and additions, respectively, per symbol.

Our CW-DA-ML CE is a D-symbol forward linear predictor, where the filter predicts ej(Δω(k+D)+θ(k+D)) as a weighted linear combination of past observed samples {r(l),lk} and symbol decisions {m^(l),lk} which are fed back. All of the above simulations used an ideal feedback delay of D = 1, which is the minimum in a feedback system. Practical coherent optical receivers supporting 100 Gb/s or higher data rates use parallel and pipelined processing to overcome the limited complementary metal-oxide-semiconductor processor speed of several Gb/s. In parallel processing, the effect of the observed sample on the feedback value is delayed by D > 1 [15]. Parallel implementation of CW-DA-ML CE, with a constant frequency offset, zero laser phase noise, and a feedback delay of D-symbols, is expected to incur marginal performance deterioration as it can estimate D-symbols ahead the linearly increasing angular frequency offset ∆ω(k + D). However, in the additional presence of laser phase noise, performance penalty originating from feedback delay is expected as pointed out in [15]. This is due to the difficulty in estimating the accumulated Gaussian random-walk of the carrier phase over the previous D-symbols, for which the estimated carrier phase information is unavailable because of feedback delay. A detailed implementation and performance analysis of CW-DA-ML CE in parallel processing is left as part of our future work.

5. Conclusions

If the frequency-offset-symbol-duration product is small, separate frequency offset estimation (FOE) may not be necessary as block Mth power and decision-aided, maximum-likelihood carrier estimation (DA-ML CE) can cope with small frequency offsets. To combat larger frequency offsets, an FOE is necessary. Hence, we proposed a symbol-by-symbol receiver employing a novel complex-weighted (CW) DA-ML CE for optical channels with both carrier phase and frequency offset uncertainty. A near-ideal FOE over a complete modulo-R reduced frequency range of ±R/2 (R = symbol rate) is demonstrated such that the bit-error rate performance is only limited by the carrier phase variance. CW-DA-ML CE is applicable to both M-ary phase-shift keying (MPSK) and quadrature amplitude modulation formats. The FOE for 4-PSK signals is more than 5 times faster than that of the differential FOE at signal-to-noise ratio values above 7 dB.

Acknowledgments

The financial support of the Singapore MoE AcRF Tier 2 Grant MOE2010-T2-1-101 is gratefully acknowledged.

References and links

1. F. Derr, “Coherent optical QPSK intradyne system: concept and digital receiver realization,” J. Lightwave Technol. 10(9), 1290–1296 (1992). [CrossRef]  

2. D. S. Ly-Gagnon, S. Tsukamoto, K. Katoh, and K. Kikuchi, “Coherent detection of optical quadrature phase-shift keying signals with carrier phase estimation,” J. Lightwave Technol. 24(1), 12–21 (2006). [CrossRef]  

3. P. Y. Kam, “Maximum likelihood carrier phase recovery for linear suppressed-carrier digital data modulations,” IEEE Trans. Commun. 34(6), 522–527 (1986). [CrossRef]  

4. S. Zhang, P. Y. Kam, J. Chen, and C. Yu, “Decision-aided maximum likelihood detection in coherent optical phase-shift-keying system,” Opt. Express 17(2), 703–715 (2009). [CrossRef]   [PubMed]  

5. “Integrable Tunable Laser Assembly Multi Source Agreement,” OIF-ITLA-MSA-01.1 (Optical Internetworking Forum, Los Angeles, 2005).

6. A. Meiyappan, P. Y. Kam, and H. Kim, “Performance of decision-aided maximum-likelihood carrier phase estimation with frequency offset,” in Proc. OFC/NFOEC, Los Angeles, CA, 2012, paper OTu2G.6.

7. A. Leven, N. Kaneda, U.-V. Koc, and Y.-K. Chen, “Frequency estimation in intradyne reception,” IEEE Photon. Technol. Lett. 19(6), 366–368 (2007). [CrossRef]  

8. W. J. Weber, “Differential encoding for multiple amplitude and phase shift keying systems,” IEEE Trans. Commun. 26(3), 385–391 (1978). [CrossRef]  

9. K.-P. Ho, Phase-modulated Optical Communication (Springer, 2005).

10. P. Y. Kam, S. S. Ng, and T. S. Ng, “Optimum symbol-by-symbol detection of uncoded digital data over the Gaussian channel with unknown carrier phase,” IEEE Trans. Commun. 42(8), 2543–2552 (1994). [CrossRef]  

11. P. Y. Kam, K. H. Chua, and X. Yu, “Adaptive symbol-by-symbol reception of MPSK on the Gaussian channel with unknown carrier phase characteristics,” IEEE Trans. Commun. 46(10), 1275–1279 (1998). [CrossRef]  

12. J. G. Proakis, Digital Communications (McGraw-Hill, 2008).

13. U. Mengali and A. N. D' Andrea, Synchronization Techniques for Digital Receivers (Plenum Press, 1997).

14. H. Meyr, M. Moeneclaey, and S. Fechtel, Digital Communication Receivers (John Wiley, 1997).

15. T. Pfau, S. Hoffmann, and R. Noe, “Hardware-efficient coherent digital receiver concept with feedforward carrier recovery for M-QAM constellations,” J. Lightwave Technol. 27(8), 989–999 (2009). [CrossRef]  

References

  • View by:
  • |
  • |
  • |

  1. F. Derr, “Coherent optical QPSK intradyne system: concept and digital receiver realization,” J. Lightwave Technol.10(9), 1290–1296 (1992).
    [CrossRef]
  2. D. S. Ly-Gagnon, S. Tsukamoto, K. Katoh, and K. Kikuchi, “Coherent detection of optical quadrature phase-shift keying signals with carrier phase estimation,” J. Lightwave Technol.24(1), 12–21 (2006).
    [CrossRef]
  3. P. Y. Kam, “Maximum likelihood carrier phase recovery for linear suppressed-carrier digital data modulations,” IEEE Trans. Commun.34(6), 522–527 (1986).
    [CrossRef]
  4. S. Zhang, P. Y. Kam, J. Chen, and C. Yu, “Decision-aided maximum likelihood detection in coherent optical phase-shift-keying system,” Opt. Express17(2), 703–715 (2009).
    [CrossRef] [PubMed]
  5. “Integrable Tunable Laser Assembly Multi Source Agreement,” OIF-ITLA-MSA-01.1 (Optical Internetworking Forum, Los Angeles, 2005).
  6. A. Meiyappan, P. Y. Kam, and H. Kim, “Performance of decision-aided maximum-likelihood carrier phase estimation with frequency offset,” in Proc. OFC/NFOEC, Los Angeles, CA, 2012, paper OTu2G.6.
  7. A. Leven, N. Kaneda, U.-V. Koc, and Y.-K. Chen, “Frequency estimation in intradyne reception,” IEEE Photon. Technol. Lett.19(6), 366–368 (2007).
    [CrossRef]
  8. W. J. Weber, “Differential encoding for multiple amplitude and phase shift keying systems,” IEEE Trans. Commun.26(3), 385–391 (1978).
    [CrossRef]
  9. K.-P. Ho, Phase-modulated Optical Communication (Springer, 2005).
  10. P. Y. Kam, S. S. Ng, and T. S. Ng, “Optimum symbol-by-symbol detection of uncoded digital data over the Gaussian channel with unknown carrier phase,” IEEE Trans. Commun.42(8), 2543–2552 (1994).
    [CrossRef]
  11. P. Y. Kam, K. H. Chua, and X. Yu, “Adaptive symbol-by-symbol reception of MPSK on the Gaussian channel with unknown carrier phase characteristics,” IEEE Trans. Commun.46(10), 1275–1279 (1998).
    [CrossRef]
  12. J. G. Proakis, Digital Communications (McGraw-Hill, 2008).
  13. U. Mengali and A. N. D' Andrea, Synchronization Techniques for Digital Receivers (Plenum Press, 1997).
  14. H. Meyr, M. Moeneclaey, and S. Fechtel, Digital Communication Receivers (John Wiley, 1997).
  15. T. Pfau, S. Hoffmann, and R. Noe, “Hardware-efficient coherent digital receiver concept with feedforward carrier recovery for M-QAM constellations,” J. Lightwave Technol.27(8), 989–999 (2009).
    [CrossRef]

2009 (2)

2007 (1)

A. Leven, N. Kaneda, U.-V. Koc, and Y.-K. Chen, “Frequency estimation in intradyne reception,” IEEE Photon. Technol. Lett.19(6), 366–368 (2007).
[CrossRef]

2006 (1)

1998 (1)

P. Y. Kam, K. H. Chua, and X. Yu, “Adaptive symbol-by-symbol reception of MPSK on the Gaussian channel with unknown carrier phase characteristics,” IEEE Trans. Commun.46(10), 1275–1279 (1998).
[CrossRef]

1994 (1)

P. Y. Kam, S. S. Ng, and T. S. Ng, “Optimum symbol-by-symbol detection of uncoded digital data over the Gaussian channel with unknown carrier phase,” IEEE Trans. Commun.42(8), 2543–2552 (1994).
[CrossRef]

1992 (1)

F. Derr, “Coherent optical QPSK intradyne system: concept and digital receiver realization,” J. Lightwave Technol.10(9), 1290–1296 (1992).
[CrossRef]

1986 (1)

P. Y. Kam, “Maximum likelihood carrier phase recovery for linear suppressed-carrier digital data modulations,” IEEE Trans. Commun.34(6), 522–527 (1986).
[CrossRef]

1978 (1)

W. J. Weber, “Differential encoding for multiple amplitude and phase shift keying systems,” IEEE Trans. Commun.26(3), 385–391 (1978).
[CrossRef]

Chen, J.

Chen, Y.-K.

A. Leven, N. Kaneda, U.-V. Koc, and Y.-K. Chen, “Frequency estimation in intradyne reception,” IEEE Photon. Technol. Lett.19(6), 366–368 (2007).
[CrossRef]

Chua, K. H.

P. Y. Kam, K. H. Chua, and X. Yu, “Adaptive symbol-by-symbol reception of MPSK on the Gaussian channel with unknown carrier phase characteristics,” IEEE Trans. Commun.46(10), 1275–1279 (1998).
[CrossRef]

Derr, F.

F. Derr, “Coherent optical QPSK intradyne system: concept and digital receiver realization,” J. Lightwave Technol.10(9), 1290–1296 (1992).
[CrossRef]

Hoffmann, S.

Kam, P. Y.

S. Zhang, P. Y. Kam, J. Chen, and C. Yu, “Decision-aided maximum likelihood detection in coherent optical phase-shift-keying system,” Opt. Express17(2), 703–715 (2009).
[CrossRef] [PubMed]

P. Y. Kam, K. H. Chua, and X. Yu, “Adaptive symbol-by-symbol reception of MPSK on the Gaussian channel with unknown carrier phase characteristics,” IEEE Trans. Commun.46(10), 1275–1279 (1998).
[CrossRef]

P. Y. Kam, S. S. Ng, and T. S. Ng, “Optimum symbol-by-symbol detection of uncoded digital data over the Gaussian channel with unknown carrier phase,” IEEE Trans. Commun.42(8), 2543–2552 (1994).
[CrossRef]

P. Y. Kam, “Maximum likelihood carrier phase recovery for linear suppressed-carrier digital data modulations,” IEEE Trans. Commun.34(6), 522–527 (1986).
[CrossRef]

Kaneda, N.

A. Leven, N. Kaneda, U.-V. Koc, and Y.-K. Chen, “Frequency estimation in intradyne reception,” IEEE Photon. Technol. Lett.19(6), 366–368 (2007).
[CrossRef]

Katoh, K.

Kikuchi, K.

Koc, U.-V.

A. Leven, N. Kaneda, U.-V. Koc, and Y.-K. Chen, “Frequency estimation in intradyne reception,” IEEE Photon. Technol. Lett.19(6), 366–368 (2007).
[CrossRef]

Leven, A.

A. Leven, N. Kaneda, U.-V. Koc, and Y.-K. Chen, “Frequency estimation in intradyne reception,” IEEE Photon. Technol. Lett.19(6), 366–368 (2007).
[CrossRef]

Ly-Gagnon, D. S.

Ng, S. S.

P. Y. Kam, S. S. Ng, and T. S. Ng, “Optimum symbol-by-symbol detection of uncoded digital data over the Gaussian channel with unknown carrier phase,” IEEE Trans. Commun.42(8), 2543–2552 (1994).
[CrossRef]

Ng, T. S.

P. Y. Kam, S. S. Ng, and T. S. Ng, “Optimum symbol-by-symbol detection of uncoded digital data over the Gaussian channel with unknown carrier phase,” IEEE Trans. Commun.42(8), 2543–2552 (1994).
[CrossRef]

Noe, R.

Pfau, T.

Tsukamoto, S.

Weber, W. J.

W. J. Weber, “Differential encoding for multiple amplitude and phase shift keying systems,” IEEE Trans. Commun.26(3), 385–391 (1978).
[CrossRef]

Yu, C.

Yu, X.

P. Y. Kam, K. H. Chua, and X. Yu, “Adaptive symbol-by-symbol reception of MPSK on the Gaussian channel with unknown carrier phase characteristics,” IEEE Trans. Commun.46(10), 1275–1279 (1998).
[CrossRef]

Zhang, S.

IEEE Photon. Technol. Lett. (1)

A. Leven, N. Kaneda, U.-V. Koc, and Y.-K. Chen, “Frequency estimation in intradyne reception,” IEEE Photon. Technol. Lett.19(6), 366–368 (2007).
[CrossRef]

IEEE Trans. Commun. (4)

W. J. Weber, “Differential encoding for multiple amplitude and phase shift keying systems,” IEEE Trans. Commun.26(3), 385–391 (1978).
[CrossRef]

P. Y. Kam, “Maximum likelihood carrier phase recovery for linear suppressed-carrier digital data modulations,” IEEE Trans. Commun.34(6), 522–527 (1986).
[CrossRef]

P. Y. Kam, S. S. Ng, and T. S. Ng, “Optimum symbol-by-symbol detection of uncoded digital data over the Gaussian channel with unknown carrier phase,” IEEE Trans. Commun.42(8), 2543–2552 (1994).
[CrossRef]

P. Y. Kam, K. H. Chua, and X. Yu, “Adaptive symbol-by-symbol reception of MPSK on the Gaussian channel with unknown carrier phase characteristics,” IEEE Trans. Commun.46(10), 1275–1279 (1998).
[CrossRef]

J. Lightwave Technol. (3)

Opt. Express (1)

Other (6)

“Integrable Tunable Laser Assembly Multi Source Agreement,” OIF-ITLA-MSA-01.1 (Optical Internetworking Forum, Los Angeles, 2005).

A. Meiyappan, P. Y. Kam, and H. Kim, “Performance of decision-aided maximum-likelihood carrier phase estimation with frequency offset,” in Proc. OFC/NFOEC, Los Angeles, CA, 2012, paper OTu2G.6.

K.-P. Ho, Phase-modulated Optical Communication (Springer, 2005).

J. G. Proakis, Digital Communications (McGraw-Hill, 2008).

U. Mengali and A. N. D' Andrea, Synchronization Techniques for Digital Receivers (Plenum Press, 1997).

H. Meyr, M. Moeneclaey, and S. Fechtel, Digital Communication Receivers (John Wiley, 1997).

Cited By

OSA participates in CrossRef's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1
Fig. 1

Phase and polarization diversity intradyne coherent receiver. LO: local oscillator, PBS: polarization beam splitter

Fig. 2
Fig. 2

Adaptation of arg( w ^ l ( k ) ) from Eq. (6) in one sample simulation run.

Fig. 3
Fig. 3

Ensemble-average squared error curves with three different ∆f and SNR values.

Fig. 4
Fig. 4

BER curves of block Mth power, DA-ML, and CW-DA-ML CE with different frequency offsets for (a) 4-PSK, and (b) 8-PSK.

Fig. 5
Fig. 5

FOE range of block Mth power CE, DA-ML CE, differential FOE, and CW-DA-ML CE for (a) 4-PSK, and (b) 8-PSK. Insets show enlarged sections of the frequency offset axis.

Tables (1)

Tables Icon

Table 1 Symbol-by-symbol receiver employing CW-DA-ML CE

Equations (9)

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

r( k )= i I ( k )+j i Q ( k ) =χ P r P LO e j( 2πΔfTk+ϕ( k )+θ( k ) ) +χ P LO n ( k ) e j( 2πΔfTk θ LO ( k ) ) + i sh
r( k )=m( k ) e j( 2πΔfTk+θ( k ) ) +n( k ),k=0,1,2,
V( k+1 )=C( k ) l=1 L r( kl+1 ) m ^ *( kl+1 )
V ( k+1 )=C( k ) l=1 L w l ( k )r( kl+1 ) m ^ *( kl+1 )
J( k )= l=1 k | e( l ) | 2 = l=1 k | r( l ) m ^ ( l ) C( l1 ) w T ( k )y( l1 ) | 2
w ^ ( k )= Φ 1 ( k )z( k ),k1
m ^ ( k )= argmin 0iM1 | r( k ) V *( k ) S i |.
m ^ ( k )= argmin 0iM1 ( | r( k ) V *( k ) | 2 + | S i | 2 Re[ 2r( k ) V *( k ) S i * ] ) = argmin 0iM1 ( | S i | 2 Re[ 2r( k ) V *( k ) S i * ] ) = argmax 0iM1 Re[ r( k ) V *( k ) S i * 1 2 | S i | 2 ].
J ( k )=Ε [ | n( k ) m( k ) | 2 + | e j( Δωk+θ( k ) ) V ( k ) | 2 +( e j( Δωk+θ( k ) ) V ( k ) ) ( n( k ) m( k ) ) * + ( e j( Δωk+θ( k ) ) V ( k ) ) * ( n( k ) m( k ) ) ].

Metrics