WO2022022521A1

WO2022022521A1 - Phase noise estimation method and apparatus

Info

Publication number: WO2022022521A1
Application number: PCT/CN2021/108706
Authority: WO
Inventors: 李峰; 杜恒; 薛江; 孟德宇; 段启宏; 每媛
Original assignee: 华为技术有限公司
Priority date: 2020-07-31
Filing date: 2021-07-27
Publication date: 2022-02-03
Also published as: CN114070682A; CN114070682B

Abstract

Provided are a phase noise estimation method and apparatus. The method comprises: firstly, receiving an analog waveform signal which comes from a sending end and has been affected by phase noise; performing demodulation processing on the analog waveform signal to obtain a time-frequency signal corresponding to the analog waveform signal; and after the time-frequency signal corresponding to the analog waveform signal is obtained, constructing a cyclic matrix according to a phase noise estimation order P and the number of sub-carriers, such that the cyclic matrix can be used to estimate inter-carrier interference caused by the phase noise, thereby estimating the inter-carrier interference caused by the phase noise. In addition, since the elements of a cyclic matrix comprise a phase tracking reference signal (PTRS) and user data, when the cyclic matrix is used to estimate inter-carrier interference caused by phase noise, the accuracy of inter-carrier interference estimation can be improved to a certain extent.

Description

Phase noise estimation method and apparatus

This application claims the priority of the Chinese patent application with the application number 202010761910.7 and the application title "Phase Noise Estimation Method and Apparatus" filed with the China Patent Office on July 31, 2020, the entire contents of which are incorporated into this application by reference.

technical field

The present application relates to the field of communication technologies, and in particular, to a phase noise estimation method and apparatus.

Background technique

In the process of signal transmission, due to the influence of frequency processing equipment, such as modulators, local oscillators, modulators, etc., and other factors, phase noise will be generated on the transmission system. The performance of phase noise is to form random phase disturbances to stable sine wave signals. , reflected in the time domain is the random jitter of the signal, that is, clock jitter; reflected in the frequency domain is the reduction of spectral purity, that is, spectrum broadening.

Especially in orthogonal frequency division multiplexing (OFDM) systems, the impact of phase noise is mainly reflected in two aspects: one is that phase noise will cause common phase errors, and the other is that phase noise will cause inter-carrier interference (Inter-Carrier Interference). -Carrier Interference, ICI). Since the common phase error caused by phase noise has the same effect on each subcarrier in the same OFDM symbol, the phase tracking pilot (Phase Tracking Reference Signal, PTRS) can be used to track and suppress the common phase error caused by phase noise. Thus, the estimation and compensation of the common phase error are realized.

However, since the interference between sub-carriers is random and the ICIs of different sub-carriers are different, it is currently impossible to estimate the inter-carrier interference caused by phase noise.

SUMMARY OF THE INVENTION

The phase noise estimation method and apparatus provided by the embodiments of the present application realize the estimation of the inter-carrier interference caused by the phase noise.

In a first aspect, an embodiment of the present application provides a phase noise estimation method, and the phase noise estimation method may include:

Receives phase noise-affected analog waveform signals from the transmitter.

The analog waveform signal is demodulated to obtain the time-frequency signal corresponding to the analog waveform signal.

Determine the phase noise estimation order P according to the time-frequency signal; P is an integer greater than or equal to 1 and less than or equal to num, where num is used to represent the number of subcarriers used when receiving the analog waveform signal.

A circulant matrix is constructed according to the phase noise estimation order P and the number of subcarriers; wherein, the elements of the circulant matrix include the phase tracking pilot PTRS and user data.

The inter-carrier interference caused by phase noise is estimated using a circulant matrix.

It can be seen that in the embodiment of the present application, the receiving end first receives the analog waveform signal affected by the phase noise from the transmitting end; demodulates the analog waveform signal to obtain the time-frequency signal corresponding to the analog waveform signal, and then obtains the analog waveform signal after obtaining the analog waveform signal. After the time-frequency signal corresponding to the waveform signal, a cyclic matrix is constructed by estimating the order P of the phase noise and the number of sub-carriers, so that the cyclic matrix can be used to estimate the inter-carrier interference caused by the phase noise, so as to realize the carrier interference caused by the phase noise. Interference is estimated. In addition, since the elements of the circulant matrix include the phase tracking pilot PTRS and user data, when the circulant matrix is used to estimate the ICI caused by the phase noise, the accuracy of the ICI estimation can be improved to a certain extent.

In a possible implementation manner, using a cyclic matrix to estimate the inter-carrier interference caused by phase noise may include:

The circulant matrix is separated, and the real part circulant matrix and the imaginary part circulant matrix corresponding to the circulant matrix are obtained.

Taking the real circulant matrix and the imaginary circulant matrix as matrix sub-blocks, the block circulant matrix corresponding to the circulant matrix is constructed.

Inter-carrier interference caused by phase noise is estimated using a block circulant matrix.

It can be understood that, in the embodiment of the present application, by separating the real part and the imaginary part corresponding to the circulant matrix, the advantage is that the number and precision of row selection can be increased, which can better estimate the inter-carrier interference.

In a possible implementation manner, the circulant matrix is a matrix with num rows and 2P+1 columns, and the block circulant matrix is a matrix with 2num rows and 4P+2 columns.

Among them, the first num rows and the first 2P+1 columns of the block circulant matrix are the real part circulant matrix, the first num rows and the last 2P+1 columns of the block circulant matrix are the imaginary part circulant matrix, The last num rows and first 2P+1 columns of the block cyclic matrix are imaginary cyclic matrices, and the last num rows and last 2P+1 columns of the block cyclic matrix are real cyclic matrices.

In a possible implementation manner, using a block cyclic matrix to estimate the inter-carrier interference caused by phase noise may include:

Screen the row vectors in the block circulant matrix to obtain the target circulant matrix.

The inter-carrier interference caused by phase noise is estimated using the target circulant matrix.

It can be understood that, in the embodiment of the present application, by screening the row vectors in the block circulant matrix, the advantage is that the accuracy of the ICI estimation can be improved to a certain extent.

In a possible implementation manner, the row vectors in the block circulant matrix are screened to obtain the target circulant matrix, which may include:

Determine whether the k-th element in the P-th column of the block cyclic matrix is a PTRS; k is any value of an integer greater than or equal to 1 and less than or equal to 2num.

If the k-th element is PTRS, the k-th row element is reserved; if the k-th element is non-PTRS, it is determined whether the value of the k-th element is greater than the preset threshold.

If the value of the kth element is greater than the preset threshold, keep the kth row element; if the value of the kth element is less than or equal to the preset threshold, delete the kth row element in the block circular matrix to obtain the target circular matrix .

In a possible implementation manner, the target cyclic matrix is used to estimate the inter-carrier interference caused by the phase noise, which may include:

Construct a column vector, the column vector is a num-dimensional column vector, and each element in the column vector is the judgment estimate of the time-frequency signal on the constellation point.

The column vector is separated to obtain the real part column vector and the imaginary part column vector corresponding to the column vector.

Taking the real part column vector and the imaginary part column vector as matrix sub-blocks, construct the block column vector corresponding to the column vector.

