WO2022201763A1 - 光波形歪み補正方法及び装置、並びに光信号受信装置 - Google Patents

光波形歪み補正方法及び装置、並びに光信号受信装置 Download PDF

Info

Publication number
WO2022201763A1
WO2022201763A1 PCT/JP2022/000917 JP2022000917W WO2022201763A1 WO 2022201763 A1 WO2022201763 A1 WO 2022201763A1 JP 2022000917 W JP2022000917 W JP 2022000917W WO 2022201763 A1 WO2022201763 A1 WO 2022201763A1
Authority
WO
WIPO (PCT)
Prior art keywords
waveform
nonlinear
waveform distortion
phase modulation
parameters
Prior art date
Application number
PCT/JP2022/000917
Other languages
English (en)
French (fr)
Inventor
崇 井上
Original Assignee
国立研究開発法人産業技術総合研究所
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 国立研究開発法人産業技術総合研究所 filed Critical 国立研究開発法人産業技術総合研究所
Priority to US18/283,208 priority Critical patent/US20240178919A1/en
Priority to JP2023508663A priority patent/JPWO2022201763A1/ja
Publication of WO2022201763A1 publication Critical patent/WO2022201763A1/ja

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/25—Arrangements specific to fibre transmission
    • H04B10/2507—Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/25—Arrangements specific to fibre transmission
    • H04B10/2507—Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion
    • H04B10/2543—Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion due to fibre non-linearities, e.g. Kerr effect
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60—Receivers
    • H04B10/61—Coherent receivers
    • H04B10/616—Details of the electronic signal processing in coherent optical receivers
    • H04B10/6163—Compensation of non-linear effects in the fiber optic link, e.g. self-phase modulation [SPM], cross-phase modulation [XPM], four wave mixing [FWM]
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04J—MULTIPLEX COMMUNICATION
    • H04J14/00—Optical multiplex systems
    • H04J14/02—Wavelength-division multiplex systems
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04J—MULTIPLEX COMMUNICATION
    • H04J14/00—Optical multiplex systems
    • H04J14/02—Wavelength-division multiplex systems
    • H04J14/0201—Add-and-drop multiplexing

