Estimation apparatus and method for nonlinear distortion and receiver

Abstract

Embodiments of the present disclosure provide an estimation apparatus and method for nonlinear distortion and a receiver. The estimation method for nonlinear distortion includes: sampling a band-limited analog signal to obtain a sampling sequence; calculating a nonlinear perturbation coefficient in nonlinear distortion estimation based on a Nyquist pulse; calculating a nonlinear perturbation term superimposed on a signal by using the nonlinear perturbation coefficient and the sampling sequence; and calculating a nonlinear distortion waveform by using the nonlinear perturbation term. With the embodiments of the present disclosure, not only any modulation formats are compatible, but also advantages of high precision and good adaptability may be achieved.

Claims

1. An estimation apparatus for nonlinear distortion, configured in a receiving device of a communication system, the estimation apparatus comprising: a digital signal processor, the digital signal processor is configured to; sample a band-limited analog signal to obtain a sampling sequence; calculate a nonlinear perturbation coefficient in a nonlinear distortion estimation based on a Nyquist pulse; calculate a nonlinear perturbation term superimposed on a signal by using the nonlinear perturbation coefficient and the sampling sequence; and calculate a nonlinear distortion waveform by using the nonlinear perturbation term; wherein the processor is further configured to; perform integral processing on a nonlinear perturbation coefficient of each optical fiber span based on link parameters; perform large dispersion approximation on the integral processed nonlinear perturbation coefficient; simplify the large dispersion approximated nonlinear perturbation coefficient by using the Nyquist pulse, so as to obtain a nonlinear perturbation coefficient of each optical fiber span in a form of a single integral; and sum up respectively obtained nonlinear perturbation coefficients of different optical fiber spans, so as to obtain a nonlinear perturbation coefficient of a whole transmission link in a form of a single integral.

2. The estimation apparatus according to claim 1, wherein the processor is further configured to: perform approximation processing on one of a link loss and a gain function in an intra-channel nonlinear distortion estimation by using a rational function; and calculate the nonlinear perturbation coefficient in the form of the single integral by using the one of the link loss and gain function and calculate the nonlinear perturbation coefficient in the form of the single integral after ignoring the one of the ink loss and gain function to obtain the nonlinear perturbation coefficient of each optical fiber span in an analytic form, and sum up respectively the obtained nonlinear perturbation coefficients of different optical fiber spans to obtain the link nonlinear perturbation coefficient of a transmission link in the analytic form.

3. The estimation apparatus according to claim 1, wherein the processor is further configured to add an effect of nonlinear inter-symbol interference in calculating the nonlinear perturbation term.