The ICI due to phase noise is estimated using the target circulant matrix and the block column vector.

It can be seen that in the embodiment of the present application, by constructing a column vector and constructing a target column vector based on the column vector, the target cyclic matrix and the target column vector can be used to estimate the inter-carrier interference caused by the phase noise, which not only realizes the Estimating the inter-carrier interference caused by phase noise; and can improve the accuracy of inter-carrier interference estimation to a certain extent. In addition, in the embodiment of the present application, by separating the real part and the imaginary part corresponding to the column vector, the advantage is that the number and precision of row selection can be increased, which can better estimate the inter-carrier interference.

In one possible implementation, the block column vector is a 2num-dimensional column vector.

The matrix sub-blocks of the first num rows of the block column vector are real part column vectors, and the matrix sub-blocks of the last num rows of the block column vector are imaginary part column vectors.

In a possible implementation manner, the target cyclic matrix and the block column vector are used to estimate the inter-carrier interference caused by the phase noise, which may include:

According to the target circulant matrix, the row vector in the block column vector is screened to obtain the target column vector; The number of rows occupied by the vector in the block column vector is the same; wherein, s=t, and s is an integer greater than or equal to 1 and less than or equal to 2num.

The ICI due to phase noise is estimated using the target circulant matrix and the target column vector.

It can be understood that, in the embodiment of the present application, by screening the row vector in the block column vector, the advantage is that the accuracy of the ICI estimation can be improved to a certain extent.

In a possible implementation manner, the target circulant matrix and the target column vector are used to estimate the inter-carrier interference caused by the phase noise, which may include:

The phase noise tap at the first moment is determined according to I=(Q ^T Q) ^-1 Q ^T U; wherein, I represents the phase noise tap, Q represents the target circulant matrix, and U represents the target column vector.

The ICI caused by the phase noise is estimated according to the phase noise tap, thereby realizing the estimation of the ICI caused by the phase noise.

In a possible implementation manner, the phase noise estimation method may further include:

A convolution operation is performed on the preset sequence and the phase noise taps to obtain a result of the convolution operation; wherein the preset sequence is a sequence of all subcarriers used at the first moment when the time-frequency signal is received.

Compensating the time-frequency signal according to the result of the convolution operation not only realizes the estimation of the inter-carrier interference caused by the phase noise, but also further realizes the signal compensation for the received time-frequency signal.

In a possible implementation, the circulant matrix is a matrix with num rows and 2P+1 columns.

Among them, the first P elements of the first column of the cyclic matrix are all 0, and the P+1th element to the numth element are the first num-P elements of the preset sequence in sequence; the first P elements of the second column of the cyclic matrix are -1 elements are all 0, the P-th element to the num-th element are the first num-P+1 elements of the preset sequence, and so on, until the first element of the second column of the cyclic matrix is 0, the second element to the num-th element are the first num-1 elements of the preset sequence.

The first element to the numth element of the P+1th of the circulant matrix are sequentially the elements in the preset sequence.

The first num-1 elements of the P+2th column of the circulant matrix are sequentially the second to numth elements of the preset sequence, and the numth element is 0; the first num-1 elements of the P+3th column of the circulant matrix The 2 elements are sequentially the 3rd element to the numth element of the preset sequence, the num-1th element and the numth element are both 0, and so on, until the first num of the 2P+1 column of the circular matrix The -P elements are sequentially the num-P+1-th element to the num-th element of the preset sequence, and the P-th element to the num-th element are all 0.

The preset sequence is R(i, 0), R(i, 1), R(i, 2), R(i, 3), ..., R(i, num-2), R(i, num -1).

In a second aspect, an embodiment of the present application provides a phase noise estimation apparatus, and the phase noise estimation apparatus may include:

The receiving unit is used for receiving the analog waveform signal affected by the phase noise from the transmitting end.

The processing unit is used to demodulate the analog waveform signal to obtain the time-frequency signal corresponding to the analog waveform signal; and determine the phase noise estimation order P according to the time-frequency signal; P is greater than or equal to 1 and less than or equal to num Integer, num is used to indicate the number of subcarriers used when receiving the analog waveform signal.

The processing unit is further configured to construct a circulant matrix according to the phase noise estimation order P and the number of subcarriers; wherein, the elements of the circulant matrix include the phase tracking pilot PTRS and user data.

The estimation unit is used for estimating the inter-carrier interference caused by the phase noise by using the cyclic matrix.

In a possible implementation manner, the processing unit is further configured to perform separation processing on the circulant matrix to obtain a real part circulant matrix and an imaginary part circulant matrix corresponding to the circulant matrix; and use the real part circulant matrix and the imaginary part circulant matrix as matrices Subblock, construct the block circulant matrix corresponding to the circulant matrix.

The estimation unit is specifically used for estimating the inter-carrier interference caused by the phase noise by using the block cyclic matrix.

In a possible implementation manner, the processing unit is further configured to perform screening processing on the row vectors in the block circulant matrix to obtain the target circulant matrix.

The estimation unit is specifically used for estimating the inter-carrier interference caused by the phase noise by using the target cyclic matrix.

In a possible implementation manner, the processing unit is specifically configured to determine whether the kth element in the Pth column of the block cyclic matrix is a PTRS; k is an integer greater than or equal to 1 and less than or equal to 2num Any value of ; if the kth element is PTRS, then keep the kth row element; if the kth element is non-PTRS, then judge whether the value of the kth element is greater than the preset threshold; if the value of the kth element is greater than the preset threshold; If the value is greater than the preset threshold, the element in the kth row is retained; if the value of the kth element is less than or equal to the preset threshold, the element in the kth row is deleted from the block cyclic matrix to obtain the target cyclic matrix.

In a possible implementation manner, the processing unit is further configured to construct a column vector, the column vector is a column vector of num dimensions, and each element in the column vector is a judgment estimate of the constellation point by the time-frequency signal; for the column vector Perform separation processing to obtain the real and imaginary column vectors corresponding to the column vectors; and use the real and imaginary column vectors as matrix sub-blocks to construct block column vectors corresponding to the column vectors.

The estimation unit is specifically used for estimating the inter-carrier interference caused by the phase noise by using the target cyclic matrix and the block column vector.

In a possible implementation manner, the processing unit is further configured to screen the row vectors in the block column vector according to the target circulant matrix to obtain the target column vector; The number of rows occupied by the matrix is the same as the number of rows occupied by the t-th row vector in the target column vector in the block column vector; where s=t, and s is an integer greater than or equal to 1 and less than or equal to 2num.

The estimation unit is specifically used for estimating the inter-carrier interference caused by the phase noise by using the target cyclic matrix and the target column vector.

In a possible implementation manner, the estimation unit is specifically configured to determine the phase noise tap at the first moment according to I=(Q ^T Q) ^-1 Q ^T U; and according to the phase noise tap, the inter-carrier interference caused by the phase noise is estimate; where I represents the phase noise tap, Q represents the target circulant matrix, and U represents the target column vector.

