Compensation coefficient calculation method

Abstract

In a method in which a compensation coefficient calculating portion (6) calculates a compensation coefficient of a compensation portion (5) which compensates transmission characteristics of a signal, a known signal is extracted from the signal. Next, a pseudo-random number is added to the extracted known signal. Next, the compensation coefficient is calculated by comparing a true value of the known signal with the known signal to which the pseudo-random number is added.

Claims

1. A compensation coefficient calculation method performed by a receiver receiving and converting an optical signal to a received signal which is an electrical signal, comprising: receiving the optical signal to which a known signal is added, converting the optical signal into the received signal, and extracting the known signal from the received signal; adding a pseudo-random number to the extracted known signal; and calculating a compensation coefficient of compensation circuitry configured to equalize transmission characteristics of the received signal by comparing a true value of the known signal with the known signal to which the pseudo-random number is added.

2. The compensation coefficient calculation method according to claim 1, wherein the known signal is a long cycle/continuous known pattern, the compensation circuitry further comprising an FIR filter, and the calculating the compensation coefficient further comprising calculating a tap coefficient of the FIR filter based on the long cycle/continuous known pattern by using a minimal mean square error method algorithm.

3. The compensation coefficient calculation method according to claim 2, further comprising: generating a pseudo-random number; and standardizing a magnitude of the pseudo-random number so that a ratio between average power of the known signal of the signal and average power of the pseudo-random number becomes a predetermined value, wherein the standardized pseudo-random number is added to the known signal in the adding the pseudo-random number to the extracted known signal.

4. The compensation coefficient calculation method according to claim 2, wherein the pseudo-random number has Gaussian distribution.

5. The compensation coefficient calculation method according to claim 2, wherein the compensation circuitry is configured to comprise a butterfly type filter, and the compensation coefficient is a tap coefficient of the butterfly type filter.

6. A compensation coefficient calculation method performed by a transmitter converting a transmission signal which is an electrical signal to an optical signal and transmitting the optical signal, comprising: extracting a known signal from an output signal of compensation circuitry configured to equalize transmission characteristics of the transmission signal to which the known signal is added; adding a pseudo-random number to the extracted known signal; and calculating a compensation coefficient of the compensation circuitry by comparing a true value of the known signal with the known signal to which the pseudo-random number is added.

7. The compensation coefficient calculation method according to claim 6, wherein the known signal is a long cycle/continuous known pattern, the compensation circuitry further comprising an FIR filter, and the calculating the compensation coefficient further comprising calculating a tap coefficient of the FIR filter based on the long cycle/continuous known pattern by using a minimal mean square error method algorithm.

8. The compensation coefficient calculation method according to claim 7, further comprising: generating a pseudo-random number; and standardizing a magnitude of the pseudo-random number so that a ratio between average power of the known signal of the signal and average power of the pseudo-random number becomes a predetermined value, wherein the standardized pseudo-random number is added to the known signal in the adding the pseudo-random number to the extracted known signal.

9. The compensation coefficient calculation method according to claim 7, wherein the pseudo-random number has Gaussian distribution.

10. The compensation coefficient calculation method according to claim 7, wherein the compensation portion is configured to comprise a butterfly type filter, and the compensation coefficient is a tap coefficient of the butterfly type filter.

11. A compensation coefficient calculation method performed by a receiver receiving and converting an optical signal to a received signal which is an electrical signal, comprising: a step of receiving the optical signal to which a known signal is added, converting the optical signal into the received signal, and extracting the known signal from the received signal; a pseudo-random number adding step of adding a pseudo-random number to the extracted known signal; and a coefficient calculating step of calculating a compensation coefficient of a compensation portion configured to equalize transmission characteristics of the received signal by comparing a true value of the known signal with the known signal to which the pseudo-random number is added.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a diagram showing a communication apparatus according to a first embodiment of the present invention.

(2) FIG. 2 is a diagram showing the adaptive equalization apparatus according to the first embodiment of the present invention.

(3) FIG. 3 is a diagram showing distribution of Z1.

(4) FIG. 4 is a diagram showing a case where pseudo-random numbers are not added to a received signal (a comparison example) and a case where pseudo-random numbers are added (the present embodiment).

