WO2020135408A1

WO2020135408A1 - Carrier phase estimation method, apparatus, device and computer-readable storage medium

Info

Publication number: WO2020135408A1
Application number: PCT/CN2019/127871
Authority: WO
Inventors: 张良俊
Original assignee: 中兴通讯股份有限公司
Priority date: 2018-12-25
Filing date: 2019-12-24
Publication date: 2020-07-02
Also published as: CN111371502A; CN111371502B

Abstract

Disclosed in the present invention are a carrier phase estimation method and apparatus, a device and a computer-readable storage medium, the method comprising: performing digital signal processing of a signal to be detected acquired in a current symbol, obtaining a phase noise signal only carrying phase noise; on the basis of various pre-configured test angles, rotating the phase noise signal so as to obtain various rotation signals, and on the basis of an error function, acquiring corresponding error values for the various rotation signals; acquiring a minimum error value from among the various error values, and determining a phase angle corresponding to the minimum error value and two phase angles adjacent to said phase angle to the left and right; determining an optimal phase angle on the basis of the minimum error value, the phase angle corresponding to the minimum error value and the two phase angles adjacent to said phase angle to the left and right; acquiring various estimated phase angles before the current symbol, and determining a target phase angle on the basis of the optimal phase angle and the various estimated phase angles.

Description

Carrier phase estimation method, device, equipment and computer readable storage medium

This disclosure claims the priority of the Chinese patent application CN201811596862.X entitled "Carrier Phase Estimation Method, Device, Equipment, and Computer-readable Storage Medium" filed on December 25, 2018, the entire contents of which are incorporated herein by reference .

Technical field

The present disclosure relates to the field of communication technologies, and in particular, to a carrier phase estimation method, device, device, and computer-readable storage medium.

Background technique

In a coherent optical communication system, since the laser at the transmitting end and the laser at the receiving end both have a certain line width, there is inevitably phase noise in the received signal. Phase noise will cause serious degradation of system performance, so we need to correctly estimate and compensate for phase noise.

Coherent digital receivers are commonly used in coherent optical communication system receivers. Coherent digital receivers can compensate for transmission impairments in the received signal in the digital domain, such as chromatic dispersion compensation, polarization mode dispersion compensation, clock recovery, frequency offset compensation, and phase Compensation, etc. The current mainstream phase estimation and compensation algorithms include Vertebi-Vertebi algorithm and blind phase search algorithm. The V-V algorithm is suitable for QPSK (Quadrature Phase Shift Keying) modulation format, but for high-order QAM (Quadrature Amplitude Modulation, quadrature amplitude modulation) compensation capacity is limited. The traditional blind phase search algorithm is transparent to the modulation format, but the algorithm complexity is too high, which is not conducive to hardware implementation

Summary of the invention

The main purpose of the present disclosure is to provide a carrier phase estimation method, device, equipment, and computer storage medium, aimed at solving the technical problem that the current algorithm is too complicated to estimate the phase noise and compensate, which is not conducive to hardware implementation.

To achieve the above objective, the present disclosure provides a carrier phase estimation method, which includes: performing digital signal processing on the signal to be measured obtained in the current symbol to obtain a phase noise signal with only phase noise; Rotate the phase noise signal based on each preset test angle to obtain each rotation signal, and obtain an error value corresponding to each rotation signal based on an error function; obtain a minimum error value among each of the error values, and determine the minimum value The phase angle corresponding to the error value and the two phase angles adjacent to the left and right of the phase angle; based on the phase angle corresponding to the minimum error value and the minimum error value and the two phases adjacent to the left and right of the phase angle The angle determines the optimal phase angle; each estimated phase angle before the current symbol is acquired, and a target phase angle is determined based on the optimal phase angle and each estimated phase angle.

In addition, in order to achieve the above object, the present disclosure also provides a carrier phase estimation device. The carrier phase estimation device includes: a signal processing unit configured to perform digital signal processing on the signal under test acquired in the current symbol to obtain only A phase noise signal with phase noise; a rotation unit for rotating the phase noise signal based on each preset test angle to obtain each rotation signal, and obtaining an error value corresponding to each rotation signal based on an error function; an obtaining unit, It is used to obtain the minimum error value of each of the error values, and determine the phase angle corresponding to the minimum error value and the two phase angles adjacent to the left and right of the phase angle; The phase angle corresponding to the value and the minimum error value and the two phase angles adjacent to the left and right of the phase angle determine the optimal phase angle; the phase angle unit is used to obtain each estimated phase angle before the current symbol, And determine the target phase angle based on the optimal phase angle and each of the estimated phase angles.

In addition, to achieve the above object, the present disclosure also provides a carrier phase estimation device; the carrier phase estimation device includes: a memory, a processor, and a computer program stored on the memory and executable on the processor, Wherein, when the computer program is executed by the processor, the steps of the carrier phase estimation method described above are implemented.

In addition, to achieve the above object, the present disclosure also provides a computer-readable storage medium; the computer-readable storage medium stores a computer program, and when the computer program is executed by the processor, the carrier phase estimation method as described above is implemented step.