In a possible implementation manner, the processing unit is further configured to perform a convolution operation on the preset sequence and the phase noise taps to obtain a result of the convolution operation; and perform signal compensation on the time-frequency signal according to the result of the convolution operation; wherein, The preset sequence is the sequence of all subcarriers used at the first moment when the time-frequency signal is received.

In a third aspect, an embodiment of the present application further provides a communication apparatus, the apparatus may include a processor and a memory, where a computer program is stored in the memory, and the processor executes the computer program stored in the memory, to The apparatus is caused to perform the phase noise estimation method described in any possible implementation manner of the first aspect above.

In a fourth aspect, an embodiment of the present application further provides a communication apparatus, where the communication apparatus may include: a processor and an interface circuit.

The interface circuit is used to receive code instructions and transmit them to the processor.

The processor is configured to run the code instructions to execute the phase noise estimation method described in any possible implementation manner of the first aspect.

In a fifth aspect, an embodiment of the present application further provides a chip, where a computer program is stored on the chip, and when the computer program is executed by a processor, the phase noise estimation described in any possible implementation manner of the first aspect is performed. method.

In a sixth aspect, an embodiment of the present application further provides a computer-readable storage medium, where an instruction is stored in the computer-readable storage medium, and when the instruction is executed on a communication device, the communication device is made to execute the above-mentioned first step. In one aspect, the phase noise estimation method described in any possible implementation manner.

It can be seen that the phase noise estimation method and device provided by the embodiments of the present application first receive the analog waveform signal from the transmitting end that has been affected by the phase noise, and then demodulate the analog waveform signal to obtain the time-frequency signal corresponding to the analog waveform signal. , after the time-frequency signal corresponding to the analog waveform signal is obtained, a cyclic matrix is constructed by estimating the order P of the phase noise and the number of sub-carriers, so that the cyclic matrix can be used to estimate the inter-carrier interference caused by the phase noise, so as to realize the phase noise. Inter-carrier interference caused by noise is estimated. In addition, since the elements of the circulant matrix include the phase tracking pilot PTRS and user data, when the circulant matrix is used to estimate the ICI caused by the phase noise, the accuracy of the ICI estimation can be improved to a certain extent.

Description of drawings

FIG. 1 is a schematic diagram of an application scenario provided by an embodiment of the present application;

FIG. 2 is a schematic flowchart of a phase noise estimation method provided by an embodiment of the present application;

FIG. 3 is a schematic structural diagram of a time-frequency signal provided by an embodiment of the present application;

FIG. 4 is a schematic flowchart of another phase noise estimation method provided by an embodiment of the present application;

FIG. 5 is a schematic structural diagram of a phase noise estimation apparatus provided by an embodiment of the present application;

FIG. 6 is a schematic structural diagram of a communication device according to an embodiment of the present application.

detailed description

The phase noise estimation method provided in the embodiments of the present application can be applied to quadrature modulation and demodulation, or a communication system similar to quadrature modulation and demodulation. In a quadrature modulation and demodulation, or a communication system similar to quadrature modulation and demodulation, the phase noise estimation method provided by the embodiment of the present application can be used to estimate the inter-carrier interference caused by the phase noise.

For example, please refer to FIG. 1 , which is a schematic diagram of an application scenario provided by an embodiment of the present application. Taking an orthogonal modulation and demodulation communication system as an example, the system includes a transmitter and a receiver. When the transmitter sends a time-frequency signal, its modulation module will first convert the time-frequency signal to be sent into an analog waveform signal, and then send the analog waveform to the receiver. However, when sending, some devices in the transmitter, such as a modulator , mixers, and frequency dividers will introduce phase noise into the analog waveform signal. In addition, due to other factors, additive Gaussian phase noise will also be introduced, so that the analog waveform signal received at the receiving end is phase noise. For the analog waveform signal affected by noise, after receiving the analog waveform signal affected by the phase noise, the receiving end demodulates the analog waveform signal, and obtains the time-frequency signal after demodulation processing. Noise-affected time-frequency signals. Phase noise will not only cause common phase error, but also cause ICI. In order to eliminate the common phase error and ICI caused by phase noise in time-frequency signals, it is necessary to estimate the common phase error and ICI. In order to realize the compensation of the time-frequency signal. When estimating the common phase error, the phase tracking pilot (Phase Tracking Reference Signal, PTRS) can be used to track and suppress the common phase error caused by the phase noise, so as to realize the estimation of the common phase error; however, due to the sub-carriers The inter-carrier interference is random, and the ICIs of different sub-carriers are different. Therefore, it is currently impossible to estimate the inter-carrier interference caused by phase noise.

In order to estimate the inter-carrier interference caused by phase noise, an embodiment of the present application provides a method for estimating phase noise. The time-frequency signal corresponding to the analog waveform signal is obtained; the phase noise estimation order P is determined according to the time-frequency signal; the circulant matrix is constructed according to the phase noise estimation order P and the number of subcarriers; wherein, the elements of the circulant matrix include the phase Tracks pilot PTRS and user data; and uses a cyclic matrix to estimate inter-carrier interference caused by phase noise. Wherein, P is an integer greater than or equal to 1 and less than or equal to num, where num is used to represent the number of subcarriers used when receiving the analog waveform signal.

It can be seen that, in the embodiment of the present application, by determining the phase noise estimation order P according to the time-frequency signal, and constructing a cyclic matrix according to the phase noise estimation order P and the number of subcarriers, the cyclic matrix can be used to cause phase noise. to estimate the inter-carrier interference caused by phase noise. In addition, since the elements of the circulant matrix include the phase tracking pilot PTRS and user data, when the circulant matrix is used to estimate the ICI caused by the phase noise, the accuracy of the ICI estimation can be improved to a certain extent.

For example, the receiving end may be a terminal device or a network device. The terminal equipment, also known as user equipment (UE), mobile station (MS), mobile terminal (MT), etc., is a device that provides voice/data connectivity to users. , for example, handheld devices with wireless connectivity, in-vehicle devices, etc. At present, some examples of terminal devices are: mobile phone (mobile phone), tablet computer, notebook computer, PDA, mobile internet device (MID), wearable device, virtual reality (VR) device, augmented Augmented reality (AR) equipment, wireless terminals in industrial control, wireless terminals in self-driving, wireless terminals in remote medical surgery, smart grid wireless terminals in transportation safety, wireless terminals in smart cities, wireless terminals in smart homes, etc. The method executed by the terminal device in the embodiment of the present application may be specifically executed by at least one chip in the terminal device.

A network device is an entity on the network side that transmits or receives signals, such as a new generation base station (generation Node B, gNodeB). A network device may be a device for communicating with terminal devices. The network equipment can be an evolved base station (evolutional Node B, eNB or eNodeB) in Long Term Evolution (LTE), or a relay station or an access point, or a vehicle-mounted device, a wearable device, and a network device in the future 5G network Or the network equipment in the future evolved public land mobile network (Public Land Mobile Network, PLMN) network, or the gNodeB in the NR system, etc. In addition, in this embodiment of the present application, a network device provides services for a cell, and a terminal device communicates with the network device through transmission resources (for example, frequency domain resources, or spectrum resources) used by the cell, and the cell may be a network device A cell corresponding to a cell (such as a base station), the cell may belong to a macro base station, or it may belong to a base station corresponding to a small cell. The small cell here may include: urban cell (metro cell), micro cell (micro cell), pico cell (pico cell), femto cell (femto cell), etc. These small cells have the characteristics of small coverage and low transmission power, and are suitable for providing high-speed data transmission services. In addition, in other possible cases, the network device may be other devices that provide wireless communication functions for the terminal device. The embodiments of the present application do not limit the specific technology and specific device form adopted by the network device. For convenience of description, in this embodiment of the present application, a device that provides a wireless communication function for a terminal device is referred to as a network device. The method performed by the network device in the embodiment of the present application may be specifically performed by at least one chip in the network device.

