# Patent

## Title
Recovering data from quadrature phase shift keying modulated optical signals

## Abstract

Systems, devices and techniques for processing received QPSK modulated optical signals include sampling the received signal at twice the baud rate, thereby producing samples that are then processed as 9-QAM symbols using a decision directed least squares optimization method. Data bits are then recovered from the resulting symbol estimates. The received optical signal may also include dual polarized signals for increased bandwidth capacity.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This patent document claims the benefit of priority of U.S. Provisional Patent Application No. 61/740,840, filed on Dec. 21, 2012. The entire content of the before-mentioned patent application is incorporated by reference as part of the disclosure of this application.


TECHNICAL FIELD
This patent document relates to optical communications and techniques for recovering data from optical signals.
BACKGROUND
With the advent of transmission of multimedia content such as video over the Internet and other communications networks, there is any ever-increasing need for increased data rate capacity on communication networks. Often, optical communication signals are used in backbones of communications network, where the increased traffic at the edges of the network aggregates to several gigabit of network traffic. Therefore, there is an increasing need for optical communication techniques to meet the demand on communications network data capacity. However, laying down additional optical transmission media such as fiber optics often requires large amount of capital expenditure and may not always be a suitable options due to the expenses involved and other regulatory issues.
Techniques for improved optical transceiver operation are needed.
SUMMARY
This patent document provides, among others, techniques for recovering data from polarization modulated quadrature phase shift keying (PM-QPSK) signals, e.g., PM-QPSK signals that have undergone strong filtering during transfer from a transmitter to a receiver.
In one aspect, methods, apparatus and computer program products for recovering data in an optical communication system are disclosed. One example method includes receiving an optical signal that is modulated using polarized modulation quadrature phase shift keying (PM-QPSK) modulation scheme, sampling the received optical signal to generate a group of samples comprising samples at time instances that are odd integer multiples of T/2, where T represents a symbol time interval of the QPSK modulation scheme, estimating received constellation symbols from the group of samples using a decision-directed least radius distance (DD-LRD) algorithm for blind equalization, wherein the DD-LRD algorithm includes iteratively refining symbol estimates using an adaptive filter, and recovering data bits transmitted in the optical signal from the estimated constellation symbols.
This and other aspects and their implementations are described in greater detail in the drawings, the description and the claims.



BRIEF DESCRIPTION OF DRAWINGS
FIG. 1A is a block diagram of an optical communications system.
FIG. 1B is a block diagram example of an optical transmitter/receiver chain.
FIG. 2 shows a signal waveform that is sampled at multiples of T and T/2.
FIG. 3 depicts example constellations of a sampled optical QPSK signal with the filter of (a) 28 GHz and (b) 24 GHz bandwidth.
FIG. 4 depicts example constellations made up from (a) T/2 samples, (b) DD-LRD processing, (c) T samples, (d) constant modulus adaptation CMA processing, (e) CMA+post-filter processing.
FIG. 5 is a graph depicting curves for measured MSE of the signals after T and T/2 sampling, DD-LRD, CMA, and CMA+post-filter processing.
FIG. 6 is a block diagram representation of a wireless communications test bed.
FIG. 7 is a graph depicting curves for measured BER as a function of OSNR (0.1 nm).
FIG. 8 depicts example constellations of the recovered signal using the algorithm of (a) 9-QAM digital processing, (b) standard CMA, (c) CMA+post-filter.
FIG. 9 is a flow chart representation of a process of optical communications.
FIG. 10 is a block diagram representation of an apparatus for optical communications.
FIG. 11 is a flowchart representation of a digital signal processing technique.



