Reconstruction of two polarization components of an optical signal field

<u>Cross-Reference To Related Applications</u>

0001This application claims benefit of United States provisional patent application serial number <patcit id="pcit0001" dnum="US60993823B"><text>60/993823, filed September 14, 2007</text></patcit>.

<u>Technical Field</u>

0002This invention relates to the reconstruction of the two polarization components of an optical signal field and compensation for polarization-mode dispersion.

<u>Background of the Invention</u>

0003As is well known, an optical signal may have two orthogonal polarization states, each of which may have different properties. Sometimes such polarization states are intentionally introduced, such as in creating a polarization-multiplexed signal in which the two orthogonal polarization states of the optical carrier are arranged so that each carries different data in order to double the spectral efficiency. Such a polarization-multiplexed signal has two so-called "generic" polarization components, each of which carries a single data modulation. Note that by a generic polarization component it is generally intended the signal at the point at which the modulation of that polarization component is completed. It should be appreciated that each generic polarization component may initially, or otherwise, exist separate from the other generic polarization component with which it is later combined.

0004The polarization orientations of the generic signal components are generally changed by the birefringence of the fiber, and possibly other fiber properties, during the passage of the signal over the optical path. Such changes may be time varying because at least the fiber birefringence is typically a function of various factors such as ambient temperature, mechanical stress, and so forth, which may vary over time and be different at various points of the transmission path. As a result, the polarization orientation of each of the generic signal components is generally unknown at the receiver.

0005Sometimes, undesirably, the fiber birefringence is so large that polarization-mode dispersion (PMD) is caused, i.e., a generic optical signal component is decomposed into two orthogonal polarization components along the two principle state of polarization (PSP) axes of the fiber, along one of which the light travels at its fastest speed through the fiber and along the other of which the light travels at its slowest speed through the fiber. In such a case, not only may the phase relationship between the two polarization components be time varying, but also each of the two orthogonal polarization components may arrive at the receiver at different times due to the PMD-induced differential group delay (DGD) between the two PSP axes. Note that, actually, as suggested above, each small section of the fiber behaves as if it is its own mini fiber that introduces its own DGD between the two PSP axes. However, for simplification purposes, one may treat the fiber as a single DGD element that introduces a certain DGD between the two axes, based on a first order approximation of the PMD. Thus, for a particular fiber or optical link, PMD is a stochastic effect, and the PMD-induced DGD may also be time varying.

0006Other linear effects distort optical signals transmitted over optical fibers. Such effects include chromatic dispersion (CD). Optical compensation methods are typically employed to reduce signal distortion that arises due to CD or PMD. Electronic chromatic dispersion compensation (EDC) has recently emerged as a technique that can flexibly reduce the distortion induced by CD in a cost effective manner. As explained by <nplcit id="ncit0001" npl-type="s"><text>M. S. O'Sullivan, K. Roberts, and C. Bontu, in "Electronic dispersion compensation techniques for optical communication systems," ECOC'05, paper Tu3.2.1, 2005</text></nplcit>, EDC can be performed at the transmitter. Alternatively, EDC can be performed at the receiver. As described by <nplcit id="ncit0002" npl-type="s"><text>S. Tsukamoto, K. Katoh, and K. Kikuchi, in "Unrepeated Transmission of 20-Gb/s Optical Quadrature Phase-Shift- Keying Signal Over 200-km Standard Single-Mode Fiber Based on Digital Processing of Homodyne-Detected Signal for Group-Velocity Dispersion Compensation," IEEE Photonics Technology Letters, Volume 18, Issue 9, 1 May 2006, pp. 1016 - 1018</text></nplcit>, EDC is implemented with a coherent-detection receiver. In addition, EDC can be implemented with a special direct differential detection receiver as explained by X. Liu and X. Wei, in United States Patent Application Serial No. <patcit id="pcit0002" dnum="US11525786B"><text>11/525786</text></patcit> entitled "Reconstruction and Restoration Of Optical Signal Field", filed on September 22, 2006 and assigned to Lucent Technologies, which shall be referred to hereinafter as Liu-Wei.

0007Unlike CD, PMD in a fiber link may change very rapidly and PMD compensation usually has to be done in the receiver. Electronic PMD compensation (EPMDC) has also attracted attention recently for its potential cost effectiveness. As explained by <nplcit id="ncit0003" npl-type="s"><text>J. Hong, R. Saunders, and S. Colaco, in "SiGe equalizer IC for PMD Mitigation and Signal Optimization of 40 Gbits/s Transmission", published in Optical Fiber Communication Conference 2005, paper OW02</text></nplcit>. However, the capability of the EPMDC with a conventional direct-detection receiver is quite limited in that the improvement in PMD tolerance is usually only about 50%.

0008<nplcit id="ncit0004" npl-type="s"><text>T. Pfau et al, "PDL-Tolerant Real-time Polarization-Multiplexed QPSK Transmission with Digital Coherent Polarization Diversity Receiver", 2007 Digest of the Leos Summer Topical Meetings, July 2007, pages 17-18</text></nplcit>, is concerned with an implementation of a real-time electronic polarization tracking algorithm.

0009<nplcit id="ncit0005" npl-type="s"><text>I. Bar-David, "Direct differential detection of phase-shift-keyed signals: a local-oscillatorless DPSK receiver", IEE Proceedings - Optoelectronics, Vol. 141, No. 1, February 1994, pages 38-42</text></nplcit>, is concerned with an application of a phase-diversity mixing concept to a direct differential detection of phase-shift-keyed (PSK) transmissions.

0010<patcit id="pcit0003" dnum="WO2006034975A1"><text>WO 2006/034975 A1</text></patcit> discloses a method and arrangement for optical transmission of data signals by means of differential phase modulation in a polarization multiplex method.

<u>Summary of the Invention</u>

0011It is an object of the invention to provide an improved optical receiver, and an improved method for use in an optical receiver.

0012This object is achieved by the optical receiver and the method for use in an optical receiver according to the independent claims. Further embodiments of the invention are part of the dependent claims.

0013In accordance with the principles of the invention, a digital version of the complex field, i.e., both amplitude and phase, e.g., with respect to a reference point, of each of two orthogonal polarization components of a received optical signal are developed at a receiver by employing a dual-polarization direct differential receiver portion that uses direct differential detection to develop a digital representation of optical signals derived from each of two orthogonal polarization components of a received optical signal which are then processed using digital signal processing (PSP) to develop a digital representation of an intensity and a phase profile representing the polarizations as received at the receiver. The reconstructed digital versions of the complex field of each of the two orthogonal polarization components of the optical signal as received at the receiver are then be further processed jointly to develop at least one so-called "generic" polarization component of the received optical. Note that due to fiber birefringence or PMD, the two orthogonal polarization components of the optical signal as received at the receiver after fiber transmission are generally not the generic polarization components of the signal. In accordance with an aspect of the invention, during the joint processing the relative phase difference between the two reference points used in the two reconstructed optical fields is determined. This may be achieved by employing a searching technique.

0014In one embodiment of the invention, a polarization beam splitter (PBS) is first used to separate the received optical signal into two arbitrarily orthogonal polarization components, <i>E_x</i> and <i>E_y'</i>. Each of the orthogonal polarization components is supplied to a special direct differential detection receiver, which employs a special pair of optical delay interferometers (ODIs) with a phase delay difference of about π/2, such as are described in Liu-Wei and which are herein referred to as an I/Q ODI pair. At least the four outputs of each I/Q ODI pair are then detected by two balanced detectors, whose two outputs are sampled by respective analog to digital converters (ADCs), are then processed to obtain a digital representation of the received signal optical field along the corresponding polarization axis, i.e., <i>x'</i> or <i>y'</i>, according to Liu-Wei.

0015In a second embodiment of the invention, the received optical signal is supplied directly into a single polarization-independent I/Q ODI pair and the resulting four outputs are each connected to a respective associated one of four PBSs, all of which have the same polarization orientation. Each of the PBS produces two outputs, so that in total there are eight outputs from the four PBSs, consisting of four outputs derived from the first polarization, e.g., <i>x</i>'-polarized outputs, and four outputs derived from the second polarization, e.g., y'-polarized outputs. Each pair of outputs of the PBSs that corresponds to a single optical delay interferometer and a single polarization are supplied to a respective one of four balanced detectors, whose outputs are sampled by a respective one of four corresponding ADCs. Each of the resulting sampled waveforms are then processed to obtain a digital representation of the received signal's optical field along each of the polarization axes <i>x'</i> and <i>y'</i> in the manner described in Liu-Wei.

0016Even though this second embodiment requires the use of three additional PBSs as compared to the first embodiment, due to the relative cost of I/Q ODIs themselves and the control electronics associated therewith, as compared to the cost of PBSs, advantageously, because only one I/Q ODI pair is employed, significant cost savings can be achieved. Furthermore, the second embodiment may be more compactly implemented.