In the embodiments of the present application, "at least one" refers to one or more, and "a plurality" refers to two or more. "And/or", which describes the association relationship of the associated objects, means that there can be three kinds of relationships, for example, A and/or B, which can mean that A exists alone, A and B exist at the same time, and B exists alone, where A and B can be singular or plural. In the text description of this application, the character "/" generally indicates that the contextual object is an "or" relationship.

Hereinafter, the phase noise estimation method provided by the embodiments of the present application will be described in detail through the following embodiments. It can be understood that the following embodiments may be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments.

FIG. 2 is a schematic flowchart of a phase noise estimation method provided by an embodiment of the present application. The phase noise estimation method may be executed by software and/or a hardware device. For example, the hardware device may be a phase noise estimation device, and the phase noise estimation The apparatus may be provided in an electronic device. For example, referring to FIG. 2 , the phase noise estimation method provided by the embodiment of the present application may include:

S201. Receive an analog waveform signal affected by phase noise from a transmitting end, and perform demodulation processing on the analog waveform signal to obtain a time-frequency signal corresponding to the analog waveform signal.

When the sending end sends time-frequency signals, its modulation module will first convert the time-frequency signals to be sent into analog waveform signals, and then send the converted analog waveform signals to the receiving end. Introduce phase noise, so that the analog waveform signal received by the receiving end is an analog waveform signal affected by phase noise; after receiving the analog waveform signal affected by phase noise, first demodulate the analog waveform signal affected by phase noise process to obtain a time-frequency signal after demodulation processing, the time-frequency signal is a two-dimensional time-frequency signal, and the structure of the two-dimensional time-frequency signal can be referred to as shown in FIG. A schematic diagram of the structure of a time-frequency signal. It can be seen that the horizontal direction of the time-frequency signal is the time dimension, the width of each grid in the horizontal direction is the duration of a symbol, the vertical direction is the frequency dimension, and the vertical direction is the frequency dimension. The width of each grid in the direction is one subcarrier interval, and each grid in FIG. 3 is defined as one resource element RE.

Exemplarily, if the time-frequency signal is received by the receiver on the j-th subcarrier of the i-th symbol, the time-frequency signal can be defined as R(i, j). When estimating the inter-carrier interference caused by the included phase noise, the phase noise estimation order P may be determined first according to the time-frequency signal, that is, the following S202 is performed:

S202. Determine the phase noise estimation order P according to the time-frequency signal.

Wherein, P is an integer greater than or equal to 1 and less than or equal to num, where num is used to represent the number of subcarriers used when receiving the analog waveform signal. For example, when determining the number of subcarriers used when receiving the analog waveform signal, if the receiving end uses subcarriers on one resource block to receive the time-frequency signal, the number of subcarriers is 12; if the receiving end uses two subcarriers When the subcarriers on the resource block receive time-frequency signals, the number of subcarriers is 24, which can be set according to actual needs. Here, the number of subcarriers is not further limited in this embodiment of the present application.

Exemplarily, when the phase estimation order P is determined according to the time-frequency signal, the selection of P is related to the characteristics of the phase noise of the physical device. Generally, the larger the phase noise width of the physical device is, the larger the value of the corresponding phase estimation order P is; on the contrary, the smaller the phase noise width of the physical device is, the higher the corresponding phase estimation order P is. It can be set according to actual needs. Here, the value of the phase noise estimation order P is not further limited in this embodiment of the present application.

After the phase noise estimation order P is determined, a circulant matrix can be constructed according to the phase noise estimation order P and the number of subcarriers, so that the circulant matrix can be used to estimate the inter-carrier interference caused by the phase noise, that is, the following S203-S204 are performed. :

S203. Construct a circulant matrix according to the phase noise estimation order P and the number of subcarriers.

The elements of the circulant matrix include phase tracking pilot PTRS and user data. The user data here can be understood as the non-phase tracking pilots in the transmission process, and the user data can be all the non-phase tracking pilots in the transmission process, or part of the non-phase tracking pilots in the transmission process. Can be set according to actual needs.

Exemplarily, when constructing a circulant matrix according to the phase noise estimation order P and the number of subcarriers, the constructed circulant matrix may be a matrix with num rows and 2P+1 columns; wherein, the first P elements of the first column of the circulant matrix are all is 0, the P+1th element to the numth element are the first num-P elements of the preset sequence in sequence; the first P-1 elements of the second column of the cyclic matrix are all 0, and the Pth element to The num-th element is the first num-P+1 elements of the preset sequence, and so on, until the first element of the second column of the cyclic matrix is 0, and the second element to the num-th element are the preset sequence. Let the first num-1 elements of the sequence; the first element to the num-th element of the P+1th of the circulant matrix are the elements in the preset sequence in turn; the first num-1 elements of the P+2th column of the circulant matrix The elements are the 2nd element to the numth element of the preset sequence, and the numth element is 0; the first num-2 elements of the P+3th column of the cyclic matrix are the 3rd element to the numth element of the preset sequence. The num-th element, the num-1-th element and the num-th element are all 0, and so on, until the first num-P elements of the 2P+1 column of the circular matrix are in turn the num-Pth elements of the preset sequence +1 element to num-th element, P-th element to num-th element are all 0; the preset sequence is R(i, 0), R(i, 1), R(i, 2), R(i, 3), ..., R(i, num-2), R(i, num-1).

Combined with the text description about the circulant matrix in S203, the circulant matrix is converted into a mathematical formula, and the following formula 1 can be referred to:

Among them, the matrix M is the constructed circulant matrix.

After the circulant matrix is obtained by constructing, the constructed circulant matrix can be used to estimate the inter-carrier interference caused by the phase noise, that is, the following S204 is performed:

S204, using a cyclic matrix to estimate the inter-carrier interference caused by the phase noise.

For example, when using a circulant matrix to estimate the inter-carrier interference caused by phase noise, since the signal is usually a complex signal, the circulant matrix can be separated first to obtain the real part circulant matrix and the imaginary part circulant matrix corresponding to the circulant matrix. The real circulant matrix and the imaginary circulant matrix are used as matrix sub-blocks to construct the block circulant matrix corresponding to the circulant matrix; and the block circulant matrix is used to estimate the inter-carrier interference caused by the phase noise. The real circulant matrix corresponding to the circulant matrix may be represented by Real(M), and the imaginary circulant matrix corresponding to the circulant matrix may be represented by Imag(M). It can be understood that, in the embodiment of the present application, by separating the real part and the imaginary part corresponding to the circulant matrix, the advantage is that the number and precision of row selection can be increased, which can better estimate the inter-carrier interference.

