CHARACTERISTIC MEASUREMENT APPARATUS, CHARACTERISTIC MEASUREMENT METHOD AND COMPUTER PROGRAM

Abstract

A characteristic measurement device that measures a characteristic between lanes of a transmitter and a receiver connected via an optical fiber transmission path, the characteristic measurement device including an adaptive equalization unit that uses a polarization multiplexed reception signal and a phase conjugate signal of the polarization multiplexed reception signal or a plurality of signals mathematically equivalent to the reception signal and the phase conjugate signal of the reception signal as input signals and performs equalization processing on the input signals, and a characteristic function derivation unit that calculates a first inverse characteristic representing an inverse characteristic of the transmitter and a second inverse characteristic representing an inverse characteristic of the receiver on the basis of a filter coefficient obtained at the time of the equalization processing performed by the adaptive equalization unit and a frequency offset.

Claims

1. A characteristic measurement device that measures a characteristic between lanes of a transmitter and a receiver connected via an optical fiber transmission path, the characteristic measurement device comprising: an adaptive equalizer configured to use a polarization multiplexed reception signal and a phase conjugate signal of the polarization multiplexed reception signal or a plurality of signals mathematically equivalent to the reception signal and the phase conjugate signal of the reception signal as input signals and performs equalization processing on the input signals; and a characteristic function deriver configured to calculate a first inverse characteristic representing an inverse characteristic of the transmitter and a second inverse characteristic representing an inverse characteristic of the receiver on the basis of a filter coefficient obtained at a time of the equalization processing performed by the adaptive equalizer and a frequency offset.

2. The characteristic measurement device according to claim 1, wherein the characteristic function deriver calculates the first inverse characteristic and the second inverse characteristic including IQ crosstalk across a plurality of spatially multiplexed signals.

3. The characteristic measurement device according to claim 1, wherein the characteristic function deriver calculates the first inverse characteristic and the second inverse characteristic including IQ crosstalk across a plurality of wavelength multiplexed signals.

4. The characteristic measurement device according to claim 1, wherein the plurality of signals mathematically equivalent to the reception signal and the phase conjugate signal of the reception signal are signals obtained by performing wavelength dispersion compensation on each of a complex signal and a phase conjugate signal of the complex signal, and phase conjugate signals of the signals on which the wavelength dispersion compensation is performed.

5. The characteristic measurement device according to claim 1, having a function of multiplying the first inverse characteristic and the second inverse characteristic by a matrix representing a degree of freedom of polarization rotation and phase rotation to transform the first inverse characteristic and the second inverse characteristic into a format capable of being handled by a combination of filters that compensate for a characteristic between some lanes of the transmitter and the receiver.

6. The characteristic measurement device according to claim 1, having a function of causing a forward characteristic obtained from the first inverse characteristic and the second inverse characteristic to act on a known signal and determining a filter coefficient so as to cancel the forward characteristic to transform the forward characteristic into a format capable of being handled by a combination of filters that compensate for a characteristic between some lanes of the transmitter and the receiver.

7. An optical transmission system comprising: the characteristic measurement device according to claim 1; and the transmitter and the receiver having a filter function of compensating for signal waveform distortion on the basis of the first inverse characteristic and the second inverse characteristic obtained by the characteristic measurement device.

8. A characteristic measurement method performed by a characteristic measurement device that measures a characteristic between lanes of a transmitter and a receiver connected via an optical fiber transmission path, the characteristic measurement method comprising: using a polarization multiplexed reception signal and a phase conjugate signal of the polarization multiplexed reception signal or a plurality of signals mathematically equivalent to the reception signal and the phase conjugate signal of the reception signal as input signals and performing equalization processing on the input signals; and calculating a first inverse characteristic representing an inverse characteristic of the transmitter and a second inverse characteristic representing an inverse characteristic of the receiver on the basis of a filter coefficient obtained at a time of the equalization processing and a frequency offset.

9. A non-transitory storage medium that stores a program for making a computer perform processes as a characteristic measurement device that measures a characteristic between lanes of a transmitter and a receiver connected via an optical fiber transmission path, the processes comprising: using a polarization multiplexed reception signal and a phase conjugate signal of the polarization multiplexed reception signal or a plurality of signals mathematically equivalent to the reception signal and the phase conjugate signal of the reception signal and performing equalization processing on the input signals; and calculating a first inverse characteristic representing an inverse characteristic of the transmitter and a second inverse characteristic representing an inverse characteristic of the receiver on the basis of a filter coefficient obtained at a time of the equalization processing performed and a frequency offset.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0018] FIG. 1 is a diagram illustrating a configuration example of a digital coherent optical transmission system in a first embodiment.

[0019] FIG. 2 is a diagram for explaining an outline of processing for deriving inverse characteristics of a transmitter and a receiver in the first embodiment.

[0020] FIG. 3 is a diagram illustrating a configuration example of a demodulation digital signal processing unit including an adaptive equalization unit in the first embodiment.

[0021] FIG. 4 is a diagram illustrating a configuration example of a characteristic function derivation unit in the first embodiment.

[0022] FIG. 5 is a flowchart illustrating a flow of processing of a receiver in the first embodiment.

[0023] FIG. 6 is a diagram for explaining an outline of processing for deriving inverse characteristics of a transmitter and a receiver in a second embodiment.

[0024] FIG. 7 is a diagram illustrating a configuration example of a demodulation digital signal processing unit including an adaptive equalization unit in the second embodiment.

[0025] FIG. 8 is a diagram illustrating a configuration example of a characteristic function derivation unit in the second embodiment.

[0026] FIG. 9 is a diagram illustrating a configuration example of a digital coherent optical transmission system 1 in a third embodiment.

[0027] FIG. 10 is a diagram illustrating a (first) example in which a plurality of filters of a 22 format is combined in the third embodiment.

[0028] FIG. 11 is a diagram illustrating a (second) example in which a plurality of filters of a 22 format is combined in the third embodiment.

[0029] FIG. 12 is a diagram illustrating an example of a method of determining a value to be input to a filter in the third embodiment.

DESCRIPTION OF EMBODIMENTS

[0030] Hereinafter, an embodiment of the present invention will be described with reference to the drawings.

First Embodiment

[0031] FIG. 1 is a diagram illustrating a configuration example of a digital coherent optical transmission system 1 in a first embodiment. The digital coherent optical transmission system 1 includes a transmitter 10 and a receiver 50. The receiver 50 receives a polarization multiplexed signal from the transmitter 10.

[0032] The transmitter 10 includes one or more transmission units 100. The transmission unit 100 outputs an optical signal of a specified wavelength to an optical fiber transmission path 30. The optical fiber transmission path 30 includes any number of optical amplifiers 31. Each optical amplifier 31 receives an input of the optical signal from the optical fiber transmission path 30 on the transmitter 10 side to amplify, and outputs the same to the optical fiber transmission path 30 on the receiver 50 side. The receiver 50 includes one or more reception units 500. The reception unit 500 receives the optical signal.

[0033] The transmission unit 100 includes a digital signal processing unit 110, a modulator driver 120, a light source 130, and an integration module 140. The digital signal processing unit 110 includes an encoding unit 111, a mapping unit 112, a training signal insertion unit 113, a frequency change unit 114, a waveform shaping unit 115, a pre-equalization unit 116, and digital-to-analog converters (DACs) 117-1 to 117-4.

[0034] The encoding unit 111 outputs a transmission signal obtained by performing forward error correction (FEC) encoding on a transmission bit string. The mapping unit 112 maps the transmission signal output from the encoding unit 111 to a symbol. The training signal insertion unit 113 inserts a known training signal into the transmission signal mapped to the symbol by the mapping unit 112. The frequency change unit 114 performs up-sampling by changing a sampling frequency for the transmission signal into which the training signal is inserted. The waveform shaping unit 115 limits a band of the transmission signal subjected to sampling.

[0035] The pre-equalization unit 116 compensates for distortion of a waveform of the transmission signal, the band of which is limited by the waveform shaping unit 115, and outputs the same to the DACs 117-1 to 117-4. The DAC 117-1 converts an I (in-phase) component of an X-polarized wave of the transmission signal input from the pre-equalization unit 116 from a digital signal to an analog signal, and outputs the same to the modulator driver 120. The DAC 117-2 converts a Q (orthogonal) component of the X-polarized wave of the transmission signal input from the pre-equalization unit 116 from a digital signal to an analog signal, and outputs the same to the modulator driver 120. The DAC 117-3 converts an I component of a Y-polarized wave of the transmission signal input from the pre-equalization unit 116 from a digital signal to an analog signal, and outputs the same to the modulator driver 120. The DAC 117-4 converts a Q component of the Y-polarized wave of the transmission signal input from the pre-equalization unit 116 from a digital signal to an analog signal, and outputs the same to the modulator driver 120.

[0036] The modulator driver 120 includes amplifiers 121-1 to 121-4. The amplifier 121-i (i is an integer of 1 or larger and 4 or smaller) amplifies the analog signal output from the DAC 117-i, and drives a modulator of the integration module 140 by the amplified analog signal. The light source 130 is, for example, an LD (semiconductor laser). The light source 130 outputs light of a specified wavelength.

[0037] The integration module 140 includes IQ modulators 141-1 and 141-2 and a polarization synthesis unit 142. The IQ modulator 141-1 outputs an optical signal of the X-polarized wave generated by modulating the optical signal output by the light source 130 with the I component of the X-polarized wave output from the amplifier 121-1 and the Q component of the X-polarized wave output from the amplifier 121-2. The IQ modulator 141-2 outputs an optical signal of the Y-polarized wave generated by modulating the optical signal output by the light source 130 with the I component of the Y-polarized wave output from the amplifier 121-3 and the Q component of the Y-polarized wave output from the amplifier 121-4. The polarization synthesis unit 142 performs polarization multiplexing on the optical signal of the X-polarized wave output by the IQ modulator 141-1 and the optical signal of the Y-polarized wave output by the IQ modulator 141-2 to output.