0017Either embodiment of the invention may be implemented with free space or fiber based optics, or any combination thereof.

0018Although to save cost it is expected that implementers generally will approximate the intensity profile of one or more of the polarization components of the received signal from the absolute value of their respective complex waveforms, they may instead employ direct intensity detection to obtain a more accurate measurement of the intensity profile.

0019The techniques of the instant invention are suitable to be employed with various types of optical differential phase-shift keying (DPSK) signals, such as differential binary phase-shift keying (DBPSK) and differential quadrature phase-shift keying (DQPSK) signals. They may also be employed with amplitude-shift keying (ASK), combined DPSK/ASK, and quadrature amplitude modulation (QAM).

<u>Brief Description of the Drawing</u>

0020In the drawing: <ul id="ul0001" list-style="none" compact="compact"><li><figref idref="f0001">FIG. 1</figref> shows an exemplary apparatus, typically in a receiver, arranged in accordance with the principles of the invention, for developing at least one so-called "generic" polarization component of a received optical signal;</li><li><figref idref="f0002">FIG. 2</figref> shows an embodiment of the invention similar to that shown in <figref idref="f0001">FIG. 1</figref> but in which the intensity detection branches are omitted;</li><li><figref idref="f0003">FIG. 3</figref> shows another embodiment of the invention similar to that shown in <figref idref="f0001">FIG. 1</figref> but which employs only a single I/Q ODI pair;</li><li><figref idref="f0004">FIG. 4</figref> shows an exemplary polarization "evolution" as an optical signal passes over a typical fiber transmission link that causes PMD;</li><li><figref idref="f0005">FIG. 5</figref> shows an exemplary arrangement to perform the digital signal processing needed to recover the generic polarization components from the reconstructed optical fields of the two orthogonal polarization components of the received optical signal, in accordance with an aspect of the invention;</li><li><figref idref="f0005">FIG. 6</figref> shows an exemplary high level block diagram of an arrangement which is suitable to be used to make up each of the processing units of <figref idref="f0005">FIG. 5</figref>;</li><li><figref idref="f0006">FIG. 7</figref> shows an exemplary process expressed in flow-chart form for finding the best guesses of the parameters needed to recover E_x(t) and E_y(t), in accordance with an aspect of the invention;</li><li><figref idref="f0006">FIG. 8</figref> shows an exemplary process, expressed in flow-chart form, which is performed by the real-time subprocessor 602 of <figref idref="f0005">FIG. 6</figref> in one embodiment of the invention;</li><li><figref idref="f0007">FIG. 9</figref> shows an exemplary high level block diagram arrangement that is suitable to be used to make up each of the processing units of <figref idref="f0005">FIG. 5</figref> but which is arranged to speed up the processing as compared to the arrangement shown in <figref idref="f0005">FIG. 6</figref>;</li><li><figref idref="f0007">FIG. 10</figref> shows exemplary high level block diagram of an arrangement suitable to be used to make up each of the processing units of <figref idref="f0005">FIG. 5</figref> but which is arranged to treat the fiber as if it was made up of multiple segments, so as to achieve better compensation for PMD, including higher order PMD, than can be achieved than using the arrangements of <figref idref="f0005">FIGs. 6</figref> or <figref idref="f0007">9</figref>;</li><li><figref idref="f0008">FIG. 11</figref> shows high level block diagram of an arrangement suitable to be used to make up each of the processing units of <figref idref="f0005">FIG. 5</figref> but which is arranged to treat the fiber as if it was made up of 2 segments each having DGD values each fixed to 0.4<i>T_S</i></li><li><figref idref="f0008">FIG. 12</figref> shows an exemplary high level block diagram of an arrangement which is suitable to be used to make up each of the processing units of <figref idref="f0005">FIG. 5</figref> but which is arranged for use with a coherent-detection receiver, in accordance with an aspect of the invention; and</li><li><figref idref="f0009">FIG. 13</figref>, shows an exemplary process, expressed in flow-chart form, that is performed in the real-time subprocessor of <figref idref="f0008">FIG. 12</figref>, in one embodiment of the invention.</li></ul>

<u>Detailed Description</u>

0021The following merely illustrates the principles of the invention.

0022Thus, for example, it will be appreciated by those skilled in the art that any block diagrams herein represent conceptual views of illustrative circuitry embodying the principles of the invention. Similarly, it will be appreciated that any flow charts, flow diagrams, state transition diagrams, pseudocode, and the like represent various processes which may be substantially represented in computer readable medium and so executed by a computer or processor, whether or not such computer or processor is explicitly shown.

0023The functions of the various elements shown in the FIGs., including any functional blocks labeled as "processors", may be provided through the use of dedicated hardware as well as hardware capable of executing software in association with appropriate software. When provided by a processor, the functions may be provided by a single dedicated processor, by a single shared processor, or by a plurality of individual processors, some of which may be shared. Moreover, explicit use of the term "processor" or "controller" should not be construed to refer exclusively to hardware capable of executing software, and may implicitly include, without limitation, digital signal processor (DSP) hardware, network processor, application specific integrated circuit (ASIC), field programmable gate array (FPGA), read-only memory (ROM) for storing software, random access memory (RAM), and non-volatile storage. Other hardware, conventional and/or custom, may also be included. Similarly, any switches shown in the FIGS, are conceptual only. Their function may be carried out through the operation of program logic, through dedicated logic, through the interaction of program control and dedicated logic, or even manually, the particular technique being selectable by the implementor as more specifically understood from the context.

0024In the claims hereof any element expressed as a means for performing a specified function is intended to encompass any way of performing that function. This may include, for example, a) a combination of electrical or mechanical elements which performs that function or b) software in any form, including, therefore, firmware, microcode or the like, combined with appropriate circuitry for executing that software to perform the function, as well as mechanical elements coupled to software controlled circuitry, if any. The invention as defined by such claims resides in the fact that the functionalities provided by the various recited means are combined and brought together in the manner which the claims call for. Applicant thus regards any means which can provide those functionalities as equivalent as those shown herein.

0025Software modules, or simply modules which are implied to be software, may be represented herein as any combination of flowchart elements or other elements indicating performance of process steps and/or textual description. Such modules may be executed by hardware that is expressly or implicitly shown.

0026Unless otherwise explicitly specified herein, the drawings are not drawn to scale.

0027In the description, identically numbered components within different ones of the FIGs. refer to the same components.

0028<figref idref="f0001">FIG. 1</figref> shows an exemplary apparatus, typically in a receiver, arranged in accordance with the principles of the invention, for developing at least one so-called "generic" polarization component of a received optical signal. This is achieved by developing an electronic version of the entire complex optical field of a received optical signal by employing direct differential detection in conjunction with digital signal processing. It is further possible to compensate for various impairments that were inflicted upon the optical signal as it traveled from its source. <figref idref="f0001">FIG. 1</figref> shows a) polarization beam splitter (PBS) 101; b) 1x3 optical splitters 102 and 103; c) optical delay interferometers (ODIs) 105, 106, 107, and 108; d) balanced intensity detectors 111, 112, 113, and 114; e) single intensity detectors 115 and 116; f) optional amplifiers 121, 122, 123, 124, 125, and 126; g) optional automatic-gain controllers (AGCs) 131, 132, 133, 134, 135, and 136; h) analog-to-digital converters (ADCs) 141, 142, 143, 144, 145, and 146; and i) digital signal processing (DSP) unit 150.

0029More specifically, polarization beam splitter 101 separates the received optical signal to produce two orthogonal polarization components, <i>E_x'</i> and <i>E_y'</i> therefrom. <i>E_x'</i> is supplied to 1x3 optical splitter 102 while <i>E</i>_<i>y</i>' is supplied to 1x3 optical splitter 103. However, the orthogonal polarization components, <i>E_x'</i> and <i>E_y'</i> are highly unlikely to correspond to the polarization of the generic components that were originally transmitted as they are currently manifest in the received signal.

00301x3 optical splitter 102 replicates the received optical signal so as to produce three copies. One of the three beams produced by 1 x3 optical splitter 102 is supplied to optical delay interferometer (ODI) 105, another of the three beams produced by 1 x3 optical splitter 102 is supplied to ODI 106, and the last beam is supplied to photodiode 115. The optical power allotted to each of the copies from the originally input optical signal is at the discretion of the implementer. In one embodiment of the invention, the power is divided up so that about between 40 to 45 percent of the input power is supplied as output to each of ODIs 105 and 106 and the remaining power, e.g., between 10 and 20 percent, is supplied to photodiode 115.