Exemplarily, when constructing a block circulant matrix corresponding to a circulant matrix, in combination with the description in S203 above, when the circulant matrix is a matrix with num rows and 2P+1 columns, the corresponding block circulant matrix may be 2num rows, 4P+ 2-column matrix. Among them, the first num rows and the first 2P+1 columns of the block circulant matrix are the real part circulant matrix, the first num rows and the last 2P+1 columns of the block circulant matrix are the imaginary part circulant matrix, The last num rows and first 2P+1 columns of the block cyclic matrix are imaginary cyclic matrices, and the last num rows and last 2P+1 columns of the block cyclic matrix are real cyclic matrices.

Combined with the text description about the block cyclic matrix in S204, the block cyclic matrix is converted into a mathematical formula, and the following formula 2 can be referred to:

in,

Represents the constructed block circulant matrix, Real(M) represents the real part circulant matrix corresponding to the circulant matrix, and Imag(M) represents the imaginary part circulant matrix corresponding to the circulant matrix.

After the real circulant matrix and the imaginary circulant matrix are used as matrix sub-blocks to construct the block circulant matrix shown in the above formula 2, the row vectors in the block circulant matrix can be screened to obtain the target circulant matrix; The inter-carrier interference caused by phase noise is estimated using the target circulant matrix. It can be understood that, in the embodiment of the present application, by screening the row vectors in the block circulant matrix, the advantage is that the accuracy of the ICI estimation can be improved to a certain extent.

For example, when screening the row vectors in the block circulant matrix, it can be first judged whether the kth element in the Pth column in the block circulant matrix is PTRS; k is greater than or equal to 1, and less than or Any value of an integer equal to 2num; that is, each element in the P-th column in the block cyclic matrix is judged in turn, if the k-th element is PTRS, the k-th row element is reserved; if the k-th element is If it is non-PTRS, then judge whether the value of the kth element is greater than the preset threshold; if the value of the kth element is greater than the preset threshold, keep the kth row element; if the value of the kth element is less than or equal to the preset threshold , then delete the k-th row element in the block circulant matrix to obtain the target circulant matrix.

The value of the preset threshold is related to the modulation method. Exemplarily, the value of the preset threshold is related to the real part and the imaginary part of the outermost constellation point of the constellation point in the modulation mode, and can be set according to actual needs. The application examples do not make specific limitations.

It can be seen that the phase noise estimation method provided by the embodiment of the present application determines the phase noise estimation order P according to the received time-frequency signal, and constructs a cyclic matrix according to the phase noise estimation order P and the number of subcarriers, and based on the cyclic The matrix is constructed to obtain the target cyclic matrix, so that the target cyclic matrix can be used to estimate the ICI caused by the phase noise, so as to realize the estimation of the ICI caused by the phase noise. In addition, since the elements of the target cyclic matrix include the phase tracking pilot PTRS and user data, when the target cyclic matrix is used to estimate the ICI caused by the phase noise, the accuracy of the ICI estimation can be improved to a certain extent.

Based on the above embodiment shown in FIG. 2 , in the above S204 , when using the target cyclic matrix to estimate the ICI caused by the phase noise, in order to further improve the accuracy of the ICI estimation, it is also possible to further combine the construction of the num dimension Column vector, where each element in the column vector is the decision estimation of the time-frequency signal for the constellation points, and the target cyclic matrix and the constructed column vector are used to jointly estimate the inter-carrier interference caused by the phase noise. 4 is a schematic flowchart of another phase noise estimation method provided by an embodiment of the present application, and the phase noise estimation method may further include:

S401. Construct a column vector.

Among them, the column vector is a num-dimensional column vector, and each element in the column vector is the judgment estimate of the time-frequency signal on the constellation point. By way of example, the decision estimation may be a least square estimation or a hard decision estimation, which may be specifically set according to actual needs.

Combined with the text description about the column vector in S401, to convert the column vector into a mathematical formula, please refer to the following formula 3:

Among them, Ls represents the decision estimate of R to the constellation point.

S402. Perform separation processing on the column vectors to obtain real column vectors and imaginary column vectors corresponding to the column vectors, and use the real column vectors and imaginary column vectors as matrix sub-blocks to construct block column vectors corresponding to the column vectors.

Exemplarily, when constructing a block column vector corresponding to a column vector corresponding to a column vector, in combination with the description in S401 above, when the column vector is a num-dimensional column vector, the corresponding block column vector is a 2num-dimensional column vector. ; wherein, the matrix sub-blocks of the first num rows of the block column vector are real part column vectors, and the matrix sub-blocks of the last num rows of the block column vector are imaginary part column vectors. The real part column vector corresponding to the block column vector may be represented by Real(Y), and the imaginary part column vector corresponding to the block column vector may be represented by Imag(Y).

Combined with the text description about the block column vector in S402, the block column vector is converted into a mathematical formula, as shown in the following formula 4:

in,

Represents the constructed block column vector, Real(Y) represents the real part column vector corresponding to the block column vector, and Imag(Y) represents the imaginary part column vector corresponding to the block column vector. It can be understood that, in the embodiment of the present application, by separating the real part and the imaginary part corresponding to the column vector, the advantage is that the number and precision of row selection can be increased, which can better estimate the ICI.

get the chunked column vector in the constructor

After that, the block column vector can also be

to filter the rows in that chunked column vector

When screening the rows of , the row vectors in the block column vector can be screened according to the target cyclic matrix in the above S204 to obtain the target column vector, that is, the following S403 is performed:

S403. Perform screening processing on the row vectors in the block column vector according to the target cyclic matrix to obtain the target column vector.

Among them, the number of rows occupied by the s-th row vector in the target circulant matrix in the block cyclic matrix is the same as the number of rows occupied by the t-th row vector in the target column vector in the block column vector; where, s=t, and s is an integer greater than or equal to 1 and less than or equal to 2num.

Combined with the description in S204 above, when screening the row vectors in the block column vector, it is not necessary to screen the row vectors in the block circulant matrix, but directly determine each row vector in the target circulant matrix. The number of rows occupied in the block circulant matrix. Assuming that the first row vector in the target circulant matrix is the third row vector in the block circulant matrix, the third row vector in the block column vector is determined. is the first row vector of the target column vector; assuming that the second row vector in the target circulant matrix is the fourth row vector in the block circulant matrix, the fourth row vector in the block column vector is determined as The second row vector of the target column vector; and so on, the filtered target column vector. It can be seen that the number of rows occupied by the sth row vector in the target circulant matrix in the block circulant matrix is the same as the number of rows occupied by the tth row vector in the target column vector in the block column vector; where, s=t. It can be understood that, in the embodiment of the present application, by performing the screening process on the row vector in the block column vector, the advantage is that the accuracy of the ICI estimation can be improved to a certain extent.

After filtering to obtain the target column vector, the following S404 can be performed:

S404 , using the target circulant matrix and the target column vector to estimate the inter-carrier interference caused by the phase noise.

