Signal processing method, signal processing device and communication system

Abstract

A signal processing method for executing: a first equalization process for converting the real and imaginary components of each polarized wave of a polarized-and-multiplexed reception signal to frequency domain signals, inputting the frequency domain signal of the real component and the frequency domain signal of the imaginary component of each polarized wave, as well as a signal derived by inverting the frequencies of each of the frequency domain signal of the real component and the frequency domain signal of the imaginary component of each polarized wave and taking a complex conjugate, as input signals, multiplying a complex transfer function by each of the real component and the imaginary component of each polarized wave, subsequently adding the resultant products together for each polarized wave, and inverse-converting from the frequency domain signal to a time domain signal; and a second equalization process for multiplying a complex transfer function by each of a signal of the real component derived by inverting the frequency of the real component of each polarized wave included in the input signals and taking a complex conjugate and a signal of the imaginary component derived by inverting the frequency of the imaginary component and taking a complex conjugate, subsequently adding the resultant products together, and inverse-converting from the frequency domain signal to a time domain signal, the signal processing method also including adding or subtracting a transmission data bias correction signal to or from a signal derived by adding together the first and second added signals.　

Description

[Technical Field]

[0001] The present invention relates to a signal processing method, a signal processing device, and a communication system. [Background technology]

[0002] In digital coherent transmission, it is necessary not only to compensate for waveform distortion occurring in the optical fiber transmission line, but also to adaptively compensate for device imperfections in the optical transmitter and receiver. In adaptive equivalent circuits used in general signal processing, compensation is mainly performed for waveform distortion occurring in the transmission line, and compensation for device imperfections in the transmitter and receiver must be performed separately in subsequent signal processing. Therefore, there is a technology that simultaneously compensates for device imperfections in the transmitter and receiver (see, for example, Patent Document 1 and Non-Patent Document 1). [Prior art documents] [Patent documents]

[0003] [Patent Document 1] Japanese Patent Application Publication No. 2020-141294 [Non-patent literature]

[0004] [Non-Patent Document 1] Takayuki Kobayashi, et. al., “35-Tb / s C-Band Transmission Over 800 km Employing 1-Tb / s PS-64QAM Signals Enhanced by Complex 8 × 2 MIMO Equalizer”, Optical Fiber Communication Conference Postdeadline Papers 2019, Th4B.2 Summary of the Invention [Problem to be solved by the invention]

[0005] The adaptive equalization circuits of Patent Document 1 and Non-Patent Document 1 have a different configuration from a 2x2 MIMO (Multiple Input Multiple Output) adaptive equalization circuit with complex inputs and complex outputs that is commonly used in conventional optical communications. In the adaptive equalization circuits of Patent Document 1 and Non-Patent Document 1, the generated tap coefficients and the like have no commonality or compatibility, and the total number of taps increases. As a result, there is a problem in that the amount of calculation increases exponentially as the number of taps increases.

[0006] In view of the above circumstances, an object of the present invention is to provide a technique that can perform equalization processing while reducing the amount of calculation in digital coherent optical transmission. [Means for solving the problem]

[0007] One aspect of the present invention includes a conversion step of converting real and imaginary components of each polarization of a polarization-multiplexed received signal into frequency domain signals; a signal input step of inputting as input signals the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization, and the converted frequency domain signals obtained by performing frequency inversion on the frequency axis for each of the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization and taking complex conjugates; a first equalization step of multiplying the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization included in the input signal by a complex transfer function, and then adding them up, thereby inversely converting the frequency domain signals into time domain signals; the frequency domain signal after the conversion of the real component of the polarization and the frequency domain signal after the conversion of the imaginary component of the polarization are multiplied by a complex transfer function, respectively, and then added together to perform second equalization processing for inversely converting the frequency domain signal into a time domain signal; and the compensation steps of, for each polarization, performing a phase rotation for frequency offset compensation on the time domain signal converted by the first equalization processing to generate a first sum signal, performing a phase rotation opposite to the phase rotation for frequency offset compensation on the time domain signal converted by the second equalization processing to generate a second sum signal, and adding or subtracting a transmit data bias correction signal to or from a signal obtained by adding the first sum signal and the second sum signal.

[0008] One aspect of the present invention includes an addition step of performing an imaginary unit multiplication process of multiplying an imaginary component of each polarization of a polarization-multiplexed received signal by an imaginary unit j, and then performing an addition process of adding the imaginary component multiplied by the imaginary unit j to a real component of each polarization of the polarization-multiplexed received signal; a conversion step of converting the signal obtained by the addition process of the imaginary component multiplied by the imaginary unit j to a frequency domain signal; and and for the frequency domain signal The frequency domain signal after frequency inversion and complex conjugation a first computed frequency domain signal obtained by adding the frequency domain signal after compensation and the frequency domain signal after compensation of the converted frequency domain signal and then multiplying by 1 / 2; a computed frequency domain signal; From the frequency domain signal, The frequency domain signal after frequency inversion and complex conjugation a second computed frequency domain signal obtained by subtracting the frequency domain signal after compensation from the frequency domain signal after compensation has been performed on the converted frequency domain signal and then multiplying by 1 / 2j;a signal input step of inputting the calculated frequency domain signal as an input signal, and No. 1 of the computed frequency domain signal and the imaginary component No. 1 A first equalization process in which each of the calculated frequency domain signals is multiplied by a complex transfer function and then added, and the frequency domain signal is inversely converted to a time domain signal; No. 2 of the computed frequency domain signal and the imaginary component No. 2 The signal processing method includes: an equalization step of multiplying each of the calculated frequency domain signals by a complex transfer function, adding the multiplied signals, and performing a second equalization process of inversely converting the frequency domain signals into time domain signals; and a compensation step of, for each polarization, performing a phase rotation for frequency offset compensation on the time domain signals converted by the first equalization process to generate a first sum signal, performing a phase rotation opposite to the phase rotation for frequency offset compensation on the time domain signals converted by the second equalization process to generate a second sum signal, and adding or subtracting a transmit data bias correction signal to or from a signal obtained by adding the first sum signal and the second sum signal.

[0009] One aspect of the present invention includes a frequency conversion unit that converts real and imaginary components of each polarization of a polarization-multiplexed received signal into frequency domain signals; a signal input unit that inputs as input signals the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization, and frequency domain signals after conversion in which frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization are frequency-inverted on the frequency axis and complex conjugates are taken for the frequency domain signals; a first equalization process that multiplies the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization included in the input signal by a complex transfer function, and then adds them together to inversely convert the frequency domain signals into time domain signals; and a compensation unit that, for each polarization, performs a phase rotation for frequency offset compensation on the time domain signal converted by the first equalization processing to generate a first sum signal, performs a phase rotation opposite to the phase rotation for frequency offset compensation on the time domain signal converted by the second equalization processing to generate a second sum signal, and adds a transmission data bias correction signal to or from a signal obtained by adding the first sum signal and the second sum signal.

[0010] According to one aspect of the present invention, there is provided an adder that performs an imaginary unit multiplication process of multiplying an imaginary component of each polarization of a polarization-multiplexed received signal by an imaginary unit j, and then performs an addition process of adding the imaginary component multiplied by the imaginary unit j to a real component of each polarization of the polarization-multiplexed received signal; a frequency converter that converts the signal obtained by the addition process of the imaginary component multiplied by the imaginary unit j to a frequency domain signal; and a frequency converter that converts the frequency domain signal of each polarization into a frequency domain signal. and for the frequency domain signal The frequency domain signal after frequency inversion and complex conjugation a first computed frequency domain signal obtained by adding the frequency domain signal after compensation and the frequency domain signal after compensation of the converted frequency domain signal and then multiplying by 1 / 2; a computed frequency domain signal; From the frequency domain signal, The frequency domain signal after frequency inversion and complex conjugation a second computed frequency domain signal obtained by subtracting the frequency domain signal after compensation from the frequency domain signal after compensation has been performed on the converted frequency domain signal and then multiplying by 1 / 2j;a signal input unit that inputs the calculated frequency domain signal as an input signal, and a signal input unit that inputs the calculated frequency domain signal as an input signal, and No. 1 of the computed frequency domain signal and the imaginary component No. 1 A first equalization process in which each of the calculated frequency domain signals is multiplied by a complex transfer function and then added, and the frequency domain signal is inversely converted to a time domain signal; No. 2 of the computed frequency domain signal and the imaginary component No. 2 the signal processing device includes: an equalization unit that performs second equalization processing by multiplying each of the calculated frequency domain signals by a complex transfer function, adding the multiplied signals, and inversely converting the frequency domain signals into time domain signals; and a compensation unit that, for each polarization, performs a phase rotation for frequency offset compensation on the time domain signals converted by the first equalization processing to generate a first sum signal, performs a phase rotation opposite to the phase rotation for frequency offset compensation on the time domain signals converted by the second equalization processing to generate a second sum signal, and adds or subtracts a transmit data bias correction signal to or from a signal obtained by adding the first sum signal and the second sum signal.