BRIEF DESCRIPTION

FIG. 1 is a schematic diagram of a terminal\device structure of a hardware operating environment involved in a solution of an embodiment of the present disclosure;

2 is a schematic flowchart of a first embodiment of a carrier phase estimation method of the present disclosure;

3 is a schematic flowchart of a second embodiment of a carrier phase estimation method of the present disclosure;

4 is a schematic diagram of functional modules of the disclosed carrier phase estimation device;

5 is a schematic diagram of a phase noise estimation device of the present disclosure;

6 is a comparison of the performance of the traditional blind phase search scheme and the disclosed phase noise estimation scheme;

FIG. 7 is a comparison of the performance of the conventional unwinding scheme and the scheme of the present disclosure.

The implementation, functional characteristics and advantages of the present disclosure will be further described in conjunction with the embodiments and with reference to the drawings.

detailed description

It should be understood that the embodiments described herein are only used to explain the present disclosure and are not intended to limit the present disclosure.

As shown in FIG. 1, FIG. 1 is a schematic diagram of a terminal structure of a hardware operating environment involved in a solution of an embodiment of the present disclosure.

The terminal of the embodiment of the present disclosure is a carrier phase estimation device.

As shown in FIG. 1, the terminal may include: a processor 1001, such as a CPU, a network interface 1004, a user interface 1003, a memory 1005, and a communication bus 1002. Among them, the communication bus 1002 is used to implement connection communication between these components. The user interface 1003 may include a display (Display), an input unit such as a keyboard (Keyboard), and the optional user interface 1003 may also include a standard wired interface and a wireless interface. The network interface 1004 may optionally include a standard wired interface and a wireless interface (such as a WI-FI interface). The memory 1005 may be a high-speed RAM memory, or may be a non-volatile memory (non-volatile memory), such as a disk memory. The memory 1005 may optionally be a storage device independent of the foregoing processor 1001.

In one embodiment, the terminal may further include a camera, an RF (Radio Frequency) circuit, a sensor, an audio circuit, a WiFi module, and so on. Among them, sensors such as light sensors, motion sensors and other sensors. In one embodiment, the light sensor may include an ambient light sensor and a proximity sensor, wherein the ambient light sensor may adjust the brightness of the display screen according to the brightness of the ambient light, and the proximity sensor may turn off the display screen when the terminal device moves to the ear And/or backlight. Of course, the terminal device can also be configured with other sensors such as gyroscopes, barometers, hygrometers, thermometers, and infrared sensors, which will not be repeated here.

A person skilled in the art may understand that the terminal structure shown in FIG. 1 does not constitute a limitation on the terminal, and may include more or less components than those illustrated, or combine certain components, or have different component arrangements.

As shown in FIG. 1, the memory 1005 as a computer storage medium may include an operating system, a network communication module, a user interface module, and a carrier phase estimation program.

In the terminal shown in FIG. 1, the network interface 1004 is mainly used to connect to the background server and perform data communication with the background server; the user interface 1003 is mainly used to connect to the client (user) and perform data communication with the client; and the processor 1001 can be used to call the carrier phase estimation program stored in the memory 1005 and perform the following operations: perform digital signal processing on the signal under test acquired in the current symbol to obtain a phase noise signal with only phase noise; Set the test angle to rotate the phase noise signal to obtain each rotation signal, and obtain the error value corresponding to each rotation signal based on the error function; obtain the minimum error value among each of the error values, and determine the corresponding minimum error value Phase angle and two phase angles adjacent to the left and right of the phase angle; based on the minimum error value and the phase angle corresponding to the minimum error value and two phase angles adjacent to the left and right of the phase angle to determine the most Optimal phase angle; acquiring each estimated phase angle before the current symbol, and determining a target phase angle based on the optimal phase angle and each of the estimated phase angles.

Referring to FIG. 2, the present disclosure provides a carrier phase estimation method. In an embodiment of the carrier phase estimation method, the carrier phase estimation method includes the following steps:

Step S10: Perform digital signal processing on the signal under test acquired in the current symbol to obtain a phase noise signal with only phase noise;

A symbol can be a binary number represented by a symbol with the same time interval in digital communication. Phase noise can refer to the random changes in the phase of the system output signal caused by the system (such as various RF devices) under the action of various noises, and is an important indicator to measure the quality of the frequency standard source frequency stability. In the current symbol, the front-end digital signal processing is performed by the signal under test received by the receiving end, including signal processing such as delay adjustment, DC removal, dispersion compensation, clock synchronization, polarization demultiplexing, frequency offset compensation, etc. Phase noise signal with only phase noise.

Step S20: Rotate the phase noise signal based on each preset test angle to obtain each rotation signal, and obtain an error value corresponding to each rotation signal based on an error function;