DETAILED DESCRIPTION
The techniques disclosed in the present document can be used to recover data from a Quadrature Phase Shift Keying (QPSK) modulated signal, even under severe filtering of the modulated signal. In the description below, the subject technology is disclosed using optical communication as an example. However, the disclosed techniques may be used in data recovery from any QPSK modulated signals. In the context of optical communications, the disclosed techniques can be used in embodiments of optical receivers such as optical trunk equipment used in Gigabit Ethernet, Synchronous Optical Networks (SONET), and other networks.
In some disclosed embodiments, a 9-QAM (Quadrature Amplitude Modulation) data recovery technique is used to process a portion of a received Polarization Modulated QPSK (PM-QPSK) signal in the presence of strong filtering. Strong filtering may occur, e.g., due to bandwidth limitation of the channel or due to bandwidth limitations in an intervening optical equipment such as a switch, or an optical front end used to receive the optically modulated data. As further discussed below, a decision-directed least radius distance (DD-LRD) algorithm for blind equalization can be used for 9-QAM recovery and inter-symbol interference (ISI) compensation. Simulation and tests have shown that the disclosed techniques are robust under strong filtering to recover a 9-QAM signal, instead of using traditional QPSK techniques only. To highlight certain features, results obtained after running an experiment of 112 Gb/s spectrum shaped PM-QPSK signal, corrupted by Wavelength Selective Switching (WSS) in a 25 GHz channel spacing Nyquist Wavelength Division Multiplexing (WDM). The final equalized signal is detected by Maximum Likelihood Sequence Detection (MLSD) for data bit error rate (BER) measurement. As will be discussed below, a 0.5 dB optical Signal to Noise Ratio (OSNR) tolerance improvement is obtained at a BER of 10−3 comparing to a traditional CMA plus post-filter algorithm.
In optical communication, signals occupying different polarization planes can be used to increase spectral efficiency. For example, a dual polarized signal could double the spectral efficiency (number of data bits transferred per second per Hertz). In some implementations, a Nyquist pulse is generated to achieve the Nyquist limit of spectral efficiency for a given baud rate. However the operating rate is limited by the speed of digital/analog converter (DAC). A partial response system gives the same high spectral efficiency at a cost of OSNR penalty due to its multilevel detection. Recently, a quadrature duobinary (QDB) has been proposed that approaches a spectral efficiency of 4 bit/s/Hz. A simpler and more realistic means is spectrum shaping based on the wavelength selective switch. A post-filter by a receiver-side duobinary shaping is used in order to apply conventional digital signal processing (DSP) schemes.
To recover data from received symbols, receivers often use timing synchronization to generate samples at integer multiples of time T (T represents the symbol period). The symbol samples are thus taken at the center of the symbol. However, when the received symbols have gone through aggressive or strong filtering, the data samples may suffer from intersymbol interference (ISI).
One approach to combat the ISI is to recover data by sampling at times that are integer multiples of T/2. Thus, even integer multiples of the sampling time instances fall on multiples of T and odd integer multiples of the sampling instances are in between adjacent symbols. Experiments and simulations performed by the present inventors have shown that the samples thus obtained show higher tolerance and robustness to the strong filtering. The constellation of the T/2 samples appears to be 9-QAM like for spectrum shaped QPSK signal. In some embodiments, a blind radius-directed equalizer (RDE) may be used to recover data from the 9-QAM like symbol spectrum. However this scheme is complicated because three constellation circles are considered. Also when the OSNR is lower (e.g., below 20 dB), it is difficult to separate three circles at acceptable bit error rate (BER), as explained in greater detail below.
In some embodiments, a decision-directed least radius distance (DD-LRD) algorithm for blind equalization is used for 9-QAM recovery and ISI suppression. Experimental results with 112 Gb/s spectrum shaped polarization multiplexing-quadrature phase shift keying (PM-QPSK) signal with a 25 GHz bandwidth WSS in the Nyquist WDM (NWDM) channel are further discussed below. As described below, the final equalized signal is detected by maximum likelihood sequence detection (MLSD) for data BER measurement.
FIG. 1A is a block diagram representation of an optical communication system 100 where the subject technology of this document can be implemented. An optical transmitter 102 transmits optical signals through an optical network 104 to one or more optical transceivers 106. The transmitted optical signals may go through intermediate optical equipment such as amplifiers, repeaters, switch, etc., which are not shown in FIG. 1A for clarity.
FIG. 1B illustrates an optical network 150 in which performance of the subject technology can be measured. A 28 Gbaud QPSK signal is generated by an I/Q modulator 152. Following this, an optical band-pass filter 154 of 4th order Gaussian type is utilized to shape the spectrum of QPSK signal. The 3 dB filter bandwidth is emulated from 22 to 30 GHz in simulation, so that the spectrum is significantly compressed to approach Nyquist bandwidth. The amplified spontaneous emission (ASE) noise (emulated as Gaussian-distributed optical white noise) is added, at 156, before optical homodyne coherent detection 158. The signal OSNR is 30 dB which is defined in a 0.1 nm noise bandwidth. Both continuous wave (CW) laser sources in the transmitter and for local oscillator in coherent receiver are of practically zero Hz linewidth. The signal is sampled at the receiver side at twice the baud rate (symbol rate).
FIG. 2 illustrates two groups of samples. The sampled signal is divided into two groups: one is the signal sampled at the time of T, i.e., the timing phase is 0 radians (202), and the other is the signal sampled at the time of T/2, i.e., the timing phase is π/2 offset (204). The T samples are called ST and T/2 samples as ST/2.
FIG. 3 shows the constellations of the received signal with the optical Gaussian filter of 28 GHz (302) and 24 GHz (306), plotted with horizontal axis 310 and vertical axis 308 representing possible data bit values. It may be seen that the constellations for the T samples are like 4-QAM (314) and T/2 samples are like 9-QAM (312). The effect by ISI after 4th order Gaussian filter can be seen in graph 306. It can be seen that each constellation point 312, 314 becomes a square-like distribution. Comparing 28 and 24 GHz filtering (graphs 302 and 306 respectively), the T samples 314 have much larger ISI when 24 GHz filter is applied, while the T/2 samples are nearly unchanged (312). In other words, a technique can recover the 9-QAM like ST/2 samples instead of conventional processing to ST samples in the presence of strong filtering.
In some embodiments, when ST is given, ST/2 can be estimated approximately by linear interpolation which is expressed as




  
    
      
        
          S
          
            T
            /
            2
          
        
        =
        
          
            
              
                
                  S
                  T
                
                ⁡
                
                  (
                  k
                  )
                
              
              +
              
                
                  S
                  T
                
                ⁡
                
                  (
                  
                    k
                    +
                    1
                  
                  )
                
              
            
            2
          
          ⁢
          
            
          
          ⁢
          
            (
            
              
                k
                =
                1
              
              ,
              2
              ,
              3
              ,
              ...
            
            )
          
        
      
    
    
      
        (
        1
        )
      
    
  




