A frequency deviation estimation method and device, electronic device, and computer readable medium

Abstract

The present disclosure provides a frequency offset estimation method, which includes: extracting sampling points from an input signal at a set interval to obtain multiple sampling point groups, and the set interval of each sampling point group is different; processing the current sampling point and the extracted sampling points in each sampling point group separately to obtain complex angle data corresponding to the set interval; and determining the frequency offset estimation value of the current input signal using the complex angle data corresponding to the set interval. The frequency offset estimation method does not require switching the frequency offset estimation during the frequency offset estimation process, and can achieve large-range and high-precision frequency offset estimation. The present disclosure also provides a frequency offset estimation device, an electronic device, and a computer-readable medium.

Description

Technical Field

[0001] The embodiments of the present disclosure relate to the field of communication technology, and in particular to a frequency offset estimation method and device, an electronic device, and a computer-readable medium. Background Art

[0002] In an optical communication system, the frequencies of the lasers in a coherent optical transmitter and a coherent optical receiver have frequency deviations and are not completely consistent. It is necessary to estimate the frequency deviation in the coherent optical receiver in order to control or compensate for it.

[0003] Figure 1 It is a common frequency offset estimation device using correlation method. The input complex signal is divided into two paths, one of which is input into register 101, and the other is directly input into multiplier 102. Register 101 delays and conjugates the complex signal to obtain a delayed conjugate signal, which is input into multiplier 102. Multiplier 102 multiplies the complex signal with the conjugate signal to obtain a multiplication result signal, and the multiplication result signal is transmitted to quarticator 103. Quarticator 103 performs a fourth power operation on the multiplication result signal, and then outputs the result to summer 104. Summer 104 is used to reduce the influence of noise on frequency offset estimation. Finally, 1 / 4 angle taker 105 performs a 1 / 4 angle taker operation on the output result of the summer, and outputs a frequency offset estimation signal.

[0004] In practical applications, the delay amount D of register 101 is an integer greater than or equal to 1. When the delay amount D is equal to 1, the frequency offset estimation range is large, but the estimation accuracy is low; when the delay amount D increases, the frequency offset estimation range becomes smaller, and the estimation accuracy is improved. In order to make the frequency offset estimation take into account both a large range and high accuracy, it is usually started with the delay amount D=1. When the frequency offset estimation converges to a certain accuracy, the value of the delay amount D is increased. This single adjustment method often results in the frequency offset estimation not yet reaching the accuracy, and the delay amount D has been switched to a large value, resulting in fuzzy frequency offset estimation and output of erroneous frequency offset estimation results.

[0005] Public Content

[0006] The embodiments of the present disclosure provide a frequency offset estimation method and device, an electronic device, and a computer-readable medium to avoid frequency offset estimation ambiguity caused by frequency offset estimation aliasing during the frequency offset estimation process.

[0007] In a first aspect, an embodiment of the present disclosure provides a frequency offset estimation method, which includes:

[0008] Extracting sampling points from the input signal at set intervals to obtain a plurality of sampling point groups, wherein the set intervals of each sampling point group are different;

[0009] Processing the current sample point and each sample point in the sample point group respectively to obtain complex angle data corresponding to the set interval;

[0010] The frequency offset estimation value of the current input signal is determined using the complex angle data corresponding to the set interval.

[0011] In a second aspect, an embodiment of the present disclosure provides a frequency offset estimation device, comprising:

[0012] A sampling module, used for extracting sampling points from the input signal according to a set interval to obtain a plurality of sampling point groups, and the set interval of each sampling point group is different;

[0013] An angle data acquisition module, used for processing the current sample point and the extracted sample points in each sample point group respectively to obtain the complex angle data corresponding to the set interval;

[0014] The frequency offset estimation module is used to determine the frequency offset estimation value of the current input signal based on the complex angle data corresponding to the set interval.

[0015] In a third aspect, an embodiment of the present disclosure provides an electronic device, comprising:

[0016] one or more processors;

[0017] A memory having one or more programs stored thereon, wherein when the one or more programs are executed by the one or more processors, the one or more processors use any one of the above-mentioned frequency offset estimation methods;

