Reconstruction and restoration of an optical signal field

RECONSTRUCTION AND RESTORATION OF AN OPTICAL SIGNAL FIELD

Technical Field

This invention relates to the reconstruction and restoration of an optical signal field.

Background of the Invention

Linear and nonlinear effects distort optical signals transmitted over optical fibers. Such effects include chromatic dispersion (CD) and self-phase modulation (SPM). Optical dispersion compensation is typically employed to reduce signal distortion that arises as a result of CD.

Electronic 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 M. S. O'Sullivan, K. Roberts, and C. Bontu, in *Electronic dispersion compensation techniques for optical communication systems,” ECOC'05, paper T113.2.1,2005, EDC can be performed at the transmitter. Doing so is referred to herein as pre-EDC. Alternatively, as described by S. Tsukamoto, K. Katoh, and K. Kikuchi, in “Unrepeated Transmission of 20-Gb/s Optical Quadrature Phase-ShiftKeying 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, EDC can be performed at the receiver, which is referred to herein as post-EDC.

Post-EDC has an advantage over pre-EDC, in that post-EDC does not require that performance feedback be supplied from the receiver to the transmitter. Unfortunately, direct intensity detection, also known as square-law detection, which is the commonly used optical detection technique of today’s optical fiber communications systems, e.g., the optical to electronic conversion performed by photodiodes, only recovers the optical signal amplitude and cannot recover the optical signal phase information, thus making the performance of post-EDC much poorer than that of pre-EDC.

To overcome this shortcoming, and hence enhance the performance of the post-EDC, the Tsukamoto et al. article suggests employing coherent detection to fully reconstruct the optical signal's complex field, i.e., both amplitude and phase. However, disadvantageous^, as compared to direct intensity detection, coherent detection is much more sophisticated, and hence more expensive and difficult to perform. Further disadvantageous^, coherent detection requires the use of an optical local oscillator (OLO), as well as phase and polarization tracking between the OLO and the signal carrier.

Summary cf the Invention

In accordance with the principles of the invention, a digital version of the 5 complex optical field, i.e., both amplitude and phase, e.g., with respect to a reference point, of a received optical is developed at a receiver by employing direct differential detection in conjunction with digital signal processing.

More specifically, as is well known, the complex optical field of any signal can be reconstructed by knowing its intensity and phase profiles. ^,The intensity profile 10 may be obtained by conventional direct intensity detection. As to obtaining the phase, in accordance with an aspect of the invention, first an electronic analog representation of a complex waveform that contains information about the phase differences between adjacent locations that are separated by a prescribed time difference ΔΤ in the received signal is obtained by employing a pair of optical delay interferometers that have 15 orthogonal phase offsets, i.e., the difference between the phase offsets is π/2, followed by two balanced intensity detectors. The output of the first interferometer after the balanced intensity detection is the real part of the complex waveform, while the output of the second interferometer after the balanced intensity detection is the imaginary part of the complex waveform. The output of each of the balanced intensity detectors, 20 and the intensity profile if obtained by direct intensity detection, are. converted to a digital representation using analog to digital conversion. The sample period for the analog to digital conversion may be shorter than ΔΤ, so that multiple samples may exist within a period of ΔΤ. From the digital representation of the complex waveform, the phase difference between the adjacent locations that are separated by ΔΤ may be 25 obtained. Then, based on the obtained phase differences, and optionally, a search for an initial phase offset among the multiple samples within a period of ΔΤ, the phase relationship among all the samples is obtained. Essentially the absolute phase profile for the received signal is thus derived with the only uncertainty being that of a constant phase shift, which is insignificant.

To simplify the hardware necessary, optionally, the intensity profile may be approximated from the absolute value of the complex waveform rather than obtaining it by direct intensity detection. Furthermore, optionally, once the intensity profile and the phase profile of the optical signal as received are recovered, digital signal processing may be employed to compensate for the distortions in the received signal, 35 e.g., signal distortions due to chromatic dispersion and SPM, so that an accurate representation of the originally transmitted optical signal waveform may be reconstructed electronically.

PCI7US2007/020367.

The 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).

Brief Description of the Drawing

In the drawing:

FIG. 1 shows an exemplary apparatus for reconstructing and restoring an optical signal field in accordance with the principles of the invention; and

FIG. 2 shows an embodiment of the invention similar to that shown in FIG. 1 but in which the intensity profile is approximated rather than directly recovered from the received optical signal.

Detailed Description

The following merely illustrates the principles of the invention. It will thus be appreciated that those skilled in the art will be able to devise various arrangements that, although not explicitly described or shown herein, embody the principles of the invention and are included within its spirit and scope. Furthermore, all examples and conditional language recited herein are principally intended expressly to be only for pedagogical purposes to aid the reader in understanding the principles of the invention and the concepts , contributed by the inventor(s) to furthering the art, and are to be construed as being without limitation to such specifically recited examples and conditions. Moreover, all statements herein reciting principles, aspects, and embodiments of the invention, as well as specific examples thereof, are intended to encompass both structural and functional equivalents thereof. Additionally, it is intended that such equivalents include both currently known equivalents as well as equivalents developed in the future, i.e., any elements developed that perform the same function, regardless of structure.

Thus, 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.

The 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.

In 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.

Software 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.

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

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

FIG. 1 shows an exemplary apparatus, typically in a receiver, arranged in accordance with the principles of the invention, for developing the entire complex optical field of a received optical signal by employing direct differential detection in conjunction with digital signal processing and for compensating for various impairments that were inflicted upon the optical signal as it traveled from its source. FIG. 1 shows a) 1x3 optical splitter 1001; b) optical delay interferometers (GDIs)

1002 and 1003; c) balanced intensity detectors 1011 and 1013; d) photodiode 1015; e) amplifiers 1021, 1022, and 1023; !)optional automatic-gain controllers (AGCs) 1031, 1032, and 1033; g) analog-to-digital converters (ADCs) 1041, 1042, and 1043; and h) digital signal processing unit 1050.

More specifically, 1x3 optical splitter 1001 replicates the incoming optical signal so as to produce three copies. The optical power allotted to each of the copies from the originally input optical signal is,at the discretion of the implemented 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 GDIs 1002 and 1003 and the remaining power, e.g., between 10 and 20 percent, is supplied to photodiode 1015.

As will be readily recognized by those of ordinary skill in the art, optical delay interferometers (ODIs) 1002 and 1003 may be any type of interferometer having the required characteristics. For example, ODIs 1002 and 1003 may be based on the well-known, so-called Mach־Zehnder interferometer. Alternatively, ODIs 1002 and

1003 may be based on the well-known, so-called Michaelson interferometer.

OD1 1002 has a delay of about ΔΤ in the optical path between its respective two arms and a phase difference, i.e., offset, of $>, where

T ״ zn

Δ?־־-*-—_t 1 <m<, sps, and mis an integer, (1) sps and where T$ is the symbol period of the signal, sps is the number of samples per symbol taken by analog to digital converters 1041, 1042, and 1043, m is an integer between 1 and sps, and is an arbitrarily selected number. If so, the free spectral range (FSR), i.e., 1/ΔΤ, of ODIs .1002 and 1003 is related to the signal symbol rate (SR) as fsr _a SR *. Note that, based on numeric simulations, it has been found m that, preferably, sps be set to a value of 4, and m can be a value of 1, 2, 3, or 4. This is because an sps value of less than 4 tends to not be sufficient to accurately represent the signai waveform sufficiently given the procedures described hereinbelow, while sps greater than 4 provides only negligible improvement.

The 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 AT*C/w, where C is the speed of light in vacuum and n 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 not precisely the phase change is effectively

The total length change used to achieve the effective length change of fa may be some percentage of the length ΔΤ-C/n. 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 5 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 ODI 1002.

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

ODI 1003 is similar to ODI 1002, in that it has delay of about ΔΤ in the optical path between its respective two arms, but between its arms it has a phase offset of 15 fa-11/2. Thus, the difference between the phase offsets of ODIs 1002 and 1003 is π/2, so ODI 1002 and 1003 are said to have orthogonal phase offsets.