In some implementations, a higher order lowpass filter may also be used for the above-discussed interpolation.
The symbol samples ST/2 is a duobinary shaped signal and to some extent exhibits properties similar to that of a narrow bandwidth of the duobinary coded signal. Therefore T/2 samples have a greater filtering tolerance (compared to T samples) and a lower ISI due to the partial response nature of the system. However, because the constellation of the T/2 samples contains points closer to each other, higher OSNR is preferred than the corresponding T sample constellation, since multilevel detection is applied to the 9-QAM signal. In some implementations, the signal can be detected with MLSD algorithm which makes use of the inherent intersymbol memory and minimizes the number of error by selecting the most probable Trellis path.
In order to decrease the signal ISI, an adaptive finite impulse response (FIR) filter can be applied (e.g., a 9-tap filter). The CMA is normally used to blindly update the FIR tap weights for the QPSK signal which has constant modulus. However, performance of CMA deteriorates with increase in ISI. A post-filter with 2 taps (i.e. one symbol delay and add), which performs the function of duobinary shaping, can be used to further combat the strong filtering limitation.
With reference to FIG. 4, graph 406 shows example T samples suffering from severe ISI. Processing the constellation shown in graph 406 with CMA (408) and post-filter (410) can reduce the ISI to a certain degree.
Alternatively, the 9-QAM T/2 samples can be recovered directly. CMA is not applicable in a straightforward manner to 9-QAM signals because 9-QAM signals do not have constant symbol amplitude and therefore error signal obtained by CMA cannot approach to zero. Instead, the DD-LRD algorithm can be used to update the filter tap weights, which is much more accurate than CMA. The error function e(n) is given by