[0018] One or more I / O interfaces are connected between the processor and the memory and are configured to implement information interaction between the processor and the memory.

[0019] In a fourth aspect, an embodiment of the present disclosure provides a computer-readable medium having a computer program stored thereon, wherein the program implements any of the above-mentioned frequency offset estimation methods when executed by a processor.

[0020] The frequency offset estimation method provided by the embodiment of the present disclosure extracts sampling points from the input signal according to the set interval to obtain multiple sampling point groups, and the set interval of each sampling point group is different; the current sampling point and the extracted sampling points in each sampling point group are processed separately to obtain the complex angle data corresponding to the set interval; the frequency offset estimation value of the current input signal is determined based on the complex angle data corresponding to the set interval. The frequency offset estimation method uses sampling point groups with different set intervals to participate in the calculation of the frequency offset estimation value at the same time, and the estimation range and accuracy of the frequency offset can be pre-set. In addition, the calculation structure is fully forward, and there is no need to switch the frequency offset estimation during the frequency offset estimation process, which avoids the frequency offset estimation ambiguity caused by the loop, thereby realizing a large-range and high-precision frequency offset estimation. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] The accompanying drawings are used to provide a further understanding of the embodiments of the present disclosure and constitute a part of the specification. Together with the embodiments of the present disclosure, they are used to explain the present disclosure and do not constitute a limitation of the present disclosure. The above and other features and advantages will become more apparent to those skilled in the art by describing detailed example embodiments with reference to the accompanying drawings, in which:

[0022] Figure 1 It is a common correlation method frequency offset estimation device;

[0023] Figure 2 A flow chart of a frequency offset estimation method provided by an embodiment of the present disclosure;

[0024] Figure 3 A specific flow chart of step 202 in the frequency offset estimation method provided in an embodiment of the present disclosure;

[0025] Figure 4 A specific flow chart of step 203 in the frequency offset estimation method provided in an embodiment of the present disclosure;

[0026] Figure 5 A schematic diagram of frequency offset calculation in the frequency offset estimation method provided in an embodiment of the present disclosure;

[0027] Figure 6 A flowchart of another frequency offset estimation method provided by an embodiment of the present disclosure;

[0028] Figure 7 A schematic diagram of a frequency offset estimation device provided in an embodiment of the present disclosure;

[0029] Figure 8 A schematic diagram of an angle data acquisition module in a frequency deviation estimation device provided in an embodiment of the present disclosure;

[0030] Fig. 9 A schematic diagram of a portion of the structure of a frequency offset estimation device provided in an embodiment of the present disclosure;

[0031] Fig.10 A block diagram of a composition of an electronic device provided in an embodiment of the present disclosure. DETAILED DESCRIPTION

[0032] In order to enable those skilled in the art to better understand the technical solution of the present disclosure, a frequency offset estimation method and device, an electronic device, and a computer-readable medium provided by the present disclosure are described in detail below in conjunction with the accompanying drawings.

[0033] Example embodiments will be described more fully below with reference to the accompanying drawings, but the example embodiments may be embodied in different forms and should not be construed as limited to the embodiments set forth herein. On the contrary, the purpose of providing these embodiments is to make the present disclosure thorough and complete and to enable those skilled in the art to fully understand the scope of the present disclosure.

[0034] In the absence of conflict, the various embodiments of the present disclosure and the various features therein may be combined with each other.

[0035] As used herein, the term "and / or" includes any and all combinations of one or more of the associated listed items.

[0036] The terms used herein are only used to describe specific embodiments and are not intended to limit the present disclosure. As used herein, the singular forms "a", "an" and "the" are also intended to include the plural forms, unless the context clearly indicates otherwise. It will also be understood that when the terms "comprising" and / or "made of" are used in this specification, the presence of the features, wholes, steps, operations, elements and / or components is specified, but the presence or addition of one or more other features, wholes, steps, operations, elements, components and / or groups thereof is not excluded.

[0037] Unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by those of ordinary skill in the art. It will also be understood that terms such as those defined in commonly used dictionaries should be interpreted as having a meaning consistent with their meaning in the context of the relevant art and the present disclosure, and will not be interpreted as having an idealized or overly formal meaning unless explicitly defined as such herein.