(5) FIG. 5 is a diagram in which frequency characteristics are compared between a case where a tap coefficient is calculated in the method of the first embodiment of the present invention and a case where a tap coefficient is calculated in the method of the comparison example.

(6) FIG. 6 is a diagram showing an adaptive equalization apparatus according to a second embodiment of the present invention.

(7) FIG. 7 is a diagram showing a transmission apparatus according to a third embodiment of the present invention.

(8) FIG. 8 is a diagram showing a conventional adaptive equalization apparatus.

DESCRIPTION OF EMBODIMENTS

(9) A compensation coefficient calculation method according to the embodiments of the present invention will be described with reference to the drawings. The same components will be denoted by the same symbols, and the repeated description thereof may be omitted.

First Embodiment

(10) FIG. 1 is a diagram showing a communication apparatus according to a first embodiment of the present invention. This communication apparatus is a digital coherent optical receiver which receives an optical signal transmitted from the transmission side.

(11) A polarized wave separating portion 1 separates an optical signal into a horizontally polarized wave and a vertical polarized wave. A photoelectric conversion portion 2 converts the set of optical signals to analog electrical signals. An AD (Analog to Digital) conversion portion 3 converts the set of analog electrical signals to a set of digital electrical signals by sampling the set of analog electrical signals at a predetermined frequency. A digital signal processing portion 4 restores (demodulates) transmit data by performing digital signal processing for received signals which are the digital electrical signals outputted from the AD conversion portion 3. The digital signal processing portion 4 has an adaptive equalization apparatus described below.

(12) FIG. 2 is a diagram showing the adaptive equalization apparatus according to the first embodiment of the present invention. An adaptive filter 5 is a compensation portion which compensates for transmission characteristics of a received signal and includes an FIR filter. A compensation coefficient of this adaptive filter 5 is calculated by a compensation coefficient calculating portion 6. In the compensation coefficient calculating portion 6, a known signal extracting portion 7 extracts a known signal from a received signal. As the known signal, a continuous known pattern configured with a plurality of symbols in a relatively long cycle (long cycle/continuous known pattern; hereinafter referred to as an LP) is used. For example, as the LP, several hundred known symbols (for example, 128 symbols, 256 symbols, 512 symbols or the like) are added to a packet or an OTU (Optical-channel Transport Unit; tens of thousands of symbols) frame of a transmission signal at a rate of once or several times.

(13) A pseudo-random number generating portion 8 generates pseudo-random numbers Z(n) according to Gaussian distribution. Specifically, uniformly distributed random numbers are generated and then Gaussian-distributed using a Box-Muller method. For example, if and are assumed to be uniformly distributed random numbers, random numbers Z1 and Z2 according to Gaussian distribution are obtained from the following equations.

(14) Z1={square root over (2 log )} cos(2)
Z2={square root over (2 log )} sin(2)[Math. 1]
The uniformly distributed random numbers can be easily obtained by using random numbers corresponding to a predetermined number of bits from a pseudo-random bit sequence PRBS (for example, the 23rd or 31st row). Further, in the case of configuring the random numbers by software, a function of generating random data is included among tools of the software, and the function can be used.

(15) FIG. 3 is a diagram showing distribution of Z1. It is seen that Z1 is according to Gaussian distribution. Probability density distribution p(z) of Z1 and Z2 can be expressed by the following equation.

(16) $\begin{array}{cc}p\left(Z1,Z2\right)=\frac{1}{\sqrt{2}}\exp \left[-\frac{Z{1}^2+Z{2}^2}{2{}^2}\right]& \left[\mathrm{Math}.2\right]\end{array}$
Here, .sup.2 indicates variance of the Gaussian distribution. Here, when 0 and 1 are satisfied, p(Z1, Z2) is standard Gaussian distribution with a mean of 0 and a variance of 1.

(17) A standardization portion 9 regards pseudo-random numbers as noise and standardizes magnitude of the pseudo-random numbers Z(n) so that a set SN is obtained. That is, the standardization portion 9 specifies the variance .sup.2 of Gaussian distribution so that a ratio between average power of the known signal of the received signal and average power of the pseudo-random numbers Z(n) becomes an S/N (a signal-to-noise ratio) given from outside to standardize the magnitude of the pseudo-random numbers Z(n) as shown by the following equation.