[0011] One aspect of the present invention is a communication system including a transmitter that transmits a polarization multiplexed signal, and a receiver that includes the signal processing device described above. [Effects of the Invention]

[0012] According to the present invention, it is possible to perform equalization processing while reducing the amount of calculation in digital coherent optical transmission. [Brief explanation of the drawings]

[0013] [Figure 1] 1 is a diagram illustrating an example of the configuration of a digital coherent optical transmission system according to a first embodiment. [Figure 2] FIG. 2 is a diagram illustrating an example of the configuration of a demodulation digital signal processing unit according to the first embodiment. [Figure 3] FIG. 10 illustrates an example of a coefficient calculation unit. [Figure 4]FIG. 10 illustrates an example of a coefficient calculation unit. [Figure 5] FIG. 10 illustrates an example of a coefficient calculation unit. [Figure 6] FIG. 10 illustrates an example of a coefficient calculation unit. [Figure 7] FIG. 10 is a diagram illustrating an example of the configuration of a demodulation digital signal processing unit in a modified example of the first embodiment. [Figure 8] FIG. 10 is a diagram for explaining the effects of the present invention. [Figure 9] FIG. 10 is a diagram for explaining the effects of the present invention. [Figure 10] FIG. 10 is a diagram illustrating an example of the configuration of a demodulation digital signal processing unit in the second embodiment. [Figure 11] FIG. 10 is a diagram illustrating an example of the configuration of a demodulation digital signal processing unit in a modified example of the second embodiment. [Figure 12] FIG. 10 is a diagram illustrating an example of the configuration of a demodulation digital signal processing unit according to a third embodiment. DETAILED DESCRIPTION OF THE INVENTION

[0014] Hereinafter, an embodiment of the present invention will be described with reference to the drawings. (First embodiment) 1 is a diagram showing an example of the configuration of a digital coherent optical transmission system 1 according to the first embodiment. The digital coherent optical transmission system 1 includes a transmitter 10 and a receiver 50. The transmitter 10 transmits a polarization multiplexed signal. The receiver 50 receives the polarization multiplexed signal from the transmitter 10.

[0015] The transmitter 10 has at least one transmitting unit 100. The transmitting unit 100 outputs a polarization multiplexed signal of a specified wavelength to an optical fiber transmission line 30. The optical fiber transmission line 30 is provided with any number of optical amplifiers 31. Each optical amplifier 31 receives a polarization multiplexed signal from the optical fiber transmission line 30 on the transmitter 10 side, amplifies the signal, and outputs it to the optical fiber transmission line 30 on the receiver 50 side. The receiver 50 has at least one receiving unit 500. The receiving unit 500 receives the polarization multiplexed signal.

[0016] First, the configuration of the transmitter 10 will be described. The transmitting unit 100 includes a digital signal processing unit 110, a modulator driver 120, a light source 130, and an integrated module 140. The digital signal processing unit 110 includes an encoding unit 111, a mapping unit 112, a training signal inserting unit 113, a frequency converting unit 114, a waveform shaping unit 115, a pre-equalization unit 116, and digital-to-analog converters (DACs) 117-1 to 117-4.

[0017] The encoding unit 111 performs FEC (forward error correction) encoding on the transmission bit string and outputs the resulting transmission signal.

[0018] Mapping section 112 maps the transmission signal output from encoding section 111 to symbols.

[0019] The training signal inserting unit 113 inserts a known training signal into the transmission signal symbol-mapped by the mapping unit 112 .

[0020] The frequency conversion unit 114 performs upsampling by changing the sampling frequency of the transmission signal into which the training signal has been inserted.

[0021] The waveform shaping unit 115 limits the band of the sampled transmission signal.

[0022] Pre-equalization section 116 compensates for distortion in the waveform of the transmission signal band-limited by waveform shaping section 115, and outputs the signal to DACs 117-1 to 117-4.

[0023] DAC 117-1 converts the I (in-phase) component of the X polarization of the transmission signal input from pre-equalization unit 116 from a digital signal to an analog signal and outputs the signal to modulator driver 120. DAC 117-2 converts the Q (quadrature) component of the X polarization of the transmission signal input from pre-equalization unit 116 from a digital signal to an analog signal and outputs the signal to modulator driver 120. DAC 117-3 converts the I component of the Y polarization of the transmission signal input from pre-equalization unit 116 from a digital signal to an analog signal and outputs the signal to modulator driver 120. DAC 117-4 converts the Q component of the Y polarization of the transmission signal input from pre-equalization unit 116 from a digital signal to an analog signal and outputs the signal to modulator driver 120.

[0024] The modulator driver 120 has amplifiers 121-1 to 121-4. The amplifier 121-i (i is an integer between 1 and 4) amplifies the analog signal output from the DAC 117-i, and drives the modulator of the integrated module 140 with the amplified analog signal.

[0025] The light source 130 is, for example, a semiconductor laser (LD), and outputs light of a specified wavelength.

[0026] The integrated module 140 includes IQ modulators 141-1 and 141-2 and a polarization combining unit 142. The IQ modulator 141-1 generates an X-polarized optical signal by modulating the optical signal output by the light source 130 based on the I component of the X-polarized wave output by the amplifier 121-1 and the Q component of the X-polarized wave output by the amplifier 121-2. The IQ modulator 141-2 generates a Y-polarized optical signal by modulating the optical signal output by the light source 130 based on the I component of the Y-polarized wave output by the amplifier 121-3 and the Q component of the Y-polarized wave output by the amplifier 121-4. The polarization combining unit 142 generates a polarization multiplexed signal by polarization multiplexing the X-polarized optical signal generated by the IQ modulator 141-1 and the Y-polarized optical signal generated by the IQ modulator 141-2. The polarization multiplexing unit 142 outputs the generated polarization multiplexed signal to the optical fiber transmission line 30 .

[0027] Next, the configuration of the receiver 50 will be described. The receiving unit 500 includes a local oscillation light source 510, an optical front end 520, and a digital signal processing unit 530. The local oscillation light source 510 is, for example, an LD. The local oscillation light source 510 outputs local oscillation light (LO).

[0028] The optical front end 520 converts the optical signal into an electrical signal while maintaining the phase and amplitude of the polarization-multiplexed phase-modulated signal. The optical front end 520 includes a polarization splitter 521, optical 90-degree hybrid couplers 522-1 and 522-2, BPDs (Balanced Photo Diodes) 523-1 to 523-4, and amplifiers 524-1 to 524-4.

[0029] The polarization separator 521 separates the input optical signal into an X-polarized optical signal and a Y-polarized optical signal. The polarization separator 521 outputs the X-polarized optical signal to the optical 90-degree hybrid coupler 522-1 and outputs the Y-polarized optical signal to the optical 90-degree hybrid coupler 522-2.

[0030] The optical 90-degree hybrid coupler 522-1 causes interference between the X-polarized optical signal and the local oscillator light output from the local oscillator light source 510, and extracts an I-component optical signal and a Q-component optical signal from the received optical field. The optical 90-degree hybrid coupler 522-1 outputs the extracted X-polarized I-component optical signal and Q-component optical signal to the BPDs 523-1 and 523-2.

[0031] The optical 90-degree hybrid coupler 522-2 causes interference between the Y-polarized optical signal and the local oscillation light output from the local oscillation light source 510, and extracts the I and Q components of the received optical field. The optical 90-degree hybrid coupler 522-2 outputs the extracted I and Q components of the Y-polarized wave to the BPDs 523-3 and 523-4.

[0032] BPDs 523-1 to 523-4 are differential input photoelectric converters. BPD 523-i outputs the difference between the photocurrents generated in two photodiodes with matching characteristics to amplifier 524-i. BPD 523-1 converts the I component of the X-polarized received signal into an electric signal and outputs it to amplifier 524-1. BPD 523-2 converts the Q component of the X-polarized received signal into an electric signal and outputs it to amplifier 524-2. BPD 523-3 converts the I component of the Y-polarized received signal into an electric signal and outputs it to amplifier 524-3. BPD 523-4 converts the Q component of the Y-polarized received signal into an electric signal and outputs it to amplifier 524-4. Amplifier 524-i (i is an integer between 1 and 4) amplifies the electric signal output from BPD 523-i and outputs it to digital signal processing unit 530.

[0033] The digital signal processing unit 530 includes analog-to-digital converters (ADCs) 531-1 to 531-4, a demodulation digital signal processing unit 532, a demapping unit 533, and a decoding unit 534.