e(n)=y(n)(|{circumflex over (d)}(n)|2−|y(n)|2),  (2)

where y(n) is the equalized signal, {circumflex over (d)}(n) is the decided symbol. The filter tap weights updating equalization is

w(n)=w(n−1)+μe(n)x(n)*,  (3)

where w(n) is the adaptive FIR filter, μ is the convergence parameter and * is the complex conjugation operation. The DD-LRD has a superior tolerance to ISI because of it is phase independent. Besides, DD-LRD has rapid convergence speed which is robust to time varying situation. The graph 402 shows the T/2 samples suffering from severe ISI. The DD-LRD algorithm presents excellent performance to reduce the effect of ISI, as shown in graph 404.
FIG. 5 is a graph 500 that shows measured mean squared error (MSE), which is defined as:




  
    
      1
      n
    
    ⁢
    
      
        ∑
        
          i
          =
          1
        
        n
      
      ⁢
      
        
          (
          
            
              
                d
                ^
              
              ⁡
              
                (
                i
                )
              
            
            -
            
              y
              ⁡
              
                (
                i
                )
              
            
          
          )
        
        2
      
    
  
  ,




where n is the number of samples, obtained from the different digital signal processing algorithms discussed before. The curve 510 shows performance of T sample processing. The curve 512 performance of T sample processing when CMA is used. The curve 514 shows performance of T sample processing when both CMA and post-filtering are used. It can be seen that approximately 1 dB MSE improvement is obtained using CMA (512) comparing to T sampling signal (510) and further 0.5 dB is gained using post-filter with 26 GHz optical filtering (514). However filter bandwidth below 26 GHz (strong filtering) results in even more serious ISI. Therefore, MSE degrades rapidly as a function of narrower bandwidth (strong filtering) and the benefit of ISI compression from CMA decreases to less than 0.5 dB when the filter bandwidth is 22 GHz. The curves 506 and 508 present the processing to T/2 samples of the 9-QAM constellation. A 1.4 dB MSE improvement is obtained using DD-LRD (508) comparing to T/2 sampling signal (506). In one advantageous aspect, there is negligible change in the MSE performance as a function of filter bandwidth varying from 22 to 30 GHz. At wideband filtering (e.g., above 28 GHz), the performance of the DD-LRD I (508) is quite comparable to the CMA+post-filter algorithm (514). In other words, recovering data samples using 9-QAM signal instead of QPSK T samples is comparable is more robust to strong filtering and shows comparable performance at weaker filtering.