0031As will be readily recognized by those of ordinary skill in the art, optical delay interferometers (ODIs) 105, 106, 107 and 108 may be any type of interferometer having the required characteristics. For example, the ODIs may be based on the well-known, so-called Mach-Zehnder interferometer. Alternatively, the ODIs may be based on the well-known, so-called Michaelson interferometer. Preferably, ODIs 105 and 106 are made in a pair so that their phase orthogonality (or π/2 offset in their differential phases between interfering arms) is automatically guaranteed, e.g., using techniques such as disclosed in United States Patent Application Serial Number <patcit id="pcit0004" dnum="US10875016B"><text>10/ 875016 applied for on June 23, 2004 by Christopher R. Doerr and Douglas M. Gill</text></patcit>, entitled "Apparatus and Method for Receiving a Quadrature Differential Phase Shift Key Modulated Optical Pulsetrain" published as <patcit id="pcit0005" dnum="US20050286911A"><text>2005/0286911 on December 29, 2005</text></patcit> and United States Patent Application Serial Number <patcit id="pcit0006" dnum="US11163190B"><text>11/163190 applied for on October 8, 2005 by Xiang Liu</text></patcit>, entitled "Optical Demodulating Apparatus and Method" published as <patcit id="pcit0007" dnum="US20070081826A"><text>2007/0081826 on April 12, 2007</text></patcit>. Further preferably, the two ODI pairs are either monolithically integrated on a same substrate so that their characteristic polarization orientations are the same. Note that the characteristic polarization orientations of an ODI is analogous to the the PSP of a fiber.

0032Each of the ODIs 105 and 107 has a delay of about ΔT in the optical path between its respective two arms and a phase difference, i.e., offset, of φ₀, where <maths id="math0001" num="(1)"><math display="block"><mi mathvariant="normal">Δ</mi><mo>⁢</mo><mi>T</mi><mo>=</mo><mfrac><msub><mi>T</mi><mi>S</mi></msub><mi mathvariant="italic">sps</mi></mfrac><mo>,</mo></math><img file="EP2201706B1_D0001.tif" /></maths> where <i>T_S</i> is the symbol period of the signal, <i>sps</i> is the number of samples per symbol taken by the analog to digital converters and φ₀ is an arbitrarily selected number, which is prefereably set at π/4. If so, the free spectral range (FSR), i.e., 1/ΔT, of the ODIs is related to the signal symbol rate (SR) as <i>FSR = SR · sps.</i> Note that, based on numerical simulations, it has been found that, preferably, <i>sps</i> be set to a value of 4. This is because an <i>sps</i> value of less than 4 tends to not be sufficient to accurately represent the signal waveform sufficiently given the procedures described hereinbelow, while <i>sps</i> greater than 4 provides only negligible improvement.

0033The delay difference may be achieved, in one embodiment of the invention, by adjusting one arm of the interferometer to have a gross length difference of ΔT*C/<i>n</i>, where C is the speed of light in vacuum and <i>n</i> is the index of refraction of the medium of the arm, and then adjusting the length further to cause a phase shift of φ₀. Note that in practice, because a phase shift of φ₀ corresponds to a very small length difference, the phase shift portion may actually be somewhat longer or shorter, so that the total length is φ₀ plus or minus a multiple of 2π. That way, even thought the length is hot precisely φ₀, the phase change is effectively φ₀.

0034The total length change used to achieve the effective length change of φ₀ may be some percentage of the length ΔT·C/<i>n</i>. While even up to 25 percent can work, preferably, the percentage is less than 10 percent, and of course, the more accurate the length can be made to match the actual desired length the better the performance will be. In other embodiments of the invention, the delay required may be divided between the arms, so long as the required delay and phase difference is achieved. Those of ordinary skill in the art will readily recognize how to develop an appropriate arrangement to implement ODIs 105 and 107.

0035While any value may be employed as the value of phase offset φ₀, for compatibility with conventional receivers, as will be seen hereinbelow, certain values of φ₀ may be advantageously employed. For example, a good value of φ₀ is π/4 for DQPSK and 0 for DBPSK.

0036Each of ODIs 106 and 108 are similar to ODIs 105 and 107, in that each has delay of about ΔT in the optical path between their respective two arms, but between their arms they each have a phase offset of φ₀ - π/2. Thus, the difference between the phase offsets of ODIs 105 and 106 is π/2, so ODIs 105 and 106 are said to have orthogonal phase offsets. Similarly, the difference between the phase offsets of ODIs 107 and 108 is π/2, so ODIs 107 and 108 are said to have orthogonal phase offsets.

0037Together, ODI 105 and 106 make up a so-called "I/Q ODI pair". The four outputs of I/Q ODI pair made up of ODIs 105 and 106 are then detected by two balanced detectors 111 and 112, respectively, in the manner shown in <figref idref="f0001">FIG. 1</figref>. The outputs of balanced detectors 111 and 112 are amplified by a respective one of amplifiers 121 and 122, and they may then be normalized by one of optional automatic-gain controllers (AGCs) 131 and 132.

0038Balanced intensity detectors 111 and 112 are conventional. Typically, each of balanced intensity detectors 111 and 112 is made up of a pair of well-matched photodiodes. Balanced intensity detectors 111 and 112 convert the output of each of the arms of ODIs 105 and 106 to an electrical representation. Thus, balanced intensity detectors 111 and 112 obtain an electrical version of the real and imaginary parts of the complex waveform that contains the information about the phase differences between two time locations separated by ΔT in the polarization component of the received optical signal supplied from PBS 101 to 1x3 optical splitter 102.

0039Photodiode 115 performs conventional direct intensity detection, and thus obtains the intensity profile of <i>E_x'</i> in electronic form.

0040Amplifiers 121, 122, and 125 amplify the signals supplied as outputs by balanced intensity detector 111, balanced intensity detector 112, and photodiode 115, respectively. Typically, amplifiers 121, 122, and 125 convert the current which is output by the various photodiodes of balanced intensity detector 111, balanced intensity detector 112, and photodiode 115 to respective corresponding voltages. To this end, amplifiers 121, 122, and 125 may be trans-impedance amplifiers. Furthermore, amplifiers 121 and 122 may be differential amplifiers. After amplification, each of the outputs is typically single ended. Optional automatic-gain controllers (AGCs) 131, 132, and 135 may be employed to normalize the electronic waveforms prior to digitization.

0041Analog-to-digital converters (ADCs) 141, 142, and 143 perform "digital sampling" of the amplified signals to develop a digital representation of the amplified signals. ADCs 141, 142, and 145 typically have the same resolution, e.g., 8 bits.

00421x3 optical splitter 103, similar to 1x3 optical splitter 102, replicates the received optical signal so as to produce three copies. One of the three beams produced by 1x3 optical splitter 103 is supplied to optical delay interferometer (ODI) 107, another of the three beams produced by 1x3 optical splitter 103 is supplied to ODI 108, and the last beam is supplied to photodiode 116.

0043Together, ODIs 107 and 108 make up an I/Q ODI pair. The four outputs of I/Q ODI pair made up of ODIs 107 and 108 are then detected by two balanced detectors 113 and 114, respectively, in the manner shown in <figref idref="f0001">FIG. 1</figref>. The outputs of balanced detectors 113 and 114 are amplified by a respective one of amplifiers 123 and 124, and they may then be normalized by one of optional automatic-gain controllers (AGCs) 133 and 134.

0044Balanced intensity detectors 113 and 114 are conventional. Typically, each of balanced intensity detectors 113 and 114 is made up of a pair of well-matched photodiodes. Balanced intensity detectors 113 and 114 convert the output of each of the arms of ODIs 107 and 108 to an electrical representation. Thus, balanced intensity detectors 113 and 114 obtain an electrical version of the real and imaginary parts of the complex waveform that contains the information about the phase differences between two time locations separated by ΔT in the polarization component of the received optical signal supplied from PBS 101 to 1x3 optical splitters 103.

0045Photodiode 116 performs conventional direct intensity detection, and thus obtains the intensity profile of <i>E_y'</i> in electronic form.

0046Amplifiers 123, 124, and 126 amplify the signals supplied as outputs by balanced intensity detector 113, balanced intensity detector 114, and photodiode 116, respectively. Typically, amplifiers 123, 124, and 126 convert the current which is output by the various photodiodes of balanced intensity detector 113, balanced intensity detector 114, and photodiode 116 to respective corresponding voltages. To this end, amplifiers 123, 124, and 126 may be trans-impedance amplifiers. Furthermore, amplifiers 123 and 124 may be differential amplifiers. After amplification, each of the outputs is typically single ended. Optional automatic-gain controllers (AGCs) 133, 134, and 136 may be employed to normalize the electronic waveforms prior to digitization.

0047Analog-to-digital converters (ADCs) 143, 144, and 146 perform "digital sampling" of the amplified signals to develop a digital representation of the amplified signals. ADCs 143, 144, and 146 typically have the same resolution, e.g., 8 bits.