[0034] The ADC 531 - i (i is an integer between 1 and 4) converts the electrical signal output from the amplifier 524 - i from an analog signal to a digital signal, and outputs the digital signal to the demodulation digital signal processing unit 532 .

[0035] The demodulation digital signal processing unit 532 receives the I component of the X-polarized received signal from the ADC 531-1, the Q component of the X-polarized received signal from the ADC 531-2, the I component of the Y-polarized received signal from the ADC 531-3, and the Q component of the Y-polarized received signal from the ADC 531-4. The demodulation digital signal processing unit 532 performs signal processing such as at least equalization processing and compensation for frequency offset and phase noise on each of the input signals. The demodulation digital signal processing unit 532 also performs signal processing such as compensation for frequency characteristics and chromatic dispersion as necessary.

[0036] Whether or not the demodulation digital signal processing unit 532 performs signal processing such as compensation for frequency characteristics and compensation for chromatic dispersion depends on the configuration of the demodulation digital signal processing unit 532. Therefore, this will be explained in detail when explaining the configuration of the demodulation digital signal processing unit 532. The demodulation digital signal processing unit 532 is one aspect of a signal processing device.

[0037] The demapping unit 533 determines the symbols of the received signal output by the demodulation digital signal processing unit 532, and converts the determined symbols into binary data.

[0038] The decoding unit 534 performs error correction decoding processing such as FEC on the binary data demapped by the demapping unit 533 to obtain a received bit string.

[0039] Although the above embodiment describes an example of a single optical fiber transmission line, the same applies to a spatially multiplexed transmission system (for example, a multi-core fiber, a multi-mode fiber, and free space transmission).

[0040] Next, the configuration of the demodulation digital signal processing unit 532 will be described. Fig. 2 is a diagram showing an example of the configuration of the demodulation digital signal processing unit 532 in the first embodiment. The demodulation digital signal processing unit 532 shown in Fig. 2 performs signal processing such as equalization processing and compensation for frequency offset and phase noise. Note that the demodulation digital signal processing unit 532 shown in Fig. 2 does not perform signal processing such as compensation for frequency characteristics and compensation for chromatic dispersion.

[0041] The demodulation digital signal processing unit 532 includes an adaptive equalization unit 54 and a frequency / phase offset compensation unit 55. The adaptive equalization unit 54 adaptively performs equalization processing on each input signal. The frequency / phase offset compensation unit 55 performs processing such as compensation for frequency offset and phase noise on the received signal that has been equalized by the adaptive equalization unit 54.

[0042] Next, the operation of demodulation digital signal processing unit 532 will be described. Adaptive equalization unit 54 of demodulation digital signal processing unit 532 receives real component XI and imaginary component XQ of the X-polarized received signal, which have been converted into digital signals by ADCs 531-1 to 531-4, and real component YI and imaginary component YQ of the Y-polarized received signal. Adaptive equalization unit 54 stores the received real component XI, imaginary component XQ, real component YI, and imaginary component YQ in the corresponding buffers. These buffers correspond to those used in the overlap save method described in Reference 1 below. (Reference 1: John J. Shynk, “Frequency-Domain and Multirate Adaptive Filtering”, January 1992.)

[0043] The adaptive equalizer 54 performs an N-point discrete Fourier transform or fast Fourier transform (corresponding to "N-DFT" shown in FIG. 2) on each of the real component XI, imaginary component XQ, real component YI, and imaginary component YQ stored in the buffer. This converts the real and imaginary components of each polarization into frequency-domain signals. That is, the adaptive equalizer 54 generates a frequency-domain signal for the real component XI, a frequency-domain signal for the imaginary component XQ, a frequency-domain signal for the real component YI, and a frequency-domain signal for the imaginary component YQ. In the following explanation, when using numerical values, a case where the number of N points is 256 will be described as an example.

[0044] The frequency domain signal of real component XI, the frequency domain signal of imaginary component XQ, the frequency domain signal of real component YI, and the frequency domain signal of imaginary component YQ generated by the adaptive equalization unit 54 are each branched into four by a branching unit. Of the four branched frequency domain signals, two frequency domain signals are input directly to the coefficient calculation unit, and the remaining two frequency domain signals are converted into inverted and complex conjugated frequency domain signals and input to the coefficient calculation unit.

[0045] For example, in the case of a frequency domain signal of real component XI, the frequency domain signal of real component XI is branched into four by the branching unit, and two of the four branched frequency domain signals of real component XI are input directly to the coefficient calculation unit, while the remaining two frequency domain signals are converted into inverted and complex conjugated frequency domain signals by the inversion and complex conjugation unit and then input to the coefficient calculation unit.

[0046] Here, the inverted and complex conjugated frequency domain signal is a signal obtained by inverting the frequency domain signal around DC (Direct current: direct current component, frequency 0) in the frequency domain and taking the complex conjugate in order to generate a complex conjugate signal in the time domain and realize the equivalent operation in the frequency domain. Consider a signal X(f) in the frequency domain, and the inverted and complex conjugated signal is X † Hereinafter, the frequency domain signal of the real component converted by the inversion / complex conjugation unit will be referred to as the "real component inverted complex conjugate signal," and the frequency domain signal of the imaginary component will be referred to as the "imaginary component inverted complex conjugate signal."

[0047] The coefficient calculation section calculates impulse responses H1 to H 16 In FIG. 2, the coefficient calculation unit performs multiplication of the complex transfer function of impulse responses H1 to H 16 Although only the values ​​are shown, the specific configuration of the coefficient calculation unit will be described with reference to FIGS.

[0048] The adaptive equalization unit 54 generates a real component XI(f) multiplied by the complex transfer function of the impulse response H1, an imaginary component XQ(f) multiplied by the complex transfer function of the impulse response H5, a real component YI(f) multiplied by the complex transfer function of the impulse response H9, and an imaginary component YI(f) multiplied by the complex transfer function of the impulse response H 13The adaptive equalization unit 54 then adds the resulting signal to the imaginary component YQ(f) multiplied by the complex transfer function YQ(f) to generate a sum signal. The sum signal generated by the adaptive equalization unit 54 is then folded in the frequency domain. The folding process involves folding back frequency components with absolute values ​​greater than half the symbol rate (the Nyquist frequency) in a line-symmetric manner around the Nyquist frequency and adding them together. This process corresponds to downsampling in the time domain.

[0049] The adaptive equalization unit 54 performs an M-point (M is a natural number, N≧M) inverse discrete Fourier transform or inverse fast Fourier transform on the sum signal that has undergone folding processing (corresponding to "M-IDFT" shown in FIG. 2). In this way, the adaptive equalization unit 54 converts the frequency domain signal into a time domain signal. Thereafter, the adaptive equalization unit 54 performs signal cutout processing on the time domain signal using the overlap-save method (corresponding to "Cut" shown in FIG. 2). In the following explanation, when using numerical values, an example will be given in which the number of M points is 128. When the number of M points is 128, 64 symbol points will be obtained by signal cutout processing using the overlap-save method.

[0050] In order to realize the above processing, the adaptive equalization unit 54 has a buffer, a Fourier transform unit, a branching unit, a coefficient calculation unit, an adding unit, a folding processing unit, an inverse Fourier transform unit, and a cutting unit.

[0051] Although the above describes a configuration in which folding, M-IDFT, and Cut processing are performed in that order, the processing may be performed in the order of M-IDFT, Cut, and downsampling.

[0052] The frequency / phase offset compensator 55 applies a frequency offset exp(jφ) to the added signal extracted by the adaptive equalizer 54 as described above. x (n)), where n represents the symbol interval.

[0053] The adaptive equalizer 54 outputs the real component inverted complex conjugate signal XI obtained by multiplying the complex transfer function of the impulse response H2.† (-f) and the complex transfer function of the impulse response H6 are multiplied to obtain the imaginary component inverted complex conjugate signal XQ. † (-f) and impulse response H 10 The real component inverted complex conjugate signal YI is multiplied by the complex transfer function † (-f) and impulse response H 14 The imaginary component inverted complex conjugate signal YQ is multiplied by the complex transfer function † (-f) to generate an added signal. The added signal generated by the adaptive equalization unit 54 is then subjected to folding, M-IDFT, and Cut processing.

[0054] The frequency / phase offset compensator 55 applies a frequency offset exp(-jφ) to the added signal extracted by the adaptive equalizer 54 as described above. x The frequency / phase offset compensation unit 55 multiplies the frequency offset exp(jφ x (n)) and the sum signal multiplied by the frequency offset exp(-jφ x (n)) to obtain the received signal of the X polarization component.