Balanced intensity detectors 1011 and 1013 are conventional. Typically, each of balanced intensity detectors 1011 and 1013 is made up of a pair of well-matched ־ photodiodes. Balanced intensity detectors 1011 and 1013 convert the output of each 20 of the arms of GDIs 1002 and 1003 to an electrical representation. Thus, balanced intensity detectors 1011 and 1013 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 ΔΤ in the received optical signal.

Photodiode 1015 performs conventional direct intensity detection, and thus obtains the intensity profile of the received optical signal in electronic form.

Amplifiers 1021, 1022, and 1023 amplify the signals supplied as outputs by balanced intensity detector 1011, balanced intensity detector 1013, and photodiode 1015, respectively. Typically, amplifiers 1021, 1022, and 1023 convert the current 30 . which is output by the various photodiodes of balanced intensity detector 1011, balanced intensity detector 1013, and photodiode 1015 to respective corresponding voltages. To this end, amplifiers 1021, 1022, and 1023 may be trans-impedance amplifiers. Furthermore, amplifiers 1021 and 1022 may be differential amplifiers. After amplification, each of the outputs is typically single ended. Optional automatic35 gain controllers (AGCs) 1031, 1032, and 1033 may be employed to normalize the electronic waveforms prior to digitization.

Analog-to-digital converters (ADCs) 1041, 1042, and 1043 perform “digital sampling” of the amplified signals to develop a digital representation of the amplified signals. ADCs 1041,1042, and 1043 typically have the same resolution, e.g'., 8 bits.

Digital signal processing unit 1050 receives the digital representation of the 5 amplified signals and develops a digital representation of the amplitude and phase profiles of the received optical signal, in accordance with an aspect of the invention. In particular, reconstruction unit 1051 performs such development. Furthermore, in accordance with another aspect of the invention, digital signal processing unit 1050 may develop a digital representation of the original waveform of the optical signal as 10 it was transmitted prior to undergoing impairments in the channel over which it passed by digitally compensating for various ones of the transmission impairments experienced by the optical signal, e.g., chromatic dispersion aiid/or self-phase modulation. Restoration unit 1052 performs such restoration. Lastly, demodulation and data recovery unit 1053 performs demodulation and conversion to actual bits.

An exemplary process for recovering the entire complex optical signal field by coupling direct differential detection with digital signal processing in accordance with the principles of the invention, using the arrangement of FIG. 1, is as follows. First, the intensity profile of the received optical field is obtained by direct intensity detection using photodiode 1015. The intensity profile, represented by I(t), is 20 computed as

W -y(M1)* (2) where y(t) is the received complex optical field as it arrives at coupler 1001 and * denotes complex conjugate.

The outputs of balanced detectors 1011 and 1013 are analog representations 25 of, respectively, the real, u_rc0!(t), and imaginary, u_lmag(t), parts of the following complex waveform that contains information about the phase differences between two time locations separated by AT = y(/)y((-Ar)*exp[;^₀] (3) = |y(i)y(f - AT)|expUW(/) -- φ{ί - AT)+0_O]} using the following definitions:

7(0 =17(01 exp[/^(i)], ³⁰ X/-AT)־|X/-AT)|expW(/-AT)]. ^U

After the analog representations of the real and imaginary parts of the complex waveform u(t) are amplified, they are converted into digital representations by sampling, e.g, by ADCs 1041 and 1042. Likewise, after the intensity profile is amplified, it too is converted to a digital representation by sampling, e.g, by ADC 35 1043. ADCs 1041 and 1042 may be consider an ADC unit, which may also include

PCI7US2007/020367

ADC 1043. Sampling of the complex waveform and the intensity profile is performed at the following time locations (t₅):

'! י...1) .¾+,/st-bit) ,,+7^,/,+7^+^7^,/1+7^+^7,,..,/,+75+^7^, (2nd-bit) /,+λΓ₅, t1+nT_s + -^T_s, / + nT_s +-^T_S,.., t_f+nT_s+^-T_s, (nth-bit) where t, is an initial, arbitrary, time position and n is an arbitrarily selected number for use in showing bow the equation is generalized to any bit position.

