The disclosure relates generally to communication systems, and more specifically, to cycle slip compensation in a coherent receiver.
The most recent generation of high-speed optical transport network systems has widely adopted receiver technologies with electronic dispersion compensation (EDC). These receivers often use phase encoded modulation formats such as phase shift keying (PSK) or quadrature PSK (QPSK). Carrier phase recovery (CPR) is a key function of coherent optical receivers. CPR algorithms compensate for effects such as laser phase noise and carrier frequency fluctuations. However, since various modulation schemes such as QPSK have rotational symmetry, errors in the carrier phase estimation may cause cycle slips. After a cycle slip occurs, all detected symbols are erroneous and they generally cannot be corrected by forward error correction (FEC) codes. To combat this catastrophic effect, differential modulation is typically used. In differential modulation schemes the information is transmitted as the phase difference between two consecutive symbols. Therefore, the effects of a cycle slip do not translate into catastrophic bit errors when differential modulation is used. However, differential modulation is prone to introduce signal-to-noise ratio (SNR) penalty compared to non-differential schemes, and is therefore an undesirable solution to the cycle slip problem.
A coherent receiver detects and corrects for cycle slips in a received data stream. A sequence of frames are received with each frame comprising a data symbol block and a pilot symbol block having one or more pilot symbols. A comparison is performed between a sequence of received pilot symbol blocks from the sequence of frames and an expected sequence of pilot symbol blocks. A phase rotation of the sequence of received pilot symbol blocks relative to the expected sequence of pilot symbol blocks is estimated based on the comparison. The estimated phase rotation is applied to a target data symbol block in a target frame of the sequence of frames to generate a compensated data symbol block. The compensated data symbol block is outputted.
In a second embodiment a sequence of frames is received with each frame comprising a data symbol block and a pilot symbol block. The pilot symbol block comprises one or more pilot symbols. The sequence of frames including at least a target frame and a subsequent frame. A first phase rotation of the target frame is estimated based on pilot symbols in the target frame and an expected sequence of pilot symbols. A second phase rotation of the subsequent frame is estimated based on pilot symbols in the subsequent frame and the expected sequence of pilot symbols. Responsive to the first phase rotation not matching the second phase rotation, a location of a cycle slip in the target frame is estimated. The cycle slip is then corrected based on the detected location.
The invention has other advantages and features which will be more readily apparent from the following detailed description of the invention and the appended claims, when taken in conjunction with the accompanying drawings, in which:
A receiver architecture and method for data recovery is described for transmissions received over an optical fiber channel. Cycle slip detection and correction methods detect and correct errors in the recovered data symbols caused by cycle slips. The described methods can beneficially be implemented with a low-complexity forward architecture suitable for high speed parallel implementation. In particular, the described methods can be implemented as a simple post-processing stage after a Viterbi & Viterbi (VV) based carrier phase estimator or any other carrier phase estimator. In one embodiment, a 100 Gb/s DP-QPSK optical system can achieve eliminate or nearly eliminate the performance degradation caused by cycle slips at a post-error correcting bit error rate of 10−6 with an 1% overhead.
The transmitter 110 comprises an encoder 112, a modulator 114, a transmitter (Tx) digital signal processor (DSP) 116, and Tx analog front end (AFE) 118. The encoder 112 receives input data 105 and encodes the data for transmission over the optical network. For example, in one embodiment, the encoder 112 encodes the input data 105 using forward error correction (FEC) codes that will enable the receiver 160 to detect, and in many cases, correct errors in the data received over the channel 130. The encoder 112 includes a pilot symbol insertion block (PSIB) 113 to insert pilot symbols in the data stream, which will be used for cycle slip detection and correction in the receiver 160, as will be described in further detail below. The modulator 114 modulates the encoded data via one or more carrier signals for transmission over the channel 130. For example, in one embodiment, the modulator 114 applies phase-shift keying (PSK) to the encoded data. For example, the encoded data may be modulated according to binary phase shift keying (BPSK), quadrature phase shift keying (QPSK), or another modulation scheme. The Tx DSP 116 adapts (by filtering, etc.) the modulator's output signal according to the channel characteristics in order to improve the overall performance of the transmitter 110. The Tx AFE 118 further processes and converts the Tx DSP's digital output signal to the analog domain before it is passed to the optical transmitter (Optical Tx) 120 where it is converted to an optical signal and transmitted via the channel 130. One example of the optical transmitter 120 transmits independent modulations on both polarizations of the optical carrier. An example modulation is QPSK, though other modulations can be used, and the choice can be made to transmit on either one or both polarizations. Although the Tx DSP 116 is illustrated as a separate block from the encoder 112 and modulator 114, in some embodiments all or parts of the encoder 112 and/or modulator 114 may be implemented in the Tx DSP 116 or a separate DSP.
In addition to the illustrated components, the transmitter 110 may comprise other conventional features of a transmitter 110 which are omitted from
The channel 130 may have a limited frequency bandwidth and may act as a filter on the transmitted data. Transmission over the channel 130 may add noise to the transmitted signal including various types of random disturbances arising from outside or within the communication system 100. Furthermore, the channel 130 may introduce fading and/or attenuation effects to the transmitted data. Additionally, the channel 130 may introduce chromatic dispersion (CD) and polarization mode dispersion (PMD) effects that cause a spreading of pulses in the channel 130. Based on these imperfections in the channel 130, the receiver 160 is designed to process the received data and recover the input data 105.
In general, the optical front end 150 receives the optical signal, converts the optical signal to an electrical signal, and passes the electrical signal to the receiver 160. The receiver 160 receives the encoded and modulated data from the transmitter 110 via the optical transmitter 120, communication channel 130, and optical front end 150, and produces recovered data 175 representative of the input data 105. The receiver 160 includes a receiver (Rx) analog front end (AFE) 168, an RX DSP 166, a data recovery system 170, a demodulator 164, and a decoder 162. The Rx AFE 168 samples the analog signal from the optical front end 150 to convert the signal to the digital domain. The Rx DSP 166 further processes the digital signal by applying one or more filters to improve signal quality. A data recovery system 170 compensates for phase and frequency errors in the received data symbols, including detection and correction of cycle slips. Embodiments of a data recovery system 170 are described in further detail below. The demodulator 164 receives the compensated signal from the data recovery system 170 and demodulates the signal. The decoder 162 decodes the demodulated signal (e.g., using error correction codes) to recover the original input data 105. Although the Rx DSP 166 is illustrated as a separate block from the decoder 162, demodulator 164, and data recovery system 170 in some embodiments all or parts of the decoder 162, demodulator 164, and/or data recovery system 170 may be implemented in the Rx DSP 166 or a separate DSP.
In addition to the illustrated components, the receiver 160 may comprise other conventional features of a receiver 160 which are omitted from
Components of the transmitter 110 and the receiver 160 described herein may be implemented, for example, as an integrated circuit (e.g., an Application-Specific Integrated Circuit (ASIC) or using a field-programmable gate array (FPGA), in software (e.g., loading program instructions to a processor (e.g., a digital signal processor (DSP)) from a non-transitory computer-readable storage medium and executing the instructions by the processor), or by a combination of hardware and software.
The carrier phase recovery block 210 comprises a phase locked loop (PLL) 212 (e.g., a low latency PLL) and carrier phase estimator (CPE) 214 (e.g., a Viterbi-Viterbi carrier phase estimator). The PLL 212 and CPE 214 operate to detect and correct for frequency and phase noise in the received data signal. While the PLL 212 and CPE 214 are generally effective at compensating for phase errors less than the constellation rotational-symmetry (for example, 90 degrees in QPSK and 180 degrees in BPSK), the PLL 212 and CPE 214 generally do not correct for full phase errors that map to a different symbol in the constellation based on the rotational symmetry (i.e., cycle slips). Instead, when a cycle slip occurs, the PLL 212 and CPE 214 may map the symbols to the wrong constellation points, thus causing corruption of the data.
The discrete time PLL input signal xk can be modeled as
x
k
=s
k
·e
jθ
+n
k, (1)
where sk is the complex-valued transmitted symbol at the k-th time instant, θk is the cumulative phase effect due to carrier offset, frequency fluctuations, and laser noise phase; nk are independent identically-distributed (iid) complex Gaussian random variables with zero mean and variance N0/2 per dimension. The output {tilde over (x)}k of the PLL 212 is given by
{tilde over (x)}
k
=x
k
·e
−j{circumflex over (θ)}
=s
k
·e
jψk
+ñ
k, (2)
where {circumflex over (θ)}k is the PLL phase correction, ψk=θk−{circumflex over (θ)}k is the residual phase error, and ñk=nk·e−j{circumflex over (θ)}
{circumflex over (x)}
k
=s
k
·e
−j(ψ
−{circumflex over (ψ)}
)
+{circumflex over (n)}
k (3)
where ψk is the error phase estimated by carrier phase estimator, applying, for example, a Viterbi-Viterbi algorithm, and {circumflex over (n)}k=ñk·e−j{circumflex over (ψ)}k.
The cycle slip correction block 230 detects and corrects for cycles slips based on a comparison of a sequence of pilot symbols (e.g., inserted by the pilot symbol insertion block 113 of
and the throughput expansion rate is given by
The sequence of transmitted pilot symbols can be comprises a predefined sequence (i.e., known by the receiver 160) which may be operator configurable. For example, in one embodiment, the predefined sequence of pilot symbols is generated based on a pseudo random bit sequence (PRBS). In various embodiments, the pilot symbol block may be inserted before or after the data symbol block, or pilot symbols may be interleaved with data symbols.
In one embodiment, P=1 and L=G*S where G is a granularity factor (generally constant) and S is programmable between a range of values and may be received, for example, via a configuration control signal. In one embodiment, S is programmable between 5 and 40. The upper end of the programmable range for S may be configured such that increasing S further ceases to show significant performance improvements in a particular architecture. The lower bound of the range of S may be selected such that decreasing S further may cause the maximum throughput of the chip to be exceeded. In one embodiment, a granularity factor G=4 is used as a factor that provides sufficient flexibility with manageable complexity. Other granularity factors may be used depending on the expected operating conditions and desired design tradeoffs.
The pilot symbol synchronization block 220 finds the position of the inserted pilot symbols in the received sequence of frames. In one embodiment, the pilot symbol synchronization block 220 applies a standard, optimum or sub-optimum, maximum likelihood searching techniques based on the knowledge of the transmitted pilot symbol sequence. In the cycle slip correction and detection block 230, the sequence of data frames are provided to a first-in-first-out (FIFO) buffer 232 and to the Pilot Aided-Cycle Slip Detector (PA-CSD) 234. To process the ith data frame, the PA-CSD 234 observes the sequence of M consecutive received pilot symbol blocks centered (or approximately centered) around the ith frame (e.g., frames i−M/2, . . . , i, . . . , i+M/2) and choose the rotation phase {circumflex over (φ)}i that maximizes the probability of these M received pilot symbol blocks being rotated by{circumflex over (φ)}i. For example, in the case where P=1, the value of {circumflex over (φ)}i at the ith sequence of M received pilot symbol blocks is determined as:
where, Prob(A|B) denotes the conditional probability of A given B, ki is the time index of the ith pilot symbol, di is the ith noiseless pilot symbol, {circumflex over (x)}k
because the constellations have a symmetry of 90 degrees. However, different sets of Φ are possible depending on the constellation for the applicable modulation scheme. It will be apparent to those of ordinary skill that {circumflex over (φ)}i can be similarly determined in the case where P≧2.
In one embodiment, a minimum Euclidean distance approach is used to estimate the rotation phase {circumflex over (φ)}i. Here, the rotation phase {circumflex over (φ)}i is chosen as the rotation phase φi that minimizes the Euclidean distance between the expected sequence of pilot symbols transmitted by the transmitter 110 and the corresponding pilot symbols received at the input of the cycle slip correction block 230 rotated by φi. For example, in the case where P=1, the rotation phase {circumflex over (φ)}i is determined as:
In another embodiment, a slicer-based approach is used. In this embodiment, received symbols are first sliced to map the symbols to the closest point on the symbol constellation. Then the rotation phase {circumflex over (φ)}i is determined based on the sliced symbols as described in the equations below:
where Q is the set of possible transmitted constellation points. It will be apparent to those of ordinary skill that {circumflex over (φ)}i can be similarly determined in the case where P≧2.
In some situations, the slicer-based approach of equations (6)-(7) achieves better performance for high values of M (e.g., above a threshold value of M) while the approach of equation (5) has better performance for low values of M (below the threshold value of M). Thus, the different approaches may be selected depending on the operating conditions and desired design tradeoffs.
The phase correction block 236 receives {circumflex over (φ)}i and the symbols of the ith frame (temporarily stored in FIFO 232) and rotates the symbols of the ith frame by {circumflex over (φ)}i. For example, in one embodiment, the phase correction block 236 rotates each of the L symbols located at the center (or approximate center) of the M received pilot symbol blocks by the phase angle {circumflex over (φ)}i. The new sequence of symbols, denoted
k
={circumflex over (x)}
k· for k=ki−1+1, . . . , ki and all i (8)
Finally, once the cycle slips are corrected by the cycle slip correction block 230, the pilot symbol deletion block 240 remove the pilot symbol blocks from {circumflex over (x)}k generating a new sequence of symbols {hacek over (x)}k. In one embodiment, the phase correction block 236 and pilot symbol deletion block 240 operate in the opposite order or the blocks 236, 240 are combined such that the rotation is performed only on the D data symbols of the ith frame, and not necessarily on the P pilot symbols of the ith frame that will be discarded.
Although a value of M=5 is shown in the example of
where
is the possible offset produced by a cycle slip, dk is the kth pilot symbol value, and {circumflex over (φ)}i indicates the estimated phase offset corresponding to the ith block. When {circumflex over (φ)}i≠{circumflex over (φ)}i+1, a cycle slip in ith frame is detected. Equation (9) assumes independent identical distributed Gaussian noise. However, depending on the equivalent channel model at the input of the PA-CSD 634, different criteria can be used. The above equation can be easily modified in the case where the pilot symbols are positioned differently within the frame.
Once a cycle slip is detected by the PA-CSD 634, the PFE 652 determines the exact position inside the frame. In one embodiment, the PFE 652 compares the phase estimation {circumflex over (ψ)}k from the carrier phase estimator 214 and a delayed version of the phase estimate, {circumflex over (ψ)}k−N. A cycle slip caused by additive noise generates a fast jump of the phase estimated by the carrier phase estimator 214 in the direction of cycle slip. Therefore, the phase difference Δ{circumflex over (ψ)}k will have a peak value where the fast jump (i.e., the cycle slip) takes places. On the other hand, a cycle slip caused by a fast change in the phase noise (which generally cannot be tracked by the carrier phase estimator 214) generates a peak value in Δ{circumflex over (ψ)}k in the direction of cycle slip. Based on these principles, the position of the cycle slip is then estimated as
where Δ{circumflex over (ψ)}k={circumflex over (ψ)}k−N−{circumflex over (ψ)}k, Δ{circumflex over (φ)}i={circumflex over (φ)}i−{circumflex over (φ)}i+1. Parameter N may be determined empirically based on operating conditions and desired tradeoffs. After the cycle slip position is estimated, the cycle slip effect is canceled for the affected symbols.
Although the detailed description contains many specifics, these should not be construed as limiting the scope but merely as illustrating different examples and aspects of the described embodiments. It should be appreciated that the scope of the described embodiments includes other embodiments not discussed in detail above. For example, the functionality of the various components and the processes described above can be performed by hardware, firmware, software, and/or combinations thereof.
Various other modifications, changes and variations which will be apparent to those skilled in the art may be made in the arrangement, operation and details of the method and apparatus of the described embodiments disclosed herein without departing from the spirit and scope of the invention as defined in the appended claims. Therefore, the scope of the invention should be determined by the appended claims and their legal equivalents.
This application is a continuation of U.S. patent application Ser. No. 15/351,173, filed on Nov. 14, 2016, which is a continuation of U.S. patent application Ser. No. 14/226,387, now U.S. Pat. No. 9,496,967, entitled “Cycle Slip Compensation in a Coherent Receiver” filed on Mar. 26, 2014, which is a divisional of U.S. patent application Ser. No. 14/209,867 now U.S. Pat. No. 9,306,676, entitled “Cycle Slip Compensation in a Coherent Receiver” filed on Mar. 13, 2014, which claims the benefit of U.S. Provisional Application No. 61/784,305 entitled “Cycle Slip Compensation” filed on Mar. 14, 2013 to Mario Alejandro Castrillon, et al., the contents of which are each incorporated by reference herein.
Number | Date | Country | |
---|---|---|---|
61784305 | Mar 2013 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 14209867 | Mar 2014 | US |
Child | 14226387 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 15351173 | Nov 2016 | US |
Child | 15714990 | US | |
Parent | 14226387 | Mar 2014 | US |
Child | 15351173 | US |