[0055] The demodulation digital signal processing unit 532 adds a transmission data bias correction signal C to the obtained reception signal of the X polarization component to cancel the bias deviation of the X polarization component. X The distortion-corrected received signal X of the X polarization component is obtained by adding (or subtracting) Rsig (n) is obtained.

[0056] On the other hand, the adaptive equalizer 54 outputs a real component XI(f) multiplied by the complex transfer function of the impulse response H3, an imaginary component XQ(f) multiplied by the complex transfer function of the impulse response H7, and an imaginary component XQ(f) multiplied by the complex transfer function of the impulse response H 11 The real component YI(f) multiplied by the complex transfer function of 15The sum signal generated by the adaptive equalization unit 54 is then subjected to folding, M-IDFT, and Cut processing. The frequency / phase offset compensation unit 55 applies a frequency offset exp(jφ y Multiply by (n).

[0057] The adaptive equalizer 54 outputs the real component inverted complex conjugate signal XI obtained by multiplying the complex transfer function of the impulse response H4. † (-f) and impulse response H 12 The imaginary component is inverted and the complex conjugate signal XQ is multiplied by the complex transfer function † (-f) and impulse response H 16 The real component inverted complex conjugate signal YI is multiplied by the complex transfer function † (-f) and impulse response H 14 The imaginary component inverted complex conjugate signal YQ is multiplied by the complex transfer function † (-f) to generate an added signal. The added signal generated by the adaptive equalization unit 54 is then subjected to folding, M-IDFT, and Cut processing.

[0058] The frequency / phase offset compensator 55 applies a frequency offset exp(-jφ) to the added signal extracted by the adaptive equalizer 54 as described above. y The frequency / phase offset compensation unit 55 multiplies the frequency offset exp(jφ y (n)) and the sum signal multiplied by the frequency offset exp(-jφ y (n)) is added to obtain the received signal of the Y polarization component.

[0059] The demodulation digital signal processing unit 532 adds a transmission data bias correction signal C to the obtained Y polarization component reception signal to cancel the bias deviation of the Y polarization component. Y The distortion-corrected received signal Y of the X polarization component is obtained by adding (or subtracting) Rsig (n) is obtained.

[0060] The values ​​of N, M, and impulse responses H1 to H 16 , and frequency offset exp(jφ x (n)), exp(-jφ x (n)), exp(jφ y (n)), exp(-jφ y (n)) are adaptively and dynamically changed. The receiver 50 obtains these values ​​by any method.

[0061] Next, the configuration and operation of the coefficient calculation unit will be described. FIGS. 3 to 6 are diagrams showing an example of the configuration of the coefficient calculation unit. As shown in FIGS. 3 to 6, the coefficient calculation unit provided in the demodulation digital signal processing unit 532 includes four coefficient calculation units. The coefficient calculation unit shown in FIG. 3 is a functional unit that calculates impulse responses H1, H3, H5, and H7. The coefficient calculation unit shown in FIG. 4 is a functional unit that calculates impulse responses H2, H4, H6, and H8. The coefficient calculation unit shown in FIG. 5 is a functional unit that calculates impulse responses H9, H10, H11, H12, H13, H14, H15, H16, H17, H18, H19, H20, H21, H22, H23, H24, H25, H26, H27, H28, H29, H30, H31, H32, H33, H34, H35, H36, H37, H38, H39, H40, H41, H42, H43, H44, H45, H46, H47, H48, H49, H50, H51, H52, H53, H54, H55, H56, H57, H58, H59, H60, H61, H62, H63, H64, H65, H66, H67, H68, H69, H70, H71, H72, H73, H74, H75, H76, H77, H78, H79, H80, H81, H81, H82, H83, H84, H85, H86, H87, H88, H89, H81, H81, H82, H83, H84, H85 11 ,H 13 ,H 15 The coefficient calculation unit shown in Figure 6 is a functional unit that calculates the impulse response H 10 ,H 12 ,H 14 ,H 16 The coefficient calculation unit includes a coefficient update unit that updates the value of the impulse response.

[0062] In the following description, the coefficient calculation unit shown in FIG. 3 will be referred to as the "first coefficient calculation unit," the coefficient calculation unit shown in FIG. 4 will be referred to as the "second coefficient calculation unit," the coefficient calculation unit shown in FIG. 5 will be referred to as the "third coefficient calculation unit," and the coefficient calculation unit shown in FIG. 6 will be referred to as the "fourth coefficient calculation unit." When there is no need to particularly distinguish between the first coefficient calculation unit to the fourth coefficient calculation unit, they will simply be referred to as the coefficient calculation unit. The operation of the coefficient calculation unit will be described below.

[0063] (Operation of the first coefficient calculation unit) The frequency domain signal of real component XI and the frequency domain signal of imaginary component XQ are input to the first coefficient calculation unit. The frequency domain signal of real component XI and the frequency domain signal of imaginary component XQ input to the first coefficient calculation unit are branched into a first path and a second path, respectively. In the first path, the frequency domain signal of real component XI and the frequency domain signal of imaginary component XQ are multiplied by a complex transfer function updated by the coefficient update unit.

[0064] In the second path, the frequency domain signal of the real component XI and the frequency domain signal of the imaginary component XQ are converted into inverted and complex conjugated frequency domain signals by the inversion and complex conjugation unit. As a result, the frequency domain signal of the real component XI input to the first coefficient calculation unit is converted into a real component inverted complex conjugate signal, and the frequency domain signal of the imaginary component XQ is converted into an imaginary component inverted complex conjugate signal.

[0065] In the first coefficient calculation unit, the real component inverted complex conjugate signal and the imaginary component inverted complex conjugate signal are multiplied by a signal based on the received signal. Here, the signal based on the received signal is a signal obtained based on the following processes (1) to (5).

[0066] (1): Reference signal (e.g., d x (n)) to the received signal (e.g., X Rsig Subtract (n) (2): A frequency offset (e.g., exp(-jφ)) is applied to the signal obtained by the process in (1). x Multiply (n))) (3): Add zeros to the signal obtained by processing (2) (corresponding to the "zero addition" shown in Figure 3) (4): The signal obtained in (3) is subjected to an M-point inverse discrete Fourier transform or inverse fast Fourier transform (corresponding to the "M-DFT" shown in Figure 3). (5): Copy the frequency domain signal obtained by processing (4) in the frequency domain (corresponding to the "foldback copy" shown in Figure 3).