For example, for sps=4, the sampling time locations are as follows: ,5^4+ן, ין» ,ι + i^s» ,!+4־^i> (Ist-bit) ,,+7^, /, + 7^ + ^-75, ,,+7^+£¾. /, + 7_S +-J7_S, (2nd-bit) ,,+75״. /, + nT_s +i7_s, /!+7₅|+75מ, /,+n75+|75, (nth-bit)

After the digital representations of the real and imaginary parts of the complex 10 waveform, and Uimagfts). are obtained, they are supplied to digital signal processing unit 1050. Likewise, after the digital representation of the intensity waveform, 7(/^, is obtained, it too is supplied to digital signal processing unit 1050.

The digital- samples are first used to reconstruct the amplitude and phase profiles of the received optical signal by reconstruction unit 1051. This reconstruction IS step may include the following procedures.

First, a group of samples from each sampled waveform Unaifts)» and are selected as a ״frame” to be processed together. The size of the frame, i.e., the number of symbols for which samples are taken, is chosen to be larger than the maximum number of optical symbols that interact during optical transmission as 20 result of chromatic dispersion or other effects during optical transmission. Note that by interacting it is meant that the pulses that make up the symbols overlap each other due to the broadening of the pulses caused by the dispersion property of the fiber. For example, for a 20-Gb/s DQPSK signal experiencing a chromatic dispersion of 17,000 ps/nm, which corresponds to the same dispersion that would be produced by 1,000 km 25 standard single-mode fiber (SSMF), the maximum number of interacting optical . symbols is about 30. For such an exemplary situation, a suitable frame size may be 64' symbols, or 64 sps samples.

Secondly, the filtering effect due to the bandwidth limitations of photodetectors 1011, 1013, and 1015 and ADCs 1041 1042, and 1043 may need to be 30 compensated for by inversely filtering the digital waveforms. In other words, the inverse of the filter transfer function caused by the superposition of the photo-detector response and the ADC response is digitally applied to the digital waveform.

Thirdly, the phase factor that represents the optical phase differences between samples that are separated by ΔΤ in time, Δ^(/,) = ^(/,)-^(/,-ΔΓ), as given by 5 equation 3, may be obtained as follows expUA^,)]־ expOWJ-^ζ -Δ7⁵) . = {[(־)

Notice that although one still needs to know the value of ¢0, the computation of equation 5 effectively removes the impact of ¢0 in obtaining the phase factor, so that ¢) can be of any arbitrary value. The finding of the value of ¢0 could be achieved 10 by a real. world search, e.g., an automated search, which varies the value of ¢0 until an optimum guess is found. The guess that yields the least bit error rate is selected as the optimum guess. Alternatively, the guess that provides the best optical signal spectrum of the signal as reconstructed as described hereinbelow can be selected as the optimum guess. Another possibility is that instead of performing a search, all the 15 results using different values of ¢0 over the range from 0 to 2n can be computed and the value of ¢0 that gives the best result is selected as the optimum guess. Doing so allows moving directly, and hence possibly more quickly, to the value of ¢0. For example, the computation could be performed for 40 possible candidate values of ¢0 with a spacing between each candidate value of 0.05 π.

Fourthly, in theory, the signal phase profile of each respective “subgroup” of samples within a frame, each subgroup consisting of those samples of the frame that have a spacing between them of ΔΤ or a whole number multiple thereof, may be obtained, based on the optical phase differences of adjacent samples within the subgroup, by determining ¢(/, + מ ΔΓ) = φ(ΐ ן )++ p ΔΤ), (subgroup 1) ¢(/, + n ΔΤ+—T, ) = +—T_s ) + Χ Δ^(Ζ, +—T_f 4 p .AT), (subgroup 2) sps sps sps ς t

^(/,+η-ΔΤ +^—^T_t) = ^(t, +—-Τ₅)+^ Δ<4(/, -T, + p-AT) (subgroup״׳) sps sps Τ., sps

6) where n is the position of a particular sample within the subgroup and when n=0 the summation is not computed at all.

Practically, rather than obtaining the phase directly, it is sufficient to obtain 30 merely the phase factors, as follows.