The preset test angle may be each test angle set by the user in advance. After obtaining the preset test angle, it is necessary to use the preset test angle to rotate the phase noise signal to obtain its corresponding rotation signal, and calculate the error value corresponding to the rotation function by means of the error function calculation. It should be noted that each preset test angle has a corresponding rotation signal, that is, an error value. The way to calculate the error value by using the error function may be to first select a modulation format, and then determine the constellation map after the rotation signal is mapped according to the modulation format, and determine the mapping value (D1, D2, D3) corresponding to the rotation signal in the constellation map ), and it is also necessary to obtain the real and imaginary values in the rotation signal at this time. Then subtract the mapping value D1 from the absolute value of the real part value to get the first real part difference, and use the absolute value of the first real part difference to subtract the mapping value D2 to get the second real part Difference, and use the absolute value of the second real difference to subtract the mapping value D3 to get the third real difference, and determine the absolute real difference of the third real difference, the same reason, The imaginary part value in the rotation signal is also subtracted from the mapping value (D1, D2, D3) to obtain the third imaginary part difference, and the absolute imaginary part difference of this third imaginary part difference is determined, and then The absolute real part difference and the absolute real part difference are added to obtain the difference calculated this time. In addition, in order to reduce the influence of Gaussian white noise, the error calculated by each symbol and the errors calculated by several symbols before and after it are summed and averaged to obtain the error value corresponding to the rotation signal.

To assist in understanding the method of error calculation in the present disclosure, an example will be described below.

For example, first use a set of test angles to rotate the signal to be restored, ie

Test corner

Select at equal intervals of (-π/4, π/4), the expression is as follows:

b ∈ (-B/2, -B/2+1,...0,1...B/2-1). Where B is the number of selected test angles. Then use the error function to calculate the error of each rotation signal Sk, b. The calculation formula is as follows: I1=|real(Sk, b)|-D1, I2=|I1|-D2, I3=|I2|-D3; Q1=|Imag(Sk, b)|-D1, Q2=|Q1 |-D2, Q3=|Q2|-D3; e=|I3|+|Q3|. Where real(Sk, b) represents the real part of the signal to be measured, Imag(Sk, b) represents the imaginary part of the signal to be tested, and e is the error value. D1, D2, and D3 are mapping values, and the values of D1, D2, and D3 are related to the modulation format. For square QPSK (Quadrature Phase Shift Keying), D1, D2, and D3 take the

values

1, 0, and 0, respectively. For square 16QAM (Quadrature Amplitude Modulation, quadrature amplitude modulation), D1, D2, and D3 take the

values

2, 1, 0, respectively. For square 32QAM and 64QAM, D1, D2, and D3 take the

values

4, 2, and 1, respectively.

Step S30: Obtain the minimum error value among each of the error values, and determine the phase angle corresponding to the minimum error value and the two phase angles adjacent to the left and right of the phase angle;

When each error value is obtained, the minimum error value of each error value needs to be determined, and then the phase angle with the smallest error is determined according to the minimum error value

And the two phase angles adjacent to it

with

And the relationship between these three phase angles should satisfy

The corresponding error values are e1, e2, and e3, which satisfy the relationship e1﹤e2﹤e3. It should also be noted that if the phase angle with the smallest error

Is the first phase angle, then

Take the last phase angle, e1 takes the corresponding error value. in case

Is the last phase angle, then

Take the first phase angle, e3 takes the corresponding error value. The phase angle may be phase noise.

Step S40: Determine an optimal phase angle based on the minimum error value and the phase angle corresponding to the minimum error value and two phase angles adjacent to the left and right of the phase angle;

When the minimum error value and the phase angle with the smallest error and the two phase angles adjacent to the left and right of the phase angle are obtained, interpolation calculation can be performed to obtain the optimal phase angle μ. For example, taking the quadratic curve interpolation as an example, the calculation formula is as follows:

Where B is the number of test corners.

Step S50: Obtain each estimated phase angle before the current symbol, and determine a target phase angle based on the optimal phase angle and each estimated phase angle.

After the optimal phase angle in the current symbol is obtained, the optimal phase angle needs to be unwound. The method is to save the phase estimation angle of N symbols before the current symbol and sum them to get the average

For μ and

Make judgments if

Then μ′=μ-π/2; if

Then μ′=μ+π/2. The final μ′ obtained is the phase noise estimated by this scheme, that is, the target phase angle. Among them, the effect of the unwinding operation is exemplified. For example, as shown in FIG. 7, taking the 16QAM modulation format as an example, comparing the effect of the traditional unwinding scheme with this embodiment, it can be seen from the figure that the traditional Unwrapping schemes are prone to frequently generate phase cycle slips, and adopting the disclosed scheme can effectively reduce the probability of phase cycle slips.

To assist in understanding the phase noise estimation method in the present disclosure, an example will be described below.

For example, as shown in Fig. 5, when the signal to be measured is acquired, the signal to be measured needs to be rotated according to the test angle (test angle 1, test angle 2...test angle B), and the various functions are calculated according to the error function calculation module The error value corresponding to the test angle, and then select the phase angle corresponding to the minimum error value and the two adjacent phase angles through the phase angle selection module, and calculate the optimal phase angle according to the interpolation calculation module, and finally unwind In the module, the optimal phase angle estimated by the current symbol and the average value of the previous N symbols are unwrapped to obtain the final phase angle, which is the target phase angle.