Definitions

  • the present invention is a technique for correcting signal waveform distortion caused by the nonlinearity of optical fibers in optical signal transmission using optical fibers.
  • Non-Patent Document 1 There are various factors that cause signal waveform distortion in optical fiber transmission lines. Of these, it is known that waveform distortion caused by a linear frequency response of a transmission line, such as group velocity dispersion, can be effectively corrected by a linear adaptive equalizer (see, for example, Non-Patent Document 1). ing. On the other hand, as a property of optical fiber, there is a phenomenon called the Kerr effect, in which the refractive index increases in proportion to the intensity of the incident light. Signal waveform is distorted. This distortion is hereinafter referred to as nonlinear waveform distortion.
  • WDM wavelength-division multiplexing
  • XPM Cross-Phase Modulation
  • DBP Digital Back Propagation
  • This method uses the following nonlinear Schrödinger equation to describe the propagation of light waves in optical fibers.
  • DSP digital signal processor
  • the Fourier transform is defined as follows as a premise of the formula (1).
  • one span of the transmission line meaning the section sandwiched between repeaters by optical amplification
  • ⁇ is divided into finite steps (sections), and at each step, the formula ( By alternately repeating the calculation (linear step) with only the linear term in 1) and the calculation with only the nonlinear term (nonlinear step), A p ( z, t).
  • This calculation method is called the split-step Fourier method (see, for example, Non-Patent Document 3). Increasing the number of steps improves the accuracy of calculation and improves the ability to correct nonlinear waveform distortion, but at the same time, the amount of calculation increases. increase. For this reason, it is preferable to find a method of improving the performance of correction while reducing the amount of calculation when implementing it in a DSP with limited computational resources.
  • Non-Patent Document 2 In the DBP proposed in Non-Patent Document 2, usually, only a single channel of the received signal waveform is cut out by a filter and only that waveform is back-propagated using equation (1). Although the distortion caused by it can be corrected, the distortion caused by XPM cannot be corrected for WDM signals.
  • Non-Patent Document 4 proposes a DBP that approximates the amount of phase shift caused by XPM for WDM signals and corrects not only SPM but also XPM distortion. However, in this method, if the number of steps per one span of the transmission line is set to less than 2, the amount of calculation becomes realistic, but the accuracy of the calculation is lowered, so there is a problem that the correction capability is significantly degraded.
  • the values of ⁇ 2 , ⁇ 3 , ⁇ , and ⁇ which are physical parameters of the optical fiber used as the transmission line, can be estimated with high accuracy and used for calculation. Desired. Of these, the values of ⁇ 2 , ⁇ 3 and ⁇ representing group velocity dispersion and propagation loss, which are linear responses, can be known with relatively high accuracy. In particular, a technique called OTDR (Optical Time-Domain Reflectometry) can be used as a method for measuring the distribution of ⁇ in the longitudinal direction. Also, although it is difficult to know the local values of ⁇ 2 and ⁇ 3 , it is possible to measure the integrated value of the measurement interval with high accuracy.
  • OTDR Optical Time-Domain Reflectometry
  • the nonlinear coefficient ⁇ it is not easy to directly measure the value of the nonlinear coefficient ⁇ , and the value can only be known indirectly by estimating it from the effective core cross-sectional area of the optical fiber. If the nonlinear waveform distortion is corrected by DBP using wrong parameters for the fiber, the correction effect will not be improved, and the waveform distortion may even increase.
  • Non-Patent Document 5 in a situation where the true values of the fiber parameters are unknown, DBP trials are repeated using the steepest descent method so that the correction result of the nonlinear waveform distortion is the best. A method for learning the optimal value of is proposed.
  • Non-Patent Document 6 using a neural network, the linear step of the split-step Fourier method is replaced with a time-domain finite impulse response (FIR: Finite Impulse Response) filter, assigned to the Affine transform, and the non-linear step is activated.
  • FIR Finite Impulse Response
  • the DBP structure using the nonlinear Schrödinger equation is applied to the neural network.
  • the FIR tap coefficients calculated from the group velocity dispersion are used as the initial values of the Affine transform coupling coefficients, and the neural network is configured to maximize the performance of the nonlinear waveform distortion correction.
  • Non-Patent Document 6 In neural networks, which are generally used for purposes such as image recognition, the coefficients of the Affine transform are initialized with random numbers, and a normalized linear function (ReLU: Rectified Linear Unit) is often used as the activation function.
  • ReLU Rectified Linear Unit
  • physical parameters are used as initial values, and a nonlinear term included in the physical evolution equation is applied to the activation function, so it is a neural network specialized for physical phenomena It can be said.
  • the performance of nonlinear waveform distortion correction is not degraded even if the number of steps per span is reduced.
  • Non-Patent Document 6 specifically reports that neural networks with 1 and 2 steps per span outperform conventional DBPs with 2 and 3 steps, respectively.
  • Non-Patent Documents 7 and 8 also report the results of nonlinear waveform distortion correction by a physical phenomenon-specific neural network under various conditions.
  • any method using a physical phenomenon-specific neural network only targets the correction of waveform distortion caused by SPM, and does not consider the distortion caused by XPM, so the distortion caused by XPM is the main cause of waveform deterioration.
  • effective correction performance cannot be exhibited.
  • Non-Patent Document 4 proposes a technique for correcting distortion due to XPM in addition to SPM. In addition, these methods have the problem that the original performance cannot be obtained unless the parameters of the transmission path are accurately estimated and input.
  • Non-Patent Document 6 by learning parameters, it is possible to exhibit optimal correction performance for nonlinear waveform distortion caused by SPM, but the distortion caused by XPM No effectiveness was found for
  • an object of the present invention is to provide a new technique for suppressing the amount of calculation for appropriately correcting waveform distortion caused by both SPM and XPM even if the parameters of the transmission path are unknown. technology.
  • the method for correcting optical waveform distortion according to the present invention includes (A) receiving an optical signal whose waveform shape has changed in a transmission line and digitizing the waveform, and then alternately calculating linear terms and nonlinear terms of the nonlinear Schrödinger equation.
  • Back propagation processing for estimating a waveform at the time of transmission by performing waveform distortion due to self-phase modulation occurring in each of a plurality of channels in the transmission path during wavelength division multiplexing, and A first parameter related to cross-phase modulation used in back-propagation processing for correcting waveform distortion caused by cross-phase modulation occurring between channels other than the channel, and a first parameter used in the back-propagation processing and self optimizing a second parameter related to phase modulation and cross-phase modulation by gradient descent; and (B) using the optimized first and second parameters to perform the back propagation process. and the step of
  • FIG. 1 is a diagram schematically showing the spectrum of a WDM signal.
  • FIG. 2 is a diagram for explaining the relationship between spans and steps.
  • FIG. 3 is a schematic diagram showing linear and non-linear steps and paths of each channel signal component.
  • FIG. 4A shows the exact path that the gradient by D 2 (0) takes from the linear step L 1,0 (0) to the subsequent stages.
  • FIG. 4B shows the approximate path that the gradient by D 2 (0) propagates from the linear step L 1,0 (0) to subsequent stages.
  • FIG. 5A is a diagram showing a flow of processing for calculating gradients of waveforms according to each transmission path parameter in parameter optimization processing.
  • FIG. 5B is a diagram showing a flow of processing for updating each transmission path parameter in parameter optimization processing.
  • FIG. 5A is a diagram showing a flow of processing for calculating gradients of waveforms according to each transmission path parameter in parameter optimization processing.
  • FIG. 5B is a diagram showing a flow of processing for updating each transmission path parameter in parameter optimization
  • FIG. 5C is a schematic diagram of an optical transmission system according to this embodiment.
  • FIG. 5D is a diagram showing a processing flow according to this embodiment.
  • FIG. 6 shows group velocity dispersion values and nonlinear constants in each span of a realistic transmission line.
  • FIG. 7A is a diagram showing the second-order group velocity dispersion parameter D 2 (j) with respect to the number of iterations of learning in the learning process in the case of an ideal 6-span transmission line and the span incident power of +2 dBm/ch.
  • FIG. 7B is a diagram showing the nonlinear coefficient parameter g (j) with respect to the number of repetitions of learning in the learning process in the case of an ideal 6-span transmission line and the span incident power of +2 dBm/ch.
  • FIG. 7A is a diagram showing the second-order group velocity dispersion parameter D 2 (j) with respect to the number of iterations of learning in the learning process in the case of an ideal 6-span transmission line and the span incident power of
  • FIG. 7C is a diagram showing d ⁇ 4 (j) with respect to the number of iterations of learning in the learning process in the case of an ideal 6-span transmission line and a span incident power of +2 dBm/ch.
  • FIG. 7D is a diagram showing d ⁇ 1 (j) with respect to the number of iterations of learning in the learning process in the case of an ideal 6-span transmission line and the span incident power of +2 dBm/ch.
  • FIG. 7E is a graph showing the mean squared error with respect to the number of repetitions of learning in the learning process in the case of an ideal 6-span transmission line and the span incident power of +2 dBm/ch.
  • FIG. 8A is a diagram showing the results of the Q value with respect to signal incident power when nonlinear waveform distortion is corrected in consideration of only SPM in ideal transmission lines of 4 spans and 6 spans.
  • FIG. 8B is a diagram showing the results of the Q value with respect to signal incident power when nonlinear waveform distortion is corrected in consideration of only SPM in ideal transmission lines of 8 spans and 10 spans.
  • FIG. 9A is a diagram showing the results of the Q value with respect to signal incident power when nonlinear waveform distortion is corrected including XPM in ideal transmission lines of 4 spans and 6 spans.
  • FIG. 9B is a diagram showing the result of the Q value with respect to the signal incident power when nonlinear waveform distortion is corrected including XPM in ideal transmission lines with 8 spans and 10 spans.
  • FIG. 10 is a diagram showing the results of the Q factor versus the number of transmission spans with and without nonlinear waveform distortion correction applied in various ways in an ideal transmission path.
  • FIG. 11A is a diagram showing the second-order group velocity dispersion parameter D 2 (j) with respect to the number of iterations of learning in the learning process when the span incident power is +2 dBm/ch in a practical 6-span transmission line.
  • FIG. 11B is a diagram showing the second-order group velocity dispersion parameter g (j) with respect to the number of iterative learning in the learning process when the span incident power is +2 dBm/ch in a practical 6-span transmission line.
  • FIG. 11A is a diagram showing the second-order group velocity dispersion parameter D 2 (j) with respect to the number of iterations of learning in the learning process when the span incident power is +2 dBm/ch in a practical 6-span transmission line.
  • FIG. 11B is a diagram showing the second-order group velocity dis
  • FIG. 11C is a diagram showing the second-order group velocity dispersion parameter d ⁇ 4 (j) with respect to the number of iterations of learning in the learning process when the span incident power is +2 dBm/ch in a practical 6-span transmission line.
  • FIG. 11D is a diagram showing the second-order group velocity dispersion parameter d ⁇ 1 (j) with respect to the number of iterations of learning in the learning process when the span incident power is +2 dBm/ch in a practical 6-span transmission line.
  • FIG. 11E is a graph showing the mean squared error with respect to the number of iterations of learning in the learning process when the span incident power is +2 dBm/ch in a realistic 6-span transmission line.
  • FIG. 12A is a diagram showing the results of the Q value with respect to signal incident power when nonlinear waveform distortion is corrected in consideration of only SPM in practical transmission paths of 4 spans and 6 spans.
  • FIG. 12B is a diagram showing the results of the Q value with respect to signal incident power when nonlinear waveform distortion is corrected in consideration of only SPM in realistic transmission paths of 8 spans and 10 spans.
  • FIG. 13A is a diagram showing the results of the Q value with respect to the signal incident power when the nonlinear waveform distortion is corrected including XPM in the realistic transmission lines of 4 spans and 6 spans.
  • FIG. 13B is a diagram showing the results of the Q value with respect to signal incident power when nonlinear waveform distortion is corrected including XPM in realistic transmission lines with 8 spans and 10 spans.
  • FIG. 14 is a diagram showing the results of the Q factor versus the number of transmission spans with and without applying nonlinear waveform distortion correction in various ways in a realistic transmission path.
  • FIG. 15 is a schematic diagram showing the transmission system used in the experiment.
  • FIG. 16A is a diagram showing changes in the second-order group velocity dispersion parameter D 2 (j) with respect to the number of parameter updates.
  • FIG. 16B is a diagram showing changes in the nonlinear coefficient parameter g (j) with respect to the number of parameter updates.
  • FIG. 16C is a diagram showing changes in the inter-polarization phase modulation parameter ⁇ (j) with respect to the number of parameter updates.
  • FIG. 16D is a diagram showing changes in the walk-off parameter d ⁇ 5 (j) with respect to the number of parameter updates.
  • FIG. 16E is a diagram showing changes in the walk-off parameter d ⁇ 1 (j) with respect to the number of parameter updates.
  • FIG. 16F is a diagram showing changes in the walk-off parameter d n (j) with respect to the number of parameter updates.
  • FIG. 16G is a diagram showing changes in Mean Square Error (MSE) and its moving average with respect to the number of parameter updates.
  • MSE Mean Square Error
  • FIG. 17 is a schematic diagram showing an optical fiber connection situation in a 6-span transmission line and a reverse propagation calculation order.
  • FIG. 18A is a diagram showing signal quality with respect to incident power when nonlinear waveform distortion correction is performed in consideration of only SPM in 6-span transmission lines and 8-span transmission lines.
  • FIG. 18B is a diagram showing signal quality with respect to incident power when nonlinear waveform distortion correction is performed in consideration of only SPM in 10-span transmission lines and 12-span transmission lines.
  • FIG. 18C is a diagram showing signal quality with respect to incident power when nonlinear waveform distortion correction is performed in consideration of only SPM in 14-span transmission lines and 16-span transmission lines.
  • FIG. 18A is a diagram showing signal quality with respect to incident power when nonlinear waveform distortion correction is performed in consideration of only SPM in 6-span transmission lines and 8-span transmission lines.
  • FIG. 18B is a diagram showing signal quality with respect to incident power when nonlinear waveform distortion correction is performed in consideration of
  • FIG. 19A is a diagram showing signal quality with respect to incident power when nonlinear waveform distortion correction is performed in consideration of SPM and XPM in 6-span transmission lines and 8-span transmission lines.
  • FIG. 19B is a diagram showing signal quality with respect to incident power when nonlinear waveform distortion correction is performed in consideration of SPM and XPM in 10-span transmission lines and 12-span transmission lines.
  • FIG. 19C is a diagram showing signal quality with respect to incident power when nonlinear waveform distortion correction is performed in consideration of SPM and XPM in 14-span transmission lines and 16-span transmission lines.
  • FIG. 20 is a diagram showing changes in signal quality with respect to transmission span.
  • FIG. 21 is a diagram showing the number of product operations with respect to the data length N in the cases of 5 channels, 11 channels and 21 channels.
  • FIG. 1 shows a schematic diagram of the spectrum of a WDM signal.
  • the envelope amplitude of Equation (1) is expressed as follows as the sum of the envelope amplitudes for each channel.
  • a p,n (t) is the baseband (meaning the center frequency is 0) envelope amplitude for the polarization component p of channel n.
  • Equation (8) expresses that the nonlinear effect in the nonlinear Schrödinger equation can be described in the form that the original nonlinear coefficient ⁇ 0 is attenuated in the longitudinal direction when the amplitude of the signal is attenuated by propagation loss.
  • the amplitude of channel number 0 is expressed as follows.
  • Equation (6) describes the linear step of the split-step Fourier method
  • Equation (11) describes the nonlinear step that considers the effects of both SPM and XPM for the signal waveform of channel number 0.
  • FIG. 2 A schematic diagram is shown. It should be noted that FIG. 2 can be easily generalized for different numbers of spans and steps.
  • the forward propagation direction is from left to right, and spans 1 and 2 are defined from the start point to the end point. Conversely, the direction from right to left is the back-propagation direction, and within each span, steps 1 and 2 are defined in the back-propagation direction.
  • the present embodiment can be applied as it is. Furthermore, the difference in the method of partitioning the steps within the span affects the results of waveform distortion correction in the methods shown in Non-Patent Documents 2 and 4 that do not perform learning, but the method that does not perform parameter learning affects the result of waveform distortion correction.
  • the effect is limited to the setting of the initial value of learning, and after the learning is converged and optimized, the difference disappears. Also, when the number of steps per span is 1, the entire 1 span is calculated as 1 step.
  • the waveform obtained as the output of the linear step is input, and the calculation of formula (11) is performed as a nonlinear step for the distance h of the entire step.
  • the waveform obtained as the output of the nonlinear step is used as an input, and the linear step calculation of Equation (6) is performed for the distance h/2 in the latter half of the step. This operation is repeated step by step, and the linear step in the latter half of one step is combined with the linear step in the first half of the next step to perform the calculation.
  • step #1 of span #2 is divided into two, and after performing calculation of the first half linear step L (0) on the input waveform, calculation of nonlinear step N (1) is performed, and then L (1) is calculated by combining the linear steps in the latter half of step #1 and the linear steps in the first half of step #2 of span #2, and similar calculations are performed thereafter.
  • FIG. 3 shows a backpropagation calculation procedure representing this situation.
  • linear and non-linear steps are alternately calculated.
  • the input waveform of channel number n and polarization component p is x p,n
  • the output waveform is y p,n .
  • x p,n is the transmitted waveform received by the receiver
  • y p,n corresponds to the estimated result of the transmitted waveform.
  • the input and output waveforms of the j-th linear step L (j) are zp ,n (j) and yp ,n ( j), respectively, the input and output waveforms of the nonlinear step N (j) are yp ,n (j ⁇ 1) and z p,n (j) .
  • the error function J( ⁇ ) is defined as follows.
  • ⁇ is a vector of parameters used in backpropagation calculation
  • y t ( ⁇ ) is the value of the signal waveform after backpropagation calculation at time t, taken as a function of the parameter ⁇
  • d t is the value of the desired signal at time t.
  • the waveform of the transmission signal is used as the desired signal.
  • the form of the error function in equation (16) is called the Mean Square Error (MSE).
  • MSE Mean Square Error
  • a data set is composed of multiple sets of the waveform d t of the transmission signal and the waveform y t ( ⁇ ) obtained after calculating the back propagation channel of the parameter ⁇ , and by the stochastic gradient descent method based on iterative calculation, Optimize the parameter ⁇ to minimize the error function J( ⁇ ).
  • an update formula for the parameter ⁇ based on the stochastic gradient descent method is obtained as follows. where ⁇ i is the parameter vector obtained as a result of the i -th update, ⁇ is a small positive number that determines the speed of learning, and ⁇ J is the gradient of the error function with respect to the parameter, and is expressed as follows: calculated to
  • the desired signal that is, the transmitted signal does not depend on the parameter ⁇
  • the parameter update formula is finally obtained as follows.
  • equation (20) when updating all the parameters in the backpropagation calculation that are elements of ⁇ , for the output signal yt, the gradient for each parameter, ⁇ yt / ⁇ , is use.
  • a formula for calculating this is derived from the differential chain rule.
  • Equation (6) the relationship between the input and output waveforms z p,n (j) and y p,n (j) at the linear step L p,n (j) is as follows.
  • the slopes of the output waveform y p,n (j) due to D 2 (j) and D 3 (j) can be calculated using the input waveform z p,n (j) .
  • the slope of y p,n (j) for any parameter ⁇ included in steps prior to L p,n (j) is obtained as follows.
  • the slope of the output waveform y p,n (j) with respect to the parameter ⁇ can be calculated using the slope ⁇ z p,n (j) / ⁇ output from the immediately preceding nonlinear step.
  • the walk-off value T n (j) is the walk - off parameter d n Using (j) and dn( j +1) , do the following: where h j is the interval width of the nonlinear step N p,n (j) .
  • Equation (11) the relationship between the input and output waveforms y p,n (j) and z p,n (j) at the nonlinear step N p,n (j) for the signal waveform of channel 0 is obtained from equation (11) as follows: become.
  • P p,n (j ⁇ 1)
  • 2 is the intensity of the input waveform y p,0 (j ⁇ 1) , and the following relationship is also be.
  • the output waveform z p , n ( _ The slope of j) is obtained by directly differentiating equation (26) as follows.
  • the output waveform z p,n (j) , the intensity of the input waveform P p,n (j ⁇ 1) , and its frequency waveform Pp,n (j ⁇ 1) (P is above P ) can be calculated using
  • equations (22) and (23) for linear steps and equations (30) to (34) for nonlinear steps are used to calculate the slope of the output waveform for that step due to the parameters used in that step, respectively. , passing the result to the next step.
  • dispersion parameters D 2 (j) and D 3 (j) are Also when calculating the gradient during learning, the gradients of both polarization components are averaged together.
  • parameters g (j) , ⁇ (j) , ⁇ 0 (j) , ⁇ n (j) , d n (j) shall be used in common.
  • nonlinear step calculations are not performed, and only linear step calculations are performed. If the calculation in the non-linear step is ignored, neither the XPM nor the SPM will be corrected. It has been confirmed that the nonlinear waveform distortion can be ignored.
  • the transmission path parameters for correcting nonlinear waveform distortion including XPM can be optimized by the gradient descent method, and the correction effect can be maximized.
  • Equation (23) shows a method of calculating the slope for the third-order group velocity dispersion of the transmission path. is so small that it can be ignored, it is possible to fix it without learning after setting it to a certain initial value, or to ignore the third-order group velocity dispersion effect itself.
  • the third-order group velocity dispersion is fixed without learning after setting the initial value.
  • the walk-off between channels changes quadratically with respect to the frequency difference between channels. will be taken into account when setting the initial values of the off-parameters d n (j) .
  • Equations (32) and (33) show how to calculate the slope for each channel's loss factor, but since the loss factor is a parameter that can be easily measured, it is fixed from the initial value without calculating the slope. Leave as is.
  • the output waveform z p,n ( 1) is calculated from equations (26) to (28), and the slope ⁇ z p,n (1) / ⁇ is calculated and sent to further steps.
  • the solid line in FIG. 4A indicates the gradient propagation path for D 2 (0) .
  • the propagation path of the gradient with a certain parameter ⁇ can be limited within the channel, resulting in a path like the solid line shown in FIG. 4B.
  • terms including H n (j) ( ⁇ ) in equation (36) are ignored, resulting in the following equation.
  • Equation (24) can be used as a gradient update formula for an arbitrary parameter ⁇ that has been propagated.
  • the waveform evolution is represented by equations (26) to (28), and the gradient calculation equation for the parameters g (j) , ⁇ (j) and d n (j) used in the step is represented by equation (30).
  • Equations (31) and (34) and equations (35) and (37) can be used as gradient update equations for any parameter ⁇ sent from the previous step.
  • FIG. 5A shows the flow of processing for calculating the slope of the waveform according to each transmission path parameter in the parameter optimization processing.
  • L (j) the evolution of the waveform from z p,n (j) to y p,n (j) is done in equation (21), but by the operations in equations (22) and (23) , L (j) for each transmission parameter (D 2 (j) and D 3 (j) ) of the waveform ( ⁇ y p,n (j) / ⁇ D 2 (j) and ⁇ y p,n (j) / ⁇ D 3 (j) ) is calculated.
  • the gradient of the waveform for all the transmission parameters ⁇ included in the steps before L (j) is calculated by the calculation according to Equation (24), and ⁇ z p,n (j) / ⁇ becomes ⁇ y p,n (j) is updated to / ⁇ .
  • ⁇ y p,n (j) / ⁇ D 2 (j) and ⁇ y p,n (j) / ⁇ D 3 (j) are included in ⁇ y p,n (j) / ⁇ and collectively It is sent to the next step, N (j+1) .
  • step N (j+1) the waveform evolution from y p,n (j) to z p,n (j+1) is done in equations (26)-(28), while equation (30) to ( 34 ) , the waveform _ _ gradient of ( ⁇ z p,n (j+1) / ⁇ g (j+1) , ⁇ z p,n (j+1) / ⁇ (j+1) , ⁇ z p,n (j+1) / ⁇ 0 (j+1) , ⁇ z p,n (j+1) / ⁇ n (j+1) , ⁇ z p,n (j+1) / ⁇ d n (j+1) ) are calculated.
  • equation (35) the slopes of the waveforms for all transmission parameters ⁇ included in the steps prior to N (j+1) are calculated by calculations according to equations (35) and (37), and ⁇ y p,n (j) / ⁇ is ⁇ zp ,n (j+1) / ⁇ .
  • equation (36) may be used instead of equation (37).
  • ⁇ zp ,n (j+1) / ⁇ g (j+1) , ⁇ zp ,n (j+1) / ⁇ (j+1) , ⁇ zp, n ( j+1)/ ⁇ 0 (j+1 ) , ⁇ zp , n (j+1) / ⁇ n (j+1) , ⁇ z p,n (j+1) / ⁇ d n (j+1) are included in ⁇ z p,n (j+1) / ⁇ and collectively in the next step sent to some L (j+1) .
  • ⁇ yt / ⁇ which is the gradient of the final output yt with respect to all parameters ⁇ in the transmission path, is calculated.
  • FIG. 5B shows a flow of processing for updating each transmission path parameter in parameter optimization processing. That is, using the output waveform y t , the gradient ⁇ y t / ⁇ ⁇ i by each transmission parameter ⁇ i , and the desired signal d t , the transmission path parameter ⁇ i at the i-th step is calculated using Equation (20) It is to be updated. As shown in FIG. 5B, in the linear step, the channel parameters D2 and D3 are updated, and in the nonlinear step, g, ⁇ , ⁇ 0 , ⁇ n , and dn are updated as the channel parameter ⁇ . .
  • FIG. 5C shows a schematic diagram of an optical transmission system according to this embodiment.
  • a WDM signal generated by a transmitter propagates through a transmission line composed of an optical fiber and an optical amplifier and reaches a receiver.
  • the receiver demultiplexes the WDM signal for each channel using a demultiplexer such as an arrayed waveguide grating.
  • the demultiplexed optical signal waveform of each channel is converted into an electrical signal waveform by a coherent receiver, and then input to a digital signal processor (DSP) for signal processing.
  • DSP digital signal processor
  • the electrical signal waveform is converted into numerical data by analog-to-digital (AD) conversion, and calculations for various demodulation processes are performed to finally convert it into a received bit string and output it.
  • AD analog-to-digital
  • Demodulation processing includes temporal sampling timing control, sampling rate conversion, clock synchronization, filtering, polarization rotation, carrier recovery, linear waveform distortion correction such as an adaptive equalizer, and nonlinear waveform distortion according to this embodiment. It includes correction, symbol determination, error correction processing, and the like.
  • the DSP can be equipped with a function for executing the parameter optimization process described above, but the parameter optimization process is performed by an external computer without being performed by the DSP, and the optimized It is also possible to download the parameters and apply them to the nonlinear waveform distortion correction according to this embodiment.
  • the AD conversion function may be separated from the DSP and installed between the DSP and the coherent receiver.
  • a DSP is provided independently for each channel, and basically the DSPs do not operate in cooperation with each other.
  • the above-mentioned XPM-based phase shift calculation is performed by collectively handling the waveform data of a plurality of channels by the processor of , or by mutually exchanging data between a plurality of independent processors.
  • FIG. 5D shows the processing flow performed by the DSP and related to the present embodiment.
  • a signal with a known waveform shape is transmitted, and based on the received signal waveform demodulated by the DSP in the receiver, a computer inside the DSP or an external computer executes parameter optimization processing to optimize the transmission path parameters (step S1).
  • the parameter optimization process is executed by a parameter optimization section configured within the DSP or an external computer.
  • the waveform distortion correction unit configured in the DSP executes back propagation processing using the optimized transmission path parameters (step S3). In this reverse propagation process, as described with reference to FIGS.
  • linear steps and nonlinear steps each have optimized transmission path parameters, and these
  • waveform distortion caused by self-phase modulation occurring within the channel and waveform distortion caused by cross-phase modulation occurring between channels other than the channel concerned are corrected, and the waveform at the time of transmission is corrected.
  • This backpropagation processing is performed by the DSP. By performing such processing, it is possible to improve accuracy while suppressing the amount of calculation.
  • the optical signal to be examined is a signal obtained by wavelength division multiplexing 9-channel wavelength division multiplexed 64-level quadrature amplitude modulation (QAM) signals with a symbol rate of 32 Gbaud (Dual-Polarization (DP)) at a frequency interval of 50 GHz. , to which random noise is added and the SNR is set to 25 dB. It is assumed that a root Nyquist filter with a roll-off coefficient of 0.05 is applied to the spectrum of the transmission signal. For 9-channel WDM signals, channel numbers are assigned from -4 to +4 from the lowest frequency side, attention is paid to the signal quality of channel number 0 in the center, and the operation of nonlinear waveform distortion correction is verified.
  • QAM quadrature amplitude modulation
  • a transmission line consisting of a standard single-mode fiber (SSMF) with a span length of 80 km and an optical amplifier that amplifies the propagation loss. Calculate the transmission of optical signals using paths.
  • the 2nd and 3rd order group velocity dispersion values are 16.641 ps/nm/km and 0.06 ps/nm 2 /km, respectively, and the nonlinear Let the coefficient be 1.3 W ⁇ 1 km ⁇ 1 and the propagation loss coefficient be 0.192 dB/km.
  • various parameters differ from span to span or from point to point, and the power of the optical signal fluctuates from the ideal state, resulting in the effects of group velocity dispersion and nonlinear effects. The size will vary.
  • a transmission line in which the second-order group velocity dispersion value and the nonlinear coefficient are varied step by step as shown in Fig. 6 is considered separately from the ideal transmission line.
  • the solid gray lines indicate the average values of the second-order group velocity dispersion value (a) and the nonlinear coefficient (b) for all parameters of 10 spans.
  • the average value within the existing spans when considering a short transmission line of 4 to 8 spans, consider the average value within the existing spans.
  • the numerical calculations for optical transmission simulation handle waveform data sampled at 32 times the symbol rate on the time axis.
  • the number of spans is set to 4, 6, 8, or 10, and the number of steps per span is set to 800 in an ideal transmission line with the same parameters in each span or a transmission line with the parameters shown in FIG.
  • the propagation of the optical signal in the transmission line is calculated according to the nonlinear Schrödinger equation of equation (1) to obtain output waveform data.
  • the value of the coefficient ⁇ corresponding to the inter-polarization cross-phase modulation is set to 1, and for the other parameters, values obtained by normalizing the parameters of the optical fiber constituting the transmission line are used.
  • the noise figure of the optical amplifier is assumed to be 6 dB, and random Gaussian distributed noise having noise power determined by the gain and noise figure is added to the signal as spontaneous emission noise during amplification.
  • the optical signal power at the start of propagation in each span is called "span incident power", and transmission simulation is performed for several span incident powers.
  • span incident power is small, although significant nonlinear waveform distortion does not occur, the SN ratio deterioration due to noise becomes significant, resulting in deterioration of signal quality.
  • the span incident power is increased, the SN ratio improves, but the nonlinear waveform distortion becomes conspicuous. The signal quality deteriorates.
  • transmit signal data with 16384 symbols is generated from a random bit pattern, a transmission simulation is performed, and the output waveform is saved.
  • a total of 200 data sets are used as data sets of pairs of input waveforms and output waveforms of the transmission line, and the gradient descent method using the above formula is used to optimize the correction of nonlinear waveform distortion.
  • Learn parameters In the backpropagation calculation involving parameter learning according to the present embodiment, the number of steps per span is set to one. On the other hand, in the backpropagation calculation without learning, both the case where the number of steps per span is 1 and the case where the number is 2 are considered, and the effect is compared with the case where learning is not performed.
  • ⁇ J i is the gradient of the error function with respect to the parameter ⁇ at the ith update
  • ⁇ in equation (38) is a learning coefficient, and a suitable value is set for each parameter.
  • Learning is performed by repeating the updating by the formula (38) 30,000 times. Since the number of data sets is 200, the order of the data sets is randomly changed every 200 updates, and learning is repeated. Note that this operation does not cause overlearning.
  • learning when performing similar learning while changing conditions such as the signal modulation format and the number of channels, learning may be performed 30,000 times or more until the parameters converge, if necessary.
  • y t and d t in equation (20) all 16384 symbols forming the data set are used.
  • another waveform with a different bit pattern and 262144 symbols is used to evaluate the signal quality after nonlinear waveform distortion correction.
  • FIGS. 7C and 7D show how the walk-off parameters d ⁇ 4 (j) and d ⁇ 1 (j) change, respectively. The values at all steps are almost the same, and almost no parameter change due to learning is observed. This is probably because in an ideal transmission path, correct parameters common to each span are given as initial values at the start of learning, so there is no need to adjust the walk-off by learning. Although not shown in the figure, other walk-off parameters converge similarly to d ⁇ 4 (j) and d ⁇ 1 (j) .
  • the received signal waveform and the desired signal that is, the MSE value calculated from the transmitted signal waveform are plotted. It can be seen that the MSE rapidly decreases immediately after the start of learning, converges when the number of iterations is about 1000, and then remains stable.
  • the MSE begins to oscillate near the minimum value, even if various parameters have not converged, updating is terminated, and steady waveform correction after learning is performed. may be started.
  • 8A and 8B show calculation of the Q value for each span incident power when the number of transmission spans is 4, 6, 8, and 10, and when nonlinear waveform distortion correction by back propagation calculation is performed and when it is not performed. shows the results.
  • the backpropagation calculation there are cases where the backpropagation calculation is performed without learning the parameters at 1step/span and 2steps/span, and after learning the parameters according to the present embodiment at 1step/span. , and all show the results when only the distortion due to SPM is corrected without correcting the non-linear waveform distortion due to XPM.
  • the Q value is obtained by converting the bit error rate (BER) value according to the following formula.
  • FIGS. 8A and 8B show that the Q value deteriorates due to the deterioration of the SN ratio in the region where the incident power is low, and conversely, the Q value deteriorates due to the nonlinear waveform distortion in the region where the incident power is large.
  • Focusing on the result of the backpropagation calculation without learning in the case of 1 step/span, there is only a slight correction effect, and in the case of 2 steps/span, a reasonable correction effect is obtained.
  • the number of steps per span is 1, signal quality slightly higher than the result of 2 steps/span without learning is obtained.
  • FIGS. 9A and 9B show the same results as FIGS. 8A and 8B, but show results when backpropagation calculations are performed taking into consideration XPM in addition to SPM.
  • the signal quality is improved by correcting the distortion by XPM.
  • the signal quality corrected by the 1 step/span back propagation calculation after the parameter learning according to the present embodiment is , surpasses the result of backpropagation calculation of 2 steps/span without learning, and is the best among them.
  • the method according to the present embodiment is sufficiently effective at 1 step/span, and can be realized with a realistic amount of calculation.
  • FIG. 4 shows the results of Q value versus number of spans under incident power. If the Q value threshold for error-free data reception is set to 7.5 dB, the maximum transmission distance without correction is 6 spans. It can be seen that by using the method of learning the parameters in consideration of , transmission over 10 spans or more is possible, and the transmission distance can be greatly extended.
  • FIG. 11A to FIG. 14 show results for a realistic transmission line whose parameters are shown in FIG.
  • FIGS. 11A to 11E show that parameters are updated 30,000 times in backpropagation calculations in which XPM is also taken into consideration for a data set with a realistic transmission line of 6 spans and a span incident power of +2 dBm/ch. It shows the fluctuation of various numerical values when However, as the initial values of the transmission path parameters at the start of learning, the average values of 6 spans are given assuming that the true parameters of the realistic transmission path shown in FIG. 6 are unknown.
  • FIG. 11A shows how the second-order group velocity dispersion parameter D 2 (j) changes, and FIG.
  • FIG. 11B shows how the nonlinear parameter g (j) changes, converging to different values for each span.
  • This result qualitatively agrees with the fact that the parameters of the original transmission line vary for each span, but the obtained value is the value that maximizes the nonlinear waveform distortion correction, and is not necessarily given in FIG. However, it does not match the transmission path parameters.
  • the walk-off parameters d ⁇ 4 (j) and d ⁇ 1 (j) respectively, converge to different values corresponding to different variance values for each span.
  • other walk-off parameters converge similarly to d ⁇ 4 (j) and d ⁇ 1 (j) . Similar to FIG. 7E, FIG.
  • 11E shows MSE values for the number of times of learning, but compared to the ideal transmission line shown in FIG. 7E, the results for the realistic transmission line in FIG. , the number of iterations converges at about 10,000. Similar to the result shown in FIG. 7E, in this case also, when the MSE starts to oscillate near the lowest value, even if various parameters do not converge, the update is terminated, and steady waveform correction after learning is performed. may be started.
  • FIGS. 12A and 12B similarly to FIGS. 8A and 8B, respectively, the nonlinear waveform distortion correction by back propagation calculation is performed and not performed. It shows the calculation result of the Q value for each span incident power.
  • the backpropagation calculation there are cases where the backpropagation calculation is performed without learning the parameters at 1 step/span and 2 steps/span, and after learning the parameters according to the present embodiment at 1 step/span. Calculations are performed, and the results are shown when only the distortion due to SPM is corrected without correcting the non-linear waveform distortion due to XPM.
  • the backpropagation calculation without learning the results are shown both when the true values shown in FIG.
  • FIGS. 13A and 13B show results similar to FIGS. 12A and 12B when nonlinear waveform distortion correction is performed in consideration of not only SPM but also XPM. Even in a realistic transmission path, the signal quality is improved by correcting the distortion by XPM. The corrected signal quality outperforms the 2steps/span backpropagation calculation with no learning and true channel parameters, and is the best among them. 13A and 13B, it should be noted that in the backpropagation calculation without learning, the result when the average value of the true values is given as the parameter of the transmission path is It is that it is remarkably deteriorated compared with the result.
  • FIG. 14 shows the results of nonlinear waveform distortion correction shown in FIGS. 12A and 12B and FIGS. 13A and 13B. , the results of the Q value versus the number of spans under the optimum incident power are shown. Similar results to the ideal transmission path were obtained for a realistic transmission path. can be extended significantly.
  • FIG. 15 shows a circular transmission experimental system, and the details of the experimental procedure using this system and the experimental results will be described below.
  • 11 channels of continuous light with different wavelengths output from the wavelength tunable light source are multiplexed by a 16 ⁇ 1 polarization maintaining coupler and then lithium niobate (LN) polarization multiplexing. It is input to the IQ modulator, modulated into an optical signal waveform by an electrical signal applied to the modulator, and output.
  • the 4-channel electrical signals applied to the modulator are generated by an arbitrary waveform generator with a sampling rate of 64 GSample/s. It corresponds to a polarization Q-channel component, a Y-polarization I-channel component, and a Y-polarization Q-channel component. These 4-channel electrical signals are amplified by driver amplifiers and applied to the modulator.
  • the optical signal output from the modulator is a polarization multiplexed QAM signal having a symbol rate of 32 Gbaud and a root Nyquist waveform with a roll-off coefficient of 0.1.
  • a uniformly distributed 16QAM signal in which the number of bits per single polarization single symbol is 4 bits, and a probability distribution (PS: Probabilistically Shaping) signal in which the number of bits per single polarization single symbol is 5 bits.
  • PS Probabilistically Shaping
  • a 64QAM signal is used.
  • the 11-channel WDM signal obtained by modulation is modulated to the same waveform in all channels, but by long-distance transmission through a transmission line in which group velocity dispersion exists, walk-off (group delay between channels) , resulting in random inter-waveform XPM after long-distance transmission.
  • the 11-channel WDM signal generated by Tx has its optical power amplified by an optical amplifier, adjusted by a variable optical attenuator (VOA), and passed through a band-pass filter (BPF). Filter) removes optical noise outside the signal band, and the noise is input to an acoustic optical modulator (AOM), which is a switch for switching cyclic transmission operations. Signal input/output to/from the loop transmission line is performed through a 3 dB coupler.
  • VOA variable optical attenuator
  • BPF band-pass filter
  • AOM acoustic optical modulator
  • the loop transmission line consists of an optical amplifier, BPF, VOA, 84.1 km long SSMF, optical amplifier, BPF, VOA, 80.5 km long SSMF, optical amplifier, isolator (right arrow), It consists of a polarization scrambler (Pol. Scrambler) and an AOM, and includes an SSMF with two spans per turn.
  • the transmission distance is set to any one of 6, 8, 10, 12, 14, and 16 spans, and signals are transmitted for each distance to verify the effect of nonlinear waveform distortion correction.
  • the SSMF with a length of 84.1 km and the SSMF with a length of 80.5 km have slightly different dispersion characteristics.
  • the average characteristic of one round of the circuit transmission line is the group velocity dispersion A value of 16.85 ps/nm/km and a dispersion slope of 0.060 ps/nm 2 /km were estimated.
  • Rx is a digital coherent receiver consisting of a 4-channel real-time oscilloscope with an electrical bandwidth of 33 GHz and 80 GSample/s, a local oscillator (LO) and an optical front end.
  • LO local oscillator
  • Normal demodulation processing up to signal quality evaluation without performing back propagation calculation includes dispersion compensation, application of the same root Nyquist filter as that applied as a matched filter during transmission, polarization rotation and polarization multiplexing components. Separation, resampling to 2 samples/symbol, retiming, carrier frequency estimation and compensation, carrier phase recovery and 3-tap feedforward linear adaptive equalization, symbol decision and bit pattern acquisition.
  • the 3-tap adaptive equalization processing is butterfly-type 2 ⁇ 2 MIMO processing that supports polarization multiplexed signals, and compensates for polarization crosstalk caused by birefringence and XPM in the transmission line. It is effective for
  • the waveform of the signal demodulated in this way has changed greatly from the waveform at the time of reception due to the application of the root Nyquist filter, and the back propagation calculation cannot be applied as it is. Therefore, without applying the root Nyquist filter to the waveform after application of dispersion compensation in the normal demodulation processing, the processing performed in the subsequent demodulation process is performed in the same way, and the waveform immediately after reception is reproduced as much as possible. faithfully reproduced.
  • This walk-off value corresponds to approximately 1400 symbols for a 32 Gbaud signal, so if it is measured under conditions that include more symbols than this, the backpropagation using a data set with a finite time width will can correctly calculate the effects of walk-off. Therefore, a data set is constructed by adding 5000 symbols to each end of the time axis for a waveform having one period of 65536 symbols.
  • one of the four data sets of the DP-16QAM signal is randomly selected, back propagation calculation is performed, and consecutive among 65536 symbols
  • the process of learning is randomized even for a limited amount of data sets, and it becomes possible to proceed with learning without causing over-learning.
  • the above-mentioned AdaBelief is used as the gradient descent algorithm used for learning.
  • FIGS. 16A to 16G show, when the number of transmission spans is 6, under the back propagation condition of 1 step/span, in the process of learning transmission path parameters for nonlinear waveform distortion correction considering XPM compensation according to the embodiment. , shows how various parameters are updated.
  • FIG. 17 shows the state of the connection of the optical fibers constituting the 6-span transmission line used in the experiment, the fiber length of each span, and the linear step L (j) (0 ⁇ j ⁇ 6) and the nonlinear step in the backward propagation calculation.
  • the calculation order of N (j) (1 ⁇ j ⁇ 6) is represented as a schematic diagram. From the results of FIGS.
  • the channel parameters and the signal quality converge after updating the parameters a certain number of times.
  • FIG. 20 plots the results of signal quality versus the number of transmission spans for each type of nonlinear waveform distortion correction, and experimental results similar to simulation results are obtained. As described above, the experimental results also revealed that the effects of the embodiment were obtained.
  • the amount of calculation required when performing nonlinear waveform distortion correction according to the embodiment is compared with the amount of calculation required when using the conventional technique.
  • the amount of computation required for the method related to the form is comparable to that required for the conventional method.
  • the nonlinear waveform distortion correction according to the embodiment is to consider both the distortion due to SPM and XPM, perform back propagation calculation at 1 step/span, and use the gradient descent method to determine the optimum value of the transmission path parameter. After learning, the waveform is corrected by back-propagation calculation.
  • attention is focused on the amount of calculation when the back-propagation calculation is performed with the optimized parameters fixed after the learning is completed.
  • the amount of calculation at this time is the same as in the case of backpropagation calculation in which XPM is considered at 1 step/span and learning is not performed. Further, the conventional technique is assumed to be one in which only correction of distortion by SPM is taken into account, and correction of distortion by XPM is not performed, and backward propagation calculation is performed at 2 steps/span.
  • numerical parameters that can be fixed regardless of the input waveform are calculated in advance and the results are saved in a look-up table (LUT). However, it is assumed that each different waveform input is read out and used, and the number of calculations required for the calculation itself is not considered.
  • the number of times of real number integration is 4N when calculating B p,0 (0, t) ⁇ ei ⁇ , which is the product of complex numbers.
  • P p,0 (t) Re[B p,0 (t)] 2 + Im [B p,0 (t )] 2
  • the number of integration times of the real numbers is 2N
  • the total number of integration times of the real numbers is 3N.
  • Equation (13) the amount of calculation required to calculate the phase shift ⁇ XPM (t) due to XPM is estimated from equation (13).
  • the intensity waveform P p,0 (t) is already obtained when computing ⁇ SPM (t), and the FFT is applied to it to compute P p,0 ( ⁇ ).
  • the number of times is 2N(log 2 N ⁇ 2).
  • P ⁇ p,0 ( ⁇ ) does not appear in equation (13), it takes computational complexity into account since it needs to be supplied to perform nonlinear waveform distortion correction for other channels.
  • P ⁇ p,n ( ⁇ ) for n ⁇ 0 appearing in Equation (13) is assumed to be separately calculated and supplied, and the necessary computational complexity is not considered.
  • this embodiment can be implemented with substantially the same amount of calculation as the conventional method.
  • the present embodiment is slightly larger than the conventional method, which is 1.6 times, but the order of calculation is almost the same. As shown in the experimental result of the case, a large nonlinear waveform distortion correction effect is obtained.
  • the number of channels is 21, the amount of calculation required for this embodiment is increased by about 2.4 times that of the conventional method, but it can be said that it remains on the same order.
  • a DSP includes a calculator and a memory. Further, the processing described above is not limited to the DSP, and may be another processor. Furthermore, a program that causes the processor to execute the processes described above is recorded in a non-volatile memory, and when executed, instructions included in the program are read out and executed by the processor. Furthermore, a dedicated circuit or a combination of a dedicated circuit and a DSP or the like may be used to execute the processing described above.
  • An optical waveform distortion correction method receives an optical signal whose waveform shape has changed due to the nonlinear optical effect and group velocity dispersion effect of an optical fiber as a transmission line, and digitizes the waveform. After that, the linear term and the nonlinear term of the nonlinear Schrödinger equation are alternately calculated to estimate the waveform at the time of transmission, thereby correcting the distortion of the optical waveform. Not only waveform distortion caused by modulation but also waveform distortion caused by cross-phase modulation between channels during wavelength division multiplexing transmission is taken into account in the calculations. It is characterized in that the number of steps per span is set to less than two. It is possible to improve the accuracy while suppressing the calculation load.
  • the number of steps per one span of the transmission line described above may be 1 or less. Even if the number of steps is reduced in this way, the accuracy can be improved.
  • parameters mentioned above may include second-order group velocity dispersion, nonlinear coefficients and walk-off.
  • the method for correcting optical waveform distortion includes (A) receiving an optical signal whose waveform shape has changed in a transmission line, digitizing the waveform, and calculating the linear term of the nonlinear Schrodinger equation.
  • Back-propagation processing for estimating a waveform at the time of transmission by alternately calculating nonlinear terms and self-phase modulation occurring in each of a plurality of channels in the transmission path during wavelength division multiplexing transmission (SPM) and waveform distortion caused by cross phase modulation (XPM) occurring between channels other than the channel used in back propagation processing and related to cross phase modulation and a second parameter used in the backpropagation process and related to self-phase modulation and cross-phase modulation, by gradient descent method; (B) optimized first and first and performing the backpropagation process using 2 parameters.
  • SPM wavelength division multiplexing transmission
  • XPM cross phase modulation
  • the second parameters eg D 2 , D 3 , g, ⁇ , ⁇ 0
  • the first parameters eg ⁇ n , d n
  • the optimization By executing back-propagation processing that corrects waveform distortion caused by SPM and XPM using the first and second parameters that have been obtained, the number of steps per transmission line span can be reduced and the calculation load can be suppressed. , it is possible to obtain sufficient calculation accuracy.
  • the second parameter may include group velocity dispersion D2 and nonlinear coefficient g, and the first parameter may include walk-off parameter dn .

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Nonlinear Science (AREA)
  • Optical Communication System (AREA)