0048Digital signal processing unit 150 receives the digital representation of all of the digitized signals supplied from ADCs 141-146 and develops at least one so-called "generic" polarization component of the received optical. Note that by a generic polarization component it is generally intended the original signal that corresponds to the received signal at the point at which the modulation of that polarization component for transmission is completed.

0049In accordance with an aspect of the invention, reconstruction unit 151-1 receives the digitized signals supplied from ADCs 141, 142, and 145 and develops a digital representation of the amplitude and phase profiles of one of the polarizations of the received optical signal, e.g., <i>x'</i>. Similarly, in accordance with an aspect of the invention, reconstruction unit 151-2 receives the digitized signals supplied from ADCs 143, 144, and 146 and develops a digital representation of the received optical signal field, i.e., the amplitude and phase profiles, of one of the other polarization of the received optical signal, e.g., <i>y'</i>. To this end, reconstruction unit 151-1 treats its inputs as if they were the entirety of the optical signal and processes those inputs according to Liu-Wei, e.g., using <i>m</i>=1, prior to any compensation for distortions, e.g., according to the processing described in Liu-Wei in connection with reconstruction unit 151 thereof. The resulting output for this reconstruction, referred to in Liu-Wei as <i>E_R</i>(<i>t_s</i>), is referred to herein as <i>E</i>_<i>x</i>'<i>(t)</i>. Similarly, reconstruction unit 151-2 treats its inputs as if they were the entirety of the optical signal and processes those inputs according to according to Liu-Wei, e.g., using <i>m</i>=1, prior to any compensation for distortions, e.g., according to the processing described in Liu-Wei in connection with reconstruction unit 151 thereof, to develop received optical signal field, i.e., the amplitude and phase profiles. The resulting output for this reconstruction, referred to in Liu-Wei as <i>E_R</i>(<i>t_s</i>), is referred to herein as <i>E_y'(t)</i>.

0050Due to fiber birefringence or PMD, the two orthogonal polarization components of the optical signal as received at the receiver after fiber transmission are generally not the generic polarization components of the signal. Therefore, in accordance with an aspect of the invention, the reconstructed digital versions of the complex field of each of the two orthogonal polarization components of the optical signal as received at the receiver <i>E_x'(t)</i> and <i>E_y'(t)</i> need to be further processed jointly to develop at least one "generic" polarization component of the received optical, as to be described below.

0051<figref idref="f0004">FIG. 4</figref> shows an exemplary polarization "evolution" as an optical signal passes over a typical fiber transmission link 402 that causes PMD. The two polarization components of the signal after fiber transmission along the two PSP axes of the fiber, defined hereafter as <maths id="math0002" num=""><math display="inline"><msubsup><mi>E</mi><mo>∥</mo><mi mathvariant="italic">out</mi></msubsup></math><img file="EP2201706B1_D0002.tif" /></maths> and <maths id="math0003" num=""><math display="inline"><msubsup><mi>E</mi><mo>⊥</mo><mi mathvariant="italic">out</mi></msubsup></math><img file="EP2201706B1_D0003.tif" /></maths> -treating the fiber as a single DGD element that introduces a certain DGD between the two axes--can be linked to the two orthogonal components of the received signal field <i>E_x'</i> and <i>E_y'</i> as follows, <maths id="math0004" num="(2)"><math display="block"><mtable columnalign="left"><mtr><mtd><mfenced open="[" close="]"><mtable><mtr><mtd><msubsup><mrow><mi>E</mi></mrow><mo>‖</mo><mi mathvariant="italic">out</mi></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><msubsup><mrow><mi>E</mi></mrow><mo>⊥</mo><mi mathvariant="italic">out</mi></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr></mtable></mfenced></mtd><mtd><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced></mtd><mtd><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced></mtd><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced></mtd></mtr></mtable></mfenced><mo>⋅</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mo>+</mo><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><msubsup><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow><mi>ʹ</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>⋅</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><msub><mi>φ</mi><mn>2</mn></msub></mrow></msup></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mspace width="1em" /></mtd><mtd><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced><mo>+</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced><mo>⋅</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><msub><mi>φ</mi><mn>2</mn></msub></mrow></msup></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced><mo>+</mo><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced><mo>⋅</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><msub><mi>φ</mi><mn>2</mn></msub></mrow></msup></mtd></mtr></mtable></mfenced></mtd></mtr></mtable><mo>,</mo></math><img file="EP2201706B1_D0004.tif" /></maths> where θ₂ is the angle between the two characteristic orientations of PBS 101 and the two PSP axes of fiber 402, and δφ₂ is the additional phase difference between the two reconstructed signal fields <i>E_x'</i> and <i>E_y'</i> as compared to the phase difference of the two received polarization components right after PBS 403. The additional phase difference includes an initially unknown relative phase difference between the two reference points used in the two reconstructed optical fields. In accordance with an aspect of the invention, during the joint processing the phase difference δφ₂ is determined by employing a searching technique, such as is described hereinbelow. <i>E_x'</i> and <i>E_y'</i> are the reconstructed optical fields of the two orthogonal polarization components of the received optical signal as separated by polarization beam splitter 403.

0052The two polarization components of the signal along the two PSP axes of fiber 402 at the input of fiber 402, defined hereafter as <maths id="math0005" num=""><math display="inline"><msubsup><mi>E</mi><mo>∥</mo><mi mathvariant="italic">in</mi></msubsup></math><img file="EP2201706B1_D0005.tif" /></maths> and <maths id="math0006" num=""><math display="inline"><msubsup><mi>E</mi><mo>⊥</mo><mi mathvariant="italic">in</mi></msubsup><mo>,</mo></math><img file="EP2201706B1_D0006.tif" /></maths> can be related to <maths id="math0007" num=""><math display="inline"><msubsup><mi>E</mi><mo>∥</mo><mi mathvariant="italic">out</mi></msubsup></math><img file="EP2201706B1_D0007.tif" /></maths> and <maths id="math0008" num=""><math display="inline"><msubsup><mi>E</mi><mo>⊥</mo><mi mathvariant="italic">out</mi></msubsup></math><img file="EP2201706B1_D0008.tif" /></maths> as <maths id="math0009" num="(3)"><math display="block"><mtable columnalign="left"><mtr><mtd><msubsup><mrow><mi>E</mi></mrow><mo>‖</mo><mi mathvariant="italic">in</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>=</mo><msubsup><mrow><mi>E</mi></mrow><mo>‖</mo><mi mathvariant="italic">out</mi></msubsup><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub></mfenced><mo>⋅</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><msub><mi>φ</mi><mn>2</mn></msub></mrow></msup><mo>,</mo></mtd></mtr><mtr><mtd><msubsup><mrow><mi>E</mi></mrow><mo>⊥</mo><mi mathvariant="italic">in</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>=</mo><msubsup><mi>E</mi><mo>⊥</mo><mi mathvariant="italic">out</mi></msubsup><mfenced><mi>t</mi></mfenced><mo>,</mo></mtd></mtr></mtable></math><img file="EP2201706B1_D0009.tif" /></maths> where <i>τ_DGD</i> is the PMD-induced DGD, and δφ is the PMD-induced or birefringence-induced phase difference between the two PSPs, which may be time varying, e.g., due to environmental, e.g., mechanical or temperature changes. Conventionally, the ∥ and ⊥ axes are called the fast PMD axis and the slow PMD axis, respectively. In the case that PMD is sufficiently small, τ<i>_DGD</i> can be approximated as 0 in Eq. (3), but the PMD-induced or birefringence-induced phase difference δφ cannot be neglected.

0053When the original signal emitted from transmitter 401 is polarization multiplexed to carry two generic polarization components, <i>E_x</i> and <i>E_y</i>, the two generic components can be linked to <maths id="math0010" num=""><math display="inline"><msubsup><mi>E</mi><mo>∥</mo><mi mathvariant="italic">in</mi></msubsup></math><img file="EP2201706B1_D0010.tif" /></maths> and <maths id="math0011" num=""><math display="inline"><msubsup><mi>E</mi><mo>⊥</mo><mi mathvariant="italic">in</mi></msubsup></math><img file="EP2201706B1_D0011.tif" /></maths> as <maths id="math0012" num="(4)"><math display="block"><mtable><mtr><mtd><mfenced open="[" close="]"><mtable><mtr><mtd><msub><mi>E</mi><mi>x</mi></msub><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><msub><mi>E</mi><mi>y</mi></msub><mfenced><mi>t</mi></mfenced></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd><mtd><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd></mtr><mtr><mtd><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd></mtr></mtable></mfenced><mo>⋅</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msubsup><mrow><mi>E</mi></mrow><mo>‖</mo><mi mathvariant="italic">in</mi></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><msubsup><mrow><mi>E</mi></mrow><mo>⊥</mo><mi mathvariant="italic">in</mi></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr></mtable></mfenced></mtd></mtr></mtable><mo>,</mo></math><img file="EP2201706B1_D0012.tif" /></maths> where θ₁ is the angle between the two orthogonal polarization components of the original signal from the transmitter and the two PSP axes of the fiber at its input.