_β/Λ׳_ι+π.ΔΠ (subgroup!) • ^P¹ ,Λ׳ι+η.Δ7+—7,) Λ״ (,7—+,׳ /Δ«׳,+αΔΓ+-!-Γ,) e ^J״׳ =e ״ .ΓΤε *״ , (subgroup 2)

17) 1־) ^,_1+״._Δτ<-—7,) «»,-U^lr,) ״ /MCf.+p-AT+^lr,) _e '״׳ =_e pje V. (subgroup״!) p־l where n is the position of a particular sample within the subgroup and when n=0 the multiplication is not computed at all. The phase factors give the phase correlation among the samples within each subgroup. However, the phase relationship among the subgroups is not yet known. Thus, it is necessary to determine m-1 phase differences. Once the phase relationship among the like-spaced samples of the subgroups, e.g., the first samples, i.e., the samples for which n=0, is known, then the phase relationship among all of the samples will be completely specified. For example, for n=0, the differences between each adjacent pairing of the terms before the summation symbol in equation 6, for example *(^,״)*>-( ן!), should be determined.

The phase relationship among the like-spaced samples of these subgroups may be estimated as follows. An initial phase difference which can be any value between 0 and 2π is selected as a candidate phase difference for candidate pairing of any two of these samples out of all possible phase offsets to obtain a “trial phase relationship” among all the samples in the frame. The initial phase difference may be between 0 and 2n, as this is the range of the actual phase difference. It has been found that a good initial candidate phase difference is 0.1 π. Furthermore, since various candidate phase offsets will be tried in order to determine the best one using a searching process, it is necessary to select a resolution for which the candidate phase offsets will be selected. A good value for the resolution has been found to be 0.1 π. Thereafter, the optical signal field is reconstructed to produce a trial reconstructed optical signal based on the selected phase difference and the known intensity profile I(t_s). This may be achieved by determining (8) where £,(/,) is the reconstructed signal for the current set of values that is the current estimate of the received optical signal.

The optical power spectrum of the trial reconstructed signal is then obtained by performing a Fourier transformation on the trial reconstructed signal. The power for that portion of the trial reconstructed signal that falls within the frequency range of [-SR, +SR] about the signal center frequency is obtained. This process is repeated by selecting a new candidate phase offset, e.g., by increasing the previous candidate phase. The set of “trial” phase offsets among the like-spaced samples of these subgroups that gives the maximal spectral power within [-SR, +SR] about the signal center frequency is selected as the best estimate. The phase relationship among all samples in the frame can then be determined based on the best estimate.

Alternatively, the set of trial phase offsets among the like-spaced samples of these subgroups that gives the minimal spectral power outside [-SR, +SR] about the signal center frequency is selected as the best estimate. The phase relationship among all samples in the frame can then be determined based on the best estimate.

In one embodiment of the invention, it may be desirable to set ΔΤ =T_t/ sps. As a result, the delay ΔΤ is equal to the sampling resolution, m=l, and hence there is only one subgroup in a frame, so all of the samples have a phase relationship with their immediately adjacent samples. In such an embodiment of the invention, the phases of all the samples in a frame can be obtained straightforwardly, in theory, by determining ¢(/, = ij + η ΔΤ) = ¢(/,) + Δ&lt;/1 + ΔΤ) + Δ^(ί1 + 2 . ΔΓ).. + Δφ(ί,), (9) which is a special case of equation 6, i.e., only the first subgroup, which is the only subgroup, is computed.

Practically, rather than obtaining the phases directly, it is sufficient to obtain merely the phase factors for each of the samples, as follows:

n

P=* (10) which is a special case of equation 7, i.e., only the first subgroup, which is the only subgroup, is computed.