In addition, in order to assist in understanding the technical effects of the present disclosure, examples will be described below.

As shown in Fig. 6, taking the 16QAM modulation format as an example, the performance of the traditional blind phase search scheme is compared with this scheme. It can be seen from the figure that compared to the traditional blind phase search scheme, which requires 32 test angles, this solution only requires 8 test angles to achieve the same performance.

In this embodiment, by performing digital signal processing on the signal under test acquired in the current symbol, a phase noise signal with only phase noise is obtained; the phase noise signal is rotated based on each preset test angle to obtain each rotation signal, and Obtaining an error value corresponding to each of the rotation signals based on an error function; obtaining a minimum error value among each of the error values, and determining a phase angle corresponding to the minimum error value and two phases adjacent to the phase angle Angle; determine the optimal phase angle based on the minimum error value and the phase angle corresponding to the minimum error value and two phase angles adjacent to the left and right of the phase angle; obtain each estimated phase angle before the current symbol And determine the target phase angle based on the optimal phase angle and each of the estimated phase angles. By using the error function to calculate the corresponding error value for the obtained signal under test, compared with the traditional calculation method using a large number of multipliers, the number of test angles is reduced, and the calculation complexity of the error function is greatly reduced , And in determining the target phase angle, that is, the estimated phase noise, in addition to estimating the current optimal phase angle, it is also related to the previously obtained estimated phase angle, thereby effectively reducing the probability of phase slip, and improving the system The performance has achieved the technical effect of greatly reducing the complexity of the algorithm without reducing the accuracy of phase noise compensation, which is conducive to hardware implementation.

Based on the first embodiment of the present disclosure, a second embodiment of the carrier phase estimation method of the present disclosure is proposed. This embodiment is a refinement of step S50 of the first embodiment of the present disclosure. Referring to FIG. 3, including:

Step S51: Obtain historical symbols before the current symbol, and acquire estimated phase angles corresponding to the historical symbols;

The historical symbol may be a symbol whose carrier phase has been estimated before the current symbol. The estimated phase angle may be the estimated phase noise in the historical symbols, and each historical symbol has an estimated phase angle corresponding to it. After obtaining the optimal phase angle of the current symbol, it is also necessary to obtain each historical symbol before the current symbol, and then obtain the estimated phase angle corresponding to each historical symbol.

Step S52: Obtain an average value of the phase angle between the optimal phase angle and each of the estimated phase angles, and determine a target phase angle based on the optimal phase angle and the average phase angle.

The average value of the phase angle may be an average value obtained by averaging the optimal phase angle and each estimated phase angle. After obtaining the optimal phase angle and each estimated phase angle, determine the total number of the optimal phase angle and each estimated phase angle, and calculate the average value of the optimal phase angle and each estimated phase angle according to the total number. The phases of the N symbols before the symbols estimate the angles and perform a summation and averaging process to obtain an average value. In addition, in order to reduce the influence of Gaussian white noise, the error calculated by each symbol and the errors calculated by several symbols before and after it are summed and averaged. To get the target phase angle.

In this embodiment, the target phase angle is determined by calculating the average of the phase angle between the optimal phase angle and each estimated phase angle to ensure the accuracy of carrier phase estimation.

In one embodiment, the step of determining the target phase angle based on the optimal phase angle and the average of the phase angles includes:

Step S521: Obtain the maximum preset test angle and the minimum preset test angle in each of the preset test angles, and compare the optimal phase angle with the average value of the phase angles;

In each preset test angle input by the user, the maximum preset test angle and the minimum preset test angle need to be determined, and the optimal phase angle and the average value of the phase angles need to be compared and judged, that is, the optimal phase angle is judged Whether it is greater than the sum of the average value of the phase angle and the maximum preset test angle, if it is greater, you can obtain the difference between the optimal phase angle minus the preset value, and use this difference as the target phase angle; determine the optimal Whether the phase angle is less than the sum of the average value of the phase angle and the minimum preset test angle, if it is less, the sum of the optimal phase angle and the preset value can be obtained, and the sum is used as the target phase angle. For example, when the optimal phase angle is μ, the average phase angle is

And test corner

If (-π/4, π/4) is selected at equal intervals, then for μ and

Make judgments if

Then μ′=μ-π/2; if

Then μ′=μ+π/2. The resulting μ′ is the phase noise estimated by this scheme.

Step S522: If the optimal phase angle is greater than the first sum value between the maximum preset test angle and the average value of the phase angles, obtain the first between the optimal phase angle and the preset angle The difference, and use the first difference as the target phase angle;

The first sum value is the sum of the maximum preset test angle and the average value of the phase angles. The first difference is the difference between the optimal phase angle and the preset angle. The preset angle can be an angle set by the user. When it is judged that the optimal phase angle is greater than the first sum value between the maximum preset test angle and the average value of the phase angles, the first difference between the optimal phase angle and the preset angle can be obtained, and the The first difference is taken as the target phase angle.