[0038] The reception 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 oscillator (LO).

[0039] The optical front end 520 converts an optical signal into an electric signal while keeping a phase and amplitude of a phase modulation signal subjected to polarization multiplexing. The optical front end 520 includes a polarization separation unit 521, optical 90-degree hybrid couplers 522-1 and 522-2, balanced photo diodes (BPDs) 523-1 to 523-4, and amplifiers 524-1 to 524-4.

[0040] The polarization separation unit 521 separates the input optical signal into the X-polarized wave and the Y-polarized wave. The polarization separation unit 521 outputs the optical signal of the X-polarized wave to the optical 90-degree hybrid coupler 522-1, and outputs the optical signal of the Y-polarized wave to the optical 90-degree hybrid coupler 522-2.

[0041] The optical 90-degree hybrid coupler 522-1 causes the optical signal of the X-polarized wave and the local oscillator output from the local oscillation light source 510 to interfere with each other, and extracts an I component and a Q component of a reception photoelectric field. The optical 90-degree hybrid coupler 522-1 outputs the extracted I component and Q component of the X-polarized wave to the BPDs 523-1 and 523-2, respectively.

[0042] The optical 90-degree hybrid coupler 522-2 causes the optical signal of the Y-polarized wave and the local oscillator output from the local oscillation light source 510 to interfere with each other, and extracts an I component and a Q component of a reception photoelectric field. The optical 90-degree hybrid coupler 522-2 outputs the extracted I component and Q component of the Y-polarized wave to the BPDs 523-3 and BPD 523-4, respectively.

[0043] The BPDs 523-1 to 523-4 are differential input type photoelectric converters. The BPD 523-i outputs, to the amplifier 524-i, a difference value between photocurrents generated in two photodiodes having even characteristics. The BPD 523-1 converts the I component of the reception signal of the X-polarized wave into an electric signal, and outputs the same to the amplifier 524-1. The BPD 523-2 converts the Q component of the reception signal of the X-polarized wave into an electric signal, and outputs the same to the amplifier 524-2. The BPD 523-3 converts the I component of the reception signal of the Y-polarized wave into an electric signal, and outputs the same to the amplifier 524-3. The BPD 523-4 converts the Q component of the reception signal of the Y-polarized wave into an electric signal, and outputs the same to the amplifier 524-4. The amplifier 524-i (i is an integer of 1 or larger and 4 or smaller) amplifies the electric signal output from the BPD 523-i, and outputs the same to the digital signal processing unit 530.

[0044] The digital signal processing unit 530 includes analog-to-digital converters (ADCs) 531-1 to 531-4, a front end correction unit 532, a wavelength dispersion compensation unit 533, an adaptive equalization unit 534, a frequency-and-phase offset compensation unit 535, a demapping unit 536, a decoding unit 537, and a characteristic function derivation unit 538. The ADC 531-i (i is an integer of 1 or larger and 4 or smaller) converts the electric signal output from the amplifier 524-i from an analog signal to a digital signal, and outputs the same to the front end correction unit 532.

[0045] The front end correction unit 532 receives an input of the I component of the reception signal of the X-polarized wave from the ADC 531-1, receives an input of the Q component of the reception signal of the X-polarized wave from the ADC 531-2, receives an input of the I component of the reception signal of the Y-polarized wave from the ADC 531-3, and receives an input of the Q component of the reception signal of the Y-polarized wave from the ADC 531-4. The front end correction unit 532 uses each input signal to generate a reception signal subjected to compensation for frequency characteristic in the optical front end 520, and outputs the same to the wavelength dispersion compensation unit 533.

[0046] The wavelength dispersion compensation unit 533 estimates wavelength dispersion performed in the optical fiber transmission path 30, compensates for the estimated wavelength dispersion for the electric signal output from the front end correction unit 532, and outputs the same to the adaptive equalization unit 534. The adaptive equalization unit 534 adaptively performs equalization processing on the reception signal output from the wavelength dispersion compensation unit 533. The adaptive equalization unit 534 outputs a filter coefficient obtained at the time of the equalization processing and a frequency offset to the characteristic function derivation unit 538, and outputs the reception signal subjected to the equalization processing to the frequency-and-phase offset compensation unit 535. The frequency-and-phase offset compensation unit 535 performs processing such as compensation for the frequency offset and phase noise on the reception signal subjected to the equalization processing by the adaptive equalization unit 534.

[0047] The demapping unit 536 determines the symbol of the reception signal output by the frequency-and-phase offset compensation unit 535, and converts the determined symbol into binary data. The decoding unit 537 performs error correction decoding processing such as FEC on the binary data demapped by the demapping unit 536, thereby obtaining a reception bit string.

[0048] The characteristic function derivation unit 538 derives an inverse characteristic of the transmitter 10 and an inverse characteristic of the receiver 50 on the basis of the filter coefficient obtained from the adaptive equalization unit 534 and the frequency offset. In a (pre-) equalization filter circuit of the transmitter or the receiver, waveform distortion occurring between IQ lanes of the transmitter 10 and the receiver 50 can be compensated for by inputting the inverse characteristic derived by the characteristic function derivation unit 538 as a fixed value.

[0049] Note that, although an example of one optical fiber transmission path is described in the above-described embodiment, the same applies in a spatially multiplexed transmission system (for example, a multi-core fiber, a multi-mode fiber, and free space transmission).

[0050] FIG. 2 is a diagram for explaining an outline of processing for deriving the inverse characteristics of the transmitter 10 and the receiver 50 in the first embodiment. In the first embodiment, the fact is used that filter coefficients (h.sub.1, . . . , h.sub.16) obtained by the adaptive equalization unit 534 using an IQ signal and a phase conjugate signal thereof as input signals, and frequency offsets (exp(j.sub.x(n/T)) and (jy(n/T))) include complete information of an IQ characteristic (corresponding to a frequency-dependent 44 matrix) of a reception system including crosstalk between lanes. A symbol interval is represented by n, and a period of the symbol is represented by T. Since this MIMO configuration is not a multistage configuration, a convergence problem that might be a conventional problem is also improved, and once an IQ characteristic function is obtained, there is no need to change the filter coefficient in the adaptive equalization unit 534 thereafter, so that calculation efficiency is improved as compared with a case where compensation is always performed dynamically.

[0051] First, signals (XI, XQ, YI, YQ) on which operations (h.sub.RXI, h.sub.RXQ, h.sub.RY1, h.sub.RYQ) for compensating for a frequency characteristic on the reception side and an operation (h.sub.CD.sup.1) for compensating for wavelength dispersion of a transmission path are performed, and phase conjugate signals of the signals (XI, XQ, YI, YQ) are input to the adaptive equalization unit 534. The adaptive equalization unit 534 performs the adaptive equalization processing. The characteristic function derivation unit 538 receives an input of the filter coefficients (h.sub.1, . . . , h.sub.16) obtained in the process of the adaptive equalization processing by the adaptive equalization unit 534 and the frequency offsets (exp(j.sub.x(n/T)) and (j.sub.y(n/T))). An inverse characteristic H.sub.T.sup.1 () of the transmitter 10 and an inverse characteristic H.sub.R.sup.1() of the receiver 50 are calculated by calculation processing performed by the characteristic function derivation unit 538. Note that, the inverse characteristic H.sub.T.sup.1 () is an aspect of a first inverse characteristic, and the inverse characteristic H.sub.R.sup.1() is an aspect of a second inverse characteristic.

[0052] FIG. 3 is a diagram illustrating a configuration example of a demodulation digital signal processing unit including the adaptive equalization unit 534 in the first embodiment. The demodulation digital signal processing unit includes the front end correction unit 532, the wavelength dispersion compensation unit 533, the adaptive equalization unit 534, and the frequency-and-phase offset compensation unit 535.

[0053] The demodulation digital signal processing unit receives an input of a real component XI and an imaginary component XQ of a reception complex signal of the X-polarized wave and a real component YI and an imaginary component YQ of a reception complex signal of the Y-polarized wave converted into the digital signal by the ADCs 531-1 to 531-4. The demodulation digital signal processing unit convolutes impulse responses (h.sub.RXI, h.sub.RXQ, h.sub.RYI, h.sub.RYQ) for compensating for the frequency characteristic of the receiver 50 and a complex impulse response h.sub.CD.sup.1 for wavelength dispersion compensation into each of the real component XI, imaginary component XQ, real component YI, and imaginary component YQ. Therefore, two complex signals are output for each of an X-polarization component and a Y-polarization component.

[0054] Subsequently, the demodulation digital signal processing unit generates a phase conjugate of each of the two complex signals, and receives an input of eight signals of the real component XI, imaginary component XQ, real component YI, and imaginary component YQ, and the phase conjugates of them for each of the X-polarization component and the Y-polarization component. Therefore, in the adaptive equalization unit 534 of the receiver 50, IQ imbalance, skew between IQ lanes, bias deviations of the IQ modulators 141-1 and 141-2 and the like caused in the transmitter 10 can be dynamically compensated for in addition to impairments caused in the optical fiber transmission path 30 and the receiver 50, and the reception signal is improved in quality.

[0055] Specifically, the demodulation digital signal processing unit applies the impulse response h.sub.RXI for compensating for the frequency characteristic of the receiver 50 and the impulse response h.sub.CD.sup.1 for wavelength dispersion compensation to the real component XI of the reception complex signal of the X-polarization component, and applies the impulse response h.sub.RXQ for compensating for the frequency characteristic of the receiver 50 and the impulse response h.sub.CD.sup.1 for wavelength dispersion compensation to the imaginary component XQ of the reception complex signal of the X-polarization component.