Finally, the digital representation of the received optical signal field, £_ff(/,), can be obtained based on the obtained phase factor and the intensity profile I(t₅) by £»(',) = /7(0 (11)

In one embodiment of the invention, shown in FIG. 2, when ΔΤ is sufficiently small as compared to the symbol period T_s, the intensity profile may be approximated by |u(t₅)|, and so «»(־&lt;׳ i/Ro[׳״׳־’.

(12) or preferably i 1V4

Note that ΔΓ may be considered to be sufficiently small when it is a least a factor of 2 smaller than the symbol period, i.e., AT s T_s /2 . For Sps=4, preferably, AT = T_s 14 .

Using this approximation means that photodiode 1015, amplifier 1023, optional automatic-gain controller 1033, and analog-to-digital converter 1043, are not 5 required, and hence they are not shown in FIG. 2. Also, 1x3 optical splitter 1001 is replaced by simpler 1x2 optical splitter 2001, since there is no need for a branch to determine the intensity, and hence only two copies are required.

After the received optical signal field is developed in the digital domain by reconstruction unit 1051, the digital representation of the optical signal field as 10 originally launched from a transmitter, £7(1 J, can then be derived by restoration unit 1052. To this end, restoration unit 1052 electronically compensates for various distortions, such as distortion caused by chromatic dispersion and, in accordance with an aspect of the invention, a) self-phase modulation (SPM) and b) combinations of chromatic dispersion and SPM that the transmitted signal suffered as it traveled to the 15 receiver.

When the signal is primarily distorted by chromatic dispersion, restoration unit 1052 may restore die original optical signal field by determining £r(G) » (13) where F(x) and are, respectively, the Fourier and inverse-Fourier 20 transformations of signals x and y, /(O_Wo/) represents the frequency-dependent modification of the optical phase of the signal due to the dispersive effect resulting from a dispersion with value D, and the sign indicates the removal of the dispersive effect. More simply, this may be approximated using conventional techniques employing finite impulse response (FIR) filters.

When the signal is distorted essentially only by SPM, such SPM may be compensated for by an embodiment of the invention in which restoration unit 1052 determines (14) where F(x) and are, respectively, the Fourier and inverse-Fourier transformations of signals x and y, as before, represents the total nonlinear phase due to the SPM, and the minus sign indicates the removal of the dispersive effect.

When, the signal is distorted by both chromatic dispersion and SPM such combined chromatic dispersion and SPM may be compensated for by an embodiment 35 of the invention in which restoration unit 1052 treats the fiber link connecting the transmitter and the receiver as being made up of N segments, each having the same dispersion and SPM effects, where the segment that is closest to the transmitter is considered to be the first segment and the segment that is closest to the receiver is considered to be the N⁰¹ segment. Restoration unit 1052 then obtains the digital representation of the original optical field by performing the iterative process embodied by the following pseudocode:

E(t״N + l) = EM, for n = N to 1

E_D(t_t) = ?{F[£(r״n +1)].

ER, n) = EM . , end ^R)־£R.l) where ER,jV)is the restored optical field at the beginning of the N-th segment, represents the total nonlinear phase due to the SPM.

After the original optical field is restored in the digital domain, it is further processed by demodulation and data recovery unit 1053. For example, when the optical signal is modulated by DQPSK format, the conventional optical DQPSK demodulation process obtains the decision variables for the in-phase (I) and quadrature (Q) data tributaries by determining

M,)*׳w/ £_rR) £_rR-r,)*.exp(/^) , w_e R ) = imag E_T (t,) E_T R - T_t) * ־ exp(y ך) .

(15)

Once the decision variables are obtained, a decision can be made to recover the original data I- and Q- tributaries transmitted at the transmitter through ״),.» ׳>! ״/»»ί».

10.

_r x — (' ״oto)²'׳* ^cQV&lt;tJ - (ο. ״_β&lt;ω&lt;ι׳_ή where l_a is the decision time and V_lh is the decision threshold, which is usually about (16) zero.

As will be readily understood by those of ordinary skill in the art, optional receiver performance monitoring can be used to provide information on how well the reconstruction and restoration processes are succeeding in recovering the original optical signal. Moreover, a feedback control may be applied to optimize each step in the reconstruction and restoration processes. For example, in the case where is slowly changing with time, e.g., due to the drift of the frequency of the optical signal carrier at the transmitter or temperature-induced path length changes in the GDIs, equation (5) can be dynamically adjusted, with a feedback control, to always find a best guess for the time-varying fa so as to thereby accurately obtain the phase factor.

As 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 5 commonly used for DPSK. detection. Furthermore, this invention may also be applied to amplitude-shift keying (ASK), combined DPSK/ASK, and differential QAM.