Step S523: If the optimal phase angle is less than the second sum between the average value of the phase angles and the minimum preset test angle, obtain the difference between the optimal phase angle and the preset angle The third sum value, and use the third sum value as the target phase angle.

The second sum value may be the sum value between the average value of the phase angle and the minimum preset test angle. The third sum value may be the sum value between the optimal phase angle and the preset angle. When it is judged that the optimal phase angle is less than the second sum between the average phase angle and the minimum preset test angle, the third sum between the optimal phase angle and the preset angle can be obtained, and this The third sum value serves as the target phase angle.

In this embodiment, the target phase angle is determined by comparing the optimal phase angle with the average value of the phase angle, thereby ensuring the accuracy of obtaining the target phase angle and effectively reducing the probability of phase slip.

Based on any one of the first to second embodiments of the present disclosure, a third embodiment of the carrier phase estimation method of the present disclosure is proposed. This embodiment is step S20 of the first embodiment of the present disclosure. The refinement of the steps of the error value corresponding to the rotation signal includes:

Step S21: Obtain the modulation format corresponding to the rotation signal, and determine the mapping value in the modulation format;

Modulation format can choose square QPSK, square 16QAM, square 32QAM, square 64QAM and square 8QAM. The mapping value may be a value obtained by mapping the rotation signal through a certain modulation format, and this value is associated with the real and imaginary parts of the rotation signal. Obtain the modulation format corresponding to the rotation signal (the modulation format can be selected by the user), and then determine each mapping value according to this modulation format.

Step S22: Obtain the real part value and the imaginary part value in the rotation signal, and calculate the real part difference between the real part value and the map value according to an error function, the imaginary part value and the map Imaginary difference between values;

Obtain the real and imaginary values in this rotation signal, and calculate the real difference between the real and mapped values through the error function calculation, and calculate the imaginary difference between the imaginary and mapped values Value, and after each calculation, you need to take the calculated absolute value.

Step S23: Obtain an error sum value between the absolute real part difference corresponding to the real part difference and the absolute imaginary part difference corresponding to the imaginary part difference, and determine the rotation signal based on the error sum The corresponding error value.

The error sum value may be the error value calculated by this error function. Obtain the absolute real part difference corresponding to the real part difference and the absolute imaginary part difference corresponding to the imaginary part difference, and calculate the error sum value between the absolute real part difference and the absolute imaginary part difference, and then pass this error sum Value to determine the error value corresponding to the rotation signal.

In this embodiment, the error value corresponding to the rotation signal is determined by determining the modulation format, thereby effectively ensuring the accuracy of the obtained error value of the signal, and improving the user's sense of experience.

In one embodiment, the step of determining the error value corresponding to the rotation signal based on the sum value includes:

Step S231: Acquire the primary error value corresponding to the error sum value and the preset number of intermediate error values before the current symbol;

The primary error value may be an error value calculated by an error function, and the primary error value and the error sum value are equal. The intermediate error value may be an error value calculated before the current symbol. The preset quantity may be a quantity set by the user. After the error sum value is calculated by the error function, the primary error value corresponding to this error sum value needs to be determined, and the intermediate error value calculated by the symbol before the current symbol should be obtained. The number of the intermediate error value and the value set by the user in advance The preset number is the same.

Step S232: Obtain an average value of errors between the primary error value and each of the intermediate error values, and use the average value of the errors as an error value corresponding to the rotation signal.

After the primary error value and the preset number of intermediate error values are obtained, the average value of the error between the primary error value and each intermediate error value needs to be calculated, and the average error value is used as the error value corresponding to the rotation signal. It should be noted that every time the rotation signal is obtained, the corresponding error average value needs to be obtained.

In this embodiment, by obtaining the average value of the error between the primary error value and each intermediate error value, and using this average value as the error value corresponding to the rotation signal, the influence of the Gaussian white noise is reduced and relative to The traditional error calculation function is to calculate the Euclidean distance between the rotated signal and the corresponding constellation point. This embodiment can effectively reduce the calculation complexity.

Based on any one of the first to third embodiments of the present disclosure, a fourth embodiment of the carrier phase estimation method of the present disclosure is proposed. This embodiment is step S30 of the first embodiment of the present disclosure to obtain each of the error values Refinement of the minimum error value in the step of determining the phase angle corresponding to the minimum error value and the two phase angles adjacent to the left and right of the phase angle, including:

Step S31: Obtain the minimum error value among each of the error values, and determine the phase angle corresponding to the minimum error value;

After the error value corresponding to each rotation signal is obtained through calculation, the minimum minimum error value among each error value needs to be obtained, and the phase angle corresponding to the minimum error value is determined according to the minimum error value.

Step S32, judging whether the phase angle corresponding to the minimum error value in the current symbol is the primary phase angle for the first test;

The primary phase angle may be the preset test phase angle used when the current symbol is tested for the first time. After obtaining the phase angle corresponding to the minimum error value, it is necessary to determine whether this phase angle is the primary phase angle for the first test.