(18) $\begin{array}{cc}\frac{S}{N}=\frac{\mathrm{average}\mathrm{power}\mathrm{of}\mathrm{known}\mathrm{signal}}{\mathrm{power}\mathrm{of}\mathrm{pseudo}\text{-}\mathrm{random}\mathrm{number}}=\frac{\mathrm{Ave}\left(\mathrm{strength}\mathrm{of}\mathrm{known}{\mathrm{signal}}^2\right)}{\mathrm{variance}\mathrm{of}\mathrm{Gaussian}\mathrm{distribution}{}^2}& \left[\mathrm{Math}.3\right]\end{array}$
Here, Ave (strength of known signal.sup.2) indicates a mean square value of the strength (voltage) of the known signal. Further, the average power of the Gaussian-distributed pseudo-random numbers is indicated by a square of standard deviation .

(19) In the case of the standard Gaussian distribution (the mean of 0 and the variance of .sup.2), standardized pseudo-random numbers NZ are expressed by the following equations in which each of Z1 and Z2 is multiplied by .
NZ1={square root over (2 log )} cos(2)
NZ2={square root over (2 log )} sin(2)[Math. 4]
At this time, power of each of the above pseudo-random numbers NZ is .sup.2. Here, for standardizing the pseudo-random numbers Z to be NZ is determined from the mean square value of the strength of the known signal and the S/N given from outside as shown by the following equation.
=[Ave(square of strength of known signal)/(S/N)].sup.(1/2)

(20) A pseudo-random number adding portion 10 adds the standardized pseudo-random numbers NZ(n) to the known signal which has been extracted. Either NZ1 or NZ2 can be used. Though values which have been converted from uniform distribution to Gaussian distribution are used as pseudo-random numbers in the above description, the distribution is not limited to Gaussian distribution, but similar distribution or the uniform distribution may be used.

(21) In a reference signal storing portion 11, a true value of an LP inserted into a data signal on the transmission side is stored as a reference signal in advance. The reference signal has a phase which the LP should originally take.

(22) A coefficient calculating portion 12 calculates the compensation coefficient of the adaptive filter 5 by comparing the known signal to which the pseudo-random numbers are added with the reference signal. Specifically, the coefficient calculating portion 12 calculates such a tap coefficient of the adaptive filter 5 that a mean square of an error e(n) between an output of the adaptive filter 5 and the reference signal becomes the smallest, using a minimal mean square error method algorithm. Similarly to a conventional MSE algorithm, the tap coefficient can be determined by the following equation.
W(n)=((X(n)+NZ).sup.T(X(n)+NZ)).sup.1(X(n)+NZ).sup.Td(n)
Here, W(n) indicates the tap coefficient of the adaptive filter 5; X(n) indicates a received signal; T indicates a transposed matrix; and d(n) indicates the reference signal.

(23) FIG. 4 is a diagram showing a case where pseudo-random numbers are not added to a received signal (a comparison example) and a case where pseudo-random numbers are added (the present embodiment). This figure shows a picture of adding pseudo-random numbers to a signal without noise for the purpose of description and does not necessarily accurately show a state where actual pseudo-random numbers are added. In the case of calculating a tap coefficient of a filter in an MSE algorithm as described above, it is possible to prevent divergence at the time of determining an inverse matrix, by adding pseudo-random numbers to a received signal. It is thought that this is because, in internal calculation, a value of a denominator is prevented from being very small at the time of performing division. Though other reasons are, of course, included, its effectiveness is confirmed in actual calculation.

(24) As described above, in the present embodiment, by adding pseudo-random numbers to a known signal, a minimum error from a true value of the known signal becomes finite, and divergence in calculation can be prevented. Therefore, it is possible to certainly calculate a compensation coefficient for compensating for transmission characteristics of a signal and improve equalization accuracy.

(25) FIG. 5 is a diagram in which frequency characteristics are compared between a case where a tap coefficient is calculated in the method of the first embodiment of the present invention and a case where a tap coefficient is calculated in the method of the comparison example. In standardization of pseudo-random numbers of the present embodiment, SNR=10 dB is assumed. Though an unnecessary high-frequency component equal to or above 20 GHz increases in the conventional method, the unnecessary high-frequency component equal to or above 20 GHz can be significantly reduced in the method of the present embodiment. A rate of reducing such a high-frequency component can be controlled by changing an SNR set for the standardization portion 9.