[0067] A reference signal (e.g., d x (n) or d y(n)) is a pilot signal inserted in advance on the transmitting side or a received signal (e.g., X Rsig (n) or Y Rsig The value provisionally determined as (n) is used. The zero adding process shown in (3) is a process of adding M / N times the signal length to be cut in the overlap save method described in Reference 1 to the input signal. In the zero adding process, M / N times the signal length to be cut is continuously added to the input signal. The frequency domain copy shown in (5) is a process of copying the frequency domain signal axisymmetrically with the Nyquist frequency as the reference. The frequency domain copy shown in (5) corresponds to upsampling in the time domain.

[0068] Although the above describes a configuration in which zero addition, M-DFT, and foldback copying are performed, upsampling and N-DFT may be performed instead.

[0069] The real component inverted complex conjugate signal and the imaginary component inverted complex conjugate signal multiplied by the signal based on the received signal are input to a coefficient update unit. The coefficient update unit performs N-IDFT, Cut, zero addition, N-DFT, multiplication by a step size μ, and addition of the value of the previous impulse response on the real component inverted complex conjugate signal and the imaginary component inverted complex conjugate signal multiplied by the signal based on the received signal. As the step size μ, a normalized LMS (Reference 1) that normalizes the step size for each frequency bin by the input signal power may be used.

[0070] Taking the process of updating the impulse response H1 as an example of the process of the first coefficient calculation unit, the coefficient update unit first performs an N (e.g., N=256)-point inverse discrete Fourier transform or inverse fast Fourier transform on a real component inverted complex conjugate signal (here, signal A1) multiplied by a signal based on the received signal. This allows the coefficient update unit to convert the frequency-domain signal A1 into the time-domain signal A1. Next, the coefficient update unit performs a signal extraction process using the overlap-save method on the time-domain signal A1. Next, the coefficient update unit adds zeros to the extracted time-domain signal A1. Next, the coefficient update unit multiplies the zero-added time-domain signal A1 by a step size μ1. Next, the coefficient update unit updates the value of the impulse response H1 by adding the value of the previously obtained impulse response H1 to the time-domain signal A1 multiplied by the step size μ1.

[0071] The process of updating the impulse response H3 in the first coefficient calculation unit is the same as the process described above except that the step size value is different. Furthermore, the process of updating the impulse responses H5 and H7 in the first coefficient calculation unit is the same as the process described above except that the imaginary component inverted complex conjugate signal multiplied by the signal based on the received signal is input to the coefficient update unit and the step size value is different.

[0072] (Operation of the second coefficient calculation unit) The second coefficient calculation unit receives as input the inverted real component complex conjugate signal of the real component XI and the inverted imaginary component complex conjugate signal of the imaginary component XQ. The inverted real component complex conjugate signal of the real component XI and the inverted imaginary component complex conjugate signal of the imaginary component XQ input to the second coefficient calculation unit are branched into a first path and a second path, respectively. In the first path, the inverted real component complex conjugate signal of the real component XI and the inverted imaginary component complex conjugate signal of the imaginary component XQ are multiplied by the complex transfer function updated by the coefficient update unit.

[0073] In the second path, the inverted complex conjugate signal of the real component XI and the inverted complex conjugate signal of the imaginary component XQ are converted into inverted and complex conjugated frequency domain signals by the inverting and complex conjugating unit. As a result, the inverted complex conjugate signal of the real component XI input to the second coefficient calculation unit is converted into the frequency domain signal of the real component XI, and the inverted complex conjugate signal of the imaginary component XQ is converted into the frequency domain signal of the imaginary component XQ.

[0074] In the second coefficient calculation unit, the frequency signal of the real component XI and the frequency domain signal of the imaginary component XQ are multiplied by the signal based on the received signal. However, in the signal based on the received signal in the second coefficient calculation unit, a frequency offset exp(jφ x The signal obtained in the process of (1) is multiplied by (n). The frequency signal of real component XI and the frequency domain signal of imaginary component XQ multiplied by the signal based on the received signal are input to the coefficient update unit. The coefficient update unit performs N-IDFT, Cut, zero addition, N-DFT, multiplication by step size μ, and addition of the value of the previous impulse response on the frequency signal of real component XI and the frequency domain signal of imaginary component XQ multiplied by the signal based on the received signal. The process performed by the coefficient update unit is the same as the process described in Figure 3, so a description thereof will be omitted.

[0075] (Operation of the third coefficient calculation unit) The processing performed by the third coefficient calculation unit is based on the following points: the input signal is a Y-polarized signal, the step size used in the coefficient update unit is different, and in generating a signal based on the received signal, a frequency offset exp(jφ y (n)) is a reference signal (e.g., d y (n)) to the received signal (e.g., Y Rsig The process is the same as that performed by the first coefficient calculation unit, except that the signal obtained by subtracting (n) is multiplied by the coefficient.

[0076] (Operation of the fourth coefficient calculation unit) The processing performed by the fourth coefficient calculation unit is based on the following points: the input signal is a Y-polarized signal, the step size used in the coefficient update unit is different, and in generating a signal based on the received signal, a frequency offset exp(-jφ y (n)) is a reference signal (e.g., d y (n)) to the received signal (e.g., Y Rsig The process is the same as that performed by the second coefficient calculation unit, except that the signal obtained by subtracting (n) is multiplied by the coefficient.

[0077] The cut and zero-add processing in the coefficient update unit corresponds to multiplication of a rectangular window function in the time domain. By changing the window function in the time domain to a Cosine window and processing it as a convolution in the frequency domain, it is possible to omit the N-IDFT and N-DFT, and thus simplify the process.

[0078] The demodulation digital signal processing unit 532 configured as described above can perform convolution operations in the frequency domain, thereby reducing the amount of calculations required. As a result, it is possible to achieve power savings in the receiver of a digital coherent optical transmission system.

[0079] (Modification of the first embodiment) The demodulation digital signal processing unit 532 may be configured to perform signal processing such as compensation for frequency characteristics and chromatic dispersion. FIG. 7 is a diagram showing an example of the configuration of a demodulation digital signal processing unit 532a in a modification of the first embodiment. The demodulation digital signal processing unit 532a includes an adaptive equalization unit 54, a frequency / phase offset compensating unit 55, and a front-end correction and chromatic dispersion estimating unit 56. The demodulation digital signal processing unit 532a differs in configuration from the demodulation digital signal processing unit 532 in that the front-end correction and chromatic dispersion estimating unit 56 is provided between the functional units that make up the adaptive equalization unit 54. The rest of the configuration of the demodulation digital signal processing unit 532a is similar to that of the demodulation digital signal processing unit 532. The differences will be described below.

[0080] The front-end correction and chromatic dispersion estimation unit 56 multiplies the frequency domain signal by a receiving device characteristic and a chromatic dispersion compensation coefficient. For example, the front-end correction and chromatic dispersion estimation unit 56 multiplies the frequency domain signal of the real component XI by a receiving device characteristic H RXI and the chromatic dispersion compensation coefficient H CD Multiply by the receiving device characteristic H RXI and the chromatic dispersion compensation coefficient H CD The frequency domain signal of the real component XI multiplied by is split into four, two of which are input to the coefficient calculation unit as is, and the remaining two are inverted and converted into complex conjugate frequency domain signals and input to the coefficient calculation unit. The subsequent processing is the same as that described above.

[0081] Similarly, the front-end correction and chromatic dispersion estimation unit 56 calculates the receiver device characteristic H RXQ and the chromatic dispersion compensation coefficient H CD Multiply by the receiving device characteristic H RXQ and the chromatic dispersion compensation coefficient H CD The frequency domain signal of the imaginary component XQ multiplied by is split into four, two of which are input to the coefficient calculation unit as is, and the remaining two are inverted and converted to complex conjugate frequency domain signals and then input to the coefficient calculation unit.

[0082] Similarly, the front-end correction and chromatic dispersion estimation unit 56 calculates the receiver device characteristic H RYI and the chromatic dispersion compensation coefficient H CD Multiply by the receiving device characteristic H RYI and the chromatic dispersion compensation coefficient H CD The frequency domain signal of the real component YI multiplied by is split into four, two of which are input to the coefficient calculation unit as is, and the remaining two are inverted and converted to complex conjugate frequency domain signals and then input to the coefficient calculation unit.

[0083] Similarly, the front-end correction and chromatic dispersion estimation unit 56 calculates the receiver device characteristic H RYQ and the chromatic dispersion compensation coefficient H CD Multiply by the receiving device characteristic H RYQ and the chromatic dispersion compensation coefficient H CD The frequency domain signal of the imaginary component YQ multiplied by is split into four, two of which are input to the coefficient calculation unit as is, and the remaining two are inverted and converted to complex conjugate frequency domain signals and then input to the coefficient calculation unit.

[0084] In the front-end correction and chromatic dispersion estimation unit 56, a value obtained by multiplying the receiving-side device characteristics and a chromatic dispersion compensation coefficient in advance may be set, or frequency offset compensation may be performed by shifting the frequency bins of the main signal and the coefficient.

[0085] In the configuration of the demodulation digital signal processing unit 532a, the front-end correction and chromatic dispersion estimation unit 56 may be provided in a stage preceding the buffer.

[0086] The demodulation digital signal processing unit 532a may not include the front-end correction and chromatic dispersion estimation unit 56, but may perform frequency offset compensation by frequency-shifting the main signal and coefficients.

[0087] Figure 8 is a diagram showing the N-DFT size dependency (DFT is calculated using FFT) of the received SNR (Signal-to-Noise Ratio) for 128 GBaud, 256 QAM (Quadrature Amplitude Modulation) using the configuration shown in Figure 7. As shown in Figure 8, increasing the DFT size increases the time response (frequency resolution) that can be compensated, and therefore the received SNR (Signal-to-Noise Ratio) is improved.

[0088] Fig. 9 shows a comparison of the number of multiplications between a conventional configuration (e.g., the configuration described in Patent Document 1) and the configuration shown in Fig. 7. In Fig. 9, the input sampling rate is 256GSample / a, the symbol rate is 128GBaud, the DFT size is N, the IDFT block size is M=N / 2 (DFT and IDFT are assumed to be calculated using fast Fourier transform (FFT) and inverse fast Fourier transform (IFFT)), and the overlap amount of the Overlap Save method is 1 / 2 (in this case, the time response length within the compensable range is (N / 2 × sampling interval), resulting in the same compensation performance as a conventional configuration with a tap length of N / 2). Note that Fig. 9 does not include the amount of calculation for the imperfection coefficients and dispersion compensation coefficients of the receiving device (only the calculation of the adaptive filter coefficients is considered).

[0089] The number of multiplications in the fast Fourier transform is 4 × (N / 2) × log2(N), the number of multiplications in the inverse fast Fourier transform is 4 × (N / 4) × log2(N / 2), and the number of multiplications for the adaptive filter coefficients is 16 × N. Under these conditions, the number of symbols that can be output from one block is N / 4, so the number of multiplications per symbol is 2 × log2(N) + 4 × log2(N / 2) + 64. With the conventional configuration, it is sufficient to consider the number of multiplications for the convolution operation per symbol, so this is 16L where L is the number of taps in the adaptive filter.

[0090] (Second embodiment) In the second embodiment, a configuration that can reduce the number of discrete Fourier transforms or fast Fourier transforms compared to the first embodiment will be described. Note that in the second embodiment, the configuration of the adaptive equalization unit included in the demodulation digital signal processing unit is different from that of the first embodiment. Therefore, only the differences from the first embodiment will be described.

[0091] Fig. 10 is a diagram showing an example of the configuration of a demodulation digital signal processing unit 532b in the second embodiment. Note that Fig. 10 omits the configuration of the frequency / phase offset compensating unit 55 and subsequent units, which are similar to those in the first embodiment. The adaptive equalizer 54b of the demodulation digital signal processing unit 532b shown in Fig. 10 differs from the adaptive equalizer 54 in the configuration of the stage preceding the branching unit. Note that the demodulation digital signal processing unit 532b does not perform signal processing such as compensation for frequency characteristics and chromatic dispersion.

[0092] The adaptive equalizer 54b receives the real component XI and imaginary component XQ of the X-polarized received signal, which have been converted into digital signals by the ADCs 531-1 to 531-4, and the real component YI and imaginary component YQ of the Y-polarized received signal. The adaptive equalizer 54b multiplies the input imaginary component XQ by an imaginary unit j to generate an imaginary component jXQ. The adaptive equalizer 54b adds the real component XI and the imaginary component jXQ together. As a result, the adaptive equalizer 54b generates an addition signal of XI + jXQ. The adaptive equalizer 54b stores the generated addition signal in a buffer.

[0093] The adaptive equalizer 54b performs an N-point discrete Fourier transform or fast Fourier transform on the sum signal stored in the buffer (corresponding to "N-DFT" shown in FIG. 10), thereby converting the X-polarized sum signal into a frequency domain signal.

[0094] The frequency-domain sum signal generated by the adaptive equalization unit 54b is branched into two. One of the branched frequency-domain sum signals is converted into an inverted and complex-conjugated frequency-domain signal. In the following description, the frequency-domain sum signal that is branched and converted into an inverted and complex-conjugated frequency-domain signal before the branching unit is referred to as the "frequency-domain converted sum signal," and the frequency-domain sum signal that is not branched and converted into an inverted and complex-conjugated frequency-domain signal is referred to as the "frequency-domain unconverted sum signal."

[0095] The frequency domain pre-conversion sum signal and the frequency domain post-conversion sum signal are each branched into two, and the adaptive equalization unit 54b adds the frequency domain pre-conversion sum signal and the frequency domain post-conversion sum signal together, and then multiplies the result by 1 / 2. This signal is equivalent to the frequency domain signal of real component XI in the first embodiment. The sum signal multiplied by 1 / 2 (the frequency domain signal of real component XI) is then branched into four by a branching unit, and two of the four branched signals are input directly to the coefficient calculation unit, while the remaining two signals are inverted and converted into complex conjugate frequency domain signals and then input to the coefficient calculation unit.

[0096] Furthermore, the adaptive equalization unit 54b subtracts the frequency-domain post-conversion sum signal from the frequency-domain pre-conversion sum signal, and then multiplies the result by 1 / 2j. This signal is equivalent to the frequency-domain signal of the imaginary component XQ in the first embodiment. Thereafter, the signal multiplied by 1 / 2j (the frequency-domain signal of the imaginary component XQ) is branched into four by a branching unit, and two of the four branched signals are input directly to the coefficient calculation unit, while the remaining two signals are inverted and converted into complex conjugate frequency-domain signals and then input to the coefficient calculation unit. The above is the processing related to the X polarization.

[0097] Similarly, the adaptive equalization unit 54b multiplies the input imaginary component YQ by the imaginary unit j to generate an imaginary component jYQ. The adaptive equalization unit 54b adds the real component YI and the imaginary component jYQ. As a result, the adaptive equalization unit 54b generates an added signal of YI + jYQ. The adaptive equalization unit 54b stores the generated added signal in a buffer.

[0098] The adaptive equalizer 54b performs an N-point discrete Fourier transform or fast Fourier transform on the sum signal stored in the buffer (corresponding to "N-DFT" shown in FIG. 10), thereby converting the Y-polarized sum signal into a frequency domain signal.

[0099] The frequency-domain sum signal generated by the adaptive equalization unit 54b is branched into two. One of the branched frequency-domain sum signals is converted into an inverted and complex-conjugated frequency-domain signal. The frequency-domain pre-conversion sum signal and the frequency-domain post-conversion sum signal are each branched into two, and the adaptive equalization unit 54b adds the frequency-domain pre-conversion sum signal and the frequency-domain post-conversion sum signal together and then multiplies the result by 1 / 2. This signal is equivalent to the frequency-domain signal of the real component YI in the first embodiment. The sum signal multiplied by 1 / 2 (the frequency-domain signal of the real component YI) is then branched into four by the branching unit, and two of the four branched signals are input directly to the coefficient calculation unit, while the remaining two signals are converted into frequency-domain signals that are inverted and complex-conjugated and then input to the coefficient calculation unit.

[0100] Furthermore, the adaptive equalization unit 54b subtracts the frequency-domain post-conversion addition signal from the frequency-domain pre-conversion addition signal, and then multiplies the result by 1 / 2j. This signal is equivalent to the frequency-domain signal of the imaginary component YQ in the first embodiment. Thereafter, the signal multiplied by 1 / 2j (the frequency-domain signal of the imaginary component YQ) is branched into four by a branching unit, and two of the four branched signals are input directly to the coefficient calculation unit, while the remaining two signals are inverted and converted into complex conjugate frequency-domain signals and then input to the coefficient calculation unit. The above is the processing related to Y polarization.

[0101] In the adaptive equalization unit 54b, the processing after the coefficient calculation unit is the same as in the first embodiment.

[0102] The demodulation digital signal processing unit 532b in the second embodiment configured as described above can reduce the number of discrete Fourier transforms or fast Fourier transforms compared to the first embodiment. Specifically, the demodulation digital signal processing unit 532 in the second embodiment performs a discrete Fourier transform or fast Fourier transform after adding the real component XI and the imaginary component XQ. This eliminates the need to perform a discrete Fourier transform or fast Fourier transform on each of the real component XI and the imaginary component XQ. Therefore, the number of discrete Fourier transforms or fast Fourier transforms can be reduced compared to the first embodiment.

[0103] (Modification of the second embodiment) The adaptive equalizer 54b may have a configuration for performing signal processing such as compensation for frequency characteristics and chromatic dispersion, as in the first embodiment. Fig. 11 is a diagram showing an example of the configuration of a demodulation digital signal processor 532c in a modified example of the second embodiment. Note that Fig. 11 omits the configuration from the frequency / phase offset compensator 55 onwards, which has the same configuration as in the first embodiment.

[0104] The demodulation digital signal processing unit 532c includes an adaptive equalization unit 54b, a frequency / phase offset compensation unit 55 (omitted in FIG. 11), and a front-end correction and chromatic dispersion estimation unit 56. The demodulation digital signal processing unit 532c differs from the configuration shown in FIG. 10 in that the front-end correction and chromatic dispersion estimation unit 56 is provided between the functional units that make up the adaptive equalization unit 54b. The rest of the configuration of the demodulation digital signal processing unit 532c is the same as the configuration shown in FIG. 10. The differences will be explained below.

[0105] The front-end correction and chromatic dispersion estimation unit 56 calculates the receiver device characteristic H RXI and the chromatic dispersion compensation coefficient H CD Multiply by the receiving device characteristic H RXI and the chromatic dispersion compensation coefficient H CDThe frequency domain signal of the real component XI multiplied by is branched into four by the branching unit, two of which are input directly to the coefficient calculation unit, and the remaining two are inverted and converted into complex conjugate frequency domain signals and then input to the coefficient calculation unit. The subsequent processing is the same as that described above.

[0106] Similarly, the front-end correction and chromatic dispersion estimation unit 56 calculates the receiving device characteristic H RXQ and the chromatic dispersion compensation coefficient H CD Multiply by the receiving device characteristic H RXQ and the chromatic dispersion compensation coefficient H CD The frequency domain signal of the imaginary component XQ multiplied by is branched into four by the branching unit, and two of the four branched signals are input directly to the coefficient calculation unit, while the remaining two signals are inverted and converted into complex conjugate frequency domain signals and then input to the coefficient calculation unit.

[0107] Similarly, the front-end correction and chromatic dispersion estimation unit 56 calculates the receiver device characteristic H RYI and the chromatic dispersion compensation coefficient H CD Multiply by the receiving device characteristic H RYI and the chromatic dispersion compensation coefficient H CD The frequency domain signal of the real component YI multiplied by is branched into four by the branching unit, and two of the four branched signals are input directly to the coefficient calculation unit, while the remaining two signals are inverted and converted into complex conjugate frequency domain signals and then input to the coefficient calculation unit.

[0108] Similarly, the front-end correction and chromatic dispersion estimation unit 56 calculates the receiving device characteristic H RYQ and the chromatic dispersion compensation coefficient H CD Multiply by the receiving device characteristic H RYQ and the chromatic dispersion compensation coefficient H CDThe frequency domain signal of the imaginary component YQ multiplied by is branched into four by the branching unit, and two of the four branched signals are input directly to the coefficient calculation unit, while the remaining two signals are inverted and converted into complex conjugate frequency domain signals and then input to the coefficient calculation unit.

[0109] The demodulation digital signal processing unit 532c may not include the front-end correction and chromatic dispersion estimation unit 56, but may perform frequency offset compensation by frequency-shifting the main signal and coefficients.

[0110] (Third embodiment) In the third embodiment, the configuration of the adaptive equalizer included in the demodulation digital signal processing unit is different from that of the second embodiment, so the differences from the second embodiment will be described.

[0111] Fig. 12 is a diagram showing an example of the configuration of a demodulation digital signal processing unit 532d in the third embodiment. Note that Fig. 12 omits the configuration from the frequency / phase offset compensating unit 55 onwards, which has the same configuration as in the second embodiment (particularly the demodulation digital signal processing unit 532c shown in Fig. 11). The demodulation digital signal processing unit 532d includes an adaptive equalizer 54d and a frequency / phase offset compensating unit 55 (omitted in Fig. 12).

[0112] The adaptive equalizer 54d applies the pre-conversion sum signal in the frequency domain of the X-polarized wave to the receiving device characteristic H RXI and the receiving device characteristics H RXQ The sum of (1 / 2 x H CD * )

[0113] Similarly, the adaptive equalizer 54d applies the reception-side device characteristic H RXI From the receiving device characteristics H RXQ The value obtained by subtracting (1 / 2 × H CD * ) multiply by 1 / 2×H CD *The frequency domain pre-conversion sum signal of X polarization multiplied by 1 / 2×H CD * Each of the converted sum signals of the X-polarized frequency domain multiplied by is split into two.

[0114] The adaptive equalizer 54d is 1 / 2×H CD * The pre-conversion sum signal of the X polarization frequency domain multiplied by 1 / 2×H CD * The converted sum signal of the frequency domain of the X-polarized wave multiplied by is then added to this. This sum signal is then split into four by a splitter, and two of the four split signals are input directly to the coefficient calculator, while the remaining two signals are inverted and converted into complex conjugate frequency domain signals and then input to the coefficient calculator.

[0115] Furthermore, the adaptive equalizer 54d is 1 / 2×H CD * The converted sum signal of the X polarization frequency domain multiplied by 1 / 2×H CD * The pre-conversion sum signal in the frequency domain of the X-polarized wave multiplied by is subtracted from the pre-conversion sum signal in the frequency domain of the X-polarized wave. Then, this subtracted signal is split into four by a splitter, and two of the four split signals are input directly to the coefficient calculator, while the remaining two signals are inverted and converted into complex conjugate frequency domain signals and then input to the coefficient calculator. The above is the processing related to the X polarization.

[0116] The adaptive equalizer 54d applies the receiving device characteristic H RYI and the receiving device characteristics H RYQ The sum of (1 / 2 x H CD * Similarly, the adaptive equalizer 54d multiplies the converted sum signal in the frequency domain of the Y-polarized wave by the receiving device characteristic H RYI From the receiving device characteristics H RYQ The value obtained by subtracting (1 / 2 × H CD * ) multiply by 1 / 2×H CD *The frequency domain pre-conversion sum signal of Y polarization multiplied by 1 / 2×H CD * Each of the converted sum signals of the Y-polarized wave in the frequency domain multiplied by is split into two.

[0117] The adaptive equalizer 54d is 1 / 2×H CD * The frequency domain pre-conversion sum signal of Y polarization multiplied by 1 / 2×H CD * The frequency domain converted sum signal of the Y-polarized wave multiplied by is then added to this. This sum signal is then split into four by a splitter, and two of the four split signals are input directly to the coefficient calculator, while the remaining two signals are inverted and converted into complex conjugate frequency domain signals, which are then input to the coefficient calculator.

[0118] Furthermore, the adaptive equalizer 54d is 1 / 2×H CD * The sum signal after the frequency domain transformation of the Y polarization multiplied by is 1 / 2×H CD * The pre-conversion sum signal in the frequency domain of the Y-polarized wave multiplied by is subtracted from the signal. The subtracted signal is then split into four by a splitter, and two of the four split signals are input directly to the coefficient calculator, while the remaining two signals are inverted and converted to complex conjugate frequency domain signals and then input to the coefficient calculator. The above is the processing related to Y polarization.

[0119] In the adaptive equalization unit 54d, the processing after the coefficient calculation unit is the same as in the second embodiment.

[0120] The demodulation digital signal processing unit 532d in the third embodiment configured as described above is different from the second embodiment in that it is possible to reduce the number of discrete Fourier transforms or fast Fourier transforms compared to the first embodiment. RXI -H RXQ ,H RYI -H RYQIf is small, it is possible to reduce the bit precision.

[0121] (Modification of the third embodiment) The demodulation digital signal processing unit 532d may perform frequency offset compensation by frequency-shifting the main signal and coefficients at a stage after the N-DFT and before the branching unit.

[0122] (Modifications common to the first to third embodiments) In each of the above embodiments, a configuration for performing wavelength division multiplexing in addition to polarization division multiplexing may be combined. In such a configuration, the digital coherent optical transmission system 1 shown in FIG. 1 differs from the system in the following respects. The transmitter 10 further includes transmitting units 100 equal to the number of WDM (Wavelength Division Multiplexing) channels. For example, if the number of WDM channels is 10, the transmitter 10 will include 10 transmitting units 100. Each transmitting unit 100 outputs an optical signal of a different wavelength. A WDM multiplexer, an optical fiber transmission line 30, and a WDM demultiplexer are provided between the transmitter 10 and the receiver 50. The WDM multiplexer multiplexes the optical signals output by each transmitting unit 100 and outputs the multiplexed optical signal to the optical fiber transmission line 30. The WDM demultiplexer demultiplexes the optical signals transmitted through the optical fiber transmission line 30 by wavelength. The receiver 50 further includes receiving units 500 equal to the number of WDM channels. For example, if the number of WDM channels is 10, the receiver 50 will include 10 receiving units 500. Each receiving unit 500 receives the optical signal demultiplexed by the WDM demultiplexer 40. The wavelengths of the optical signals received by each receiver 500 are different from each other. The processing executed by the receiver 500 is the same as the processing described above.

[0123] In each of the above embodiments, when N=M, folding processing in the adaptive equalizers 54, 54b, and 54d does not need to be performed.

[0124] Some of the functional units of the receiver 50 in the above-described embodiment may be implemented by a computer. In this case, a program for implementing the functions may be recorded on a computer-readable recording medium, and the program recorded on the recording medium may be read into a computer system and executed. Note that the term "computer system" here includes hardware such as an OS and peripheral devices.

[0125] Furthermore, "computer-readable recording media" refers to portable media such as flexible disks, optical magnetic disks, ROMs (Read Only Memory), and CD-ROMs, as well as storage devices such as hard disks built into computer systems. Furthermore, "computer-readable recording media" may also include devices that dynamically store programs for a short period of time, such as communication lines used when transmitting programs over networks like the Internet or communication lines like telephone lines, or devices that store programs for a fixed period of time, such as volatile memory within a computer system that serves as a server or client in such cases. Furthermore, the program may be one that implements some of the aforementioned functions, or one that can realize the aforementioned functions in combination with a program already stored in the computer system, or one that can be implemented using a programmable logic device such as an FPGA (Field-Programmable Gate Array).

[0126] Although an embodiment of the present invention has been described above in detail with reference to the drawings, the specific configuration is not limited to this embodiment, and includes designs within the scope of the gist of the present invention. [Industrial Applicability]

[0127] The present invention can be applied to a technique for receiving a single-carrier polarization multiplexed signal in digital coherent optical transmission. [Explanation of symbols]

[0128] 1...Digital coherent optical transmission system 10...Transmitter 30...Optical fiber transmission line 31...Optical amplifier 50...Receiver 54, 54b, 54d...Adaptive equalization section 55...Frequency / phase offset compensation section 56...Front-end correction and chromatic dispersion estimation unit 100...Transmitter 110...Digital signal processing unit 111...encoding section 112...Mapping section 113...Training signal insertion unit 114...Frequency conversion unit 115...Waveform shaping section 116...Pre-equalization section 117-1 to 117-4...Digital-to-analog converters 120...Modulator driver 121-1~121-4...Amplifier 130...Light source 140...Integrated module 141-1, 141-2...IQ modulator 142...Polarization combining unit 500...Receiver 510...Local oscillator light source 520...Optical front end 521...Polarization separation unit 522-1, 522-2... Optical 90-degree hybrid coupler 523-1~523-4...BPD 524-1~524-4...Amplifier 530...Digital signal processing unit 531-1 to 531-4...Analog-to-digital converters 532, 532a, 532b, 532c, 532d...Demodulation digital signal processing section 533…Demapping section 534...Decoding unit

Claims

1. a conversion step of converting the real and imaginary components of each polarization of the polarization-multiplexed received signal into frequency domain signals; a signal input step of inputting, as input signals, the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization, and frequency domain signals obtained by performing frequency inversion on the frequency axis and complex conjugate transformation on the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization; an equalization step of performing a first equalization process for each polarization, in which the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization included in the input signal are multiplied by a complex transfer function, and then added together to inversely convert the frequency domain signals into time domain signals; and a second equalization process for multiplying the frequency domain signals of the real components of each polarization included in the input signal after conversion by a complex transfer function, and then added together to inversely convert the frequency domain signals into time domain signals. a compensation step of performing, for each polarization, a phase rotation for frequency offset compensation on the time domain signal converted by the first equalization processing to generate a first sum signal, performing a phase rotation opposite to the phase rotation for frequency offset compensation on the time domain signal converted by the second equalization processing to generate a second sum signal, and adding or subtracting a transmit data bias correction signal to or from a signal obtained by adding the first sum signal and the second sum signal; A signal processing method comprising:

2. an addition processing step of performing an imaginary unit multiplication process of multiplying an imaginary component of each polarization of the polarization-multiplexed received signal by an imaginary unit j, and then performing an addition processing of adding the imaginary component multiplied by the imaginary unit j and a real component of each polarization of the polarization-multiplexed received signal; a transform step of transforming the signal obtained by adding the imaginary component multiplied by the imaginary unit j and the real component into a frequency domain signal; a signal input step of inputting as input signals: a first computed frequency domain signal obtained by adding the frequency domain signal of each polarization to a converted frequency domain signal obtained by frequency inverting the frequency domain signal on the frequency axis and taking a complex conjugate, and then multiplying by ½; or a first computed frequency domain signal obtained by adding the frequency domain signal after compensation has been performed on the frequency domain signal to a converted frequency domain signal after compensation has been performed on the converted frequency domain signal, and then multiplying by ½; a second computed frequency domain signal obtained by subtracting from the frequency domain signal a converted frequency domain signal obtained by frequency inverting the frequency domain signal on the frequency axis and taking a complex conjugate, and then multiplying by ½; or a second computed frequency domain signal obtained by subtracting from the compensated frequency domain signal the converted frequency domain signal a compensated frequency domain signal, and then multiplying by ½; an equalization step for performing, for each polarization, a first equalization process in which a first calculated frequency domain signal of the real component and a first calculated frequency domain signal of the imaginary component of each polarization included in the input signal are respectively multiplied by a complex transfer function, and then added together to inversely convert from a frequency domain signal to a time domain signal; and a second equalization process in which a second calculated frequency domain signal of the real component and a second calculated frequency domain signal of the imaginary component of each polarization included in the input signal are respectively multiplied by a complex transfer function, and then added together to inversely convert from a frequency domain signal to a time domain signal; a compensation step of performing, for each polarization, a phase rotation for frequency offset compensation on the time domain signal converted by the first equalization processing to generate a first sum signal, performing a phase rotation opposite to the phase rotation for frequency offset compensation on the time domain signal converted by the second equalization processing to generate a second sum signal, and adding or subtracting a transmit data bias correction signal to or from a signal obtained by adding the first sum signal and the second sum signal; A signal processing method comprising:

3. In the signal input step, the frequency domain signal is branched into a first path and a second path, and a first signal processing is performed for each polarization, in which the frequency domain signal branched into the first path is added to a frequency domain signal branched into the second path that has been frequency-inverted and complex-conjugated, and the resulting signal is multiplied by 1 / 2; and a second signal processing is performed for each polarization, in which the frequency domain signal branched into the first path is subtracted from the frequency domain signal branched into the second path that has been frequency-inverted and complex-conjugated, and the resulting signal is multiplied by 1 / 2j.

3. The signal processing method according to claim 2.

4. In the signal input step, the frequency domain signal is branched into a first path and a second path, and a first signal processing is performed for each polarization, the first signal processing being a frequency domain signal branched into the first path after compensation for frequency characteristics and chromatic dispersion compensation has been performed, and a second signal processing being a frequency domain signal branched into the second path after frequency inversion and complex conjugation has been performed, compensation for frequency characteristics and chromatic dispersion compensation has been performed, is added; and the second signal processing being a frequency domain signal branched into the second path after frequency inversion and complex conjugation has been performed, compensation for frequency characteristics and chromatic dispersion compensation has been performed, is subtracted from the frequency domain signal branched into the first path after compensation for frequency characteristics and chromatic dispersion compensation has been performed, and the input signal is then generated.

3. The signal processing method according to claim 2.

5. a frequency conversion unit that converts the real and imaginary components of each polarization of the polarization-multiplexed received signal into frequency domain signals; a signal input unit that receives, as input signals, the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization, and frequency domain signals obtained by performing frequency inversion on the frequency axis and complex conjugation on the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization; an equalization unit that performs, for each polarization, a first equalization process in which the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization included in the input signal are multiplied by a complex transfer function, and then summed, thereby performing an inverse conversion from the frequency domain signals to time domain signals; and a second equalization process in which the frequency domain signals of the real components and the frequency domain signals of the imaginary components of each polarization included in the input signal are multiplied by a complex transfer function, and then summed, thereby performing an inverse conversion from the frequency domain signals to time domain signals; a compensator that performs, for each polarization, a phase rotation for frequency offset compensation on the time domain signal converted by the first equalization processing to generate a first sum signal, performs a phase rotation opposite to the phase rotation for frequency offset compensation on the time domain signal converted by the second equalization processing to generate a second sum signal, and adds or subtracts a transmit data bias correction signal to or from a signal obtained by adding the first sum signal and the second sum signal; A signal processing device comprising:

6. an adder that performs an imaginary unit multiplication process of multiplying an imaginary component of each polarization of the polarization-multiplexed received signal by an imaginary unit j, and then performs an addition process of adding the imaginary component multiplied by the imaginary unit j and a real component of each polarization of the polarization-multiplexed received signal; a frequency transform unit that transforms the signal obtained by adding the imaginary component multiplied by the imaginary unit j and the real component into a frequency domain signal; a signal input unit that receives as input signals: a first computed frequency domain signal obtained by adding the frequency domain signal of each polarization to a converted frequency domain signal obtained by frequency inverting the frequency domain signal on the frequency axis and taking a complex conjugate, and then multiplying by ½; or a first computed frequency domain signal obtained by adding the frequency domain signal after compensation has been performed on the frequency domain signal to a converted frequency domain signal after compensation has been performed on the converted frequency domain signal, and then multiplying by ½; a second computed frequency domain signal obtained by subtracting from the frequency domain signal a converted frequency domain signal obtained by frequency inverting the frequency domain signal on the frequency axis and taking a complex conjugate, and then multiplying by ½; or a second computed frequency domain signal obtained by subtracting from the compensated frequency domain signal the converted frequency domain signal a compensated frequency domain signal, and then multiplying by ½; an equalization unit that performs, for each polarization, a first equalization process in which a first calculated frequency domain signal of the real component and a first calculated frequency domain signal of the imaginary component of each polarization included in the input signal are multiplied by a complex transfer function, and then summed to inversely convert the frequency domain signal into a time domain signal; and a second equalization process in which a second calculated frequency domain signal of the real component and a second calculated frequency domain signal of the imaginary component of each polarization included in the input signal are multiplied by a complex transfer function, and then summed to inversely convert the frequency domain signal into a time domain signal; a compensator that performs, for each polarization, a phase rotation for frequency offset compensation on the time domain signal converted by the first equalization processing to generate a first sum signal, performs a phase rotation opposite to the phase rotation for frequency offset compensation on the time domain signal converted by the second equalization processing to generate a second sum signal, and adds or subtracts a transmit data bias correction signal to or from a signal obtained by adding the first sum signal and the second sum signal; A signal processing device comprising:

7. 7. A communication system comprising: a transmitter for transmitting a polarization multiplexed signal; and a receiver having the signal processing device according to claim 5.

Images