Abstract

伝送路のパラメータが未知であっても、ＳＰＭ及びＸＰＭの両方に起因する波形歪みを適切に補正するための計算量を抑制するため、本方法は、非線形シュレディンガー方程式の線形項と非線形項とを交互に演算することで送信時の波形を推定する逆伝搬処理であって、波長分割多重伝送時における上記伝送路内の複数のチャネルの各々について、当該チャネル内で生ずるＳＰＭに起因する波形歪みと、当該チャネル以外のチャネルとの間で生じるＸＰＭに起因する波形歪みとを補正する逆伝搬処理において用いられ且つＸＰＭに関連する第１のパラメータと、ＳＰＭ及びＸＰＭに関連する第２のパラメータとを、勾配降下法によって最適化する工程と、第１及び第２のパラメータを用いて、上記逆伝搬処理を実行する工程とを含む。

Description

光波形歪み補正方法及び装置、並びに光信号受信装置
　本発明は、光ファイバを用いた光信号伝送において、光ファイバの非線形性に起因する信号波形の歪みを補正する技術である。
　光ファイバ伝送路には、信号の波形歪みを引き起こす様々な要因がある。このうち、群速度分散などの伝送路の線形周波数応答によって生じる波形歪みについては、線形の適応等化器（例えば非特許文献１を参照のこと）によって効果的に波形歪みを補正できることが知られている。一方で、光ファイバの性質として、屈折率が入射光の強度に比例して大きくなるという、カー効果と呼ばれる現象があって、この効果と群速度分散効果の組み合わせによって、光ファイバを伝搬する光信号の波形が歪む。以下ではこの歪みを非線形波形歪みという。波長分割多重（ＷＤＭ：Wavelength-Division Multiplexing）によって周波数あるいは波長が異なる複数の光信号を多重して光ファイバを伝送する場合、伝搬中にカー効果が発生することによって、単一チャネルの波形内で生じる位相シフトを自己位相変調（ＳＰＭ：Self-Phase Modulation）といい、異なるチャネルの波形同士の相互作用によって引き起こされる位相シフトを相互位相変調（ＸＰＭ：Cross-Phase Modulation）という。これらの位相シフトと群速度分散の効果とが組み合わさって生じる非線形波形歪みを、線形の適応等化器で効果的に補正することはできない。
　非線形波形歪みを補正する手法として、デジタル逆伝搬法（ＤＢＰ：Digital Back Propagation）が提案された（例えば非特許文献２を参照のこと）。この方法は、光ファイバ中の光波伝搬を記述する以下の非線形シュレディンガー方程式を用いる。
Figure JPOXMLDOC01-appb-M000001
 　そして、受信機内のデジタル信号処理プロセッサ（ＤＳＰ：Digital Signal Processor）によって、光信号受信波形を伝送路の逆方向に伝搬させる計算を行うことで、送信時の波形を推定して波形歪みを補正するものである。ここでＡｐ（ｚ，ｔ）は、直交する二つの偏波成分（ｐ＝１，２）に対応する光信号波形の複素包絡線振幅であり、ファイバの長手方向についての距離ｚと、時間ｔの関数である。なお、式（１）の前提として、フーリエ変換は、以下のように定義している。
