WO2018054053A1 - Phase correction method and device, and computer storage medium - Google Patents

Phase correction method and device, and computer storage medium Download PDF

Info

Publication number
WO2018054053A1
WO2018054053A1 PCT/CN2017/081739 CN2017081739W WO2018054053A1 WO 2018054053 A1 WO2018054053 A1 WO 2018054053A1 CN 2017081739 W CN2017081739 W CN 2017081739W WO 2018054053 A1 WO2018054053 A1 WO 2018054053A1
Authority
WO
WIPO (PCT)
Prior art keywords
value
average
phase
blocks
phase offset
Prior art date
Application number
PCT/CN2017/081739
Other languages
French (fr)
Chinese (zh)
Inventor
邵晨峰
周晏
芦秋雁
Original Assignee
深圳市中兴微电子技术有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 深圳市中兴微电子技术有限公司 filed Critical 深圳市中兴微电子技术有限公司
Publication of WO2018054053A1 publication Critical patent/WO2018054053A1/en

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/25Arrangements specific to fibre transmission
    • H04B10/2507Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion
    • H04B10/2543Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion due to fibre non-linearities, e.g. Kerr effect
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • H04B10/613Coherent receivers including phase diversity, e.g., having in-phase and quadrature branches, as in QPSK coherent receivers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60Receivers
    • H04B10/61Coherent receivers
    • H04B10/616Details of the electronic signal processing in coherent optical receivers
    • H04B10/6165Estimation of the phase of the received optical signal, phase error estimation or phase error correction
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J14/00Optical multiplex systems
    • H04J14/02Wavelength-division multiplex systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04JMULTIPLEX COMMUNICATION
    • H04J14/00Optical multiplex systems
    • H04J14/02Wavelength-division multiplex systems
    • H04J14/0227Operation, administration, maintenance or provisioning [OAMP] of WDM networks, e.g. media access, routing or wavelength allocation
    • H04J14/0254Optical medium access
    • H04J14/0272Transmission of OAMP information
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes

Definitions

  • the present invention relates to the field of optical communications, and in particular, to a phase correction method and apparatus, and a computer storage medium.
  • optical fiber is an indispensable medium. How to transmit the most abundant information with the least amount of optical fiber. For this kind of exploration, the development of optical transmission has basically gone through the following stages: air separation The (SDM, Space Division Multiplexing) phase, the Time Division Multiplexing (TDM) phase, and the Wavelength Division Multiplexing (WDM) phase.
  • SDM Space Division Multiplexing
  • TDM Time Division Multiplexing
  • WDM Wavelength Division Multiplexing
  • the cable transmission is dominated by the WDM system.
  • the commercial 40G WDM transmission has gradually evolved to 100G and 400G transmission.
  • the data transmission distance is also constantly expanding.
  • the sub-system will bring about many problems such as chromatic dispersion, polarization film dispersion, strong filtering effect, frequency offset and other deviations in the transmission process.
  • phase deviation makes the constellation of the signal extend in the direction of the arc. If it is too large, it will cause the signal points to overlap in the middle, resulting in bit errors.
  • the mainstream phase-off phase correction method is four times, but the industry Most of the methods are calculated by using the data in the plane rectangular coordinate system. For example, the amplitude value of the sample is obtained by the data in the Cartesian coordinate system, and the sample is determined according to the amplitude value. According to the classification result, the sample is subjected to the quadratic phase deviation estimation in the Cartesian coordinate system according to the classification result, and then the phase deviation angle is calculated according to the phase deviation estimation value; although this calculation method is widely used, the method is applicable to each sample point. All calculations are carried out, which leads to a large amount of data calculation. The circuit structure used is complicated and fixed, and the circuit area is large and the power consumption is large. When the phase correction is performed, the processing efficiency is low.
  • embodiments of the present invention are expected to provide a phase correction method and apparatus, and a computer storage medium, which reduce the amount of calculation in performing phase correction and improve the processing efficiency of phase correction.
  • an embodiment of the present invention provides a phase correction method, including: performing M-th power processing on signal values of respective samples, and dividing signal values of each sample processed by the M-th power into at least one parallel segment.
  • M is a preset value; dividing a signal value of the sample point in each parallel segment to obtain a preset number of average blocks, and combining signal values of the sample points in each average block to obtain each of the average blocks a signal value, performing low-pass filtering processing on the signal value of each of the average blocks to obtain a low-pass filtered signal value of each of the average blocks; and according to the low-pass filtered signal value of each of the average blocks, Determining a phase offset value of each of the average blocks; obtaining, according to a phase offset value of each of the average blocks and a phase value of each of the samples in each of the average blocks, The corrected phase value of each sample.
  • the dividing the signal value of each sample processed by the Mth power into at least one parallel segment comprises: determining the size of the sliding window according to one or more of the following: transmission of signals of each sample point a laser line width of the signal of each sample; dividing the signal values of the samples processed by the M-th power to obtain the at least one parallel segment, wherein the number of samples of each parallel segment The size of the sliding window determined.
  • the combining the signal values of the samples in each average block to obtain the signal values of each of the average blocks includes: summing the signal values of the samples in each average block, The signal values of each of the average blocks are obtained.
  • the summation of the signal values of the samples in each average block is performed to obtain the signal value of each of the average blocks, including: the signal value of the sample points in each of the average blocks
  • the abscissa and the ordinate are respectively subjected to a summation calculation to obtain an abscissa of the signal value of each of the average blocks and an ordinate of the signal value of each of the average blocks; correspondingly, the pair of the averages
  • the signal value of the block is subjected to low-pass filtering processing to obtain a low-pass filtered signal value of each of the average blocks, including: an abscissa of the signal value of each of the average blocks and a signal of each of the average blocks
  • the ordinate of the value is subjected to low-pass filtering processing to obtain the abscissa of the low-pass filtered signal value of each of the average blocks and the ordinate of the low-pass filtered signal value of each of the average blocks.
  • determining, according to the low-pass filtered signal value of each average block, determining a phase offset value of each of the corresponding average blocks comprising: low-pass filtering each of the average blocks The subsequent signal value is converted into a low-pass filtered signal value of each average block in polar coordinates; the phase value corresponding to the low-pass filtered signal value of each average block in the polar coordinates is divided by M, The phase offset value of each of the average blocks is derived.
  • the method further comprising: in each of the parallel segments Determining, as a first difference, a difference between a phase offset value of each of the average blocks and a phase offset value of the corresponding previous average block; a preset according to the first difference A numerical interval that adjusts a phase offset value of each of the average blocks.
  • the method further includes: each parallel segment a difference between a phase offset value of the first average block and a phase offset value of a last average block in the corresponding previous parallel segment is determined as a second difference; falling according to the second difference
  • the preset value interval determines a correction value corresponding to the second difference; and the correction value is parallel to each of the Phase offset values of the first average block in the segment are summed to obtain a phase offset value of the adjusted first parallel block; and according to the modified value and the phase offset of the adjusted first parallel block
  • the obtaining, according to the phase offset value of each of the average blocks and the phase value of each sample in each of the average blocks, obtaining a correction for each sample in each of the average blocks a post-phase value comprising: a phase offset value for each of the average blocks, a phase value of each of the samples in each of the average blocks, and a - ⁇ /4 summation calculation to obtain each of the average blocks The corrected phase value for each sample.
  • an embodiment of the present invention provides a phase correction apparatus, including: a dividing module configured to perform a power processing on a signal value of each sample point, and divide a signal value of each sample processed by the M-th power processing To at least one parallel segment, wherein M is a preset value; the processing module is configured to divide the signal values of the samples in each parallel segment to obtain a preset number of average blocks, and the signals of the samples in each average block The values are combined to obtain the signal values of the average blocks, and the signal values of the average blocks are subjected to low-pass filtering processing to obtain low-pass filtered signal values of each of the average blocks; Determining, according to the low-pass filtered signal value of each of the average blocks, a phase offset value of each of the average blocks; and a correction module configured to: according to the phase offset value of each of the average blocks and each of the The phase values of each sample in the average block yield the corrected phase values for each of the samples in each of the average blocks.
  • the dividing module is further configured to determine a size of the sliding window according to one or more of: a transmission distance of a signal of each sample point, a laser line width of a signal of each sample point; and the M The signal values of the samples processed by the power are divided to obtain the at least one parallel segment, wherein the number of samples of each parallel segment is the size of the determined sliding window.
  • the processing module is further configured to perform a summation calculation on the signal values of the samples in each average block to obtain a signal value of each of the average blocks.
  • the processing module is further configured to sample the average block Calculating the abscissa and the ordinate of the signal value of the point respectively, obtaining the abscissa of the signal value of each of the average blocks and the ordinate of the signal value of each of the average blocks; respectively for each of the averages Performing low-pass filtering processing on the abscissa of the signal value of the block and the ordinate of the signal value of each of the average blocks, obtaining the abscissa of the low-pass filtered signal value of each of the average blocks and each of the average blocks The ordinate of the low pass filtered signal value.
  • the determining module is further configured to convert the low-pass filtered signal value of each average block into a low-pass filtered signal value of each average block in polar coordinates;
  • the phase value corresponding to the low-pass filtered signal value of each average block in polar coordinates is divided by M to obtain the phase offset value of each of the average blocks.
  • the apparatus further includes: a first adjustment module, configured to determine a phase offset value of each of the corresponding average blocks after the low-pass filtered signal value according to each of the average blocks Before, according to the phase offset value of each of the average blocks and the phase value of each sample in each of the average blocks, before obtaining the corrected phase values of each of the samples in each of the average blocks, In each parallel segment, determining a difference between a phase offset value of each of the average blocks and a phase offset value of the corresponding previous average block as a first difference; according to the first difference The preset value interval in which the value falls, adjusting the phase offset value of each of the average blocks.
  • a first adjustment module configured to determine a phase offset value of each of the corresponding average blocks after the low-pass filtered signal value according to each of the average blocks Before, according to the phase offset value of each of the average blocks and the phase value of each sample in each of the average blocks, before obtaining the corrected phase values of each of the samples in each of the average blocks, In each parallel segment, determining
  • the device further includes: a second adjustment module, configured to adjust a phase offset value of each of the average blocks after a preset value interval that falls according to the first difference, Determining, by the difference between the phase offset value of the first average block in each parallel segment and the phase offset value of the last average block in the corresponding previous parallel segment as a second difference; a preset value interval in which the second difference falls to determine a correction value corresponding to the second difference; summing the correction value and a phase offset value of the first average block in each parallel segment Obtaining a phase offset value of the adjusted first parallel block; adjusting the jth parallel in each parallel segment according to the modified value and the adjusted phase offset value of the first parallel block The phase offset value of the block, where j is greater than or equal to two.
  • the correction module is further configured to calculate a phase offset value of each of the average blocks, a phase value of each sample in each of the average blocks, and a - ⁇ /4 summation calculation, A corrected phase value for each of the samples in each of the average blocks is obtained.
  • an embodiment of the present invention provides a computer storage medium storing a computer program configured to perform the phase correction method.
  • the phase correction method and device and the computer storage medium provided by the embodiments of the present invention first perform M-th power processing on the signal values of the samples, and divide the signal values of the samples processed by the M-th power into at least one parallel segment. Then, the signal values of the samples in each parallel segment are divided to obtain a preset number of average blocks, and then the signal values of the samples in each average block are combined to obtain the signal value of each average block, thus avoiding For the defect of calculating the signal value of each sample, only a preset number of average blocks need to be low-pass filtered to determine the phase offset value of each average block, and finally, according to the phase of each average block
  • the offset value and the phase value of each sample in each average block yield a corrected phase value for each sample in each average block, that is, in the embodiment of the present invention, by sampling each parallel segment After the signal values of the points are divided and combined, the phase offset value is determined for each average block, which greatly reduces the calculation amount in determining the phase offset value, thereby reducing the phase correction process
  • FIG. 1 is a schematic flow chart of a phase correction method according to an embodiment of the present invention.
  • FIG. 3 is another optional flow chart of a phase correction method according to an embodiment of the present invention.
  • FIG. 4 is a schematic diagram of an optional addition tree structure according to an embodiment of the present invention.
  • FIG. 5 is an optional flow chart of adjusting a phase offset value according to an embodiment of the present invention.
  • FIG. 6 is another optional flow chart of adjusting a phase offset value according to an embodiment of the present invention.
  • FIG. 7 is a schematic structural diagram of a phase correcting apparatus according to an embodiment of the present invention.
  • the embodiment of the present invention provides a phase correction method, which can be applied to the correction of phase offset phase in an optical transmission network, and can also be applied to phase correction in other transmission networks, and the present invention is not specifically limited herein;
  • FIG. 1 is a schematic flowchart of a phase correction method according to an embodiment of the present invention. As shown in FIG. 1, the method includes:
  • S101 Perform M-th power processing on the signal values of the samples, and divide the signal values of the samples processed by the M-th power into at least one parallel segment;
  • M is a preset value, and the value of M may be 4, 8, and 12, etc., and the present invention is not specifically limited.
  • the parallel segment includes a fixed number of samples processed by the M-th power. Signal value, the fixed number is equal to the size of the sliding window described below;
  • FIG. 2 is an optional flow chart of the phase correction method in the embodiment of the present invention, such as As shown in Figure 2, the phase of the received signal can be expressed as follows:
  • ⁇ s (k) is the signal phase of the kth sample
  • ⁇ wkT i is the frequency offset phase
  • ⁇ n is the phase deviation phase caused by the laser line width
  • ⁇ ASE is the noise phase
  • V 4 (k) exp ⁇ j4 ⁇ s (k) ⁇ exp ⁇ j4 ⁇ n ⁇ exp ⁇ j4 ⁇ ASE ⁇ (2)
  • the phase value V 4 (k) after the fourth power processing further includes ⁇ s (k), ⁇ n and ⁇ ASE ; wherein, since ⁇ s (k) There are only a few fixed values: 0, + ⁇ /2, - ⁇ /2, and ⁇ , so V 4 (k) can remove ⁇ s (k), then the phase value V 4 (k) after the fourth power processing ) also includes ⁇ n and ⁇ ASE .
  • FIG. 3 is another optional flow chart of the phase correction method according to the embodiment of the present invention.
  • FIG. 3 after performing M-th power processing on the signal values of the samples, there are three values, which are rectangular coordinates.
  • S101 Can include:
  • the size of the sliding window is determined according to one or more of the following: the transmission distance of the signals of each sample point, and the laser line width of the signal of each sample point; the signal values of the samples after the M-th power processing are divided to obtain at least A parallel segment in which the number of samples of each parallel segment is the size of the determined sliding window.
  • the transmission distance of the signals of the above samples refers to the distance between the signals from the transmitting end to the receiving end, and the laser line width of the signals of the above various points is the laser line that generates the phase-biased phase. width;
  • the phase deviation phase of the symbol is regarded as the same premise is destroyed.
  • the phase deviation of the first sample and the last sample in the group is actually quite different, which affects the accuracy of the phase estimation, so it needs to be configured differently according to the scene.
  • the size of the sliding window is actually quite different, which affects the accuracy of the phase estimation, so it needs to be configured differently according to the scene. The size of the sliding window.
  • the method for adaptively adjusting the size of the sliding window is: after the system is powered on, the default value of the size of the sliding window (for example, the default value of the size of the sliding window is 4) is used for phase offset estimation and Compensate, and count the bit error rate under the size of this sliding window and record it in the software, then change the size of the sliding window through software, such as changing the size of the sliding window to 12, and counting the bit error rate under 12, if 12 If the error rate is lower than 4, the size of the sliding window is increased to 36. If the ratio of 36 is smaller than 12, the size of the sliding window is increased to 84.
  • the default value of the size of the sliding window for example, the default value of the size of the sliding window is 4
  • the error rate of 84 is larger than 36, then Select 36 as the size of the final sliding window; if the bit error rate under 84 is smaller than 36, continue to increase the size of the sliding window until the maximum supported by the hardware, if the error rate at the maximum is the smallest, select the maximum
  • the size of the sliding window the above principle of adjusting the size of the sliding window is to select the size of the sliding window with the smallest bit error rate by changing the size of the sliding window.
  • the size of the sliding window can be configured as four different values of 84, 36, 12, and 4.
  • the size of the sliding window is an integer multiple of 4, and the signal values of each sample are divided according to the size of the sliding window to obtain at least A parallel segment, the number of samples of each parallel segment is equal to the size of the sliding window, and is processed for each parallel segment instead of processing each sample, and the size of the sliding window is flexibly configured according to different application scenarios. Then, according to different scenarios, the size of different sliding windows is adaptively selected, so that the phase correction performance is optimal, and different addition tree structures are adopted according to the size of the sliding window, which reduces the amount of calculation, reduces resource consumption, and simplifies. design.
  • S102 Dividing the signal values of the samples in each parallel segment to obtain a preset number of average blocks, and combining the signal values of the samples in each average block to obtain signal values of each average block, for each average block.
  • the signal value is subjected to low-pass filtering processing to obtain a low-pass filtered signal value of each average block;
  • the two polarization states (in_i, in_q) in the Cartesian coordinate system are independently subjected to low-pass filtering processing, and the two do not affect each other, and the result after low-pass filtering processing is obtained.
  • each parallel segment dividing the signal values of the samples in each parallel segment to obtain a preset number of average blocks, where each average
  • the number of samples of the block is the sliding displacement of the window, and the sliding displacement of the window can be configured as 4, 8, 12, etc., wherein the preset number of values is the size of the sliding window divided by the sliding displacement of the window;
  • processing is performed for each average block.
  • the signal for each of the average block samples in S102 The values are combined to obtain the signal value of each average block, which may include: summing the signal values of the samples in each average block to obtain the signal value of each average block.
  • FIG. 4 is a schematic diagram of an optional addition tree structure in the embodiment of the present invention, as shown in FIG.
  • the addition branch of 12 is multiplexed with a window sliding displacement of 4, and the signal values of 12 samples are divided into 3 average blocks, and each 4 samples are divided into an average block, for 0-3.
  • the signal value of the sample is summed to obtain sum4_0, and the signal values of the 4th to 7th samples are summed to obtain sum4_1, and the signal values of the 8th to 11th samples are summed to obtain sum4_2;
  • the method is to perform multiplexing on different addition calculation paths.
  • the size of the sliding window can be configured as K1, K2, K3, K4, etc., wherein K1, K2, K3, K4, etc. are positive integers, for example, a sliding window.
  • the size can be configured as four different values of 84, 36, 12, and 4.
  • the size of the sliding window is an integer multiple of 4.
  • the multiplexed window has a sliding displacement of 4.
  • Structure, sliding window size 36, addition branch multiplexing window sliding The structure with a displacement of 12 and the addition branch with a sliding window size of 84 also multiplex the structure with a sliding displacement of 12; here, the addition branch of 84, the addition branch of 36 is similar to the addition branch of 12, and only the slip is given in the figure.
  • the addition branch of the window size is 12, and the phase offset value can be estimated by the addition tree structure, which further reduces the hardware implementation complexity and resource consumption.
  • the size of the sliding window is L, and the sliding displacement of the window is fixed by a value N, which can ensure the accuracy of the circuit is sufficient and can be simplified.
  • the parallel N-path summation is performed on the parallel samples, so that the number of sample points to be processed is reduced.
  • the processing method achieves the compression of the computational quantity under the condition that the processing precision requirement is satisfied, and the subsequent processing is performed based on the combined average block, which greatly reduces the complexity of the subsequent processing;
  • the scheme of merging the signal values of adjacent road samples has a generalization significance in similar designs.
  • L can take values of 64, 128, 256
  • N can take values of 4 and 8.
  • the summation of the signal values of the samples in each average block to obtain the signal value of each average block may include: summing the abscissa and the ordinate of the signal values of the samples in each average block respectively Obtaining the abscissa of the signal value of each average block and the ordinate of the signal value of each average block;
  • the signal value of each average block is subjected to low-pass filtering processing in S102 to obtain a low-pass filtered signal value of each average block, which may include: respectively, the abscissa and the signal value of each average block signal value.
  • the ordinate of the signal value of the average block is subjected to low-pass filtering processing to obtain the abscissa of the low-pass filtered signal value of each average block and the ordinate of the low-pass filtered signal value of each average block.
  • the received signal values are values in a Cartesian coordinate system, so the signal values of the samples include the abscissa and the ordinate.
  • the abscissa and the longitudinal direction of the signal values of the respective samples are respectively The coordinates are summed to obtain the abscissa of the signal value of each average block and the ordinate of the signal value of each average block;
  • the signal value of each average block is subjected to low-pass filtering processing, and the signal values obtained by low-pass filtering of each average block include denoising processing of the signal value of each average block and processing for eliminating residual frequency offset, That is to say, the low-pass filtering process can not only correct the phase deviation caused by the line width of the laser, but also eliminate a part of the residual frequency offset;
  • phase offset estimation and phase offset compensation are generally after frequency offset compensation and frequency offset estimation, that is, phase offset estimation and phase offset compensation.
  • the signal received in the reimbursement is the data after the majority of the frequency offset is removed, but there is still a partial residual frequency offset; since the frequency offset and the phase offset are both operated on the phase angle, in the embodiment of the present invention, the fourth power is passed.
  • the phase offset value estimated by the variable-pass low-pass filtering operation of the processing and sliding window is a real-time value. If there is a residual frequency offset, the estimated phase offset value will also change, and the phase offset value compensation is performed. It is equivalent to doing a frequency offset compensation, thereby eliminating the damage of the residual frequency offset to the signal.
  • the size of the sliding window can be configured to different sizes such as K1, K2, K3, K4, etc., and different sliding window sizes are adaptively adapted according to different residual frequency offsets.
  • the size of the larger sliding window is used to eliminate the residual frequency offset
  • the size of the smaller sliding window is used to eliminate the residual frequency offset.
  • the process of denoising by low-pass filtering is still taking the quadratic processing as an example.
  • the formula (2) is obtained.
  • V 4 (k) can remove the signal.
  • the phase is then added by multiple V 4 (k), that is, the noise phase is removed by low-pass filtering, and then the angle is extracted to obtain the phase estimation result.
  • the obtained phase estimation value ⁇ n is expressed as follows:
  • phase offset value ⁇ e is optimized to obtain the phase offset value ⁇ e as:
  • a variable-size sliding window is used to average the signal values of the samples in the window, which is done by summing the signal values of the samples in the sliding window, in order to balance the processing.
  • the invention provides a fixed window sliding displacement and supports the accuracy, the adaptability of different application scenarios, the complexity of the circuit, the reduction of the area and the power consumption, and the speed of the processing.
  • Several different sliding window sizes, maximizing the degree of logical multiplexing, compared to traditional unoptimized practices The method can ensure the accuracy of the data, and significantly reduce the circuit area by more than 50%, improve the processing speed of the circuit, and significantly reduce the power consumption.
  • S103 Determine a phase offset value of each average block according to the low-pass filtered signal value of each average block.
  • the obtained low-pass filtered signal value of each average block includes a phase offset value which is a value in a Cartesian coordinate system.
  • the method may include:
  • the low-pass filtered signal value of each average block is converted into a low-pass filtered signal value of each average block in polar coordinates; the phase corresponding to the low-pass filtered signal value of each average block in polar coordinates Dividing the value by M yields the phase offset value for each average block.
  • CORDIC CoDIC Rotation Digital Computer
  • CORDIC can convert the Cartesian coordinate system into a polar coordinate system, so that each average block low pass in the polar coordinate system can be obtained.
  • the phase value corresponding to the filtered signal value since the M-th power processing is performed for each sample point, the phase value corresponding to the low-pass filtered signal value of each average block is divided by M, and the phase offset value of each average block is obtained.
  • the low-pass filter is used to denoise and eliminate the residual frequency offset, and the amplitude is obtained by CORDIC, and the obtained amplitude value is divided by 4.
  • the phase offset value can be obtained.
  • the CORDIC processing is added to make the multiplication operation in the original Cartesian coordinate system become the addition operation in the polar coordinate system, and the calculation amount is greatly reduced.
  • the method may include:
  • a difference between a phase offset value of each average block and a phase offset value of the corresponding previous average block is determined as a first difference; according to the first difference
  • the preset value interval adjusts the phase offset value of each average block.
  • theta1 is first adjusted intra-segment; in the specific implementation process, the preset numerical interval can be divided into five preset numerical intervals by several thresholds; for each average block, two before and after The difference between the phase offset values of the average blocks is compared with several thresholds, and different adjustment operations are performed according to the comparison results to achieve fine adjustment, and the adjustments in the parallel segments are as follows:
  • the first average block in the segment is not adjusted.
  • the phase offset value of each average block is represented by PE k
  • the adjusted phase offset value in each average block segment is PE k ' ⁇ [- ⁇ /4, ⁇ /4) indicates that in order to make the adjusted phase offset value in each average block segment use PE k ' and the previous average block adjusted phase offset estimate PE k-1 ' ⁇
  • the deviation between [- ⁇ , ⁇ ) does not exceed ⁇ /4.
  • FIG. 5 is an optional flow chart of adjusting the phase offset value in the embodiment of the present invention, as shown in FIG. 5, and the specific adjustment method is as follows: :
  • PE' k PE k + ⁇
  • PE' k PE k + ⁇ /2;
  • PE' k PE k ;
  • PE' k PE k - ⁇ /2;
  • PE' k PE k - ⁇ .
  • the difference between PE k and PE k-1 ' is firstly performed, and the method branch of the intra-segment adjustment is obtained according to the range of the difference and the range of PE k , and then the corresponding intra-segment adjustment result PE is selected according to the obtained method branch.
  • k ' the input of the structure diagram in Fig. 5 is the phase offset value, and the output is the phase offset value PE k ' after the segment adjustment.
  • phase offset value In order to further adjust the phase offset value to obtain a more accurate phase offset value, After the parallel segment adjustment is completed, inter-segment adjustment is also required. Then, in an optional embodiment, the phase of each average block is adjusted according to a preset value interval falling according to the first difference.
  • the method may include: determining a difference between a phase offset value of the first average block in each parallel segment and a phase offset value of a last average block in the corresponding previous parallel segment as a second difference; determining a correction value according to a preset value interval in which the second difference falls; summing the correction value and a phase offset value of the first average block in each parallel segment to obtain an adjusted a phase offset value of a parallel block; adjusting a phase offset value of the jth parallel block in each parallel segment according to the modified value and the adjusted phase offset value of the first parallel block, where j is greater than or equal to 2.
  • FIG. 3 after the intra-segment adjustment of theta1, the inter-segment adjustment is performed to obtain theta2; in the specific implementation process, the inter-segment adjustment can be refined into two steps, and FIG. 6 is a phase shift in the embodiment of the present invention.
  • Another optional flow diagram for adjusting the values, as shown in Figure 6, is the first step of the correction of the first average block of the parallel segment: as with the intra-segment adjustment, it will not be described here, and is adjusted from the above segment. It can be seen that the range of the angle value of the first average block correction of the parallel segment (corresponding to the above correction value) is as follows:
  • the second step is the quadratic correction in the parallel segment: according to the correction of the first average block, the average block of the latter average block from the second one is corrected again, wherein PE k ' is the intra-segment adjustment After the sample, PE kf ' is the result of the adjustment between the segments, the specific method is as follows:
  • the above method branches 1) and 7) are combined, and 2) and 8) are combined to obtain the first step and the inter-segment adjustment.
  • the difference between PE k and PE k-1 ' is first obtained to obtain a second difference, and the corrected angle value of the first average block is obtained according to the range of the second difference.
  • after obtaining ⁇ , the range of the adjusted sample PE k ' in the combined segment is obtained as a branch of the second correction, and the corresponding corrected second result PE kf ' is selected according to the calculated branch.
  • the input of the structure diagram in Fig. 6 is the data after the adjustment in the segment, and the output is the data after the inter-segment adjustment.
  • the above method for adjusting the phase offset value adopts a combination of parallel segment adjustment and inter-stage adjustment two-step calculation, and the adjustment value is calculated between the partitions during the adjustment process, and has the characteristics of high adjustment precision.
  • S104 Obtain a corrected phase value of each sample in each average block according to a phase offset value of each average block and a phase value of each sample in each average block.
  • S104 may include: a phase offset value for each average block, a phase value of each sample point in each average block, and a - ⁇ /4 summation calculation to obtain each average The corrected phase value for each sample in the block.
  • the phase offset value theta2, the phase offset value theta2, and the signal angle value in_theta after the frequency offset correction are obtained, and in the compensation operation, each of the compensation operations
  • the phase value of the sample is compensated for the phase offset value of the average block where each sample is located, that is, in the actual application, the adjusted phase offset value plus the signal angle value after the frequency offset correction is further reduced.
  • phase offset correction To fix the phase ⁇ /4, complete the phase offset correction; complete the phase offset correction phase angle value (theta3) through the cosine cos and sine sin lookup table, get the corresponding cos(theta3) and sin(theta3)
  • the value of cos(theta3) is multiplied by the input amplitude value (in_mod) to obtain the output out_i in the Cartesian coordinate system.
  • the sin(theta3) is multiplied by the input amplitude value (in_mod) to obtain the output out_q in the Cartesian coordinate system.
  • the conversion from polar coordinates to Cartesian coordinates yields the final output.
  • the phase correction method provided by the embodiment of the present invention first performs M-th power processing on the signal values of the samples, and divides the signal values of the samples processed by the M-th power into at least one parallel segment, and then parallelizes each of the samples.
  • the signal values of the samples in the segment are divided to obtain a preset number of average blocks, and then the signal values of the samples in each average block are combined to obtain the signal value of each average block, thus avoiding the
  • the defect that the signal values are calculated is only required to perform low-pass filtering processing on the preset number of average blocks to determine the phase offset value of each average block, and finally, according to the phase offset value of each average block and each
  • the phase value of each sample in the average block is obtained as a corrected phase value for each sample in each average block, that is, in the embodiment of the present invention, by dividing the signal value of the sample points in each parallel segment After the combination, the phase offset value is determined for each average block, which greatly reduces the calculation amount in determining the phase offset value, thereby
  • FIG. 7 is a schematic structural diagram of a phase correction device according to an embodiment of the present invention. As shown in FIG. 7, the device includes: a division module 71, a processing module 72, Determining module 73 and correction module 74;
  • the dividing module 71 is configured to perform M-th power processing on the signal values of the samples, and divide the signal values of the samples processed by the M-th power into at least one parallel segment, where M is a preset value;
  • the module 72 is configured to divide the signal values of the samples in each parallel segment to obtain a preset number of average blocks, and combine the signal values of the samples in each average block to obtain the signal values of each average block, for each
  • the signal values of the average blocks are subjected to low-pass filtering processing to obtain low-pass filtered signal values for each average block; and the determining module 73 is configured to determine the average block for each low-pass filtered signal value of each average block.
  • a correction module 74 configured to obtain each sample in each average block based on a phase offset value of each average block and a phase value of each sample in each average block The corrected phase value.
  • the at least one parallel segment is obtained by dividing the signal values of the samples after the M-th power processing, in an optional embodiment.
  • the dividing module 71 is further configured to determine the size of the sliding window according to one or more of the following: the transmission distance of the signals of the samples, the laser line width of the signals of the samples, and the processing after the M-th power processing
  • the signal values of the samples are divided to obtain at least one parallel segment, wherein the number of samples of each parallel segment is the size of the determined sliding window.
  • the phase correction means after obtaining a predetermined number of average blocks, is processed for each average block.
  • the processing module 72 is further configured to The summation of the signal values of the samples in each average block is performed to obtain the signal value of each average block.
  • the processing module 72 is further configured to perform a summation calculation on the abscissa and the ordinate of the signal values of the samples in each average block to obtain the abscissa of the signal value of each average block and the signal value of each average block.
  • the ordinate of the signal; the abscissa of the signal value of each average block and the ordinate of the signal value of each average block are respectively subjected to low-pass filtering processing, and the abscissa and each of the low-pass filtered signal values of each average block are obtained.
  • the ordinate of the low-pass filtered signal value of the average block is a summation calculation on the abscissa and the ordinate of the signal values of the samples in each average block.
  • the obtained low-pass filtered signal value of each average block includes a phase offset value which is a value in a Cartesian coordinate system.
  • the determining module 73 is further configured to convert the low-pass filtered signal value of each average block into a low-pass filtered signal value of each average block in polar coordinates; each average block in polar coordinates The phase value corresponding to the low-pass filtered signal value is divided by M to obtain the phase offset value of each average block.
  • the device further includes: a first adjustment module (not shown) configured to determine a phase offset value of each corresponding average block after the low-pass filtered signal value according to each average block, The phase offset value of each average block and the phase value of each sample in each average block, before obtaining the corrected phase value of each sample in each average block, in each parallel segment, each will The difference between the phase offset value of the average block and the phase offset value of the corresponding previous average block is determined as a first difference; each average block is adjusted according to a preset value interval in which the first difference falls Phase offset value.
  • the above device further includes: a second adjustment module (not shown) configured to adjust the phase offset value of each average block after the preset value interval falling according to the first difference, in each parallel segment
  • the difference between the phase offset value of the first average block and the phase offset value of the last average block in the corresponding previous parallel segment is determined as the second difference; the preset according to the second difference
  • the value interval determines a correction value corresponding to the second difference; summing the correction value and the phase offset value of the first average block in each parallel segment to obtain a phase offset value of the adjusted first parallel block; Adjusting the phase offset value of the jth parallel block in each parallel segment according to the modified value and the adjusted phase offset value of the first parallel block, where j is greater than or equal to 2.
  • the correction module 74 in order to determine the corrected phase value, is further configured as a phase offset value for each average block, and a phase value of each sample in each average block. And the - ⁇ /4 summation calculation yields the corrected phase values for each sample in each average block.
  • the dividing module 71, the processing module 72, the determining module 73, and the correcting module 74 may all be provided by a central processing unit (CPU), a microprocessor (MPU, a microprocessor unit), an application specific integrated circuit ( ASIC, Application Specific Integrated Circuit) or Field-Programmable Gate Array (FPGA).
  • CPU central processing unit
  • MPU microprocessor
  • ASIC Application Specific Integrated Circuit
  • FPGA Field-Programmable Gate Array
  • This embodiment describes a computer storage medium, which may be a ROM (for example, a read only memory, a FLASH memory, a transfer device, etc.), a magnetic storage medium (for example, a magnetic tape, a magnetic disk drive, etc.), an optical storage medium (for example, a CD-ROM). , a DVD-ROM, a paper card, a paper tape, etc., and other well-known types of program memory; computer-readable medium storing computer-executable instructions that, when executed, cause at least one processor to perform operations comprising:
  • Performing M-th power processing on the signal values of the respective points dividing the signal values of the samples processed by the M-th power into at least one parallel segment; dividing the signal values of the samples in each parallel segment to obtain a preset number
  • the average block, the signal values of the samples in each average block are combined to obtain the signal value of each average block, and the signal values of each average block are subjected to low-pass filtering processing to obtain low-pass filtered each average block.
  • a signal value; determining a phase offset value of each average block according to a low-pass filtered signal value of each average block; a phase offset value according to each average block and a phase value of each sample point in each average block Get the corrected phase value for each sample in each average block.
  • the phase correction method provided by the embodiment of the present invention first performs M-th power processing on the signal values of the samples, and divides the signal values of the samples processed by the M-th power into at least one parallel segment, and then parallelizes each of the samples.
  • the signal values of the samples in the segment are divided to obtain a preset number of average blocks, and then the signal values of the samples in each average block are combined to obtain the signal value of each average block, thus avoiding the
  • the defect that the signal values are calculated is only required to perform low-pass filtering processing on the preset number of average blocks to determine the phase offset value of each average block, and finally, according to the phase offset value of each average block and each
  • the phase value of each sample in the average block is obtained as a corrected phase value for each sample in each average block, that is, in the embodiment of the present invention, by dividing the signal value of the sample points in each parallel segment After the combination, the phase offset value is determined for each average block, which greatly reduces the calculation amount in determining the phase offset value, thereby
  • the disclosed apparatus and method may be implemented in other manners.
  • the device embodiments described above are merely illustrative.
  • the division of the unit is only a logical function division.
  • there may be another division manner such as: multiple units or components may be combined, or Can be integrated into another system, or some features can be ignored or not executed.
  • the coupling, or direct coupling, or communication connection of the various components shown or discussed may be through some interface, device or unit.
  • the indirect coupling or communication connection can be electrical, mechanical or other form.
  • the units described above as separate components may or may not be physically separated, and the components displayed as the unit may or may not be physical units; they may be located in one place or distributed on multiple network units; Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of the embodiment.
  • each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may be separately used as one unit, or two or more units may be integrated into one unit;
  • the unit can be implemented in the form of hardware or in the form of hardware plus software functional units.
  • the foregoing program may be stored in a computer readable storage medium, and when executed, the program includes The foregoing steps of the method embodiment; and the foregoing storage medium includes: a removable storage device, a read only memory (ROM), a magnetic disk or an optical disk, and the like, which can store program codes.
  • ROM read only memory
  • the above-described integrated unit of the present invention may be stored in a computer readable storage medium if it is implemented in the form of a software function module and sold or used as a standalone product.
  • the technical solution of the embodiments of the present invention may be embodied in the form of a software product in essence or in the form of a software product stored in a storage medium, including a plurality of instructions.
  • a computer device (which may be a personal computer, server, or network device, etc.) is caused to perform all or part of the methods described in various embodiments of the present invention.
  • the foregoing storage medium includes various media that can store program codes, such as a mobile storage device, a ROM, a magnetic disk, or an optical disk.
  • the signal value of each sample is processed by M-th power, and the signal value of each sample processed by the M-th power is divided into at least one parallel segment, and then the sample is sampled for each parallel segment.
  • the signal value of the point is divided to obtain a preset number of average blocks, and then the signal values of the samples in each average block are combined to obtain the signal value of each average block, thus avoiding the signal values for each sample point.
  • the defect of calculation is that only a preset number of average blocks need to be low-pass filtered to determine the phase offset value of each average block, and finally, according to the phase offset value of each average block and each average block The phase value of each sample is obtained as the corrected phase value of each sample in each average block.
  • the signal values of the samples in each parallel segment are divided and combined, and then averaged for each average.
  • the block is used to determine the phase offset value, which greatly reduces the calculation amount in determining the phase offset value, thereby reducing the amount of calculation in the phase correction process, and further improving the processing efficiency of the phase correction.

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Nonlinear Science (AREA)
  • Optical Communication System (AREA)
  • Image Analysis (AREA)