Exemplarily, when using the target circulant matrix and the target column vector to estimate the inter-carrier interference caused by the phase noise, the phase noise tap at the first moment can be determined according to I=(Q ^T Q) ^-1 Q ^T U; and according to the phase noise Noise taps estimate inter-carrier interference due to phase noise. where I is the phase noise tap, Q is the target circulant matrix, and U is the target column vector.

It can be seen that in the embodiment of the present application, by constructing a column vector and constructing a target column vector based on the column vector, the target cyclic matrix and the target column vector can be used to estimate the inter-carrier interference caused by the phase noise, which not only realizes the Estimating the inter-carrier interference caused by phase noise; and can improve the accuracy of inter-carrier interference estimation to a certain extent.

Based on the above embodiment shown in FIG. 2 or FIG. 4 , after estimating the inter-carrier interference caused by the phase noise, signal compensation may also be performed on the received time-frequency signal. For example, when performing signal compensation on the received time-frequency signal, a convolution operation may be performed on the preset sequence and the phase noise taps to obtain the result of the convolution operation; and the time-frequency signal may be signal compensated according to the result of the convolution operation. , which not only realizes the estimation of the inter-carrier interference caused by the phase noise, but also further realizes the signal compensation of the received time-frequency signal.

The preset sequence is a sequence of all subcarriers used at the first moment when the time-frequency signal is received. Combined with the description in Figure 2 above, the preset sequence can be R(i, 0), R(i, 1), R(i, 2), R(i, 3), ..., R(i, num-2 ), R(i, num-1).

FIG. 5 is a schematic structural diagram of a phase noise estimation apparatus 50 provided by an embodiment of the present application. For example, referring to FIG. 5 , the phase noise estimation apparatus 50 may include:

The receiving unit 501 is configured to receive an analog waveform signal affected by phase noise from a transmitting end.

The processing unit 502 is configured to demodulate the analog waveform signal to obtain a time-frequency signal corresponding to the analog waveform signal; and determine the phase noise estimation order P according to the time-frequency signal; P is greater than or equal to 1 and less than or equal to num The integer of num is used to indicate the number of subcarriers used when receiving the analog waveform signal.

The processing unit 502 is further configured to construct a circulant matrix according to the phase noise estimation order P and the number of subcarriers; wherein, the elements of the circulant matrix include the phase tracking pilot PTRS and user data.

The estimation unit 503 is used for estimating the inter-carrier interference caused by the phase noise by using a cyclic matrix.

Optionally, the processing unit 502 is further configured to perform separation processing on the circulant matrix, to obtain the real part circulant matrix and the imaginary part circulant matrix corresponding to the circulant matrix; and use the real part circulant matrix and the imaginary part circulant matrix as matrix sub-blocks to construct The block circulant matrix corresponding to the circulant matrix.

The estimation unit 503 is specifically configured to estimate the inter-carrier interference caused by the phase noise by using the block cyclic matrix.

Optionally, the circulant matrix is a matrix with num rows and 2P+1 columns, and the block circulant matrix is a matrix with 2num rows and 4P+2 columns.

Optionally, the processing unit 502 is further configured to perform screening processing on the row vectors in the block circulant matrix to obtain the target circulant matrix.

The estimation unit 503 is specifically configured to use the target cyclic matrix to estimate the inter-carrier interference caused by the phase noise.

Optionally, the processing unit 502 is specifically configured to determine whether the k-th element in the P-th column of the block cyclic matrix is a PTRS; k is any value in an integer greater than or equal to 1 and less than or equal to 2num If the kth element is PTRS, then keep the kth row element; if the kth element is non-PTRS, then judge whether the value of the kth element is greater than the preset threshold; if the value of the kth element is greater than the preset threshold , then keep the k-th row element; if the value of the k-th element is less than or equal to the preset threshold, delete the k-th row element in the block circular matrix to obtain the target circular matrix.

Optionally, the processing unit 502 is further configured to construct a column vector, where the column vector is a num-dimensional column vector, and each element in the column vector is a judgment estimate of the constellation point by the time-frequency signal; the column vector is separated and processed, Obtain the real part column vector and the imaginary part column vector corresponding to the column vector; take the real part column vector and the imaginary part column vector as matrix sub-blocks, and construct the block column vector corresponding to the column vector.

The estimation unit 503 is specifically configured to estimate the inter-carrier interference caused by the phase noise by using the target cyclic matrix and the block column vector.

Optionally, the block column vector is a 2num-dimensional column vector.

Optionally, the processing unit 502 is further configured to perform screening processing on the row vectors in the block column vector according to the target circulant matrix to obtain the target column vector; The number of rows is the same as the number of rows occupied by the t-th row vector in the target column vector in the block column vector; where s=t, and s is an integer greater than or equal to 1 and less than or equal to 2num.

The estimation unit 503 is specifically configured to use the target cyclic matrix and the target column vector to estimate the inter-carrier interference caused by the phase noise.

Optionally, the estimation unit 503 is specifically configured to determine the phase noise tap at the first moment according to I=(Q ^T Q) ⁻¹ Q ^T U; and estimate the inter-carrier interference caused by the phase noise according to the phase noise tap; wherein , I denotes the phase noise tap, Q denotes the target circulant matrix, and U denotes the target column vector.

Optionally, the processing unit 502 is further configured to perform a convolution operation on the preset sequence and the phase noise taps to obtain a result of the convolution operation; and perform signal compensation on the time-frequency signal according to the result of the convolution operation; wherein the preset sequence is: The sequence of all sub-carriers used at the first moment when the time-frequency signal is received.

Optionally, the circulant matrix is a matrix with num rows and 2P+1 columns.

The phase noise estimation apparatus 50 shown in the embodiment of the present application can execute the phase noise estimation method in the embodiment shown in any of the above drawings, and its implementation principle and beneficial effects are similar to those of the phase noise estimation method, No further description is given here.

FIG. 6 is a schematic structural diagram of a communication apparatus 60 provided by an embodiment of the present application. For an example, please refer to FIG. 6 . The communication apparatus 60 may include a processor 601 and a memory 602, where a computer program is stored in the memory 602 , the processor 601 executes the computer program stored in the memory 602, so that the communication device 60 executes the phase noise estimation method in the embodiment shown in any of the above drawings, its realization principle, beneficial effects and phase noise The realization principles and beneficial effects of the estimation method are similar, and are not repeated here.

An embodiment of the present application further provides a communication apparatus, and the communication apparatus may include: a processor and an interface circuit.

The processor is used to run the code instructions to execute the phase noise estimation method in the embodiment shown in any of the above drawings, and its implementation principle and beneficial effects are similar to those of the phase noise estimation method. It will not be repeated here.

The embodiment of the present application also provides a chip, a computer program is stored on the chip, when the computer program is executed by the processor, the phase noise estimation method in the embodiment shown in any of the above drawings is executed, and the implementation principle and beneficial effects thereof are implemented. The implementation principle and beneficial effects of the phase noise estimation method are similar, and are not repeated here.

Embodiments of the present application further provide a computer-readable storage medium, where instructions are stored in the computer-readable storage medium, and when the instructions are executed on a communication device, the communication device is made to execute any of the above-mentioned drawings. The implementation principle and beneficial effects of the phase noise estimation method in the illustrated embodiment are similar to those of the phase noise estimation method, which will not be repeated here.