FIG. 6 shows a test system 600 of 28 Gbaud Narrow Wavelength Division Multiplexing (NWDM) PM-QPSK with WSS spectrum shaping. The NWDM subchannels 602 are generated from a comb generator based on phase and intensity modulators with 25 GHz carrier spacing and equal tone power. The odd and even channels are separated by using a 25/50 GHz optical interleaver (IL) 604. The 28 Gb/s binary electrical signals are generated from the two channels pulse pattern generator (PPG) with a pseudorandom binary sequence (PRBS) length of 211−1. The optical QPSK signal is generated using a commercial I/Q modulator 606. The polarization multiplexing is realized by the polarization multiplexer 608, which comprises a polarization maintaining coupler (PMC) to halve the signal, an optical delay line (ODL) to provide a delay of 150 symbols, and a polarization beam combiner (PBC) to recombine the signal. The even and odd channels are modulated and polarization multiplexed individually. After that, they are combined with 25 GHz channel spacing and the spectrum of each channel is shaped by a waveshaper (i.e., WSS) 612. A gain adjuster 614 is optionally used. The WSS was measured with 3-dB bandwidth of 21.6 GHz, 10-dB bandwidth of 30 GHz, and 20-dB bandwidth of 37.1 GHz. At the receiver, one tunable optical band-pass filter (OBPF) 616 with 3 dB bandwidth of 0.4 nm is employed to choose the measured subchannel. Polarization diversity homodyne detection is utilized at the receiver. The linewidth of external cavity lasers (ECL) at the transmitter and for local oscillator (LO) 618 at the receiver are both smaller than 100 kHz. A 50-GSa/s Analog/Digital conversion (ADC) 620 sampling is used in the Oscilloscope. The received data is then recovered by offline digital processed by a computer.
The received signal is resampled to 4 times symbol rate in order to process the timing phase estimation with square timing method (626). 2 samples/symbol is processed by cubic interpolation according to the right extracted clock. Four 17-taps T/2-spaced adaptive butterfly FIR filters are applied for polarization demultiplexing 628. The filters' weights are first updated by CMA for pre-convergence. The final adaptation is switched to DD-LRD 634 for precise feedback control. The adaptive FIR filter and DD-LRD in the meanwhile play a role of reducing the ISI and interchannel crosstalk. The frequency offset estimation (FOE) 630 is based on the fast Fourier transform (FFT) method and carrier phase estimation (CPE) 632 is based on blind phase search (BPS) algorithm. Finally the signal is detected by MLSD in use of intersymbol memory for data BER measurement (636). As a comparison, CMA plus post-filter scheme as well as MLSD detection is also evaluated.
The measured BERs based on different algorithms as a function of OSNR are shown in graph 700 of FIG. 7. The OSNR (horizontal axis 702) is measured in a 0.1 nm noise bandwidth. The data that is only processed with CMA (without MLSD decision) has poor performance due to the severe ISI (curve 708). This technique could not achieve 10−3 BER in our measurements. We can observe the obvious ISI noise in graphs 806, 808 (X and Y polarization) of FIG. 8. The OSNR for BER at 10−3 based on the above-discussed 9-QAM recovery algorithm is 18 dB (curve 706). The OSNR tolerance is 0.5 dB better than using CMA plus post-filter (curve 710). The constellations recovered from 9-QAM and CMA+post-filter methods are shown in FIG. 8, graphs 802, 804 (X and Y polarization) and 810, 812 (X and Y polarization) respectively.
The robustness to strong filtering and superiorly rapid convergence speed make the above-discussed 9-QAM digital recovery scheme more attractive and valuable for practical applications in optical networks.
FIG. 9 is a flowchart representation of a process 900 of optical communication by which data can be recovered in an optical receiver. At 902, an optical signal modulated using a quadrature phase shift keying (QPSK) modulation scheme is received. At 904, the received optical signal is sampled to generate a group of samples comprising samples taken at odd integer multiples of T/2, where T represents a symbol time interval of the QPSK modulation scheme. At 906, the received constellation symbols are estimated from the group of samples using a decision directed least radius distance (DD-LRD) algorithm for blind equalization. As previously described, the DD-LRD algorithm includes iteratively refining symbol estimates using an adaptive filter. At 908, data bits transmitted in the optical signal are recovered from the estimated constellation symbols.
As previously described, in some implementations, the constellation symbols are processed as 9-QAM symbols. In some implementations, the optical signal may include dual polarized QPSK signals. The dual polarizations may be orthogonal to each other (e.g., X and Y planes). In some implementations, the adaptive filter may comprise an FIR filter (e.g., the previously described 9-tap filter). In some implementations, the adaptive filter is iteratively refined using an LSE error criterion.
FIG. 10 is a block diagram representation of an optical receiver apparatus 1000. The module 1002 is for receiving an optical signal that is modulated using a quadrature phase shift keying (QPSK) modulation scheme. The module 1004 is for sampling the received optical signal to generate a group of samples comprising samples at time instances that are odd integer multiples of T/2, where T represents a symbol time interval of the QPSK modulation scheme. The module 1006 is for estimating received constellation symbols from the group of samples using a decision-directed least radius distance (DD-LRD) algorithm for blind equalization, wherein the DD-LRD algorithm includes iteratively refining symbol estimates using an adaptive filter. The module 1008 is for recovering data bits transmitted in the optical signal from the estimated constellation symbols. The apparatus 1000 and modules 1002, 1004, 1006, 1008 may further be configured to implement some of the techniques disclosed in this document.
In some embodiments, an optical communication receiver includes an optical front end that receives a polarized modulation quadrature phase shift keying (PM-QPSK) signal, a timing synchronizer that produces a sampling signal indicative of sampling instances at odd multiples of half of a baud rate of the received PM-QPSK, a signal sampler that produces sample values of the received PM-QPSK at the sampling instances, a symbol generator that estimates received symbols from the sample values using a decision directed adaptive filtering method in which decisions about received symbols are iteratively improved using a feedback signal indicative of an error of calculation of previous decisions, and a data recoverer that recovers data bits transmitted in the optical signal from the estimated received symbols.
FIG. 11 is a flowchart representation of a process 1100 of processing digital signals during the above-discussed data recovery procedures. At 1102, the previously described retiming may be performed to obtain sampling locations at odd multiples of T/2, using an estimated or a known baud rate of incoming signal. At 1104, the above-described constant modulus algorithm may be applied. Frequency offset compensation may be performed at 1106. Carrier phase estimation may be performed at 1106. The above-described DD-LRD estimation symbol estimation may be performed at 1108. Based on results obtained in 1108, a feedback may be provided to the CMA pre-convergence algorithm to refine the CMA step. Finally, at 1110, MLSE may be used to finalize symbol estimates from the 9 QAM constellation. In addition to, or in place of the algorithms described above, other methods such as the Gardner timing method for the retiming part may be used.
One of skill in the art will appreciate that a 9-QAM data recovery for PM-QPSK signal in presence of strong filtering to achieve Nyquist bandwidth is disclosed. In some embodiments, the DD-LRD algorithm for blind equalization is used for 9-QAM recovery and ISI compression. It will further be appreciated that the proposed techniques are robust under strong filtering to recover 9-QAM signal compared to the traditional QPSK data recovery techniques. It will further be appreciated that an experiment that includes a 112 Gb/s spectrum shaped PM-QPSK signal by WSS in a 25 GHz channel spacing Nyquist WDM shows the robustness of the disclosed technology. The final equalized signal is detected by MLSD for data BER measurement. 0.5 dB OSNR tolerance is improved at a BER of 10-3 comparing to CMA plus post-filter algorithm.
One of skill in the art will appreciate that the disclosed techniques are application to any communication system, and especially optical communication system, in which bandwidth efficiency approaching Nyquist limit (i.e., bandwidth of electrical signal close to 1⁄2 symbol rate and bandwidth of optical signal close to symbol rate) can be achieved by recovering transmitted symbols using T/2 located samples.
Further features and aspects are disclosed in the attached Appendix A which forms an integral part of the present document.
The disclosed and other embodiments, modules and the functional operations described in this document (e.g., an optical front end, a timing synchronizer, a signal sampler, a symbol generator and a data recoverer) can be implemented in digital electronic circuitry, or in computer software, firmware, or hardware, including the structures disclosed in this document and their structural equivalents, or in combinations of one or more of them. The disclosed and other embodiments can be implemented as one or more computer program products, i.e., one or more modules of computer program instructions encoded on a computer readable medium for execution by, or to control the operation of, data processing apparatus. The computer readable medium can be a machine-readable storage device, a machine-readable storage substrate, a memory device, a composition of matter effecting a machine-readable propagated signal, or a combination of one or more them. The term “data processing apparatus” encompasses all apparatus, devices, and machines for processing data, including by way of example a programmable processor, a computer, or multiple processors or computers. The apparatus can include, in addition to hardware, code that creates an execution environment for the computer program in question, e.g., code that constitutes processor firmware, a protocol stack, a database management system, an operating system, or a combination of one or more of them. A propagated signal is an artificially generated signal, e.g., a machine-generated electrical, optical, or electromagnetic signal, that is generated to encode information for transmission to suitable receiver apparatus.
A computer program (also known as a program, software, software application, script, or code) can be written in any form of programming language, including compiled or interpreted languages, and it can be deployed in any form, including as a stand alone program or as a module, component, subroutine, or other unit suitable for use in a computing environment. A computer program does not necessarily correspond to a file in a file system. A program can be stored in a portion of a file that holds other programs or data (e.g., one or more scripts stored in a markup language document), in a single file dedicated to the program in question, or in multiple coordinated files (e.g., files that store one or more modules, sub programs, or portions of code). A computer program can be deployed to be executed on one computer or on multiple computers that are located at one site or distributed across multiple sites and interconnected by a communication network.
The processes and logic flows described in this document can be performed by one or more programmable processors executing one or more computer programs to perform functions by operating on input data and generating output. The processes and logic flows can also be performed by, and apparatus can also be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application specific integrated circuit).
Processors suitable for the execution of a computer program include, by way of example, both general and special purpose microprocessors, and any one or more processors of any kind of digital computer. Generally, a processor will receive instructions and data from a read only memory or a random access memory or both. The essential elements of a computer are a processor for performing instructions and one or more memory devices for storing instructions and data. Generally, a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto optical disks, or optical disks. However, a computer need not have such devices. Computer readable media suitable for storing computer program instructions and data include all forms of non volatile memory, media and memory devices, including by way of example semiconductor memory devices, e.g., EPROM, EEPROM, and flash memory devices; magnetic disks, e.g., internal hard disks or removable disks; magneto optical disks; and CD ROM and DVD-ROM disks. The processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.
While this patent document contains many specifics, these should not be construed as limitations on the scope of an invention that is claimed or of what may be claimed, but rather as descriptions of features specific to particular embodiments. Certain features that are described in this document in the context of separate embodiments can also be implemented in combination in a single embodiment. Conversely, various features that are described in the context of a single embodiment can also be implemented in multiple embodiments separately or in any suitable sub-combination. Moreover, although features may be described above as acting in certain combinations and even initially claimed as such, one or more features from a claimed combination can in some cases be excised from the combination, and the claimed combination may be directed to a sub-combination or a variation of a sub-combination. Similarly, while operations are depicted in the drawings in a particular order, this should not be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed, to achieve desirable results.
Only a few examples and implementations are disclosed. Variations, modifications, and enhancements to the described examples and implementations and other implementations can be made based on what is disclosed.

