Efficient backward-propagation using wavelet-based filtering for fiber backward-propagation

With the goal of reducing the number of operations required for digital backward-propagation used for fiber impairment compensation, wavelet-based filtering is presented. The wavelet-based design relies on signal decomposition using time-limited basis functions and hence is more compatible with the dispersion operator, which is also time-limited. This is in comparison with inverse-Fourier filter design which by definition is not time-limited due to the use of harmonic basis functions for signal decomposition. Artificial, after-the-fact windowing may be employed in this case; however only a limited amount of saving in the number of operations can be achieved, compared to the wavelets-base filter design. Waveletbased filter design procedure and numerical simulations which validate this approach are presented in this paper. 2009 Optical Society of America OCIS codes: (060.1660) Coherent communications; (060.4370) Nonlinear optics, fibers. References and links 1. A. D. Ellis and F. C. G. Gunning, "Spectral density enhancement using coherent WDM," IEEE Photon. Technol. Lett. 17, 504-506 (2005). 2. G. Goldfarb, G. Li, and M. G. Taylor, "Orthogonal Wavelength-Division Multiplexing Using Coherent Detection," IEEE Photon. Technol. Lett. 19, 2015-2017 (2007). 3. S. J. Savory, G. Gavioli, R. I. Killey, and P. Bayvel, "Electronic compensation of chromatic dispersion using a digital coherent receiver," Opt. Express 15, 2120-2126 (2007). 4. G. P. Agrawal, Nonlinear Fiber Optics, 4th. ed., (Academic Press, 2007). 5. A. Yariv, D. Fekete, and D. M. Pepper, "Compensation for channel dispersion by nonlinear optical phase conjugation," Opt. Lett. 4, 52–54 (1979). 6. R. A. Fisher, B. R. Suydam, and D. Yevick, "Optical phase conjugation for time-domain undoing of dispersive self-phase-modulation effects," Opt. Lett. 12, 611-613 (1983). 7. C. Pare, A. Villeneuve, P. A. Belanger, and N. J. Doran, "Compensating for dispersion and the nonlinear Kerr effect without phase conjugation," Opt. Lett. 21, 459-461 (1996). 8. R.-J. Essiambre, P. J. Winzer, X. Q. Wang, W. Lee, C. A. White, and E. C. Burrows, "Electronic predistortion and fiber nonlinearity," IEEE Photon. Technol. Lett. 18, 1804-1806 (2006). 9. R. I. Killey, P. M. Watts, V. Mikhailov, M. Glick, and P. Bayvel, "Electronic dispersion compensation by signal predistortion using digital Processing and a dual-drive Mach-Zehnder Modulator," IEEE Photon. Technol. Lett. 17, 714-716 (2005). 10. K. Roberts, L. Chuandong, L. Strawczynski, M. O'Sullivan, and I. Hardcastle, "Electronic precompensation of optical nonlinearity," IEEE Photon. Technol. Lett. 18, 403-405 (2006). 11. E. Yamazaki, F. Inuzuka, K. Yonenaga, A. Takada, and M. Koga, "Compensation of Interchannel Crosstalk Induced by Optical Fiber Nonlinearity in Carrier Phase-Locked WDM System," IEEE Photon. Technol. Lett. 19, 9-11 (2007). 12. L. B. Y. Du and A. J. Lowery, "Fiber nonlinearity precompensation for long-haul links using direct-detection optical OFDM," Opt. Express 16, 6209-6215 (2008). 13. X. Li, X. Chen, G. Goldfarb, E. Mateo, I. Kim, F. Yaman, and G. Li, "Electronic post-compensation of WDM transmission impairments using coherent detection and digital signal processing," Opt. Express 16, 880-888 (2008). 14. G. Goldfarb, M. G. Taylor, and G. Li, "Experimental Demonstration of Fiber Impairment Compensation Using the Split-Step Finite-Impulse-Response Filtering Method," IEEE Photon. Technol. Lett. 20, 18871889 (2008). 15. G. Goldfarb and G. Li, "Demonstration of fibre impairment compensation using split-step infinite-impulseresponse filtering method," Electron. Lett. 44, 814-816 (2008). 16. X. Li, C. Xingzhong, and M. Qasmi, "A broad-band digital filtering approach for time-domain Simulation of pulse propagation in optical fiber," J. Lightwave Technol. 23, 864-875 (2005). #107052 $15.00 USD Received 2 Feb 2009; revised 4 Apr 2009; accepted 6 Apr 2009; published 12 May 2009 (C) 2009 OSA 25 May 2009 / Vol. 17, No. 11 / OPTICS EXPRESS 8815 17. T. Kremp and W. Freude, "Fast split-step wavelet collocation method for WDM system parameter optimization," J. Lightwave Technol. 23, 1491-1502 (2005). 18. G. Beylkin, "On the Representation of Operators in Bases of Compactly Supported Wavelets," SIAM J. on Numerical Analysis 29, 1716-1740 (1992).

This Article is brought to you for free and open access by the Faculty Bibliography at STARS. It has been accepted for inclusion in Faculty Bibliography 2000s by an authorized administrator of STARS. For more information, please contact lee.dotson@ucf.edu.