In the above-mentioned various embodiments, the processor may be a general-purpose processor, a digital signal processor (DSP), an application specific integrated circuit (ASIC), an off-the-shelf programmable gate array (field programmable gate array, FPGA) Or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components. The methods, steps, and logic block diagrams disclosed in the embodiments of this application can be implemented or executed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in conjunction with the embodiments of the present application may be directly embodied as executed by a hardware decoding processor, or executed by a combination of hardware and software modules in the decoding processor. Software modules can be located in random access memory (RAM), flash memory, read-only memory (ROM), programmable read-only memory or electrically erasable programmable memory, registers, etc. in the storage medium. The storage medium is located in the memory, and the processor reads the instructions in the memory, and completes the steps of the above method in combination with its hardware.

In the several embodiments provided in this application, it should be understood that the disclosed apparatus and method may be implemented in other manners. For example, the apparatus embodiments described above are only illustrative. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components may be combined or Can be integrated into another system, or some features can be ignored, or not implemented. On the other hand, the shown or discussed mutual coupling or direct coupling or communication connection may be through some interfaces, indirect coupling or communication connection of devices or units, and may be in electrical, mechanical or other forms.

The units described as separate components may or may not be physically separated, and components displayed as units may or may not be physical units, that is, may be located in one place, or may be distributed to multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution in this embodiment.

In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically alone, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware, or may be implemented in the form of hardware plus software functional units.

Claims

A phase noise estimation method, comprising:

Receive the phase noise-affected analog waveform signal from the transmitter;

demodulating the analog waveform signal to obtain a time-frequency signal corresponding to the analog waveform signal;

Determine the phase noise estimation order P according to the time-frequency signal; P is an integer greater than or equal to 1 and less than or equal to num, where num is used to represent the number of subcarriers used when receiving the analog waveform signal;

Construct a circulant matrix according to the phase noise estimation order P and the number of subcarriers; wherein, the elements of the circulant matrix include phase tracking pilot PTRS and user data;

The inter-carrier interference caused by the phase noise is estimated using the circulant matrix.
The method according to claim 1, wherein the using the cyclic matrix to estimate the inter-carrier interference caused by the phase noise comprises:

Separating the circulant matrix to obtain a real part circulant matrix and an imaginary part circulant matrix corresponding to the circulant matrix;

Using the real part circulant matrix and the imaginary part circulant matrix as matrix sub-blocks, construct the block circulant matrix corresponding to the circulant matrix;

The inter-carrier interference caused by the phase noise is estimated using the block cyclic matrix.
The method of claim 2, wherein:

The cyclic matrix is a matrix of num rows and 2P+1 columns, and the block cyclic matrix is a matrix of 2num rows and 4P+2 columns;

Wherein, the first num rows and the first 2P+1 columns of the block cyclic matrix are the real part circulant matrix, the first num rows of the block circulant matrix, and the last 2P+1 columns of the matrix subblock is the imaginary cyclic matrix, the last num rows of the block cyclic matrix, and the matrix sub-blocks of the first 2P+1 columns are the imaginary cyclic matrix, the last num rows of the block cyclic matrix, the last 2P+ The one-column matrix sub-block is the real circulant matrix.
The method according to claim 2 or 3, wherein the using the block cyclic matrix to estimate the inter-carrier interference caused by phase noise comprises:

Screening the row vectors in the block cyclic matrix to obtain a target cyclic matrix;

The inter-carrier interference caused by the phase noise is estimated using the target cyclic matrix.
The method according to claim 4, characterized in that, performing a screening process on the row vectors in the block circulant matrix to obtain a target circulant matrix, comprising:

Determine whether the k-th element in the P-th column in the block cyclic matrix is a PTRS; k is any value of an integer greater than or equal to 1 and less than or equal to 2num;

If the kth element is a PTRS, then keep the kth row element; if the kth element is a non-PTRS, then judge whether the value of the kth element is greater than a preset threshold;

If the value of the k-th element is greater than the preset threshold, keep the k-th row element; if the value of the k-th element is less than or equal to the preset threshold, execute the block loop Delete the k-th row element from the matrix to obtain the target cyclic matrix.
The method according to claim 4, wherein the using the target cyclic matrix to estimate the inter-carrier interference caused by the phase noise comprises:

Constructing a column vector, the column vector is a num-dimensional column vector, and each element in the column vector is a judgment estimate of the constellation point by the time-frequency signal;

Separating the column vector to obtain a real part column vector and an imaginary part column vector corresponding to the column vector;

Using the real part column vector and the imaginary part column vector as matrix sub-blocks, construct the block column vector corresponding to the column vector;

The inter-carrier interference caused by phase noise is estimated using the target circulant matrix and the block column vector.
The method of claim 6, wherein:

The block column vector is a 2num-dimensional column vector;

The matrix sub-blocks of the first num rows of the block column vector are the real part column vectors, and the matrix sub-blocks of the last num rows of the block column vector are the imaginary part column vectors.
The method according to claim 6 or 7, wherein the using the target circulant matrix and the block column vector to estimate the inter-carrier interference caused by phase noise comprises:

According to the target cyclic matrix, the row vectors in the block column vector are screened to obtain a target column vector; the number of rows occupied by the s-th row vector in the target cyclic matrix in the block cyclic matrix, is the same as the number of rows occupied by the t-th row vector in the target column vector in the block column vector; wherein, s=t, and s is an integer greater than or equal to 1 and less than or equal to 2num;

The inter-carrier interference caused by phase noise is estimated using the target circulant matrix and the target column vector.
The method according to claim 8, wherein the using the target cyclic matrix and the target column vector to estimate the inter-carrier interference caused by phase noise comprises:

Determine the phase noise tap at the first moment according to I=(Q T Q) -1 Q T U; wherein, I represents the phase noise tap, Q represents the target circulant matrix, and U represents the target column vector;

The phase noise induced inter-carrier interference is estimated from the phase noise taps.
The method according to claim 9, wherein the method further comprises:

Perform a convolution operation on the preset sequence and the phase noise taps to obtain a result of the convolution operation; wherein the preset sequence is a sequence of all subcarriers used at the first moment when receiving the time-frequency signal;

Signal compensation is performed on the time-frequency signal according to the result of the convolution operation.
The method according to any one of claims 1-10, wherein,

The circulant matrix is a matrix with num rows and 2P+1 columns;

Wherein, the first P elements of the first column of the circulant matrix are all 0, and the P+1th element to the numth element are sequentially the first num-P elements of the preset sequence; the second element of the circulant matrix The first P-1 elements of the column are all 0, and the Pth element to the numth element are sequentially the first num-P+1 elements of the preset sequence, and so on, until the second column of the cyclic matrix The first element of is 0, and the second element to the numth element are the first num-1 elements of the preset sequence;

The first element to the numth element of the P+1th of the circulant matrix are sequentially elements in the preset sequence;