[0038] The 400Gb / s PM-16QAM (Polarization Multiplexing-Quadrature Phase Shift Keying) coherent optical communication receiver obtains two signals, namely the X polarization signal and the Y polarization signal, after equalization and polarization demultiplexing. Since the frequency of the laser in the coherent optical communication transmitter is not completely consistent with the frequency of the laser in the coherent optical communication receiver, there is a frequency deviation in the two polarization signals output by the equalizer. It is necessary to estimate the frequency deviation to obtain a frequency deviation estimation signal to compensate for the frequency of the coherent optical communication receiver, thereby improving the signal quality of the optical communication system.

[0039] The frequency offset estimation method of the disclosed embodiment can perform frequency offset estimation on any one of the X polarization signal and the Y polarization signal, and then apply the frequency offset estimation result to the other polarization signal; or, perform frequency offset estimation on both the X polarization signal and the Y polarization signal at the same time. For the convenience of description, the following embodiment introduces the simultaneous frequency offset estimation of the X polarization signal and the Y polarization signal.

[0040] In a first aspect, an embodiment of the present disclosure provides a frequency offset estimation method, which is applied to a coherent optical communication receiver to estimate the frequency offset of a laser in a coherent optical communication transmitter and a laser in a coherent optical communication receiver to compensate or correct the frequency of the coherent optical communication receiver.

[0041] Figure 2 A flowchart of a frequency offset estimation method provided by an embodiment of the present disclosure. Figure 2 , the frequency offset estimation method includes:

[0042] Step 201 , extracting sampling points from an input signal according to a set interval to obtain a plurality of sampling point groups, and the set interval of each sampling point group is different.

[0043] The input signal is a time domain signal obtained by equalizing and filtering the digital baseband electrical signal, and the signal after equalization and filtering is a signal with a parallelism of 64.

[0044] Sampling points are extracted from the input signal, and the input signal can be stored first, and the storage length at least satisfies the number of extracted sampling points, plus the set interval of the extracted sampling points. The set interval is also called the delay interval, which can be set arbitrarily as needed. If the set interval is 4, then one sampling point is extracted every 4 sampling points. In addition, the number of extracted sampling points can also be set arbitrarily as needed, such as extracting 20, 50, or 80 sampling points. More specifically, when one extracted sampling point is extracted every 4 sampling points, and each sampling point group includes 5 extracted sampling points, the storage length should at least include 21 sampling points.

[0045] In this embodiment, the set interval of each sample point group is different. In other words, each sample point group is extracted using a different set interval. For example, for the first sample point group, the set interval may be 1, including 32 extracted sample points. For the second sample point group, the set interval may be 4, including 32 extracted sample points. For the third sample point group, the set interval may be 16, including 32 extracted sample points.

[0046] It should be noted that the sampling objects of each sampling point group in the embodiment of the present disclosure are the same, that is, different sampling point groups are extracted from input signals within the same sampling period.

[0047] Step 202: Process the current sample point and the extracted sample points in each sample point group respectively to obtain complex angle data corresponding to the set interval.

[0048] The current sampling point is any sampling point in the input signal, and the input signal and the input signal of the extracted sampling point group are input signals with the same sampling period.

[0049] In step 202, the current sample point and the sample points in the sample point group are processed respectively, such as fourth power operation, summation, filtering and angle calculation, to obtain complex angle data.

[0050] The current sample point is processed with the extracted sample points in each sample point group to obtain a plurality of complex angle data, that is, each sample point group corresponds to a complex angle data. For example, when there are three sample point groups, three complex angle data can be obtained. Of course, when there are a large number of sample point groups, complex angle data less than the number of sample point groups can also be obtained.

[0051] Step 203: Determine a frequency offset estimation value of the current input signal based on the complex angle data corresponding to the set interval.

[0052] In step 203, the frequency offset estimation value of the current input signal is determined using the complex angle data.

[0053] In some embodiments, the frequency offset estimation value may be estimated step by step by setting the complex angle data corresponding to the interval, and the estimation range of each level is 1 / 4 of the previous level.