Introduction
One of the benefits associated with coherent optical transmission systems is the possibility to obtain higher spectral-efficiency compared to direct or differential detection. High spectralefficiency may be achieved in several ways: tight channel spacing (e.g., [1,2]), advanced modulation formats (e.g., M-ary phase shift keying or quadrature-amplitude modulation), polarization multiplexing (e.g., [3]) or a combination of these techniques. The increase in spectral-efficiency is usually accompanied by higher sensitivity to degrading fiber propagation effects which stem from the interaction between dispersion and nonlinearity. Mitigation of nonlinear effects can be achieved by using dispersion-managed systems [4], but this solution introduces added noise and is quite rigid.
Alternatively, digital signal processing (DSP) may be used to achieve backwardpropagation, as explained in the following. The concept of backward-propagating an optical signal to mitigate transmission effect stems from Yariv et al. which in 1979 suggested the use of phase-conjugation for dispersion compensation [5]. Fisher et al. extended this notion in 1983 to compensate for both dispersion and Kerr nonlinearity [6]. The use of a medium with negative nonlinear index (e.g., semiconductors) to reverse the effects of transmission without phase conjugation was suggested in [7]. An examination of the nonlinear Schrödinger equation (NLSE), which governs the evolution of optical signal propagation in fiber with Kerr nonlinearity, reveals that all these techniques amount to reversing the fiber parameters (either dispersion, nonlinear index, or both). This is mathematically equivalent to reversing the spatial variable sign in the NLSE and hence the technique is referred to as backwardpropagation.
Electronic pre-distortion of the signal was considered in [8][9][10] and experimental demonstrations using dispersion-managed links with subsequent nonlinear phase fluctuation mitigation (lumped post-compensation) using DSP are presented in [11,12]. Postcompensation via backward-propagation takes into account all deterministic effects in fiber, allows the use of conventional transmitter architecture and offers flexibility through adaptive processing. Numerical simulations of a wavelength-division multiplexing (WDM) system with post-compensation of both linear and nonlinear effects were presented in [13]. An experimental demonstration of the feasibility of receiver-side backward-propagation in an orthogonal-WDM system environment soon followed [14]. Efficient real-time implementation of any DSP algorithm in a communication system context is crucial in order to reduce transmission latency and power consumption. Finite impulse response (FIR) filtering is highly compatible with real-time DSP implementation, as compared to other filtering techniques. It does, however, suffer from limited efficiency. In this paper, the authors suggest a wavelet-based FIR filter design used for backward-propagation, which is shown to achieve higher efficiency compared to alternative filter design procedures. In section 2 the method of achieving digital backward-propagation is presented. Section 3 explains the wavelet-based design procedure and in section 4 simulation results validating the efficiency of this approach are presented. Conclusions are presented in section 5.