## Claims

1. A method of recovering data in an optical communication system, the method comprising:
receiving an optical signal that is modulated using a quadrature phase shift keying (QPSK) modulation scheme;
sampling the received optical signal only once per symbol to generate a group of samples comprising samples at time instances that are odd integer multiples of T/2 instead of samples at time instances that are multiples of T, where T represents a symbol time interval of the QPSK modulation scheme;
estimating received constellation symbols from the group of samples using a decision-directed least radius distance (DD-LRD) algorithm for blind equalization, wherein the DD-LRD algorithm includes iteratively refining symbol estimates using an adaptive filter; and
recovering data bits transmitted in the optical signal from the estimated constellation symbols.

2. The method of claim 1, wherein the constellation symbols are processed as 9-Quadrature Amplitude Modulation (9-QAM) symbols.

3. The method of claim 1, wherein the optical signal comprises at least two polarized components.

4. The method of claim 1, wherein:
the adaptive filter comprises a finite impulse response (FIR) filter.

5. The method of claim 4, wherein the adaptive filter is refined using a least square error (LSE) optimization criterion.

6. The method of claim 4, wherein the DD-LRD algorithm includes a step of updating tap weights of the FIR filter.

7. The method of claim 6, wherein the step of updating tap weights of the FIR filter includes updating the tap weights using:


w(n)=w(n−1)+μe(n)x(n)*,

where w(n) represents the FIR filter, μ is a convergence parameter, * is a complex conjugation operation, e(n) is an error function, x(n) is an equalized signal, and n is an increment number.

8. The method of claim 1, wherein the recovering data bits comprises performing a frequency offset compensation using a fast Fourier transform (FFT) computation.

9. The method of claim 1, wherein the recovering data bits comprises performing a carrier phase estimation based on a blind phase search algorithm.

10. The method of claim 1, wherein the recovering data bits comprises detecting final equalized signal by performing a Maximum Likelihood Sequence Detection (MLSD) for data bit error rate (BER) measurement.

11. An optical communication receiver, comprising:
an optical front end that receives a polarized modulation quadrature phase shift keying (PM-QPSK) signal;
a timing synchronizer that produces a sampling signal indicative of sampling instances once a symbol period at odd multiples of half of a baud rate of the received PM-QPSK instead of sampling instances at multiples of a baud rate of the received PM-QPSK;
a signal sampler that produces sample values of the received PM-QPSK at the sampling instances, thereby generating one sample value per received symbol;
a symbol generator that estimates received symbols from the sample values using a decision directed adaptive filtering method in which decisions about received symbols are iteratively improved using a feedback signal indicative of an error of calculation of previous decisions;
a data recoverer that recovers data bits transmitted in the optical signal from the estimated received symbols.