4. The estimation apparatus according to claim 3, wherein an analytic expression of the nonlinear perturbation coefficient taking the nonlinear inter-symbol interference into account is: $\begin{array}{cc}{C}_{m,n}(\mathrm{lT})=\underset{i}{\overset{}{.Math.}}{C}_{m,n}^i(\mathrm{lT})=\frac{8}{9}\frac{L_{\mathrm{eff}}}{L_s}\{\begin{array}{c}\frac{1}{4{}^2.Math.{}_2.Math.{\mathrm{lTf}}_s^3}\underset{i}{\overset{}{.Math.}}[{}^{-\frac{jt}{T}}{E}_i\left(-j\frac{T\left[\mathrm{mnT}-\min .Math.\mathrm{lT}\right]}{C\left(z_i\right)}\right)-\\ {}^{\frac{jt}{T}}{E}_i\left(-j\frac{T\left[\mathrm{mnT}-\max .Math.\mathrm{lT}\right]}{C\left(z_i\right)}\right)]{|}_{z_i^0}^{L_S},\left(l0\right)\\ \underset{i}{\overset{}{j.Math.}}\left[\frac{1}{2.Math.{}_2.Math.f_s^2}{E}_i\left(-j\frac{T^2\mathrm{mn}}{C\left(z_i\right)}\right)+j\frac{\min -\max }{4{}^2.Math.{}_2.Math.f_s^4\mathrm{mn}}{}^{j\frac{{\mathrm{mnT}}^2}{C\left(z_i\right)}}\right]{|}_{z_i^0}^{L_S},\left(l=0\right)\end{array}& \end{array}$ where I is an integer, T is a symbol cycle, L.sub.eff=1/ is an effective transmission distance, is an attenuation coefficient of an optical fiber, m,n denote that three pulse locations of a nonlinear function are t=mT, nT and (m+n)T, respectively, min and max respectively denote min(m,n,m+n) and max(m,n,m+n), z.sub.i denotes a transmission distance of a signal at an i-th span, L.sub.s is a length of an optical fiber span, .sub.2 is a dispersion coefficient, is a residual dispersion rate, C(z.sub.i)=.sub.2(i1)L.sub.s+.sub.2z.sub.i is accumulative dispersion of the i-th span, and Ei denotes an exponential integral function.

5. A receiver, comprising the estimation apparatus for nonlinear distortion as claimed in claim 1.

6. An estimation method for nonlinear distortion, comprising the following steps performed by a receiving device of a communication system: sampling a band-limited analog signal to obtain a sampling sequence; calculating a nonlinear perturbation coefficient in a nonlinear distortion estimation based on a Nyquist pulse; calculating a nonlinear perturbation term superimposed on a signal by using the nonlinear perturbation coefficient and the sampling sequence; and calculating a nonlinear distortion waveform by using the nonlinear perturbation term; wherein the calculating a nonlinear perturbation coefficient in nonlinear distortion estimation based on a Nyquist pulse comprises the following steps performed by the receiving device; performing integral processing on a nonlinear perturbation coefficient of each optical fiber span based on link parameters; performing large dispersion approximation on the integral processed nonlinear perturbation coefficient; simplifying the large dispersion approximated nonlinear perturbation coefficient by using the Nyquist pulse, so as to obtain a nonlinear perturbation coefficient of each optical fiber span in a form of a single integral; and summing up respectively obtained nonlinear perturbation coefficients of different optical fiber spans, so as to obtain a nonlinear perturbation coefficient of a whole transmission link in a form of a single integral.

7. The method according to claim 6, wherein the calculating the nonlinear perturbation coefficient in the nonlinear distortion estimation based on the Nyquist pulse further comprises the following steps performed by the receiving device: performing approximation processing on one of link loss and a gain function in an intra-channel nonlinear distortion estimation by using a rational function; and one of calculating the nonlinear perturbation coefficient in the form of the single integral by using one of the link loss and gain function, and calculating the nonlinear perturbation coefficient in the form of the single integral after ignoring the one of the link loss and gain function to obtain the nonlinear perturbation coefficient of each optical fiber span in an analytic form, and summing up the respectively obtained nonlinear perturbation coefficients of different optical fiber spans to obtain the link nonlinear perturbation coefficient of a transmission link in the analytic form.

8. The method according to claim 6, wherein the method further comprises the following steps performed by the receiving device: adding an effect of nonlinear inter-symbol interference in calculating the nonlinear perturbation term.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The drawings are included to provide further understanding of the present disclosure, which constitute a part of the specification and illustrate the preferred embodiments of the present disclosure, and are used for setting forth the principles of the present disclosure together with the description. It is obvious that the accompanying drawings in the following description are some embodiments of the present disclosure only, and a person of ordinary skill in the art may obtain other drawings according to these accompanying drawings without making an inventive effort. In the drawings:

(2) FIG. 1 is a schematic diagram of the structure of a long-haul optical fiber transmission system;

(3) FIG. 2 is a flowchart of the estimation method of Embodiment 1 of the present disclosure;

(4) FIG. 3 is a schematic diagram of the structure of the nonlinear estimation model of any waveform of Embodiment 1 of the present disclosure;

(5) FIG. 4 is a schematic diagram of calculating a nonlinear perturbation coefficient of a Nyquist pulse shape of Embodiment 1 of the present disclosure;

(6) FIG. 5 is a schematic diagram of comparison of a Nyquist shadow pulse and a Gaussian shadow pulse of Embodiment 1 of the present disclosure;

(7) FIG. 6 is a schematic diagram of the structure of the estimation apparatus of Embodiment 2 of the present disclosure;

(8) FIG. 7 is a schematic diagram of the structure of the coefficient calculating unit of Embodiment 2 of the present disclosure; and

(9) FIG. 8 is a schematic diagram of the receiver of Embodiment 3 of the present disclosure.

DETAILED DESCRIPTION

(10) These and further aspects and features of the present disclosure will be apparent with reference to the following description and attached drawings. In the description and drawings, particular embodiments of the disclosure have been disclosed in detail as being indicative of some of the ways in which the principles of the disclosure may be employed, but it is understood that the disclosure is not limited correspondingly in scope. Rather, the disclosure includes all changes, modifications and equivalents coming within the terms of the appended claims.

(11) The embodiments of the present disclosure are applicable to a long-haul optical fiber communication system, and quantitative estimation is performed to nonlinear distortion of data signals in a transmission process. FIG. 1 is a schematic diagram of the structure of a long-haul optical fiber transmission system. As shown in FIG. 1, a signal transmitted by a transmitter reaches a receiver after passing through different devices in a transmission link (such as an optical fiber, an optical amplifier, and a dispersion compensation optical fiber, etc.). In the system shown in FIG. 1, nonlinear distortion estimation needs to be performed at the receiver end.

(12) Under a one-order perturbation framework, the embodiments of the present disclosure provide an estimation method and apparatus for nonlinear distortion applicable to any band-limited waveform. The method has an advantage of high precision and good adaptability, is applicable to a non-dispersion compensation link (NDM) and a dispersion management link (DM), in which any modulation formats are compatible; and it is also applicable to a single polarization or polarization multiplexing system. The estimation method is based on calculation of a nonlinear perturbation coefficient (which may also be referred to as a weighting coefficient) of a Nyquist pulse, and gives a method of calculation of nonlinear perturbation of a sampling sequence, thereby forming a set of nonlinear distortion estimation models of any waveform based on sampling, coefficient calculation, perturbation and nonlinear distortion calculation.

(13) Embodiment 1

(14) An embodiment of the present disclosure provides an estimation method for nonlinear distortion. FIG. 2 is a flowchart of the estimation method of the embodiment of the present disclosure. As shown in FIG. 2, the estimation method includes:

(15) step 201: sampling a band-limited analog signal to obtain a sampling sequence;

(16) step 202: calculating a nonlinear perturbation coefficient in nonlinear distortion estimation based on a Nyquist pulse;

(17) step 203: calculating a nonlinear perturbation term superimposed on a signal by using the nonlinear perturbation coefficient and the sampling sequence; and

(18) step 204: calculating a nonlinear distortion waveform by using the nonlinear perturbation term.

(19) In this embodiment, a nonlinear estimation model of any band-limited signal takes a sampling value of a signal waveform and link coefficients as input, and a nonlinear distortion waveform as output. It should be noted that the order of executing the steps in FIG. 2 is not limited thereto; for example, the order of executing steps 201 and 202 may be adjusted, or steps 201 and 202 may be executed at the same time, etc., according to an actual situation.

(20) FIG. 3 is a schematic diagram of the structure of the nonlinear estimation model of any waveform of the embodiment of the present disclosure. As shown in FIG. 3, the embodiment of the present disclosure is carried out by four functions; those are completed respectively: sampling of an analog signal, calculation of a nonlinear perturbation coefficient, calculation of a nonlinear perturbation term, and calculation of distorted waveform, which shall be described below in detail with reference to FIGS. 2 and 3.

(21) In step 201, any band-limited analog signal may be sampled. It is pointed out in the Nyquist sampling theorem that a band-limited analog signal u(t) (for example, spectrum occupies B+B) may be uniquely expressed by a sampling value of a sampling interval not greater than B. On the premise that the sampling theorem is satisfied, u(t) may be expressed as infinite series of an orthogonal sampling function (a sin c function), as shown in equation (1); where, T is a sampling period. Roots of the infinite series are sin c functions (Nyquist pulses) orthogonal to each other, and the coefficient is a sapling value u(kT) of u(t). Equation (1) shows that the band-limited signal may be recovered without distortion from its sampling value sequence by a matrix filter.
u(t)=.sub.k=.sup.+2B.Math.u(kTs).Math.sinc [2B(tkTs)](1)

(22) On the premise that the sampling theorem is not satisfied, i.e. a sampling frequency f.sub.s<2B, u(t) cannot be expressed by a sampling sequence without distortion. At this moment, expressing an analog signal by using a sampling sequence will introduce inevitable frequency domain aliasing distortion, the size of the distortion being related to a signal frequency and a sampling rate.

(23) In summary, no matter whether the aliasing is introduced, the band-limited analog signal may be expressed by a Nyquist pulse sequence. Considering from a perspective of applicability of a nonlinear model, if a theoretical model is capable of efficiently describing a nonlinear evolvement feature of the Nyquist pulse sequence, it may be applied to analysis of nonlinear distortion of any waveform based on the sampling theorem. A sampling rate of a pulse sequence is not limited in the embodiment of the present disclosure, and the theoretical analysis is applicable to cases with aliasing distortion or without aliasing distortion.

(24) In step 202, the nonlinear perturbation coefficient in the nonlinear distortion estimation may be calculated based on the Nyquist pulse; wherein, an optical fiber in the whole transmission link may be divided into multiple optical fiber spans, and for each optical fiber span, following processing may be performed: performing integral processing on a nonlinear perturbation coefficient of each optical fiber span based on link parameters; performing large dispersion approximation on the integral processed nonlinear perturbation coefficient; simplifying the large dispersion approximated nonlinear perturbation coefficient by using the Nyquist pulse, so as to obtain a nonlinear perturbation coefficient of each optical fiber span in a form of a single integral.

(25) In this embodiment, the respectively obtained nonlinear perturbation coefficients of different optical fiber spans may be summed up, so as to obtain a nonlinear perturbation coefficient of the whole transmission link in a form of a single integral.

(26) Furthermore, approximation processing may be performed on link loss or a gain function in intra-channel nonlinear distortion estimation by using a rational function, and the nonlinear perturbation coefficient in a form of a single integral is calculated by using the approximated link loss or gain function, so as to obtain a nonlinear perturbation coefficient of each optical fiber span in an analytic form; and the respectively obtained nonlinear perturbation coefficients of different optical fiber spans are summed up, so as to obtain a nonlinear perturbation coefficient of a whole transmission link in an analytic form.

(27) Furthermore, the nonlinear perturbation coefficient in a form of a single integral may be calculated after ignoring the link loss or gain function, so as to obtain a nonlinear perturbation coefficient of each optical fiber span in an analytic form; and the respectively obtained nonlinear perturbation coefficients of different optical fiber spans are summed up, so as to obtain a nonlinear perturbation coefficient of a whole transmission link in an analytic form.

(28) FIG. 4 is a schematic diagram of calculating a nonlinear perturbation coefficient of a Nyquist pulse shape of the embodiment of the present disclosure. Calculation of the nonlinear perturbation coefficient of the embodiment of the present disclosure shall be described below in detail with reference to FIG. 4.

(29) In a case where a transmission sequence pulse shape and link parameters are know and no approximation is taken into account, for an i-th optical fiber span, a one-order nonlinear perturbation weighting coefficient may be expressed as

(30) $\begin{array}{cc}{C}_{m,n}^i\left(t\right)=F^{-1}\left[{\tilde{C}}_{m,n}^i\left(f\right)\right]=j\frac{8}{9}{F}^{-1}\left\{{}_0^{L_s}{}^{-z_i-jC\left(z_i\right){\left(2f\right)}^2/2}F\left[g\left(t-\mathrm{mT},z_i\right).Math.g\left(t-\mathrm{nT},z_i\right).Math.g^{*}\left(t-\left(m+n\right)T,z_i\right)\right]z_i\right\}=j\frac{8}{9}{F}^{-1}\{{}_0^{L_s}{}^{-z_i-jC\left(z_i\right){\left(2f\right)}^2/2}\left[\tilde{g}\left(f,z_i\right){}^{j2\mathrm{fmT}}\right]*[\tilde{g}(f,z_i){}^{j2\mathrm{fnT}}*[\tilde{g}(-f,z_i){}^{j2f\left(m+n\right)T}]z_i\}& \left(2\right)\end{array}$

(31) Where, , are an attenuation coefficient and a nonlinear coefficient of an optical fiber, m, n denote that three pulse positions of a nonlinear action are t=mT, nT, (m+n)T, g(t, z) and {tilde over (g)}(f, z.sub.i) are a time domain expression and a frequency domain expression of a symbol pulse when it is transmitted to z.sub.i respectively, z.sub.i denoting a transmission distance of the signal at the i-th span, L.sub.s is a length of the optical fiber span, .sub.2 is a dispersion coefficient, is a residual dispersion rate, C(z.sub.i)=.sub.2(i1)L.sub.s+.sub.2z.sub.i is accumulated dispersion of the i-th span, F(.Math.) is Fourier transform, F.sup.1(.Math.) is inverse Fourier transform, and .sup.* denotes a convolution operation.

(32) It can be seen from Equation (2) that in a case of no approximation, the nonlinear perturbation coefficient C.sub.m,n.sup.i(t) may be uniformly expressed in a quadruple integral form (where, F.sup.1(.Math.) is a single integer, and the frequency domain two-dimensional convolution is a double integer), thereby performing calculation by using a numerical value method.

(33) As the calculation of the perturbation coefficient is the core of the estimation of the nonlinear distortion, the current study uniformly tends to how to simplify the quadruple integral shown in Equation (2), so that achieve nonlinear estimation of high precision by hardware of low complexity. A process of approximating Equation (2) into a single integral or an analytic solution in a specific case shall be described below.

(34) Large dispersion approximation is introduced first, and evolution of a symbol pulse under only an action of dispersion can be described approximately by using a far field.

(35) $\begin{array}{cc}\begin{array}{c}g\left(t,z\right)={}_{-}^{+}{}^{j\left[2\mathrm{ft}+2{}^{2}C\left(z_i\right)f^2\right]}\tilde{g}\left(f,0\right)f\\ \sqrt{\frac{1}{j2C\left(z_i\right)}}.Math.\tilde{g}\left(-\frac{t}{C\left(z_i\right)},0\right){}^{-j\frac{t^2}{2C\left(z_i\right)}}\end{array}& \left(3\right)\end{array}$

(36) The far field approximation may be physically interpreted as that in a case where dispersion is sufficiently large, a shape of a pulse is evolved into a Fourier transform of a transmission pulse, and long-haul transmission maps a pulse frequency domain shape to a time domain, with its principle being similar to the Fraunhofer far field diffractive integral, hence, it is referred to as far field approximation.

(37) The quadruple integral may be simplified into an expression of double integral by bringing Equation (3) into Equation (2):

(38) $\begin{array}{cc}{C}_{m,n}^i\left(t\right)=j\frac{8}{9}.Math.{}_{-}^{+}{}^{\mathrm{j2}\mathrm{ft}}f{}_0^{L_s}\frac{{}^{-z_i}}{2C\left(z_i\right)}\tilde{g}\left(f-\frac{\mathrm{mT}}{2C\left(z_i\right)},0\right)\tilde{g}\left(f-\frac{\mathrm{nT}}{2C\left(z_i\right)},0\right)\tilde{g}\left(-f+\frac{\left(m+n\right)T}{2C\left(z_i\right)},0\right){}^{j\frac{{\mathrm{mnT}}^2}{C\left(z_i\right)}}{z}_i& \left(4\right)\end{array}$

(39) Equation (4) shows that the two-dimensional convolution in large dispersion approximation may be approximated by multiplication.

(40) In taking a Nyquist shaping pulse (a roll-off coefficient is 0) of sampling points into account, Equation (4) may further be simplified. Taking the spectrum of the Nyquist pulse into account,

(41) $\begin{array}{cc}\tilde{g}\left(f,0\right)=\{\begin{array}{cc}1\text{/}{f}_s,& -f_s\text{/}2f{f}_s\text{/}2\\ 0,& \mathrm{others}\end{array}& \left(5\right)\end{array}$

(42) where, f.sub.s=1/T is a sampling rate. Bringing Equation (5) into Equation (4), an integral factor becomes into:

(43) $\begin{array}{cc}{\tilde{X}}_{m,n}\left(f,z_i\right)=\frac{{}^{-z_i}}{2C\left(z_i\right)}{{\tilde{g}}^{}}_{m,n}\left(f\right){}^{j\frac{{\mathrm{mnT}}^2}{C\left(z_i\right)}}& \left(6\right)\end{array}$

(44) where, overlapped spectrum {tilde over (g)}.sub.m,n(f) is denoted as:

(45) $\begin{array}{cc}{{\tilde{g}}^{}}_{m,n}\left(f\right)=\{\begin{array}{cc}1\text{/}{f}_s^3,& l_{m,n}f{r}_{m,n}\\ 0,& \mathrm{others}\end{array}& \left(7\right)\end{array}$
l.sub.m,n and r.sub.m,n decide upper and lower limits of the overlapped spectrum, l.sub.m,n=max(m, n, m+nT2Czifs2

(46) $\begin{array}{cc}{r}_{m,n}=\min \left(m,n,m+n\right)\frac{T}{2C\left(z_i\right)}+\frac{f_s}{2}& \left(8\right)\end{array}$

(47) As l.sub.m,n r.sub.m,n needs to be ensured, for determined m and n, the integral range is further constrained as:

(48) $\begin{array}{cc}{z}_i{z}_i^0\left(m,n\right)=\left[\max \left(m,n,m+n\right)-\min \left(m,n,m+n\right)\right]\frac{T^2}{2{}_2}-\left(-1\right)L_s& \left(9\right)\end{array}$

(49) Bringing Equation (6) into Equation (4) and performing inverse Fourier transform on rectangular overlapped pulses, a time domain nonlinear perturbation coefficient of the i-th span may be expressed as:

(50) $\begin{array}{cc}{C}_{m,n}^i\left(t\right)=j\frac{8}{9}{}_{z_i^0}^{L_s}\frac{{}^{-z_i}}{2C\left(z_i\right)f_s^3}{}^{j\frac{{\mathrm{mnT}}^2}{C\left(z_i\right)}}\left(r_{m,n}-l_{m,n}\right){}^{-j\left(r_{m,n}+l_{m,n}\right)t}\mathrm{sinc}\left[\left(r_{m,n}-l_{m,n}\right)t\right]z_i& \left(10\right)\end{array}$

(51) Equation (10) denotes that the nonlinear perturbation coefficient of the Nyquist shaping pulse may represent a one-dimensional integral form for the transmission distance. Such a one-dimensional integral may be calculated by using a numerical value method, and for the quadruple integral shown in Equation (2), the calculation complex may be outstandingly lowered. In the numerical value calculation, a perturbation coefficient of the whole transmission link may perform direct summation on the optical fiber spans, so as to obtain an expression of a perturbation coefficient in a one-dimensional form.

(52) 0$\begin{array}{cc}{C}_{m,n}\left(t\right)=\underset{i}{.Math.}{C}_{m,n}^i\left(t\right)=j\frac{8}{9}\underset{i}{\overset{}{.Math.}}{}_{z_i^0}^{L_s}\frac{{}^{-z_i}}{2C\left(z_i\right)f_s^3}{}^{j\frac{{\mathrm{mnT}}^2}{C\left(z_i\right)}}\left(r_{m,n}-l_{m,n}\right){}^{-j\left(r_{m,n}+l_{m,n}\right)t}\mathrm{sinc}\left[\left(r_{m,n}-l_{m,n}\right)t\right]z_i& \left(11\right)\end{array}$

(53) For Equation (11), a reason that there exists no analytic form in the one-dimensional integral is that there exists an exponential attenuation factor in an integral factor.

(54) For further obtaining an analytic solution of the perturbation coefficient, it may be taken into account to ignore the exponential attenuation, which is similar to a method for processing an attenuation factor by a Mecozzi model in a Gaussian pulse shape; wherein, let .sub.0.sup.Lse.sup.azi. f(zi)dzi1.Math.0Lsf(zi)dzi, and an analytic expression of Equation (12) may be obtained in such approximation:

(55) $\begin{array}{cc}{C}_{m,n}^i\left(t\right)=\{\begin{array}{c}\frac{8}{9}\frac{L_{\mathrm{eff}}}{L_s}\frac{1}{4{}^2.Math.{}_2.Math.{\mathrm{tf}}_s^3}\left[\begin{array}{c}{}^{-\frac{jt}{T}}{E}_i\left(-j\frac{T\left[\mathrm{mnT}-\min .Math.t\right]}{C\left(z_i\right)}\right)-\\ {}^{\frac{jt}{T}}{E}_i\left(-j\frac{T\left[\mathrm{mnT}-\max .Math.t\right]}{C\left(z_i\right)}\right)\end{array}\right]{|}_{z_i^0}^{L_S},\left(t0\right)\\ j\frac{8}{9}\frac{L_{\mathrm{eff}}}{L_s}\left[\frac{1}{2.Math.{}_2.Math.f_s^2}{E}_i\left(-j\frac{T^2\mathrm{mn}}{C\left(z_i\right)}\right)+j\frac{\min -\max }{4{}^2.Math.{}_2.Math.f_s^4\mathrm{mn}}{}^{j\frac{{\mathrm{mnT}}^2}{C\left(z_i\right)}}\right]{|}_{z_i^0}^{L_S},\left(t=0\right)\end{array}& \left(12\right)\end{array}$

(56) where, L.sub.eff=1/ is defined as an effective transmission distance, min, max denote min(m, n, m+n) and max(m, n, m+n), respectively. An analytic expression of a total perturbation coefficient may be obtained by performing summation on exponents of the spans in Equation (12).

(57) Another manner for processing attenuation factor is: replacing the exponential attenuation function with a specific rational function, so as to achieve a goal of making Equation (11) integrable. Such a processing method is better in estimation precision over the method of ignoring the attenuation, and is one of effective means for processing exponential attenuation.

(58) The calculation method of a perturbation coefficient of the Nyquist shaping pulse is given above, which is applicable to describing an evolvement process of nonlinear distortion of any waveform transmitted in any modulation format in the optical fiber link, does not limit the transmission link, and is a universal estimation method for nonlinear distortion.

(59) In step 203, the nonlinear perturbation coefficient and the sampling sequence may be used to calculate a nonlinear perturbation term superimposed on a signal. Wherein, an effect of nonlinear inter-symbol crosstalk is added in calculating the nonlinear perturbation term.

(60) In this embodiment, a physical meaning of C.sub.m,n(t) is time domain waveforms of shadow pulses produced at t=0 by Nyquist pulses located at three positions of t=mT, nT, (m+n)T, and a nonlinear distortion waveform of the whole pulse sequence is superimposition of shadow pulses at every time. Therefore, a characteristic of a shadow pulse decides a characteristic of nonlinear distortion to a certain extent. And it is necessary to further discuss a shadow pulse and a method for calculating a nonlinear distortion term.

(61) A nonlinear perturbation term is conventionally expressed in a form of a weighted sum of a product of three items of a pulse sequence, which is expressed as follows in a polarization multiplexing system:
u.sup.H/V(kT)=.sub.m,n[A.sub.m+k.sup.H/VA.sub.n+k.sup.H/V(A.sub.m+n+k.sup.H/V)*+A.sub.m+k.sup.H/VA.sub.n+k.sup.V/H(A.sub.m+n+k.sup.V/H)*]C.sub.m,n(0)(13);

(62) where, u.sup.H/V(kT) denotes nonlinear distortion at the H or V polarization state at a time of kT, and A.sub.m.sup.H/V denotes digital information of a transmission sequence at the H or V polarization state at a time of mT. This expression only uses information on the nonlinear perturbation coefficient at a time of t=0, with an implied assumption that a width of the nonlinear shadow pulse is sufficiently narrow, so that an effect of the shadow pulse on other sampling points may be ignored. However, such an assumption is applicable a Gaussian pulse sequence only, and for any waveform, an ideal sampling pulse is a Nyquist pulse, of which a characteristic of the nonlinear shadow pulse being different from a characteristic of a Gaussian shadow pulse, which is represented as that an expression of the nonlinear perturbation coefficient is different.

(63) FIG. 5 is a schematic diagram of comparison of a Nyquist shadow pulse and a Gaussian shadow pulse of the embodiment of the present disclosure. Without loss of generality, FIG. 5 shows shapes of shadow pulses produced at t=0 by two pairs of Gaussian pulses or Nyquist pulses with m=20T and n=40T. Simulation conditions are that a pulse period T=31.25 ps, pulse average power is 3 dBm, and the transmission link is a 100 km single-mode optical fiber with 50017 ps/nm dispersion.

(64) As shown in FIG. 5, the width of the Nyquist shadow pulse is obvious greater than that of the Gaussian shadow pulse, and at a time of t=T, the effect of the Gaussian shadow pulse may be ignored, but the effect of the Nyquist shadow pulse is remarkable. Therefore, the Nyquist pulse at t=(m+k)T, t=(n+k)T and t=(m+n+k)T not only produces nonlinear perturbation at t=kT, but also produces perturbation at t=(k1)T, (k2)T, . . . .

(65) In summary, it can be called that contribution of the Nyquist pulse at tkT located at t=(m+k)T, (n+k)T, (m+n+k)T is nonlinear symbol crosstalk (NL-ISI). Taking nonlinear distortion of any waveform constituted by the Nyquist sampling pulse into account, the embodiment of the present disclosure proposes that a nonlinear symbol crosstalk term should be taken into account in calculating the nonlinear perturbation term. In this case, the expression of the nonlinear perturbation term is changed into:
u.sup.H/V(kT)=.sub.l,m,n[A.sub.m+k+l.sup.H/VA.sub.n+k+l.sup.H/V(A.sub.m+n+k+l.sup.H/V)*+A.sub.m+k+l.sup.H/VA.sub.n+k+l.sup.V/H(A.sub.m+n+k+l.sup.V/H)*]C.sub.m,n(IT)(14)

(66) that is, in the dual-polarization state,

(67) $u^H\left(\mathrm{KT}\right)=\underset{l,m,n}{\overset{}{.Math.}}\left[A_{m+k+l}^H{A_{n+k+l}^H\left(A_{m+n+k+l}^H\right)}^{*}+A_{m+k+l}^H{A_{n+k+l}^V\left(A_{m+n+k+l}^V\right)}^{*}\right]C_{m,n}\left(\mathrm{lT}\right),u^V\left(\mathrm{KT}\right)=\underset{l,m,n}{\overset{}{.Math.}}\left[A_{m+k+l}^V{A_{n+k+l}^V\left(A_{m+n+k+l}^V\right)}^{*}+A_{m+k+l}^V{A_{n+k+l}^H\left(A_{m+n+k+l}^H\right)}^{*}\right]C_{m,n}\left(\mathrm{lT}\right);$

(68) where, I is an integer, T is a symbol period, A.sub.m+k+l.sup.H and A.sub.m+k+l.sup.V denote symbol information of a pulse at the horizontal polarization state and the vertical polarization state at a time of m+k+l, respectively, A.sub.n+k+l.sup.H and A.sub.n+k+l.sup.V denote symbol information of the pulse at the horizontal polarization state and the vertical polarization state at a time of n+k+l, respectively, and (A.sub.m+n+k+l.sup.H)* and (A.sub.m+n+k+l.sup.V)* denote conjugates of the symbol information of the pulse at the horizontal polarization state and the vertical polarization state at a time of m+n+k+l;

(69) or, in a single-polarization state, the expression of the nonlinear perturbation term is:

(70) $u\left(\mathrm{KT}\right)=\underset{l,m,n}{\overset{}{.Math.}}\left[A_{m+k+l}{A_{n+k+l}\left(A_{m+n+k+l}\right)}^{*}\right]C_{m,n}\left(\mathrm{lT}\right);$

(71) where, I is an integer, T is a symbol period, A.sub.m+k+l denotes symbol information of a pulse at the single-polarization state at a time of m+k+l, A.sub.n+k+l denotes symbol information of the pulse at the single-polarization state at a time of n+k+l, and (A.sub.m+n+k+l)* denotes conjugates of the symbol information of the pulse at the single-polarization state at a time of m+n+k+l.

(72) In this embodiment, an analytic expression of a perturbation coefficient taking nonlinear symbol crosstalk into account may be obtained by using Equation (12) to perform direction summation on the optical fiber spans:

(73) $\begin{array}{cc}{C}_{m,n}(\mathrm{lT})=\underset{i}{\overset{}{.Math.}}{C}_{m,n}^i(\mathrm{lT})=\frac{8}{9}\frac{L_{\mathrm{eff}}}{L_s}\{\begin{array}{c}\frac{1}{4{}^2.Math.{}_2.Math.{\mathrm{lTf}}_s^3}\underset{i}{\overset{}{.Math.}}[{}^{-\frac{jt}{T}}{E}_i\left(-j\frac{T\left[\mathrm{mnT}-\min .Math.\mathrm{lT}\right]}{C\left(z_i\right)}\right)-\\ {}^{\frac{jt}{T}}{E}_i\left(-j\frac{T\left[\mathrm{mnT}-\max .Math.\mathrm{lT}\right]}{C\left(z_i\right)}\right)]{|}_{z_i^0}^{L_S},\left(l0\right)\\ \underset{i}{\overset{}{j.Math.}}\left[\frac{1}{2.Math.{}_2.Math.f_s^2}{E}_i\left(-j\frac{T^2\mathrm{mn}}{C\left(z_i\right)}\right)+j\frac{\min -\max }{4{}^2.Math.{}_2.Math.f_s^4\mathrm{mn}}{}^{j\frac{{\mathrm{mnT}}^2}{C\left(z_i\right)}}\right]{|}_{z_i^0}^{L_S},\left(l=0\right)\end{array}& \left(14a\right)\end{array}$

(74) In Equation (14), the information on the nonlinear perturbation coefficient at other sampling points than 0 is used to calculate a nonlinear perturbation term, and adaptability is better and accuracy of the result is higher in using this equation for calculating the nonlinear perturbation. In using Equation (14) for calculating the nonlinear perturbation, as a range of values of I is limited, a truncation operation needs to be performed on I, and as shown in FIG. 5, the nonlinear inter-symbol crosstalk is relatively remarkable at T and 2T to which the shadow pulses are close, and its effect is relatively small at other locations. Therefore, a range of truncation of I is usually taken as 3I3.

(75) Calculation of the nonlinear distortion waveform in step 204 shall be described below. Actually, in the one-order nonlinear perturbation theorem, the nonlinear perturbation term shown in Equation (14) is completely deemed as additive perturbation, hence, superimposing the nonlinear perturbation term directly on the basis of the analytic solution may be deemed as a method for calculating a distortion waveform, that is,
u.sup.H/V(kT)=u.sub.0.sup.H/V(kT)+u.sup.H/V(kT)(15)
where, u.sub.0.sup.H/V(kT) is a linear solution of the nonlinear Schrodinger equation.

(76) Equation (15) is an expression of conventional calculation of nonlinear distortion waveform. Although Equation (15) is a theoretical solution of a one-order perturbation method, uniformly deeming nonlinear perturbation terms as additive terms is irrational, and such irrationality is actually originated from theoretical deviation of quasi-linear approximation of the one-order perturbation method when a nonlinear action is relatively intense.

(77) Further studies show that properly introducing high-order nonlinear terms may improve the precision of the estimation method and increase power adaptability range of the estimation method. A relatively typical result of this aspect is the previously-proposed nonlinear additive multiplicative (AM) model. In the additive multiplicative model, starting from the result produced by nonlinear pulse interaction, a difference between nonlinear perturbation terms related to three types of pulse interaction and two types of pulse interaction is discussed, pointing out that the nonlinear perturbation terms related to two types of pulse interaction are actually represented as nonlinear phase noises, thereby obtaining a method for calculating a nonlinear distortion waveform of higher application value, which is expressed as follows taking u.sup.H(kT) as an example:

(78) $\begin{array}{cc}{u}^H\left(\mathrm{kT}\right)=\left[u_0^H\left(\mathrm{kT}\right)+{}_{\mathrm{IFWM}}\left(\mathrm{kT}\right)\right]\exp \left[j_{\mathrm{IXPM}}\left(\mathrm{kT}\right)\right]& \left(16\right)\\ \mathrm{where},& \\ {}_{\mathrm{IFWM}}\left(\mathrm{kT}\right)=\underset{m0,n0}{\overset{}{.Math.}}\left[A_{m+k}^H{A_{n+k}^H\left(A_{m+n+k}^H\right)}^{*}+A_{m+k}^H{A_{m+k}^V\left(A_{m+n+k}^V\right)}^{*}\right]C_{m,n}\left(0\right)+\underset{m0,n}{\overset{}{.Math.}}{A}_{m+k}^H{A_{n+k}^V\left(A_{m+n+k}^V\right)}^{*}{C}_{m,n}\left(0\right)& \left(16a\right)\\ {}_{\mathrm{IXPM}}\left(\mathrm{kT}\right)={}_{m0}\left(2{.Math.A_{m+k}^H.Math.}^2+{.Math.A_{m+k}^V.Math.}^2\right)C_{m,0}\left(0\right)+\left({.Math.A_k^H.Math.}^2+{.Math.A_k^V.Math.}^2\right)C_{0,0}\left(0\right)& \left(16b\right)\end{array}$

(79) .sub.IFWM(kT) is an intra-channel quadruple mixing term, and .sub.IXPM(kT) is an intra-channel cross phase modulation term, corresponding respectively to nonlinear distortion introduced by three pieces of pulse interaction and two pieces of pulse interaction; wherein, the distortion of the three pieces of pulse interaction is still deemed as an additive term, and the distortion of the two pieces of pulse interaction is deemed as phase distortion, represented as a form of a product of the distortion and a nonlinear term. Equation (16) exhibits good precision of calculation in calculating a nonlinear distortion waveform of a Gaussian pulse.

(80) In a case of Nyquist pulse shape, as nonlinear intra-symbol crosstalk needs to be taken into account, the number of terms of nonlinear triple term product related to equations (16a) and (16b) is increased. Under a condition of accepting the physical meaning described by the additive multiplicative model, the expression of the additive multiplicative model under the Nyquist pulse needs to be modified, with a theoretical basis for the modification being that the terms of triple term product without participation of a local symbol are concluded as a .sub.IFWM(kT) term, and the terms of triple term product with participation of the local symbol are concluded as a .sub.IXPM(kT) term.

(81) Based on the above rule, the expressions of .sub.IFWM(kT) and .sub.IXPM(kT) in Equation (16) are modified as:
.sub.IFWM(kT)=.sub.m+k0,n+k0[A.sub.m+k.sup.HA.sub.n+k.sup.H(A.sub.m+n+k.sup.H)*+A.sub.m+k.sup.HA.sub.n+k.sup.V(A.sub.m+n+k.sup.V)*]C.sub.m,n(IT)+m+k0,nAm+kHAn+kVAm+n+kV*Cm,n(IT) (17a)
.sub.IXPM(kT)=.sub.m+k0(2|A.sub.m+k.sup.H|.sup.2+|A.sub.m+k.sup.V|.sup.2)C.sub.m,0(0)+(|A.sub.k.sup.H|.sup.2+|A.sub.k.sup.V|.sup.2)C.sub.0,0(0) (17b)

(82) The last term in the above two equations is a term of triple term product added taking the nonlinear inter-symbol crosstalk into account. Hence, in the Nyquist pulse shape, the nonlinear distortion waveform based on the additive multiplicative model is commonly described by equations (16), (17a) and (17b). When the nonlinear inter-symbol crosstalk is not taken into account, equations (17a) and (17b) are degenerated into equations (16a) and (16b).

(83) As calculating the perturbation coefficient by using Equation (11) or (14a) is relatively complex, in a general system (a dispersion compensation or dispersion management link), introduction of the method for calculating a C.sub.0,0(0) term under the Gaussian pulse assumption (a single integral expression) into Equation (17b) will not bring in many errors, hence, in actual perturbation calculation, the Gaussian pulse assumption is often used for simplification of the C.sub.0,0(0) term.

(84) It can be seen from the above embodiment that sampling is performed on a band-limited analog signal and the nonlinear distortion waveform is estimated by calculating the nonlinear perturbation coefficient based on the Nyquist pulse, in which not only any modulation formats are compatible, but also advantages of high precision and good adaptability may be achieved.

(85) Embodiment 2

(86) An embodiment of the present disclosure provides an estimation apparatus for nonlinear distortion, corresponding to the estimation method for nonlinear distortion described in Embodiment 1, with the identical contents being not going to be described any further.

(87) FIG. 6 is a schematic diagram of the structure of the estimation apparatus of the embodiment of the present disclosure. As shown in FIG. 6, the estimation apparatus 600 includes: a signal sampling unit 601, a coefficient calculating unit 602, a perturbation term calculating unit 603 and a waveform estimating unit 604.

(88) Wherein, the signal sampling unit 601 is configured to sample a band-limited analog signal to obtain a sampling sequence; the coefficient calculating unit 602 is configured to calculate a nonlinear perturbation coefficient in nonlinear distortion estimation based on a Nyquist pulse; the perturbation term calculating unit 603 is configured to calculate a nonlinear perturbation term superimposed on a signal by using the nonlinear perturbation coefficient and the sampling sequence; and the waveform estimating unit 604 is configured to calculate a nonlinear distortion waveform by using the nonlinear perturbation term.

(89) FIG. 7 is a schematic diagram of the structure of the coefficient calculating unit of the embodiment of the present disclosure. As shown in FIG. 7, the coefficient calculating unit 602 may include: an integral processing unit 701, a first approximating unit 702, a second approximating unit 703 and a coefficient summating unit 704.

(90) Wherein, the integral processing unit 701 is configured to perform integral processing on a nonlinear perturbation coefficient of each optical fiber span based on link parameters; the first approximating unit 702 is configured to perform large dispersion approximation on the integral processed nonlinear perturbation coefficient; the second approximating unit 703 is configured to simplify the large dispersion approximated nonlinear perturbation coefficient by using the Nyquist pulse, so as to obtain a nonlinear perturbation coefficient of each optical fiber span in a form of a single integral; and the coefficient summating unit 704 is configured to sum up respectively obtained nonlinear perturbation coefficients of different optical fiber spans, so as to obtain a nonlinear perturbation coefficient of a whole transmission link in a form of a single integral.

(91) As shown in FIG. 7, the coefficient calculating unit 602 may further include: a third approximating unit 705 and a fourth approximating unit 706. Wherein, the third approximating unit 705 is configured to perform approximation processing on link loss or a gain function in intra-channel nonlinear distortion estimation by using a rational function; and the fourth approximating unit 706 is configured to calculate the nonlinear perturbation coefficient in a form of a single integral by using the approximated link loss or gain function, or to calculate the nonlinear perturbation coefficient in a form of a single integral after ignoring the link loss or gain function, so as to obtain a nonlinear perturbation coefficient of each optical fiber span in an analytic form; and the coefficient summating unit 704 is further configured to sum up respectively obtained nonlinear perturbation coefficients of different optical fiber spans, so as to obtain a nonlinear perturbation coefficient of a whole transmission link in an analytic form.

(92) In this embodiment, the perturbation term calculating unit 603 adds an effect of nonlinear inter-symbol interference in calculating the nonlinear perturbation term. Wherein an analytic expression of the nonlinear perturbation coefficient taking the nonlinear inter-symbol interference into account is:

(93) $\begin{array}{cc}{C}_{m,n}(\mathrm{lT})=\underset{i}{\overset{}{.Math.}}{C}_{m,n}^i(\mathrm{lT})=\frac{8}{9}\frac{L_{\mathrm{eff}}}{L_s}\{\begin{array}{c}\frac{1}{4{}^2.Math.{}_2.Math.{\mathrm{lTf}}_s^3}\underset{i}{\overset{}{.Math.}}[{}^{-\frac{jt}{T}}{E}_i\left(-j\frac{T\left[\mathrm{mnT}-\min .Math.\mathrm{lT}\right]}{C\left(z_i\right)}\right)-\\ {}^{\frac{jt}{T}}{E}_i\left(-j\frac{T\left[\mathrm{mnT}-\max .Math.\mathrm{lT}\right]}{C\left(z_i\right)}\right)]{|}_{z_i^0}^{L_S},\left(l0\right)\\ \underset{i}{\overset{}{j.Math.}}\left[\frac{1}{2.Math.{}_2.Math.f_s^2}{E}_i\left(-j\frac{T^2\mathrm{mn}}{C\left(z_i\right)}\right)+j\frac{\min -\max }{4{}^2.Math.{}_2.Math.f_s^4\mathrm{mn}}{}^{j\frac{{\mathrm{mnT}}^2}{C\left(z_i\right)}}\right]{|}_{z_i^0}^{L_S},\left(l=0\right)\end{array}& \end{array}$

(94) where, I is an integer, T is a symbol cycle, L.sub.eff=1/ is an effective transmission distance, is an attenuation coefficient of an optical fiber, m, n denote that three pulse locations of a nonlinear function are t=mT, nT and (m+n)T, respectively, min and max respectively denote min(m, n, m+n) and max(m, n, m+n), z.sub.i denotes a transmission distance of a signal at an i-th span, L.sub.s is a length of an optical fiber span, .sub.2 is a dispersion coefficient, is a residual dispersion rate, C(z.sub.i)=.sub.2(i1)L.sub.s+.sub.2z.sub.i is accumulative dispersion of the i-th span, and Ei denotes an exponential integral function.

(95) It can be seen from the above embodiment that sampling is performed on a band-limited analog signal and the nonlinear distortion waveform is estimated by calculating the nonlinear perturbation coefficient based on the Nyquist pulse, in which not only any modulation formats are compatible, but also advantages of high precision and good adaptability may be achieved.

(96) Embodiment 3

(97) An embodiment of the present disclosure provides a receiver, including the estimation apparatus 600 for nonlinear distortion described in Embodiment 2.

(98) FIG. 8 is a schematic diagram of the receiver of the embodiment of the present disclosure. The receiver 800 may include a central processing unit (CPU) 200 and a memory 210, the memory 210 being coupled to the central processing unit 200. Wherein, the memory 210 may store various data; furthermore, it stores programs for information processing, and executes the programs under the control of the central processing unit 200.

(99) In an implementation, functions of the estimation apparatus 600 may be integrated into the central processing unit 200. Wherein, the central processing unit 200 may be configured to carry out the estimation method described in Embodiment 1.

(100) In another implementation, the estimation apparatus 600 and the central processing unit 200 may be configured separately. For example, the estimation apparatus 600 may be configured as a chip connected to the central processing unit 200, with the functions of the estimation apparatus 600 being realized under control of the central processing unit 200.

(101) Furthermore, as shown in FIG. 8, the receiver 800 may include a transceiver 220, and an antenna 230, etc. Wherein, the functions of the above components are similar to those of the prior art, and shall not be described herein any further.

(102) It should be noted that the receiver 800 does not necessarily include all the components shown in FIG. 8; furthermore, the receiver 800 may include components not shown in FIG. 8, and the prior art may be referred to for details.

(103) An embodiment of the present disclosure further provides a computer-readable program, wherein when the program is executed in an electronic device, the program enables a computer to carry out the estimation method for nonlinear distortion as described in Embodiment 1 in the electronic device.

(104) An embodiment of the present disclosure further provides a storage medium in which a computer-readable program is stored, wherein the computer-readable program enables a computer to carry out the estimation method for nonlinear distortion as described in Embodiment 1 in an electronic device.

(105) The above apparatuses and methods of the present disclosure may be implemented by hardware, or by hardware in combination with software. The present disclosure relates to such a computer-readable program that when the program is executed by a logic device, the logic device is enabled to carry out the apparatus or components as described above, or to carry out the methods or steps as described above. The present disclosure also relates to a storage medium for storing the above program, such as a hard disk, a floppy disk, a CD, a DVD, and a flash memory, etc.

(106) The present disclosure is described above with reference to particular embodiments. However, it should be understood by those skilled in the art that such a description is illustrative only, and not intended to limit the protection scope of the present disclosure. Various variants and modifications may be made by those skilled in the art according to the spirits and principle of the present disclosure, and such variants and modifications fall within the scope of the present disclosure.