[0054] Figure 3 This is a specific flow chart of step 202 in the frequency offset estimation method provided in the embodiment of the present disclosure. Figure 3 As shown, the current sample point and the extracted sample points in each sample point group are processed as follows:

[0055] Step 301, conjugate multiplication, amplitude multiplication and phase subtraction are performed on the current sample point and the extracted sample point in the sample point group, and M-th power operation is performed to obtain an M-th power correlation value.

[0056] In step 301, for the X polarization signal and the Y polarization signal, the current sample point and the sample points in the sample point group are conjugate multiplied, amplitude multiplied and phase subtracted, and an M-th power operation is performed to obtain an M-th power correlation value, where M is an integer greater than 2.

[0057] The current sample point is conjugate multiplied, amplitude multiplied, and phase subtracted with the three sample point groups, and the results are processed to the fourth power to obtain fourth-power correlation values ​​R1, R2, and R3.

[0058] For example, the fourth-power correlation values ​​R1, R2, and R3 are calculated using formula (1).

[0059]

[0060] In formula (1), R1 represents the fourth power correlation value with a setting interval of 1, R2 represents the fourth power correlation value with a setting interval of 4, and R3 represents the fourth power correlation value with a setting interval of 16. represents the signal amplitude, represents the phase of the sampling point, t is the sampling point time, e j4 represents the fourth power processing, x represents the X polarization signal, y represents the Y polarization signal, Indicates the multiplication of the signal amplitude of the X-polarization signal of the current sampling point and the sampling point t+1. Indicates the phase subtraction between the X-polarization signal at the current sample point and the t+1 sample point.

[0061] It should be noted that this embodiment is described by taking the fourth power processing as an example. In fact, the present invention can also perform other power processing.

[0062] Step 302, every N beats, the M-th power correlation values ​​are added once to obtain accumulated data.

[0063] Wherein, N is an integer greater than 2, and N beats represent N M-th power correlation values, that is, N consecutive M-th power correlation values ​​are accumulated. In this embodiment, N can be selected as 32 beats, that is, every 32 consecutive M-th power correlation values ​​are accumulated once.

[0064] For example, 32 fourth-power correlation values ​​are added and summed using formula (2) to obtain accumulated data.

[0065]

[0066] In formula (2), R sum,1 represents the accumulated data corresponding to the sample group with a set interval of 1, R1[count] represents the fourth power correlation value with a set interval of 1, R sum,2 represents the accumulated data corresponding to the sample group with a set interval of 4, R2[count] represents the fourth power correlation value with a set interval of 4, R sum,3 It represents the accumulated data corresponding to the sample group with a set interval of 4, and R3[count] represents the fourth-power correlation value with a set interval of 16.

[0067] Step 303: filter the accumulated data to obtain filtered data.

[0068] In step 303, a first-order IIR digital filter is used to perform low-pass filtering to smooth the accumulated data.

[0069] For example: use formula (3) to calculate the filtered data.

[0070]

[0071] In formula (3), α1, α2, and α3 are the configured filter coefficients, and R avg,1 Indicates the filtered data corresponding to the sample group with a set interval of 1, R avg,2 Indicates the filtered data corresponding to the sample group with a set interval of 4, Ravg,16 Indicates the filtered data corresponding to the sample group with a set interval of 16.

[0072] Step 304: Obtain complex angle data corresponding to the set interval according to the filtered data and the set interval.

[0073] In step 304, a complex angle is calculated according to the filtered data and the set interval.

[0074] For example, the complex angle is calculated using formula (4).

[0075]

[0076] In formula (4), F1, F2, and F3 represent complex angle data, Arg() represents the complex angle to be calculated, the unit is turns, and the value range is usually [-0.5, 0.5] turns. 4, 16, and 64 are values ​​related to the set interval.

[0077] In some embodiments, the frequency offset value is estimated step by step for the complex angle data corresponding to the set interval to obtain the frequency offset estimation value of the current input signal. When estimating the frequency offset value step by step, the estimation range of each step is 1 / 4 of the previous step.

[0078] Figure 4 This is a specific flow chart of step 203 in the frequency offset estimation method provided in the embodiment of the present disclosure. Figure 4 As shown, determining the frequency deviation estimation value of the current input signal based on the complex angle data corresponding to the set interval includes:

[0079] Step 401 , classifying the complex angle data according to a set interval, and the smaller the set interval is, the lower the level of the complex angle data is.

[0080] For example, when the intervals are set to 1, 4, and 16, the complex data angle F1 with an interval of 1 is the first level, the complex data angle F2 with an interval of 4 is the second level, and the complex data angle F3 with an interval of 16 is the third level.

[0081] Step 402, estimating the frequency offset value of the complex angle data step by step to obtain the frequency offset estimation value of the current input signal.

[0082] In some embodiments, a first-level intermediate frequency offset estimation value is obtained based on two adjacent levels of complex angle data; a next-level intermediate frequency offset estimation value is obtained based on the previous-level intermediate frequency offset estimation value and the next-level complex angle data, until a final frequency offset estimation value of the current input signal is obtained.

[0083] The complex angle data of two adjacent levels are directly intercepted at a high position to obtain the angle direct intercept data; based on the angle direct intercept data and the complex angle data, the first-level intermediate frequency deviation estimation value is obtained.

[0084] like Figure 5 As shown, the complex angle data is used to determine the frequency deviation estimate of the current input signal through two adders. Among them, the output bit width of the subtractor is the same as the bit width of its input subtrahend, and the redundant high-end bits are truncated. Assume that the fixed-point bit width of the complex angle is represented as (sign, bit_eff, bit_least), where sign represents the number of signs, bit_eff represents the effective bit width, and bit_least represents the weight of the least significant bit. For example, in the fixed-point bit width (1, 12, 18), "1" represents the number of signs, "12" represents the effective bit width, and "18" represents the least significant bit is 2. -18 .

[0085] In some embodiments, determining the frequency deviation estimation value of the current input signal using the complex angle data includes: n-1 and nth level F n The complex angle value is directly intercepted at the high level to obtain the intercept angle value; the intercept angle value is combined with the n-1th level complex angle value F n-1 Add together to get the intermediate value M of the n-1th level angle n-1 ; The n+1th level complex angle value and the n-1th level angle intermediate value M n-1 High-level direct cut, get a new direct cut angle value; the new direct cut angle value and the middle value M of the n-1th level angle are combined. n-1 The frequency deviation estimation value of the current input signal is finally obtained.

[0086] The following is an introduction using three complex angle values ​​F1, F2, and F3 as examples, wherein the first complex angle value F1 is (1, 12, 14), the second complex angle value F2 is (1, 12, 16), and the third complex angle value F1 is (1, 12, 18).

[0087] After the high-order truncation of F2-F1, the first-level truncation angle value (1, 12, 16) is obtained; the first-level truncation angle value (1, 12, 16) is added to F1 to obtain the intermediate angle value M1, and the fixed-point bit width of the intermediate angle value M1 is (1, 14, 16); after the high-order truncation of F3-M1, the second-level truncation angle value (1, 12, 18) is obtained; after the second-level truncation angle value (1, 12, 18) is added to M1, the intermediate angle value M2 is obtained, and the fixed-point bit width of the intermediate angle value M2 is (1, 16, 18). When there are no other complex angle values, M2 is the final frequency offset estimate. If there are other complex angle values, follow Figure 3 The frequency deviation calculation process continues.

[0088] Figure 6 Flow chart of another frequency offset estimation method provided by an embodiment of the present disclosure. Figure 6As shown, the frequency offset estimation method includes:

[0089] Step 601: pre-process the input signal according to a set threshold to obtain a valid input signal.

[0090] The threshold is set according to the characteristics of the sample points, that is, it is set according to the discrete degree of the input signal.

[0091] In some embodiments, the preprocessing performs circle processing on the input signal. Specifically, a threshold for circle processing is set according to the characteristics of the sample points, and sample points suitable for quartic estimation are found. For example, for a 16QAM input signal, the inner circle 4 sample points and the outermost circle 4 sample points in the constellation diagram are extracted according to the distance of the sample points from the coordinate axis, and the other sample points are set to 0.