Step S33, if the phase angle corresponding to the minimum error value is the primary phase angle for the first test, obtain the final phase angle of the last test in the current symbol, and divide the second phase in the current symbol The second phase angle and the final phase angle of the second test are taken as two phase angles adjacent to the left and right of the phase angle.

When it is judged that the phase angle corresponding to the minimum error value is the primary phase angle for the first test, you need to obtain the final phase angle used in the last test in the current symbol and the second phase angle for the second test. The second phase angle, and the second phase angle and the final phase angle are taken as two adjacent phase angles. Wherein, the ultimate phase angle may be the phase angle obtained during the last test in the current symbol. The second phase angle may be the phase angle obtained during the second test in the current symbol.

In this embodiment, by judging whether the phase angle corresponding to the minimum error value is the primary phase angle obtained during the first test, and when the phase angle corresponding to the minimum error value is the primary phase angle, the second phase angle The phase angle corresponding to the final phase angle as the minimum error value is adjacent to the left and right, thereby effectively ensuring that the phase angle corresponding to the minimum error value can be obtained at any time, and the carrier phase is improved. Estimated accuracy.

In one embodiment, after the step of determining whether the phase angle corresponding to the minimum error value is the primary phase angle for the first test, it includes:

Step S34, if the phase angle corresponding to the minimum error value is not the primary phase angle for the first test, determine whether the phase angle corresponding to the minimum error value is the final phase angle for the last test;

When it is judged that the phase angle corresponding to the minimum error value is not the primary phase angle obtained during the first test, it is necessary to judge again whether the phase angle corresponding to the minimum error value is the final phase angle obtained during the last test .

Step S35, if the phase angle corresponding to the minimum error value is the final phase angle of the last test, obtain the primary phase angle of the first test, and perform the primary phase angle and the penultimate time The tested third phase angle is taken as two phase angles adjacent to the left and right of the phase angle.

When it is judged that the phase angle corresponding to the minimum error value is the final phase angle obtained during the last test, you need to obtain the primary phase angle obtained during the first test and the penultimate phase angle obtained during the second test. For the third phase angle, the primary phase angle and the third phase angle are regarded as the phase angle corresponding to the minimum error value, and the two adjacent phase angles are left and right. However, when it is judged that the phase angle corresponding to the minimum error value is not the final phase angle obtained during the last test, the two phase angles adjacent to the left and right of the phase angle are directly obtained. Among them, the third phase angle is the symbol obtained during the penultimate test of the current symbol.

In this embodiment, by judging whether the phase angle corresponding to the minimum error value is the final phase angle during the last test, the accuracy of the carrier phase estimation is effectively ensured, and the user experience is improved.

In addition, referring to FIG. 4, an embodiment of the present disclosure also proposes a carrier phase estimation device. The carrier phase estimation device includes: a signal processing unit configured to perform digital signal processing on the signal to be measured obtained in the current symbol to obtain A phase noise signal with only phase noise; a rotation unit for rotating the phase noise signal based on each preset test angle to obtain each rotation signal, and obtaining an error value corresponding to each rotation signal based on an error function; an obtaining unit , Used to obtain the minimum error value of each of the error values, and to determine the phase angle corresponding to the minimum error value and the two phase angles adjacent to the phase angle to the left and right; The phase angle corresponding to the error value and the minimum error value and the two phase angles adjacent to the left and right of the phase angle determine an optimal phase angle; a phase angle unit is used to obtain each estimated phase angle before the current symbol And determine the target phase angle based on the optimal phase angle and each of the estimated phase angles.

In one embodiment, the phase angle unit is configured to: obtain each historical symbol before the current symbol, and obtain each estimated phase angle corresponding to each historical symbol; and obtain the optimal phase angle And an average value of phase angles between each of the estimated phase angles, and a target phase angle is determined based on the optimal phase angle and the average value of phase angles.

In one embodiment, the phase angle unit is configured to: obtain a maximum preset test angle and a minimum preset test angle in each of the preset test angles, and combine the optimal phase angle and the phase angle Compare the average values; if the optimal phase angle is greater than the first sum between the maximum preset test angle and the average value of the phase angles, then obtain the difference between the optimal phase angle and the preset angle A first difference, and use the first difference as a target phase angle; if the optimal phase angle is less than the second sum between the average phase angle and the minimum preset test angle, obtain A third sum value between the optimal phase angle and the preset angle, and using the third sum value as the target phase angle.

In one embodiment, the rotation unit is further configured to: acquire a modulation format corresponding to the rotation signal, and determine a mapping value in the modulation format; and acquire real value and imaginary value in the rotation signal And calculate the real difference between the real value and the mapped value according to the error function, the imaginary difference between the imaginary value and the mapped value; get the real difference corresponding to An error sum value between the absolute real part difference value and the absolute imaginary part difference value corresponding to the imaginary part difference value, and the error value corresponding to the rotation signal is determined based on the error sum value.

In one embodiment, the rotating unit is further configured to: obtain a primary error value corresponding to the error sum value and a preset number of intermediate error values before the current symbol; obtain the primary error value and each An average value of errors between the intermediate error values, and the average value of the errors is used as an error value corresponding to the rotation signal.