Figure JPOXMLDOC01-appb-M000002
　β２及びβ３は、それぞれ光ファイバの二次及び三次の群速度分散、αは伝搬損失係数である。また、γ０は非線形係数、そしてδは偏波間位相変調の程度を表す係数で、通常はδ＝１とおかれる。式（１）の逆伝搬計算を行う際には、伝送路の１スパン（光増幅による中継器に挟まれた区間を意味する）を有限のステップ（区間）に区切り、各ステップでは、式（１）で線形項のみとした場合の計算（線形ステップ）と、非線形項のみとした場合の計算（非線形ステップ）を交互に繰り返すことで、近似的に式（１）の解としてのＡｐ（ｚ，ｔ）を計算する。
　この計算方法はスプリットステップフーリエ法（例えば非特許文献３を参照のこと）と呼ばれ、ステップ数を増やせば計算の精度が高まって非線形波形歪みを補正する能力も向上するが、同時に計算量が増大する。このため、計算資源の限られたＤＳＰに実装するにあたっては、計算量を低減しつつ、補正の性能を高める手法を見い出すことが好ましい。
　非特許文献２で提案されたＤＢＰでは、通常、受信信号波形の単一チャネルのみをフィルタで切り出し、式（１）を用いてその波形だけを逆伝搬させるため、非線形波形歪みのうちのＳＰＭに起因する歪みを補正することはできるが、ＷＤＭ信号に対してＸＰＭに起因する歪みを補正することはできない。一方、非特許文献４では、ＷＤＭ信号を対象として、ＸＰＭによって生じる位相シフトの量を近似的に計算し、ＳＰＭだけでなくＸＰＭによる歪みも補正するＤＢＰが提案されている。しかしながら、この方法では、伝送路１スパンあたりのステップ数を２未満にすると、計算量は現実的なものとなるが、計算精度が落ちるため、補正能力が著しく劣化するという問題がある。
　一方で、ＤＢＰによる非線形波形歪み補正を実施するにあたっては、伝送路である光ファイバが有する物理パラメータであるβ２、β３、α、γの値を高い精度で推定し、計算に用いることが求められる。このうち、線形応答である群速度分散と伝搬損失を表すβ２、β３及びαについては、比較的精度よく値を知ることができる。特に、αの長手方向の分布を測定する方法として、ＯＴＤＲ（Optical Time-Domain Reflectometry）と呼ばれる手法が利用可能である。また、β２及びβ３については局所的な値を知ることは難しいものの、測定区間の積分値を高精度に測定することは可能である。ところが、非線形係数γの値を直接測定することは容易ではなく、光ファイバの実効コア断面積から推定するなどで間接的にしか値を知ることができない。ファイバのパラメータとして誤ったパラメータを用いてＤＢＰによる非線形波形歪み補正を実施した場合、補正効果があがらないどころか、かえって波形歪みが増大することすらある。
　そこで非特許文献５では、ファイバパラメータの真値が不明な状況において、非線形波形歪みの補正結果が最良となるように、最急降下法を用いてＤＢＰの試行を繰り返すことで、ＤＢＰで用いるファイバパラメータの最適値を学習する手法が提案されている。
　また非特許文献６では、ニューラルネットワークを用い、スプリットステップフーリエ法の線形ステップを時間領域の有限インパルス応答（ＦＩＲ：Finite Impulse Response）フィルタに置き換えたうえでAffine変換に割り当て、さらに非線形ステップを活性化関数に割り当てることで、非線形シュレディンガー方程式を用いたＤＢＰの構造をニューラルネットワークに適用している。伝送路の構成から推定されるファイバパラメータのうち、群速度分散から計算されるＦＩＲタップ係数をAffine変換結合係数の初期値に用い、非線形波形歪み補正の性能が最大化されるようにニューラルネットワークを学習することで、光ファイバのパラメータが未知の状況で、群速度分散に限らない任意の線形応答が伝送路に含まれる場合においても、効果的な非線形波形歪み補正が可能となる。
　一般に画像認識等を目的として使用されるニューラルネットワークでは、Affine変換の係数は乱数で初期化され、さらに活性化関数には正規化線形関数（ＲｅＬＵ：Rectified Linear Unit）がよく用いられるのに対して、非特許文献６の方法では、物理的なパラメータを初期値として用い、さらに物理的な発展方程式に含まれる非線形項を活性化関数に適用しているため、言わば物理現象に特化したニューラルネットワークと言えるものである。この物理現象特化型ニューラルネットワークで得られる付加的な効果として、スパンあたりのステップ数を削減しても、非線形波形歪み補正の性能が劣化しないことが挙げられる。非特許文献６では具体的に、スパンあたりのステップ数が１及び２のニューラルネットワークが、それぞれステップ数が２及び３の従来型ＤＢＰの補正能力を上回るという結果が報告されている。
　このような結果を受けて、非特許文献７及び８でも、様々な条件で物理現象特化型ニューラルネットワークによる非線形波形歪み補正の結果が報告されている。しかしながら、物理現象特化型ニューラルネットワークを用いたいずれの手法もＳＰＭによって生じる波形歪みの補正のみを対象としており、ＸＰＭによる歪みを考慮していないため、ＸＰＭによる歪みが主な波形劣化要因となるＷＤＭ伝送系においては、有効な補正性能を発揮することができていない。
　非特許文献２で提案されたＤＢＰでは、ＳＰＭによる波形歪みのみを補正の対象としているため、ＸＰＭによる歪みの補正を行うことができない。非特許文献４ではＳＰＭに加えてＸＰＭによる歪みを補正する手法が提案されたが、１スパンを区切るステップ数を多くしないと性能が発揮できず、計算量が増大するという問題がある。また、これらの手法では、伝送路のパラメータを精度よく推定して入力しなければ、本来の性能が得られないという問題がある。
　さらに、非特許文献６で提案された物理現象特化型ニューラルネットワークでは、パラメータの学習を行うことで、ＳＰＭに起因する非線形波形歪みに対しては最適な補正性能を発揮できるが、ＸＰＭによる歪みに対しては有効性を見いだせない。
特開２０２０－１４５５６１号公報
S.Haykin, "Adaptive Filter Theory," Pearson (2013) E. Ip and J. M. Kahn, "Compensation of Dispersion and Nonlinear Impairments Using Digital Backpropagation," J. Lightw. Technol., vol.26, no.20, pp.3416-3425 (2008) G. P. Agrawal, "Nonlinear Fiber Optics," Academic Press (2001) E. F. Mateo, F. Yaman, and G. Li, "Efficient compensation of inter-channel nonlinear effects via digital backward propagation in WDM optical transmission," Opt. Express, vol.18, no.14, pp.15144-15154 (2010) T. Tanimura, T. Hoshida, T. Tanaka, L. Li, S. Oda, H. Nakashima, Z. Tao, and J. C. Rasmussen, "Semi-blind Nonlinear Equalization in Coherent Multi-Span Transmission System with Inhomogeneous Span Parameters," Proceedings of OFC/NFOEC2010, Paper OMR6 C. Hager and H. D. Pfister, "Nonlinear Interference Mitigation via Deep Neural Networks," Proceedings of OFC2018, Paper W3A.4 Q. Fan, G. Zhou, T. Gui, C. Lu, A. P. T. Lau, "Advancing theoretical understanding and practical performance of signal processing for nonlinear optical communications through machine learning," Nat. Commun, vol.11, 3694 (2020) B. Bitachon, A. Ghazisaeidi, M. Eppenberger, B. Baeuerle, M. Ayata, and J. Leuthold, "Deep learning based digital backpropagation demonstrating SNR gain at low complexity in a 1200km transmission link," Opt. Express, vol.28, no.20, pp.29318-29334 (2020) J. Zhuang, T. Tang, Y. Ding, S. Tatikonda, N. Dvornek, X. Papademetris, and J. S. Duncan, "AdaBelief Optimizer: Adapting Stepsizes by the Belief in Observed Gradients," arXiv, 2010.07468 (2020)
　従って、本発明の目的は、一側面によれば、伝送路のパラメータが未知であっても、ＳＰＭ及びＸＰＭの両方に起因する波形歪みを適切に補正するための計算量を抑制するための新たな技術を提供することである。
　本発明に係る光波形歪み補正方法は、（Ａ）伝送路において波形の形状が変化した光信号を受信して波形を数値化した後、非線形シュレディンガー方程式の線形項と非線形項とを交互に演算することで送信時の波形を推定する逆伝搬処理であって、波長分割多重伝送時における上記伝送路内の複数のチャネルの各々について、当該チャネル内で生ずる自己位相変調に起因する波形歪みと、当該チャネル以外のチャネルとの間で生じる相互位相変調に起因する波形歪みとを補正する逆伝搬処理において用いられ且つ相互位相変調に関連する第１のパラメータと、上記逆伝搬処理において用いられ且つ自己位相変調及び相互位相変調に関連する第２のパラメータとを、勾配降下法によって最適化するステップと、（Ｂ）最適化された第１及び第２のパラメータを用いて、上記逆伝搬処理を実行するステップとを含む。
図１は、ＷＤＭ信号のスペクトルを模式的に示す図である。 図２は、スパンとステップとの関係を説明するための図である。 図３は、線形ステップ及び非線形ステップ、そして各チャネル信号成分の経路を示す模式図である。 図４Ａは、Ｄ２ （０）による勾配が線形ステップＬ１，０ （０）を発してから後段に波及する厳密な経路を示す図である。 図４Ｂは、Ｄ２ （０）による勾配が線形ステップＬ１，０ （０）を発してから後段に波及する近似経路を示す図である。 図５Ａは、パラメータ最適化処理において、各伝送路パラメータによる波形の勾配を計算する処理のフローを示す図である。 図５Ｂは、パラメータ最適処理において、各伝送路パラメータを更新する処理のフローを示す図である。 図５Ｃは、本実施の形態に係る光伝送システムの模式図である。 図５Ｄは、本実施の形態に係る処理フローを示す図である。 図６は、現実的な伝送路の各スパンにおける群速度分散値と非線形定数とを示すものである。 図７Ａは、理想的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対する２次の群速度分散パラメータＤ２ （ｊ）を表す図である。 図７Ｂは、理想的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対する非線形係数パラメータｇ（ｊ）を表す図である。 図７Ｃは、理想的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対するｄ－４ （ｊ）を表す図である。 図７Ｄは、理想的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対するｄ－１ （ｊ）を表す図である。 図７Ｅは、理想的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対する平均二乗誤差を表す図である。 図８Ａは、４スパン及び６スパンの理想的な伝送路で、ＳＰＭのみを考慮して非線形波形歪みを補正した場合に、信号入射パワーに対するＱ値の結果を表す図である。 図８Ｂは、８スパン及び１０スパンの理想的な伝送路で、ＳＰＭのみを考慮して非線形波形歪みを補正した場合に、信号入射パワーに対するＱ値の結果を表す図である。 図９Ａは、４スパン及び６スパンの理想的な伝送路で、ＸＰＭを含めて非線形波形歪みを補正した場合に、信号入射パワーに対するＱ値の結果を表す図である。 図９Ｂは、８スパン及び１０スパンの理想的な伝送路で、ＸＰＭを含めて非線形波形歪みを補正した場合に、信号入射パワーに対するＱ値の結果を表す図である。 図１０は、理想的な伝送路において様々な方法で非線形波形歪み補正を適用した場合としない場合に、伝送スパン数に対するＱ値の結果を表す図である。 図１１Ａは、現実的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対する２次の群速度分散パラメータＤ２ （ｊ）を表す図である。 図１１Ｂは、現実的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対する２次の群速度分散パラメータｇ（ｊ）を表す図である。 図１１Ｃは、現実的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対する２次の群速度分散パラメータｄ－４ （ｊ）を表す図である。 図１１Ｄは、現実的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対する２次の群速度分散パラメータｄ－１ （ｊ）を表す図である。 図１１Ｅは、現実的な６スパン伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合の学習過程において、繰り返し学習回数に対する平均二乗誤差を表す図である。 図１２Ａは、４スパン及び６スパンの現実的な伝送路で、ＳＰＭのみを考慮して非線形波形歪みを補正した場合に、信号入射パワーに対するＱ値の結果を表す図である。 図１２Ｂは、８スパン及び１０スパンの現実的な伝送路で、ＳＰＭのみを考慮して非線形波形歪みを補正した場合に、信号入射パワーに対するＱ値の結果を表す図である。 図１３Ａは、４スパン及び６スパンの現実的な伝送路で、ＸＰＭを含めて非線形波形歪みを補正した場合に、信号入射パワーに対するＱ値の結果を表す図である。 図１３Ｂは、８スパン及び１０スパンの現実的な伝送路で、ＸＰＭを含めて非線形波形歪みを補正した場合に、信号入射パワーに対するＱ値の結果を表す図である。 図１４は、現実的な伝送路において様々な方法で非線形波形歪み補正を適用した場合としない場合に、伝送スパン数に対するＱ値の結果を表す図である。 図１５は、実験に用いた伝送系を示す模式図である。 図１６Ａは、パラメータ更新回数に対する、２次の群速度分散パラメータＤ２ （ｊ）の変化を表す図である。 図１６Ｂは、パラメータ更新回数に対する、非線形係数パラメータｇ（ｊ）の変化を表す図である。 図１６Ｃは、パラメータ更新回数に対する、偏波間位相変調パラメータδ（ｊ）の変化を表す図である。 図１６Ｄは、パラメータ更新回数に対する、ウォークオフパラメータｄ－５ （ｊ）の変化を表す図である。 図１６Ｅは、パラメータ更新回数に対する、ウォークオフパラメータｄ－１ （ｊ）の変化を表す図である。 図１６Ｆは、パラメータ更新回数に対する、ウォークオフパラメータｄｎ （ｊ）の変化を表す図である。 図１６Ｇは、パラメータ更新回数に対する、Mean Square Error（ＭＳＥ）とその移動平均の変化を表す図である。 図１７は、６スパン伝送路における光ファイバ接続状況と、逆伝搬計算順序を示す模式図である。 図１８Ａは、６スパン伝送路及び８スパン伝送路においてＳＰＭのみを考慮して非線形波形歪み補正を実施した場合における入射パワーに対する信号品質を表す図である。 図１８Ｂは、１０スパン伝送路及び１２スパン伝送路においてＳＰＭのみを考慮して非線形波形歪み補正を実施した場合における入射パワーに対する信号品質を表す図である。 図１８Ｃは、１４スパン伝送路及び１６スパン伝送路においてＳＰＭのみを考慮して非線形波形歪み補正を実施した場合における入射パワーに対する信号品質を表す図である。 図１９Ａは、６スパン伝送路及び８スパン伝送路においてＳＰＭ及びＸＰＭを考慮して非線形波形歪み補正を実施した場合における入射パワーに対する信号品質を表す図である。 図１９Ｂは、１０スパン伝送路及び１２スパン伝送路においてＳＰＭ及びＸＰＭを考慮して非線形波形歪み補正を実施した場合における入射パワーに対する信号品質を表す図である。 図１９Ｃは、１４スパン伝送路及び１６スパン伝送路においてＳＰＭ及びＸＰＭを考慮して非線形波形歪み補正を実施した場合における入射パワーに対する信号品質を表す図である。 図２０は、伝送スパンに対する信号品質の変化を表す図である。 図２１は、５チャネル、１１チャネル及び２１チャネルの場合における、データ長Ｎに対する、積演算回数を示す図である。
［本発明の実施の形態における基本的な考え方］
　周波数間隔がΔωであるＷＤＭ信号を対象に、本実施の形態に係る非線形波形歪みの補正技術を説明する。図１にＷＤＭ信号のスペクトルの模式図を示す。ここで、補正対象の信号をチャネル番号０とし、この信号の中心周波数を基準値ω０＝０とし、その他のチャネル番号ｎの信号の中心周波数をωｎ＝ｎΔω （ｎ＝±１，±２，．．．）とする。このとき、式（１）の包絡線振幅は、チャネル毎の包絡線振幅の和として以下のように表される。
Figure JPOXMLDOC01-appb-M000003
 　但し、Ａｐ，ｎ（ｔ）は、チャネルｎの偏波成分ｐに関するベースバンド（中心周波数が０を意味する）の包絡線振幅である。
　また、式（１）の非線形シュレディンガー方程式の線形項のみを書き出すと以下のようになる。
Figure JPOXMLDOC01-appb-M000004
 　光ファイバの長手方向の座標ｚ＝０における波形Ａｐ，ｎ（ｚ＝０，ｔ）を入力とし、距離ｈ伝搬後の波形についての式（４）の解は、周波数領域でチャネル毎に分離することができて、以下のようになる。
Figure JPOXMLDOC01-appb-M000005
　ここで、Ｄ２（ｈ）＝－β２ｈ／２及びＤ３＝－β３ｈ／６は、それぞれ２次及び３次の群速度分散の累積値である。また、Ｔｎ（ｈ、ωｎ）＝（β２ωｎ＋β３ωｎ ２／２）ｈは、チャネル番号ｎの信号に生ずる群遅延（ウォークオフ）を表している。式（５）の波形を、フーリエ変換を表す演算子Ｆを用いて時間領域で表すと、以下のようになる。
Figure JPOXMLDOC01-appb-M000006
　一方で、非線形シュレディンガー方程式(1)の非線形項のみを書き出すと以下のようになる。
Figure JPOXMLDOC01-appb-M000007
 　ここで、右辺の伝搬損失の項は本来線形項であるが、信号の強度が伝搬損失によって減衰することに伴う非線形性の変化を考慮するために含められている。包絡線振幅ＢｐをＡｐ＝Ｂｐ（ｚ）ｅｘｐ（－αｚ／２）と定義すると、Ｂｐは伝搬損失による減衰量を分離した振幅であり、式（７）はＢｐに対して次のように表される。
Figure JPOXMLDOC01-appb-M000008
 　ここで、γ（ｚ）＝γ０ｅｘｐ（－αｚ）である。
　式（８）は、伝搬損失によって信号の振幅が減衰する際に、非線形シュレディンガー方程式における非線形効果は、元の非線形係数γ０が長手方向に減衰していく形で記述できることを表している。式（８）をチャネル毎に分解すると、チャネル番号０の振幅については以下のように表される。
Figure JPOXMLDOC01-appb-M000009
 　式（９）の右辺において、チャネル番号０の信号強度の時間波形である｜Ｂｐ，０｜（ｐ＝１，２）を含む項はＳＰＭ、そしてチャネル番号ｎ≠０の信号強度の時間波形である｜Ｂｐ，ｎ｜（ｐ＝１，２）を含む項はＸＰＭを生じるものである。
　非特許文献４で提案されているように、式（９）において、伝搬時の波形変化に関する仮定を導入することで、近似解を得られる。第１の仮定は、チャネル番号０の信号の時間波形強度は距離ｚに対して不変であるとして、｜Ｂｐ，０（ｚ，ｔ）｜２＝｜Ｂｐ，０（０，ｔ）｜２とおくことであり、この仮定はスプリットステップフーリエ法の計算において通常用いられているものである。
　第２の仮定は、チャネル番号ｎ≠０の信号強度の時間波形は距離ｚに対して形状が不変である一方、ウォークオフによる群遅延の発生を考慮して｜Ｂｐ，ｎ（ｚ，ｔ）｜２＝｜Ｂｐ，ｎ（０，ｔ－ｄｎｚ）｜２とおくことである。この式は、強度の形状が伝搬距離に依存せず初期波形のまま保たれる一方で、伝搬距離に比例した遅延ｄｎｚが時間軸上で発生することを表している。ここで、以下で表すｄｎは、ウォークオフを表すパラメータであり、群速度の逆数に相当する。
Figure JPOXMLDOC01-appb-M000010
 　この仮定を導入した結果、周波数領域において式（９）をｚで積分することができて、距離ｈ伝搬後の波形について以下の解が得られる。