0054Combing equations (2), (3), and (4), two generic polarization components, <i>E_x</i> and <i>E_y</i>, can then be expressed in terms of the received polarization components as <maths id="math0013" num="(5)"><math display="block"><mtable columnalign="left"><mtr><mtd><mfenced open="[" close="]"><mtable><mtr><mtd><msub><mi>E</mi><mi>x</mi></msub><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><msub><mi>E</mi><mi>y</mi></msub><mfenced><mi>t</mi></mfenced></mtd></mtr></mtable></mfenced></mtd><mtd><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd><mtd><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd></mtr><mtr><mtd><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd></mtr></mtable></mfenced><mo>⋅</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msubsup><mrow><mi>E</mi></mrow><mo>‖</mo><mi mathvariant="italic">in</mi></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><msubsup><mrow><mi>E</mi></mrow><mo>‖</mo><mi mathvariant="italic">in</mi></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mspace width="1em" /></mtd><mtd><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd><mtd><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd></mtr><mtr><mtd><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced></mtd></mtr></mtable></mfenced><mo>⋅</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msubsup><mi>E</mi><mo>‖</mo><mi mathvariant="italic">out</mi></msubsup><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub></mfenced><mo>⋅</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><mi>φ</mi></mrow></msup></mtd></mtr><mtr><mtd><msubsup><mi>E</mi><mo>⊥</mo><mi mathvariant="italic">out</mi></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced><mo>⁢</mo><msubsup><mi>E</mi><mo>‖</mo><mi mathvariant="italic">out</mi></msubsup><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub></mfenced><mo>⋅</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><mi>φ</mi></mrow></msup><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced><mo>⁢</mo><msubsup><mi>E</mi><mo>⊥</mo><mi mathvariant="italic">out</mi></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced><mo>⁢</mo><msubsup><mi>E</mi><mo>‖</mo><mi mathvariant="italic">out</mi></msubsup><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub></mfenced><mo>⋅</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><mi>φ</mi></mrow></msup><mo>+</mo><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced><mo>⁢</mo><msubsup><mi>E</mi><mo>⊥</mo><mi mathvariant="italic">out</mi></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr></mtable></mfenced></mtd></mtr><mtr><mtd><mspace width="1em" /></mtd><mtd><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced><mo>⁢</mo><mfenced open="[" close="]" separators=""><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub></mfenced><mo>+</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub></mfenced><mo>⋅</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><msub><mi>φ</mi><mn>2</mn></msub></mrow></msup></mfenced><mo>⁢</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><mi>φ</mi></mrow></msup><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced><mo>⁢</mo><mfenced open="[" close="]" separators=""><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced><mo>+</mo><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced><mo>⁢</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><msub><mi>φ</mi><mn>2</mn></msub></mrow></msup></mfenced></mtd></mtr><mtr><mtd><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced><mo>⁢</mo><mfenced open="[" close="]" separators=""><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub></mfenced><mo>+</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub></mfenced><mo>⋅</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><msub><mi>φ</mi><mn>3</mn></msub></mrow></msup></mfenced><mo>⁢</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><mi>φ</mi></mrow></msup><mo>+</mo><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>1</mn></msub></mfenced><mo>⁢</mo><mfenced open="[" close="]" separators=""><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced><mo>+</mo><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced><mo>⁢</mo><msub><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced><mo>⁢</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><msub><mi>φ</mi><mn>3</mn></msub></mrow></msup></mfenced></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></math><img file="EP2201706B1_D0013.tif" /></maths> In the case that the original signal is singly polarized, i.e., it has only one generic polarization component at the transmitter, e.g., <i>E_x</i>, only half of the computation in Eq. (5) is needed.

0055As shown in Eq. (5), five parameters, θ₁, <i>θ</i>₂, δφ, δφ₂, and τ<i>_DGD</i>, are generally needed to recover the original optical signal field, which can be either single polarized or polarization multiplexed. When PMD is sufficiently small, e.g., the PMD induced DGD is much smaller than the signal symbol period, τ<i>_DGD</i> may be safely set to zero in deriving the original signal field, leaving four parameters, θ₁, θ₂, δφ, and δφ₂, to be determined. Since these parameters are generally time varying, it is needed to find the values of these parameters dynamically.

0056The digital signal processing needed to recover the generic polarization components from the reconstructed optical fields of the two orthogonal polarization components of the received optical signal can be performed on a block by block basis, with each block having multiple samples. <figref idref="f0005">FIG. 5</figref> shows an exemplary arrangement to perform the digital signal processing needed to recover the generic polarization components from the reconstructed optical fields of the two orthogonal polarization components of the received optical signal, in accordance with an aspect of the invention. This circuit consists of demultiplexers 501 and 502, M processing units (PUs) 505, and multiplexers 503 and 504.

0057The inputs to the arrangement of <figref idref="f0005">FIG. 5</figref> are <i>E_x'(t)</i> and <i>E_y'(t)</i> and the outputs therefrom are <i>E_x(t)</i> and <i>E_y(t)</i>. Each of demultiplexers 501 and 502 divides the samples it receives over M parallel paths, thereby reducing the processing speed requirement of PUs 505. Eventually multiplexers 503 and 504 multiplex the processed samples to construct <i>E_x(t)</i> and <i>E_y(t)</i>. Note that at any given time, the blocks of samples supplied to one of PUs 505 may have samples overlapping with those of its adjacent PUs. Note also that the multiplexers and demultiplexers can be shared with the field reconstruction process, e.g., as described in Liu-Wei.

0058<figref idref="f0005">FIG. 6</figref> shows an exemplary high level block diagram of arrangement 600 which is suitable to be used to make up each of PUs 505. Arrangement 600, when employed as a PU 505, receives at one time, for each time period, two corresponding blocks of samples, each block being of length N. N is typically greater than 4 and less than 40, with a suitable value being about 10. The value selected for N represents a tradeoff between accuracy achievable and the speed of computation needed to process the samples. The received samples are supplied to feed-forward subprocessor 601 and real-time subprocessor 602. Feed-forward subprocessor 601 finds the best guesses of the parameters needed to recover <i>E_x(t)</i> and <i>E_y(t)</i>, and feeds these parameters, except for δφ₂ to real-time subprocessor 602. Real-time subprocessor 602 receives N pairs of <i>E_x'(t)</i> and <i>E_y'(t)</i> samples, as well as the best guesses of parameters θ₁, θ₂, δφ, and τ_DGD, which were determined by feed-forward path 601, and supplies as an output the N pairs of <i>E_x(t)</i> and <i>E_y(t)</i> samples after processing the received inputs as described hereinbelow.

0059<figref idref="f0006">FIG. 7</figref> shows an exemplary process expressed in flow-chart form for finding the best guesses of the parameters needed to recover E_x(t) and E_y(t), in accordance with an aspect of the invention. This process may be performed in feed-forward subprocessor 601. The process begins in step 702, when N pairs of E_x'(t) and E_y'(t) samples are received. Next, in step 703, using Eq. (5), candidate values of E_x(t) and E_y(t) are calculated for each of the received N pairs of E_x'(t) and E_y'(t) samples. To do so, for each E_x'(t) and E_y'(t) pair, a candidate value is calculated for each possible combination of values for each of the five parameters over their respective physically allowable ranges. For example, the physically allowable ranges may be <i>θ</i>₁ ∈ [0, π), θ₂ ∈ [0, π), δφ ∈ [0,2π), and δφ₂ ∈ [0,2π). For τ_DGD the range employed may be from 0 to the symbol period of the signal, although it is recognized that τ_DGD may actually be larger.

0060Preferably, this is performed by selecting a combination of values for each of the five parameters and using them to compute N candidate values of E_x(t) and E_y(t). In one embodiment of the invention, the guess values for each of the parameters are uniformly distributed over within its allowed range. Typically, 10 to 20 guess values for each parameter should be sufficient.

