Detailed Description
For the convenience of understanding the embodiments of the present invention, the following description will be further explained by taking several specific embodiments as examples in conjunction with the drawings, and the embodiments are not to be construed as limiting the embodiments of the present invention.
Example one
The processing flow of the method for compensating chromatic dispersion in an optical communication system provided by this embodiment is shown in fig. 1, and includes the following processing steps:
step 11, expressing the CD value of the signal to be compensated in the optical communication system by using a CD vector with a set length, where the value of the set length is the number of frequency points obtained after FFT conversion of the signal to be compensated, and the initial value of the CD data in the CD vector is the group velocity dispersion value of the signal to be compensated.
In a communication system, the same signal may be characterized by a signal in a time domain or a signal in a frequency domain, and thus, the processing of the signal in the communication system may be divided into time domain processing and frequency domain processing. In the present invention, dispersion compensation may be performed based on a signal in the frequency domain, or dispersion compensation may be performed based on a signal in the time domain.
When performing dispersion compensation based on the frequency domain signal, the step is performed based on the frequency domain digital signal obtained by transforming after FFT in the optical communication system, and specifically may include:
acquiring the number of frequency points obtained after FFT, expressing the CD value of the digital signal as a CD vector containing a plurality of CD data, wherein the number of the CD data is the number of the frequency points, the initial value of each CD data is the same, the initial value of each CD data is the group velocity dispersion value corresponding to the signal to be compensated, and the initial value of the CD vector is formed by utilizing the initial values of the CD data.
And dividing the CD data in the initial value of the CD vector into a plurality of groups, wherein the CD data in the same group are always kept equal.
When the dispersion compensation is performed on the time domain-based signal, the step is performed on the basis of the time domain signal in the optical communication system, and compared with the method for performing the dispersion compensation on the frequency domain-based signal, the method only needs to convert the time domain signal to be compensated into the frequency domain digital signal after serial-to-parallel conversion and FFT conversion, and the subsequent process is the same as the process described above.
And step 12, calculating an initial value of an error function of the CD vector according to the initial value of the CD data in the CD vector, changing the value of the CD vector for multiple times, calculating the value of the error function for multiple times respectively, and determining a final value of the CD vector according to the maximum value or the minimum value of the error function after obtaining the maximum value or the minimum value of the error function.
Calculating an initial value of a transfer function corresponding to the CD vector according to the initial value of the CD vector of the signal to be compensated, multiplying the digital signal by the initial value of the transfer function corresponding to the CD vector, and calculating an initial value of an error function corresponding to the CD vector according to the multiplied signal.
Selecting a first group of CD data in the CD vectors, changing the values of the first group of CD data, constructing a new CD vector according to the changed values of the first group of CD data, and recalculating the value of an error function corresponding to the new CD vector according to the value of the new CD vector; repeatedly executing the calculation process of the value of the error function, and after the maximum value or the minimum value of the error function is obtained, determining the final value of the first group of CD data according to the CD vector corresponding to the maximum value or the minimum value of the error function;
selecting CD data of other groups in the CD vector, determining the final values of the CD data of the other groups according to the determination process of the final value of the first group of CD data, and determining the final value of the CD vector according to the final values of all the groups of CD data.
And step 13, calculating a dispersion compensation function of the signal to be compensated according to the final value of the CD vector, and performing dispersion compensation on the signal to be compensated by using the compensation function.
Calculating a final value of a transfer function corresponding to the digital signal according to the final value of the CD vector, and taking the final value of the transfer function as a dispersion compensation function of the signal to be compensated;
and compensating the chromatic dispersion of the signal to be compensated according to the chromatic dispersion compensation function. The specific compensation process has two schemes: a scheme for performing compensation based on a signal in a frequency domain and a scheme for performing compensation based on a signal in a time domain.
In the scheme of compensating for a signal based on a frequency domain, compensating for the dispersion of the signal to be compensated according to the dispersion compensation function may specifically include: and multiplying the digital signal of the frequency domain obtained after the FFT by the dispersion compensation function in the optical communication system, and performing Inverse Fast Fourier Transform (IFFT) and parallel-serial conversion on the obtained digital signal to obtain a time domain signal and outputting the time domain signal. The output time domain signal is a signal with the chromatic dispersion compensated.
In the scheme of compensating for a signal based on a time domain, compensating for the dispersion of the signal to be compensated according to the dispersion compensation function may specifically include: and transforming the dispersion compensation function into a finite long single-bit impulse response filter FIR filter of a time domain through IFFT, performing online filtering compensation on a signal to be compensated of the time domain in an optical communication system through the FIR filter, and outputting the obtained time domain signal.
It can be seen from the technical solution provided in the above embodiment that, in the embodiment, the CD value of the signal to be compensated is represented by CD vectors with the same number of frequency points obtained after FFT conversion corresponding to the signal to be compensated, which is equivalent to setting one CD data for each frequency component of the signal to be compensated, so that the CD value compensation method can be adapted to scenes in which the CD values of the frequency components of the signal to be compensated are not equal or equal in an optical communication system, and particularly can effectively compensate dispersion ripples caused by various optical devices in the optical communication system.
Example two
The processing flow of the method for compensating chromatic dispersion in an optical communication system provided by this embodiment is shown in fig. 2, and includes the following processing steps:
and step 21, calculating the group velocity dispersion value of the signal to be compensated.
Inputting a time domain signal to be compensated, and converting the time domain signal into a frequency domain digital signal after serial-parallel conversion and FFT conversion in an optical communication system.
And then, performing group velocity dispersion estimation based on the digital signal to obtain a group velocity dispersion value of the signal to be compensated. The process of performing group velocity dispersion estimation on the signal to be compensated mainly comprises the following steps: converting the time domain signal to be compensated into the frequency domain, and multiplying the frequency domain signal by a dispersion compensation transfer function HCDH aboveCDThe CD value is calculated according to the trial value of the CD value of the digital signal, and the specific calculation formula is as follows:
HCD=exp(-if2CDλ2π/c)
in the above formula, f is the frequency of the digital signal, c is the speed of light, λ is the center wavelength of the digital signal, j is an imaginary unit, and CD is the trial value of the CD value of the digital signal.
Then according to the frequency domain signal and HCDMultiplying and calculating a feedback quantity according to the noise level of the multiplied signal. The trial value of the CD value is changed within a certain dispersion range, and the feedback amount is calculated. For example, the dispersion is about-1000 ps/nm to 30000ps/nm in the application scenario of the system, the dispersion range can be set to [ -5000, 40000 [)]. And acquiring a corresponding CD value according to the minimum value of the feedback quantity, and taking the CD value as the group velocity dispersion value of the digital signal.
And step 22, expressing the CD value of the signal to be compensated by a CD vector with a set length, wherein the set length is the number of frequency points obtained after FFT conversion is carried out on the signal to be compensated, the initial value of CD data in the CD vector is the group velocity dispersion value of the signal to be compensated, and the initial value of an error function is calculated according to the initial value of the CD vector.
In an actual optical communication system with dispersion ripple, the CD values of the frequency components of the signal to be compensated are not the same, and there is no correlation between the CD values of the frequency components, but they are approximately random, that is, there are CD values equal to the number of the frequency points.
And acquiring the number of frequency points obtained after the signal to be compensated is subjected to FFT. And expressing the CD value of the signal to be compensated by using a CD vector with a set length, wherein the value of the set length is the number of the frequency points. And if the number of the frequency points is N, the CD vector comprises N CD data, each CD data is equivalent to a frequency component corresponding to the signal to be compensated, the initial value of each CD data is the same, and the initial value of each CD data is the group velocity dispersion value.
The initial value of the CD vector can be determined according to the initial value of each piece of CD data, and the initial value of the transfer function corresponding to the CD vector is calculated according to the initial value of the CD vector.
The CD data in the CD vector are divided into a plurality of groups, the CD data in the same group have equal values, and the CD data in different groups may have equal or unequal values.
For example, if the number of frequency points is N equal to 8 points and the estimated group velocity dispersion value is 10, the initial value of the CD vector is: [ 1010101010101010 ] the length of the CD vector is 8.
The initial values of the CD vectors may be divided into 4 groups, which are: [10, 10], [10, 10] and [10, 10 ].
First, an error function corresponding to the CD vector at an initial value is calculated. Calculating an initial value of a transfer function of the CD vector using the CD vector [10, 10, 10, 10, 10, 10, 10] as an initial value according to the following formula:
HCD=exp(-if2CDλ2π/c)
in the above formula, f is the frequency of the digital signal, c is the speed of light, λ is the center wavelength of the digital signal, j is an imaginary unit, and CD is the initial value of the CD vector.
Then, transfer function H of the digital signal and the CD vector is carried outCDThe initial value of (2) is multiplied, and an initial error function error (1) corresponding to the digital signal is calculated from the multiplied signal, and the calculation principle of the error (1) may be calculated based on the noise level of the signal or the strength of the clock component.
And step 23, changing the values of a group of CD data in the CD vector for multiple times, calculating the values of the error function for multiple times, and determining the final value of the group of CD data according to the maximum value or the minimum value of the error function after the maximum value or the minimum value of the error function is obtained.
Selecting a set of CD data in the CD vector, and changing the value of the set of CD data. And constructing a new CD vector according to the changed values of the set of CD data. And according to the calculation process of the error function of the CD vector, recalculating the value of the error function corresponding to the new CD vector. And repeating the processing process, and determining the final value of the group of CD data according to the CD vector corresponding to the maximum value or the minimum value of the error function after the maximum value or the minimum value of the error function is obtained.
The process of determining the final value of the set of CD data will be described below by taking the minimum value of the error function as an example.
For example, the CD data of the first group of the CD vector [ 1010101010101010 ] is increased by 10, so that the CD vector becomes [20, 20, 10, 10, 10, 10, 10, 10], and the error function error (2) is calculated again according to the calculation process of the error function corresponding to the CD vector.
When error (1) > error (2), the CD data of the first group will continue to increase to 30, i.e. the CD vector is updated to [ 3030101010101010 ], and the error function error (3) is calculated again according to the calculation process of the error function corresponding to the CD vector.
Until the error function error (k +1) is obtained when the CD data of the first group is updated to [ (k +1) 10(k +1) 10101010101010 ], when error (k +1) > error (k), that is, when error (k-1) ≧ error (k) ≦ error (k +1), that is, when error (k) is a minimum value, the final value of the CD data of the first group is determined to be [ k 10 ].
Then, the next set of CD data in the CD vector is selected, and the same operation is performed on the next set of CD data until the final value of the next set of CD data is determined. After the final values of all sets of CD data in the CD vector are determined, the final values of the CD vector are determined.
And 24, calculating a dispersion compensation function corresponding to the digital signal according to the final value of the CD vector, and compensating the dispersion of the signal to be compensated by using the dispersion compensation function.
Calculating the transfer function H corresponding to the digital signal according to the final value of the CD vectorCDThe specific calculation formula of the final value of (a) is as follows:
HCD=exp(-if2CDλ2π/c)
in the above formula, f is the frequency of the digital signal, c is the speed of light, λ is the center wavelength of the digital signal, j is an imaginary unit, and CD is the final value of the CD vector.
The above transfer function HCDAs a dispersion compensation function for said digital signal.
And compensating the dispersion of the digital signal according to the dispersion compensation function of the digital signal. The digital signal is multiplied by a dispersion compensation function of the digital signal, and the obtained digital signal is subjected to IFFT (Inverse Fast Fourier Transform) conversion and parallel-to-serial conversion. And finally, outputting the obtained time domain signal.
It can be seen from the technical solution provided in the above embodiment that, in the embodiment, the CD value of the digital signal to be compensated is represented by a CD vector having the same length as the number of FFT points corresponding to the digital signal, which is equivalent to setting a CD data for each frequency component of the digital signal, so that the method can adapt to a scenario in which the CD values of the frequency components of the digital signal to be compensated in the optical communication system are not equal or equal, and particularly can effectively compensate dispersion ripples caused by various optical devices in the optical communication system, thereby improving the transmission performance and distance of the optical transmission system.
EXAMPLE III
In practical applications, for a scheme of performing dispersion on a signal in a time domain, after the dispersion compensation function h (f) of the signal to be compensated is calculated, the dispersion compensation function h (f) of the signal to be compensated may be transformed into a FIR (Finite Impulse Response filter) filter in the time domain by IFFT, where the FI R filter may be represented as h (t) ([ b0, b1, b 2., bN-1] (h (f)) and the FIR filter is used as a dispersion compensation function of the time domain signal as an input of the signal to be compensated. Then, the input time domain signal to be compensated is filtered and compensated by the FIR filter in the optical communication system, so that y [ N ] ═ b0x [ N ] + b1x [ N-1] + b2x [ N-2] +. + bN-1x [ N-N +1], where x [ N ], x [ N-1],. once, x [ N-N +1] are input time domain signals, and the obtained y [ N ] is taken as an output signal.
It can be seen from the technical solution provided in the above embodiment that the dispersion compensation part is not used for IFFT transformation and parallel-to-serial transformation of digital signals, and the speed of performing dispersion compensation on time-domain signals can be increased.
Example four
In practical applications, after the group velocity dispersion value of the digital signal is calculated, a group velocity dispersion compensation function may be calculated based on the group velocity dispersion value, and the digital signal may be subjected to dispersion compensation using the group velocity dispersion compensation function.
Then, after the dispersion compensation function of the digital signal is calculated according to the processing procedure of the first embodiment, the above-mentioned process of performing online dispersion compensation on the digital signal by using the group velocity color compensation function is stopped. And then the dispersion compensation function of the digital signal is utilized to carry out online dispersion compensation on the digital signal.
As can be seen from the technical solutions provided by the above embodiments, in the embodiments, before performing dispersion compensation on a digital signal by using a dispersion compensation function of the digital signal, the group velocity color compensation function is used to perform rapid dispersion compensation on the digital signal, so that the compensation speed and efficiency of the digital signal can be provided.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), or the like.
The embodiment of the present invention further provides a device for compensating chromatic dispersion in an optical communication system, the specific structure of which is shown in fig. 3, and the device comprises the following modules:
the CD vector establishing module 31 is configured to represent a CD value of a signal to be compensated by a CD vector with a set length, where the set length is a number of frequency points obtained after performing Fast Fourier Transform (FFT) on the signal to be compensated, initial values of CD data in the CD vectors are the same, and the initial value of each CD data is a group velocity dispersion value corresponding to the signal to be compensated;
a CD vector final value obtaining module 32, configured to calculate an initial value of an error function corresponding to the CD vector according to an initial value of CD data in the CD vector, change the value of the CD vector multiple times, calculate the value of the error function corresponding to the CD vector multiple times respectively, and determine a final value of the CD vector according to a maximum value or a minimum value of the error function;
and the dispersion compensation processing module 33 is configured to calculate a dispersion compensation function corresponding to the signal to be compensated according to the final value of the CD vector, and compensate the dispersion of the signal to be compensated by using the dispersion compensation function.
It can be seen from the technical solution provided in the above embodiment that, in the embodiment, the CD value of the signal to be compensated is represented by CD vectors with the same number of frequency points obtained after FFT conversion corresponding to the digital signal, which is equivalent to setting one CD data for each frequency component of the signal to be compensated, so that the CD value of each frequency component of the signal to be compensated in the optical communication system can be adapted to scenes with unequal or equal CD values, and dispersion ripples caused by various optical devices in the optical communication system can be compensated effectively.
The embodiment also provides another apparatus for compensating chromatic dispersion in an optical communication system, and the specific structure of the apparatus is shown in fig. 4, and the apparatus includes the following modules:
a CD vector establishing module 41, configured to represent a CD value of a signal to be compensated by a CD vector with a set length, where the set length is a number of frequency points obtained after performing Fast Fourier Transform (FFT) on the signal to be compensated, initial values of CD data in the CD vector are the same, and the initial value of each CD data is a group velocity dispersion value corresponding to the signal to be compensated;
a CD vector final value obtaining module 42, configured to calculate an initial value of an error function corresponding to the CD vector according to an initial value of CD data in the CD vector, change the value of the CD vector multiple times, calculate a value of the error function corresponding to the CD vector multiple times respectively, and determine a final value of the CD vector according to a maximum value or a minimum value of the error function;
and a dispersion compensation processing module 43, configured to calculate a dispersion compensation function corresponding to the signal to be compensated according to the final value of the CD vector, and compensate the dispersion of the signal to be compensated by using the dispersion compensation function.
Specifically, the CD vector establishing module 41 may specifically include:
a group velocity dispersion value obtaining module 411, configured to input a signal to be compensated, convert the signal to be compensated into a frequency domain digital signal through serial-to-parallel conversion and FFT, and perform group velocity dispersion estimation on the digital signal to obtain a group velocity dispersion value corresponding to the signal to be compensated;
a CD vector initial value obtaining module 412, configured to obtain the number of frequency points obtained after FFT, represent the CD value corresponding to the signal to be compensated as a CD vector including a plurality of CD data, where the number of CD data is the number of frequency points, an initial value of each CD data is a group velocity dispersion value corresponding to the signal to be compensated, and the initial values of the CD data are used to form the initial value of the CD vector;
and the packet processing module 413 is configured to divide the CD data in the initial value of the CD vector into a plurality of groups, where the CD data in the same group always remains equal.
Specifically, the CD vector final value obtaining module 42 may specifically include:
an error function calculation module 421, configured to calculate an initial value of a transfer function corresponding to the CD vector according to the initial value of the CD vector of the signal to be compensated, multiply the digital signal with the initial value of the transfer function corresponding to the CD vector, and calculate an initial value of an error function corresponding to the CD vector according to the multiplied signal, where a calculation principle of the error function may be based on a noise level of the signal or an intensity of a clock component;
the error function calculation module 421 selects a first set of CD data in the CD vectors, changes the values of the first set of CD data, constructs a new CD vector according to the changed values of the first set of CD data, and recalculates the value of the error function corresponding to the new CD vector according to the new values of the CD vector; repeatedly performing a calculation process of the value of the error function;
a final value determining module 422 for grouping CD data, configured to determine a final value of the first group of CD data according to a CD vector corresponding to a maximum value or a minimum value of the error function after the maximum value or the minimum value of the error function is obtained; selecting CD data of other groups in the CD vector, determining the final values of the CD data of the other groups according to the determination process of the final value of the first group of CD data, and determining the final value of the CD vector according to the final values of all the groups of CD data.
Specifically, the dispersion compensation processing module 43 may specifically include:
dispersion compensation function calculation modelAnd a block 431, configured to calculate a final value of a transfer function corresponding to the digital signal according to the final value of the CD vector, and use the final value of the transfer function as a dispersion compensation function of the signal to be compensated. Final value H of the transfer functionCDThe specific calculation formula of (A) is as follows:
HCD=exp(-if2CDλ2π/c)
in the above formula, f is the frequency of the digital signal, c is the speed of light, λ is the center wavelength of the digital signal, j is an imaginary unit, and CD is the final value of the CD vector.
And the compensation processing module 432 is configured to compensate the chromatic dispersion of the signal to be compensated according to the chromatic dispersion compensation function. The specific compensation process of the compensation processing module 432 includes two schemes: a scheme for performing compensation based on a signal in a frequency domain and a scheme for performing compensation based on a signal in a time domain.
In the scheme of performing compensation based on the frequency domain signal, the compensation process of the compensation processing module 432 may include: and multiplying the digital signal of the frequency domain obtained after the FFT by the dispersion compensation function in the optical communication system, and performing Inverse Fast Fourier Transform (IFFT) and parallel-serial conversion on the obtained digital signal to obtain a time domain signal and outputting the time domain signal. The output time domain signal is a signal with the chromatic dispersion compensated.
In the scheme of performing compensation based on the time domain signal, the compensation process of the compensation processing module 432 may include: and transforming the dispersion compensation function into a finite long single-bit impulse response filter FIR filter of a time domain through IFFT, performing online filtering compensation on a signal to be compensated of the time domain in an optical communication system through the FIR filter, and outputting the obtained time domain signal.
The dispersion compensation processing module 43 is further configured to calculate a group velocity color compensation function according to the group velocity dispersion value after calculating the group velocity dispersion value of the digital signal, and perform online dispersion compensation on the digital signal by using the group velocity color compensation function;
after the dispersion compensation function of the digital signal is calculated, the process of performing online dispersion compensation on the digital signal by using the group velocity color compensation function is terminated, and the online dispersion compensation on the digital signal is performed by using the dispersion compensation function of the digital signal.
In summary, it can be seen from the technical solution provided in the above embodiment that, in the embodiment of the present invention, a CD value of a signal to be compensated is represented by a CD vector having the same number as frequency points obtained after FFT conversion corresponding to the signal to be compensated, which is equivalent to setting a CD data for each frequency component of the signal to be compensated, so that the method can adapt to scenes in which CD values of frequency components of the signal to be compensated are not equal or equal in an optical communication system, effectively compensate dispersion ripples caused by various optical devices in the optical communication system, and improve transmission performance and distance of the optical transmission system.
When the dispersion compensation part of the embodiment of the invention is carried out in the time domain, the dispersion compensation part is not used for IFFT conversion and parallel-serial conversion of digital signals, and the speed of carrying out dispersion compensation on signals to be compensated can be improved.
In the embodiment of the invention, the group velocity color compensation function is firstly utilized to carry out rapid dispersion compensation on the digital signal before the online dispersion compensation is carried out on the digital signal by utilizing the dispersion compensation function of the digital signal, so that the compensation speed and efficiency of the signal to be compensated can be improved.
According to the embodiment of the invention, the transfer function of the digital signal is calculated off-line, and the transfer function of the digital signal can be calculated by adopting a high-complexity estimation algorithm, so that the calculation accuracy of the transfer function of the digital signal can be improved. The algorithm for estimating the transfer function of the Digital Signal can be implemented in a low-speed DSP (Digital Signal Processor) or a CPU, without increasing the logic scale of an ASIC (Application-Specific Integrated Circuit).
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.