Figure JPOXMLDOC01-appb-M000011
Figure JPOXMLDOC01-appb-M000012
Figure JPOXMLDOC01-appb-M000013
　但し、式（１１）乃至（１３）を導出するために距離ｚに対して長さｈの積分を行う際、積分区間を［－ｈ／２，ｈ／２］とした。また、これらの式で用いられる記号は、以下のように定義される。
Figure JPOXMLDOC01-appb-M000014
　式（６）によりスプリットステップフーリエ法の線形ステップが記述され、式（１１）によってチャネル番号０の信号波形に対して、ＳＰＭとＸＰＭの両方の効果を考慮する非線形ステップが記述された。以下では、これらの式に基づいて、ある光ファイバ伝送路で受信した光信号波形を逆伝搬させる計算を行い、送信波形を推定する方法を述べる。
　逆伝搬計算の手順を説明するために、図２に、スパン数が２の伝送路で、スパンあたりのステップ数を２とした場合のステップの区割りと、線形ステップ及び非線形ステップの計算順序を表す模式図を示す。なお、異なるスパン数やステップ数の場合でも、図２を容易に一般化して考えることが可能である。図２において、左から右の向きを順伝搬方向とし、始点から終点にかけてスパン１及び２を定義する。逆に、右から左の向きを逆伝搬方向とし、各スパン内で逆伝搬方向にステップ１及び２を定義する。スパン内で各ステップの長さの比は等分配ではなく、ファイバの損失による信号パワーの減衰を考慮した非線形係数γ（ｚ）＝γ０ｅｘｐ（－αｚ）の積分値が等しくなるように設定する。例えば、ステップ数を２とした場合に、スパンの始点をｚ＝０、終点をｚ＝ｚｓ、ステップの区切りをｚ＝ｚ１とおいた場合に、以下の式が成立するようにｚ１を定める。
Figure JPOXMLDOC01-appb-M000015
　なお、スパン内で各ステップの長さの比を等分配した場合でも、本実施の形態はそのまま適用可能である。さらに言えば、スパン内でステップを区割りするやり方の違いは、学習を行わない非特許文献２及び４に示された方法では波形歪み補正の結果に影響を与えるが、パラメータの学習を実施する非特許文献５乃至８及び本実施の形態においては、学習の初期値設定に影響を与えるにとどまり、学習が収束して最適化された後は、差が出なくなる。また、スパンあたりのステップ数を１とした場合は、１スパン全体を１ステップとして計算する。
　次に、各ステップの計算を行うにあたっては、「対称型」のスプリットステップフーリエ法（非特許文献３を参照のこと）に基づく計算を行う。すなわち、距離がｈである一つのステップを前半と後半に等分割し、まずステップ前半の距離ｈ／２に対して、線形ステップとして式（６）を用いた計算を行う。
　次に、線形ステップの出力として得られた波形を入力とし、ステップ全体の距離ｈに対して、非線形ステップとして式（１１）の計算を行う。最後に、非線形ステップの出力として得られた波形を入力として、ステップ後半の距離ｈ／２に対して、式（６）の線形ステップの計算を実施する。この操作をステップ毎に繰り返すが、あるステップの後半の線形ステップは、次のステップの前半の線形ステップと合わせて計算を実施する。
　図２の例では、スパン＃２のステップ＃１を２分割し、入力波形に対して前半の線形ステップＬ（０）の計算を実施した後、非線形ステップＮ（１）の計算を行い、続いてステップ＃１後半の線形ステップとスパン＃２のステップ＃２前半の線形ステップを結合したＬ（１）の計算を行い、以下同様の計算を行う。スパン数がＮｓｐａｎで、スパンあたりのステップ数がＮｓｔｅｐの伝送路では、ｍ＝Ｎｓｐａｎ×Ｎｓｔｅｐとおくと、合計で（ｍ＋１）回の線形ステップと、ｍ回の非線形ステップが存在することになる。図２の例では、Ｎｓｐａｎ＝２及びＮｓｔｅｐ＝２であるから、ｍ＝４であり、線形ステップはＬ（０）からＬ（４）までのｍ＋１＝５回、非線形ステップはＮ（１）からＮ（４）までのｍ＝４回となる。
　式（６）の線形ステップでは、各チャネルの信号波形が独立して発展するが、式（１１）乃至（１３）に示した非線形ステップでは、１つのチャネルの信号波形の発展を計算する際に、ＸＰＭを考慮した計算を行うために、他のチャネルの信号波形を取り入れている。図３に、この様子を表した逆伝搬計算手順を示す。上述のとおり、線形ステップと非線形ステップを交互に計算することになる。逆伝搬計算において、チャネル番号ｎ、偏波成分ｐの入力波形をｘｐ，ｎとし、出力波形をｙｐ，ｎとする。ただし、ｘｐ，ｎは受信機で受信した伝送後の波形であり、またｙｐ，ｎは送信波形の推定結果に相当する。
　一般に、ｊ番目の線形ステップＬ（ｊ）の入出力波形をそれぞれｚｐ，ｎ （ｊ）及びｙｐ，ｎ （ｊ）とすると、非線形ステップＮ（ｊ）の入出力波形はｙｐ，ｎ （ｊ－１）及びｚｐ，ｎ （ｊ）である。また、図３に示すように、ｚｐ，ｎ （０）＝ｘｐ，ｎ及びｙｐ，ｎ （ｍ）＝ｙｐ，ｎである。Ｌ（ｊ）の出力波形ｙｐ，ｎ （ｊ）がＮ（ｊ＋１）に入力される際に、実線の矢印で示された経路は、式（１２）で示されるチャネル内のＳＰＭによる位相シフトを計算するための波形データの流れであり、点線の矢印で示された経路は、式（１３）で示されるチャネル間のＸＰＭによる位相シフトを計算するための波形データの流れである。
　次に、逆伝搬計算で用いるパラメータを勾配降下法によって最適化し、非線形波形歪み補正の性能を最大化する、本実施の形態における手法を説明する。以下では具体的な学習方法について示す。誤差関数Ｊ(θ)を次のように定義する。
Figure JPOXMLDOC01-appb-M000016
 　ここで、θは逆伝搬計算で用いるパラメータからなるベクトルであり、ｅｔ(θ)＝ｙｔ(θ)ーｄｔは時刻ｔにおける誤差である。ｙｔ(θ)は逆伝搬計算後の信号波形の時刻ｔにおける値であって、パラメータθの関数であるとみなし、ｄｔは時刻ｔにおける所望信号の値である。ここでは所望信号として、送信信号の波形を用いる。
　ここで、式（１６）の誤差関数の形は、平均二乗誤差（ＭＳＥ：Mean Square Error）と呼ばれる。送信信号の波形ｄｔと、パラメータθの逆伝搬伝送路を計算後に得られた波形ｙｔ(θ)との複数の組でデータセットを構成し、繰り返し計算に基づく確率的勾配降下法によって、誤差関数Ｊ(θ)を最小化するようにパラメータθを最適化する。非特許文献１に記載されているように、確率的勾配降下法に基づくパラメータθの更新式は次のように得られる。
Figure JPOXMLDOC01-appb-M000017
 　ただし、θｉはｉ回目の更新結果として得られたパラメータベクトルであり、ηは学習の速度を決定する微少な正の数、そして∇Ｊはパラメータに対する誤差関数の勾配を表し、次式のように計算される。
Figure JPOXMLDOC01-appb-M000018
　ここで、所望信号、すなわち送信信号はパラメータθに依存せず、
Figure JPOXMLDOC01-appb-M000019
 　であるから、結局パラメータの更新式は以下のように得られる。
Figure JPOXMLDOC01-appb-M000020
　式（２０）を用いて、θの要素である逆伝搬計算中の全てのパラメータについての更新を行う際、出力信号ｙｔに対して、それぞれのパラメータに対する勾配である∂ｙｔ／∂θを用いる。これを計算するための計算式を、微分連鎖規則により導出する。
　式（６）より、線形ステップＬｐ，ｎ （ｊ）における入出力波形ｚｐ，ｎ （ｊ）及びｙｐ，ｎ （ｊ）の関係は、以下のようになる。
Figure JPOXMLDOC01-appb-M000021
　式（２１）を直接微分することで、以下の式が得られる。
Figure JPOXMLDOC01-appb-M000022
Figure JPOXMLDOC01-appb-M000023
　これにより、出力波形ｙｐ，ｎ （ｊ）のＤ２ （ｊ）及びＤ３ （ｊ）による勾配は、入力波形ｚｐ，ｎ （ｊ）を用いて計算できる。ここで得られた勾配は、微分連鎖規則によって次のステップに送られ、最終的には、出力波形ｙｐ，ｎ＝ｙｐ，ｎ （ｍ）に対する勾配の形に改められる。Ｌｐ，ｎ （ｊ）より前のステップに含まれる任意のパラメータεに対するｙｐ，ｎ （ｊ）の勾配は、以下のように得られる。
Figure JPOXMLDOC01-appb-M000024
　これにより、出力波形ｙｐ，ｎ （ｊ）のパラメータεに対する勾配は、直前の非線形ステップより出力される勾配∂ｚｐ，ｎ （ｊ）／∂εを用いて計算できる。ただし、ウォークオフの値Ｔｎ （ｊ）は、線形ステップＬｐ，ｎ （ｊ）の前後の非線形ステップＮｐ，ｎ （ｊ）及びＮｐ，ｎ （ｊ＋１）に含まれるウォークオフパラメータｄｎ （ｊ）及びｄｎ （ｊ＋１）を用いて、以下のようにする。
Figure JPOXMLDOC01-appb-M000025
 　ここで、ｈｊは非線形ステップＮｐ，ｎ （ｊ）の区間幅である。
　次に、チャネル０の信号波形についての非線形ステップＮｐ，ｎ （ｊ）における入出力波形ｙｐ，ｎ （ｊ）及びｚｐ，ｎ （ｊ）の関係は、式（１１）より、以下のようになる。
Figure JPOXMLDOC01-appb-M000026
Figure JPOXMLDOC01-appb-M000027
Figure JPOXMLDOC01-appb-M000028
 　但し、Ｐｐ，ｎ （ｊ－１）＝｜ｙｐ，ｎ （ｊ－１）｜２は、入力波形ｙｐ，０ （ｊ－１）の強度であり、また、以下のような関係もある。
Figure JPOXMLDOC01-appb-M000029
　これから、非線形ステップＮｐ，ｎ （ｊ）で用いられるパラメータｇ（ｊ），δ（ｊ），α０ （ｊ），αｎ （ｊ）及びｄｎ （ｊ）に対する出力波形ｚｐ，ｎ （ｊ）の勾配は、式（２６）を直接微分することで以下のように得られる。
Figure JPOXMLDOC01-appb-M000030
Figure JPOXMLDOC01-appb-M000031
Figure JPOXMLDOC01-appb-M000032
Figure JPOXMLDOC01-appb-M000033
Figure JPOXMLDOC01-appb-M000034
 　このように、出力波形ｚｐ，ｎ （ｊ）と入力波形の強度Ｐｐ，ｎ （ｊ－１）とその周波数波形Ｐ~ｐ，ｎ （ｊ－１）（Ｐ~はＰの上に～）を用いて計算できる。