Second Embodiment

(26) FIG. 6 is a diagram showing an adaptive equalization apparatus according to a second embodiment of the present invention. In a case where an optical signal is separated into a horizontally polarized wave and a vertically polarized wave, and compensation is performed on each of the polarized waves by the adaptive filter 5, the adaptive filter 5 has a butterfly type filter to cope with influence by crosstalk between the polarized waves.

(27) The butterfly type filter is configured with four n-stage FIR filters. An output Y.sub.h of a horizontally polarized wave component of the adaptive filter 5 and an output Y.sub.v of a vertically polarized wave component are expressed by the following equations.

(28) $\begin{array}{cc}{Y}_h=\left[X\right]\left[H\right]=\left[\begin{array}{cc}\mathrm{H\_in}& \mathrm{V\_in}\end{array}\right]\left[\begin{array}{c}\mathrm{HH}\\ \mathrm{VH}\end{array}\right]& \left[\mathrm{Math}.5\right]\\ Y_v=\left[X\right]\left[V\right]=\left[\begin{array}{cc}\mathrm{H\_in}& \mathrm{V\_in}\end{array}\right]\left[\begin{array}{c}\mathrm{HV}\\ \mathrm{VV}\end{array}\right]& \end{array}$
Here, HH, VH, HV and VV indicate tap coefficients of the four FIR filters, respectively; H_in indicates an input of a horizontally polarized wave to the adaptive filter 5; and V_in indicates an input of a vertically polarized wave.

(29) Here, HH, VH, HV and VV are indicated by HH.sub.1 to HH, VH.sub.1 to VH.sub.n, HV.sub.1 to HV.sub.n, and VV.sub.1 to VV.sub.n, respectively; H_in is indicated by a sampling data string H_in.sub.1, H_in.sub.2, . . . , and V_in is indicated by a sampling data string V_in.sub.1, V_in.sub.2, . . . . Outputs Yh and Yv at the time when the sampling data strings are inputted to the butterfly type filter, respectively, are expressed by the following equations.

