1. Field of the Art
The disclosure relates generally to communication systems, and more specifically, to using equalizer coefficients to improve performance and diagnostics in an optical receiver system.
2. Description of the Related Art
Optical communication systems typically communicate data over a plurality of channels corresponding to different phases and/or polarizations of the optical signal. For example, in a dual-polarization coherent optical communication system, data is communicated over four channels: (1) a horizontally polarized in-phase (HI) channel; (2) a horizontally polarized quadrature (HQ) channel; (3) a vertically polarized in-phase (VI) channel; and (4) a vertically polarized quadrature (VQ) channel. While the data communicated over the different channels is typically aligned relative to a common clock when transmitted by the transmitter, delay (or skew) may be introduced into one or more of the channels based on characteristics of the transmitter, receiver, and/or the optical fiber. As a result, the relative timing of the data in the various channels may be misaligned at the receiver, causing degradation of the recovered data.
A system performs automatic skew calibration to compensate for skew between different input channels. An analog-to-digital converter converts an analog signal vector to a digital signal vector. The analog signal vector has a plurality of channels. The analog-to-digital converter samples each of the plurality of channels based on a corresponding sampling clock. A skew compensation block generates the corresponding sampling clocks for each of the plurality of channels based on skew compensation values specifying an amount of skew compensation for each of the plurality of channels. An adaptive equalizer equalizes the digital signal vector. The adaptive equalizer furthermore adaptively updates coefficients based on characteristics of the digital signal vector and generates an equalized signal vector. A skew estimator obtains the coefficients of the adaptive equalizer and determines an amount of skew between at least two of the plurality of channels. The skew estimator then programs the skew compensation block with the skew compensation values based on the determined amount of skew.
In one embodiment of a calibration process, a receiver sets a first skew compensation value for a first channel to a first value for testing and sets a second skew compensation value for a second channel to a second skew value for testing. The first and second skew compensation values control an amount of skew compensation in their respective channels. The receiver obtains coefficients of an adaptive equalizer when operating the receiver using the first and second skew compensation values. A skew metric is computed that approximates a skew between the first and second channel based on the obtained coefficients of the adaptive equalizer. At least one of the first skew compensation value for the first channel and the second skew compensation value for the second channel is adjusted based on the skew metric.
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:
Overview
The disclosed embodiments generally relate to improving performance of coherent communication systems that employ electronic equalization of communication channels in the receiver. In a first aspect, a receiver performs an automatic skew calibration to compensate for skew between different input channels during a calibration stage of the receiver. The receiver skew is estimated by observing the coefficients of an adaptive equalizer, which adjusts the coefficients based on time-varying properties of the multi-channel input signal. The receiver compensates for the estimated skew by adjusting the phases of the sampling clocks for the different channels.
In a second aspect, the receiver performs channel diagnostics to extract channel information from the equalizer coefficients. The channel information can be used to automatically estimate differential group delay introduced by channel impairments during real-time operation. The differential group delay is accurately estimated based on the equalizer coefficients using a frequency averaging approach. Framer information can furthermore be utilized to estimate differential group delay that is an integer multiple of the symbol rate. The diagnostics information may also be used to reset the device when substantial performance degradation is detected.
High Level System Architecture
The transmitter 110 comprises a transmitter (Tx) framer 111, an encoder 112, a modulator 114, a Tx digital signal processor (DSP) 116, and a Tx analog front end (AFE) 118. The Tx framer 111 receives input data 105 from a host interface and formats the data for processing by the transmitter 110. The encoder 112 encodes the framed 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 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) or different phase-shift keying (DPSK) to the encoded data. 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. Example modulations are binary phase shift keying (BPSK) and quadrature phase-shift keying (QPSK), though other modulations can be used, and the choice can be made to transmit on either one or both polarizations. Additional example modulations include differential BPSK (DBPSK) and differential QPSK (DQPSK).
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 distortion of the signals transmitted through the channel 130.
The optical front end 150 receives the optical signal from the channel 130, 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 demodulator 164, a decoder 162, and an Rx framer 163. The Rx AFE 168 samples the analog signal from the optical front end 150 based on a clock signal 181 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. The demodulator 164 receives the modulated signal from the Rx DSP 166 and demodulates the signal. The decoder 162 decodes the demodulated signal (e.g., using error correction codes). The Rx framer 163 aligns and formats the decoded data for synchronization with a host interface that receives the recovered input data 105.
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 microprocessor or a digital signal processor (DSP)) from a computer-readable storage medium and executing the instructions by the processor), or by a combination of hardware and software.
The physical characteristics of the transmitter 110 or receiver 160 (or both) may introduce skew between the various channels HI, HQ, VI, VQ in the digital signal vector 205 if left uncompensated. For example, skew between the HI and HQ channels or between the VI and VQ channels may introduce error in the quadrature phase relationships between the signals (i.e., the phase difference between the signals may vary from the expected 90° relationship). Furthermore, skew between the horizontal and vertical polarizations may be introduced that appear as a differential group delay (DGD) between the polarizations. The receiver 160 includes various components to compensate for this skew, thereby improving performance of the receiver 160.
In the AFE 168, a multi-channel analog-to-digital converter (ADC) 210 samples the analog input signal vector 203 to generate a digital signal vector 205. Each channel of the analog input signal 203 is sampled by the ADC 210 according to sampling clocks 207 which each corresponding to a different channel of the analog input signal 203 (e.g., sampling clocks CLKHI, CLKHQ, CLKVI, CLKVQ). By varying the relative phases of the sampling clocks 207, the different channels of the analog input signal 203 can be sampled with different relative timings, thereby compensating for skew.
A skew compensation block 220 generates the sampling clocks 207 for each channel based on an input clock signal 211. Based on a programmed level of skew compensation for each channel, the skew compensation block 211 adjusts the phase of the input clock signal 211 separately for each channel to generate each of the sampling clocks 207. For example, in one embodiment, the skew compensation block 220 comprises four programmable phase interpolators, one for each of the four channels. The phase interpolator for a given channel includes a programmable skew compensation register that can be programmed to introduce a particular amount of phase shift to the sampling clock for that channel. Based on the programmed value, the phase interpolator applies a phase shift to the input clock signal 211 to generate the sampling clock for the corresponding channel, thus compensating for the detected skew. Values for the skew compensation registers in the skew compensation block 220 may be programmed via a control signal 213 during a calibration phase described below. In one embodiment, this calibration phase is performed during an initial set up of the receiver at the beginning of its lifetime. Additionally, the receiver 160 may be recalibrated at later times (e.g., during a service call) by re-performing the calibration. Typically, the calibration is performed under ideal or near ideal laboratory conditions to avoid impairments effects introduced by a noisy channel. In alternative embodiments, however, calibration may also be performed periodically or continuously during normal operation of the receiver 160 in the field.
The DSP 166 receives the digital signal vector 205 and performs various processing operations to generate the digital output signal 209 representing the originally transmitted data. In one embodiment, the DSP 166 includes an equalizer 230 to equalize the channels of the digital signal vector 205.
In the DSP 166, the equalizer 230 compensates for channel impairments present in the digital signal vector 205 including, for example, chromatic dispersion, polarization mode dispersion, nonlinear optic effects, and polarization scrambling. In one embodiment, the equalizer 230 comprises a four channel feedforward adaptive equalizer that implements the following equation:
where rHI, rHQ, rVI, and rVQ are the inputs to the equalizer 230 representing the horizontally polarized in-phase channel (HI), the horizontally polarized quadrature channel (HQ), the vertically polarized in-phase channel (VI), and the vertically polarized quadrature channel respectively (VQ). qHI, qHQ, qVI, and qVQ are the outputs out of the equalizer 230. cab, (for a,bε[0, 1, 2, 3]) are filters (e.g., finite impulse response (FIR) filters), each having L taps (e.g., 16 taps). Here, the term “filter” is used in the mathematical sense to describe each of the individual vectors cab in equation (1). These vectors cab do not necessarily operate to filter signals independently of each other, but may instead be applied collectively in the equalizer 230 to achieve an overall filtering effect between an input signal and an output signal as described in equation (1).
During normal operation, the adaptive equalizer 230 dynamically adapts the taps of each of the filters cab based on the time-varying properties of the communication channel. The matrix C converges to values that compensate for the impairments and continuously changes as the impairment conditions change. For example, in one embodiment a least-mean-square (LMS) algorithm is used to adapt the taps of the filters cab to produce the least mean squares of an error signal. The equalizer 230 may alternatively use a different adaptive algorithm. As a result of the adaptive nature of the equalizer 230, the filters cab will adapt if the relative timing of an input channel changes with respect to the other channels. Thus, the equalizer coefficients can provide useful information about the operation of the transmitter 110, receiver 160 and channel 130, which can be used for skew calibration and channel diagnostics according to various techniques described below.
Other components of the DSP 166 are omitted for simplicity in the illustration. In various embodiments, the DSP 166 may furthermore include processing and filtering components such as, for example, a bulk chromatic dispersion compensation block, a polarization mode dispersion compensation block, a timing recovery block, etc.
In one embodiment, a separate processor 180 is coupled to the DSP 166 and AFE 168 includes a skew estimator 240 and a channel diagnostics block 270. For example, the processor 180 may be embodied on the same chip as the receiver 160 or in a physically separate device (e.g., on a host board). In one embodiment, the functions attributed to processor 180 described herein are carried out by executing program instructions stored to a computer-readable storage medium. Alternatively, the processor 180 can be implemented as an application specific integrated circuit or a field programmable gate array (FPGA)).
The skew estimator 240 is utilized during calibration of the receiver 160. For example, in one embodiment, the skew estimator 240 receives the equalizer coefficients from the equalizer 230 and generates the control signal 213 to program the skew compensation registers of the skew compensation block 220 during calibration. Further details regarding operation of the skew estimator are described below.
In another embodiment, the control signal 213 can be provided by a different external device. For example, in one embodiment, the equalizer coefficients are outputted external to the DSP 166 where they can be to manually determine calibration settings. The skew compensation values can then be programmed into the skew compensation block 220 via firmware, a field programmable gate array (FPGA), or erasable programmable read-only memory (EPROM) that controls the skew compensation block 220.
The channel diagnostics block 270 also receives the equalizer coefficients from the equalizer 230 and generates various types of channel information such as a differential group delay (DGD) estimate. In one embodiment, information from the Rx framer 163 is also used in performing the channel diagnostics computations. The channel diagnostics information is used to provide a channel diagnostics output 260 useful for calibration and reporting. Additionally, the DGD estimate could optionally be used by the skew estimator 240 to estimate H/V skew during the calibration phase.
Skew Compensation Based on Equalizer Coefficients
As described above, equalizer coefficients are used during calibration of the receiver 160 to estimate the skew between the channels and program the skew compensation block 220 to compensate for the estimated skew. It is generally desirable to perform calibration under ideal or near ideal operating conditions. Thus, in one embodiment, calibration is performed with the transmitter 110 and receiver 160 connected back-to-back (effectively bypassing the optical fiber). For example, in one embodiment, the back-to-back setup involves connecting the transmitter 110 to the receiver 160 using a channel 130 of negligible length. Furthermore, in one embodiment, certain features of the receiver 160 such as coarse carrier recovery may be bypassed during calibration. Additionally, it may be desirable to reduce the dynamic response of the equalizer 230 using available programmable options.
The calibration process is based on the adaptive operation of the equalizer 230 with respect to the incoming signal. Specifically, if the timing of rHQ changes with respect to the other inputs, the corresponding filters c01(k), c11(k), c21(k), c31(k) will adapt to reflect the change in timing. The inputs rHI, rHQ, rVI, and rVQ are not independent, but rather will maintain certain relationships when the channels are properly aligned. Because the equalizer 230 adapts to the channel data, the equalizer filters will maintain the channel relationships when the signals are aligned. Particularly, when the channels are properly aligned, the filters will reflect the appropriate in-phase and quadrature relationships between the input signals. Additionally, when the channels are properly aligned, the filters will maintain time reversal symmetry between the horizontally polarized inputs and the vertically polarized inputs. In an equalizer 230 having the characteristics of equation (1), the filter coefficients will have the following relationships when the HI and HQ channels are aligned: c11≈c00 and −c01≈c10. Furthermore, when the VI and VQ channels are aligned, the following relationships will be maintained: c33≈c22 and −c23≈c32.
When the channels are misaligned (e.g., due to skew between HI and HQ or between VI and VQ), the filter coefficients will adapt in response and will reflect the misalignment. In order to determine the amount of skew in the HQ channel relative to the HI channel, c11 is compared to c00 and −c01 is compared to c10. Similarly, to determine the amount of skew in the VQ channel relative to the VI channel, c33 is compared to c22 and −c23 is compared to c32.
Correlation Method
The skew estimator 240 obtains 404 the coefficients of the relevant equalizer filters from the equalizer 230 while operating the receiver 160 with the skew compensation block 220 set based on the current skew parameters. For example, in one embodiment, a test data pattern is transmitted to the receiver 160 until the equalizer 230 converges (i.e., the filter coefficients adapt until they reach substantially stable values). The relevant filter coefficients are then read from the equalizer 230. To determine the skew between HI and HQ channels, coefficients of filters c11, c00, c01, and c10 are obtained. To determine the amount of skew between VQ and VI channels, coefficients of filters c33, c22, c23, and c32 are obtained. The skew estimator 240 computes 406 a skew metric between the in-phase and quadrature channels based on the obtained filter coefficients. For example, in one embodiment the skew metric is determined based on a correlation between the relevant filters using the following equation:
R(ab,cd)(τI,τQ)=Σn=1Lcab(n,τI)·ccd(n,τQ);a,bε[0,1,2,3] (2)
where R(ab,cd)(τI, τQ) represents a correlation metric between filters cab and ccd. For example, to determine skew between HI and HQ channels, R(00,11)(τ1, τQ) and R(01,10) (τI, τQ) are computed, which represent the correlation between filters c00 and c11 and between filters c01 and c10 respectively. Then an overall correlation metric RHIHQ for the horizontally polarized channels can be determined by scaling each term by the energy of the corresponding filters:
Similarly, to determine skew between VI and VQ channels, R(33,22) (τI, τQ) and R(23,32) (τI, τQ) are computed, which represent the correlation between filters c33 and c22 and between filters c23 and c32 respectively. Then an overall correlation metric RVIVQ for the vertically polarized channels can be determined as follows:
In step 408, the skew estimator 240 determines whether there are additional skew compensation value to test. If there are additional skew compensation values to test, the skew compensation block 220 is programmed with new values and the process repeats from step 402. Once the full range of skew compensation values are tested, the optimal skew values are determined 410 based on the skew metrics (e.g., the metric indicating the highest correlation). The skew compensation block 220 is then configured with these optimal values for use during operation of the receiver 160. For example, in one embodiment, skew compensation registers associated with phase interpolators are configured to apply an appropriate amount of phase shift to the sampling clocks associated with each channel based on the determined skew.
Various strategies may be employed in selecting the range of skew compensation values for testing. For example, in one embodiment, the skew compensation registers for the in phase channel is set to the middle of its range (e.g., τI=0) and the skew compensation register for the quadrature channel is swept across a predefined testing range (e.g., ±2 unit intervals). Alternatively, the skew compensation register for the quadrature channel is fixed and the skew parameter for the in phase channel is swept across its range. In other embodiments, parametric sweeps may be performed to test different combinations of values for τI, τQ.
Center of Gravity Method
In another embodiment, a “center of gravity” method is used to determine the relative skew between the channels based on the equalizer coefficients.
where CGH is the center of gravity of the horizontally polarized channels and CGV is the center of gravity for the vertically polarized channels. In another embodiment, additional coefficients may be considered in the center of gravity equations. For example, in another embodiment, the centers of gravity may be determined as follows:
In yet another alternative embodiment, the centers of gravity may be determined as follows:
In yet another alternative embodiment, the centers of gravity may be determined as follows:
A sweep of the skew compensation values then begins. For example, in one embodiment, the skew compensation values for the in-phase channels (HI and VI) are swept across a predefined testing range while the quadrature channels (HQ and VQ) remain constant at their centered value. Alternatively, the skew compensation values for the quadrature channels may be swept while the in-phase channels are held constant. To begin the sweep, the skew compensation value for the channel under calibration is set 458 to a first test value in a predefined testing range. The equalizer coefficients are then determined 460 once the equalizer 230 converges with the current skew compensation values. Error values for the current skew parameter are then calculated 512 as follows:
where x represents a tap range parameter (e.g., x=3) that defines a range around the tap closest to the center of gravity to be used in the error calculation, round ( ) is a function that returns the nearest integer of the argument, and where:
HIn=c00(n)−c10(n)
HQn=c11(n)+c01(n)
VIn=c22(n)−c32(n)
VQn=c33(n)+c23(n) (11)
If in step 464, there are additional skew compensation values to test in the defined testing range, the next skew compensation value is applied and the process loops back to step 468.
Otherwise, once correlation metrics are determined for all skew compensation values in the desired testing range, the minimum error values ΔHmin and ΔVmin are determined 466 from all of the tested skew compensation values. The skew compensation values associated with the minimum error values are written to the skew compensation registers of the skew compensation block 220.
Group Delay Method
In another embodiment, the skew metric is determined by computing a difference in group delay between the relevant filters being compared. Let the discrete Fourier transform (DFT) of a filter be H(w)=A(ω)e−iΘ(ω). Where A(ω) is the magnitude of the frequency response and Θ(ω) is the phase. Then the group delay is defined as:
For H polarization IQ skew, the group delay is computed as follows:
The group delays can be combined by to compute an overall skew value by scaling the group delays by the magnitude response at DC. In other words,
τHIHQ=([A002(0)+A112(0)]τ(00,11)+[A102(0)+A012(0)]τ(10,01))/[A002(0)+A112(0)+A102(0)+A012(0)] (13b)
For V polarization IQ skew, the group delay is computed as follows:
The group delays can be combined to compute an overall skew value by scaling the group delays by the magnitude response at DC. In other words,
τVIVQ=([A222(0)+A332(0)]τ(22,33)+[A322(0)+A232(0)]τ(32,23))/[A222(0)+A332(0)+A322(0)+A232(0)] (14b)
The skew compensation block 220 can then be set to compensate for the overall skew determined in Eqs. (13b) and (14b) above.
HV Skew Calibration
In another embodiment, a calibration process detects skew between the horizontally polarized channels and the vertically polarized channels (referred to herein as “HV skew”) and programs the skew compensation block 220 to compensate for the detected skew. To extract the HV skew information from the equalizer coefficients, IQ and time reversal symmetry relationships are applied. This relationship is embodied in the calculation of differential group delay (DGD) between the horizontally and vertically polarized channels.
The above described process may be used to generate skew metrics while varying the HI and HQ channels together across a testing range while maintaining the determined offset between HI and HQ. Alternatively, the above described process may be used to generate skew metrics while varying the VI and VQ channels together across a testing range while maintaining the determined offset between VI and VQ.
Channel Diagnostics Using Equalizer Coefficients
In one embodiment, the receiver 160 can also use the equalizer coefficients to perform channel diagnostics. In contrast to the calibration procedure discussed above, channel diagnostics is typically performed dynamically during real-time operation of the receiver 160 under field conditions in order to obtain real-time diagnostic information. Channel diagnostics estimate various time-varying characteristics of the incoming data signal, including, for example, the differential group delay (DGD) in the optical channel 130. In one embodiment, the obtained information can then be used to compensate for impairments introduced in the data channel. Additionally, the channel diagnostics information may be externally reported where it may be logged or analyzed to extract various types of performance information. The channel diagnostics process will produce the most accurate results when the receiver 160 is correctly calibrated to compensate for I/Q skew, which can be accomplished using any of the techniques mentioned above.
As with the skew compensation process described above, channel diagnostics information can be obtained based on the adaptive operation of the equalizer 230 with respect to the incoming signal. Eq. (1) can be rewritten as filtering electrical fields of the input signals to the equalizer 230 as follows:
where:
hHHI=0.5*(c00+c11);
hHHQ=0.5*(c10−c02);
hHVI=0.5*(c02+c13);
hHVQ=0.5*(c12−c03);
hVHI=0.5*(c20+c31);
hVHQ=0.5*(c30−c21);
hVVI=0.5*(c22+c33); and
hVVQ=0.5*(c32−c23). (16)
If Eqn. (15) is represented in the frequency domain, the operation will be matrix multiplication expressed as:
Q(ω)=H(ω)R(ω). (17)
where ω is the angular frequency. If H(ω) is assumed to be a perfect inverse of the linear optical channel distortions, it is then possible to derive the channel parameters from H(ω):
where D(ω) represents the frequency response of chromatic dispersion, and Un(ω) and En are matrices representing local birefringence and polarization dependent Loss (PD), respectively.
The chromatic dispersion (CD) portion of the channel distortions can be isolated from H(ω) by taking the determinant of H(ω):
Therefore, the amount of CD can be determined by quadratic fitting of the unwrapped phase of det(H(ω))1/2.
To isolate the effect of differential group delay (DGD), H(ω) is normalized by its determinant. The normalization will remove the effect of CD and minimize the effect of polarization-dependent loss (PDL):
T(ω)=H(ω)/det(H(ω)). (20)
Where T(ω) represents the normalized frequency response of an ideal equalizer. The DGD can be derived from T(ω) using a frequency based method. In general, it is reasonable to assume that T(ω) for |ω−ω0|<B/(4π) can be approximated by:
where B is the symbol rate, ω0 is the offset between the signal and a local oscillator, and V and U are unitary transformation matrices. From Eqn. (19), DGD can be determined using the following equation:
τ(ωo)=√{square root over (det(dT(ω)/dω))}|ω
The DGD varies with frequency. An accurate estimate of the DGD can be determined, in one embodiment, by averaging DGD over frequency. In one embodiment, DGD can be estimated by based on the fractional group delay of the frequency responses of each equalizer filter.
The equalizer can be seen as an approximation to an ideal fractional delay filter hu(n)=sin(π(n−μ))/(π(n−μ)) over the symbol rate. The DGD varies with frequency due to truncation of these ideal fractional delay filters. The approximation is the result of windowing the filter response. Its effect can be approximated by the convolution of the frequency response of the window function:
Hw(ω)=sin(ω(N/2))/(N sin(ω/2)) (23)
where N is number of taps.
The main lobe of the window function is 2/N. Given that the group delay and magnitude of frequency response of the ideal fractional delay filter are both constant in frequency, the variation is expected to have periods of 2/N. If group delay is averaged over the symbol rate, N/4 cycles of variations are expected. Assuming that the number of taps is a multiple of 4, the overall group delay can be accurately estimated by averaging the group delay over a range of frequencies:
where B is the symbol rate, Ts, is the sampling rate, and
In some implementations, the equalizer 230 itself may naturally add DGD in the receiver 160 based on the operating characteristics of the equalizer 230. This DGD is added after the data has been decoded by the receiver 160, so it does not directly affect performance and need not be compensated for by the skew compensation block 220. Nevertheless, it is desirable for the channel diagnostics information to accurately report the DGD added by the channel 130 and transmitter 110. Because this particular DGD is instead added by the receiver 160, it is desirable to mathematically remove this component from the DGD estimation included in the channel diagnostics output 260.
The DGD introduced by the equalizer 230 may include DGD that is an integer number of symbol periods. Blind adaptation of the equalizer 230 generally cannot resolve DGD that is an integer number of symbol period. Instead, information embedded in the Rx framer 163 can be used to detect this DGD. The DGD estimate can then be adjusted 608 based on the framer information. For example, in one embodiment, the Rx framer 163 monitors the H and V lane skew at the Rx framer 163 and provides this information to the channel diagnostics block 270. In one embodiment, this information can be obtained by reading a specific framer register after the DSP has converged.
The frequency response of the equalizer 230 can be represented as follows:
where G is the frequency response of the equalizer, T is the equalizer induced skew that is an integer multiple of a symbol period determined by the Rx framer 163, and H is an inverse of the linear optical channel distortions. If the raw equalizer frequency response is normalized by its determinant, as described above, it is a reasonable to estimate the skews as being evenly distributed between H and V channels.
The channel diagnostics block 270 then applies an inverse matrix Ito G to adjust the DGD estimate, where:
The estimated DGD, as adjusted, is then outputted 612 in the channel diagnostics output 260.
Steps 602, 604, 606, 608 of
In another embodiment, the channel diagnostics block 270 can use the equalizer coefficients to detect performance degradation resulting in high bit error rate and synchronization loss, and reset the DSP 166 in response. The DSP 166 applies various iterative algorithms that should converge on a solution during start up conditions. However, under certain channel conditions, these algorithms may not properly converge to best solution, thus resulting in high bit error rate and synchronization loss. By resetting the DSP 166 when such conditions are detected, the algorithms can start over under different channel conditions, and may thus converge to a new solution that provides improved performance.
Polarization Scrambling
In an alternative embodiment, an accurate estimate of DGD can instead be determined by computing an average DGD over a range of polarization states instead of averaging over a range of frequencies.
In another embodiment, polarization scrambling takes place in the DSP 166 instead of in the optical front end 850. In this embodiment, a matrix rotator in the DSP 166 applies a matrix rotation R to the input vector:
where and R (θ, φ) is the rotation matrix given by
and where θ the polarization angle and φ is a random phase angle. More general forms of the rotation matrix may be used. Here, θ and φ are varied over time to achieve the different polarization angles. In one embodiment, θ and φ are varied such that all combinations of θ and φ are covered within a predefined time period. Multiple group delays τt (or τθ, φ) are calculated for different polarization orientations at a single frequency ω0, and averaged over polarization orientations to calculate a mean group delay
In another embodiment, the matrix rotation R is applied to the equalizer coefficients outside the DSP 166, such as in processor 180, for the purpose of estimating the DGD.
Optionally, the process described above can also be performed during calibration of the receiver 160. Here, the optional step 910 is added in which the channel diagnostics block 270 adjusts 910 the skew compensation values of the skew compensation block 220 to compensate for the detected DGD (e.g., by adjusting the phases of the sampling clocks of the horizontally polarized channels relative to the vertically polarized channels or vice versa).
In one embodiment, where polarization scrambling is achieved in by matrix rotation, the receiver 160 includes a separate DSP path for channel estimation that is independent from the data path. An example architecture is illustrated in
Alternatively, in the embodiment of
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 claims the benefit of the following U.S. Provisional Applications: U.S. Provisional Patent Application Ser. No. 61/587,572 entitled “Automatic Skew Compensation Scheme Based on FFE Coefficients,” filed on Jan. 17, 2012 to Shih Cheng Wang; U.S. Provisional Application No. 61/677,460 entitled “Feed-Forward Equalizer Center-of-Gravity Monitoring in High-Speed Optical Communications Networks” filed on Jul. 30, 2012 to Mario R. Hueda, et al.; U.S. Provisional Patent Application Ser. No. 61/683,169 entitled “Automatic Receiver Deskew Procedure,” filed on Aug. 14, 2012 to Diego Crivelli, et al; U.S. Provisional Patent Application Ser. No. 61/704,490 entitled “Method Using Frequency Domain Averaging to Calculate Impairments in a Fiber-Optic Channel to Invert an FFE Matrix,” filed on Sep. 22, 2012 to Shih Cheng Wang, et al.; and U.S. Provisional Application No. 61/745,389 entitled “Method of Estimating Channel Performance Using Equalizer Tap Coefficients and Framer Information” filed on Dec. 21, 2012 to Shih Cheng Wang, et al. The contents of each of the above referenced applications are incorporated by reference herein in their entirety. This application is related to the following applications: U.S. Utility application Ser. No. 13/744,211 entitled “Channel Diagnostics Based on Equalizer Coefficients” filed on Jan. 17, 2013 to Shih Cheng Wang, et al.; and U.S. Utility application Ser. No. 13/744,235 entitled “Reset in a Receiver Using Center of Gravity of Equalizer Coefficients” filed on Jan. 17, 2013 to Mario Rafael Hueda, et al. The contents of each of the above referenced applications are incorporated by reference herein in their entirety.
Number | Name | Date | Kind |
---|---|---|---|
4815103 | Cupo et al. | Mar 1989 | A |
5353312 | Cupo et al. | Oct 1994 | A |
5870372 | Kuribayashi | Feb 1999 | A |
6791995 | Azenkot et al. | Sep 2004 | B1 |
7693214 | Shida | Apr 2010 | B2 |
8077762 | Agazzi | Dec 2011 | B2 |
8565621 | Ibragimov et al. | Oct 2013 | B2 |
8731413 | Dave et al. | May 2014 | B1 |
8873358 | Saito et al. | Oct 2014 | B2 |
20070206963 | Koc | Sep 2007 | A1 |
20080175590 | Perkins et al. | Jul 2008 | A1 |
20090110046 | Aziz et al. | Apr 2009 | A1 |
20100008413 | Yamazaki et al. | Jan 2010 | A1 |
20110064421 | Zhang et al. | Mar 2011 | A1 |
20110150506 | Tanimura et al. | Jun 2011 | A1 |
20110229127 | Sakamoto et al. | Sep 2011 | A1 |
20120189320 | Zelensky et al. | Jul 2012 | A1 |
20120230676 | Mo et al. | Sep 2012 | A1 |
20120257652 | Malipatil et al. | Oct 2012 | A1 |
20120307025 | Werner Lerche et al. | Dec 2012 | A1 |
20130051442 | Cooper et al. | Feb 2013 | A1 |
20130084080 | Shibutani | Apr 2013 | A1 |
20130209089 | Harley et al. | Aug 2013 | A1 |
20130243127 | Chmelar et al. | Sep 2013 | A1 |
20130251082 | Abe, Junichi | Sep 2013 | A1 |
20130259490 | Malouin et al. | Oct 2013 | A1 |
20130343490 | Wertz et al. | Dec 2013 | A1 |
20140161470 | Zelensky et al. | Jun 2014 | A1 |
20140254644 | Gotman et al. | Sep 2014 | A1 |
Entry |
---|
El-Chammas, M. et al., “A 12-Gs/s 81-mW 5-bit Time-Interleaved Flash ADC With Background Timing Skew Calibration,” IEEE Journal Solid-State Circuits, Apr. 2011, pp. 838-847, vol. 46, No. 4. |
Geyer, J.C. et al., “Optical Performance Monitoring using a 43Gb/s Realtime Coherent Receiver (Invited),” IEEE, 2009, pp. 93-94. |
Geyer, J.C. et al., “Performance Monitoring Using Coherent Receivers,” IEEE, 2009, three pages. |
Hauske, F. et al., “Optical Performance Monitoring in Digital Coherent Receivers,” Journal of Lightwave Technology, Aug. 15, 2009, pp. 3623-3631, vol. 27, No. 16. |
Leven, A. et al., “Real-Time Implementation of Digital Signal Processing for Coherent Optical Digital Communication Systems,” IEEE Journal of Selected Topics in Quantum Electronics, Sep./Oct. 2010, pp. 1227-1234, vol. 16, No. 5. |
Nelson, L. et al., “Performance of 46-Gbps Dual-Polarization QPSK Transceiver With Real-Time Coherent Equalization Over High PMD Fiber,” J. Lightwave Tech., Feb. 2009, pp. 158-167, vol. 27, No. 3. |
Viterbi, A. et al., “Nonlinear Estimation of PSK-Modulated Carrier Phase with Application to Burst Digital Transmission,” IEEE Transactions on Information Theory, Jul. 1983, pp. 543-551, vol. 29, No. 4. |
Greshishchev, Y. et al., “A 40GS/s 6b ADC in 65nm CMOS,” ISSCC Dig. IEEE International Solid-State Circuits Conference, 2010, pp. 390-392. |
Number | Date | Country | |
---|---|---|---|
61587572 | Jan 2012 | US | |
61677460 | Jul 2012 | US | |
61683169 | Aug 2012 | US | |
61704490 | Sep 2012 | US | |
61745389 | Dec 2012 | US |