Abstract

Disclosed is a phase correction method, comprising: performing Mth-power processing on signal values of sample points, and dividing the signal values of the sample points having undergone Mth-power processing into at least one parallel section; dividing signal values of sample points in each parallel section into a preset number of average blocks, combining the signal values of the sample points in each average block to obtain the signal value of each average block, and low-pass filtering the signal value of each average block to obtain the low-pass filtered signal value of each average block; determining the phase offset value of each average block according to the low-pass filtered signal value of each average block; and obtaining a corrected phase value for each sample point in each average block according to the phase offset value of each average block and the phase value of each sample point in each average block. Also disclosed are a phase correction device and a computer storage medium.

Description

一种相位校正方法和装置、计算机存储介质Phase correction method and device, computer storage medium
相关申请的交叉引用Cross-reference to related applications
本申请基于申请号为201610841514.9、申请日为2016年09月22日的中国专利申请提出,并要求该中国专利申请的优先权,该中国专利申请的全部内容在此引入本申请作为参考。The present application is filed on the basis of the Chinese Patent Application No. PCT Application No.
技术领域Technical field
本发明涉及光通信领域,尤其涉及一种相位校正方法和装置、计算机存储介质。The present invention relates to the field of optical communications, and in particular, to a phase correction method and apparatus, and a computer storage medium.
背景技术Background technique
在传输技术的发展中,光纤是一种不可或缺的媒介,如何用最少量的光纤传输最丰富的信息,出于这种探索,光传输的发展基本经历了以下几个阶段:空分复用(SDM,Space Division Multiplexing)阶段、时分复用(TDM,Time Division)阶段和波分复用(WDM,Wavelength Division Multiplexing)阶段。In the development of transmission technology, optical fiber is an indispensable medium. How to transmit the most abundant information with the least amount of optical fiber. For this kind of exploration, the development of optical transmission has basically gone through the following stages: air separation The (SDM, Space Division Multiplexing) phase, the Time Division Multiplexing (TDM) phase, and the Wavelength Division Multiplexing (WDM) phase.
目前,有线传输以波分***为主,随着通信技术的发展,商用的40G波分传输逐渐演变到100G、400G传输,与此同时,在数据传输距离上也在不断的拓展,这样,波分***在传输过程中会带来色度色散、偏振膜色散、强滤波效应、频偏相偏等诸多问题需要解决。At present, the cable transmission is dominated by the WDM system. With the development of communication technology, the commercial 40G WDM transmission has gradually evolved to 100G and 400G transmission. At the same time, the data transmission distance is also constantly expanding. The sub-system will bring about many problems such as chromatic dispersion, polarization film dispersion, strong filtering effect, frequency offset and other deviations in the transmission process.
其中,相偏的存在使得信号的星座图在弧度方向发生延展,如果过大,会导致信号点中间发生交叠,从而产生误码;目前主流的相偏相位校正方法是四次方法,但业界的大多数做法是采用平面直角坐标系下的数据进行计算,例如通过直角坐标系下的数据求出样点的幅度值,依据幅度值对样 点分类,根据分类结果对样点在直角坐标系下进行四次方相偏估计,再根据相偏估计值计算相偏角度;虽然这一计算方法被广泛采用,但是该方法对每个样点均进行计算,进而导致数据运算量大,采用的电路结构复杂、固定,实现起来电路面积大、功耗大,进行相位校正时处理效率较低。Among them, the existence of phase deviation makes the constellation of the signal extend in the direction of the arc. If it is too large, it will cause the signal points to overlap in the middle, resulting in bit errors. Currently, the mainstream phase-off phase correction method is four times, but the industry Most of the methods are calculated by using the data in the plane rectangular coordinate system. For example, the amplitude value of the sample is obtained by the data in the Cartesian coordinate system, and the sample is determined according to the amplitude value. According to the classification result, the sample is subjected to the quadratic phase deviation estimation in the Cartesian coordinate system according to the classification result, and then the phase deviation angle is calculated according to the phase deviation estimation value; although this calculation method is widely used, the method is applicable to each sample point. All calculations are carried out, which leads to a large amount of data calculation. The circuit structure used is complicated and fixed, and the circuit area is large and the power consumption is large. When the phase correction is performed, the processing efficiency is low.
发明内容Summary of the invention
有鉴于此,本发明实施例期望提供一种相位校正方法和装置、计算机存储介质,减小了在进行相位校正过程中的运算量,提高了相位校正的处理效率。In view of this, embodiments of the present invention are expected to provide a phase correction method and apparatus, and a computer storage medium, which reduce the amount of calculation in performing phase correction and improve the processing efficiency of phase correction.
第一方面,本发明实施例提供了一种相位校正方法,包括:对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分至至少一个并行段,其中,M为预设值;对每个并行段中样点的信号值进行划分得到预设数目的平均块,对每个平均块中样点的信号值进行合并得到所述每个平均块的信号值,对所述每个平均块的信号值进行低通滤波处理,得到所述每个平均块低通滤波后的信号值;根据所述每个平均块低通滤波后的信号值,确定所述每个平均块的相位偏移值;根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相位值。In a first aspect, an embodiment of the present invention provides a phase correction method, including: performing M-th power processing on signal values of respective samples, and dividing signal values of each sample processed by the M-th power into at least one parallel segment. Wherein, M is a preset value; dividing a signal value of the sample point in each parallel segment to obtain a preset number of average blocks, and combining signal values of the sample points in each average block to obtain each of the average blocks a signal value, performing low-pass filtering processing on the signal value of each of the average blocks to obtain a low-pass filtered signal value of each of the average blocks; and according to the low-pass filtered signal value of each of the average blocks, Determining a phase offset value of each of the average blocks; obtaining, according to a phase offset value of each of the average blocks and a phase value of each of the samples in each of the average blocks, The corrected phase value of each sample.
在上述技术方案中,所述将M次方处理后的各样点的信号值划分至至少一个并行段,包括:根据以下一项或多项确定滑动窗口的大小:各样点的信号的传输距离、各样点的信号的激光器线宽;对所述M次方处理后的各样点的信号值进行划分,得到所述至少一个并行段,其中,每个并行段的样点的个数为所确定的滑动窗口的大小。In the above technical solution, the dividing the signal value of each sample processed by the Mth power into at least one parallel segment comprises: determining the size of the sliding window according to one or more of the following: transmission of signals of each sample point a laser line width of the signal of each sample; dividing the signal values of the samples processed by the M-th power to obtain the at least one parallel segment, wherein the number of samples of each parallel segment The size of the sliding window determined.
在上述技术方案中,所述对每个平均块中样点的信号值进行合并得到所述每个平均块的信号值,包括:对每个平均块中样点的信号值进行求和计算,得到所述每个平均块的信号值。 In the above technical solution, the combining the signal values of the samples in each average block to obtain the signal values of each of the average blocks includes: summing the signal values of the samples in each average block, The signal values of each of the average blocks are obtained.
在上述技术方案中,所述对每个平均块中样点的信号值进行求和计算,得到所述每个平均块的信号值,包括:对所述每个平均块中样点的信号值的横坐标和纵坐标分别进行求和计算,得到所述每个平均块的信号值的横坐标和所述每个平均块的信号值的纵坐标;相应地,所述对所述每个平均块的信号值进行低通滤波处理,得到所述每个平均块低通滤波后的信号值,包括:分别对所述每个平均块的信号值的横坐标和所述每个平均块的信号值的纵坐标进行低通滤波处理,得到所述每个平均块低通滤波后的信号值的横坐标和所述每个平均块低通滤波后的信号值的纵坐标。In the above technical solution, the summation of the signal values of the samples in each average block is performed to obtain the signal value of each of the average blocks, including: the signal value of the sample points in each of the average blocks The abscissa and the ordinate are respectively subjected to a summation calculation to obtain an abscissa of the signal value of each of the average blocks and an ordinate of the signal value of each of the average blocks; correspondingly, the pair of the averages The signal value of the block is subjected to low-pass filtering processing to obtain a low-pass filtered signal value of each of the average blocks, including: an abscissa of the signal value of each of the average blocks and a signal of each of the average blocks The ordinate of the value is subjected to low-pass filtering processing to obtain the abscissa of the low-pass filtered signal value of each of the average blocks and the ordinate of the low-pass filtered signal value of each of the average blocks.
在上述技术方案中,所述根据所述每个平均块低通滤波后的信号值,确定对应的所述每个平均块的相位偏移值,包括:将所述每个平均块低通滤波后的信号值,转换成极坐标下的每个平均块低通滤波后的信号值;将所述极坐标下的每个平均块低通滤波后的信号值对应的相位值除以M,得出所述每个平均块的相位偏移值。In the above technical solution, determining, according to the low-pass filtered signal value of each average block, determining a phase offset value of each of the corresponding average blocks, comprising: low-pass filtering each of the average blocks The subsequent signal value is converted into a low-pass filtered signal value of each average block in polar coordinates; the phase value corresponding to the low-pass filtered signal value of each average block in the polar coordinates is divided by M, The phase offset value of each of the average blocks is derived.
在上述技术方案中,在根据所述每个平均块低通滤波后的信号值,确定对应的所述每个平均块的相位偏移值之后,在根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相位值之前,所述方法还包括:在所述每个并行段中,将所述每个平均块的相位偏移值与对应的上一平均块的相位偏移值之间的差值确定为第一差值;根据所述第一差值所落入的预设数值区间,调整所述每个平均块的相位偏移值。In the above technical solution, after determining the phase offset value of the corresponding each average block according to the signal value of the low-pass filtering of each of the average blocks, after the phase offset according to each of the average blocks Before the value and the phase value of each of the samples in each of the average blocks, the corrected phase values of each of the samples in each of the average blocks are obtained, the method further comprising: in each of the parallel segments Determining, as a first difference, a difference between a phase offset value of each of the average blocks and a phase offset value of the corresponding previous average block; a preset according to the first difference A numerical interval that adjusts a phase offset value of each of the average blocks.
在上述技术方案中,在根据所述第一差值所落入的预设数值区间,调整所述每个平均块的相位偏移值之后,所述方法还包括:将所述每个并行段中第一个平均块的相位偏移值与对应的上一并行段中最后一个平均块的相位偏移值之间的差值确定为第二差值;根据所述第二差值所落入的预设数值区间确定所述第二差值对应的修正值;对所述修正值与所述每个并行 段中第一个平均块的相位偏移值进行求和,得到调整后的第一个并行块的相位偏移值;根据所述修正值和所述调整后的第一个并行块的相位偏移值,调整所述每个并行段中第j个并行块的相位偏移值,其中,j大于等于2。In the above technical solution, after adjusting the phase offset value of each average block according to a preset value interval that falls within the first difference, the method further includes: each parallel segment a difference between a phase offset value of the first average block and a phase offset value of a last average block in the corresponding previous parallel segment is determined as a second difference; falling according to the second difference The preset value interval determines a correction value corresponding to the second difference; and the correction value is parallel to each of the Phase offset values of the first average block in the segment are summed to obtain a phase offset value of the adjusted first parallel block; and according to the modified value and the phase offset of the adjusted first parallel block And shifting, adjusting a phase offset value of the jth parallel block in each parallel segment, where j is greater than or equal to 2.
在上述技术方案中,所述根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相位值,包括:对所述每个平均块的相位偏移值、所述每个平均块中每个样点的相位值以及-π/4求和计算,得到所述每个平均块中每个样点的校正后相位值。In the above technical solution, the obtaining, according to the phase offset value of each of the average blocks and the phase value of each sample in each of the average blocks, obtaining a correction for each sample in each of the average blocks a post-phase value, comprising: a phase offset value for each of the average blocks, a phase value of each of the samples in each of the average blocks, and a -π/4 summation calculation to obtain each of the average blocks The corrected phase value for each sample.
第二方面,本发明实施例提供了一种相位校正装置,包括:划分模块,配置为对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分至至少一个并行段,其中,M为预设值;处理模块,配置为对每个并行段中样点的信号值进行划分得到预设数目的平均块,对每个平均块中样点的信号值进行合并得到所述每个平均块的信号值,对所述每个平均块的信号值进行低通滤波处理,得到所述每个平均块低通滤波后的信号值;确定模块,配置为根据所述每个平均块低通滤波后的信号值,确定所述每个平均块的相位偏移值;校正模块,配置为根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相位值。In a second aspect, an embodiment of the present invention provides a phase correction apparatus, including: a dividing module configured to perform a power processing on a signal value of each sample point, and divide a signal value of each sample processed by the M-th power processing To at least one parallel segment, wherein M is a preset value; the processing module is configured to divide the signal values of the samples in each parallel segment to obtain a preset number of average blocks, and the signals of the samples in each average block The values are combined to obtain the signal values of the average blocks, and the signal values of the average blocks are subjected to low-pass filtering processing to obtain low-pass filtered signal values of each of the average blocks; Determining, according to the low-pass filtered signal value of each of the average blocks, a phase offset value of each of the average blocks; and a correction module configured to: according to the phase offset value of each of the average blocks and each of the The phase values of each sample in the average block yield the corrected phase values for each of the samples in each of the average blocks.
在上述技术方案中,所述划分模块,还配置为根据以下一项或多项确定滑动窗口的大小:各样点的信号的传输距离、各样点的信号的激光器线宽;对所述M次方处理后的各样点的信号值进行划分,得到所述至少一个并行段,其中,每个并行段的样点的个数为所确定的滑动窗口的大小。In the above technical solution, the dividing module is further configured to determine a size of the sliding window according to one or more of: a transmission distance of a signal of each sample point, a laser line width of a signal of each sample point; and the M The signal values of the samples processed by the power are divided to obtain the at least one parallel segment, wherein the number of samples of each parallel segment is the size of the determined sliding window.
在上述技术方案中,所述处理模块,还配置为对每个平均块中样点的信号值进行求和计算,得到所述每个平均块的信号值。In the above technical solution, the processing module is further configured to perform a summation calculation on the signal values of the samples in each average block to obtain a signal value of each of the average blocks.
在上述技术方案中,所述处理模块,还配置为对所述每个平均块中样 点的信号值的横坐标和纵坐标分别进行求和计算,得到所述每个平均块的信号值的横坐标和所述每个平均块的信号值的纵坐标;分别对所述每个平均块的信号值的横坐标和所述每个平均块的信号值的纵坐标进行低通滤波处理,得到所述每个平均块低通滤波后的信号值的横坐标和所述每个平均块低通滤波后的信号值的纵坐标。In the above technical solution, the processing module is further configured to sample the average block Calculating the abscissa and the ordinate of the signal value of the point respectively, obtaining the abscissa of the signal value of each of the average blocks and the ordinate of the signal value of each of the average blocks; respectively for each of the averages Performing low-pass filtering processing on the abscissa of the signal value of the block and the ordinate of the signal value of each of the average blocks, obtaining the abscissa of the low-pass filtered signal value of each of the average blocks and each of the average blocks The ordinate of the low pass filtered signal value.
在上述技术方案中,所述确定模块,还配置为将所述每个平均块低通滤波后的信号值,转换成极坐标下的每个平均块低通滤波后的信号值;将所述极坐标下的每个平均块低通滤波后的信号值对应的相位值除以M,得出所述每个平均块的相位偏移值。In the above technical solution, the determining module is further configured to convert the low-pass filtered signal value of each average block into a low-pass filtered signal value of each average block in polar coordinates; The phase value corresponding to the low-pass filtered signal value of each average block in polar coordinates is divided by M to obtain the phase offset value of each of the average blocks.
在上述技术方案中,所述装置还包括:第一调整模块,配置为在根据所述每个平均块低通滤波后的信号值,确定对应的所述每个平均块的相位偏移值之后,在根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相位值之前,在所述每个并行段中,将所述每个平均块的相位偏移值与对应的上一平均块的相位偏移值之间的差值确定为第一差值;根据所述第一差值所落入的预设数值区间,调整所述每个平均块的相位偏移值。In the above technical solution, the apparatus further includes: a first adjustment module, configured to determine a phase offset value of each of the corresponding average blocks after the low-pass filtered signal value according to each of the average blocks Before, according to the phase offset value of each of the average blocks and the phase value of each sample in each of the average blocks, before obtaining the corrected phase values of each of the samples in each of the average blocks, In each parallel segment, determining a difference between a phase offset value of each of the average blocks and a phase offset value of the corresponding previous average block as a first difference; according to the first difference The preset value interval in which the value falls, adjusting the phase offset value of each of the average blocks.
在上述技术方案中,所述装置还包括:第二调整模块,配置为在根据所述第一差值所落入的预设数值区间,调整所述每个平均块的相位偏移值之后,将所述每个并行段中第一个平均块的相位偏移值与对应的上一并行段中最后一个平均块的相位偏移值之间的差值确定为第二差值;根据所述第二差值所落入的预设数值区间确定所述第二差值对应的修正值;对所述修正值与所述每个并行段中第一个平均块的相位偏移值进行求和,得到调整后的第一个并行块的相位偏移值;根据所述修正值和所述调整后的第一个并行块的相位偏移值,调整所述每个并行段中第j个并行块的相位偏移值,其中,j大于等于2。 In the above technical solution, the device further includes: a second adjustment module, configured to adjust a phase offset value of each of the average blocks after a preset value interval that falls according to the first difference, Determining, by the difference between the phase offset value of the first average block in each parallel segment and the phase offset value of the last average block in the corresponding previous parallel segment as a second difference; a preset value interval in which the second difference falls to determine a correction value corresponding to the second difference; summing the correction value and a phase offset value of the first average block in each parallel segment Obtaining a phase offset value of the adjusted first parallel block; adjusting the jth parallel in each parallel segment according to the modified value and the adjusted phase offset value of the first parallel block The phase offset value of the block, where j is greater than or equal to two.
在上述技术方案中,所述校正模块,还配置为对所述每个平均块的相位偏移值、所述每个平均块中每个样点的相位值以及-π/4求和计算,得到所述每个平均块中每个样点的校正后相位值。In the above technical solution, the correction module is further configured to calculate a phase offset value of each of the average blocks, a phase value of each sample in each of the average blocks, and a -π/4 summation calculation, A corrected phase value for each of the samples in each of the average blocks is obtained.
第三方面,本发明实施例提供了一种计算机存储介质,该计算机存储介质存储有计算机程序,该计算机程序配置为执行上述相位校正方法。In a third aspect, an embodiment of the present invention provides a computer storage medium storing a computer program configured to perform the phase correction method.
本发明实施例所提供的相位校正方法和装置、计算机存储介质,首先对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分成至少一个并行段,然后对每个并行段中样点的信号值进行划分得到预设数目的平均块,再对每个平均块中样点的信号值进行合并得到每个平均块的信号值,这样,避免了对每个样点的信号值都进行计算的缺陷,仅仅只需要对预设数目的平均块进行低通滤波处理来确定每个平均块的相位偏移值,最后,根据每个平均块的相位偏移值和每个平均块中每个样点的相位值,得到每个平均块中每个样点的校正后相位值,也就是说,本发明实施例,通过对每个并行段中样点的信号值划分和合并之后,再针对每个平均块来进行相位偏移值的确定,大大减少了在确定相位偏移值中的计算量,进而也减少了在进行相位校正过程中的运算量,进一步地提高了相位校正的处理效率。The phase correction method and device and the computer storage medium provided by the embodiments of the present invention first perform M-th power processing on the signal values of the samples, and divide the signal values of the samples processed by the M-th power into at least one parallel segment. Then, the signal values of the samples in each parallel segment are divided to obtain a preset number of average blocks, and then the signal values of the samples in each average block are combined to obtain the signal value of each average block, thus avoiding For the defect of calculating the signal value of each sample, only a preset number of average blocks need to be low-pass filtered to determine the phase offset value of each average block, and finally, according to the phase of each average block The offset value and the phase value of each sample in each average block yield a corrected phase value for each sample in each average block, that is, in the embodiment of the present invention, by sampling each parallel segment After the signal values of the points are divided and combined, the phase offset value is determined for each average block, which greatly reduces the calculation amount in determining the phase offset value, thereby reducing the phase correction process. Amount of computation, the processing efficiency is further improved phase correction.
附图说明DRAWINGS
图1为本发明实施例中相位校正方法的流程示意图;1 is a schematic flow chart of a phase correction method according to an embodiment of the present invention;
图2为本发明实施例中相位校正方法的一种可选的流程框图;2 is an optional flow chart of a phase correction method according to an embodiment of the present invention;
图3为本发明实施例中相位校正方法的另一种可选的流程框图;3 is another optional flow chart of a phase correction method according to an embodiment of the present invention;
图4为本发明实施例中一种可选的加法树形结构的示意图;4 is a schematic diagram of an optional addition tree structure according to an embodiment of the present invention;
图5为本发明实施例中对相位偏移值进行调整的一种可选的流程框图;FIG. 5 is an optional flow chart of adjusting a phase offset value according to an embodiment of the present invention; FIG.
图6为本发明实施例中对相位偏移值进行调整的另一种可选的流程框图; 6 is another optional flow chart of adjusting a phase offset value according to an embodiment of the present invention;
图7为本发明实施例中相位校正装置的结构示意图。FIG. 7 is a schematic structural diagram of a phase correcting apparatus according to an embodiment of the present invention.
具体实施方式detailed description
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。The technical solutions in the embodiments of the present invention will be clearly and completely described in the following with reference to the accompanying drawings.
本发明实施例提供一种相位校正方法,该方法可以应用于光传输网络中对相偏相位的校正,还可以应用于其他传输网络中对相位的校正,这里本发明不做具体限定;The embodiment of the present invention provides a phase correction method, which can be applied to the correction of phase offset phase in an optical transmission network, and can also be applied to phase correction in other transmission networks, and the present invention is not specifically limited herein;
图1为本发明实施例中相位校正方法的流程示意图,如图1所示,该方法包括:1 is a schematic flowchart of a phase correction method according to an embodiment of the present invention. As shown in FIG. 1, the method includes:
S101:对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分至至少一个并行段;S101: Perform M-th power processing on the signal values of the samples, and divide the signal values of the samples processed by the M-th power into at least one parallel segment;
其中,M为预设值,M的取值可以为4、8和12等,这里,本发明不做具体限定,另外,上述并行段包括有固定数目的经过M次方处理后的各样点的信号值,该固定数目等于下述滑动窗口的大小;Wherein, M is a preset value, and the value of M may be 4, 8, and 12, etc., and the present invention is not specifically limited. In addition, the parallel segment includes a fixed number of samples processed by the M-th power. Signal value, the fixed number is equal to the size of the sliding window described below;
需要说明的是,仅仅对各样点的信号值的相位进行M次方处理,若以M=4为例,图2为本发明实施例中相位校正方法的一种可选的流程框图,如图2所示,接收信号的相位,可以表示如下:It should be noted that only the phase of the signal value of each sample is processed by M-th power. If M=4 is taken as an example, FIG. 2 is an optional flow chart of the phase correction method in the embodiment of the present invention, such as As shown in Figure 2, the phase of the received signal can be expressed as follows:
θk=θs(k)+ΔωkTinASE                 (1)θ ks (k)+ΔωkT inASE (1)
其中,θs(k)为第k个样点的信号相位,ΔwkTi为频偏相位,θn为由激光器线宽引起的相偏相位,θASE为噪声相位;其中,假设在进行四次方处理之前进行频偏估计,那么ΔwkTi经过频偏估计已经被去除,则剩下的相位进行四次方得到如下公式:Where θ s (k) is the signal phase of the kth sample, ΔwkT i is the frequency offset phase, θ n is the phase deviation phase caused by the laser line width, and θ ASE is the noise phase; wherein, it is assumed that it is performed four times. The frequency offset estimation is performed before the square processing, then the ΔwkT i has been removed by the frequency offset estimation, and the remaining phase is quadratic to obtain the following formula:
V4(k)=exp{j4θs(k)}·exp{j4θn}·exp{j4θASE}            (2)V 4 (k)=exp{j4θ s (k)}·exp{j4θ n }·exp{j4θ ASE } (2)
由上述公式(2)可以看出,经过四次方处理后的相位值V4(k)中还包括 θs(k)、θn和θASE;其中,由于θs(k)的取值只有如下几个固定值:0,+π/2,-π/2和π,所以,V4(k)可以去掉θs(k),那么,四次方处理后的相位值V4(k)中还包括θn和θASEIt can be seen from the above formula (2) that the phase value V 4 (k) after the fourth power processing further includes θ s (k), θ n and θ ASE ; wherein, since θ s (k) There are only a few fixed values: 0, +π/2, -π/2, and π, so V 4 (k) can remove θ s (k), then the phase value V 4 (k) after the fourth power processing ) also includes θ n and θ ASE .
图3为本发明实施例中相位校正方法的另一种可选的流程框图,如图3所示,在对各样点的信号值进行M次方处理之后,有三路值,分别为直角坐标系下的(in_i,in_q),信号角度值in_theta和信号幅度值in_mod;其中,对(in_i,in_q)进行低通滤波处理,对信号角度值in_theta和信号幅度值in_mod进行相位偏移值的补偿处理。FIG. 3 is another optional flow chart of the phase correction method according to the embodiment of the present invention. As shown in FIG. 3, after performing M-th power processing on the signal values of the samples, there are three values, which are rectangular coordinates. The system (in_i, in_q), the signal angle value in_theta and the signal amplitude value in_mod; wherein, (in_i, in_q) is subjected to low-pass filtering processing, and the signal angle value in_theta and the signal amplitude value in_mod are compensated for the phase offset value. deal with.
在对各样点的信号值的相位进行M次方处理之后,为了对M次方处理后的各样点的信号值进行划分得到至少一个并行段,在一种可选的实施例中,S101可以包括:After the M-th power processing of the phase of the signal values of the samples, the signal values of the samples after the M-th power processing are divided to obtain at least one parallel segment. In an optional embodiment, S101 Can include:
根据以下一项或多项确定滑动窗口的大小:各样点的信号的传输距离、各样点的信号的激光器线宽;对M次方处理后的各样点的信号值进行划分,得到至少一个并行段,其中,每个并行段的样点的个数为所确定的滑动窗口的大小。The size of the sliding window is determined according to one or more of the following: the transmission distance of the signals of each sample point, and the laser line width of the signal of each sample point; the signal values of the samples after the M-th power processing are divided to obtain at least A parallel segment in which the number of samples of each parallel segment is the size of the determined sliding window.
其中,由于应用于不同的场景,上述各样点的信号的传输距离是指信号从发射端到接收端之间的距离,上述各样点的信号的激光器线宽为产生相偏相位的激光器线宽;Wherein, since applied to different scenes, the transmission distance of the signals of the above samples refers to the distance between the signals from the transmitting end to the receiving end, and the laser line width of the signals of the above various points is the laser line that generates the phase-biased phase. width;
举例来说,传输距离越长噪声越大,激光器线宽越宽相偏相位变化的越快,滑动窗口的大小越大对噪声的抑制作用越好,但是滑动窗口的大小越大,连续若干个符号的相偏相位视为相同的前提遭到破坏,该组第一个样点和最后一个样点的相位偏差实际上相差较大,从而影响相位估计的准确性,所以需要根据场景不同配置不同的滑动窗口的大小。For example, the longer the transmission distance, the larger the noise, the wider the line width of the laser, the faster the phase deviation changes. The larger the size of the sliding window, the better the suppression of noise, but the larger the size of the sliding window, the number of consecutive The phase deviation phase of the symbol is regarded as the same premise is destroyed. The phase deviation of the first sample and the last sample in the group is actually quite different, which affects the accuracy of the phase estimation, so it needs to be configured differently according to the scene. The size of the sliding window.
具体自适应调整滑动窗口的大小的方法为:***上电后,利用该滑动窗口的大小的默认值(比如滑动窗口的大小的默认值为4)进行相偏估计和 补偿,并且统计此滑动窗口的大小下的误码率并在软件中记录下来,然后通过软件改变滑动窗口的大小,比如将滑动窗口的大小修改为12,统计12下的误码率,如果12下的误码率比4小,则继续增大滑动窗口的大小为36,如果36下的比12小,则继续增大滑动窗口的大小为84,如果84的误码率比36大,则选择36为最终的滑动窗口的大小;如果84下的误码率比36小,则继续增大滑动窗口的大小,直到硬件支持的最大值,如果最大值下的误码率最小,则选择最大的滑动窗口的大小,上述调整滑动窗口的大小的原则就是通过改变滑动窗口的大小,选择误码率最小的滑动窗口的大小。The method for adaptively adjusting the size of the sliding window is: after the system is powered on, the default value of the size of the sliding window (for example, the default value of the size of the sliding window is 4) is used for phase offset estimation and Compensate, and count the bit error rate under the size of this sliding window and record it in the software, then change the size of the sliding window through software, such as changing the size of the sliding window to 12, and counting the bit error rate under 12, if 12 If the error rate is lower than 4, the size of the sliding window is increased to 36. If the ratio of 36 is smaller than 12, the size of the sliding window is increased to 84. If the error rate of 84 is larger than 36, then Select 36 as the size of the final sliding window; if the bit error rate under 84 is smaller than 36, continue to increase the size of the sliding window until the maximum supported by the hardware, if the error rate at the maximum is the smallest, select the maximum The size of the sliding window, the above principle of adjusting the size of the sliding window is to select the size of the sliding window with the smallest bit error rate by changing the size of the sliding window.
上述滑动窗口的大小可以配置为84、36、12和4这4种不同值,滑动窗口的大小都为4的整数倍,根据滑动窗口的大小来对各样点的信号值进行划分,得到至少一个并行段,每个并行段的样点的个数等于滑动窗口的大小,针对每个并行段进行处理,而不是对每个样点进行处理,并且根据不同的应用场景灵活配置滑动窗口的大小,那么依场景不同自适应选择不同滑动窗口的大小,从而可以达到相位校正性能最优,并根据滑动窗口的大小采用了不同的加法树形结构,降低了运算量,减少了资源消耗,简化了设计。The size of the sliding window can be configured as four different values of 84, 36, 12, and 4. The size of the sliding window is an integer multiple of 4, and the signal values of each sample are divided according to the size of the sliding window to obtain at least A parallel segment, the number of samples of each parallel segment is equal to the size of the sliding window, and is processed for each parallel segment instead of processing each sample, and the size of the sliding window is flexibly configured according to different application scenarios. Then, according to different scenarios, the size of different sliding windows is adaptively selected, so that the phase correction performance is optimal, and different addition tree structures are adopted according to the size of the sliding window, which reduces the amount of calculation, reduces resource consumption, and simplifies. design.
S102:对每个并行段中样点的信号值进行划分得到预设数目的平均块,对每个平均块中样点的信号值进行合并得到每个平均块的信号值,对每个平均块的信号值进行低通滤波处理,得到每个平均块低通滤波后的信号值;S102: Dividing the signal values of the samples in each parallel segment to obtain a preset number of average blocks, and combining the signal values of the samples in each average block to obtain signal values of each average block, for each average block. The signal value is subjected to low-pass filtering processing to obtain a low-pass filtered signal value of each average block;
仍然以图3为例进行说明,对直角坐标系下的(in_i,in_q)这两个偏振态相互独立地进行低通滤波处理,两者之间互不影响,得到低通滤波处理之后的结果为(i1,q1);Still taking FIG. 3 as an example, the two polarization states (in_i, in_q) in the Cartesian coordinate system are independently subjected to low-pass filtering processing, and the two do not affect each other, and the result after low-pass filtering processing is obtained. Is (i1,q1);
具体来说,在对各样点进行划分得到至少一个并行段之后,针对每个并行段,再对每个并行段中样点的信号值进行划分得到预设数目的平均块,这里每个平均块的样点的个数为窗口滑动位移,窗口滑动位移可以配置为 4、8、12等,其中,这里预设数目的取值为滑动窗口的大小除以窗口滑动位移;Specifically, after dividing each sample into at least one parallel segment, for each parallel segment, dividing the signal values of the samples in each parallel segment to obtain a preset number of average blocks, where each average The number of samples of the block is the sliding displacement of the window, and the sliding displacement of the window can be configured as 4, 8, 12, etc., wherein the preset number of values is the size of the sliding window divided by the sliding displacement of the window;
在得到预设数目的平均块之后,针对每个平均块进行处理,为了得到每个平均块的信号值,在一种可选的实施例中,S102中对每个平均块中样点的信号值进行合并得到每个平均块的信号值,可以包括:对每个平均块中样点的信号值进行求和计算,得到每个平均块的信号值。After a predetermined number of average blocks are obtained, processing is performed for each average block. To obtain the signal value of each average block, in an alternative embodiment, the signal for each of the average block samples in S102 The values are combined to obtain the signal value of each average block, which may include: summing the signal values of the samples in each average block to obtain the signal value of each average block.
上述在对每个平均块中样点的信号值进行求和计算中采用了一种加法树形结构,图4为本发明实施例中一种可选的加法树形结构的示意图,如图4所示,为12的加法分支复用窗口滑动位移为4的结构,对12个样点的信号值划分成3个平均块,每4个样点划分为一个平均块,对第0-3个样点的信号值进行求和计算得到sum4_0,对第4-7个样点的信号值进行求和计算得到sum4_1,对第8-11个样点的信号值进行求和计算得到sum4_2;这一方法是对于不同加法计算路径上采取复用,根据应用场景不同,滑动窗口的大小可配置为K1、K2、K3、K4等,其中K1、K2、K3、K4等为正整数,例如,滑动窗口的大小可配置为84、36、12和4这4种不同值,滑动窗口的大小都为4的整数倍,当滑动窗口的大小配置为12的加法分支时,复用窗口滑动位移为4的结构,滑动窗口的大小为36的加法分支复用窗口滑动位移为12的结构、滑动窗口的大小为84的加法分支也复用窗口滑动位移为12的结构;这里,84的加法分支、36加法分支与12的加法分支类似,图中仅仅给出了滑动窗口的大小为12的加法分支,通过加法树形结构可以较为方便的推算得到相位偏移值,进一步降低硬件实现复杂度及资源消耗。In the above-mentioned summation calculation of the signal values of the samples in each average block, an additive tree structure is adopted. FIG. 4 is a schematic diagram of an optional addition tree structure in the embodiment of the present invention, as shown in FIG. As shown, the addition branch of 12 is multiplexed with a window sliding displacement of 4, and the signal values of 12 samples are divided into 3 average blocks, and each 4 samples are divided into an average block, for 0-3. The signal value of the sample is summed to obtain sum4_0, and the signal values of the 4th to 7th samples are summed to obtain sum4_1, and the signal values of the 8th to 11th samples are summed to obtain sum4_2; The method is to perform multiplexing on different addition calculation paths. According to different application scenarios, the size of the sliding window can be configured as K1, K2, K3, K4, etc., wherein K1, K2, K3, K4, etc. are positive integers, for example, a sliding window. The size can be configured as four different values of 84, 36, 12, and 4. The size of the sliding window is an integer multiple of 4. When the size of the sliding window is configured as an additive branch of 12, the multiplexed window has a sliding displacement of 4. Structure, sliding window size 36, addition branch multiplexing window sliding The structure with a displacement of 12 and the addition branch with a sliding window size of 84 also multiplex the structure with a sliding displacement of 12; here, the addition branch of 84, the addition branch of 36 is similar to the addition branch of 12, and only the slip is given in the figure. The addition branch of the window size is 12, and the phase offset value can be estimated by the addition tree structure, which further reduces the hardware implementation complexity and resource consumption.
在具体实施过程中,假设滑动窗口的大小为L,将窗口滑动位移采用固定值N,这样既能保证电路的精度足够,又能方便简化。首先,对并行的样点进行相邻N路求和合并,这样就将要处理的样点路数减少,这种预合 并处理方法在满足处理精度需求的条件下,实现了运算量的压缩,并且,此后的处理都是基于合并后的平均块进行,大大降低了后续处理的复杂度;这种在精度允许的条件下将相邻路样点的信号值合并的方案在类似设计中也具有推广意义,在实际应用中,L可以取64、128、256等值,N可以取4、8等值。In the specific implementation process, it is assumed that the size of the sliding window is L, and the sliding displacement of the window is fixed by a value N, which can ensure the accuracy of the circuit is sufficient and can be simplified. First, the parallel N-path summation is performed on the parallel samples, so that the number of sample points to be processed is reduced. The processing method achieves the compression of the computational quantity under the condition that the processing precision requirement is satisfied, and the subsequent processing is performed based on the combined average block, which greatly reduces the complexity of the subsequent processing; The scheme of merging the signal values of adjacent road samples has a generalization significance in similar designs. In practical applications, L can take values of 64, 128, 256, and N can take values of 4 and 8.
由于各样点的信号值包括横坐标值和纵坐标值,为了将每个平均块中样点的信号值进行求和得到每个平均块的信号值,在一种可选的实施例中,对每个平均块中样点的信号值进行求和计算,得到每个平均块的信号值,可以包括:对每个平均块中样点的信号值的横坐标和纵坐标分别进行求和计算,得到每个平均块的信号值的横坐标和每个平均块的信号值的纵坐标;Since the signal values of the various points include the abscissa value and the ordinate value, in order to sum the signal values of the samples in each average block to obtain the signal value of each average block, in an optional embodiment, The summation of the signal values of the samples in each average block to obtain the signal value of each average block may include: summing the abscissa and the ordinate of the signal values of the samples in each average block respectively Obtaining the abscissa of the signal value of each average block and the ordinate of the signal value of each average block;
相应地,S102中对每个平均块的信号值进行低通滤波处理,得到每个平均块低通滤波后的信号值,可以包括:分别对每个平均块的信号值的横坐标和每个平均块的信号值的纵坐标进行低通滤波处理,得到每个平均块低通滤波后的信号值的横坐标和每个平均块低通滤波后的信号值的纵坐标。Correspondingly, the signal value of each average block is subjected to low-pass filtering processing in S102 to obtain a low-pass filtered signal value of each average block, which may include: respectively, the abscissa and the signal value of each average block signal value. The ordinate of the signal value of the average block is subjected to low-pass filtering processing to obtain the abscissa of the low-pass filtered signal value of each average block and the ordinate of the low-pass filtered signal value of each average block.
其中,接收到的各样的信号值是在直角坐标系下的值,所以各样点的信号值中包含有横坐标和纵坐标,这里,分别对各样点的信号值的横坐标和纵坐标进行求和计算,得到每个平均块的信号值的横坐标和每个平均块的信号值的纵坐标;Wherein, the received signal values are values in a Cartesian coordinate system, so the signal values of the samples include the abscissa and the ordinate. Here, the abscissa and the longitudinal direction of the signal values of the respective samples are respectively The coordinates are summed to obtain the abscissa of the signal value of each average block and the ordinate of the signal value of each average block;
这里,对每个平均块的信号值进行低通滤波处理,得到每个平均块低通滤波后的信号值包括对每个平均块的信号值进行去噪处理和消除剩余频偏的处理,也就是说,通过低通滤波处理不仅能校正由于激光器线宽引起的相偏,还可以消除一部分剩余频偏;Here, the signal value of each average block is subjected to low-pass filtering processing, and the signal values obtained by low-pass filtering of each average block include denoising processing of the signal value of each average block and processing for eliminating residual frequency offset, That is to say, the low-pass filtering process can not only correct the phase deviation caused by the line width of the laser, but also eliminate a part of the residual frequency offset;
通过低通滤波处理消除残余频偏为具体原理是:在接收机中,相偏估计和相偏补偿一般是在频偏补偿与频偏估计的后面,即相偏估计和相偏补 偿中接收到的信号是将大部分的频偏去除之后的数据,但还是会有部分残余频偏;由于频偏和相偏都是对相角进行操作,本发明实施例中通过四次方处理、滑动窗口的大小可变的低通滤波操作估计出的相位偏移值是实时值,如果有残余频偏,估计出的相位偏移值也会随之变化,在进行相位偏移值补偿时也就相当于做了一个频偏补偿,从而消除了剩余频偏对信号的损伤。The specific principle of eliminating residual frequency offset by low-pass filtering is: in the receiver, phase offset estimation and phase offset compensation are generally after frequency offset compensation and frequency offset estimation, that is, phase offset estimation and phase offset compensation. The signal received in the reimbursement is the data after the majority of the frequency offset is removed, but there is still a partial residual frequency offset; since the frequency offset and the phase offset are both operated on the phase angle, in the embodiment of the present invention, the fourth power is passed. The phase offset value estimated by the variable-pass low-pass filtering operation of the processing and sliding window is a real-time value. If there is a residual frequency offset, the estimated phase offset value will also change, and the phase offset value compensation is performed. It is equivalent to doing a frequency offset compensation, thereby eliminating the damage of the residual frequency offset to the signal.
由上述加法树形结构可知,根据应用场景不同,滑动窗口的大小可配置为K1、K2、K3、K4等等不同大小,根据剩余频偏的不同自适应采用不同的滑动窗口的大小,当剩余频偏较大时采用较大的滑动窗口的大小来消除剩余频偏,剩余频偏较小时采用较小的滑动窗口的大小来消除剩余频偏。According to the above-mentioned addition tree structure, according to different application scenarios, the size of the sliding window can be configured to different sizes such as K1, K2, K3, K4, etc., and different sliding window sizes are adaptively adapted according to different residual frequency offsets. When the frequency offset is large, the size of the larger sliding window is used to eliminate the residual frequency offset, and when the residual frequency offset is small, the size of the smaller sliding window is used to eliminate the residual frequency offset.
通过低通滤波进行去噪的过程,仍然以四次方处理为例,在对相位进行四次方处理之后得到公式(2),在上述公式(2)中,V4(k)可以去掉信号相位,然后采用多个V4(k)相加,也就是经过低通滤波去除噪声相位,再提取幅角,可得到相位估计的结果,得到的相位估计值θn表示如下:The process of denoising by low-pass filtering is still taking the quadratic processing as an example. After the fourth-order processing of the phase, the formula (2) is obtained. In the above formula (2), V 4 (k) can remove the signal. The phase is then added by multiple V 4 (k), that is, the noise phase is removed by low-pass filtering, and then the angle is extracted to obtain the phase estimation result. The obtained phase estimation value θ n is expressed as follows:
Figure PCTCN2017081739-appb-000001
Figure PCTCN2017081739-appb-000001
由于四次方之后,噪声也被放大了,相角噪声的放大避免不了,为尽量减小幅度噪声被放大,将上式进行优化,得到相位偏移值θe为:Since the noise is also amplified after the fourth power, the amplification of the phase angle noise can not be avoided. To minimize the amplitude noise is amplified, the above equation is optimized to obtain the phase offset value θ e as:
Figure PCTCN2017081739-appb-000002
Figure PCTCN2017081739-appb-000002
在低通滤波处理中,采用可变大小的滑动窗口来对窗口内的各样点的信号值取平均操作,该操作通过对滑动窗口内的样点的信号值求和来完成,为了兼顾处理精度、处理不同应用场景的适应性、尽量降低电路的复杂程度、尽量减小面积和功耗、加快处理速度等多方面的优化需求,本发明实施例提出了一种具有固定窗口滑动位移、支持几种不同滑动窗口的大小、最大限度的增加逻辑复用度的方案,与传统未加优化的做法相比,此种做 法能够保证数据精度的需求,且显著减少电路面积50%以上,提高了电路的处理速度,同时显著降低了功耗。In the low-pass filtering process, a variable-size sliding window is used to average the signal values of the samples in the window, which is done by summing the signal values of the samples in the sliding window, in order to balance the processing. The invention provides a fixed window sliding displacement and supports the accuracy, the adaptability of different application scenarios, the complexity of the circuit, the reduction of the area and the power consumption, and the speed of the processing. Several different sliding window sizes, maximizing the degree of logical multiplexing, compared to traditional unoptimized practices The method can ensure the accuracy of the data, and significantly reduce the circuit area by more than 50%, improve the processing speed of the circuit, and significantly reduce the power consumption.
S103:根据每个平均块低通滤波后的信号值,确定每个平均块的相位偏移值。S103: Determine a phase offset value of each average block according to the low-pass filtered signal value of each average block.
在图3中,将直角坐标系转换成极坐标系,得到(i1,q1)对应的相位偏移值为theta1;In Figure 3, the Cartesian coordinate system is converted into a polar coordinate system, and the phase offset value corresponding to (i1, q1) is obtained as theta1;
经过上述低通滤波处理,得到的每个平均块低通滤波后的信号值中包含有相位偏移值是直角坐标系下的值,为了得到相位偏移值,在一种可选的实施例中,S103可以包括:After the low-pass filtering process, the obtained low-pass filtered signal value of each average block includes a phase offset value which is a value in a Cartesian coordinate system. In order to obtain a phase offset value, in an optional embodiment In S103, the method may include:
将每个平均块低通滤波后的信号值,转换成极坐标下的每个平均块低通滤波后的信号值;将极坐标下的每个平均块低通滤波后的信号值对应的相位值除以M,得出每个平均块的相位偏移值。The low-pass filtered signal value of each average block is converted into a low-pass filtered signal value of each average block in polar coordinates; the phase corresponding to the low-pass filtered signal value of each average block in polar coordinates Dividing the value by M yields the phase offset value for each average block.
在实际应用中,是通过坐标旋转数字计算机(CORDIC,Coordinate Rotation Digital Computer)来实现的,通过CORDIC可以将直角坐标系转换成极坐标系,这样可以得到极坐标系下的每个平均块低通滤波后的信号值对应的相位值,由于各样点进行M次方处理,每个平均块低通滤波后的信号值对应的相位值除以M,得出每个平均块的相位偏移值;如图2所示,在对样点的信号相位进行四次方处理之后,经过低通滤波去噪和消除剩余频偏,通过CORDIC求得幅角,将得到的幅角值除以4才能够得到相位偏移值;本发明实施例中加入CORDIC处理,使原来的直角坐标系下的乘法运算变为极坐标系下的加法运算,大大减少运算量。In practical applications, it is realized by CoDIC Rotation Digital Computer (CORDIC). CORDIC can convert the Cartesian coordinate system into a polar coordinate system, so that each average block low pass in the polar coordinate system can be obtained. The phase value corresponding to the filtered signal value, since the M-th power processing is performed for each sample point, the phase value corresponding to the low-pass filtered signal value of each average block is divided by M, and the phase offset value of each average block is obtained. As shown in Fig. 2, after performing the fourth power processing on the signal phase of the sample, the low-pass filter is used to denoise and eliminate the residual frequency offset, and the amplitude is obtained by CORDIC, and the obtained amplitude value is divided by 4. The phase offset value can be obtained. In the embodiment of the present invention, the CORDIC processing is added to make the multiplication operation in the original Cartesian coordinate system become the addition operation in the polar coordinate system, and the calculation amount is greatly reduced.
在得到每个平均块的相位偏移值之后,为了得到更加精确的相位偏移值,消除求辐角时角度周期性变化导致的不连续性,即避免发生90度相位模糊,需要对相位偏移值进行段内调整,那么,在一种可选的实施例中,在S103之后,在S104之前,该方法可以包括: After obtaining the phase offset value of each average block, in order to obtain a more accurate phase offset value, the discontinuity caused by the periodic variation of the angle when the angle of the argument is obtained is eliminated, that is, the 90 degree phase blur is avoided, and the phase offset is required. The shifting is performed in the intra-segment adjustment. Then, in an optional embodiment, after S103, before S104, the method may include:
在每个并行段中,将每个平均块的相位偏移值与对应的上一平均块的相位偏移值之间的差值确定为第一差值;根据第一差值所落入的预设数值区间,调整每个平均块的相位偏移值。In each parallel segment, a difference between a phase offset value of each average block and a phase offset value of the corresponding previous average block is determined as a first difference; according to the first difference The preset value interval adjusts the phase offset value of each average block.
在图3中,对theta1先进行段内调整;在具体实施过程中,上述预设数值区间可以通过几个阈值来划分成5个预设数值区间;针对每个平均块之间,将前后两个平均块的相位偏移值的差值与几个阈值进行比较,根据比较的结果进行不同的调整操作,实现精细化调整,并行段内的调整具体如下:In Figure 3, theta1 is first adjusted intra-segment; in the specific implementation process, the preset numerical interval can be divided into five preset numerical intervals by several thresholds; for each average block, two before and after The difference between the phase offset values of the average blocks is compared with several thresholds, and different adjustment operations are performed according to the comparison results to achieve fine adjustment, and the adjustments in the parallel segments are as follows:
段内第一个平均块不做调整,从第二个平均块开始,每个平均块的相位偏移值用PEk表示,每个平均块段内调整后的相位偏移值用PEk'∈[-π/4,π/4)表示,为了使得每个平均块段内调整后的相位偏移值用PEk'与前一平均块调整后的相偏估值PEk-1'∈[-π,π)之间的偏差不超过π/4,图5为本发明实施例中对相位偏移值进行调整的一种可选的流程框图,如图5所示,具体调整方法如下:The first average block in the segment is not adjusted. Starting from the second average block, the phase offset value of each average block is represented by PE k , and the adjusted phase offset value in each average block segment is PE k ' ∈[-π/4, π/4) indicates that in order to make the adjusted phase offset value in each average block segment use PE k ' and the previous average block adjusted phase offset estimate PE k-1 '∈ The deviation between [-π, π) does not exceed π/4. FIG. 5 is an optional flow chart of adjusting the phase offset value in the embodiment of the present invention, as shown in FIG. 5, and the specific adjustment method is as follows: :
当PEk<0并且|PEk-PE'k-1|>3π/4时,PE'k=PEk+π;When PE k <0 and |PE k -PE' k-1 |>3π/4, PE' k =PE k +π;
当-3π/4≤(PEk-PE'k-1)<-π/4时,PE'k=PEk+π/2;When -3π/4≤(PE k -PE' k-1 )<-π/4, PE' k =PE k +π/2;
当-π/4≤(PEk-PE'k-1)<π/4时,PE'k=PEkWhen -π/4≤(PE k -PE' k-1 )<π/4, PE' k =PE k ;
当π/4≤(PEk-PE'k-1)<3π/4时,PE'k=PEk-π/2;When π/4≤(PE k -PE' k-1 )<3π/4, PE' k =PE k -π/2;
当PEk≥0并且|PEk-PE'k-1|>3π/4时,PE'k=PEk-π。When PE k ≥ 0 and |PE k -PE' k-1 |> 3π/4, PE' k = PE k - π.
图5中先对PEk和PEk-1'做差值,根据该差值所在范围和PEk所在范围得到段内调整的方法分支,再根据得到的方法分支选择对应的段内调整结果PEk',图5中结构图的输入为相位偏移值,输出为经过段内调整后的相位偏移值PEk'。In FIG. 5, the difference between PE k and PE k-1 ' is firstly performed, and the method branch of the intra-segment adjustment is obtained according to the range of the difference and the range of PE k , and then the corresponding intra-segment adjustment result PE is selected according to the obtained method branch. k ', the input of the structure diagram in Fig. 5 is the phase offset value, and the output is the phase offset value PE k ' after the segment adjustment.
为了进一步的对相位偏移值进行调整以得到更加精确的相位偏移值, 在进行并行段内调整完成后,还需要进行段间调整,那么,在一种可选的实施例中,在根据第一差值所落入的预设数值区间,调整每个平均块的相位偏移值之后,该方法可以包括:将每个并行段中第一个平均块的相位偏移值与对应的上一并行段中最后一个平均块的相位偏移值之间的差值确定为第二差值;根据第二差值所落入的预设数值区间确定修正值;对修正值与每个并行段中第一个平均块的相位偏移值进行求和,得到调整后的第一个并行块的相位偏移值;根据修正值和调整后的第一个并行块的相位偏移值,调整每个并行段中第j个并行块的相位偏移值,其中,j大于等于2。In order to further adjust the phase offset value to obtain a more accurate phase offset value, After the parallel segment adjustment is completed, inter-segment adjustment is also required. Then, in an optional embodiment, the phase of each average block is adjusted according to a preset value interval falling according to the first difference. After the offset value, the method may include: determining a difference between a phase offset value of the first average block in each parallel segment and a phase offset value of a last average block in the corresponding previous parallel segment as a second difference; determining a correction value according to a preset value interval in which the second difference falls; summing the correction value and a phase offset value of the first average block in each parallel segment to obtain an adjusted a phase offset value of a parallel block; adjusting a phase offset value of the jth parallel block in each parallel segment according to the modified value and the adjusted phase offset value of the first parallel block, where j is greater than or equal to 2.
在图3中,对theta1进行段内调整之后,再进行段间调整得到theta2;在具体实施过程中,段间调整可以细化为两个步骤,图6为本发明实施例中对相位偏移值进行调整的另一种可选的流程框图,如图6所示,第一步是并行段的第一个平均块的修正:与段内调整一样,这里不再赘述,由上述段内调整可知,并行段的第一个平均块修正的角度值(相当于上述修正值)的范围如下所示:In FIG. 3, after the intra-segment adjustment of theta1, the inter-segment adjustment is performed to obtain theta2; in the specific implementation process, the inter-segment adjustment can be refined into two steps, and FIG. 6 is a phase shift in the embodiment of the present invention. Another optional flow diagram for adjusting the values, as shown in Figure 6, is the first step of the correction of the first average block of the parallel segment: as with the intra-segment adjustment, it will not be described here, and is adjusted from the above segment. It can be seen that the range of the angle value of the first average block correction of the parallel segment (corresponding to the above correction value) is as follows:
Figure PCTCN2017081739-appb-000003
Figure PCTCN2017081739-appb-000003
第二步是并行段内的二次修正:根据第一个平均块的修正情况,对后一平均块从第二个开始的每一个平均块再次进行修正,其中,PEk'为段内调整后的样点,PEkf'为段间调整后的结果,具体方法如下:The second step is the quadratic correction in the parallel segment: according to the correction of the first average block, the average block of the latter average block from the second one is corrected again, wherein PE k ' is the intra-segment adjustment After the sample, PE kf ' is the result of the adjustment between the segments, the specific method is as follows:
1)当Δ=-π并且PEk'∈[-π,0)时,PE'kf=PE'k+π;1) When Δ=-π and PE k '∈[-π,0), PE' kf =PE' k +π;
2)当Δ=-π并且PEk'∈[0,π)时,PE'kf=PE'k-π;2) When Δ=-π and PE k '∈[0,π), PE' kf =PE' k -π;
3)当Δ=-π/2并且PEk'∈[-π,-π/2)时,PE'kf=PE'k+3π/2;3) When Δ=-π/2 and PE k '∈[-π, -π/2), PE' kf =PE' k +3π/2;
4)当Δ=-π/2并且PEk'∈[-π/2,π)时,PE'kf=PE'k-π/2;4) When Δ=-π/2 and PE k '∈[-π/2,π), PE' kf =PE' k -π/2;
5)当Δ=π/2并且PEk'∈[-π,π/2)时,PE'kf=PE'k+π/2; 5) When Δ = π/2 and PE k '∈[-π, π/2), PE' kf = PE' k + π/2;
6)当Δ=π/2并且PEk'∈[π/2,π)时,PE'kf=PE'k-3π/2;6) When Δ=π/2 and PE k '∈[π/2,π), PE' kf =PE' k -3π/2;
7)当Δ=π并且PEk'∈[-π,0)时,PE'kf=PE'k+π;7) When Δ=π and PE k '∈[-π,0), PE' kf =PE' k +π;
8)当Δ=π并且PEk'∈[0,π)时,PE'kf=PE'k-π;8) When Δ=π and PE k '∈[0,π), PE' kf =PE' k -π;
9)当Δ=0时,PE'kf=PE'k9) When Δ = 0, PE' kf = PE'k;
结合以上9种情况和并行段的第一个平均块的修正方法进行分析,将上述方法分支1)和7)进行合并,2)和8)进行合并,可得包括段间调整第一步和第二步的整体结构图,图6中先对PEk和PEk-1'做差值,得到第二差值,根据该第二差值所在范围求出第一个平均块的修正角度值Δ,得到Δ后结合段内调整后的样点PEk'所在范围得出二次修正的分支,根据计算的分支选择对应的二次修正后的结果PEkf'。图6中结构图的输入为经过段内调整后的数据,输出为经过段间调整后的数据。Combining the above nine cases and the correction method of the first average block of the parallel segment, the above method branches 1) and 7) are combined, and 2) and 8) are combined to obtain the first step and the inter-segment adjustment. In the overall structure diagram of the second step, in FIG. 6, the difference between PE k and PE k-1 ' is first obtained to obtain a second difference, and the corrected angle value of the first average block is obtained according to the range of the second difference. Δ, after obtaining Δ, the range of the adjusted sample PE k ' in the combined segment is obtained as a branch of the second correction, and the corresponding corrected second result PE kf ' is selected according to the calculated branch. The input of the structure diagram in Fig. 6 is the data after the adjustment in the segment, and the output is the data after the inter-segment adjustment.
上述对相位偏移值进行调整的方法,采用了并行段内调整和段间调整两步计算相结合的方法,在调整过程中分区间计算调整值,具有调整精度高的特点。The above method for adjusting the phase offset value adopts a combination of parallel segment adjustment and inter-stage adjustment two-step calculation, and the adjustment value is calculated between the partitions during the adjustment process, and has the characteristics of high adjustment precision.
S104:根据每个平均块的相位偏移值和每个平均块中每个样点的相位值,得到每个平均块中每个样点的校正后相位值。S104: Obtain a corrected phase value of each sample in each average block according to a phase offset value of each average block and a phase value of each sample in each average block.
在一种可选的实施例中,S104可以包括:对每个平均块的相位偏移值、每个平均块中每个样点的相位值以及-π/4求和计算,得到每个平均块中每个样点的校正后相位值。In an optional embodiment, S104 may include: a phase offset value for each average block, a phase value of each sample point in each average block, and a -π/4 summation calculation to obtain each average The corrected phase value for each sample in the block.
在图3中,经过段内调整和段间调整后,可以获得相位偏移值theta2,相位偏移值theta2和经由频偏校正后的信号角度值in_theta做补偿运算,在补偿运算中,每个样点的相位值是与每个样点所在的平均块的相位偏移值来做补偿的,即,在实际应用中,调整后相位偏移值加上频偏校正后的信号角度值再减去固定相位π/4,完成相偏的修正;完成相偏修正后的相角值(theta3)通过余弦cos和正弦sin查找表,得到相应cos(theta3)和sin(theta3) 的值,cos(theta3)与输入的幅度值(in_mod)相乘得到直角坐标系下的输出out_i,sin(theta3)与输入的幅度值(in_mod)相乘得到直角坐标系下的输出out_q,完成由极坐标向直角坐标的转换,从而得到最终的输出。In FIG. 3, after the intra-segment adjustment and the inter-segment adjustment, the phase offset value theta2, the phase offset value theta2, and the signal angle value in_theta after the frequency offset correction are obtained, and in the compensation operation, each of the compensation operations The phase value of the sample is compensated for the phase offset value of the average block where each sample is located, that is, in the actual application, the adjusted phase offset value plus the signal angle value after the frequency offset correction is further reduced. To fix the phase π/4, complete the phase offset correction; complete the phase offset correction phase angle value (theta3) through the cosine cos and sine sin lookup table, get the corresponding cos(theta3) and sin(theta3) The value of cos(theta3) is multiplied by the input amplitude value (in_mod) to obtain the output out_i in the Cartesian coordinate system. The sin(theta3) is multiplied by the input amplitude value (in_mod) to obtain the output out_q in the Cartesian coordinate system. The conversion from polar coordinates to Cartesian coordinates yields the final output.
本发明实施例所提供的相位校正方法,首先对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分成至少一个并行段,然后对每个并行段中样点的信号值进行划分得到预设数目的平均块,再对每个平均块中样点的信号值进行合并得到每个平均块的信号值,这样,避免了对每个样点的信号值都进行计算的缺陷,仅仅只需要对预设数目的平均块进行低通滤波处理来确定每个平均块的相位偏移值,最后,根据每个平均块的相位偏移值和每个平均块中每个样点的相位值,得到每个平均块中每个样点的校正后相位值,也就是说,本发明实施例,通过对每个并行段中样点的信号值划分和合并之后,再针对每个平均块来进行相位偏移值的确定,大大减少了在确定相位偏移值中的计算量,进而也减少了在进行相位校正过程中的运算量,进一步地提高了相位校正的处理效率。The phase correction method provided by the embodiment of the present invention first performs M-th power processing on the signal values of the samples, and divides the signal values of the samples processed by the M-th power into at least one parallel segment, and then parallelizes each of the samples. The signal values of the samples in the segment are divided to obtain a preset number of average blocks, and then the signal values of the samples in each average block are combined to obtain the signal value of each average block, thus avoiding the The defect that the signal values are calculated is only required to perform low-pass filtering processing on the preset number of average blocks to determine the phase offset value of each average block, and finally, according to the phase offset value of each average block and each The phase value of each sample in the average block is obtained as a corrected phase value for each sample in each average block, that is, in the embodiment of the present invention, by dividing the signal value of the sample points in each parallel segment After the combination, the phase offset value is determined for each average block, which greatly reduces the calculation amount in determining the phase offset value, thereby reducing the amount of calculation in the phase correction process, further Phase correction processing efficiency.
基于同一发明构思,本发明实施例还提供一种相位校正装置,图7为本发明实施例中相位校正装置的结构示意图,如图7所示,该装置包括:划分模块71、处理模块72、确定模块73和校正模块74;Based on the same inventive concept, the embodiment of the present invention further provides a phase correction device. FIG. 7 is a schematic structural diagram of a phase correction device according to an embodiment of the present invention. As shown in FIG. 7, the device includes: a division module 71, a processing module 72, Determining module 73 and correction module 74;
其中,划分模块71,配置为对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分至至少一个并行段,其中,M为预设值;处理模块72,配置为对每个并行段中样点的信号值进行划分得到预设数目的平均块,对每个平均块中样点的信号值进行合并得到每个平均块的信号值,对每个平均块的信号值进行低通滤波处理,得到每个平均块低通滤波后的信号值;确定模块73,配置为根据每个平均块低通滤波后的信号值,确定每个平均块的相位偏移值;校正模块74,配置为根据每个平均块的相位偏移值和每个平均块中每个样点的相位值,得到每个平均块中每个样点 的校正后相位值。The dividing module 71 is configured to perform M-th power processing on the signal values of the samples, and divide the signal values of the samples processed by the M-th power into at least one parallel segment, where M is a preset value; The module 72 is configured to divide the signal values of the samples in each parallel segment to obtain a preset number of average blocks, and combine the signal values of the samples in each average block to obtain the signal values of each average block, for each The signal values of the average blocks are subjected to low-pass filtering processing to obtain low-pass filtered signal values for each average block; and the determining module 73 is configured to determine the average block for each low-pass filtered signal value of each average block. a phase offset value; a correction module 74 configured to obtain each sample in each average block based on a phase offset value of each average block and a phase value of each sample in each average block The corrected phase value.
相位校正装置在对各样点的信号值的相位进行M次方处理之后,为了对M次方处理后的各样点的信号值进行划分得到至少一个并行段,在一种可选的实施例中,上述划分模块71,还配置为根据以下一项或多项确定滑动窗口的大小:各样点的信号的传输距离、各样点的信号的激光器线宽;对M次方处理后的各样点的信号值进行划分,得到至少一个并行段,其中,每个并行段的样点的个数为所确定的滑动窗口的大小。After the phase correction device performs the M-th power processing on the phase of the signal value of each sample, the at least one parallel segment is obtained by dividing the signal values of the samples after the M-th power processing, in an optional embodiment. The dividing module 71 is further configured to determine the size of the sliding window according to one or more of the following: the transmission distance of the signals of the samples, the laser line width of the signals of the samples, and the processing after the M-th power processing The signal values of the samples are divided to obtain at least one parallel segment, wherein the number of samples of each parallel segment is the size of the determined sliding window.
相位校正装置在得到预设数目的平均块之后,针对每个平均块进行处理,为了得到每个平均块的信号值,,在一种可选的实施例中,上述处理模块72,还配置为对每个平均块中样点的信号值进行求和计算,得到每个平均块的信号值。The phase correction means, after obtaining a predetermined number of average blocks, is processed for each average block. To obtain the signal value of each average block, in an optional embodiment, the processing module 72 is further configured to The summation of the signal values of the samples in each average block is performed to obtain the signal value of each average block.
由于各样点的信号值包括横坐标值和纵坐标值,为了将每个平均块中样点的信号值进行求和得到每个平均块的信号值,在一种可选的实施例中,上述处理模块72,还配置为对每个平均块中样点的信号值的横坐标和纵坐标分别进行求和计算,得到每个平均块的信号值的横坐标和每个平均块的信号值的纵坐标;分别对每个平均块的信号值的横坐标和每个平均块的信号值的纵坐标进行低通滤波处理,得到每个平均块低通滤波后的信号值的横坐标和每个平均块低通滤波后的信号值的纵坐标。Since the signal values of the various points include the abscissa value and the ordinate value, in order to sum the signal values of the samples in each average block to obtain the signal value of each average block, in an optional embodiment, The processing module 72 is further configured to perform a summation calculation on the abscissa and the ordinate of the signal values of the samples in each average block to obtain the abscissa of the signal value of each average block and the signal value of each average block. The ordinate of the signal; the abscissa of the signal value of each average block and the ordinate of the signal value of each average block are respectively subjected to low-pass filtering processing, and the abscissa and each of the low-pass filtered signal values of each average block are obtained. The ordinate of the low-pass filtered signal value of the average block.
经过上述低通滤波处理,得到的每个平均块低通滤波后的信号值中包含有相位偏移值是直角坐标系下的值,为了得到相位偏移值,在一种可选的实施例中,上述确定模块73,还配置为将每个平均块低通滤波后的信号值,转换成极坐标下的每个平均块低通滤波后的信号值;将极坐标下的每个平均块低通滤波后的信号值对应的相位值除以M,得出每个平均块的相位偏移值。After the low-pass filtering process, the obtained low-pass filtered signal value of each average block includes a phase offset value which is a value in a Cartesian coordinate system. In order to obtain a phase offset value, in an optional embodiment The determining module 73 is further configured to convert the low-pass filtered signal value of each average block into a low-pass filtered signal value of each average block in polar coordinates; each average block in polar coordinates The phase value corresponding to the low-pass filtered signal value is divided by M to obtain the phase offset value of each average block.
在得到每个平均块的相位偏移值之后,为了得到更加精确的相位偏移 值,消除求辐角时角度周期性变化导致的不连续性,即避免发生90度相位模糊,相位校正装置需要对相位偏移值进行段内调整,那么,在一种可选的实施例中,上述装置还包括:第一调整模块(图中未示出),配置为在根据每个平均块低通滤波后的信号值,确定对应的每个平均块的相位偏移值之后,在根据每个平均块的相位偏移值和每个平均块中每个样点的相位值,得到每个平均块中每个样点的校正后相位值之前,在每个并行段中,将每个平均块的相位偏移值与对应的上一平均块的相位偏移值之间的差值确定为第一差值;根据第一差值所落入的预设数值区间,调整每个平均块的相位偏移值。After obtaining the phase offset value of each average block, in order to get a more accurate phase shift Value, eliminating the discontinuity caused by the periodic variation of the angle when the angle is obtained, that is, avoiding the 90 degree phase ambiguity, the phase correction device needs to perform the intra-segment adjustment of the phase offset value, then, in an alternative embodiment The device further includes: a first adjustment module (not shown) configured to determine a phase offset value of each corresponding average block after the low-pass filtered signal value according to each average block, The phase offset value of each average block and the phase value of each sample in each average block, before obtaining the corrected phase value of each sample in each average block, in each parallel segment, each will The difference between the phase offset value of the average block and the phase offset value of the corresponding previous average block is determined as a first difference; each average block is adjusted according to a preset value interval in which the first difference falls Phase offset value.
为了进一步的对相位偏移值进行调整以得到更加精确的相位偏移值,在进行并行段内调整完成后,还需要进行段间调整,那么,在一种可选的实施例中,上述装置还包括:第二调整模块(图中未示出),配置为在根据第一差值所落入的预设数值区间,调整每个平均块的相位偏移值之后,将每个并行段中第一个平均块的相位偏移值与对应的上一并行段中最后一个平均块的相位偏移值之间的差值确定为第二差值;根据第二差值所落入的预设数值区间确定第二差值对应的修正值;对修正值与每个并行段中第一个平均块的相位偏移值进行求和,得到调整后的第一个并行块的相位偏移值;根据修正值和调整后的第一个并行块的相位偏移值,调整每个并行段中第j个并行块的相位偏移值,其中,j大于等于2。In order to further adjust the phase offset value to obtain a more accurate phase offset value, inter-segment adjustment is also required after the parallel segment adjustment is completed. Then, in an optional embodiment, the above device The method further includes: a second adjustment module (not shown) configured to adjust the phase offset value of each average block after the preset value interval falling according to the first difference, in each parallel segment The difference between the phase offset value of the first average block and the phase offset value of the last average block in the corresponding previous parallel segment is determined as the second difference; the preset according to the second difference The value interval determines a correction value corresponding to the second difference; summing the correction value and the phase offset value of the first average block in each parallel segment to obtain a phase offset value of the adjusted first parallel block; Adjusting the phase offset value of the jth parallel block in each parallel segment according to the modified value and the adjusted phase offset value of the first parallel block, where j is greater than or equal to 2.
在一种可选的实施例中,为了确定出校正后的相位值,上述校正模块74,还配置为对每个平均块的相位偏移值、每个平均块中每个样点的相位值以及-π/4求和计算,得到每个平均块中每个样点的校正后相位值。In an optional embodiment, in order to determine the corrected phase value, the correction module 74 is further configured as a phase offset value for each average block, and a phase value of each sample in each average block. And the -π/4 summation calculation yields the corrected phase values for each sample in each average block.
在实际应用中,划分模块71、处理模块72、确定模块73和校正模块74均可由位于装置的中央处理器(CPU,Central Processing Unit)、微处理器(MPU,Microprocessor Unit)、专用集成电路(ASIC,Application Specific  Integrated Circuit)或现场可编程门阵列(FPGA,Field-Programmable Gate Array)等实现。In practical applications, the dividing module 71, the processing module 72, the determining module 73, and the correcting module 74 may all be provided by a central processing unit (CPU), a microprocessor (MPU, a microprocessor unit), an application specific integrated circuit ( ASIC, Application Specific Integrated Circuit) or Field-Programmable Gate Array (FPGA).
本实施例记载一种计算机存储介质,可以为ROM(例如,只读存储器、FLASH存储器、转移装置等)、磁存储介质(例如,磁带、磁盘驱动器等)、光学存储介质(例如,CD-ROM、DVD-ROM、纸卡、纸带等)以及其他熟知类型的程序存储器;计算机可读介质中存储有计算机可执行指令,当执行指令时,引起至少一个处理器执行包括以下的操作:This embodiment describes a computer storage medium, which may be a ROM (for example, a read only memory, a FLASH memory, a transfer device, etc.), a magnetic storage medium (for example, a magnetic tape, a magnetic disk drive, etc.), an optical storage medium (for example, a CD-ROM). , a DVD-ROM, a paper card, a paper tape, etc., and other well-known types of program memory; computer-readable medium storing computer-executable instructions that, when executed, cause at least one processor to perform operations comprising:
对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分至至少一个并行段;对每个并行段中样点的信号值进行划分得到预设数目的平均块,对每个平均块中样点的信号值进行合并得到每个平均块的信号值,对每个平均块的信号值进行低通滤波处理,得到每个平均块低通滤波后的信号值;根据每个平均块低通滤波后的信号值,确定每个平均块的相位偏移值;根据每个平均块的相位偏移值和每个平均块中每个样点的相位值,得到每个平均块中每个样点的校正后相位值。Performing M-th power processing on the signal values of the respective points, dividing the signal values of the samples processed by the M-th power into at least one parallel segment; dividing the signal values of the samples in each parallel segment to obtain a preset number The average block, the signal values of the samples in each average block are combined to obtain the signal value of each average block, and the signal values of each average block are subjected to low-pass filtering processing to obtain low-pass filtered each average block. a signal value; determining a phase offset value of each average block according to a low-pass filtered signal value of each average block; a phase offset value according to each average block and a phase value of each sample point in each average block Get the corrected phase value for each sample in each average block.
本发明实施例所提供的相位校正方法,首先对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分成至少一个并行段,然后对每个并行段中样点的信号值进行划分得到预设数目的平均块,再对每个平均块中样点的信号值进行合并得到每个平均块的信号值,这样,避免了对每个样点的信号值都进行计算的缺陷,仅仅只需要对预设数目的平均块进行低通滤波处理来确定每个平均块的相位偏移值,最后,根据每个平均块的相位偏移值和每个平均块中每个样点的相位值,得到每个平均块中每个样点的校正后相位值,也就是说,本发明实施例,通过对每个并行段中样点的信号值划分和合并之后,再针对每个平均块来进行相位偏移值的确定,大大减少了在确定相位偏移值中的计算量,进而也减少了在进行相位校正过程中的运算量,进一步地提高了相位校正的处理效率。 The phase correction method provided by the embodiment of the present invention first performs M-th power processing on the signal values of the samples, and divides the signal values of the samples processed by the M-th power into at least one parallel segment, and then parallelizes each of the samples. The signal values of the samples in the segment are divided to obtain a preset number of average blocks, and then the signal values of the samples in each average block are combined to obtain the signal value of each average block, thus avoiding the The defect that the signal values are calculated is only required to perform low-pass filtering processing on the preset number of average blocks to determine the phase offset value of each average block, and finally, according to the phase offset value of each average block and each The phase value of each sample in the average block is obtained as a corrected phase value for each sample in each average block, that is, in the embodiment of the present invention, by dividing the signal value of the sample points in each parallel segment After the combination, the phase offset value is determined for each average block, which greatly reduces the calculation amount in determining the phase offset value, thereby reducing the amount of calculation in the phase correction process, further Phase correction processing efficiency.
这里需要指出的是:以上装置实施例项的描述,与上述方法描述是类似的,具有同方法实施例相同的有益效果,因此不做赘述。对于本发明装置实施例中未披露的技术细节,本领域的技术人员请参照本发明方法实施例的描述而理解,为节约篇幅,这里不再赘述。It should be noted here that the description of the above device embodiment items is similar to the above method description, and has the same beneficial effects as the method embodiments, and therefore will not be described again. For the technical details that are not disclosed in the embodiments of the present invention, those skilled in the art should refer to the description of the method embodiments of the present invention, and the details are not described herein.
这里需要指出的是:What needs to be pointed out here is:
应理解,说明书通篇中提到的“一个实施例”或“一实施例”意味着与实施例有关的特定特征、结构或特性包括在本发明的至少一个实施例中。因此,在整个说明书各处出现的“在一个实施例中”或“在一实施例中”未必一定指相同的实施例。此外,这些特定的特征、结构或特性可以任意适合的方式结合在一个或多个实施例中。应理解,在本发明的各种实施例中,上述各过程的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。It is to be understood that the phrase "one embodiment" or "an embodiment" or "an" Thus, "in one embodiment" or "in an embodiment" or "an" In addition, these particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments. It should be understood that, in various embodiments of the present invention, the size of the sequence numbers of the above processes does not mean the order of execution, and the order of execution of each process should be determined by its function and internal logic, and should not be taken to the embodiments of the present invention. The implementation process constitutes any limitation. The serial numbers of the embodiments of the present invention are merely for the description, and do not represent the advantages and disadvantages of the embodiments.
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。It is to be understood that the term "comprises", "comprising", or any other variants thereof, is intended to encompass a non-exclusive inclusion, such that a process, method, article, or device comprising a series of elements includes those elements. It also includes other elements that are not explicitly listed, or elements that are inherent to such a process, method, article, or device. An element that is defined by the phrase "comprising a ..." does not exclude the presence of additional equivalent elements in the process, method, item, or device that comprises the element.
在本申请所提供的几个实施例中,应该理解到,所揭露的设备和方法,可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,如:多个单元或组件可以结合,或可以集成到另一个***,或一些特征可以忽略,或不执行。另外,所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口,设备或单元 的间接耦合或通信连接,可以是电性的、机械的或其它形式的。In the several embodiments provided by the present application, it should be understood that the disclosed apparatus and method may be implemented in other manners. The device embodiments described above are merely illustrative. For example, the division of the unit is only a logical function division. In actual implementation, there may be another division manner, such as: multiple units or components may be combined, or Can be integrated into another system, or some features can be ignored or not executed. In addition, the coupling, or direct coupling, or communication connection of the various components shown or discussed may be through some interface, device or unit. The indirect coupling or communication connection can be electrical, mechanical or other form.
上述作为分离部件说明的单元可以是、或也可以不是物理上分开的,作为单元显示的部件可以是、或也可以不是物理单元;既可以位于一个地方,也可以分布到多个网络单元上;可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。The units described above as separate components may or may not be physically separated, and the components displayed as the unit may or may not be physical units; they may be located in one place or distributed on multiple network units; Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of the embodiment.
另外,在本发明各实施例中的各功能单元可以全部集成在一个处理单元中,也可以是各单元分别单独作为一个单元,也可以两个或两个以上单元集成在一个单元中;上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may be separately used as one unit, or two or more units may be integrated into one unit; The unit can be implemented in the form of hardware or in the form of hardware plus software functional units.
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:移动存储设备、只读存储器(ROM,Read Only Memory)、磁碟或者光盘等各种可以存储程序代码的介质。It will be understood by those skilled in the art that all or part of the steps of implementing the foregoing method embodiments may be performed by hardware related to program instructions. The foregoing program may be stored in a computer readable storage medium, and when executed, the program includes The foregoing steps of the method embodiment; and the foregoing storage medium includes: a removable storage device, a read only memory (ROM), a magnetic disk or an optical disk, and the like, which can store program codes.
或者,本发明上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括:移动存储设备、ROM、磁碟或者光盘等各种可以存储程序代码的介质。Alternatively, the above-described integrated unit of the present invention may be stored in a computer readable storage medium if it is implemented in the form of a software function module and sold or used as a standalone product. Based on such understanding, the technical solution of the embodiments of the present invention may be embodied in the form of a software product in essence or in the form of a software product stored in a storage medium, including a plurality of instructions. A computer device (which may be a personal computer, server, or network device, etc.) is caused to perform all or part of the methods described in various embodiments of the present invention. The foregoing storage medium includes various media that can store program codes, such as a mobile storage device, a ROM, a magnetic disk, or an optical disk.
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明 的保护范围应以所述权利要求的保护范围为准。The above is only a specific embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily think of changes or substitutions within the technical scope of the present invention. It should be covered by the scope of the present invention. Therefore, the present invention The scope of protection shall be subject to the scope of protection of the claims.
工业实用性Industrial applicability
本发明实施例的技术方案,首先对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分成至少一个并行段,然后对每个并行段中样点的信号值进行划分得到预设数目的平均块,再对每个平均块中样点的信号值进行合并得到每个平均块的信号值,这样,避免了对每个样点的信号值都进行计算的缺陷,仅仅只需要对预设数目的平均块进行低通滤波处理来确定每个平均块的相位偏移值,最后,根据每个平均块的相位偏移值和每个平均块中每个样点的相位值,得到每个平均块中每个样点的校正后相位值,本发明实施例通过对每个并行段中样点的信号值划分和合并之后,再针对每个平均块来进行相位偏移值的确定,大大减少了在确定相位偏移值中的计算量,进而也减少了在进行相位校正过程中的运算量,进一步地提高了相位校正的处理效率。 In the technical solution of the embodiment of the present invention, first, the signal value of each sample is processed by M-th power, and the signal value of each sample processed by the M-th power is divided into at least one parallel segment, and then the sample is sampled for each parallel segment. The signal value of the point is divided to obtain a preset number of average blocks, and then the signal values of the samples in each average block are combined to obtain the signal value of each average block, thus avoiding the signal values for each sample point. The defect of calculation is that only a preset number of average blocks need to be low-pass filtered to determine the phase offset value of each average block, and finally, according to the phase offset value of each average block and each average block The phase value of each sample is obtained as the corrected phase value of each sample in each average block. In the embodiment of the present invention, the signal values of the samples in each parallel segment are divided and combined, and then averaged for each average. The block is used to determine the phase offset value, which greatly reduces the calculation amount in determining the phase offset value, thereby reducing the amount of calculation in the phase correction process, and further improving the processing efficiency of the phase correction.