0061One way to perform the calculation is to implement double loop, where the outer loop is the parameter values and the inner loop is the N sample pairs. A loop so arranged facilitates the computation of step 704, in which the particular values of the five parameters that minimizes a variance-type quantity of the N candidate <i>E_x(t)</i> and <i>E_y(t)</i> is selected. For example, the values of the five parameters that minimizes the variance of a candidate set of <i>E_x(t),</i> i.e., <maths id="math0014" num=""><math display="inline"><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><msup><mfenced open="[" close="]" separators=""><msup><mfenced open="|" close="|" separators=""><msub><mi>E</mi><mi>x</mi></msub><mfenced><mi>t</mi></mfenced></mfenced><mn>2</mn></msup><mo>-</mo><mover><msup><mfenced open="|" close="|" separators=""><msub><mi>E</mi><mi>x</mi></msub><mfenced><mi>t</mi></mfenced></mfenced><mn>2</mn></msup><mo>‾</mo></mover></mfenced><mn>2</mn></msup></math><img file="EP2201706B1_D0014.tif" /></maths> is selected as the best guess. Alternatively, the values of the five parameters that minimizes the variance of a candidate set of <i>E_y(t),</i> i.e., <maths id="math0015" num=""><math display="inline"><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><msup><mfenced open="[" close="]" separators=""><msup><mfenced open="|" close="|" separators=""><msub><mi>E</mi><mi>y</mi></msub><mfenced><mi>t</mi></mfenced></mfenced><mn>2</mn></msup><mo>-</mo><mover><msup><mfenced open="|" close="|" separators=""><msub><mi>E</mi><mi>y</mi></msub><mfenced><mi>t</mi></mfenced></mfenced><mn>2</mn></msup><mo>‾</mo></mover></mfenced><mn>2</mn></msup></math><img file="EP2201706B1_D0015.tif" /></maths> is selected as the best guess. Alternatively, some combination of the two variances may be specified as variance-type quantity to be minimized.

0062The forgoing selection of the generic polarization components of the signal, i.e., <i>E_x(t)</i> and <i>E_y(t,)</i> assumes that as originally transmitted the generic polarization components intrinsically had a constant intensity, i.e., amplitude, which is generally the case for DPSK-type formats, which include at least DBPSK and DQPSK. Alternatively, the best guesses of these parameters may be found using approaches similar to or based on the constant modulus algorithm (CMA).

0063In step 705, the best guesses for the four parameters θ₁, <i>θ</i>₂, δφ, and τ<i>_DGD</i> are supplied as an output, and control passes back to step 702 to process the next N pairs of <i>E_x'(t)</i> and <i>E_y'(t)</i> samples. Note that the values of the four parameters θ₁, θ₂, δφ, and τ_DGD typically tend to change at a rate that is much slower than the signal symbol rate. Thus, feed-forward subprocessor 601 need not process all of the blocks of N pairs of <i>E_x'(t)</i> and <i>E_y'(t)</i> samples that it receives, since doing so will yield essentially the same values for those periods of time over which the parameters remain substantially unchanged. For example, the rate of fiber PMD change is usually slower than 10 KHz, which is 10⁶ times slower than the symbol rate of a 10-Gbaud signal. Advantageously, this significantly relaxes the computation speed required feed-forward subprocessor 601. Of course, should there be a situation in which the rate of change of the values of four parameters θ₁, θ₂, δφ, and τ<i>_DGD</i> is more rapid, they may be computed more often, or even for every block.

0064In another embodiment of the invention, rather than use 10 to 20 guesses for each of the slow varying parameters, θ₁, θ₂, δφ, and τ<i>_DGD</i> only three guess values are employed for each parameter, one being the previous best guess value and the other two being its nearest neighboring guess values. For the angular parameters, the nearest neighbor guess values should be the cyclic neighbors, by which it is generally meant taking modulus with respect to the appropriate value, e.g., 2π for δφ and π for θ₁ and θ₂. The cyclic spacing between the two nearest neighboring guess values should be much smaller than the allowable range of the parameter. Preferably, the spacing is at least 10 times smaller than the allowable range of the parameter. For τ<i>_DGD</i> the nearest neighbors are those values that are one minimum step up and one minimum step down, each step being substantially smaller than the symbol period. Preferably, the spacing is at least 5 times smaller than the allowable range of the parameter. Doing so advantageously reduces the amount of computation that is required.

0065In yet another embodiment of the invention, the guess value for τ<i>_DGD</i> does not need to be searched. Rather, the guess value for τ<i>_DGD</i> can be fixed to a fraction of the symbol period, e.g., 0.4<i>T_S</i>, and useful PMD compensation still results.

0066<figref idref="f0006">FIG. 8</figref> shows an exemplary process, expressed in flow-chart form, which is performed by real-time subprocessor 602 in one embodiment of the invention. In step 802, N pairs of <i>E_x'(t)</i> and <i>E_y'(t)</i> samples, as well as the best guesses of parameters θ₁, θ₂, δφ, and τ<i>_DGD</i> obtained by feed-forward subprocessor 601, are received. Next, step 803 computes <i>E_x(t)</i> and <i>E_y(t)</i> for a set of guess values of δφ₂ ∈ [0, <i>2π</i>) using Eq. (5).

0067Thereafter, in step 804, the best guess of δφ₂ is found. In one embodiment of the invention, the best guess is the guess by which at least one of the variance of the N quantities that represents <i>E_x(t),</i> i.e., <maths id="math0016" num=""><math display="inline"><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><msup><mfenced open="[" close="]" separators=""><msup><mfenced open="|" close="|" separators=""><msub><mi>E</mi><mi>x</mi></msub><mfenced><mi>t</mi></mfenced></mfenced><mn>2</mn></msup><mo>-</mo><mover><msup><mfenced open="|" close="|" separators=""><msub><mi>E</mi><mi>x</mi></msub><mfenced><mi>t</mi></mfenced></mfenced><mn>2</mn></msup><mo>‾</mo></mover></mfenced><mn>2</mn></msup><mo>,</mo></math><img file="EP2201706B1_D0016.tif" /></maths> and that represents <i>E_y(t)</i>, i.e., <maths id="math0017" num=""><math display="inline"><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><msup><mfenced open="[" close="]" separators=""><msup><mfenced open="|" close="|" separators=""><msub><mi>E</mi><mi>y</mi></msub><mfenced><mi>t</mi></mfenced></mfenced><mn>2</mn></msup><mo>-</mo><mover><msup><mfenced open="|" close="|" separators=""><msub><mi>E</mi><mi>y</mi></msub><mfenced><mi>t</mi></mfenced></mfenced><mn>2</mn></msup><mo>‾</mo></mover></mfenced><mn>2</mn></msup><mo>,</mo></math><img file="EP2201706B1_D0017.tif" /></maths> is minimized. Note that typically minimizing one of variance of the N quantities that represents <i>E_x(t)</i> and variance of the N quantities that represents <i>E_y(t)</i> results in the other also being minimized. However, this is not always so, e.g., in the presence of noise, and the implementer may instead choose to minimize the difference between the variance of the N quantities that represents <i>E_x(t)</i> and the variance of the N quantities that represents <i>E_y(t)</i>.

0068In step 805 the N pairs of <i>E_x(t)</i> and <i>E_y(t)</i> samples that correspond to the best guess of δφ₂ are supplied as outputs and control passes back to step 802 to process the next N pairs of E_<i>x</i>'<i>(t)</i> and <i>E_y'(t)</i> samples.

0069<figref idref="f0007">FIG. 9</figref> shows an exemplary high level block diagram arrangement 900 that is suitable to be used to make up each of PUs 505 but which is arranged to speed up the processing as compared to the arrangement shown in <figref idref="f0005">FIG. 6</figref>. The arrangement of <figref idref="f0007">FIG. 9</figref> takes into consideration the fact that parameters θ₁, θ₂, and τ<i>_DGD</i> typically change much more slowly than δφ changes, and therefore they may to be computed at a slower rate. To this end feed-forward subprocessor 601 of <figref idref="f0005">FIG. 6</figref> is further split into first feed-forward subprocessor 901, which computes at a slower rate as compared to feed-forward subprocessor 601 of <figref idref="f0005">FIG. 6</figref>, and second feed-forward subprocessor 903, which computes at the same rate as did feed-forward subprocessor 601 of <figref idref="f0005">FIG. 6</figref>.

0070First feed-forward subprocessor 901 receives as input N pairs of <i>E_x'(t)</i> and <i>E_y'(t)</i> samples and supplies as outputs the best guesses of parameters θ₁, θ₂, and τ<i>_DGD</i>. Second feed-forward subprocessor 903 likewise receives as input N pairs of <i>E_x'(t)</i> and <i>E_y'(t)</i> samples and it also receives as input the best guesses of θ₁, θ₂, and τ<i>C_DGD</i> which are supplied as outputs by first feed-forward subprocessor 901. Second feed-forward subprocessor 903 supplies the best guess value of go as an output to real-time subprocessor 602 and it also passes on the best guesses of θ₁, θ₂, and τ<i>_DGD</i>. Note that, alternatively, the best guesses of θ₁, θ₂, and τ<i>_DGD</i> could be supplied directly from first feed-forward subprocessor 901 to real-time subprocessor 602. Advantageously, since the update rate of first feed-forward subprocessor 901 can be much slower than that of second feed-forward subprocessor 903, the computational speed requirement of the feed-forward path is reduced overall.