(30) $\begin{array}{cc}\left[\begin{array}{c}{Y}_{h1}\\ Y_{h2}\\ Y_{h3}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ Y_{hM}\end{array}\right]=[\begin{array}{cccccccccc}{\mathrm{H\_in}}_1& {\mathrm{H\_in}}_2& {\mathrm{H\_in}}_3& \mathrm{.Math.}& {\mathrm{H\_in}}_n& {\mathrm{V\_in}}_1& {\mathrm{V\_in}}_2& {\mathrm{V\_in}}_3& \mathrm{.Math.}& {\mathrm{V\_in}}_n\\ {\mathrm{H\_in}}_3& {\mathrm{H\_in}}_4& {\mathrm{H\_in}}_5& \mathrm{.Math.}& {\mathrm{H\_in}}_{n+2}& {\mathrm{V\_in}}_3& {\mathrm{V\_in}}_4& {\mathrm{V\_in}}_5& \mathrm{.Math.}& {\mathrm{V\_in}}_{n+2}\\ {\mathrm{H\_in}}_5& {\mathrm{H\_in}}_6& {\mathrm{H\_in}}_7& \mathrm{.Math.}& {\mathrm{H\_in}}_{n+4}& {\mathrm{V\_in}}_5& {\mathrm{V\_in}}_6& {\mathrm{V\_in}}_7& \mathrm{.Math.}& {\mathrm{V\_in}}_{n+4}\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& & \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& & \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{H\_in}}_{L-n+1}& {\mathrm{H\_in}}_{L-n+2}& {\mathrm{H\_in}}_{L-n+3}& \mathrm{.Math.}& {\mathrm{H\_in}}_L& {\mathrm{V\_in}}_{L-n+1}& {\mathrm{V\_in}}_{L-n+2}& {\mathrm{V\_in}}_{L-n+3}& \mathrm{.Math.}& {\mathrm{V\_in}}_L\end{array}]\left[\begin{array}{c}{\mathrm{HH}}_1\\ {\mathrm{HH}}_2\\ {\mathrm{HH}}_3\\ \mathrm{.Math.}\\ {\mathrm{HH}}_n\\ {\mathrm{VH}}_1\\ {\mathrm{VH}}_2\\ {\mathrm{VH}}_3\\ \mathrm{.Math.}\\ {\mathrm{VH}}_n\end{array}\right]& \left[\mathrm{Math}.6\right]\\ \left[\begin{array}{c}{Y}_{v1}\\ Y_{v2}\\ Y_{v3}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ Y_{\mathrm{vM}}\end{array}\right]=[\begin{array}{cccccccccc}{\mathrm{H\_in}}_1& {\mathrm{H\_in}}_2& {\mathrm{H\_in}}_3& \mathrm{.Math.}& {\mathrm{H\_in}}_n& {\mathrm{V\_in}}_1& {\mathrm{V\_in}}_2& {\mathrm{V\_in}}_3& \mathrm{.Math.}& {\mathrm{V\_in}}_n\\ {\mathrm{H\_in}}_3& {\mathrm{H\_in}}_4& {\mathrm{H\_in}}_5& \mathrm{.Math.}& {\mathrm{H\_in}}_{n+2}& {\mathrm{V\_in}}_3& {\mathrm{V\_in}}_4& {\mathrm{V\_in}}_5& \mathrm{.Math.}& {\mathrm{V\_in}}_{n+2}\\ {\mathrm{H\_in}}_5& {\mathrm{H\_in}}_6& {\mathrm{H\_in}}_7& \mathrm{.Math.}& {\mathrm{H\_in}}_{n+4}& {\mathrm{V\_in}}_5& {\mathrm{V\_in}}_6& {\mathrm{V\_in}}_7& \mathrm{.Math.}& {\mathrm{V\_in}}_{n+4}\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& & \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& & \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{H\_in}}_{L-n+1}& {\mathrm{H\_in}}_{L-n+2}& {\mathrm{H\_in}}_{L-n+3}& \mathrm{.Math.}& {\mathrm{H\_in}}_L& {\mathrm{V\_in}}_{L-n+1}& {\mathrm{V\_in}}_{L-n+2}& {\mathrm{V\_in}}_{L-n+3}& \mathrm{.Math.}& {\mathrm{V\_in}}_L\end{array}]\left[\begin{array}{c}{\mathrm{HV}}_1\\ {\mathrm{HV}}_2\\ {\mathrm{HV}}_3\\ \mathrm{.Math.}\\ {\mathrm{HV}}_n\\ {\mathrm{VV}}_1\\ {\mathrm{VV}}_2\\ {\mathrm{VV}}_3\\ \mathrm{.Math.}\\ {\mathrm{VV}}_n\end{array}\right]& \end{array}$

(31) In the case of performing sampling twice per symbol to acquire sampling data, data is set for the FIR filters for each sampling, but, as for timing of input to the filters, it is efficient to perform input per symbol. Therefore, at the next timing, data is inputted every two samplings.

(32) The known signal extracting portion 7 extracts an LP from each of the horizontally polarized wave input H_in and the vertically polarized wave input V_in. The pseudo-random number adding portion 10 adds standardized pseudo-random numbers NZ1(n) and NZ2(n) to the extracted LPs of the horizontally polarized wave input and the vertically polarized wave input, respectively. A reference signal LP.sub.h for an LP of a horizontally polarized wave and a reference signal LP.sub.v for an LP of a vertically polarized wave are stored in the reference signal storing portion 11 in advance.

(33) The coefficient calculating portion 12 determines such filter coefficients of the butterfly type filter that errors between received signals of the LPs to which the pseudo-random numbers are added, and the reference signals LP.sub.h and LP.sub.v, respectively, become the smallest. The filter coefficients are calculated by substituting Y.sub.h=LP.sub.h and Y.sub.v=LP.sub.v into the above equations. Therefore, [X][H]=[LP.sub.h] is calculated as shown below.