The first num-1 elements of the P+2th column of the circulant matrix are sequentially the second element to the numth element of the preset sequence, and the numth element is 0; the P+3th column of the circulant matrix The first num-2 elements are sequentially the 3rd element to the numth element of the preset sequence, the num-1th element and the numth element are both 0, and so on, until the 2Pth element of the cyclic matrix The first num-P elements of the +1 column are sequentially the num-P+1th element to the numth element of the preset sequence, and the Pth element to the numth element are all 0;

Wherein, the preset sequence is R(i, 0), R(i, 1), R(i, 2), R(i, 3), ..., R(i, num-2), R(i , num-1).
A phase noise estimation device, comprising:

The receiving unit is used for receiving the analog waveform signal affected by the phase noise from the transmitting end;

a processing unit, configured to perform demodulation processing on the analog waveform signal to obtain a time-frequency signal corresponding to the analog waveform signal; and determine the phase noise estimation order P according to the time-frequency signal; P is greater than or equal to 1, and an integer less than or equal to num, where num is used to indicate the number of subcarriers used when receiving the analog waveform signal;

The processing unit is further configured to construct a circulant matrix according to the phase noise estimation order P and the number of subcarriers; wherein, the elements of the circulant matrix include phase tracking pilot PTRS and user data;

an estimation unit, configured to use the cyclic matrix to estimate the inter-carrier interference caused by the phase noise.
The device of claim 12, wherein:

The processing unit is further configured to perform separation processing on the circulant matrix to obtain a real part circulant matrix and an imaginary part circulant matrix corresponding to the circulant matrix; and use the real part circulant matrix and the imaginary part circulant matrix as Matrix sub-block, constructs the block circulant matrix corresponding to the circulant matrix;

The estimation unit is specifically configured to estimate the inter-carrier interference caused by the phase noise by using the block cyclic matrix.
The device of claim 13, wherein:

The cyclic matrix is a matrix of num rows and 2P+1 columns, and the block cyclic matrix is a matrix of 2num rows and 4P+2 columns;

Wherein, the first num rows and the first 2P+1 columns of the block cyclic matrix are the real part circulant matrix, the first num rows of the block circulant matrix, and the last 2P+1 columns of the matrix subblock is the imaginary cyclic matrix, the last num rows of the block cyclic matrix, and the matrix sub-blocks of the first 2P+1 columns are the imaginary cyclic matrix, the last num rows of the block cyclic matrix, the last 2P+ The one-column matrix sub-block is the real circulant matrix.
The device according to claim 13 or 14, characterized in that,

The processing unit is further configured to perform screening processing on the row vectors in the block cyclic matrix to obtain a target cyclic matrix;

The estimation unit is specifically configured to use the target cyclic matrix to estimate the inter-carrier interference caused by the phase noise.
The apparatus of claim 15, wherein:

The processing unit is specifically configured to determine whether the k-th element in the P-th column of the block cyclic matrix is a PTRS; k is any value in an integer greater than or equal to 1 and less than or equal to 2num; If the k-th element is a PTRS, the k-th row element is reserved; if the k-th element is a non-PTRS, it is determined whether the value of the k-th element is greater than a preset threshold; if the k-th element is not a PTRS If the value of the element is greater than the preset threshold, the k-th row element is retained; if the value of the k-th element is less than or equal to the preset threshold, the block circulant matrix is deleted. With k row elements, the target cyclic matrix is obtained.
The apparatus of claim 15, wherein:

The processing unit is further configured to construct a column vector, where the column vector is a num-dimensional column vector, and each element in the column vector is a judgment estimate of the time-frequency signal for constellation points; The vector is separated and processed to obtain the real part column vector and the imaginary part column vector corresponding to the column vector; the real part column vector and the imaginary part column vector are used as matrix sub-blocks, and the part corresponding to the column vector is constructed. block column vector;

The estimation unit is specifically configured to estimate the inter-carrier interference caused by the phase noise by using the target cyclic matrix and the block column vector.
The apparatus of claim 17, wherein:

The block column vector is a 2num-dimensional column vector;

The matrix sub-blocks of the first num rows of the block column vector are the real part column vectors, and the matrix sub-blocks of the last num rows of the block column vector are the imaginary part column vectors.
The device according to claim 17 or 18, characterized in that,

The processing unit is further configured to perform screening processing on the row vectors in the block column vector according to the target circulant matrix to obtain a target column vector; The number of rows occupied by the circular matrix is the same as the number of rows occupied by the t-th row vector in the target column vector in the block column vector; wherein, s=t, and s is greater than or equal to 1 and less than or an integer equal to 2num;

The estimation unit is specifically configured to use the target cyclic matrix and the target column vector to estimate the inter-carrier interference caused by the phase noise.
The apparatus of claim 19, wherein:

The estimation unit is specifically configured to determine the phase noise tap at the first moment according to I=(Q T Q) -1 Q T U; and estimate the inter-carrier interference caused by the phase noise according to the phase noise tap; wherein, I denotes the phase noise tap, Q denotes the target circulant matrix, and U denotes the target column vector.
The apparatus of claim 20, wherein:

The processing unit is further configured to perform a convolution operation on the preset sequence and the phase noise taps to obtain a convolution operation result; and perform signal compensation on the time-frequency signal according to the convolution operation result; wherein, the The preset sequence is a sequence of all subcarriers used at the first moment when the time-frequency signal is received.
The device according to any one of claims 12-21, characterized in that,

The circulant matrix is a matrix with num rows and 2P+1 columns;

Wherein, the first P elements of the first column of the circulant matrix are all 0, and the P+1th element to the numth element are sequentially the first num-P elements of the preset sequence; the second element of the circulant matrix The first P-1 elements of the column are all 0, and the Pth element to the numth element are sequentially the first num-P+1 elements of the preset sequence, and so on, until the second column of the cyclic matrix The first element of is 0, and the second element to the numth element are the first num-1 elements of the preset sequence;

The first element to the numth element of the P+1th of the circulant matrix are sequentially elements in the preset sequence;

The first num-1 elements of the P+2th column of the circulant matrix are sequentially the second element to the numth element of the preset sequence, and the numth element is 0; the P+3th column of the circulant matrix The first num-2 elements are sequentially the 3rd element to the numth element of the preset sequence, the num-1th element and the numth element are both 0, and so on, until the 2Pth element of the cyclic matrix The first num-P elements of the +1 column are sequentially the num-P+1th element to the numth element of the preset sequence, and the Pth element to the numth element are all 0;

Wherein, the preset sequence is R(i, 0), R(i, 1), R(i, 2), R(i, 3), ..., R(i, num-2), R(i , num-1).
A communication device, characterized in that the device comprises a processor and a memory, wherein a computer program is stored in the memory, and the processor executes the computer program stored in the memory, so that the device executes the method as claimed in the claims The phase noise estimation method described in any one of 1 to 11.
A communication device, comprising: a processor and an interface circuit;

the interface circuit for receiving code instructions and transmitting them to the processor;

The processor is configured to execute the code instructions to perform the phase noise estimation method according to any one of claims 1 to 11.
A computer-readable storage medium storing instructions in the computer-readable storage medium, characterized in that, when the instructions are executed on a communication device, the communication device is made to execute any one of claims 1 to 11 The described phase noise estimation method.