0071Real-time subprocessor 602 receives the same N pairs of <i>E_x'(t)</i> and <i>E_y'(t)</i> samples and the best guesses of parameters θ₁, θ₂, δφ, and τ_DGD, as it did in <figref idref="f0005">FIG. 6</figref>, albeit from second feed-forward subprocessor 903 rather than feed-forward subprocessor 601, and outputs the N pairs of <i>E_x(t)</i> and <i>Ey(t)</i> samples based on the signal processing as described hereinabove.

0072The effectiveness of the electronic PMD compensation (PMDC) described hereinabove can be further improved by treating the fiber as if it was made up of multiple segments, each of which has its own "virtual" DGD parameters. More specifically, instead of the three parameters θ₁, δφ, and τ<i>_DGD</i> previously used to describe the fiber-induced PMD, one can treat the fiber link as a concatenation of <i>M</i> PMD segments, i.e., mini fibers, each described by three parameters, θ₁<i>ⁱ</i>, δφ<i>ⁱ</i>, and τ<i>_DGDⁱ</i> where i = <i>1, 2,</i> ... <i>M</i> is the index of the virtual PMD segment. The two received polarization components can be generally linked to the generic polarization components, E_x and <i>E_y</i>, as <maths id="math0018" num="(6)"><math display="block"><mfenced open="[" close="]"><mtable><mtr><mtd><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mi>T</mi><mo>⋅</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msub><mi>E</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>E</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mi>P</mi><mo>⋅</mo><msub><mi>R</mi><mn>2</mn></msub><mo>⋅</mo><mstyle displaystyle="true"><munderover><mo>∏</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover></mstyle><msup><mi mathvariant="italic">PMD</mi><mrow><mi>M</mi><mo>-</mo><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>⋅</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msub><mi>E</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>E</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></math><img file="EP2201706B1_D0018.tif" /></maths> where matrix <i>T</i> represents the polarization transformation of the fiber link, <i>R₂</i> is the rotation matrix associated with the projection of the signal components along the fiber PMD PSP axes at the fiber output on the polarization axes of the PBS used in the receiver, <i>PMDⁱ</i> is the matrix describing the PMD effect of the <i>i</i>-th segment, and <i>P</i> is a phase-delay matrix representing the additional phase delay between the two reconstructed fields after the polarization beam splitting at the receiver.

0073Using the notations shown in <figref idref="f0004">Fig. 4</figref>, the rotation matrix <i>R₂</i> can be written as <maths id="math0019" num="(7)"><math display="block"><msub><mi>R</mi><mn>2</mn></msub><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced></mtd><mtd><mo>-</mo><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced></mtd></mtr><mtr><mtd><mi>sin</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced></mtd><mtd><mi>cos</mi><mfenced><msub><mi>θ</mi><mn>2</mn></msub></mfenced></mtd></mtr></mtable></mfenced><mn>.</mn></math><img file="EP2201706B1_D0019.tif" /></maths> The phase-delay matrix P can be written as <maths id="math0020" num="(8)"><math display="block"><mi>P</mi><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><msub><mi mathvariant="italic">δφ</mi><mn>2</mn></msub></mrow></msup></mtd></mtr></mtable></mfenced><mn>.</mn></math><img file="EP2201706B1_D0020.tif" /></maths> The PMD matrix of the i-th segment, <i>PMDⁱ</i>, can be written in the frequency domain, e.g., after a Fourier transform from the time domain, as <maths id="math0021" num="(9)"><math display="block"><msup><mi mathvariant="italic">PMD</mi><mi>i</mi></msup><mfenced><mi mathvariant="italic">Δf</mi></mfenced><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mfenced separators=""><mn>2</mn><mo>⁢</mo><mi>π</mi><mo>⋅</mo><mi mathvariant="normal">Δ</mi><mo>⁢</mo><mi>f</mi><mo>⋅</mo><msup><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub><mi>i</mi></msup><mo>+</mo><mi>δ</mi><mo>⁢</mo><msup><mi>φ</mi><mi>i</mi></msup></mfenced></mrow></msup></mtd></mtr></mtable></mfenced><mo>⋅</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msup><msub><mi>θ</mi><mn>1</mn></msub><mi>i</mi></msup></mfenced></mtd><mtd><mi>sin</mi><mfenced><msup><msub><mi>θ</mi><mn>1</mn></msub><mi>i</mi></msup></mfenced></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><mfenced><msup><msub><mi>θ</mi><mn>1</mn></msub><mi>i</mi></msup></mfenced></mtd><mtd><mi>cos</mi><mfenced><msup><msub><mi>θ</mi><mn>1</mn></msub><mi>i</mi></msup></mfenced></mtd></mtr></mtable></mfenced><mo>,</mo></math><img file="EP2201706B1_D0021.tif" /></maths> where τ<i>ⁱ_DGD</i> and δφ<i>ⁱ</i> are, respectively, the differential group-delay (DGD) and the phase delay between the two PSP axes of the <i>i</i>-th PMD segment, <maths id="math0022" num=""><math display="inline"><msubsup><mi>θ</mi><mn>1</mn><mi>i</mi></msubsup></math><img file="EP2201706B1_D0022.tif" /></maths> is the angle between the orientation of the signal polarization at the input of the i-th PMD segment and one of its PSP axes, and Δ<i>f</i> is the frequency offset from the center frequency of the signal. In the time domain, the PMD matrix acts on an input signal as follows, <maths id="math0023" num="(10)"><math display="block"><msup><mi mathvariant="italic">PMD</mi><mi>i</mi></msup><mo>⋅</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msubsup><mi>E</mi><mi>x</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><msubsup><mi>E</mi><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mfenced><mi>y</mi></mfenced></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open="[" close="]"><mtable><mtr><mtd><mi>cos</mi><mfenced><msup><msub><mi>θ</mi><mn>1</mn></msub><mi>i</mi></msup></mfenced><mo>⁢</mo><msubsup><mi>E</mi><mi>x</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mfenced><mi>t</mi></mfenced><mo>+</mo><mi>sin</mi><mfenced><msup><msub><mi>θ</mi><mn>1</mn></msub><mi>i</mi></msup></mfenced><mo>⁢</mo><msubsup><mi>E</mi><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><mfenced open="[" close="]" separators=""><mo>-</mo><mi>sin</mi><mfenced><msup><msub><mi>θ</mi><mn>1</mn></msub><mi>i</mi></msup></mfenced><mo>⁢</mo><msubsup><mi>E</mi><mi>x</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msup><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub><mi>i</mi></msup></mfenced><mo>+</mo><mi>cos</mi><mfenced><msup><msub><mi>θ</mi><mn>1</mn></msub><mi>i</mi></msup></mfenced><mo>⁢</mo><msubsup><mi>E</mi><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>⁢</mo><mfenced separators=""><mi>t</mi><mo>-</mo><msup><msub><mi>τ</mi><mi mathvariant="italic">DGD</mi></msub><mi>i</mi></msup></mfenced></mfenced><mo>⁢</mo><msup><mi>e</mi><mrow><mi>j</mi><mo>⋅</mo><mi>δ</mi><mo>⁢</mo><msup><mi>φ</mi><mi>i</mi></msup></mrow></msup></mtd></mtr></mtable></mfenced><mo>,</mo></math><img file="EP2201706B1_D0023.tif" /></maths> where <maths id="math0024" num=""><math display="inline"><msubsup><mi>E</mi><mi>x</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup></math><img file="EP2201706B1_D0024.tif" /></maths> and <maths id="math0025" num=""><math display="inline"><msubsup><mi>E</mi><mi>y</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msubsup></math><img file="EP2201706B1_D0025.tif" /></maths> are the two orthogonal signal polarization components along the two PSP axes of PMD segment <i>i</i>-<i>1</i> at its output when <i>i ≥</i> 2, or the generic signal polarization components, <i>E_x</i> and <i>E_y</i>, when <i>i</i> = 1. In principle, when the polarization transformation matrix T is known, the generic signal polarization components can then be derived from the reconstructed polarization components using <maths id="math0026" num="(11),"><math display="block"><mfenced open="[" close="]"><mtable><mtr><mtd><msub><mi>E</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>E</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><msup><mi>T</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>⋅</mo><mfenced open="[" close="]"><mtable><mtr><mtd><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced></mtd></mtr><mtr><mtd><msub><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub><mfenced><mi>t</mi></mfenced></mtd></mtr></mtable></mfenced><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∏</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover></mstyle><msup><mfenced><msup><mi mathvariant="italic">PMD</mi><mi>i</mi></msup></mfenced><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>⋅</mo><msubsup><mi>R</mi><mn>2</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>⋅</mo><msup><mi>P</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mfenced open="[" close="]"><mtable><mtr><mtd><msub><mi>E</mi><mrow><mi>x</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>E</mi><mrow><mi>y</mi><mo>⁢</mo><mi>ʹ</mi></mrow></msub></mtd></mtr></mtable></mfenced><mn>.</mn></math><img file="EP2201706B1_D0026.tif" /></maths> where -1 indicates the standard matrix inverse operation, i.e., the product of a matrix and its inverse is the identity matrix I.