[0092] Step 602, conjugate multiplication, amplitude multiplication and phase subtraction are performed on the current sample point and the extracted sample point in the sample point group, and M-th power operation is performed to obtain an M-th power correlation value.

[0093] In step 602, the current sample point and the extracted sample points in the sample point group are all from the sample points after the circle processing. The other processing methods are the same as step 301 in the above embodiment, and will not be repeated here.

[0094] Step 603, every N beats, the M-th power correlation values ​​are added once to obtain accumulated data.

[0095] Step 604: filter the accumulated data to obtain filtered data.

[0096] Step 605, obtaining complex angle data corresponding to the set interval according to the filtered data and the set interval.

[0097] The specific implementation of steps 603 to 605 is the same as that of steps 302 to 304 in the above embodiment, and will not be repeated here.

[0098] Step 606: Determine a frequency offset estimation value of the current input signal based on the complex angle data corresponding to the set interval.

[0099] The specific implementation of step 606 is the same as step 203 and will not be repeated here.

[0100] The frequency offset estimation method provided by the embodiment of the present disclosure extracts sampling points from the input signal according to the set interval to obtain multiple sampling point groups, and the set interval of each sampling point group is different; the current sampling point and the extracted sampling points in each sampling point group are processed separately to obtain the complex angle data corresponding to the set interval; the frequency offset estimation value of the current input signal is determined based on the complex angle data corresponding to the set interval. The frequency offset estimation method uses sampling point groups with different set intervals to participate in the calculation of the frequency offset estimation value at the same time, and the estimation range and accuracy of the frequency offset can be pre-set. In addition, the calculation structure is fully forward, and there is no need to switch the frequency offset estimation during the frequency offset estimation process, which avoids the frequency offset estimation ambiguity caused by the loop, thereby realizing a large-range and high-precision frequency offset estimation.

[0101] In a second aspect, an embodiment of the present disclosure provides a frequency offset estimation device, which can be applied to a coherent optical communication receiver to estimate the frequency deviation of a laser in a coherent optical communication transmitter and a laser in a coherent optical communication receiver to compensate or correct the frequency of the coherent optical communication receiver.

[0102] Figure 7 FIG. 1 is a schematic diagram of a frequency offset estimation device provided by an embodiment of the present disclosure. Figure 7 As shown, the frequency offset estimation device includes:

[0103] The sampling module 701 is used to extract sampling points from the input signal according to a set interval to obtain multiple sampling point groups, and the set interval of each sampling point group is different.

[0104] The specific implementation method of the sampling module 701 is the same as step 201 in the above embodiment, and will not be repeated here.

[0105] The angle data obtaining module 702 is used to process the current sample point and the extracted sample points in each sample point group respectively to obtain complex angle data corresponding to the set interval.

[0106] The specific implementation method of the angle data obtaining module 702 is the same as step 202 in the above embodiment, and will not be repeated here.

[0107] The frequency offset estimation module 703 is used to determine the frequency offset estimation value of the current input signal based on the complex angle data corresponding to the set interval.

[0108] In some embodiments, the frequency offset estimation module 703 estimates the frequency offset value step by step for the complex angle data corresponding to the set interval to obtain the frequency offset estimation value of the current input signal.

[0109] The specific implementation of the frequency offset estimation module 703 is the same as step 203 in the above embodiment, and will not be repeated here.

[0110] like Figure 8 As shown, the angle data acquisition module includes:

[0111] A fourth power calculation unit 801 is used to perform conjugate multiplication, amplitude multiplication and phase subtraction on the current sample point and the extracted sample point in the sample point group, and perform M-th power operation to obtain an M-th power correlation value; wherein M is an integer greater than 2;

[0112] A summing unit 802 is used to add the M-th power correlation value of each N beats once to obtain accumulated data, wherein N is an integer greater than 2;

[0113] The filtering unit 803 is used to filter the accumulated data to obtain filtered data;

[0114] The complex angle calculation unit 804 is used to obtain the complex angle data corresponding to the set interval according to the filtered data and the set interval.

[0115] In some embodiments, the angle data acquisition module further includes a preprocessing unit for preprocessing the input signal according to a set threshold value to obtain a valid input signal.