In an embodiment, the acquiring unit is further configured to: acquire a minimum error value among each of the error values, and determine a phase angle corresponding to the minimum error value; determine the minimum value in the current symbol Whether the phase angle corresponding to the error value is the primary phase angle for the first test; if the phase angle corresponding to the minimum error value is the primary phase angle for the first test, obtain the last test in the current symbol The final phase angle of the current symbol, and use the second phase angle of the second test in the current symbol and the final phase angle as the two phase angles adjacent to the left and right of the phase angle.

In one embodiment, the acquiring unit is further configured to: if the phase angle corresponding to the minimum error value is not the primary phase angle for the first test, determine whether the phase angle corresponding to the minimum error value is the last The final phase angle for the test at a time; if the phase angle corresponding to the minimum error value is the final phase angle for the last test, obtain the primary phase angle for the first test and sum the primary phase angle The third phase angle for the penultimate test is used as the two phase angles adjacent to the left and right of the phase angle.

For the steps implemented by the functional modules of the carrier phase estimation device, reference may be made to the embodiments of the carrier phase estimation method of the present disclosure, which will not be repeated here.

The present disclosure also provides a terminal including: a memory, a processor, a communication bus, and a carrier phase estimation program stored on the memory: the communication bus is used to implement connection communication between the processor and the memory; The processor is used to execute the carrier phase estimation program to implement the steps of the above embodiments of the carrier phase estimation method.

The present disclosure also provides a computer-readable storage medium that stores one or more programs, and the one or more programs may also be executed by one or more processors for implementation The steps of the above embodiments of the carrier phase estimation method.

The implementation of the computer-readable storage medium of the present disclosure is basically the same as the above embodiments of the carrier phase estimation method, and details are not described herein again.

It should be noted that in this article, the terms "include", "include" or any other variant thereof are intended to cover non-exclusive inclusion, so that a process, method, article or system that includes a series of elements includes not only those elements, It also includes other elements that are not explicitly listed, or include elements inherent to this process, method, article, or system. Without more restrictions, the element defined by the sentence "include a..." does not exclude that there are other identical elements in the process, method, article or system that includes the element.

The sequence numbers of the above-mentioned embodiments of the present disclosure are for description only, and do not represent the merits of the embodiments.

Through the description of the above embodiments, those skilled in the art can clearly understand that the methods in the above embodiments can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware, but in many cases the former is better Implementation. Based on such an understanding, the technical solution of the present disclosure can be embodied in the form of a software product in essence or part that contributes to the existing technology, and the computer software product is stored in a storage medium (such as ROM/RAM) as described above , Magnetic disks, optical disks), including several instructions to enable a terminal device (which may be a mobile phone, computer, server, air conditioner, or network device, etc.) to perform the methods described in various embodiments of the present disclosure.

The present disclosure achieves the technical effect of greatly reducing the complexity of the algorithm without reducing the accuracy of the phase noise compensation, which is beneficial to hardware implementation. The present disclosure obtains a phase noise signal with only phase noise by performing digital signal processing on the signal under test acquired in the current symbol; rotating the phase noise signal based on each preset test angle to obtain each rotation signal, and based on The error function obtains the error value corresponding to each of the rotation signals; obtains the minimum error value among each of the error values, and determines the phase angle corresponding to the minimum error value and the two phase angles adjacent to the left and right of the phase angle Determine the optimal phase angle based on the minimum error value and the phase angle corresponding to the minimum error value and the two phase angles adjacent to the left and right of the phase angle; obtain each estimated phase angle before the current symbol, And determine the target phase angle based on the optimal phase angle and each of the estimated phase angles. By using the error function to calculate the corresponding error value for the obtained signal under test, compared with the traditional calculation method using a large number of multipliers, the number of test angles is reduced, and the calculation complexity of the error function is greatly reduced , And in determining the target phase angle, that is, the estimated phase noise, in addition to estimating the current optimal phase angle, it is also related to the previously obtained estimated phase angle, thereby effectively reducing the probability of the occurrence of phase cycle slips and improving the system. The performance has achieved the technical effect of greatly reducing the complexity of the algorithm without reducing the accuracy of the phase noise compensation, which is conducive to hardware implementation.

The above are only preferred embodiments of the present disclosure and do not limit the patent scope of the present disclosure. Any equivalent structure or equivalent process transformation made by using the contents of the specification and drawings of the present disclosure, or directly or indirectly used in other related technical fields , The same reason is included in the scope of patent protection of this disclosure.

Claims

A carrier phase estimation method, wherein the carrier phase estimation method includes the following steps:

Perform digital signal processing on the signal under test acquired in the current symbol to obtain a phase noise signal with only phase noise;

Rotating the phase noise signal based on each preset test angle to obtain each rotation signal, and obtaining an error value corresponding to each rotation signal based on an error function;

Acquiring a minimum error value among each of the error values, and determining a phase angle corresponding to the minimum error value and two phase angles adjacent to the left and right of the phase angle;