(34) $\begin{array}{cc}\left[\begin{array}{cccc}{\mathrm{H\_in}}_1+\mathrm{NZ}1\left(1\right)& {\mathrm{H\_in}}_2+\mathrm{NZ}2\left(1\right)& {\mathrm{H\_in}}_3+\mathrm{NZ}1\left(2\right)\mathrm{.Math.}& {\mathrm{H\_in}}_n+\mathrm{NZ}2\left(n/2\right)\\ {\mathrm{H\_in}}_3+\mathrm{NZ}1\left(2\right)& {\mathrm{H\_in}}_4+\mathrm{NZ}2\left(2\right)& {\mathrm{H\_in}}_5+\mathrm{NZ}1\left(3\right)\mathrm{.Math.}& {\mathrm{H\_in}}_{n+2}+\mathrm{NZ}2\left(n/2+1\right)\\ {\mathrm{H\_in}}_5+\mathrm{NZ}1\left(3\right)& {\mathrm{H\_in}}_6+\mathrm{NZ}2\left(3\right)& {\mathrm{H\_in}}_7+\mathrm{NZ}1\left(4\right)\mathrm{.Math.}& {\mathrm{H\_in}}_{n+4}+\mathrm{NZ}2\left(n/2+2\right)\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ \begin{array}{c}{\mathrm{H\_in}}_{L-n+1}+\\ \mathrm{NZ}1\left(\left(L-n+1\right)/2\right)\end{array}& \begin{array}{c}{\mathrm{H\_in}}_{L-n+2}+\\ \mathrm{NZ}2\left(\left(L-n+1\right)/2\right)\end{array}& \begin{array}{c}{\mathrm{H\_in}}_{L-n+3}+\\ \mathrm{NZ}1\left(\left(L-n+1\right)/2+1\right)\mathrm{.Math.}\end{array}& \begin{array}{c}{\mathrm{H\_in}}_L+\\ \mathrm{NZ}2\left(n/2+L-n\right)\end{array}\end{array}\begin{array}{cccc}{\mathrm{V\_in}}_1+\mathrm{NZ}1\left(1\right)& \begin{array}{c}{\mathrm{V\_in}}_2+\\ \mathrm{NZ}2\left(1\right)\end{array}& \begin{array}{c}{\mathrm{V\_in}}_3+\\ \mathrm{NZ}1\left(2\right)\mathrm{.Math.}\end{array}& \begin{array}{c}{\mathrm{V\_in}}_n+\\ \mathrm{NZ}2\left(n/2\right)\end{array}\\ {\mathrm{V\_in}}_3+\mathrm{NZ}1\left(2\right)& \begin{array}{c}{\mathrm{V\_in}}_4+\\ \mathrm{NZ}2\left(2\right)\end{array}& \begin{array}{c}{\mathrm{V\_in}}_5+\\ \mathrm{NZ}1\left(3\right)\mathrm{.Math.}\end{array}& \begin{array}{c}{\mathrm{V\_in}}_{n+2}+\\ \mathrm{NZ}2\left(n/2+1\right)\end{array}\\ {\mathrm{V\_in}}_5+\mathrm{NZ}1\left(3\right)& \begin{array}{c}{\mathrm{V\_in}}_6+\\ \mathrm{NZ}2\left(3\right)\end{array}& \begin{array}{c}{\mathrm{V\_in}}_7+\\ \mathrm{NZ}1\left(4\right)\mathrm{.Math.}\end{array}& \begin{array}{c}{\mathrm{V\_in}}_{n+4}+\\ \mathrm{NZ}2\left(n/2+2\right)\end{array}\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ \begin{array}{c}{\mathrm{V\_in}}_{L-n+1}+\\ \mathrm{NZ}1\left(\left(L-n+1\right)/2\right)\end{array}& \begin{array}{c}{\mathrm{V\_in}}_{L-n+2}+\\ \mathrm{NZ}2\left(\left(L-n+1\right)/2\right)\end{array}& \begin{array}{c}{\mathrm{V\_in}}_{L-n+3}+\\ \mathrm{NZ}1\left(\left(L-n+1\right)/2+1\right)\mathrm{.Math.}\end{array}& \begin{array}{c}{\mathrm{V\_in}}_L+\\ \mathrm{NZ}2\left(n/2+L-n\right)\end{array}\end{array}\right]\left[\begin{array}{c}{\mathrm{HH}}_1\\ {\mathrm{HH}}_2\\ {\mathrm{HH}}_3\\ \mathrm{.Math.}\\ {\mathrm{HH}}_n\\ {\mathrm{VH}}_1\\ {\mathrm{VH}}_2\\ {\mathrm{VH}}_3\\ \mathrm{.Math.}\\ {\mathrm{VH}}_n\end{array}\right]=\left[\begin{array}{c}{\mathrm{LP}}_{h1}\\ {\mathrm{LP}}_{h2}\\ {\mathrm{LP}}_{h3}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ \mathrm{.Math.}\\ {\mathrm{LP}}_{hM}\end{array}\right]& \left[\mathrm{Math}.7\right]\end{array}$
Here, [X][V]=[LP.sub.v] is also similarly calculated.