Claims (17)

  1. 一种相位校正方法,包括:A phase correction method comprising:
    对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分至至少一个并行段,其中,M为预设值;Performing M-th power processing on the signal values of the samples, and dividing the signal values of the samples processed by the M-th power into at least one parallel segment, where M is a preset value;
    对每个并行段中样点的信号值进行划分得到预设数目的平均块,对每个平均块中样点的信号值进行合并得到所述每个平均块的信号值,对所述每个平均块的信号值进行低通滤波处理,得到所述每个平均块低通滤波后的信号值;Dividing the signal values of the samples in each parallel segment to obtain a preset number of average blocks, and combining the signal values of the samples in each average block to obtain signal values of each of the average blocks, for each of the The signal values of the average block are subjected to low-pass filtering processing to obtain low-pass filtered signal values of each of the average blocks;
    根据所述每个平均块低通滤波后的信号值,确定所述每个平均块的相位偏移值;Determining, according to the low-pass filtered signal value of each average block, a phase offset value of each of the average blocks;
    根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相位值。A corrected phase value of each of the samples in each of the average blocks is obtained based on a phase offset value of each of the average blocks and a phase value of each of the samples in each of the average blocks.
  2. 根据权利要求1所述的方法,其中,所述将M次方处理后的各样点的信号值划分至至少一个并行段,包括:The method according to claim 1, wherein the dividing the signal value of each sample processed by the Mth power into at least one parallel segment comprises:
    根据以下一项或多项确定滑动窗口的大小:各样点的信号的传输距离、各样点的信号的激光器线宽;Determining the size of the sliding window according to one or more of the following: the transmission distance of the signals of each sample, and the laser line width of the signal of each sample;
    对所述M次方处理后的各样点的信号值进行划分,得到所述至少一个并行段,其中,每个并行段的样点的个数为所确定的滑动窗口的大小。The signal values of the samples processed by the M-th power are divided to obtain the at least one parallel segment, wherein the number of samples of each parallel segment is the size of the determined sliding window.
  3. 根据权利要求1所述的方法,其中,所述对每个平均块中样点的信号值进行合并得到所述每个平均块的信号值,包括:The method of claim 1, wherein the combining the signal values of the samples in each average block to obtain the signal values of each of the average blocks comprises:
    对每个平均块中样点的信号值进行求和计算,得到所述每个平均块的信号值。The summation of the signal values of the samples in each average block is performed to obtain the signal values of each of the average blocks.
  4. 根据权利要求3所述的方法,其中,所述对每个平均块中样点的信号值进行求和计算,得到所述每个平均块的信号值,包括:The method according to claim 3, wherein said summing a signal value of a sample point in each average block to obtain a signal value of said each average block comprises:
    对所述每个平均块中样点的信号值的横坐标和纵坐标分别进行求和 计算,得到所述每个平均块的信号值的横坐标和所述每个平均块的信号值的纵坐标;Sum of the abscissa and the ordinate of the signal values of the samples in each of the average blocks Calculating, obtaining an abscissa of signal values of each of the average blocks and an ordinate of signal values of each of the average blocks;
    相应地,所述对所述每个平均块的信号值进行低通滤波处理,得到所述每个平均块低通滤波后的信号值,包括:分别对所述每个平均块的信号值的横坐标和所述每个平均块的信号值的纵坐标进行低通滤波处理,得到所述每个平均块低通滤波后的信号值的横坐标和所述每个平均块低通滤波后的信号值的纵坐标。Correspondingly, performing low-pass filtering processing on the signal values of each of the average blocks to obtain low-pass filtered signal values of each of the average blocks, including: respectively, signal values of each of the average blocks Performing low-pass filtering processing on the abscissa and the ordinate of the signal value of each of the average blocks, obtaining the abscissa of the low-pass filtered signal value of each of the average blocks and the low-pass filtering of each of the average blocks The ordinate of the signal value.
  5. 根据权利要求1所述的方法,其中,所述根据所述每个平均块低通滤波后的信号值,确定对应的所述每个平均块的相位偏移值,包括:The method according to claim 1, wherein the determining the phase offset value of each of the average blocks according to the low-pass filtered signal value of each of the average blocks comprises:
    将所述每个平均块低通滤波后的信号值,转换成极坐标下的每个平均块低通滤波后的信号值;Converting the low-pass filtered signal value of each average block into a low-pass filtered signal value of each average block in polar coordinates;
    将所述极坐标下的每个平均块低通滤波后的信号值对应的相位值除以M,得出所述每个平均块的相位偏移值。The phase value corresponding to the low-pass filtered signal value of each average block in the polar coordinates is divided by M to obtain the phase offset value of each of the average blocks.
  6. 根据权利要求1所述的方法,其中,在根据所述每个平均块低通滤波后的信号值,确定对应的所述每个平均块的相位偏移值之后,在根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相位值之前,所述方法还包括:The method according to claim 1, wherein after determining a phase offset value of said each of said average blocks based on said low-pass filtered signal value of said average block, Before the phase offset value of the block and the phase value of each sample in each of the average blocks, before obtaining the corrected phase value of each sample in each of the average blocks, the method further includes:
    在所述每个并行段中,将所述每个平均块的相位偏移值与对应的上一平均块的相位偏移值之间的差值确定为第一差值;Determining, in each of the parallel segments, a difference between a phase offset value of each of the average blocks and a phase offset value of the corresponding previous average block as a first difference;
    根据所述第一差值所落入的预设数值区间,调整所述每个平均块的相位偏移值。And adjusting a phase offset value of each of the average blocks according to a preset value interval in which the first difference falls.
  7. 根据权利要求6所述的方法,其中,在根据所述第一差值所落入的预设数值区间,调整所述每个平均块的相位偏移值之后,所述方法还包括: The method according to claim 6, wherein after adjusting the phase offset value of each of the average blocks according to the preset value interval in which the first difference falls, the method further includes:
    将所述每个并行段中第一个平均块的相位偏移值与对应的上一并行段中最后一个平均块的相位偏移值之间的差值确定为第二差值;Determining, by the difference between the phase offset value of the first average block in each parallel segment and the phase offset value of the last average block in the corresponding previous parallel segment as a second difference;
    根据所述第二差值所落入的预设数值区间确定所述第二差值对应的修正值;Determining a correction value corresponding to the second difference according to a preset value interval in which the second difference falls;
    对所述修正值与所述每个并行段中第一个平均块的相位偏移值进行求和,得到调整后的第一个并行块的相位偏移值;And summing the correction value and the phase offset value of the first average block in each parallel segment to obtain a phase offset value of the adjusted first parallel block;
    根据所述修正值和所述调整后的第一个并行块的相位偏移值,调整所述每个并行段中第j个并行块的相位偏移值,其中,j大于等于2。And adjusting a phase offset value of the jth parallel block in each parallel segment according to the modified value and the adjusted phase offset value of the first parallel block, where j is greater than or equal to 2.
  8. 根据权利要求1所述的方法,其中,所述根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相位值,包括:The method according to claim 1, wherein said each of said average blocks is obtained based on a phase offset value of said each average block and a phase value of each of said each of said average blocks The corrected phase values of the samples, including:
    对所述每个平均块的相位偏移值、所述每个平均块中每个样点的相位值以及-π/4求和计算,得到所述每个平均块中每个样点的校正后相位值。Correcting a phase offset value of each of the average blocks, a phase value of each of the samples in each of the average blocks, and a -π/4 summation to obtain a correction for each of the samples in each of the average blocks Post phase value.
  9. 一种相位校正装置,包括:A phase correction device comprising:
    划分模块,配置为对各样点的信号值进行M次方处理,将M次方处理后的各样点的信号值划分至至少一个并行段,其中,M为预设值;The dividing module is configured to perform M-th power processing on the signal values of the samples, and divide the signal values of the samples processed by the M-th power into at least one parallel segment, where M is a preset value;
    处理模块,配置为对每个并行段中样点的信号值进行划分得到预设数目的平均块,对每个平均块中样点的信号值进行合并得到所述每个平均块的信号值,对所述每个平均块的信号值进行低通滤波处理,得到所述每个平均块低通滤波后的信号值;a processing module configured to divide a signal value of a sample point in each parallel segment to obtain a preset number of average blocks, and combine signal values of the sample points in each average block to obtain signal values of each of the average blocks, Performing a low-pass filtering process on the signal values of each of the average blocks to obtain a low-pass filtered signal value of each of the average blocks;
    确定模块,配置为根据所述每个平均块低通滤波后的信号值,确定所述每个平均块的相位偏移值;a determining module, configured to determine a phase offset value of each of the average blocks according to the low-pass filtered signal value of each of the average blocks;
    校正模块,配置为根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相 位值。a correction module configured to obtain a corrected phase of each sample in each of the average blocks according to a phase offset value of each of the average blocks and a phase value of each of the samples in each of the average blocks Bit value.
  10. 根据权利要求9所述的装置,其中,所述划分模块,还配置为根据以下一项或多项确定滑动窗口的大小:各样点的信号的传输距离、各样点的信号的激光器线宽;对所述M次方处理后的各样点的信号值进行划分,得到所述至少一个并行段,其中,每个并行段的样点的个数为所确定的滑动窗口的大小。The apparatus according to claim 9, wherein the dividing module is further configured to determine a size of the sliding window according to one or more of: a transmission distance of a signal of each sample point, and a laser line width of a signal of each sample point And dividing the signal values of the samples processed by the M-th power to obtain the at least one parallel segment, wherein the number of samples of each parallel segment is the size of the determined sliding window.
  11. 根据权利要求9所述的装置,其中,所述处理模块,还配置为对每个平均块中样点的信号值进行求和计算,得到所述每个平均块的信号值。The apparatus according to claim 9, wherein said processing module is further configured to perform a summation calculation on signal values of samples in each average block to obtain signal values of said each average block.
  12. 根据权利要求11所述的装置,其中,所述处理模块,还配置为对所述每个平均块中样点的信号值的横坐标和纵坐标分别进行求和计算,得到所述每个平均块的信号值的横坐标和所述每个平均块的信号值的纵坐标;分别对所述每个平均块的信号值的横坐标和所述每个平均块的信号值的纵坐标进行低通滤波处理,得到所述每个平均块低通滤波后的信号值的横坐标和所述每个平均块低通滤波后的信号值的纵坐标。The apparatus according to claim 11, wherein the processing module is further configured to perform a summation calculation on the abscissa and the ordinate of the signal values of the samples in each of the average blocks to obtain the average The abscissa of the signal value of the block and the ordinate of the signal value of each of the average blocks; respectively lowering the abscissa of the signal value of each of the average blocks and the ordinate of the signal value of each of the average blocks Through the filtering process, the abscissa of the low-pass filtered signal value of each average block and the ordinate of the low-pass filtered signal value of each of the average blocks are obtained.
  13. 根据权利要求9所述的装置,其中,所述确定模块,还配置为将所述每个平均块低通滤波后的信号值,转换成极坐标下的每个平均块低通滤波后的信号值;将所述极坐标下的每个平均块低通滤波后的信号值对应的相位值除以M,得出所述每个平均块的相位偏移值。The apparatus according to claim 9, wherein the determining module is further configured to convert the low-pass filtered signal value of each of the average blocks into a low-pass filtered signal of each average block in polar coordinates. a value; dividing the phase value corresponding to the low-pass filtered signal value of each of the average blocks in the polar coordinates by M to obtain a phase offset value of each of the average blocks.
  14. 根据权利要求9所述的装置,其中,所述装置还包括:The apparatus of claim 9 wherein said apparatus further comprises:
    第一调整模块,配置为在根据所述每个平均块低通滤波后的信号值,确定对应的所述每个平均块的相位偏移值之后,在根据所述每个平均块的相位偏移值和所述每个平均块中每个样点的相位值,得到所述每个平均块中每个样点的校正后相位值之前,在所述每个并行段中,将所述每个平均块的相位偏移值与对应的上一平均块的相位偏移值之间的差值确 定为第一差值;根据所述第一差值所落入的预设数值区间,调整所述每个平均块的相位偏移值。a first adjustment module, configured to determine a phase offset value according to each of the average blocks after determining a phase offset value of each of the corresponding average blocks after the low-pass filtered signal value according to each of the average blocks And shifting the value of each of the samples in each of the average blocks to obtain a corrected phase value for each of the samples in each of the average blocks, in each of the parallel segments, each of said The difference between the phase offset value of the average block and the phase offset value of the corresponding previous average block The first difference is determined; and the phase offset value of each of the average blocks is adjusted according to a preset value interval in which the first difference falls.
  15. 根据权利要求14所述的装置,其中,所述装置还包括:The apparatus of claim 14 wherein said apparatus further comprises:
    第二调整模块,配置为在根据所述第一差值所落入的预设数值区间,调整所述每个平均块的相位偏移值之后,将所述每个并行段中第一个平均块的相位偏移值与对应的上一并行段中最后一个平均块的相位偏移值之间的差值确定为第二差值;根据所述第二差值所落入的预设数值区间确定所述第二差值对应的修正值;对所述修正值与所述每个并行段中第一个平均块的相位偏移值进行求和,得到调整后的第一个并行块的相位偏移值;根据所述修正值和所述调整后的第一个并行块的相位偏移值,调整所述每个并行段中第j个并行块的相位偏移值,其中,j大于等于2。a second adjustment module, configured to adjust a first one of each parallel segment after adjusting a phase offset value of each average block according to a preset value interval that falls within the first difference a difference between a phase offset value of the block and a phase offset value of a last average block in the corresponding previous parallel segment is determined as a second difference; a preset value interval according to the second difference Determining a correction value corresponding to the second difference; summing the correction value with a phase offset value of the first average block in each parallel segment to obtain a phase of the adjusted first parallel block Offset value; adjusting a phase offset value of the jth parallel block in each parallel segment according to the modified value and the adjusted phase offset value of the first parallel block, where j is greater than or equal to 2.
  16. 根据权利要求9所述的装置,其中,所述校正模块,还配置为对所述每个平均块的相位偏移值、所述每个平均块中每个样点的相位值以及-π/4求和计算,得到所述每个平均块中每个样点的校正后相位值。The apparatus of claim 9, wherein the correction module is further configured to phase offset values for each of the average blocks, phase values of each of the samples in each of the average blocks, and -π/ 4 summation calculations to obtain corrected phase values for each of the samples in each of the average blocks.
  17. 一种计算机存储介质,所述计算机存储介质中存储有计算机可执行指令,该计算机可执行指令配置为执行权利要求1-8任一项所述的相位校正方法。 A computer storage medium having stored therein computer executable instructions configured to perform the phase correction method of any of claims 1-8.