Determine an optimal phase angle based on the minimum error value and the phase angle corresponding to the minimum error value and two phase angles adjacent to the left and right of the phase angle;

Obtain each estimated phase angle before the current symbol, and determine a target phase angle based on the optimal phase angle and each of the estimated phase angles.
The carrier phase estimation method according to claim 1, wherein the step of acquiring each estimated phase angle before the current symbol and determining a target phase angle based on the optimal phase angle and each of the estimated phase angles ,include:

Acquiring each historical symbol before the current symbol, and acquiring each estimated phase angle corresponding to each historical symbol;

Obtain the average value of the phase angle between the optimal phase angle and each of the estimated phase angles, and determine the target phase angle based on the optimal phase angle and the average phase angle.
The carrier phase estimation method according to claim 2, wherein the step of determining a target phase angle based on the optimal phase angle and the average value of the phase angles includes:

Acquiring the maximum preset test angle and the minimum preset test angle among each of the preset test angles, and comparing the optimal phase angle with the average value of the phase angles;

If the optimal phase angle is greater than the first sum between the maximum preset test angle and the average value of the phase angles, obtain a first difference between the optimal phase angle and the preset angle, And use the first difference as the target phase angle;

If the optimal phase angle is less than the second sum value between the average phase angle and the minimum preset test angle, a third sum between the optimal phase angle and the preset angle is obtained Value, and use the third sum value as the target phase angle.
The carrier phase estimation method according to claim 1, wherein the step of obtaining an error value corresponding to each of the rotation signals based on an error function includes:

Acquiring the modulation format corresponding to the rotation signal, and determining the mapping value in the modulation format;

Acquiring the real part value and the imaginary part value in the rotation signal, and calculating the real part difference between the real part value and the mapped value according to an error function, between the imaginary part value and the mapped value The difference of the imaginary part of

Acquiring an error sum value between the absolute real part difference corresponding to the real part difference and the absolute imaginary part difference corresponding to the imaginary part difference, and determining the error corresponding to the rotation signal based on the error sum value.
The carrier phase estimation method according to claim 4, wherein the step of determining the error value corresponding to the rotation signal based on the error sum value includes:

Acquiring a primary error value corresponding to the error sum value and a preset number of intermediate error values before the current symbol;

Acquiring an average value of errors between the primary error value and each of the intermediate error values, and using the average error value as an error value corresponding to the rotation signal.
The carrier phase estimation method according to claim 1, wherein the minimum error value of each of the error values is obtained, and a phase angle corresponding to the minimum error value and two adjacent left and right sides of the phase angle are determined The phase angle steps include:

Acquiring the minimum error value among each of the error values, and determining the phase angle corresponding to the minimum error value;

Judging in the current symbol whether the phase angle corresponding to the minimum error value is the primary phase angle for the first test;

If the phase angle corresponding to the minimum error value is the primary phase angle to be tested first, the final phase angle of the last test in the current symbol is obtained, and the second test in the current symbol is tested The second phase angle and the final phase angle serve as two phase angles adjacent to the left and right of the phase angle.
The carrier phase estimation method according to claim 6, wherein the step of determining whether the phase angle corresponding to the minimum error value is the primary phase angle for the first test includes:

If the phase angle corresponding to the minimum error value is not the primary phase angle for the first test, determine whether the phase angle corresponding to the minimum error value is the final phase angle for the last test;

If the phase angle corresponding to the minimum error value is the final phase angle of the last test, the primary phase angle of the first test is obtained, and the primary phase angle and the penultimate time of the second test are obtained. The three-phase angle serves as two phase angles adjacent to the left and right of the phase angle.
A carrier phase estimation device, wherein the carrier phase estimation device includes:

The signal processing unit is used to perform digital signal processing on the signal under test acquired in the current symbol to obtain a phase noise signal with only phase noise;

A rotation unit, configured to rotate the phase noise signal based on each preset test angle to obtain each rotation signal, and obtain an error value corresponding to each rotation signal based on an error function;

An obtaining unit, configured to obtain a minimum error value among each of the error values, and determine a phase angle corresponding to the minimum error value and two phase angles adjacent to the left and right of the phase angle;

A determining unit, configured to determine an optimal phase angle based on the minimum error value and a phase angle corresponding to the minimum error value and two phase angles adjacent to the left and right of the phase angle;

The phase angle unit is used to obtain each estimated phase angle before the current symbol, and determine a target phase angle based on the optimal phase angle and each estimated phase angle.
A carrier phase estimation apparatus, wherein the carrier phase estimation apparatus includes: a memory, a processor, and a carrier phase estimation program stored on the memory and operable on the processor, the carrier phase estimation program is The processor executes the steps of the carrier phase estimation method according to any one of claims 1 to 7 when executed.
A computer-readable storage medium, wherein a carrier-phase estimation program is stored on the computer-readable storage medium, and when the carrier-phase estimation program is executed by a processor, it is implemented according to any one of claims 1 to 7. Carrier phase estimation method steps.