[0056] Similarly, the demodulation digital signal processing unit applies the impulse response h.sub.RYI for compensating for the frequency characteristic of the receiver 50 and the impulse response h.sub.CD.sup.1 for wavelength dispersion compensation to the real component YI of the reception complex signal of the Y-polarization component, and applies the impulse response h.sub.RYQ for compensating for the frequency characteristic of the receiver 50 and the impulse response h.sub.CD.sup.1 for wavelength dispersion compensation to the imaginary component YQ of the reception complex signal of the Y-polarization component.

[0057] The demodulation digital signal processing unit branches each of the real component XI, the imaginary component XQ, the real component YI, and the imaginary component YQ into which the impulse response for compensating for the frequency characteristic of the receiver 50 and the impulse response for wavelength dispersion compensation are convoluted into four signals, directly inputs two signals out of the branched four signals to the adaptive equalization unit 534, and converts the remaining two signals into the phase conjugate signals to input to the adaptive equalization unit 534.

[0058] The adaptive equalization unit 534 adds the real component XI into which the impulse response h.sub.1 is convoluted, the imaginary component XQ into which the impulse response h.sub.5 is convoluted, the real component YI into which the impulse response h.sub.9 is convoluted, and the imaginary component YQ into which the impulse response h.sub.13 is convoluted. The addition signal is multiplied by the frequency offset exp(j.sub.x(n/T)). Moreover, the adaptive equalization unit 534 adds a real component phase conjugate XI* into which the impulse response h.sub.2 is convoluted, an imaginary component phase conjugate XQ into which the impulse response h.sub.6 is convoluted, a real component phase conjugate YI into which the impulse response h.sub.10 is convoluted, and an imaginary component phase conjugate YQ* into which the impulse response h.sub.14 is convoluted. The addition signal is multiplied by a frequency offset exp(j.sub.x(n/T)).

[0059] The demodulation digital signal processing unit adds the addition signal by which the frequency offset exp(j.sub.x(n/T)) is multiplied and the addition signal by which the frequency offset exp(j.sub.x(n/T)) is multiplied, thereby obtaining the reception signal of the X-polarization component. The demodulation digital signal processing unit adds (or subtracts) a transmission data bias correction signal C.sub.X for canceling the bias deviation of the X-polarization component to (or from) the obtained reception signal of the X-polarization component, thereby obtaining a reception signal X.sub.Rsig(n) of the X-polarization component subjected to distortion correction. The demapping unit 536 outputs a reception signal X{circumflex over ()}R.sub.sig(n) obtained as a result of performing symbol determination on the reception signal X.sub.Rsig(n).

[0060] In contrast, the adaptive equalization unit 534 adds the real component XI into which the impulse response h.sub.3 is convoluted, the imaginary component XQ into which the impulse response h.sub.7 is convoluted, the real component YI into which the impulse response hu is convoluted, and the imaginary component YQ into which the impulse response his is convoluted. The addition signal is multiplied by a frequency offset exp(j.sub.y(n/T)). Moreover, the adaptive equalization unit 534 adds a real component phase conjugate XI* into which the impulse response h.sub.4 is convoluted, an imaginary component phase conjugate XQ into which the impulse response h.sub.12 is convoluted, a real component phase conjugate YI into which the impulse response h.sub.16 is convoluted, and an imaginary component phase conjugate YQ* into which the impulse response h.sub.14 is convoluted. The addition signal is multiplied by a frequency offset exp(j.sub.y(n/T)).