　また、Ｎｐ，ｎ （ｊ）より前のステップに含まれる任意のパラメータεに対するｚｐ，ｎ （ｊ）の勾配は、式（２６）をεで微分することで以下の式が得られる。
Figure JPOXMLDOC01-appb-M000035
Figure JPOXMLDOC01-appb-M000036
　なお式（２６）乃至（３６）では、チャネル番号０の信号についての非線形ステップにおける計算式を記載したが、一般のチャネル番号の信号についても同様に記述することが可能である。
　まとめると、線形ステップでは式（２２）及び（２３）、そして非線形ステップでは式（３０）乃至（３４）を用いて、そのステップで用いられるパラメータによる、そのステップの出力波形の勾配をそれぞれ計算し、その結果を次のステップに渡す。また線形ステップでは式（２４）、非線形ステップでは式（３５）及び（３６）を用いて、手前のステップより渡された勾配を、そのステップの出力についての勾配に更新して、さらに次のステップに渡す。入力側から出力側に向けてこのような計算を続けていくことで、逆伝搬計算で用いられるすべてのパラメータによる、最終的な出力波形ｙｐ，ｎ＝ｙｐ，ｎ （ｍ）の勾配を計算することができ、式（２０）によってパラメータの更新が可能となる。
　なお、図３において、チャネル番号ｎの各直交偏波成分についての線形ステップＬ１，ｎ （ｊ）及びＬ２，ｎ （ｊ）では、分散パラメータＤ２ （ｊ）及びＤ３ （ｊ）を共有し、学習の際に勾配を計算する際も、両偏波成分の勾配を一括して平均をとる。また非線形ステップＮ１，ｎ （ｊ）及びＮ２，ｎ （ｊ）でも同様に、パラメータｇ（ｊ），δ（ｊ），α０ （ｊ），αｎ （ｊ），ｄｎ （ｊ）を共通して利用するものとする。
　また、以下で述べる実施例では、チャネル番号ｎ＝０以外の信号については、非線形ステップの計算を実施せず、線形ステップにおける計算のみを実施する。非線形ステップにおける計算を無視した場合、ＸＰＭはおろかＳＰＭも補正されないことになるが、チャネル番号ｎ＝０に対する、それ以外のチャネルによるＸＰＭの影響を計算する上では、それ以外のチャネルの波形に生じる非線形波形歪みは無視することができることを確認済みである。
　以上に記した方法により、ＸＰＭを含めた非線形波形歪みを補正するための伝送路パラメータを勾配降下法によって最適化し、補正効果を最大化することができる。
　以下では、最適化するパラメータを取捨選択したうえで、さらに近似を適用することで計算量を削減する手法を示す。
　式（２３）では伝送路の３次の群速度分散に対する勾配を計算する方法を示しているが、シンボルレートが数十Ｇbaud以下の信号に対しては、３次の群速度分散が単一チャネルの波形に及ぼす影響は無視できるほど小さいので、ある初期値に設定した後は学習を行わずに固定しておくか、あるいは３次の群速度分散効果そのものを無視することもできる。ここでは３次の群速度分散について、初期値を設定した後は学習を行わずに固定するものとする。ただし３次の群速度分散に起因して、チャネル間のウォークオフが、チャネル間の周波数差に対して２次で変化することになるが、この効果はウォークオフＴｎ （ｊ）や、ウォークオフパラメータｄｎ （ｊ）の初期値を設定する際に考慮することにする。
　式（３２）及び（３３）では、各チャネルの損失係数に関する勾配を計算する方法を示しているが、損失係数は容易に測定できるパラメータであるため、勾配を計算することなく初期値から固定したままとする。
　次に、勾配の伝播の仕方に制限を加える。図４Ａ及び図４Ｂは、チャネル番号がｎ＝０及び１である、２チャネルの偏波多重信号の入力波形ｘｐ，ｎ（ｐ＝１，２；ｎ＝０，１）に対して、ｍ＝２の場合における線形ステップ及び非線形ステップの経路を示すものである。ここで、線形ステップＬ１，０ （０）で用いられる２次の群速度分散Ｄ２ （０）についての勾配が後段のステップに波及する経路について考える。
　線形ステップＬ１，０ （０）の出力波形ｙ１，０ （０）は、非線形ステップＮｐ，ｎ （１）（ｐ＝１，２；ｎ＝０，１）に送られる。それぞれの非線形ステップでは、式（２６）乃至（２８）より出力波形ｚｐ，ｎ （１）が計算されるとともに、式（３５）及び（３６）よりε＝Ｄ２ （０）の場合における勾配∂ｚｐ，ｎ （１）／∂εが計算されて、さらに後のステップに送られる。Ｄ２ （０）に対する勾配は、線形ステップＬｐ，ｎ （１）（ｐ＝１，２；ｎ＝０，１）を経て、非線形ステップＮ１，０ （２）に伝播する。
　図４Ａの実線は、Ｄ２ （０）に対する勾配の伝播経路を示すものである。しかしながら、Ｌｐ，０ （０）（ｐ＝１，２）で共通して用いられるＤ２ （０）の値が変化した場合に、ｙｐ，０ （０）がＸＰＭを通じてチャネル番号ｎ＝１の波形に影響し、さらにＸＰＭを通じて非線形ステップＮｐ，０ （２）でチャネル番号ｎ＝０の信号波形ｚｐ，０ （２）に及ぶ影響は、極めて小さいものであって、無視することが可能である。このように考えると、あるパラメータεによる勾配の伝播経路をチャネル内に限定することができ、図４Ｂに示した実線のような経路となる。このとき、式（３６）においてＨｎ （ｊ）(ω)を含む項を無視することになり、結果として次式のような形となる。
Figure JPOXMLDOC01-appb-M000037
　以上に示した近似による手法をまとめると、線形ステップでは、波形の発展として式（２１）、そのステップで用いられる２次の群速度分散に関する勾配の計算式として式（２２）、そして前のステップより伝播してきた任意のパラメータεに関する勾配の更新式として式（２４）を用いて計算できる。また非線形ステップでは、波形の発展として式（２６）乃至（２８）、そのステップで用いられるパラメータｇ（ｊ），δ（ｊ），ｄｎ （ｊ）に関する勾配の計算式として式（３０）、（３１）及び（３４）、そして前のステップより送られてきた任意のパラメータεに関する勾配の更新式として、式（３５）及び（３７）を用いることができる。
　図５Ａに、パラメータ最適化処理において、各伝送路パラメータによる波形の勾配を計算する処理のフローを示す。線形ステップＬ（ｊ）においては、式（２１）でｚｐ，ｎ （ｊ）からｙｐ，ｎ （ｊ）への波形の発展がなされるが、式（２２）及び（２３）による演算によって、Ｌ（ｊ）で用いられる各伝送パラメータ（Ｄ２ （ｊ）及びＤ３ （ｊ））に関する波形の勾配（∂ｙｐ，ｎ （ｊ）／∂Ｄ２ （ｊ）及び∂ｙｐ，ｎ （ｊ）／∂Ｄ３ （ｊ））が計算される。また式（２４）による演算により、Ｌ（ｊ）より前のステップに含まれる全ての伝送パラメータεに関する波形の勾配が計算され、∂ｚｐ，ｎ （ｊ）／∂εが∂ｙｐ，ｎ （ｊ）／∂εに更新される。∂ｙｐ，ｎ （ｊ）／∂Ｄ２ （ｊ）及び∂ｙｐ，ｎ （ｊ）／∂Ｄ３ （ｊ）は、∂ｙｐ，ｎ （ｊ）／∂εに含められ、まとめて次のステップであるＮ（ｊ＋１）に送られる。
　次に、非線形ステップＮ（ｊ＋１）においては、式（２６）乃至（２８）でｙｐ，ｎ （ｊ）からｚｐ，ｎ （ｊ＋１）への波形の発展がなされるが、式（３０）乃至（３４）による演算によって、Ｎ（ｊ＋１）で用いられる各伝送パラメータ（ｇ（ｊ＋１），δ（ｊ＋１），α０ （ｊ＋１），αｎ （ｊ＋１），ｄｎ （ｊ＋１））に関する、波形の勾配（∂ｚｐ，ｎ （ｊ＋１）／∂ｇ（ｊ＋１），∂ｚｐ，ｎ （ｊ＋１）／∂δ（ｊ＋１），∂ｚｐ，ｎ （ｊ＋１）／∂α０ （ｊ＋１），∂ｚｐ，ｎ （ｊ＋１）／∂αｎ （ｊ＋１），∂ｚｐ，ｎ （ｊ＋１）／∂ｄｎ （ｊ＋１））が計算される。また、式（３５）及び（３７）による演算により、Ｎ（ｊ＋１）より前のステップに含まれる全ての伝送パラメータεに関する波形の勾配が計算され、∂ｙｐ，ｎ （ｊ）／∂εが∂ｚｐ，ｎ （ｊ＋１）／∂εに更新される。なお、式（３７）の代わりに式（３６）を用いても良い。∂ｚｐ，ｎ （ｊ＋１）／∂ｇ（ｊ＋１）、∂ｚｐ，ｎ （ｊ＋１）／∂δ（ｊ＋１）、∂ｚｐ，ｎ （ｊ＋１）／∂α０ （ｊ＋１）、∂ｚｐ，ｎ （ｊ＋１）／∂αｎ （ｊ＋１）、∂ｚｐ，ｎ （ｊ＋１）／∂ｄｎ （ｊ＋１）は、∂ｚｐ，ｎ （ｊ＋１）／∂εに含められ、まとめて次のステップであるＬ（ｊ＋１）に送られる。
　このような計算をステップ毎に繰り返すことで、伝送路中のすべてのパラメータθに関する、最終的な出力ｙｔの勾配である、∂ｙｔ／∂θが計算されるのである。
　さらに、図５Ｂに、パラメータ最適化処理において、各伝送路パラメータを更新する処理のフローを示す。すなわち、出力波形ｙｔと各伝送パラメータθｉによる勾配∂ｙｔ／∂θｉと所望信号ｄｔとを用いて、ｉ番目のステップにおける伝送路パラメータθｉを、式（２０）を用いて更新するものである。図５Ｂに示すように、線形ステップにおいては、伝送路パラメータＤ２及びＤ３が更新され、非線形ステップにおいては、伝送路パラメータθとしてｇ、δ、α０、αｎ、ｄｎが更新される。
　なお、光ファイバの損失係数α０及びαｎの値や、偏波間相互位相変調に対応する係数δの値が１であることが判明している場合は、これらの伝送路パラメータについての学習を省略しても良い。また、チャネル内の３次分散効果が無視できる場合には、Ｄ３の学習を省略しても良い。
［実施の形態に係るシステム構成］
　図５Ｃは、本実施の形態に係る光伝送システムの模式図を示している。送信機によって生成されたＷＤＭ信号は、光ファイバと光増幅器から構成される伝送路を伝搬し、受信機に到達する。受信機では、アレイ導波路格子などの分波装置を用いて、ＷＤＭ信号をチャネル毎に分波する。分波された各チャネルの光信号波形は、コヒーレント受信機によって電気信号波形に変換された後、信号処理のためのデジタルシグナルプロセッサ（ＤＳＰ：Digital Signal Processor）に入力される。ＤＳＰでは、電気信号波形がアナログデジタル（ＡＤ：Analog-to-Digital）変換によって数値データとなり、最終的に受信ビット列に変換して出力するための各種復調処理のための計算が実施される。復調処理としては、時間的なサンプリングのタイミング制御、サンプリングレートの変換、クロック同期、フィルタリング、偏波回転、搬送波再生、適応等化器などの線形波形歪み補正、本実施の形態に係る非線形波形歪み補正、シンボル判定、誤り訂正処理などが含まれる。なお、ＤＳＰには、上で述べたパラメータ最適化処理を実行する機能を搭載することが可能であるが、パラメータ最適化処理をＤＳＰでは実施せずに外部の計算機で実施し、最適化されたパラメータをダウンロードして本実施の形態に係る非線形波形歪み補正に適用することも可能である。ＡＤ変換機能については、ＤＳＰから分離して、ＤＳＰとコヒーレント受信機の間に設置してもよい。従来の光伝送システムでは、ＤＳＰはチャネル毎に独立して設けられており、ＤＳＰ同士が連係して動作することは基本的にないが、本実施の形態では、図５Ｃのように、単一のプロセッサで複数チャネルの波形データを一括して扱うか、独立した複数のプロセッサが相互にデータをやり取りすることで、上で述べたＸＰＭによる位相シフトの計算を行う。
　図５Ｄに、ＤＳＰによって実行され、且つ本実施の形態に関係する処理フローを示す。まず、波形形状が既知である信号を伝送し、受信機でＤＳＰが復調した受信信号波形をもとに、ＤＳＰ内もしくは外部の計算機において伝送路パラメータを最適化するパラメータ最適化処理を実行する（工程Ｓ１）。これは、図５Ａ及び図５Ｂに示したようなフローに従った演算を含む。パラメータ最適化処理は、ＤＳＰ内に構成されるパラメータ最適化部もしくは外部の計算機により実行される。次に、ＤＳＰ内に構成される波形歪み補正部は、最適化された伝送路パラメータを用いた逆伝搬処理を実行する（工程Ｓ３）。この逆伝搬処理では、図２及び図３で述べたように、波長分割多重伝送時における複数のチャネルの各々について、線形ステップと非線形ステップのそれぞれが最適化された伝送路パラメータを持ち、これらの計算を順番に行うことで、当該チャネル内で生ずる自己位相変調に起因する波形歪みと、当該チャネル以外のチャネルとの間で生じる相互位相変調に起因する波形歪みとを補正し、送信時の波形を推定する。この逆伝搬処理は、ＤＳＰにより実行される。このような処理を行うことで、計算量を抑えつつ、精度を向上させることが出来るようになる。
［実施例］
　数値計算を用いた光伝送シミュレーションに基づく具体的な事例を基に、本実施の形態において示した手法による、ＷＤＭ光信号の非線形波形歪み補正の効果について説明する。
　検討対象とする光信号は、シンボルレートが３２Ｇbaudの偏波多重（ＤＰ：Dual-Polarization）６４値直交振幅変調（ＱＡＭ：Quadrature Amplitude Modulation）信号を、周波数間隔５０ＧＨｚで９チャネル波長分割多重した信号であり、これにランダムな雑音を付与して、ＳＮ比を２５ｄＢに設定する。なお送信信号のスペクトルは、ロールオフ係数が０．０５のルートナイキストフィルタが適用されたものとする。９チャネルのＷＤＭ信号について、周波数が低い側からチャネル番号を－４から＋４まで割り振り、中心のチャネル番号０の信号品質に注目し、非線形波形歪み補正の動作を検証する。
　１スパンが長さ８０ｋｍの標準シングルモード光ファイバ（ＳＳＭＦ：Standard Single-Mode Fiber）と、その伝搬損失を増幅する光増幅器で構成され、４スパンから１０スパンまでの伝送路を考え、これらの伝送路を用いて光信号を伝送する計算を行う。理想的な伝送路では、全スパンのＳＳＭＦが同じパラメータを持つものとして、２次および３次の群速度分散値をそれぞれ１６．６４１ｐｓ／ｎｍ／ｋｍおよび０．０６ｐｓ／ｎｍ２／ｋｍとし、非線形係数を１．３Ｗ－１ｋｍ－１、そして伝搬損失係数を０．１９２ｄＢ／ｋｍとおく。一方、現実の伝送路では、様々なパラメータがスパンごと、あるいは地点ごとに異なっており、また光信号のパワーも理想的な状態から揺らいでおり、結果的に群速度分散の効果や非線形効果の大きさがばらつくことになる。この状況を考慮するために、理想的な伝送路とは別に、図６のように２次の群速度分散値と非線形係数をステップごとに変動させた伝送路を考え、以下ではこの伝送路を現実的な伝送路と呼ぶことにする。なお、図６において、灰色実線は、２次の群速度分散値（ａ）と非線形係数（ｂ）のそれぞれについて、１０スパンすべてのパラメータに対する平均値を示している。ただし、４乃至８スパンまでの短い伝送路を考える場合は、存在するスパン内での平均値を考えるものとする。
　光伝送シミュレーションを実施する数値計算では、時間軸でシンボルレートの３２倍にサンプリングした波形データを扱う。スパン数を４，６，８，１０のいずれかに設定し、各スパンで同じパラメータの理想的な伝送路又は図６に示したパラメータを有する伝送路において、１スパンあたりのステップ数を８００とするスプリットステップフーリエ法を用いて、式（１）の非線形シュレディンガー方程式に従って光信号が伝送路を伝搬する様子を計算し、出力波形データを得る。式（１）で偏波間相互位相変調に対応する係数δの値を１と設定し、その他のパラメータについては、伝送路を構成する光ファイバのパラメータを規格化して得られた値を用いる。また光増幅器の雑音指数を６ｄＢとし、利得と雑音指数によって定められる雑音電力を持つランダムなガウス分布雑音を、自然放出光雑音として増幅の際に信号に付加する。各スパンで伝搬開始時の光信号パワーを「スパン入射パワー」と呼び、いくつかのスパン入射パワーに対して伝送シミュレーションを実施する。スパン入射パワーが小さい場合は、顕著な非線形波形歪みは生じないものの、雑音によるＳＮ比劣化が顕著になり、信号品質が低下する。スパン入射パワーを大きくしていくと、ＳＮ比は向上していくものの、非線形波形歪みが顕著になっていき、あるスパン入射パワーにおいて信号品質がピーク値を取った後、さらに大きなスパン入射パワーに対しては信号品質が劣化する。
　検討する伝送スパン数ならびに入射パワーそれぞれの条件に対して、シンボル数が１６３８４の送信信号データをランダムビットパターンより生成し、伝送シミュレーションを実施し、出力波形を保存する。伝送路の入力波形と出力波形の組をデータセットとして、合計２００のデータセットを用い、上記数式を用いた勾配降下法によって、非線形波形歪みの補正が最適化されるように、逆伝搬計算のパラメータを学習する。本実施例に係る、パラメータの学習を伴う逆伝搬計算においては、１スパンあたりのステップ数を１に設定する。これに対して、学習を伴わない逆伝搬計算においては、１スパンあたりのステップ数が１の場合と２の場合の両方を考え、学習を行わない場合との効果の比較を行う。
　各パラメータの学習にあたっては勾配降下法を用いるが、パラメータ更新の具体的な実装方法として、非特許文献９で提案されているAdaBeliefという手法を用いる。なお、AdaBelief以外の手法を用いても学習は可能である。AdaBeliefの更新式は、μ０＝０，ν０＝０とおいて以下のように得られる。
Figure JPOXMLDOC01-appb-M000038
　但し、∇Ｊｉはｉ回目の更新時におけるパラメータθに対する誤差関数の勾配であり、またｂ１＝０．９，ｂ２＝０．９９９，ｅ＝１０－８は定数である。式（３８）のηは学習係数であり、パラメータ毎に適した値を設定する。本実施例においては、２次の分散値に対応するパラメータＤ２に対する学習係数をη＝１．０とおき、非線形係数に対応するパラメータｇに対してはη＝２．０×１０－６、偏波間相互位相変調係数に対応するパラメータδに対してはη＝２．０×１０－５、そしてウォークオフに対応するパラメータｄｎに対してはη＝２．０×１０－４を用い、式（３８）による更新を３万回繰り返して学習を行う。なおデータセットの個数は２００であるので、２００回の更新ごとにデータセットの並び順をランダムに入れ替え、繰り返し学習を行う。なお、この操作によって過学習が発生することはない。また、信号の変調フォーマットやチャネル数などの条件を変更して同様の学習を実施する際に、必要に応じてパラメータが収束するまで３万回以上の学習を実施することもあり得る。式（２０）のｙｔ及びｄｔとしては、データセットを構成する１６３８４シンボルのすべてを用いる。学習後の検証作業には、異なるビットパターンで、シンボル数が２６２１４４である別の波形を用い、非線形波形歪み補正後の信号品質を評価する。
　まず、パラメータが全スパンで同じ値を持つ理想的な伝送路において、非線形波形歪み補正の結果を示す。図７Ａ乃至図７Ｅは、６スパンの理想的な伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合のデータセットに対して、ＸＰＭも考慮した逆伝搬計算において、パラメータの更新を３万回繰り返した際の、各種数値の変動を示している。図７Ａでは、パラメータ更新の繰り返し回数に対して、線形ステップＬ（ｊ）（０≦ｊ≦６）における分散パラメータＤ２ （ｊ）の変動の様子を示している。ｊ＝０および６の線形ステップでは、距離が他のステップの半分であるため、Ｄ２ （ｊ）の初期値も半分である。学習を進めていくと、ｊ＝１，２，３，４，５のパラメータは同じ値に収束し、ｊ＝０のパラメータも近い値に収束する。これに対して、ｊ＝６の場合、すなわち最も送信機に近い線形ステップのみ、他とは異なる値に収束している。これらの結果は、いま検討の対象としている理想的な伝送路においてパラメータが全スパンで等しいという事実と必ずしも一致しないが、学習による性能の最適化が行われた結果として得られたものであり、後述するように学習を行わない場合よりも性能が向上している原因のひとつになっていると考えられる。
　図７Ｂでは、非線形ステップＮ（ｊ）（１≦ｊ≦６）の非線形パラメータｇ（ｊ）の変動の様子を示している。ｊ＝６の場合、すなわち最も送信機に近い非線形ステップのみ他とはやや異なる値に収束しているが、その他のステップではほぼ同じ値に収束している。図７Ｃ及び図７Ｄでは、それぞれウォークオフパラメータｄ－４ （ｊ）およびｄ－１ （ｊ）の変化の様子を示している。すべてのステップにおける値がほぼ重なっており、また学習によるパラメータの変化もほぼ見られない。理想的な伝送路では、学習開始時に各スパンに共通な正しいパラメータを初期値として与えているため、学習によってウォークオフを調整する必要がなかったことによるものと考えられる。図には記載していないが、その他のウォークオフパラメータもｄ－４ （ｊ）及びｄ－１ （ｊ）と同様に収束している。
　図７Ｅでは、受信信号波形と、所望信号、すなわち送信信号波形から計算されるＭＳＥの値をプロットしている。学習開始直後からＭＳＥが急激に減少し、繰り返し回数が１０００程度でほぼ収束し、その後は安定に推移していることがわかる。実際の伝送システムで本手法を用いる際には、ＭＳＥが最低値付近で振動するようになったときに、各種パラメータが収束していなくても更新を終了し、学習後の定常的な波形補正を開始してもよい。
　図８Ａ及び図８Ｂは、伝送スパン数が４，６，８，１０それぞれの場合で、逆伝搬計算による非線形波形歪み補正を行った場合と行わない場合について、各スパン入射パワーに対するＱ値の計算結果を示している。ただし逆伝搬計算の条件としては、１step/spanおよび２steps/spanでパラメータを学習しない逆伝搬計算を行った場合と、１step/spanで本実施の形態に係るパラメータの学習を行ってから逆伝搬計算を実施した場合を考えていて、いずれもＸＰＭによる非線形波形歪みを補正せずに、ＳＰＭによる歪みのみを補正したときの結果を示している。またＱ値は、ビット誤り率（ＢＥＲ：Bit Error Rate）の値から以下の式に従って変換したものである。
Figure JPOXMLDOC01-appb-M000039
　図８Ａ及び図８Ｂの結果は全体として、入射パワーが低い領域ではＳＮ比の劣化が原因でＱ値が劣化しており、逆に入射パワーが大きい領域では、非線形波形歪みによってＱ値が劣化している。学習しない場合の逆伝搬計算の結果に注目すると、１step/spanの場合は補正の効果がわずかしかなく、２steps/spanではそれなりの補正効果が得られている。また本実施の形態に係る方法でパラメータを学習した逆伝搬計算の結果では、スパンあたりのステップ数が１であるにも関わらず、学習しない場合の２steps/spanの結果をやや上回る信号品質が得られている。このことは、学習を実施することで、波形補正に真に必要なパラメータが獲得でき、本質的には２steps/spanの計算が不要である可能性を示唆している。実際に、結果は記載しないが、２steps/spanで本実施の形態に係る方法でパラメータの学習を行うと、ＭＳＥが下がりきった条件における非線形波形歪み補正の性能は、１step/spanの場合と変わらないものが得られる。
　非特許文献８で報告されている、物理現象特化型ニューラルネットワークによる非線形波形歪み補正の場合についても、１２スパンの伝送路に対して、１step/spanの結果が２steps/spanの結果とほぼ変わらない一方で、０．５step/spanの構成、すなわちステップ数が６の場合でも、性能は劣化しつつもそれなりの補正効果が得られたという結果が報告されている。このことから、本実施の形態に係る手法を用いた場合は、２steps/span未満のステップ数で十分な補正効果が得られ、さらに１step/span未満の構成でもそれなりの補正効果が得られるものと考えられる。
　図９Ａ及び図９Ｂは、図８Ａ及び図８Ｂと同様の結果であるが、ＳＰＭに加えて、ＸＰＭも考慮に入れた逆伝搬計算を実施した場合の結果を示している。ＸＰＭによる歪みの補正を行うことで信号品質が向上しており、その中でも、本実施の形態に係るパラメータの学習を実施した後の、１step/spanの逆伝搬計算によって波形補正された信号品質が、学習を行わない２steps/spanの逆伝搬計算の結果を上回り、この中では最良である。また非特許文献４の手法と比較すると、本実施の形態に係る手法は１step/spanで十分な効果が得られており、現実的な計算量で実現することが可能である。
　図１０は、図８Ａ及び図８Ｂ並びに図９Ａ及び図９Ｂに示した非線形波形歪み補正の結果から、いくつかの条件で非線形波形歪み補正をする場合と、補正を行わない場合について、それぞれ最適な入射パワーの下での、スパン数に対するＱ値の結果を示したものである。誤りなくデータを受信できるＱ値の閾値を７．５ｄＢに設定すると、補正しない場合の最大伝送距離は６スパンであるのに対して、本実施の形態に係る方式である、ＳＰＭとＸＰＭの両方を考慮してパラメータの学習を実施する手法を用いることで、１０スパン以上の伝送が可能となり、大幅に伝送距離を延伸可能であることがわかる。
　次に、図６にパラメータを示した現実的な伝送路についての結果を図１１Ａ乃至図１４に示す。図１１Ａ乃至図１１Ｅは、６スパンの現実的な伝送路で、スパン入射パワーが＋２ｄＢｍ／ｃｈの場合のデータセットに対して、ＸＰＭも考慮した逆伝搬計算において、パラメータの更新を３万回繰り返した際の、各種数値の変動を示している。ただし、学習を開始するにあたっての伝送路パラメータの初期値としては、図６に示した現実的な伝送路の真のパラメータが不明であるという状況を想定し、６スパンの平均値を与えている。図１１Ａでは２次の群速度分散パラメータＤ２ （ｊ）、図１１Ｂでは非線形パラメータｇ（ｊ）の変化の様子を示しており、それぞれスパンごとに異なる値に収束している。この結果は、元の伝送路のパラメータがスパン毎に変化している事実に定性的に一致するが、得られた値は非線形波形歪み補正を最大化する値であって、必ずしも図６で与えた伝送路パラメータに一致するわけではない。図１１Ｃ及び図１１Ｄでは、それぞれウォークオフパラメータｄ－４ （ｊ）及びｄ－１ （ｊ）が、各スパンの異なる分散値に対応して異なる値に収束していることがわかる。図には記載していないが、その他のウォークオフパラメータもｄ－４ （ｊ）及びｄ－１ （ｊ）と同様に収束している。図１１Ｅでは図７Ｅと同様に、学習の回数に対するＭＳＥの値を示しているが、図７Ｅに示した理想的な伝送路の場合と比較して、図１１Ｅの現実的な伝送路に対する結果では、繰り返し回数が約１００００で収束している。図７Ｅに示した結果と同様に、この場合もＭＳＥが最低値付近で振動するようになったときに、各種パラメータが収束していなくても更新を終了し、学習後の定常的な波形補正を開始してもよい。
　図１２Ａ及び図１２Ｂでは、図８Ａ及び図８Ｂと同様に、伝送スパン数が４，６，８，１０それぞれの場合で、逆伝搬計算による非線形波形歪み補正を行った場合と行わない場合について、各スパン入射パワーに対するＱ値の計算結果を示している。但し、逆伝搬計算の条件としては、１step/span及び２steps/spanでパラメータを学習しない逆伝搬計算を行った場合と、１step/spanで本実施の形態に係るパラメータの学習を行ってから逆伝搬計算を実施した場合で、いずれもＸＰＭによる非線形波形歪みを補正せずに、ＳＰＭによる歪みのみを補正したときの結果を示している。また学習しない逆伝搬計算については、計算に用いるパラメータとして、図６に示した真の値を与えた場合と、その平均値を与えた場合との両方の結果を示している。理想的な伝送路について得られた図８Ａ及び図８Ｂの結果と同じく、現実的な伝送路に対しても、本実施の形態に係るパラメータの学習を実施する逆伝搬計算によって波形補正された信号品質がこの中では最良であることがわかる。学習を行わない逆伝搬計算においては、伝送路パラメータの真の値を与えた場合と、平均値を与えた場合で結果に大きな差がなく、ＳＰＭのみを補正する状況においては、補正の性能は伝送路のパラメータに大きく依存していないことがわかる。
　図１３Ａ及び図１３Ｂは、ＳＰＭだけでなくＸＰＭも考慮した非線形波形歪み補正を実施した場合について、図１２Ａ及び図１２Ｂと同様の結果を示している。現実的な伝送路においても、ＸＰＭによる歪みの補正を行うことで信号品質が向上しており、その中でも、本実施の形態に係るパラメータの学習を実施する、１step/spanの逆伝搬計算によって波形補正された信号品質が、学習を行わず、真の伝送路パラメータを与えた２steps/spanの逆伝搬計算の結果を上回り、この中では最良である。図１３Ａ及び図１３Ｂの結果において注目すべき点は、学習を行わない逆伝搬計算において、伝送路のパラメータとして真の値の平均値を与えた場合の結果が、真の値を与えた場合の結果と比較して、著しく劣化していることである。これは、ＸＰＭを補正するにあたって、異なるチャネル間の時間的な位置関係を表すウォークオフの値を正しく与えなければ、ＸＰＭによる位相シフトが正しく計算できずに補正性能が劣化することが原因である。ゆえに、伝送路のパラメータがスパンごとに揺らいでおり、かつ正しい値が未知の状況においては、本実施の形態で示した手法による、伝送路パラメータを学習することでウォークオフを最適値に設定し、効果的にＸＰＭを補正できるという点が、極めて重要な意味を持つのである。
　図１４は図１０と同様に、図１２Ａ及び図１２Ｂ並びに図１３Ａ及び図１３Ｂに示した非線形波形歪み補正の結果から、いくつかの条件で非線形波形歪み補正をする場合と、補正を行わない場合について、それぞれ最適な入射パワーの下での、スパン数に対するＱ値の結果を示したものである。現実的な伝送路においても、理想的な伝送路の場合と同様の結果が得られており、本実施の形態に係るパラメータの学習を実施した後で逆伝搬計算を実施することで、伝送距離を大幅に延伸することが可能となる。
　上述のシミュレーション結果に加え、周回伝送実験によっても実施の形態の効果を検証した。図１５は周回伝送実験系を示しており、以下ではこの系を用いた実験手順の詳細と、実験結果を示す。
　送信機（Ｔｘ：Transmitter）の構成として、波長可変光源より出力される波長の異なる１１チャネルの連続光が、１６×１偏波保持カプラで合波されてからリチウムニオブ酸（ＬＮ）偏波多重ＩＱ変調器に入力され、変調器に印加される電気信号によって光信号波形に変調されて出力される。１１チャネルの連続光が持つ周波数は、中心周波数を１９３．１ＴＨｚ（波長が１５５２．５２４ｎｍ）として、１９２．８５ＴＨｚ乃至１９３．３５ＴＨｚまで５０ＧＨｚ間隔で設定されており、周波数の低い順からチャネル番号をｎ＝－５，－４，・・・，４，５と定める。変調器に印加される４チャネルの電気信号は、サンプリングレートが６４ＧSample／ｓの任意波形発生器より生成されたものであり、各チャネルは偏波多重ＩＱ変調信号のＸ偏波Ｉチャネル成分、Ｘ偏波Ｑチャネル成分、Ｙ偏波Ｉチャネル成分、そしてＹ偏波Ｑチャネル成分に相当する。この４チャネル電気信号は、それぞれドライバアンプによって増幅され、変調器に印加される。変調器から出力された光信号は、シンボルレートが３２Ｇbaudであり、ロールオフ係数が０．１のルートナイキスト波形を持つ、偏波多重ＱＡＭ信号である。変調フォーマットとしては、単一偏波単一シンボル当たりのビット数が４bitである均一分布１６ＱＡＭ信号と、単一偏波単一シンボル当たりのビット数が５bitである確率的分布（ＰＳ：Probabilistically Shaping）６４ＱＡＭ信号を用いる。１６ＱＡＭおよびＰＳ－６４ＱＡＭの各変調フォーマットについて、ランダムビットパターンより変調される、単一偏波あたり６５５３６シンボルで構成される信号波形を４パターンずつ生成する。４パターンのうちの一つを選択すると、Ｔｘからはそのパターンの波形が繰り返し送信される。なお変調によって得られる１１チャネルのＷＤＭ信号は、すべてのチャネルで同じ波形に変調されているが、群速度分散が存在する伝送路を長距離伝送することで、ウォークオフ（チャネル間の群遅延）が生じるため、長距離伝送後にランダムな波形間のＸＰＭが生じることになる。
　Ｔｘで生成された１１チャネルのＷＤＭ信号は、光増幅器によって光パワーが増幅され、可変光減衰器（ＶＯＡ：Variable Optical Attenuator）によって光パワーが調整された後、帯域通過フィルタ（ＢＰＦ：Band-Pass Filter）によって信号帯域外の光雑音が除去され、周回伝送動作を切り替えるスイッチである音響光学変調器（ＡＯＭ：Acoustic Optical Modulator）に入力される。周回伝送路に対する信号の入出力は、３ｄＢカプラを通じて行われる。周回伝送路は入力側から順に、光増幅器、ＢＰＦ、ＶＯＡ、長さが８４．１ｋｍのＳＳＭＦ、光増幅器、ＢＰＦ、ＶＯＡ、長さが８０．５ｋｍのＳＳＭＦ、光増幅器、アイソレータ（右向き矢印）、偏波スクランブラ（Pol. Scrambler）、そしてＡＯＭで構成され、一周あたり２スパンのＳＳＭＦを含む。実験では伝送距離を６，８，１０，１２，１４，１６スパンのいずれかに設定して、それぞれの距離に対して信号を伝送し、非線形波形歪み補正の効果を検証する。なお長さが８４．１ｋｍのＳＳＭＦと、８０．５ｋｍのＳＳＭＦはやや異なる分散特性を持っており、測定の結果、周波数１９３．１ＴＨｚにおける群速度分散値と分散スロープの値は、前者がそれぞれ１７．１４ｐｓ／ｎｍ／ｋｍおよび０．０６２ｐｓ／ｎｍ２／ｋｍ、後者がそれぞれ１６．５５ｐｓ／ｎｍ／ｋｍおよび０．０５８ｐｓ／ｎｍ２／ｋｍであり、周回伝送路一周の平均特性としては群速度分散値が１６．８５ｐｓ／ｎｍ／ｋｍ、分散スロープが０．０６０ｐｓ／ｎｍ２／ｋｍと見積もられた。
　周回伝送路より出力された信号は、通過帯域が５０ＧＨｚのＢＰＦによって１１チャネルのうちのいずれか１チャネルのみが切り出され、光増幅器で増幅された後、受信機（Ｒｘ：Receiver）に入力される。Ｒｘは、電気帯域が３３ＧＨｚで８０ＧSample／ｓの４チャネルリアルタイムオシロスコープと、局所発振光源（ＬＯ：Local Oscillator）および光フロントエンドからなるデジタルコヒーレント受信機であり、オシロスコープで取得したリアルタイム波形をオフラインデジタル信号処理によって復調した後、非線形波形歪み補正のための信号処理をやはりオフラインで実施する。ただし実験では、チャネル番号がｎ＝０である中心のチャネルの信号品質にのみ着目し、これを最大化するように、実施の形態に係る非線形波形歪み補正のための学習を実施して、補正後の信号品質を評価する。
　ＷＤＭ信号を受信した後に非線形波形歪み補正のための逆伝搬計算をオフラインで行うために、全１１チャネルの波形をチャネルごとに受信する。逆伝搬計算を実施せずに、信号の品質評価まで行う通常の復調処理としては、分散補償、マッチドフィルタとして送信時に適用したものと同じルートナイキストフィルタの適用、偏波回転および偏波多重成分の分離、２サンプル／シンボルへのリサンプリング、リタイミング、搬送波周波数推定および補償、搬送波位相再生および３タップのフィードフォワード型線形適応等化処理、そしてシンボル判定およびビットパターンの取得を行う。ここで、３タップの適応等化処理は、偏波多重信号に対応したバタフライ型２×２ＭＩＭＯ処理となっていて、伝送路の複屈折やＸＰＭなどに起因して発生する偏波クロストークを補償するうえで有効である。一方で、このように復調した信号波形は、ルートナイキストフィルタの適用によって、受信時の波形からは大きく変化しており、そのまま逆伝搬計算を適用することはできない。そこで、通常の復調処理の中で分散補償を適用後の波形に対して、ルートナイキストフィルタを適用せずに、以降の復調過程で実施された処理を同様に行って、受信直後の波形をできるだけ忠実に再現する。
　ところで、本来であれば、複数のトランシーバを用いて全チャネルの光信号を同時に受信し、すべてのチャネルの波形を使用して、それらを復調することなしに、本実施の形態に係る非線形波形歪み補正のための逆伝搬計算を行って、最後に復調を実施すべきである。しかしながら、ここで実施している実験では、１台のＲｘを用いて各チャネルを順に受信して復調していくため、測定された全チャネルの波形は、チャネル間の同期がとれたものではない。特に、長距離の光ファイバを伝送後のＷＤＭ信号には、群速度分散の効果によってチャネル内の線形波形歪みに加えて、チャネル間に群遅延（ウォークオフ）が生じており、逆伝搬計算を行う際には、受信時の正確な群遅延量を保ったまま計算を開始する必要があるが、非同期に受信した波形を復調せずにそのまま逆伝搬計算すると、その条件が満たされない。そこで、いったん全チャネルを独立に受信して、分散補償を含めた復調を実施し、得られた波形から既知のパイロットシンボルを検出し、全チャネルのタイミングを同期させてから、補償した群速度分散値を再度付与することで各チャネルにウォークオフを与え、一括受信時に得られるであろうＷＤＭ信号波形を再現し、逆伝搬計算を開始するという手順を採る。
　実験では、異なる伝送スパン数および入射パワーの各条件において、各変調フォーマットの信号をそれぞれ伝送、受信、且つ復調し、全チャネルに共通する送信波形と、上に述べた全チャネル一括受信波形を組み合わせて、データセットを生成する。各変調フォーマットについて、異なる４つの波形パターンに対応して４つのデータセットが生成される。なお逆伝搬計算を実施中に発生するウォークオフを正しく処理するために、受信時には送信波形の一周期を構成する６５５３６シンボルよりも十分に多くのシンボルが含まれるように測定する。チャネル番号がｎ＝０である中心の信号と、チャネル番号がｎ＝±５である端のチャネルは周波数差が２５０ＧＨｚであり、最大伝送スパン数が１６であることから、ウォークオフの最大値は約４３５００ｐｓと見積もられる。このウォークオフ値は、３２Ｇbaud信号に対しては約１４００シンボル分に相当するため、これより多くのシンボル数を含む条件で測定すれば、有限の時間幅を持つデータセットを用いて逆伝搬させる際にウォークオフの影響を正しく計算できる。そのため、一周期が６５５３６シンボルである波形に対して、時間軸の両端にそれぞれ５０００シンボルずつ付加した波形をもってデータセットを構成する。
　実験に先だって、シミュレーションに基づく事前の検討によって、以下の新たな事実を見出した。すなわち、一定の複雑さを持つ変調フォーマット（例えばＤＰ－１６ＱＡＭ）の信号を用いて、一定値以上の大きな入射パワー（例えばチャネル当たりの入射パワーが＋３ｄＢｍ）について得られた波形に対してのみ伝送路パラメータの学習を行えば、そこで得られたパラメータは、入射パワーが学習の際に使用した値以下である任意の変調フォーマットの信号に対する非線形波形歪み補正に適用することが可能になるということである。この事実に基づいて、実験では、チャネル当たりの入射パワーが＋３ｄＢｍであるＤＰ－１６ＱＡＭ信号を伝送、受信、且つ復調してデータセットを構成し、逆伝搬計算に用いる伝送路パラメータの学習を行う。次に、チャネル当たりの入射パワーが－５ｄＢｍ乃至＋２ｄＢｍであるＤＰ－ＰＳ－６４ＱＡＭ信号を伝送、受信、且つ復調してデータセットを構成し、先に得られた伝送路パラメータを用いた非線形波形歪み補正を実施して、信号品質の改善量を評価する。
　伝送路パラメータの学習に際しては、一回の学習ステップにおいて、ＤＰ－１６ＱＡＭ信号の４つのデータセットからいずれか一つをランダムで選び、逆伝搬計算を実施したうえで、６５５３６シンボルの中の連続する１０２４シンボルをランダムで選んで、その振幅波形を式（２０）の出力信号波形ｙｔとし、さらに対応する送信信号波形を所望信号ｄｔとして、誤差信号から各パラメータの勾配を計算することで、パラメータの更新を行う。このような学習を行うことで、限られた量のデータセットに対しても学習の過程がランダム化され、過学習を生じることなく学習を進めることが可能となる。なお学習に用いる勾配降下法のアルゴリズムとしては、上で述べたAdaBeliefを用いる。
　図１６Ａ乃至図１６Ｇは、伝送スパン数が６の場合に、１step／spanの逆伝搬条件で、実施の形態に係るＸＰＭ補償を考慮した非線形波形歪み補正のための伝送路パラメータを学習するプロセスにおいて、各種パラメータが更新される様子を示している。また図１７は、実験で用いた６スパン伝送路を構成する光ファイバの接続の様子と各スパンのファイバ長、そして逆伝搬計算における線形ステップＬ（ｊ）（０≦ｊ≦６）および非線形ステップＮ（ｊ）（１≦ｊ≦６）の計算順序を模式図として表したものである。図１６Ａ乃至図１６Ｇの結果から、伝送路パラメータと信号品質は一定回数のパラメータ更新後に収束していることがわかる。図１６Ｄ及び図１６Ｅでは、それぞれウォークオフパラメータｄ（ｊ） －５およびｄ（ｊ） －１について、ｊ＝１，３，５の場合には初期値より小さいほぼ同じ値に収束し、ｊ＝２，４，６の場合は初期値より大きいほぼ同じ値に収束している。このことは、図１５に示した長さ８４．１ｋｍのファイバと、長さ８０．５ｋｍのファイバがやや異なる分散特性を持つという事実を反映した結果である。つまり、非線形ステップＮ（ｊ）（１≦ｊ≦６）に含まれるウォークオフパラメータｄ（ｊ） ｎは、チャネル番号がｎである信号の、中心チャネル（番号がｎ＝０）の信号に対する、単位距離当たりの群遅延量に相当するパラメータであり、ファイバの分散スロープ値に比例する値であるが、ｊ＝１，３，５の場合は分散スロープ値の小さい長さ８０．５ｋｍのファイバのパラメータを、ｊ＝２，４，６の場合は分散スロープ値の大きい長さ８４．１kmのファイバのパラメータをそれぞれ正しく学習できたことを意味している。このように、実験で用いた現実の環境においても、実施の形態で提案する学習方法によって、非線形波形歪み補正に有効なパラメータを正しく学習できていることが示唆される。
　学習が完了し、各伝送スパン数に対して、非線形波形歪み補正を実施する伝送路パラメータが得られた。このパラメータにもとづき、チャネル当たりの入射パワー範囲が－４ｄＢｍから＋２ｄＢｍの範囲にある１１チャネルＰＳ－６４ＱＡＭ信号に対する非線形波形歪み補正の効果を検証した。図１８Ａ乃至図１８Ｃ及び図１９Ａ乃至図１９Ｃは、図８Ａ、図８Ｂ、図９Ａ及び図９Ｂと同様に、各伝送スパン数において、チャネル当たりの入射パワーに対する信号品質をプロットしたものである。前述のシミュレーション結果と同様に、実施の形態に係る非線形波形歪み補正の方法、すなわちチャネル内に生じるＳＰＭに加えて、チャネル間に生ずるＸＰＭによって発生する非線形シフトも考慮し、学習によるパラメータの最適化を行った１step／spanの逆伝搬計算を実施したときの信号品質が最良であることがわかる。図２０は図１０と同様に、伝送スパン数に対する信号品質の結果を、非線形波形歪み補正の有無あるいは種別ごとにプロットしたもので、やはりシミュレーション結果と同様の実験結果が得られている。以上のように、実験結果においても、実施の形態の効果が得られていることが明らかになった。
　以下では、実施の形態に係る非線形波形歪み補正を実施する場合に要する計算量と、従来技術を用いた場合に要する計算量との比較を行い、一定のチャネル数以下という条件においては、実施の形態に係る手法に要する計算量が従来手法に要する計算量と同程度であることを示す。なお、実施の形態に係る非線形波形歪み補正とは、ＳＰＭとＸＰＭによる歪みの両方を考慮し、１step／spanで逆伝搬計算を行う上で、勾配降下法を用いて伝送路パラメータの最適値を学習したうえで、逆伝搬計算によって波形補正するものであるが、ここでは、学習が終了して最適化されたパラメータを固定して逆伝搬計算する場合の計算量に注目する。このときの計算量は、１step／spanでＸＰＭを考慮し、学習を行わない逆伝搬計算の場合と同じである。また従来技術とは、ＳＰＭによる歪みの補正のみを考慮し、ＸＰＭによる歪みの補正を実施せず、２steps／spanで逆伝搬計算を行うものとする。計算量を導出する際の前提として、入力波形とは無関係に固定することが可能な数値パラメータについては、事前に計算を行ったうえでルックアップテーブル（ＬＵＴ：Look-Up Table）に結果を保存し、異なる波形入力のたびに読み出して使用することを想定し、それ自体の計算に必要な演算回数についてはこれを考慮しないものとする。
　スパン数をＳ、データ長をＮ、そしてＸＰＭ補正の計算で考慮するチャネル数をＣとおき、チャネルあたり、そして偏波あたりの計算回数を求める。ＤＳＰの計算負荷としては、積演算が主たるものであるため、その回数を求める。なお複素数どうしの積ａ×ｂ＝Re[a]Re[b]－Im[a]Im[b]＋ｉ（Re[a]Im[b]＋Re[b]Im[a]）に要する実数の積演算回数を４回として、実数の積演算回数の合計値を求める。また一般に、サイズがＮ＝２ｎである複素数信号に対するＦＦＴを実施する際の複素数の積算回数はＮ（ｌｏｇ２Ｎ－２）／２であるから、実数の積算回数はこの４倍の２Ｎ（ｌｏｇ２Ｎ－２）である。
　まず線形ステップである式（６）に要する計算量を見積もる。Ａ~は、Ａの上に～が載った記号であるとして、Ａ~ｐ，ｎ（０，ω）＝ＦＡｐ，ｎ（０，ｔ）を計算する際の実数の積算回数は２Ｎ（ｌｏｇ２Ｎ－２）である。ｅｘｐ（－αｈ／２＋ｉ（Ｄ２ω２＋Ｄ３ω３－Ｔｎω））の値は入力波形に無関係であるからＬＵＴに格納されているとして、これをＡ~ｐ，ｎ（０，ω）に積算する際の実数の積演算回数は４Ｎである。最後に逆ＦＦＴに要する演算を考慮すると、線形ステップ全体で実数の総積算回数は４Ｎ＋２×２Ｎ（ｌｏｇ２Ｎ－２）＝４Ｎ（ｌｏｇ２Ｎ－１）となる。
　続いて、非線形ステップである式（１１）乃至（１３）の計算に要する計算量を見積もる。実数である非線形位相シフト量φが得られたとして、ｅｉφの計算に必要な実数の積算回数は、以下の４次のTaylor展開において、φ２／２の計算に２Ｎ回、φ２の結果を再利用することで(φ２)２／２４の計算に２Ｎ回、そしてφ２×φ／６の計算に２Ｎ回必要であることから、合計６Ｎ回となる。
Figure JPOXMLDOC01-appb-M000040
　また複素数同士の積であるＢｐ，０（０，ｔ）×ｅｉφを計算する際の実数の積算回数は４Ｎ回である。次に、式（１２）でＳＰＭによる位相シフトφＳＰＭ（ｔ）を計算する際に、Ｐｐ，０（ｔ）＝Re[Ｂｐ，０（ｔ）]２＋Im[Ｂｐ，０（ｔ）]２であるから、実数の積算回数は２Ｎ回であり、さらに係数ｇＨ０をかけることを考慮すると、実数の積算回数は合計３Ｎ回である。さらに、別途供給される直交偏波成分の強度Ｐ３－ｐ，０（ｔ）に係数ｇδＨ０をかける際にもＮ回必要なので、結局φＳＰＭ（ｔ）を計算する際の実数の積算回数は４Ｎ回である。
　次に、式（１３）からＸＰＭによる位相シフトφＸＰＭ（ｔ）を計算する際に必要な計算量を見積もる。強度波形Ｐｐ，０（ｔ）は、φＳＰＭ（ｔ）を計算する際にすでに得られており、これにＦＦＴを適用してＰ~ｐ，０（ω）を計算するための実数の積算回数は２Ｎ（ｌｏｇ２Ｎ－２）である。Ｐ~ｐ，０（ω）は式（１３）には現れないが、他のチャネルの非線形波形歪み補正を実施するために供給する必要があるため、計算量を考慮に入れる。一方で、式（１３）に現れるｎ≠０のＰ~ｐ，ｎ（ω）は、別途計算されて供給されるものとし、必要な計算量を考慮しない。Ｐ~ｐ，ｎ（ω）およびＰ~３－ｐ，ｎ（ω）にそれぞれ実数の係数である２ｇおよびδｇをかける際に２Ｎ回ずつ、合計４Ｎ回の実数の積演算が必要である。式（１４）のＨｎ（ｈ，ω）はＬＵＴに保存されているものとし、２ｇＰ~ｐ，ｎ（ω）＋δｇＰ~３－ｐ，ｎ（ω）の結果にＨｎ（ｈ，ω）をかける際に必要な実数の積算回数は４Ｎであるため、Ｈｎ（ｈ，ω）｛２ｇＰ~ｐ，ｎ（ω）＋δｇＰ~３－ｐ，ｎ（ω）｝を得るために必要な実数の積算回数は合計で８Ｎである。この演算がｎ≠０であるＣ－１チャネル分必要なので、実数の積算回数は８Ｎ（Ｃ－１）である。最後に逆ＦＦＴを適用する際の計算量を考慮すると、φＸＰＭ（ｔ）を計算する際に必要な計算量は、２Ｎ（ｌｏｇ２Ｎ－２）＋８Ｎ（Ｃ－１）＋２Ｎ（ｌｏｇ２Ｎ－２）＝４Ｎ（ｌｏｇ２Ｎ＋２Ｃ－４）である。以上より、一つの非線形ステップにおいてφＸＰＭ（ｔ）を考慮しない場合の計算量は６Ｎ＋４Ｎ＋４Ｎ＝１４Ｎであり、ＸＰＭ補償を実施するためにφＸＰＭ（ｔ）を考慮した場合の計算量は、１４Ｎ＋４Ｎ（ｌｏｇ２Ｎ＋２Ｃ－４）となる。
　スパン数Ｓの伝送路において、スパンあたりのステップ数をＭとする逆伝搬計算では、線形ステップの数が合計ＭＳ＋１、非線形ステップの数が合計ＭＳである。上に示した結果から、ＸＰＭ補償を行わない場合の実数積の回数は、以下のとおりである。
４Ｎ（ｌｏｇ２Ｎ－１）×（ＭＳ＋１）＋１４Ｎ×ＭＳ
　ＸＰＭ補償を行う場合の実数席の回数は、以下のとおりである。
４Ｎ（ｌｏｇ２Ｎ－１）×（ＭＳ＋１）＋｛１４Ｎ＋４Ｎ（ｌｏｇ２Ｎ＋２Ｃ－４）｝×ＭＳ
　図２１は、スパン数をＳ＝１０、チャネル数をＣ＝５，１１，２１としたときに、本実施の形態としてＭ＝１step／spanでＸＰＭ補償ありとする場合と、従来方式としてＭ＝２steps／spanでＸＰＭ補償なしとする場合についての実数の積算回数を、データ長Ｎに対してプロットしたものである。チャネル数が５のとき、本実施の形態は従来方式とほぼ同じ計算量で実施可能である。チャネル数が１１のときは、本実施の形態が従来方式よりやや大きくなって１．６倍となるが、オーダーとしてほぼ同じ計算量にとどまり、前に９チャネルの場合のシミュレーション結果と１１チャネルの場合の実験結果を示したように、大きな非線形波形歪み補正効果が得られる。チャネル数が２１のときは本実施の形態に必要な計算量が従来方式の約２．４倍に上昇するが、依然として同じオーダーに留まると言える。
　以上、本発明の実施の形態について説明したが、本発明はこれらに限定されるものではない。例えば、確率的勾配降下法を用いる例を示したが、勾配降下法の各種バリエーションを適用可能である。また、上でも述べたように、ＳＰＭ及びＸＰＭの影響を全てのステップにおいて考慮するのではなく、一部のステップにのみ両方の影響を考慮したとしても、十分な効果が得られる場合がある。
　なお、ＤＳＰは、演算器とメモリを含んだものである。また、上で述べた処理を実行するのはＤＳＰに限らず、他のプロセッサであっても良い。さらに、上で述べた処理をプロセッサに実行させるプログラムは、不揮発性メモリに記録され、実行時には、当該プログラムに含まれる命令がプロセッサに読み出されて実行される。さらに、専用の回路又は専用の回路とＤＳＰ等との組み合わせにて、上で述べた処理を実行するようにしても良い。
　以上述べた実施の形態をまとめると以下のようになる。
　本実施の形態における第１の側面に係る光波形歪み補正方法は、伝送路である光ファイバの非線形光学効果と群速度分散効果によって波形の形状が変化した光信号を受信して波形を数値化した後、非線形シュレディンガー方程式の線形項と非線形項を交互に演算することで送信時の波形を推定することで光波形の歪みを補正する光波形歪み補正方法であって、チャネル内で生じる自己位相変調に起因する波形歪みのみならず、波長分割多重伝送時にチャネル間で生じる相互位相変調に起因する波形歪みも考慮に入れて計算し、さらに計算に用いるパラメータを勾配降下法によって最適化し、伝送路１スパンあたりのステップ数を２未満としたことを特徴とするものである。計算負荷を抑制しつつ、精度を向上させることが出来る。
　上で述べた伝送路１スパンあたりのステップ数が、１以下であってもよい。このようにステップ数を減らしても、精度を向上させることが出来る。
　さらに、上で述べたパラメータが、二次の群速度分散、非線形係数及びウォークオフを含むようにしてもよい。
　本実施の形態における第２の側面に係る光波形歪み補正方法は、（Ａ）伝送路において波形の形状が変化した光信号を受信して波形を数値化した後、非線形シュレディンガー方程式の線形項と非線形項とを交互に演算することで送信時の波形を推定する逆伝搬処理であって、波長分割多重伝送時における上記伝送路内の複数のチャネルの各々について、当該チャネル内で生ずる自己位相変調（ＳＰＭ）に起因する波形歪みと、当該チャネル以外のチャネルとの間で生じる相互位相変調（ＸＰＭ）に起因する波形歪みとを補正する逆伝搬処理において用いられ且つ相互位相変調に関連する第１のパラメータと、上記逆伝搬処理において用いられ且つ自己位相変調及び相互位相変調に関連する第２のパラメータとを、勾配降下法によって最適化するステップと、（Ｂ）最適化された第１及び第２のパラメータを用いて、上記逆伝搬処理を実行するステップとを含む。
　このように、第２のパラメータ（例えば、Ｄ２、Ｄ３、ｇ、δ、α０）だけではなく、第１のパラメータ（例えば、αｎ、ｄｎ）についても最適化すると共に、最適化された第１及び第２のパラメータを用いてＳＰＭ及びＸＰＭに起因する波形歪みを補正する逆伝搬処理を実行することで、伝送路１スパンあたりのステップ数を減少させて計算負荷を抑えても、十分な計算精度を得ることが出来るようになる。
　なお、相互位相変調に起因する上記波形歪みを、伝搬距離に依存せず強度について初期波形が維持されるが、伝搬距離に比例した遅延が時間軸上で発生するという近似の下、補正を行うようにしても良い。このような近似により補正（補償とも呼ぶ）の計算負荷がより抑制できる。
　さらに、上記第２のパラメータは、群速度分散Ｄ２及び非線形係数ｇを含み、上記第１のパラメータは、ウォークオフパラメータｄｎを含むようにしても良い。最適化するパラメータを絞り込こめば、さらに計算負荷をさらに抑制できる。