(35) Here, the standardized pseudo-random numbers NZ1 are added to odd-numbered data of the sample data string, and the standardized pseudo-random numbers NZ2 are added to even-numbered data of the sample data string. This is for the purpose of increasing randomness more and improving certainty of convergence and prevention of increase in a high-frequency range in the MSE algorithm. The above equations merely show an example of addition of the pseudo-random numbers NZ. If randomness is secured, a method is not limited to the above method, and it is not necessary to separately use NZ1 and NZ2. Further, variables may be arbitrary. Further, the pseudo-random numbers to be added may be either the same or different between a received signal of a horizontally polarized wave and a received signal of a vertically polarized wave.

(36) In the MSE algorithm of the present embodiment, mean square errors (MSE) are expressed by the following equations.
e.sub.h=X.Math.HLP.sub.h.sup.2, e.sub.v=X.Math.VLP.sub.v.sup.2[Math. 8]
Furthermore, the tap coefficients at the time when the mean square errors are the smallest are determined by the following equations.
H=(X.sup.HX).sup.1X.sup.HLP.sub.h, V=(X.sup.HX).sup.1X.sup.HLP.sub.v
Here, pseudo-random numbers are added to sample data of a received signal X.

(37) As described above, in the present embodiment, the adaptive filter 5 has a butterfly type filter, and compensation coefficients are tap coefficients of the butterfly type filter. Thus, in the case of applying the present invention to calculation of filter coefficients of an adaptive equalization filter for receiving double polarized waves also, it is possible to certainly calculate compensation coefficients for compensating for transmission characteristics and improve equalization accuracy similarly to the case of an adaptive equalization filter for receiving single polarized wave.

Third Embodiment

(38) FIG. 7 is a diagram showing a transmission apparatus according to a third embodiment of the present invention. The present embodiment is a configuration example in which the compensation coefficient calculation method of the present invention is applied to a transmission apparatus. A compensation portion 13 has a function of compensating for transmission characteristics of the transmission apparatus. An amplifier 14 amplifies an output of the compensation portion 13.

(39) An output of the amplifier 14 is fed back and inputted to the known signal extracting portion 7. Other components and operations are similar to those of the first and second embodiments in which the present invention is applied to a reception circuit. Thus, in the case of applying the present invention to a transmission circuit also, it is possible to certainly calculate a compensation coefficient for compensating for transmission characteristics and improve equalization accuracy without divergence in an MSE algorithm.

(40) A program for implementing the compensation coefficient calculation method according to the first to third embodiment may be recorded in a computer-readable recording medium so that the program recorded in the recording medium is read in and executed by a computer system or a programmable logic device to perform the compensation coefficient calculation. The computer system described here may contain hardware such as OS, peripheral equipment, etc. The computer system may contain a WWW system having a homepage providing environment (or displaying environment). The computer-readable recording medium may be a portable medium such as a flexible disc, a magneto-optic disc, ROM and CD-ROM, or a storage device such as a hard disc incorporated in a computer system. Furthermore, the computer-readable recording medium may contain an element for holding a program for a constant time like a volatile memory (RAM) in a computer system serving as a server or a client when a program is transmitted through a network such as the Internet or a communication line such as a telephone line. The program may be transmitted from a computer system having a storage device or the like in which the program is stored, to another computer system via a transmission medium or by transmission waves in a transmission medium. Here, the transmission medium for transmitting the program is a medium having a function of transmitting information like a network (communication network) such as the Internet or a communication circuit (communication line) such as a telephone line or the like. The program may be provided to perform a part of the foregoing function. Furthermore, the program may be provided to perform the foregoing function in combination with a program which has been already stored in the computer system, so-called a differential file (differentia program).

REFERENCE SIGNS LIST

(41) 5 adaptive filter; 6 compensation coefficient calculating portion; 7 known signal extracting portion; 8 pseudo-random number generating portion; 9 standardization portion; 10 pseudo-random number adding portion; 11 reference signal storing portion; 12 coefficient calculating portion; 13 compensation portion