[0061] The demodulation digital signal processing unit adds the addition signal to which the frequency offset exp(j.sub.y(n/T)) is applied and the addition signal to which the frequency offset exp(j.sub.y(n/T) is applied, thereby obtaining the reception signal of the Y-polarization component. The demodulation digital signal processing unit adds (or subtracts) a transmission data bias correction signal Cy for canceling the bias deviation of the Y-polarization component to (or from) the obtained reception signal of the Y-polarization component, thereby obtaining a reception signal Y.sub.Rsig(n) of the X-polarization component subjected to distortion correction. The demapping unit 536 outputs a reception signal Y{circumflex over ()}R.sub.sig(n) obtained as a result of performing symbol determination on the reception signal Y.sub.Rsig(n).

[0062] Note that, the complex impulse response h.sub.CD.sup.1 for wavelength dispersion compensation, the impulse responses h.sub.1 to h.sub.16, and the frequency offsets exp(j.sub.x(n/T)), exp(j.sub.x(n/T)), exp(j.sub.y(n/T)), and exp(j.sub.y(n/T)) are adaptively and dynamically changed. The receiver 50 acquires these values by any method.

[0063] Note that, the convolution of the impulse responses h.sub.RXI, h.sub.RXQ, h.sub.RYI, and h.sub.RYQ corresponds to processing of the front end correction unit 532 illustrated in FIG. 1, and the convolution of the impulse response h.sub.CD.sup.1 for wavelength dispersion compensation corresponds to processing of the wavelength dispersion compensation unit 533. Multiplication processing of the frequency offsets exp(j.sub.x(n/T)), exp(j.sub.x(n/T)), exp(j.sub.y(n/T)), and exp(j.sub.y(n/T) on the addition signals corresponds to the function of the frequency-and-phase offset compensation unit 535.

[0064] The filter coefficients (h.sub.1, . . . , h.sub.16) obtained in the processing of the adaptive equalization unit 534 in the demodulation digital signal processing unit and the frequency offsets (exp(j.sub.x(n/T)) and (j.sub.y(n/T))) are output to the characteristic function derivation unit 538.

[0065] Next, a detailed principle of the above-described demodulation digital signal processing unit will be described.

[0066] First, a general variable a(x) is set to ()a()=a*() while paying attention to following Expression (1) when FT(S.sub.i,in(t))=S.sub.i,in(). Note that, () is attached above a.

[00001]$\begin{array}{cc}\left[\mathrm{Math}.1\right]& \end{array}$$\begin{array}{cc}\mathrm{FT}\left(s_{i,\mathrm{in}}^{*}\right)=s_{i,\mathrm{in}}^{*}\left(-\right)& \mathrm{Expression}\left(1\right)\end{array}$

[0067] In order to consider skew and crosstalk between the IQ lanes, a signal waveform is represented by 41 vector representation as following Expression (2).

[00002]$\begin{array}{cc}\left[\mathrm{Math}.2\right]& \end{array}$$\begin{array}{cc}{s}_{\mathrm{in}}\left(\right)=\left(\begin{array}{c}{s}_{x,\mathrm{in}}\left(\right)\\ {\tilde{s}}_{x,\mathrm{in}}\left(\right)\\ s_{y,\mathrm{in}}\left(\right)\\ {\tilde{s}}_{x,\mathrm{in}}\left(\right)\end{array}\right),s_{\mathrm{out}}\left(\right)=\left(\begin{array}{c}{s}_{x,\mathrm{out}}\left(\right)\\ {\tilde{s}}_{x,\mathrm{out}}\left(\right)\\ s_{y,\mathrm{out}}\left(\right)\\ {\tilde{s}}_{x,\mathrm{out}}\left(\right)\end{array}\right)& \mathrm{Expression}\left(2\right)\end{array}$

[0068] Herein, each vector before and after propagation can be expressed as following Expression (3).

[00003]$\begin{array}{cc}\begin{array}{cc}\left[\mathrm{Math}.3\right]& \end{array}& \end{array}$$\begin{array}{cc}{s}_{\mathrm{out}}\left(\right)=H_R\left(\right)\mathrm{FT}\left(H_{\mathrm{fR}}\left(t\right){\mathrm{FT}}^{-1}\left(H_{\mathrm{CD}}\left(\right)H_{\mathrm{couple}}\left(\right)\mathrm{FT}\left(H_{\mathrm{fT}}\left(t\right)\left({\mathrm{FT}}^{-1}\left(H_T\left(\right)s_{\mathrm{in}}\left(\right)\right)\right)\right)\right)\right)& \mathrm{Expression}\left(3\right)\end{array}$

[0069] Each sign in Expression (3) means the following content. [0070] H.sub.R(): 44 matrix representing the IQ characteristics of the transmitter 10 and the receiver 50 including the IQ skew, imbalance, and crosstalk between the lanes on the reception side [0071] H.sub.fR(t): a 44 matrix representing a frequency of local oscillation light of each receiver [0072] H.sub.CD(): a 44 matrix representing an influence of the wavelength dispersion of the transmission path [0073] H.sub.couple(): a 44 matrix representing channel crosstalk in the transmission path [0074] H.sub.fT(t): a 44 matrix representing a frequency of carrier light of each transmitter [0075] H.sub.T(): a 44 matrix representing the IQ characteristics of the transmitter 10 and the receiver 50 including the IQ skew, imbalance, and crosstalk between the lanes on the transmission side

[0076] It is possible to calculate S.sub.in() from S.sub.out() by calculating an inverse matrix of each sign in Expression (3). Following Expression (4) expresses S.sub.in().

[00004]$\begin{array}{cc}\left[\mathrm{Math}.4\right]& \end{array}$$\begin{array}{cc}{s}_{\mathrm{in}}\left(\right)=H_T^{-1}\left(\right)\mathrm{FT}(H_{\mathrm{fT}}^{-1}\left(t\right){\mathrm{FT}}^{-1}\left(H_{\mathrm{couple}}^{-1}\left(\right)H_{\mathrm{CD}}^{-1}\left(\right)\left(H_{\mathrm{CD}}^{-1}\left(\right){\mathrm{FT}}^{}\left(H_{\mathrm{fR}}^{-1}\left(t\right){\mathrm{FT}}^{-1}\left(H_R^{-1}\left(\right)s_{\mathrm{out}}\left(\right)\right)\right)\right)\right)& \mathrm{Expression}\left(4\right)\end{array}$

[0077] Each matrix of H.sub.R.sup.1(), H.sub.fR.sup.1(t), H.sub.CD.sup.1(), H.sub.couple.sup.1(), H.sub.fT.sup.1 (t), and H.sub.T.sup.1 (t) in Expression (4) is expressed as following Expressions (5) to (10), respectively.

[00005]$\begin{array}{cc}\left[\mathrm{Math}.5\right]& \end{array}$$\begin{array}{cc}{H}_R^{-1}\left(\right):=\left(\begin{array}{cccc}{h}_{R11}\left(\right)& h_{R12}\left(\right)& h_{R13}\left(\right)& h_{R14}\left(\right)\\ h_{R21}\left(\right)& h_{R22}\left(\right)& h_{R23}\left(\right)& h_{R24}\left(\right)\\ h_{R31}\left(\right)& h_{R32}\left(\right)& h_{R33}\left(\right)& h_{R34}\left(\right)\\ h_{R41}\left(\right)& h_{R42}\left(\right)& h_{R43}\left(\right)& h_{R44}\left(\right)\end{array}\right)& \mathrm{Expression}\left(5\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.6\right]& \end{array}$$\begin{array}{cc}{H}_{\mathrm{fR}}^{-1}\left(t\right):=\left(\begin{array}{cccc}\exp \left(i{}_{R}t\right)& 0& 0& 0\\ 0& \exp \left(-i{}_{R}t\right)& 0& 0\\ 0& 0& \exp \left(i{}_{R}t\right)& 0\\ 0& 0& 0& \exp \left(-i{}_{R}t\right)\end{array}\right)& \mathrm{Expression}\left(6\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.7\right]& \end{array}$$\begin{array}{cc}{H}_{\mathrm{CD}}^{-1}\left(\right):=\left(\begin{array}{cccc}{h}_{\mathrm{CD}}^{-1}\left(\right)& 0& 0& 0\\ 0& {\tilde{h}}_{\mathrm{CD}}^{-1}\left(\right)& 0& 0\\ 0& 0& h_{\mathrm{CD}}^{-1}\left(\right)& 0\\ 0& 0& 0& {\tilde{h}}_{\mathrm{CD}}^{-1}\left(\right)\end{array}\right)& \mathrm{Expression}\left(7\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.8\right]& \end{array}$$\begin{array}{cc}{H}_{\mathrm{couple}}^{-1}\left(\right):=\left(\begin{array}{cccc}{h}_{\mathrm{xx}}^{-1}\left(\right)& 0& h_{\mathrm{xy}}^{-1}\left(\right)& 0\\ 0& {\tilde{h}}_{\mathrm{xx}}^{-1}\left(\right)& 0& {\tilde{h}}_{\mathrm{xy}}^{-1}\left(\right)\\ h_{\mathrm{yx}}^{-1}\left(\right)& 0& h_{\mathrm{yy}}^{-1}\left(\right)& 0\\ 0& {\tilde{h}}_{\mathrm{yx}}^{-1}\left(\right)& 0& {\tilde{h}}_{\mathrm{yy}}^{-1}\left(\right)\end{array}\right)& \mathrm{Expression}\left(8\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.9\right]& \end{array}$$\begin{array}{cc}{H}_{\mathrm{fT}}^{-1}\left(t\right):=\left(\begin{array}{cccc}\exp \left(-i{}_{T_x}t\right)& 0& 0& 0\\ 0& \exp \left(i{}_{T_x}t\right)& 0& 0\\ 0& 0& \exp \left(-i{}_{T_y}t\right)& 0\\ 0& 0& 0& \exp \left(i{}_{T_y}t\right)\end{array}\right)& \mathrm{Expression}\left(9\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.10\right]& \end{array}$$\begin{array}{cc}{H}_T^{-1}\left(\right)=\left(\begin{array}{cccc}{h}_{T11}\left(\right)& h_{T12}\left(\right)& h_{T13}\left(\right)& h_{T14}\left(\right)\\ h_{T21}\left(\right)& h_{T22}\left(\right)& h_{T23}\left(\right)& h_{T24}\left(\right)\\ h_{T31}\left(\right)& h_{T32}\left(\right)& h_{T33}\left(\right)& h_{T34}\left(\right)\\ h_{T41}\left(\right)& h_{T42}\left(\right)& h_{T43}\left(\right)& h_{T44}\left(\right)\end{array}\right)& \mathrm{Expression}\left(10\right)\end{array}$

[0078] Herein, it is assumed that the local light on the reception side has the same frequency and phase fluctuation between the polarized waves. Supposing that .sub.x=.sub.Rx.sub.T, and .sub.y=.sub.RyT, A(), B(), and M are defined as Expressions (11) to (13), respectively. Note that, x of .sub.Rx is a subscript of R, and y of .sub.Ry is a subscript of R.

[00006]$\begin{array}{cc}\left[\mathrm{Math}.11\right]& \end{array}$$\begin{array}{cc}A\left(\right)=\left(\begin{array}{cccc}{h}_{T11}\left(-{}_x\right)& h_{T12}\left(+{}_x\right)& h_{T13}\left(-{}_y\right)& h_{T14}\left(+{}_y\right)\\ h_{T21}\left(-{}_x\right)& h_{T22}\left(+{}_x\right)& h_{T23}\left(-{}_y\right)& h_{T24}\left(+{}_y\right)\\ h_{T31}\left(-{}_x\right)& h_{T32}\left(+{}_x\right)& h_{T33}\left(-{}_y\right)& h_{T34}\left(+{}_y\right)\\ h_{T41}\left(-{}_x\right)& h_{T42}\left(+{}_x\right)& h_{T43}\left(-{}_y\right)& h_{T44}\left(+{}_y\right)\end{array}\right)\left(\begin{array}{cccc}{h}_{\mathrm{xx}}^{-1}\left(-{}_{R_x}\right)& 0& h_{\mathrm{xy}}^{-1}\left(-{}_{R_x}\right)& 0\\ 0& 0& 0& 0\\ h_{\mathrm{yx}}^{-1}\left(-{}_{R_y}\right)& 0& h_{\mathrm{yy}}^{-1}\left(-{}_{R_y}\right)& 0\\ 0& 0& 0& 0\end{array}\right)\left(\begin{array}{cccc}{h}_{R11}\left(\right)& h_{R12}\left(\right)& h_{R13}\left(\right)& h_{R14}\left(\right)\\ h_{R21}\left(\right)& h_{R22}\left(\right)& h_{R23}\left(\right)& h_{R24}\left(\right)\\ h_{R31}\left(\right)& h_{R32}\left(\right)& h_{R33}\left(\right)& h_{R34}\left(\right)\\ h_{R41}\left(\right)& h_{R42}\left(\right)& h_{R43}\left(\right)& h_{R44}\left(\right)\end{array}\right)& \mathrm{Expression}\left(11\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.12\right]& \end{array}$$B\left(\right)=\begin{array}{cc}\left(\begin{array}{cccc}{h}_{T11}\left(-{}_x\right)& h_{T12}\left(+{}_x\right)& h_{T13}\left(-{}_y\right)& h_{T14}\left(+{}_y\right)\\ h_{T21}\left(-{}_x\right)& h_{T22}\left(+{}_x\right)& h_{T23}\left(-{}_y\right)& h_{T24}\left(+{}_y\right)\\ h_{T31}\left(-{}_x\right)& h_{T32}\left(+{}_x\right)& h_{T33}\left(-{}_y\right)& h_{T34}\left(+{}_y\right)\\ h_{T41}\left(-{}_x\right)& h_{T42}\left(+{}_x\right)& h_{T43}\left(-{}_y\right)& h_{T44}\left(+{}_y\right)\end{array}\right)\left(\begin{array}{cccc}0& 0& 0& 0\\ 0& {\tilde{h}}_{\mathrm{xx}}^{-1}\left(-{}_{R_x}\right)& 0& {\tilde{h}}_{\mathrm{xy}}^{-1}\left(-{}_{R_x}\right)\\ 0& 0& 0& 0\\ 0& {\tilde{h}}_{\mathrm{yx}}^{-1}\left(-{}_{R_y}\right)& 0& {\tilde{h}}_{\mathrm{yy}}^{-1}\left(-{}_{R_y}\right)\end{array}\right)\left(\begin{array}{cccc}{h}_{R11}\left(\right)& h_{R12}\left(\right)& h_{R13}\left(\right)& h_{R14}\left(\right)\\ h_{R21}\left(\right)& h_{R22}\left(\right)& h_{R23}\left(\right)& h_{R24}\left(\right)\\ h_{R31}\left(\right)& h_{R32}\left(\right)& h_{R33}\left(\right)& h_{R34}\left(\right)\\ h_{R41}\left(\right)& h_{R42}\left(\right)& h_{R43}\left(\right)& h_{R44}\left(\right)\end{array}\right)& \mathrm{Expression}\left(12\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.13\right]& \end{array}$$\begin{array}{cc}M=\left(\begin{array}{cccc}1& i& 0& 0\\ 1& -i& 0& 0\\ 0& 0& 1& i\\ 0& 0& 1& -i\end{array}\right)& \mathrm{Expression}\left(13\right)\end{array}$

[0079] In this case, S.sub.in() is expressed as following Expression (14).

[00007]$\begin{array}{cc}\left[\mathrm{Math}.14\right]& \end{array}$$\begin{array}{cc}{s}_{\mathrm{in}}\left(\right)=\mathrm{FT}\left(\begin{array}{cccc}\exp \left(i{}_{x}t\right)& 0& 0& 0\\ 0& \exp \left(i{}_{x}t\right)& 0& 0\\ 0& 0& \exp \left(i{}_{y}t\right)& 0\\ 0& 0& 0& \exp \left(i{}_{y}t\right)\end{array}\right){\mathrm{FT}}^{-1}A\left(\right)M\left(\begin{array}{c}{h}_{\mathrm{CD}}^{-1}\left(\right)\mathrm{Re}\left(s_{x,\mathrm{out}}\left(\right)\right)\\ h_{\mathrm{CD}}^{-1}\left(\right)\mathrm{Im}\left(s_{x,\mathrm{out}}\left(\right)\right)\\ h_{\mathrm{CD}}^{-1}\left(\right)\mathrm{Re}\left(s_{y,\mathrm{out}}\left(\right)\right)\\ h_{\mathrm{CD}}^{-1}\left(\right)\mathrm{Im}\left(s_{x,\mathrm{out}}\left(\right)\right)\end{array}\right)+\mathrm{FT}\left(\begin{array}{cccc}\exp \left(-i{}_{x}t\right)& 0& 0& 0\\ 0& \exp \left(-i{}_{x}t\right)& 0& 0\\ 0& 0& \exp \left(-i{}_{y}t\right)& 0\\ 0& 0& 0& \exp \left(-i{}_{y}t\right)\end{array}\right){\mathrm{FT}}^{-1}B\left(\right)M\left(\begin{array}{c}{\left[h_{\mathrm{CD}}^{-1}\left(\right)\mathrm{Re}\left(s_{x,\mathrm{out}}\left(\right)\right)\right]}^{*}\\ {\left[h_{\mathrm{CD}}^{-1}\left(\right)\mathrm{Im}\left(s_{x,\mathrm{out}}\left(\right)\right)\right]}^{*}\\ {\left[h_{\mathrm{CD}}^{-1}\left(\right)\mathrm{Re}\left(s_{y,\mathrm{out}}\left(\right)\right)\right]}^{*}\\ {\left[h_{\mathrm{CD}}^{-1}\left(\right)\mathrm{Im}\left(s_{x,\mathrm{out}}\left(\right)\right)\right]}^{*}\end{array}\right)& \mathrm{Expression}\left(14\right)\end{array}$

[0080] When Expression (14) is converted into a time domain, this is expressed as Expression (15).

[00008]$\begin{array}{cc}\left[\mathrm{Math}.15\right]& \end{array}$$\begin{array}{cc}{s}_{\mathrm{in}}\left(t\right)=\left(\begin{array}{cccc}\exp \left(i{}_{x}t\right)& 0& 0& 0\\ 0& \exp \left(i{}_{x}t\right)& 0& 0\\ 0& 0& \exp \left(i{}_{y}t\right)& 0\\ 0& 0& 0& \exp \left(i{}_{y}t\right)\end{array}\right)\left[A\left(\right)M*\left(\begin{array}{c}{h}_{\mathrm{CD}}^{-1}\left(t\right)*\mathrm{Re}\left(s_{x,\mathrm{out}}\left(t\right)\right)\\ h_{\mathrm{CD}}^{-1}\left(t\right)*\mathrm{Im}\left(s_{x,\mathrm{out}}\left(t\right)\right)\\ h_{\mathrm{CD}}^{-1}\left(t\right)*\mathrm{Re}\left(s_{y,\mathrm{out}}\left(t\right)\right)\\ h_{\mathrm{CD}}^{-1}\left(t\right)*\mathrm{Im}\left(s_{x,\mathrm{out}}\left(t\right)\right)\end{array}\right)\right]+\left(\begin{array}{cccc}\exp \left(-i{}_{x}t\right)& 0& 0& 0\\ 0& \exp \left(-i{}_{x}t\right)& 0& 0\\ 0& 0& \exp \left(-i{}_{y}t\right)& 0\\ 0& 0& 0& \exp \left(-i{}_{y}t\right)\end{array}\right)\left[B\left(t\right)M*\left(\begin{array}{c}{\left[h_{\mathrm{CD}}^{-1}\left(t\right)\mathrm{Re}\left(s_{x,\mathrm{out}}\left(t\right)\right)\right]}^{*}\\ {\left[h_{\mathrm{CD}}^{-1}\left(t\right)\mathrm{Im}\left(s_{x,\mathrm{out}}\left(t\right)\right)\right]}^{*}\\ {\left[h_{\mathrm{CD}}^{-1}\left(t\right)\mathrm{Re}\left(s_{y,\mathrm{out}}\left(t\right)\right)\right]}^{*}\\ {\left[h_{\mathrm{CD}}^{-1}\left(t\right)\mathrm{Im}\left(s_{x,\mathrm{out}}\left(t\right)\right)\right]}^{*}\end{array}\right)\right]& \mathrm{Expression}\left(15\right)\end{array}$

[0081] According to Expression (15), in order to compensate for distortion occurring during the transmission, it is sufficient to perform wavelength dispersion compensation on each of a real part and an imaginary part of the signal of each polarized wave, then convolute an appropriate function (for example, corresponding to A(t)M and B(t)M based on Expressions (11) to (13) above) in a 44 matrix format into the set of four signals or their complex conjugate signals, and multiply by a term for correcting the frequency offset of the local oscillation light. One of adaptive filters that perform equalization on this principle is an 82 MIMO configuration illustrated in FIG. 3.

[0082] In the 82 MIMO configuration, odd-numbered rows of A(t)M and B(t)M are adaptively obtained. That is, h.sub.1, h.sub.3, . . . , h.sub.15 correspond to elements in the odd-numbered rows of A(t)M, and h.sub.2, h.sub.4, . . . , h.sub.16 correspond to elements in the odd-numbered rows of B(t)M. Under the condition that complex conjugate of a component of a (2i+1)-th row of S.sub.in() is equal to a component of a 2i-th row, components of even-numbered rows of A() and B() are obtained from components of the odd-numbered rows of the other, and A() and B() are obtained as Expressions (16) and (17) from Fourier transform of the filter coefficients h.sub.1 to h.sub.16 of the 82 MIMO configuration.

[00009]$\begin{array}{cc}\left[\mathrm{Math}.16\right]& \end{array}$$\begin{array}{cc}A\left(\right)=M^{-1}(\begin{array}{cccc}{h}_1\left(\right)& h_5\left(\right)& h_9\left(\right)& h_{13}\left(\right)\\ {\tilde{h}}_2\left(\right)& {\stackrel{}{h}}_6\left(\right)& {\stackrel{}{h}}_{10}\left(\right)& {\tilde{h}}_{14}\left(\right)\\ h_3\left(\right)& h_7\left(\right)& h_{11}\left(\right)& h_{15}\left(\right)\\ {\tilde{h}}_4\left(\right)& {\stackrel{}{h}}_8\left(\right)& {\tilde{h}}_{12}\left(\right)& {\tilde{h}}_{16}\left(\right)\end{array})& \mathrm{Expression}\left(16\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.17\right]& \end{array}$$\begin{array}{cc}B\left(\right)=M^{-1}\left(\begin{array}{cccc}{h}_2\left(\right)& h_6\left(\right)& h_{10}\left(\right)& h_{14}\left(\right)\\ {\tilde{h}}_1\left(\right)& {\stackrel{}{h}}_5\left(\right)& {\stackrel{}{h}}_9\left(\right)& {\tilde{h}}_{13}\left(\right)\\ h_4\left(\right)& h_8\left(\right)& h_{12}\left(\right)& h_{16}\left(\right)\\ {\tilde{h}}_3\left(\right)& {\stackrel{}{h}}_7\left(\right)& {\tilde{h}}_{11}\left(\right)& {\tilde{h}}_{15}\left(\right)\end{array}\right)& \mathrm{Expression}\left(17\right)\end{array}$

[0083] In the first embodiment, the characteristic function derivation unit 538 illustrated in FIG. 4 derives the characteristic H.sub.T() of the transmitter 10 and the characteristic H.sub.R() of the receiver 50 on the basis of A(). As for H.sub.T() and H.sub.R(), a complex conjugate of a component of a (2i+1)-th row of an input/output vector is equal to a component of a 2i-th row with respect to a value that can be taken by each element, and even if a matrix representing any polarization rotation and phase rotation is multiplied on the transmission path side, this can be regarded as the same transceiver characteristic, so that H.sub.T() and H.sub.R() can reduce the degrees of freedom from 16 complex numbers to 4 complex numbers, and can be regarded as Expressions (18) and (19).

[00010]$\begin{array}{cc}\left[\mathrm{Math}.18\right]& \end{array}$$\begin{array}{cc}{H}_T^{-1}\left(\right)=\left(\begin{array}{cccc}1& {\stackrel{}{h}}_{T21}\left(\right)& 0& {\tilde{h}}_{T23}\left(\right)\\ h_{T21}\left(\right)& 1& h_{T23}\left(\right)& 0\\ 0& {\tilde{h}}_{T41}\left(\right)& 1& {\tilde{h}}_{T43}\left(\right)\\ h_{T41}\left(\right)& 0& h_{T43}\left(\right)& 1\end{array}\right)& \mathrm{Expression}\left(18\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.19\right]& \end{array}$$\begin{array}{cc}{H}_R^{-1}\left(\right)=\left(\begin{array}{cccc}1& h_{R12}\left(\right)& 0& h_{R14}\left(\right)\\ {\stackrel{}{h}}_{R12}\left(\right)& 1& {\stackrel{}{h}}_{R14}\left(\right)& 0\\ 0& h_{R32}\left(\right)& 1& h_{R34}\left(\right)\\ {\stackrel{}{h}}_{R32}\left(\right)& 0& {\stackrel{}{h}}_{R34}\left(\right)& 1\end{array}\right)& \mathrm{Expression}\left(19\right)\end{array}$

[0084] Herein, from the definition of A, A() is expressed as following Expression (20).

[00011]$\begin{array}{cc}\left[\mathrm{Math}.20\right]& \end{array}$$\begin{array}{cc}A\left(\right)=\left(\begin{array}{cccc}1& {\stackrel{}{h}}_{T21}\left(-{}_x\right)& 0& {\tilde{h}}_{T23}\left(-{}_y\right)\\ h_{T21}\left(-{}_x\right)& 1& h_{T23}\left(-{}_y\right)& 0\\ 0& {\tilde{h}}_{T41}\left(-{}_x\right)& 1& {\tilde{h}}_{T43}\left(-{}_y\right)\\ h_{T41}\left(-{}_x\right)& 0& h_{T43}\left(-{}_y\right)& 1\end{array}\right)\left(\begin{array}{cccc}{h}_{\mathrm{xx}}^{-1}\left(-{}_{R_x}\right)& 0& h_{\mathrm{xy}}^{-1}\left(-{}_{R_x}\right)& 0\\ 0& 0& 0& 0\\ h_{\mathrm{yx}}^{-1}\left(-{}_{R_y}\right)& 0& h_{\mathrm{yy}}^{-1}\left(-{}_{R_y}\right)& 0\\ 0& 0& 0& 0\end{array}\right)\left(\begin{array}{cccc}1& h_{R12}\left(\right)& 0& h_{R14}\left(\right)\\ {\stackrel{}{h}}_{R12}\left(\right)& 1& {\stackrel{}{h}}_{R14}\left(\right)& 0\\ 0& h_{R32}\left(\right)& 1& h_{R34}\left(\right)\\ {\stackrel{}{h}}_{R32}\left(\right)& 0& {\stackrel{}{h}}_{R34}\left(\right)& 1\end{array}\right)& \mathrm{Expression}\left(20\right)\end{array}$

[0085] As a result, A() is expressed as following Expression (21).

[00012]$\begin{array}{cc}\left[\mathrm{Math}.21\right]& \end{array}$$\begin{array}{cc}A\left(\right)=\left(\begin{array}{cccc}{h}_{\mathrm{xx}}^{-1}& h_{\mathrm{xx}}^{-1}{h}_{R12}+h_{\mathrm{xy}}^{-1}{h}_{R32}& h_{\mathrm{xy}}^{-1}& h_{\mathrm{xx}}^{-1}{h}_{R14}+h_{\mathrm{xy}}^{-1}{h}_{R34}\\ h_{\mathrm{xx}}^{-1}{h}_{T21}+h_{\mathrm{yx}}^{-1}{h}_{T23}& -& h_{\mathrm{xy}}^{-1}{h}_{T21}+h_{\mathrm{yy}}^{-1}{h}_{T23}& -\\ h_{\mathrm{yx}}^{-1}& h_{\mathrm{yx}}^{-1}{h}_{R12}+h_{\mathrm{yy}}^{-1}{h}_{R32}& h_{\mathrm{yy}}^{-1}& h_{\mathrm{yx}}^{-1}{h}_{R12}+h_{\mathrm{yy}}^{-1}{h}_{R32}\\ h_{\mathrm{xx}}^{-1}{h}_{T41}+h_{\mathrm{yx}}^{-1}{h}_{T43}& -& h_{\mathrm{xy}}^{-1}{h}_{T41}+h_{\mathrm{yy}}^{-1}{h}_{T43}& -\end{array}\right)& \mathrm{Expression}\left(21\right)\end{array}$

[0086] When A() is obtained as described above and A() is placed as following Expression (22), a relationship of Expression (23) is derived on the basis of Expressions (21) and (22).

[00013]$\begin{array}{cc}\left[\mathrm{Math}.22\right]& \\ A\left(\right)=\left(\begin{array}{cccc}{A}_{11}& A_{12}& A_{13}& A_{14}\\ A_{21}& A_{22}& A_{23}& A_{24}\\ A_{31}& A_{32}& A_{33}& A_{34}\\ A_{41}& A_{42}& A_{43}& A_{44}\end{array}\right)& \mathrm{Expression}\left(22\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.23\right]& \\ \left(\begin{array}{cc}{h}_{\mathrm{xx}}^{-1}& h_{\mathrm{xy}}^{-1}\\ h_{\mathrm{yx}}^{-1}& h_{\mathrm{yy}}^{-1}\end{array}\right)=\left(\begin{array}{cc}{A}_{11}& A_{13}\\ A_{31}& A_{33}\end{array}\right)& \mathrm{Expression}\left(23\right)\end{array}$

[0087] Moreover, a relationship of Expression (24) is also derived on the basis of Expressions (21) and (22).

[00014]$\begin{array}{cc}\left[\mathrm{Math}.24\right]& \\ \left(\begin{array}{cc}{A}_{21}& A_{23}\\ A_{41}& A_{43}\end{array}\right)=\left(\begin{array}{cc}{h}_{T21}& h_{T23}\\ h_{T41}& h_{T43}\end{array}\right)\left(\begin{array}{cc}{h}_{\mathrm{xx}}^{-1}& h_{\mathrm{xy}}^{-1}\\ h_{\mathrm{yx}}^{-1}& h_{\mathrm{yy}}^{-1}\end{array}\right)& \mathrm{Expression}\left(24\right)\end{array}$

[0088] Expression (25) is derived on the basis of Expressions (24) and (23).

[00015]$\begin{array}{cc}\left[\mathrm{Math}.25\right]& \\ \left(\begin{array}{cc}{h}_{T21}& h_{T23}\\ h_{T41}& h_{T43}\end{array}\right)=\left(\begin{array}{cc}{A}_{21}& A_{23}\\ A_{41}& A_{43}\end{array}\right){\left(\begin{array}{cc}{A}_{11}& A_{13}\\ A_{31}& A_{33}\end{array}\right)}^{-1}& \mathrm{Expression}\left(25\right)\end{array}$

[0089] Moreover, a relationship of Expression (26) is also derived on the basis of Expressions (21) and (22).

[00016]$\begin{array}{cc}\left[\mathrm{Math}.26\right]& \\ \left(\begin{array}{cc}{A}_{12}& A_{14}\\ A_{32}& A_{34}\end{array}\right)=\left(\begin{array}{cc}{h}_{\mathrm{xx}}^{-1}& h_{\mathrm{xy}}^{-1}\\ h_{\mathrm{yx}}^{-1}& h_{\mathrm{yy}}^{-1}\end{array}\right)\left(\begin{array}{cc}{h}_{R12}& h_{R14}\\ h_{R32}& h_{R34}\end{array}\right)& \mathrm{Expression}\left(26\right)\end{array}$

[0090] Expression (27) is derived on the basis of Expressions (26) and (23).

[00017]$\begin{array}{cc}\left[\mathrm{Math}.27\right]& \\ \left(\begin{array}{cc}{h}_{R12}& h_{R14}\\ h_{R32}& h_{R34}\end{array}\right)={\left(\begin{array}{cc}{A}_{11}& A_{13}\\ A_{31}& A_{33}\end{array}\right)}^{-1}\left(\begin{array}{cc}{A}_{12}& A_{14}\\ A_{32}& A_{34}\end{array}\right)& \mathrm{Expression}\left(27\right)\end{array}$

[0091] Each component is calculated by performing the calculation as described above. Although factors of .sub.x and .sub.y remain in each element of obtained H.sub.T.sup.1, since the frequency offset of the local oscillation light is obtained at the time of adaptive equalization, it is possible to perform correction by performing convolution operation of this value. The characteristic function derivation unit 538 derives the inverse characteristic H.sub.T.sup.1 () and the inverse characteristic H.sub.R.sup.1() using the expressions obtained as described above. Specifically, the characteristic function derivation unit 538 derives the inverse characteristic H.sub.T.sup.1 () by applying each matrix element of Expression (25) to Expression (18), and derives the inverse characteristic H.sub.R.sup.1() by applying each matrix element of Expression (27) to Expression (19).

[0092] FIG. 5 is a flowchart illustrating a flow of processing of the receiver 50 in the first embodiment.

[0093] The optical front end 520 receives the optical signal (polarization multiplexed signal) transmitted through the optical fiber transmission path 30 (step S101). Each functional unit in the optical front end 520 performs, on the received optical signal, polarization separation by the polarization separation unit 521, extraction of the I component and the Q component of the X-polarized wave by the optical 90-degree hybrid coupler 522-1, extraction of the I component and the Q component of the Y-polarized wave by the optical 90-degree hybrid coupler 522-2, conversion into the electric signal, and amplification of the electric signal.

[0094] The ADC 531-i converts the electric signal output from the amplifier 524-i from the analog signal to the digital signal (step S102). The front end correction unit 532 uses each input signal to generate the reception signal subjected to the compensation for the frequency characteristic in the optical front end 520 (step S103). The wavelength dispersion compensation unit 533 performs the wavelength dispersion compensation on the electric signal output from the front end correction unit 532 (step S104).

[0095] The adaptive equalization unit 534 performs the equalization processing on the reception signal output from the wavelength dispersion compensation unit 533 (step S105). The adaptive equalization unit 534 outputs the filter coefficient obtained at the time of the equalization processing and the frequency offset to the characteristic function derivation unit 538. Note that, in FIG. 5, the description of the frequency-and-phase offset compensation unit 535 and after that will be omitted.

[0096] The characteristic function derivation unit 538 derives the inverse characteristics of the transmitter 10 and the receiver 50 on the basis of the filter coefficient output from the adaptive equalization unit 534 and the frequency offset (step S106).

[0097] According to the receiver 50 configured as described above, the inverse characteristic H.sub.T.sup.1() representing the inverse characteristic of the transmitter 10 and the inverse characteristic H.sub.R.sup.1() representing the inverse characteristic of the receiver 50 are calculated on the basis of the adaptive equalization unit 534 that performs the equalization processing on the input signal using the polarization multiplexed reception signal and the phase conjugate signal of the polarization multiplexed reception signal as the input signals, the filter coefficient obtained at the time of the equalization processing performed by the adaptive equalization unit 534, and the frequency offset. As described above, by analyzing the polarization multiplexed reception signal in the receiver 50, it is possible to estimate the characteristic of the transmitter 10 and the characteristic of the receiver 50 including the IQ crosstalk between the lanes beyond polarization, core, and mode. Therefore, it is possible to implement a highly efficient and highly reliable optical communication system by inputting an inverse characteristic function to the equivalent filters of the transmitter 10 and the receiver 50 to correct communication distortion.

Modification Example of First Embodiment

[0098] In the above-described embodiment, the configuration in which the signal input to the adaptive equalization unit 534 is the set of the IQ signal and the phase conjugate signal of the IQ signal has been described. As the signal input to the adaptive equalization unit 534, a set of signals mathematically equivalent to the set of the IQ signal and the phase conjugate signal of the IQ signal may be used. For example, a set of a total of four signals including signals obtained by performing wavelength dispersion compensation on a complex signal and a phase conjugate signal thereof, and phase conjugate signals of the two signals is associated with the set of the IQ signal and the phase conjugate signal of the IQ signal by orthogonal linear conversion. Therefore, the set of four signals described above may be input to the adaptive equalization unit 534 as the input signal. Even in such a configuration, it is possible to measure the characteristic of the transmitter 10 and the characteristic of the receiver 50. At the time of measurement, the filter coefficient may be converted using inverse conversion of the set used as the input signal and the set of the IQ signal and the phase conjugate signal thereof.

Second Embodiment

[0099] Although only polarization division multiplexing is considered in the first embodiment, in a second embodiment, a configuration of extending to any multiplex number N (N2) obtained by combining space division multiplexing and wavelength division multiplexing in addition to the polarization division multiplexing will be described. A basic system configuration of a digital coherent optical transmission system of the second embodiment is different from the digital coherent optical transmission system 1 illustrated in FIG. 1 in a following configuration.

[0100] A transmitter 10 further includes as many transmission units 100 as the number of channels of wavelength division multiplexing (WDM). The transmission units 100 output optical signals of different wavelengths. A WDM multiplexer, an optical fiber transmission path 30, and a WDM demultiplexer are provided between the transmitter 10 and a receiver 50. The WDM multiplexer multiplexes optical signals output by the transmission units 100, and outputs the same to the optical fiber transmission path 30. The WDM demultiplexer demultiplexes the optical signal transmitted through the optical fiber transmission path 30 by wavelength. The receiver 50 further includes as many reception units 500 as the number of channels of the WDM. Each reception unit 500 receives the optical signal demultiplexed by a WDM demultiplexer 40. Wavelengths of the optical signals received by the reception units 500 are different from each other. The configuration described heretofore is that in a case where the polarization division multiplexing and the wavelength division multiplexing are combined.

[0101] In a case of further combining the spatial division multiplexing, a point that the transmitter 10 transmits a spatially N multiplexed polarization multiplexed signal, a point that a device for spatial multiplexing/demultiplexing such as a mode multiplexer/demultiplexer is inserted in addition to the WDM multiplexer and the WDM demultiplexer, and a point that, in the receiver 50, as many optical front ends 520 as spatial multiplex number are arranged, the number of inputs of a MIMO equalizer (demodulation digital signal processing unit) and the number of complex impulse responses increase to 16N.sup.2, and N sets of polarization multiplexed signals are demodulated are added. The spatially N multiplexed polarization multiplexed signal is transmitted to the receiver 50 by, for example, a multi-core fiber a multi-mode or the like.

[0102] FIG. 6 is a diagram for explaining an outline of processing for deriving inverse characteristics of the transmitter 10 and the receiver 50 in the second embodiment. An adaptive equalization unit 534 performs adaptive equalization processing on 2N inputs of X.sub.1I, X.sub.1Q, . . . , X.sub.NI, and X.sub.NQ. The characteristic function derivation unit 538 receives an input of filter coefficients h.sub.1, . . . , and h.sub.4(N.sub.2.sub.) obtained in the process of the adaptive equalization processing by the adaptive equalization unit 534 and exp(j.sub.x1(n/T)), . . . , and exp(j.sub.xN(n/T)). Note that, herein, 1 in .sub.x1 is a subscript of x, and N in .sub.xN is a subscript of x. An inverse characteristic H.sub.T.sup.1 () of the transmitter 10 and an inverse characteristic H.sub.R.sup.1() of the receiver 50 are calculated by calculation processing performed by the characteristic function derivation unit 538.

[0103] FIG. 7 is a diagram illustrating a configuration example of a demodulation digital signal processing unit including the adaptive equalization unit 534 in the second embodiment. FIG. 7 illustrates a configuration of the adaptive equalization unit 534 in a case where 82 MIMO is extended a case of the multiplex number N to obtain a (4NN) MIMO configuration.

[0104] The demodulation digital signal processing unit sets an I component signal of an X-polarization component of a k-th (k is an integer of 1 or larger and N or smaller) polarization multiplexed reception signal output by the optical front end 520 as a real component X.sub.kI, a Q component signal as an imaginary component X.sub.kQ, an I component signal of a Y-polarization component as a real component Y.sub.kI, and a Q component signal as an imaginary component Y.sub.kQ. The demodulation digital signal processing unit convolutes an impulse response for compensating for a frequency characteristic of the receiver and a complex impulse response for wavelength dispersion compensation into each of the real component X.sub.kI, imaginary component X.sub.kQ, real component Y.sub.kI, and imaginary component Y.sub.kQ of the k-th polarization multiplexed reception signal according to each component.

[0105] The demodulation digital signal processing unit branches each of the convoluted real component X.sub.kI, imaginary component X.sub.kQ, real component Y.sub.kI, and imaginary component Y.sub.kQ into 4N signals. The demodulation digital signal processing unit directly inputs 2N signals out of the branched 4N signals into the adaptive equalization unit 534, and converts the remaining 2N signals into the phase conjugate signals to input to the adaptive equalization unit 534.

[0106] Phase conjugates of the real component X.sub.kI, the imaginary component X.sub.kQ, the real component Y.sub.kI, and the imaginary component Y.sub.kQ are set as a real component phase conjugate X.sub.kI*, an imaginary component phase conjugate X.sub.kQ*, a real component phase conjugate Y.sub.kI*, and an imaginary component phase conjugate Y.sub.kQ*, respectively. Each of 2N sets including the real component X.sub.kI, the imaginary component X.sub.kQ, the real component Y.sub.kI, the imaginary component Y.sub.kQ, the real component phase conjugate X.sub.kI*, the imaginary component phase conjugate X.sub.kQ*, the real component phase conjugate Y.sub.kI*, and the imaginary component phase conjugate Y.sub.kQ* corresponds to the X-polarization component and Y-polarization component of N polarization multiplexed reception signals.

[0107] The adaptive equalization unit 534 convolutes the impulse response into each of 2N real components X.sub.1I to X.sub.NI, imaginary components X.sub.1Q to X.sub.NQ, real components Y.sub.1I to Y.sub.NI, imaginary components Y.sub.1Q to Y.sub.NQ, real component phase conjugates X.sub.1I* to X.sub.NI*, imaginary component phase conjugates X.sub.1Q* to X.sub.NQ*, real component phase conjugates Y.sub.1I* to Y.sub.NI*, and imaginary component phase conjugates Y.sub.1Q* to Y.sub.NQ*. The adaptive equalization unit 534 adds the real components X.sub.1I to X.sub.NI, the imaginary components X.sub.1Q to X.sub.NQ, the real components Y.sub.1I to Y.sub.NI, and the imaginary components Y.sub.1Q to Y.sub.NQ into which the impulse response according to the polarization and each component is convoluted, for each polarization of each polarization multiplexed reception signal. The demodulation digital signal processing unit applies phase rotation for frequency offset compensation to the addition signal to generate a first addition signal.

[0108] Similarly, the adaptive equalization unit 534 adds the real component phase conjugates X.sub.1I* to X.sub.NI*, the imaginary component phase conjugates X.sub.1Q* to X.sub.NQ*, the real component phase conjugates Y.sub.1I* to Y.sub.NI*, and the imaginary component phase conjugates Y.sub.1Q* to Y.sub.NQ* into which the impulse response according to the polarization and each phase conjugate is convoluted, for each polarization of each polarization multiplexed reception signal. The demodulation digital signal processing unit applies phase rotation opposite to the phase rotation for frequency offset compensation to the addition signal to generate a second addition signal.

[0109] When obtaining a reception signal by adding the first addition signal and the second addition signal generated for the polarization for each polarization of each polarization multiplexed reception signal, the demodulation digital signal processing unit performs distortion correction by adding (or subtracting) a transmission data bias correction signal of the polarization.

[0110] The filter coefficients (h.sub.1, . . . , h.sub.4(N.sub.2.sub.) obtained in the processing of the adaptive equalization unit 534 in the demodulation digital signal processing unit and the frequency offsets (exp(j.sub.x1(n/T)), . . . , and (j.sub.xN(n/T))) are output to the characteristic function derivation unit 538.

[0111] In the second embodiment, the characteristic function derivation unit 538 illustrated in FIG. 8 derives the characteristic H.sub.T() of the transmitter 10 and the characteristic H.sub.R() of the receiver 50 on the basis of A().

[0112] The characteristic function derivation unit 538 calculates A() on the basis of following Expression (29) using a matrix M defined by following Expression (28).

[00018]$\begin{array}{cc}\left[\mathrm{Math}.28\right]& \\ M=\left(\begin{array}{ccccccc}1& i& 0& 0& .Math.& 0& 0\\ 1& -i& 0& 0& .Math.& 0& 0\\ 0& 0& 1& i& .Math.& 0& 0\\ 0& 0& 1& -1& .Math.& 0& 0\\ .Math.& .Math.& .Math.& .Math.& & .Math.& .Math.\\ 0& 0& 0& 0& .Math.& 1& i\\ 0& 0& 0& 0& .Math.& 1& -1\end{array}\right)& \mathrm{Expression}\left(28\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.29\right]& \\ A\left(\right)=M^{-1}\left(\begin{array}{cccc}{h}_1\left(\right)& h_{2N+1}\left(\right)& .Math.& h_{4N^2-2N+1}\left(\right)\\ {\tilde{h}}_2\left(\right)& {\tilde{h}}_{2N+2}\left(\right)& .Math.& {\tilde{h}}_{4N^2-2N+2}\left(\right)\\ .Math.& .Math.& & .Math.\\ {\tilde{h}}_{2N}\left(\right)& {\tilde{h}}_{4N}\left(\right)& .Math.& {\tilde{h}}_{4N^2}\left(\right)\end{array}\right)& \mathrm{Expression}\left(29\right)\end{array}$

[0113] The characteristic function derivation unit 538 calculates the inverse characteristic H.sub.R.sup.1() of the receiver 50 on the basis of following Expression (30).

[00019]$\begin{array}{cc}\left[\mathrm{Math}.30\right]& \\ \left(\begin{array}{ccc}{A}_{21}& .Math.& A_{22N-1}\\ .Math.& & .Math.\\ A_{2N1}& .Math.& A_{2N2N-1}\end{array}\right){\left(\begin{array}{ccc}{A}_{11}& .Math.& A_{12N-1}\\ .Math.& & .Math.\\ A_{2N-11}& .Math.& A_{2N-12N-1}\end{array}\right)}^{-1}& \mathrm{Expression}\left(30\right)\end{array}$

[0114] The characteristic function derivation unit 538 calculates the inverse characteristic H.sub.T.sup.1() of the transmitter 10 by performing a convolution operation on exp(j.sub.x1(n/T)), . . . , exp(j.sub.xN(n/T)) for frequency offset compensation on following Expression (31).

[00020]$\begin{array}{cc}\left[\mathrm{Math}.31\right]& \\ {\left(\begin{array}{ccc}{A}_{12}& .Math.& A_{12N}\\ .Math.& & .Math.\\ A_{2N-12}& .Math.& A_{2N-12N}\end{array}\right)}^{-1}\left(\begin{array}{ccc}{A}_{11}& .Math.& A_{12N-1}\\ .Math.& & .Math.\\ A_{2N-11}& .Math.& A_{2N-12N-1}\end{array}\right)& \mathrm{Expression}\left(31\right)\end{array}$

[0115] According to the receiver 50 in the second embodiment configured as described above, even in a case where the multiplex number is increased to N, it is possible to calculate characteristics of a transceiver including IQ crosstalk across a plurality of polarized waves.

Modification Example of Second Embodiment

[0116] As in the first embodiment, as the signal input to the adaptive equalization unit 534, a signal equivalent to the IQ signal and the phase conjugate signal of the IQ signal may be used.

Modification Example Common to First Embodiment and Second Embodiment

[0117] The adaptive equalization unit 534 and the characteristic function derivation unit 538 may be configured as a characteristic measurement device for measuring a characteristic function between lanes of the transmitter 10 and the receiver 50. In the above-described example, the configuration in which the characteristic measurement device is provided in the receiver 50 has been described, but the characteristic measurement device may be provided in a housing different from the receiver 50.

Third Embodiment

[0118] In a third embodiment, a configuration of compensating for signal waveform distortion on the basis of an inverse characteristic obtained by a characteristic function derivation unit by the method described in each of the above-described embodiments will be described. Various modifications and calculation processing described below are performed by the characteristic function derivation unit. In order to improve a signal quality using the inverse characteristic H.sub.T.sup.1() of the transmitter and the inverse characteristic H.sub.R.sup.1() of the receiver obtained in each of the above-described embodiments, a pre-equalization unit in the transmitter and a front end correction unit in the receiver have a function of performing (static) filtering processing on the signal in a corresponding format, and each inverse characteristic is input as a filter coefficient. The characteristic function derivation unit inputs the inverse characteristic H.sub.T.sup.1 () of the transmitter in the pre-equalization unit in the transmitter and the inverse characteristic H.sub.R.sup.1() of the receiver in the front end correction unit in the receiver as the filter coefficient.

[0119] Herein, the inverse characteristic obtained in each of the above-described embodiments has a format of a 44 matrix that acts on a signal vector in a form of (s.sub.x()()s.sub.x()s.sub.y()()s.sub.y()), and in order to compensate for the signal waveform distortion using this, a static filter configuration in a 44 format that acts on the signal vector in the form described above has been required.

[0120] However, in some signal processing circuits in the transmitter and the receiver, not a complex signal vector in the above-described format but an I component (real component) and a Q component (imaginary component) of the signal are basic elements of signal processing, so that there is a case where it is better to handle the signal vector as (Re(s.sub.x())Im(s.sub.x())Re(s.sub.y()) Im(s.sub.y())) instead of the above-described format. In this case, a characteristic function derivation unit 538 may input M.sup.1H.sup.1.sub.T/RM to the filter as a characteristic instead of H.sup.1.sub.T/R.

[0121] Hereinafter, signal processing using handling of this characteristic function is assumed. Hereinafter, a single mode fiber is assumed as a transmission medium, but a spatial multiplexing fiber is also similarly applicable. In order to further reduce an amount of calculation, as illustrated in a pre-equalization unit 116 and a front end correction unit 532 in FIG. 9, there is a case where it is desired to perform distortion compensation by using a plurality of smaller, for example, 22 format filters (for example, static filters 151-1 to 151-3, static filters 550-1 to 550-3) in combination. However, in this case, the inverse characteristic obtained by the above-described method cannot be directly applied to the signal processing. FIG. 9 is a diagram illustrating a configuration example of a digital coherent optical transmission system 1 in a third embodiment. The third embodiment is different from the first embodiment in configurations of the pre-equalization unit 116 and the front end correction unit 532. It is hereinafter described focusing on differences.

[0122] As an example of using a plurality of stages of filters in a 22 format (for example, static filters 151-1 to 151-3, static filters 550-1 to 550-3) in combination, it is considered to apply a filter including three static filters in a 22 format applicable to an XI lane and an XQ lane, the XQ lane and a YI lane, and the YI lane and a YQ lane, respectively, to the signal as illustrated in FIG. 10. At that time, an action of the filter on the signal expressed in a format of following Expression (32) is in a format of following Expression (33), and some components are always zero.

[00021]$\begin{array}{cc}\left[\mathrm{Math}.32\right]& \\ \left(\begin{array}{c}\mathrm{Re}\left(s_x\left(\right)\right)\\ \mathrm{Im}\left(s_x\left(\right)\right)\\ \mathrm{Re}\left(s_y\left(\right)\right)\\ \mathrm{Im}\left(s_y\left(\right)\right)\end{array}\right)& \mathrm{Expression}\left(32\right)\end{array}$$\begin{array}{cc}\left[\mathrm{Math}.33\right]& \\ \left(\begin{array}{cccc}{h}_{11}\left(\right)& h_{12}\left(\right)& 0& 0\\ h_{21}\left(\right)& h_{22}\left(\right)& h_{23}\left(\right)& h_{24}\left(\right)\\ h_{31}\left(\right)& h_{32}\left(\right)& h_{33}\left(\right)& h_{34}\left(\right)\\ 0& 0& h_{43}\left(\right)& h_{44}\left(\right)\end{array}\right)& \mathrm{Expression}\left(33\right)\end{array}$

[0123] However, a transceiver characteristic represented by M.sup.1H.sup.1.sub.T/RM cannot be handled by the filter of the configuration of this example because the above-described component is not zero in general. More generally, as for a characteristic matrix input to a filter for characteristic compensation, there is a case where only a certain (i,j) component must necessarily be zero. In the third embodiment, as described above, in a case where a plurality of static filters for compensating for a characteristic between some lanes is provided, and compensation is wanted to be performed by a combination thereof, a value to be input to the filter can be obtained from the inverse characteristic by an appropriate operation.

[0124] A method of modification will be specifically described below. First, a transmission side will be described below. Regarding the transceiver characteristic, even if a matrix representing any polarization rotation and phase rotation is multiplied on the transmission path side, this may be regarded as the same transceiver characteristic, so that as for any H.sub.c() in the form of following Expression (34), even if the inverse characteristic H.sup.1.sub.T of the transmitter is multiplied by H.sub.c() from the left and a function to be input to the filter (for example, the pre-equalization unit 116) is set to M.sup.1H.sub.c()H.sup.1.sub.TM, there is no change in compensation effect.

[00022]$\begin{array}{cc}\left[\mathrm{Math}.34\right]& \\ H_c\left(\right):=\left(\begin{array}{cccc}{h}_{c,11}\left(\right)& 0& h_{c,13}\left(\right)& 0\\ 0& {\tilde{h}}_{c,11}\left(\right)& 0& {\tilde{h}}_{c,13}\left(\right)\\ h_{c,31}\left(\right)& 0& h_{c,33}\left(\right)& 0\\ 0& {\tilde{h}}_{c,31}\left(\right)& 0& {\tilde{h}}_{c,33}\left(\right)\end{array}\right)& \mathrm{Expression}\left(34\right)\end{array}$

[0125] Therefore, H.sub.c() may be determined so as to satisfy a condition of a format of the filter that is wanted to be obtained. For this purpose, an equation of following Expression (35) may be solved for all (i,j) in which the (i,j) component of the characteristic matrix is zero.

[00023]$\begin{array}{cc}\left[\mathrm{Math}.35\right]& \\ {\left(M^{-1}{H}_c\left(\right)H_T^{-1}M\right)}_{ij}=0& \mathrm{Expression}\left(35\right)\end{array}$