Claims (5)

  1. 　伝送路において波形の形状が変化した光信号を受信して波形を数値化した後、非線形シュレディンガー方程式の線形項と非線形項とを交互に演算することで送信時の波形を推定する逆伝搬処理であって、波長分割多重伝送時における前記伝送路内の複数のチャネルの各々について、当該チャネル内で生ずる自己位相変調に起因する波形歪みと、当該チャネル以外のチャネルとの間で生じる相互位相変調に起因する波形歪みとを補正する逆伝搬処理において用いられ且つ前記相互位相変調に関連する第１のパラメータと、前記逆伝搬処理において用いられ且つ前記自己位相変調及び前記相互位相変調に関連する第２のパラメータとを、勾配降下法によって最適化するステップと、
    　最適化された前記第１及び第２のパラメータを用いて、前記逆伝搬処理を実行するステップと、
    　を含む光波形歪み補正方法。
  2. 　前記相互位相変調に起因する波形歪みを、伝搬距離に依存せず強度について初期波形が維持されるが、伝搬距離に比例した遅延が時間軸上で発生するという近似の下、補正を行う
    　ことを特徴とする請求項１記載の波形歪み補正方法。
  3. 　前記第１のパラメータは、ウォークオフパラメータｄｎを含み、
    　前記第２のパラメータは、群速度分散Ｄ２及び非線形係数ｇを含む
    　請求項１又は２記載の波形歪み補正方法。
  4. 　伝送路において波形の形状が変化した光信号を受信して波形を数値化した後、非線形シュレディンガー方程式の線形項と非線形項とを交互に演算することで送信時の波形を推定する逆伝搬処理であって、波長分割多重伝送時における前記伝送路内の複数のチャネルの各々について、当該チャネル内で生ずる自己位相変調に起因する波形歪みと、当該チャネル以外のチャネルとの間で生じる相互位相変調に起因する波形歪みとを補正する逆伝搬処理において用いられ且つ前記相互位相変調に関連する第１のパラメータと、前記逆伝搬処理において用いられ且つ前記自己位相変調及び前記相互位相変調に関連する第２のパラメータとを、勾配降下法によって最適化する最適化部と、
    　最適化された前記第１及び第２のパラメータを用いて、前記逆伝搬処理を実行する逆伝搬処理部と、
    　を有する光波形歪み補正装置。
  5. 　請求項４記載の光波形歪み補正装置を有する光信号受信装置。