12. The receiver of claim 11, wherein the symbol generator produces 9-Quadrature Amplitude Modulation (9-QAM) symbols.

13. The receiver of claim 12, further including:
an equalizer that uses a decision-directed least radius distance (DD-LRD) algorithm.

14. The receiver of claim 11, wherein:
the adaptive filtering method uses a 9-tap finite impulse response (FIR) adaptive filter.

15. The receiver of claim 14, wherein the adaptive filtering method uses a least square error (LSE) optimization criterion.

16. A non-transitory computer-readable medium having code stored thereon, the code, when executed by a processor, causing the processor to implement a method of recovering data in an optical communication system, the method comprising:
receiving an optical signal that is modulated using a quadrature phase shift keying (QPSK) modulation scheme;
sampling the received optical signal only once per symbol to generate a group of samples comprising samples at time instances that are odd integer multiples of T/2 instead of samples at time instances that are multiples of T, where T represents a symbol time interval of the QPSK modulation scheme;
estimating received constellation symbols from the group of samples using a decision-directed least radius distance (DD-LRD) algorithm for blind equalization, wherein the DD-LRD algorithm includes iteratively refining symbol estimates using an adaptive filter; and
recovering data bits transmitted in the optical signal from the constellation symbols.

17. The non-transitory computer-readable medium of claim 16, wherein the constellation symbols are processed as 9-Quadrature Amplitude Modulation (9-QAM) symbols.

18. The non-transitory computer-readable medium of claim 16, wherein the optical signal comprises at least two polarized components.

19. The non-transitory computer-readable medium of claim 16, wherein:
the adaptive filter comprises a finite impulse response (FIR) filter.

20. The non-transitory computer-readable medium of claim 19, wherein the adaptive filter is refined using a least square error (LSE) optimization criterion.