PCT/CN2017/081739 2016-09-22 2017-04-24 Phase correction method and device, and computer storage medium WO2018054053A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201610841514.9 2016-09-22
CN201610841514.9A CN107864017B (en) 2016-09-22 2016-09-22 A kind of method for correcting phase and device

Publications (1)

Publication Number Publication Date
WO2018054053A1 true WO2018054053A1 (en) 2018-03-29

Family

ID=61689308

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2017/081739 WO2018054053A1 (en) 2016-09-22 2017-04-24 Phase correction method and device, and computer storage medium

Country Status (2)

Country Link
CN (1) CN107864017B (en)
WO (1) WO2018054053A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112165440A (en) * 2020-09-25 2021-01-01 锐迪科创微电子(北京)有限公司 Phase smoothing method and device, and computer readable storage medium
CN115102817A (en) * 2022-08-24 2022-09-23 鹏城实验室 Phase jump correction method and related equipment
CN115314349A (en) * 2022-06-29 2022-11-08 西安电子科技大学广州研究院 Frequency offset estimation method of wireless system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1303561A (en) * 1998-04-14 2001-07-11 弗兰霍菲尔运输应用研究公司 Method and equipment for multi-carrier modulation and demodulation and method and equipment for performing echo phase offset correction
US20020080898A1 (en) * 2000-07-21 2002-06-27 Broadcom Incorporated Methods and systems for DSP-based receivers
CN1549622A (en) * 2003-05-23 2004-11-24 乐金电子(中国)研究开发中心有限公 Method and apparatus for carrier deviation estimation in mobile communication system
CN102971975A (en) * 2010-07-09 2013-03-13 华为技术有限公司 Method and apparatus for carrier phase estimation and correction in a coherent optical system

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20120068337A (en) * 2010-12-17 2012-06-27 한국전자통신연구원 Method and apparatus for transmitting and receiving coherent optical ofdm
CN104219185B (en) * 2014-08-25 2018-06-15 武汉光迅科技股份有限公司 A kind of line width compensating module and its compensation method applied to DQPSK systems

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1303561A (en) * 1998-04-14 2001-07-11 弗兰霍菲尔运输应用研究公司 Method and equipment for multi-carrier modulation and demodulation and method and equipment for performing echo phase offset correction
US20020080898A1 (en) * 2000-07-21 2002-06-27 Broadcom Incorporated Methods and systems for DSP-based receivers
CN1549622A (en) * 2003-05-23 2004-11-24 乐金电子(中国)研究开发中心有限公 Method and apparatus for carrier deviation estimation in mobile communication system
CN102971975A (en) * 2010-07-09 2013-03-13 华为技术有限公司 Method and apparatus for carrier phase estimation and correction in a coherent optical system

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112165440A (en) * 2020-09-25 2021-01-01 锐迪科创微电子(北京)有限公司 Phase smoothing method and device, and computer readable storage medium
CN115314349A (en) * 2022-06-29 2022-11-08 西安电子科技大学广州研究院 Frequency offset estimation method of wireless system
CN115314349B (en) * 2022-06-29 2023-06-02 西安电子科技大学广州研究院 Frequency offset estimation method of wireless system
CN115102817A (en) * 2022-08-24 2022-09-23 鹏城实验室 Phase jump correction method and related equipment
CN115102817B (en) * 2022-08-24 2022-12-13 鹏城实验室 Phase jump correction method and related equipment