PCT/JP2022/000917 2021-03-25 2022-01-13 光波形歪み補正方法及び装置、並びに光信号受信装置 WO2022201763A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US18/283,208 US20240178919A1 (en) 2021-03-25 2022-01-13 Method and Apparatus for Correcting Optical Waveform Distortion and Optical Signal Receiving Apparatus
JP2023508663A JPWO2022201763A1 (ja) 2021-03-25 2022-01-13

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2021052107 2021-03-25
JP2021-052107 2021-03-25

Publications (1)

Publication Number Publication Date
WO2022201763A1 true WO2022201763A1 (ja) 2022-09-29

Family

ID=83396865

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2022/000917 WO2022201763A1 (ja) 2021-03-25 2022-01-13 光波形歪み補正方法及び装置、並びに光信号受信装置

Country Status (3)

Country Link
US (1) US20240178919A1 (ja)
JP (1) JPWO2022201763A1 (ja)
WO (1) WO2022201763A1 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015201726A (ja) * 2014-04-07 2015-11-12 富士通株式会社 デジタルコヒーレント受信器および受信方法
WO2018215850A1 (en) * 2017-05-24 2018-11-29 Ciena Corporation Adjustment of control parameters of section of optical fiber network and power control in an optical fiber network
JP2019118018A (ja) * 2017-12-27 2019-07-18 三菱電機株式会社 光通信装置、光通信方法、及び制御プログラム
US20190266480A1 (en) * 2018-02-26 2019-08-29 Nec Laboratories America, Inc Single-step nonlinearity compensation using artificial intelligence for digital coherent transmission systems

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015201726A (ja) * 2014-04-07 2015-11-12 富士通株式会社 デジタルコヒーレント受信器および受信方法
WO2018215850A1 (en) * 2017-05-24 2018-11-29 Ciena Corporation Adjustment of control parameters of section of optical fiber network and power control in an optical fiber network
JP2019118018A (ja) * 2017-12-27 2019-07-18 三菱電機株式会社 光通信装置、光通信方法、及び制御プログラム
US20190266480A1 (en) * 2018-02-26 2019-08-29 Nec Laboratories America, Inc Single-step nonlinearity compensation using artificial intelligence for digital coherent transmission systems

Also Published As

Publication number Publication date
JPWO2022201763A1 (ja) 2022-09-29
US20240178919A1 (en) 2024-05-30

Similar Documents

Publication Publication Date Title
Yan et al. Low complexity digital perturbation back-propagation
Gao et al. Assessment of intrachannel nonlinear compensation for 112 Gb/s dual-polarization 16QAM systems
US10439730B2 (en) Method, device and communication system for reducing optical transmission impairments
US9002210B2 (en) Method and system for non-linearity compensation in optical transmission systems
US10903912B2 (en) Method and apparatus for nonlinear compensation in coherent optical links
EP2858272B1 (en) Non-linear distortion compensator, method of compensating non-linear distortion, and optical receiver
US20220149954A1 (en) Asymmetric direct detection of optical signals
Do et al. Data-aided OSNR estimation for QPSK and 16-QAM coherent optical system
Li et al. Fiber nonlinearity mitigation in single carrier 400 G and 800 G Nyquist-WDM systems
Goroshko et al. Overcoming performance limitations of digital back propagation due to polarization mode dispersion
Sillekens et al. Time-domain learned digital back-propagation
Asif et al. Logarithmic step-size based digital backward propagation in N-channel 112Gbit/s/ch DP-QPSK transmission
Kumar et al. Enhanced regular perturbation-based nonlinearity compensation technique for optical transmission systems
Inoue et al. Learning-based digital back propagation to compensate for fiber nonlinearity considering self-phase and cross-phase modulation for wavelength-division multiplexed systems
EP2723004A1 (en) Method for non-linear equalization of the optical channel in the frequency domain
Li et al. Inter-channel fiber nonlinearity mitigation in high baud-rate optical communication systems
Mussolin et al. Polarization multiplexed 224 Gb/s 16QAM transmission employing digital back-propagation
Nakashima et al. Digital nonlinear compensation technologies in coherent optical communication systems
EP3352387A1 (en) Digital multi-channel nonlinearity compensation scheme for optical coherent communication
US11777612B2 (en) Method for nonlinear compensation of coherent high-capacity high-order qam system
WO2022201763A1 (ja) 光波形歪み補正方法及び装置、並びに光信号受信装置
CN107888291B (zh) 一种双向拉曼放大的光纤通信中信号补偿系统及方法
Nouri et al. Chromatic dispersion and nonlinear phase noise compensation based on KLMS method
Liang et al. Perturbation-assisted DBP for nonlinear compensation in polarization multiplexed systems
He et al. A layer-reduced neural network based digital backpropagation algorithm for fiber nonlinearity mitigation

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22774568

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2023508663

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 18283208

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 22774568

Country of ref document: EP

Kind code of ref document: A1