[0116] In some embodiments, Fig. 9 As shown, in the frequency offset estimation device, the angle data acquisition module includes a plurality of quartic calculation units, summing units, filtering units and complex angle calculation units, and each sample point group corresponds to a quartic calculation unit, summing unit, filtering unit and complex angle calculation unit. For example, when there are three sample point groups, the angle data acquisition module includes a preprocessing unit 901, three quartic calculation units 902, summing units 903, filtering units 904 and complex angle calculation units 905, and a frequency offset estimation module 906. According to the extracted sampling points in a sample group, a set interval of complex angle data is obtained according to the quartic calculation unit 902, summing unit 903, filtering unit 904 and complex angle calculation unit 905, and then the complex angle data is sent to the frequency offset estimation module 906 to calculate the frequency offset estimation value.

[0117] It should be noted that the frequency offset estimation method and device provided in the above embodiments do not require a training sequence, nor do they depend on the synchronous state of the system, and can also be applied in an asynchronous state.

[0118] The frequency offset estimation device provided by the embodiment of the present disclosure has a sampling module that extracts sampling points from the input signal according to a set interval to obtain multiple sampling point groups, and the set interval of each sampling point group is different; the angle data acquisition module processes the current sampling point and the extracted sampling points in each sampling point group respectively to obtain the complex angle data corresponding to the set interval; the frequency offset estimation module determines the frequency offset estimation value of the current input signal based on the complex angle data corresponding to the set interval. The frequency offset estimation device uses sampling point groups with different set intervals to simultaneously participate in the calculation of the frequency offset estimation value, and the estimation range and accuracy of the frequency offset can be pre-set. In addition, the calculation structure is fully forward, and there is no need to switch the frequency offset estimation during the frequency offset estimation process, thereby avoiding the frequency offset estimation ambiguity caused by the loop, thereby realizing a large-range and high-precision frequency offset estimation.

[0119] Thirdly, refer to Fig.10 , an embodiment of the present disclosure provides an electronic device, comprising:

[0120] One or more processors 1001;

[0121] A memory 1002, on which one or more programs are stored, and when the one or more programs are executed by one or more processors, the one or more processors implement any one of the above frequency offset estimation methods;

[0122] One or more I / O interfaces 1003 are connected between the processor and the memory and are configured to implement information exchange between the processor and the memory.

[0123] Among them, the processor 1001 is a device with data processing capabilities, including but not limited to a central processing unit (CPU); the memory 1002 is a device with data storage capabilities, including but not limited to random access memory (RAM, more specifically SDRAM, DDR, etc.), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory (FLASH); the I / O interface (read-write interface) 1003 is connected between the processor 1001 and the memory 1002, and can realize information interaction between the processor 1001 and the memory 1002, including but not limited to a data bus (Bus), etc.

[0124] In some embodiments, the processor 1001 , the memory 1002 , and the I / O interface 1003 are connected to each other through a bus, and further connected to other components of the computing device.

[0125] In a fourth aspect, an embodiment of the present disclosure provides a computer-readable medium having a computer program stored thereon, and when the program is executed by a processor, any one of the above-mentioned frequency offset estimation methods is implemented.

[0126] It will be appreciated by those skilled in the art that all or some of the steps, systems, and functional modules / units in the methods disclosed above may be implemented as software, firmware, hardware, and appropriate combinations thereof. In hardware implementations, the division between the functional modules / units mentioned in the above description does not necessarily correspond to the division of physical components; for example, a physical component may have multiple functions, or a function or step may be performed by several physical components in cooperation. Some or all physical components may be implemented as software executed by a processor, such as a central processing unit, a digital signal processor, or a microprocessor, or implemented as hardware, or implemented as an integrated circuit, such as an application-specific integrated circuit. Such software may be distributed on a computer-readable medium, which may include a computer storage medium (or non-transitory medium) and a communication medium (or temporary medium). As known to those skilled in the art, the term computer storage medium includes volatile and non-volatile, removable and non-removable media implemented in any method or technology for storing information (such as computer-readable instructions, data structures, program modules, or other data). Computer storage media include, but are not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tapes, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to store the desired information and can be accessed by a computer. In addition, it is well known to those skilled in the art that communication media typically contain computer-readable instructions, data structures, program modules, or other data in a modulated data signal such as a carrier wave or other transport mechanism, and may include any information delivery media.