Also Published As

Publication number Publication date
CN107864017B (en) 2019-10-18
CN107864017A (en) 2018-03-30

Similar Documents

Publication Publication Date Title
JP7485133B2 (en) Digital Signal Processing Device
EP3094037B1 (en) Clock recovery method, device and system and computer storage medium
JP5704077B2 (en) Phase deviation / carrier frequency deviation compensation apparatus and phase deviation / carrier frequency deviation compensation method
KR101855019B1 (en) Method for sending/receiving signal, corresponding device and system
WO2016107236A1 (en) Self-adaptive equalization filtering apparatus and method in optical fibre communications, and storage medium
JP6476716B2 (en) Compensation apparatus and compensation method for nonlinear degradation between channels
JP6310091B2 (en) Chromaticity dispersion measuring method, apparatus and digital coherent receiver
WO2017000406A1 (en) Frequency shift and phase shift compensation method and device, and storage medium
JP2017073780A (en) Channel spacing detection device, method and system
WO2015188517A1 (en) Method and device for compensating phase deviation
WO2018054053A1 (en) Phase correction method and device, and computer storage medium
WO2020052683A1 (en) Frequency offset estimation and compensation method and apparatus, communication device, and computer readable storage medium
JP6458730B2 (en) Frequency deviation compensation method, frequency deviation compensation method, and program
CN106101044B (en) A kind of frequency deviation estimating method based on multistep interpolation for coherent optical communication system
WO2018006800A1 (en) Clock balancing method, apparatus and system for qpsk system
CN111262634B (en) Dispersion estimation method, apparatus, receiver and storage medium
WO2019033242A1 (en) Method and apparatus for determining tap coefficient of filter
WO2016101668A1 (en) Carrier phase estimation method and device and storage medium
CN103109481A (en) Nonlinear compensation method, device and signal receiving system
WO2018195722A1 (en) Optical receiver and delay estimation method
JPH10117160A (en) Demodulator
WO2015103805A1 (en) Frequency-domain dispersion search method based on training sequence

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17852127

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 17852127

Country of ref document: EP

Kind code of ref document: A1