0074<figref idref="f0007">FIG. 10</figref> shows exemplary high level block diagram of arrangement 1000 suitable to be used to make up each of PUs 505 but which is arranged to treat the fiber as if it was made up of multiple segments, so as to achieve better compensation for PMD, including higher order PMD, than can be achieved than using the arrangements of <figref idref="f0005">FIGs. 6</figref> or <figref idref="f0007">9</figref>. The process for performing such electronic PMDC assumes that the fiber was made up of <i>M</i> "virtual" PMD elements. Feed-forward subprocessor 1001 receives as inputs N pairs of <i>E_x'(t)</i> and <i>E_y'(t)</i> samples and supplies as its output the best guesses of parameters θ₁¹, δφ¹, and τ<i>_DGD</i>¹ ... θ₁<i>^M</i>, δφ<i>^M</i>, and τ<i>_DGD^M.</i> Real-time subprocessor 1002 receives as inputs the N pairs of E_x'(t) and E_y'(t) samples as well as the best guesses of the parameters supplied by feed-forward subprocessor 1001, and supplies as outputs the N pairs of E_x(t) and E_y(t) samples using the methods described hereinabove.

0075As M increases, the PMDC capability also increases. However, the needed computation power to perform the PMDC calculations increases as well. Thus, there is a tradeoff between required computation power and the PMDC that is performed. Note that oftentimes setting M=2 is sufficient to provide better than first order PMD compensation without requiring a severe increase in processing power. For further simplification, the guess values for τ_DGD¹ ... τ_DGD^M can be fixed to a fraction of the symbol period, e.g., 0.4<i>T</i>_s, as noted hereinabove.

0076<figref idref="f0008">FIG. 11</figref> shows high level block diagram of arrangement 1100 suitable to be used to make up each of PUs 505 but which is arranged to treat the fiber as if it was made up of 2 segments each having DGD values each fixed to 0.4<i>T</i>_s.

0077Note that, in accordance with an aspect of the invention, the digital PMD compensation schemes described above may also be employed with so-called "dual- polarization coherent-detection" receivers, where digital representations of two orthogonal polarization components of the received optical signal are obtained. A typical dual-polarization coherent-detection receiver with DSP is shown in "Uncompensated Transmission of 86Gbit/s Polarization Multiplexed RZ-QPSK over 100km of NDSF Employing Coherent Equalisation" by Fludger et al., which was published as ECOC'06 post-deadline paper Th4.3.3. <figref idref="f0008">FIG. 12</figref> shows an exemplary high level block diagram of arrangement 1200 which is suitable to be used to make up each of PUs 505 but which is arranged for use with a coherent-detection receiver, in accordance with an aspect of the invention.

0078As is well known, such coherent detection receivers employ, an optical local oscillator (OLO) that provides an absolute phase reference for both received polarization components. As a result, there is no uncertainty in the additional phase difference δφ₂ due to the field reconstruction process. Consequently, δφ₂ does not need to be estimated in real-time subprocessor 1202, and instead can be estimated in feed-forward subprocessor 1201 so that the computational effort in the real-time subprocessor is much reduced. Typically, δφ₂ changes at a very low speed, e.g., &lt;1 KHz, so feed-forward subprocessor 1201 only needs to update δφ₂ at a much lower speed than the signal date rate.

0079<figref idref="f0009">FIG. 13</figref>, similar to <figref idref="f0006">FIG. 8</figref>, shows an exemplary process, expressed in flowchart form, that is performed in real-time subprocessor 1202, in one embodiment of the invention. In step 1302 N pairs of E_x'(t) and E_y'(t) samples, as well as the best guesses of parameters θ₁, <i>θ</i>₂, δφ, δφ₂, and τ<i>_DGD</i> for each of M virtual segments obtained by feed-forward subprocessor 1201, are received. Next, step 1303 computes <i>E_x(t)</i> and <i>E_y(t)</i> using Eq. (11). In step 1305 the N pairs of <i>E_x(t)</i> and <i>E_y(t)</i> samples are supplied as outputs and control passes back to step 1302 to process the next N pairs of <i>E_x'(t)</i> and <i>E_y'(t)</i> samples.

0080As mentioned hereinabove, to simplify, the DGD values can be fixed, e.g., each equal to about 0.4Ts. Also, M can be chosen to be 2.

0081Once the original signal field is obtained, further compensation for other impairments, e.g., chromatic dispersion and/or self-phase modulation, and data recovery following suitable demodulation, can be performed to extract the data content from at least on generic polarization, using compensation for other impairments unit 152, demodulation unit 154, and data recovery unit 155.

0082A practical issue with the use of differential detection is that it employs ODIs, which typically exhibit polarization-dependent phase shift (PDPS). In other words, generally, the phase offset between the two arms of an ODI, φ, is dependent on the polarization state of the optical signal. When the signal polarization is aligned with one of the two characteristic polarization orientations of the ODI, the phase offset φ reaches its maximum or its minimum. The PDPS is the difference between the maximum and the minimum phase offsets.

0083There are three common types of ODI: 1) fiber-based, 2) planar lightwave circuit (PLC)-based, and 3) free-space optics-based. The PDPS of a fiber-based ODI is typically due to the birefringence of the fiber resulting from mechanical stress. The PDPS of a PLC-based ODI is typically due to the birefringence of the waveguide structure of the PLC. The PDPS of a free-space optics based ODI is typically due to the polarization-dependent phase delay of the beam splitter used in forming the two optical interference paths. The PDPS can range from about 2 degrees (0.035 rad.) to about 20 degrees (0.35 rad.), the particular value for any ODI depending on its design.

0084The field reconstruction process as described herein relies on the phase difference estimation at multiple sampling locations. If the signal polarization is not aligned with one of the two characteristic polarization orientations of the ODI, the PDPS will accumulate as the number of sampling points increases and prevent accurate phase estimation. Thus, it is preferred to align PBS 101 in such a way that the polarization of each of the two split signals, i.e., <i>E_x'</i> and <i>E_y'</i> is aligned with one of the two characteristic polarization orientations of the ODIs. This can be achieved by using, e.g., 1) polarization-maintaining fibers with suitable orientations to connect the two outputs of PBS 101 with the two inputs of ODI pairs 105-106 and 107-108, or 2) polarization-maintaining free space optical connections between PBS 101 and the two inputs of ODI pairs 105-106 and 107-108 when the ODIs are free-space optics based.

0085<figref idref="f0002">FIG. 2</figref> shows an embodiment of the invention similar to that shown in <figref idref="f0001">FIG. 1</figref> but in which the intensity detection branches are omitted. In accordance with an aspect of the invention, the intensity profile for each polarization component is approximated from the absolute value of its respective one of the complex waveforms rather than directly recorvered from the received optical signal.

0086<figref idref="f0003">FIG. 3</figref> shows another embodiment of the invention similar to that shown in <figref idref="f0001">FIG. 1</figref> but which employs only a single I/Q ODI pair. This is achieved by eliminating PBS 101 of <figref idref="f0001">FIG. 1</figref> and employing in lieu thereof four PBSs 301 at the four outputs of the first I/Q ODI. Furthermore, as in <figref idref="f0002">FIG. 2</figref>, the intensity detection branches are omitted and the intensity profile for each polarization component is approximated rather than directly recovered from the received optical signal. Advantageously, the cost of the overall arrangement, as with respect to the arrangement of <figref idref="f0001">FIG. 1</figref>, is significantly reduced.

0087In the embodiment of <figref idref="f0003">FIG. 3</figref> the received signal is supplied directly into a single polarization-independent I/Q ODI pair made up of ODIs 305 and 306, whose four outputs are each connected to a respective one of PBSs 301 with the same polarization orientation. The eight outputs from PBSs 301, consisting of four <i>x</i>'-polarized outputs and four <i>y</i>'-polarized outputs, are then treated in the same manner as the outputs of the I/Q ODI pairs of <figref idref="f0001">FIG. 1</figref>, i.e., detected by balanced detectors whose outputs are sampled, after optional amplification and gain control, by a respective ADC. The sampled waveforms are then processed to obtain a digital representation of the signal optical field as described hereinabove.

0088As will be readily understood by those of ordinary skill in the art, the instant invention may be applied to optical differential phase-shift keying (DPSK) signals, such as differential binary phase-shift keying (DBPSK) and differential quadrature phase-shift keying (DQPSK) signals, since ODI(s) and balanced detection are commonly used for DPSK detection. Furthermore, this invention may also be applied to amplitude-shift keying (ASK), combined DPSK/ASK, and differential QAM.