[0127] Example embodiments have been disclosed herein, and although specific terms are employed, they are used and should be interpreted only in a general illustrative sense and not for limiting purposes. In some instances, it will be apparent to those skilled in the art that, unless otherwise expressly noted, features, characteristics, and / or elements described in conjunction with a particular embodiment may be used alone or in combination with features, characteristics, and / or elements described in conjunction with other embodiments. Therefore, those skilled in the art will appreciate that various changes in form and detail may be made without departing from the scope of the present disclosure as set forth in the appended claims.

Claims

1. A frequency offset estimation method, comprising: Extracting sampling points from the input signal according to a set interval to obtain a plurality of sampling point groups, wherein the set interval of each sampling point group is different; Processing the current sample point and the extracted sample points in each of the sample point groups to obtain complex angle data corresponding to the set interval; The current sample point is any sample point in the input signal, and the input signal and the input signal of the sample point group are input signals of the same sampling period; The frequency offset estimation value of the current input signal is determined using the complex angle data corresponding to the set interval.

2. The method according to claim 1, wherein: The processing of the current sample point and the extracted sample points in the sample point group to obtain the complex angle data corresponding to the set interval includes: The current sample point and the extracted sample points in each sample point group are processed as follows: Conjugate multiplication, amplitude multiplication and phase subtraction are performed on the current sample point and the sample points in the sample point group, and M-th power operation is performed to obtain an M-th power correlation value; wherein M is an integer greater than 2; The M-th power correlation value is added once every N beats to obtain accumulated data; wherein N is an integer greater than 2; filtering the accumulated data to obtain filtered data; Complex angle data corresponding to the set interval is obtained according to the filtered data and the set interval.

3. The method according to claim 1, wherein: The determining the frequency deviation estimation value of the current input signal by using the complex angle data corresponding to the set interval includes: classifying the complex angle data according to the set interval; The frequency offset value is estimated step by step for the complex angle data to obtain a frequency offset estimation value of the current input signal.

4. The method according to claim 3, wherein: The step of estimating the frequency offset value of the complex angle data step by step to obtain the frequency offset estimation value of the current input signal includes: Obtaining a first-stage intermediate frequency deviation estimation value based on the complex angle data of two adjacent stages; The intermediate frequency offset estimation value of the next level is obtained based on the intermediate frequency offset estimation value of the previous level and the complex angle data of the next level, until the final frequency offset estimation value of the current input signal is obtained.

5. The method according to claim 4, wherein: The obtaining of a first-level intermediate frequency offset estimation value based on the complex angle data of two adjacent levels comprises: Performing high-order direct cutting of the complex angle data of two adjacent levels to obtain angle direct cutting data; The first-stage intermediate frequency offset estimation value is obtained based on the angle direct data and the complex angle data.

6. The method according to any one of claims 1 to 5, wherein: Before extracting sampling points from the input signal at set intervals to obtain multiple sampling point groups, the method further includes: The input signal is preprocessed according to a set threshold to obtain a valid input signal; wherein the set threshold is set according to the discrete degree of the input signal.

7. A frequency offset estimation device, comprising: A sampling module, used for extracting sampling points from the input signal according to a set interval to obtain a plurality of sampling point groups, and the set interval of each sampling point group is different; An angle data acquisition module, used for processing the current sample point and the extracted sample points in each sample point group respectively to obtain the complex angle data corresponding to the set interval; The current sample point is any sample point in the input signal, and the input signal and the input signal of the sample point group are input signals of the same sampling period; The frequency offset estimation module is used to determine the frequency offset estimation value of the current input signal based on the complex angle data corresponding to the set interval.

8. An electronic device, comprising: one or more processors; A storage device having one or more programs stored thereon, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the method according to any one of claims 1 to 6; One or more I / O interfaces are connected between the processor and the memory and are configured to implement information interaction between the processor and the memory.

9. A computer readable medium having a computer program stored thereon, wherein the program, when executed by a processor, implements the method according to any one of claims 1 to 6.

Images