Backward-propagation using the split-step method
Backward propagation is achieved using DSP after coherent detection (with which complete information of the received complex-field is preserved). After sampling at the appropriate rate, the received signal is backward-propagated by solving the NLSE (for a single polarization) with the spatial variable negated: where A is the complex electric field, D is the linear operator accounting for dispersion as well as fiber loss and N is the nonlinear operator, which takes Kerr nonlinearity into account. These operators are given by: where α β β 2 3 , , and γ are the attenuation factor, first-and second-order group-velocity dispersion and the nonlinearity parameter, respectively. The symmetric-SSM is used to solve N sequential nonlinear and linear full-step compensation operations follow and the last step consists of a single nonlinear (full-step) and linear (half-step) operations. This simple modification to the asymmetric-SSM achieves higher accuracy for a given step size [4]. Fiber loss is accounted for at each step, while amplification (required after each span to compensate for fiber loss of Backward-propagation for a single span using the symmetric-SSM is illustrated in Fig. 1. Accurate backward-propagation of a multichannel optical signal requires the SSM step size to be small enough so that the assumption that dispersion and nonlinearity act independently is valid [4]. A major difficulty in implementing this technique is the high computational load associated with the dispersion operator. This operator can be implemented in the frequency-domain (via fast Fourier transform, FFT) or in the time-domain using filters. For real-time compatible implementation, FIR filtering is a practical choice as it does not require block-processing (as opposed to FFT implementation). Moreover, FIR filtering does not require a feedback path, as opposed to the use of infinite-impulse response filtering, which was shown to be highly efficient [15], but is hard to implement in real-time.
The FIR filter coefficients may be obtained by an inverse-Fourier transform (IFT) of the dispersion frequency response of each step, given in Eq. (2a). This makes the filter length identical to the IFT length, which should be made long enough to achieve high accuracy of the SSM. To reduce the SSM computational load, time-domain windowing of the dispersion compensating FIR filter may be considered; e.g., using Tukey windowing [16]. The window width should be approximately on the order of the time-frame related to the dispersioninduced pulse spreading over a single step, given by: ( ) 1, f H ω ω = ∀ . As multiple steps of dispersion compensation are required by the SSM, the accumulated error from the windowing of the FIR filter becomes significant and the ultimate quality of the backward-propagated signal is deteriorated. IFT-based design of the FIR dispersion compensation filter is suboptimal since the Fourier transform relies on a harmonic decomposition basis which is not time-limited and windowing is hence required. Wavelet-based filter design of the dispersion operator is based on time-limited decomposition. This approach is shown to achieve accurate back-propagation while reducing the number of operations required for the linear operator, as shown in the following.

Wavelet-based FIR filter design
Wavelet-based FIR filter design makes use of compact-support decomposition functions. Compact support functions are only non-zero on a finite interval, as opposed to the harmonic basis. As suggested in [17], . The interpolating scaling function is defined by a filter whose coefficients are given by [17]: where k δ is the Kroneker delta function. The compact support of ( )  [18]. Keeping the same notation as in [17] ( 2, 3 µ = for the second and third order derivative, respectively), the following set of equations is to be solved:  [17] to introduce second-order group-velocity dispersion is included: In order for f h to obtain high accuracy, (2,3) D are zero-padded (before the Fourier transform) so that the impulse response is longer than the limit given by Eq.

Numerical simulation
To investigate the benefit of wavelet-based FIR filter design, simulation of a 9 channel WDM system was performed. Quadrature phase-shift keying (QPSK) modulation at 10GBaud, with 20GHz channel spacing was employed. The simulation sampling rate used was set to 320GHz , where s D is the dispersion slope. The loss in each span was compensated using an optical amplifier with noise figure 5 F N = . Verification of the benefit of using backward-propagation as compared to linear (dispersion) compensation only is seen in Fig. 3. The total launching power is scanned to demonstrate the ability to achieve better performance using backward-propagation, compared to linear compensation only. Backward-propagation is implemented as discussed in Section 2, with frequency domain dispersion compensation (SSM/FFT). Using frequency domain compensation no time-domain effects (such as FIR windowing) deteriorate the performance. SSM/FFT is hence used as the benchmark for backward-propagation performance. Also included is the performance limit when all the fiber parameters except loss are set to zero; in this case, only amplifier noise and linear interference affect the performance. The phase standard-deviation (Phase-STD) serves as a metric for the performance since phase shift keying is considered in this numerical study. From Fig. 3, the benefit of using backwardpropagation is clear as the optimal launching power is increased and the performance improved. The optimum total launching power is found to be 5dBm and an SSM step size of 2km was determined to achieve close to zero penalty performance, compared to smaller step sizes. These values will be used in the following analysis. A comparison of three FIR filters using IFT, IFT with Tukey windowing (both in frequency and time, as in [16]) and wavelet-based design is seen in Fig. 4. On the left, the time-domain magnitude, namely the impulse response is plotted on a logarithmic scale for clarity. The plot reveals the slow decay of the IFT filter; Tukey windowing limits the extent (and hence the number of operations required for filtering) of the IFT-based filter. Wavelets-based design achieves tighter windowing limits due to the compact-support of the decomposition function it is based upon. Moreover, observing the magnitude response comparison between IFT/Tukey and wavelets-based designs (Fig. 4, top  right) shows that the all-pass characteristic of the wavelet-based filter is preserved much better than with the IFT/Tukey design. Although the IFT/Tukey filter group-delay (Fig. 4, bottom right) is more accurate (namely, closer to the required linear group-delay response) at the band edges than the wavelet-based filter, loss of the all-pass feature combined with multiple linearcompensation operations required when employing the SSM is a significantly more severe issue. The number of FIR filter coefficients for IFT/Tukey was scanned by varying the Tukey window ratios of both frequency and time windows used. For the wavelets-based filter, the parameter s was scanned increasingly by taking values starting from lin N given in Eq. From Fig. 5, it is observed that the wavelets-based FIR filtering approach achieves equivalent performance as the IFT/Tukey design (both being close to the frequency-domain linear compensation implementation). The wavelet-based filter achieves this performance using 39 taps, as opposed to 70 taps for the IFT/Tukey design. The number of filter taps required for the SSM linear operator is hence reduced by a factor of approximately 1.8, validating the wavelet-based approach to dispersion-compensation FIR filter design. Although the wavelets-based is more efficient than IFT-based ones, the implementation of backwardpropagation leads to a larger computational load and latency, compared to lumped linear compensation only. The backward-propagation scheme's additional computational load and latency should be taken into account as an engineering trade-off with its superior performance, compared to linear compensation only.

Conclusions
In this contribution, the authors suggest an alternative approach to FIR filter design used for digital backward-propagation of an optical signal propagated in a nonlinear fiber. FIR filtering is highly compatible with real-time implementation, yet suffers from a large computational load requirement. The compact support of the wavelet decomposition basis function leads to a rapid decay of the FIR filter coefficients, as opposed to inverse-FT design. This rapid decay allows truncation of filter coefficients without severely affecting the all-pass property of the FIR filter used. Simulation was conducted for a WDM system with 9 channels, each QPSK modulated at 10GBaud , with 20GHz channels spacing. The optical signal was propagated in 48 spans of 100km NZ-DSF and the SSM was used for backward propagation. The waveletbased design achieves equivalent performance as IFT-based design with Tukey windowing, with a reduction of approximately 1.8 in filter length for the WDM configuration presented.