5G non-ground network millimeter wave uplink waveform phase noise estimation and compensation method and device

Abstract

The invention discloses a 5G non-ground network millimeter wave uplink waveform phase noise estimation and compensation method and device, and the method comprises the steps: obtaining a data symbol after channel equalization, carrying out the cyclic shift of a subcarrier index in the data symbol according to a timing error estimation value, and obtaining a subcarrier index after cyclic shift, extracting a phase tracking reference signal group in the data symbol, and performing conjugate multiplication on the phase tracking reference signal group and a local phase tracking reference signal group to obtain a phase noise estimation value of the phase tracking reference signal group; selecting a data subcarrier of which the power value exceeds a power threshold value as an auxiliary interpolation point, and calculating a phase noise estimation value of the data subcarrier; and performing interpolation calculation on the phase noise estimation value of the phase tracking reference signal group and the phase noise estimation value of the auxiliary interpolation point to obtain a phase noise interpolation result of all subcarriers in the data symbol, and performing phase noise compensation on all subcarriers in the data symbol to obtain a result after phase noise compensation. According to the invention, the bit error rate performance of phase noise estimation and compensation can be improved.

Description

Technical Field

[0001] The present invention belongs to the field of communication technology, and specifically relates to a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method and device. Background Art

[0002] In recent years, with the development of 5G standards and the introduction of concepts such as broadband satellite Internet and non-terrestrial networks (NTN), broadband satellite communications based on 5G NTN have become a hot topic in current research, and the integration of satellite communications and traditional terrestrial mobile communications to build an integrated air-ground-space network has become a trend in the development of future mobile communications. At the same time, with the development of mobile communication technology, the existing low-frequency band (Sub-6GHz) spectrum resources are becoming increasingly scarce. The millimeter wave band (30GHz~300GHz) can provide a larger communication bandwidth and higher transmission rate than the traditional low-frequency band. 5G millimeter wave broadband satellite communications have become a key technology in the construction of non-terrestrial networks.

[0003] Phase noise refers to the noise interference caused by oscillation of various RF hardware, which will cause the phase of the signal to change, thus affecting the performance of the system. Generally speaking, the impact of phase noise on low frequency bands is small. However, with the increase of operating frequency, especially in high frequency bands above 6GHz, the phase noise generated by the non-ideal nature of hardware such as oscillators cannot be ignored, which has a great impact on high-order constellation demodulation and seriously affects the bit error rate performance of the system. Therefore, it is very important to design an accurate phase noise estimation and compensation algorithm at the receiving end. For this reason, the 5G protocol specifically defines the phase tracking reference signal (PTRS) for uplink waveform phase noise estimation and compensation. In the uplink of 5G non-terrestrial networks, in order to reduce the peak-to-average ratio of the signal, DFT-S-OFDM waveform is generally used for transmission. Compared with the CP-OFDM waveform, the transform precoding operation is added in the frequency domain. At this time, the phase noise is manifested as random phase rotation at each constellation point.

[0004] However, in 5G millimeter wave broadband systems, there are many data subcarriers and the phase noise changes quickly. Relying only on the phase estimation value at the PTRS for interpolation will not be effective in estimating and compensating the phase noise. In addition, timing errors will inevitably occur during uplink reception, which will cause the data phase in the symbol to cyclically shift. This phase discontinuity will further deteriorate the estimation and compensation results. Therefore, improving the accuracy of phase noise estimation and compensation has become an urgent problem to be solved. Summary of the invention

[0005] In order to solve the above problems existing in the prior art, the present invention provides a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method and device. The technical problem to be solved by the present invention is achieved by the following technical solutions:

[0006] In a first aspect, the present invention provides a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method, comprising:

[0007] Acquire the data symbol after channel equalization, and cyclically shift the subcarrier index in the data symbol according to the timing error estimation value to obtain the subcarrier index after cyclic shift;

[0008] Extracting a phase tracking reference signal group in the data symbol according to the subcarrier index after cyclic shift, and conjugate-multiplying the phase tracking reference signal group with the local phase tracking reference signal group to obtain a phase noise estimation value of the phase tracking reference signal group;

[0009] Selecting data subcarriers whose power values ​​exceed a power threshold from the data symbols as auxiliary interpolation points, and calculating phase noise estimation values ​​of the auxiliary interpolation points;

[0010] Interpolate the phase noise estimation value of the phase tracking reference signal group and the phase noise estimation value of the auxiliary interpolation point to obtain the phase noise interpolation result of all subcarriers in the data symbol;

[0011] According to the phase noise interpolation result, phase noise compensation is performed on all subcarriers in the data symbol to obtain a phase noise compensated result.

[0012] In a second aspect, the present invention further provides a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation device, comprising:

[0013] A data acquisition module is used to acquire data symbols after channel equalization, and cyclically shift the subcarrier index in the data symbol according to the timing error estimation value to obtain the subcarrier index after cyclic shift;

[0014] A data processing module is used to extract a phase tracking reference signal group in the data symbol according to the subcarrier index after cyclic shift, and conjugate-multiply the phase tracking reference signal group with the local phase tracking reference signal group to obtain a phase noise estimation value of the phase tracking reference signal group;

[0015] A data auxiliary processing module, used for selecting data subcarriers whose power values ​​exceed a power threshold from data symbols as auxiliary interpolation points, and calculating phase noise estimation values ​​of the auxiliary interpolation points;

[0016] A phase noise interpolation result calculation module is used to interpolate the phase noise estimation value of the phase tracking reference signal group and the phase noise estimation value of the auxiliary interpolation point to obtain the phase noise interpolation results of all subcarriers in the data symbol;

[0017] The phase noise compensation module is used to perform phase noise compensation on all subcarriers in the data symbol according to the phase noise interpolation result to obtain the phase noise compensated result.

[0018] Beneficial effects of the present invention:

[0019] The present invention provides a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method and device, which take into account the impact of the receiving timing error on the data phase, and use data-assisted interpolation to improve the interpolation accuracy of the phase noise estimation. It can effectively improve the bit error rate performance of the phase noise estimation and compensation, and alleviate the data phase discontinuity problem caused by the receiving timing error. It has low complexity and is easy to implement.

[0020] The present invention will be further described in detail below with reference to the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] Figure 1 It is a flow chart of a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method provided by an embodiment of the present invention;

[0022] Figure 2 It is a schematic diagram of a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method provided by an embodiment of the present invention;

[0023] Figure 3 It is a schematic diagram of a power spectral density curve of a phase noise model provided by an embodiment of the present invention under the condition that the carrier frequency is 45 GHz;

[0024] Figure 4 It is a schematic diagram of comparing the bit error rate of the method provided by the embodiment of the present invention with that of the traditional method at a signal-to-noise ratio of 17-28 dB;

[0025] Figure 5 It is a schematic diagram of a comparison between the minimum signal-to-noise ratio required for achieving a block error rate of 10% by the proposed method and the traditional method when the modulation coding scheme is MCS25 under the condition that different timing errors exist at the receiving end provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0026] The present invention is further described in detail below with reference to specific embodiments, but the embodiments of the present invention are not limited thereto.

[0027] See also Figure 1 and Figure 2 , Figure 1 is a flow chart of a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method provided by an embodiment of the present invention, Figure 21 is a schematic diagram of a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method provided by an embodiment of the present invention. The present invention provides a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method, including:

[0028] S101. Acquire data symbols after channel equalization, and cyclically shift subcarrier indexes in the data symbols according to a timing error estimation value to obtain cyclically shifted subcarrier indexes.

[0029] Specifically, in this embodiment, the cyclically shifted subcarrier index includes the cyclically shifted phase tracking reference signal subcarrier index and the cyclically shifted data subcarrier index.

[0030] Cyclic shifting of a subcarrier index in a data symbol according to a timing error estimation value to obtain a cyclically shifted subcarrier index includes:

[0031] According to the estimated value of the timing error, the cyclic shift length m is calculated, and its expression is:

[0032]

[0033] Among them, N STO represents the estimated value of the timing error, N SC Indicates the number of data symbol subcarriers, N FFT Indicates the number of Fourier transform points when demodulating data symbols, [·] indicates rounding to the nearest integer;

[0034] According to the cyclic shift length m, calculate the subcarrier index of the phase tracking reference signal after cyclic shift and the cyclically shifted data subcarrier index The expressions are:

[0035]

[0036] Among them, k ptrs represents the subcarrier index of the phase tracking reference signal before cyclic shift, k data Represents the data subcarrier index before cyclic shift, rem() N It represents the remainder modulo N.

[0037] S102. Extract a phase tracking reference signal group in the data symbol according to the subcarrier index after cyclic shift, and conjugate-multiply the phase tracking reference signal group with the local phase tracking reference signal group to obtain a phase noise estimation value of the phase tracking reference signal group.

[0038] Specifically, in this embodiment, the phase noise estimation value θ of the phase tracking reference signal group ptrs The expression is:

[0039]

[0040] in, represents the subcarrier index of the i-th phase tracking reference signal group, N S represents the number of phase tracking reference signal symbols in each phase tracking reference signal group, r(n) represents the phase tracking reference signal group in the extracted data symbol, t(n) represents the local phase tracking reference signal group, * represents conjugate, and arg represents angle operation.

[0041] In this embodiment, it also includes:

[0042] Get the phase noise estimate subcarrier index P of the phase tracking reference signal group ptrs , whose expression is:

[0043]

[0044] in, Indicates taking the minimum value of the subcarrier index in each phase tracking reference signal group, N G Indicates the number of phase tracking reference signal groups in each data symbol.

[0045] S103: Select data subcarriers whose power values ​​exceed a power threshold from the data symbols as auxiliary interpolation points, and calculate phase noise estimation values ​​of the auxiliary interpolation points.

[0046] Specifically, in this embodiment, the phase noise estimation value of the auxiliary interpolation point is calculated by the following process, including:

[0047] Calculate the power value S of the data subcarrier index in the data symbol, and its expression is:

[0048]

[0049] Among them, S( i ) represents the power value of the i-th data subcarrier index, and d(i) represents the data of the i-th data subcarrier;

[0050] Select the data subcarrier whose power value exceeds the power threshold from the data symbol as the auxiliary interpolation point, and calculate the data subcarrier index P corresponding to the auxiliary interpolation point data , whose expression is:

[0051]

[0052] Among them, S th Indicates the power threshold;

[0053] According to the positive and negative of the real and imaginary parts of the data subcarrier corresponding to the auxiliary interpolation point, the reference constellation point q of the data subcarrier index corresponding to the auxiliary interpolation point is calculated, and its expression is:

[0054]

[0055] Among them, q( i ) represents the reference constellation point of the i-th point in the data subcarrier index corresponding to the auxiliary interpolation point;

[0056] The data of the data subcarrier corresponding to the auxiliary interpolation point is conjugate-multiplied with the reference constellation point indexed by the data subcarrier corresponding to the auxiliary interpolation point to obtain the phase noise estimate θ of the auxiliary interpolation point. data , whose expression is:

[0057] θ data ( i ) = arg ( d(i)q * (i) ) ,i∈P data ;

[0058] Among them, θ data ( i ) represents the phase noise estimation value of the i-th auxiliary interpolation point.

[0059] S104, interpolating the phase noise estimation value of the phase tracking reference signal group and the phase noise estimation value of the auxiliary interpolation point to obtain the phase noise interpolation result of all subcarriers in the data symbol.

[0060] Specifically, in this embodiment, the expression of the phase noise interpolation result θ of all subcarriers in the data symbol is:

[0061]

[0062] Among them, θ( i ) represents the phase noise interpolation result of the i-th subcarrier in the data symbol, and P represents the phase noise estimation value subcarrier index P of the phase tracking reference signal group ptrs The data subcarrier index P corresponding to the auxiliary interpolation point data P1 and P2 represent the two subcarrier indices closest to the current auxiliary interpolation point in the set P, and the distance between P1 and the current auxiliary interpolation point is less than the distance between P2 and the current auxiliary interpolation point. P1 and θ P2 They represent the phase noise estimates of P1 and P2 respectively.

[0063] S105 . Perform phase noise compensation on all subcarriers in the data symbol according to the phase noise interpolation result to obtain a phase noise compensated result.

[0064] Specifically, in this embodiment, the phase noise compensation result r comp The expression is:

[0065]

[0066] Where r(i) represents the i-th received data on the data symbol, i = 0, ..., N SC -1.

[0067] In summary, the present invention provides a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method, which takes into account the impact of the receiving timing error on the data phase, and uses data-assisted interpolation to improve the interpolation accuracy of the phase noise estimation process. It can effectively improve the bit error rate performance after phase noise estimation and compensation, and alleviate the data phase discontinuity problem caused by the receiving timing error. It has low complexity and is easy to implement.

[0068] Based on the same inventive concept, the present invention also provides a 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation device, which is used to implement the 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method provided by the above embodiment of the present invention. The embodiment of the method is referred to above and will not be repeated here. The device includes:

[0069] A data acquisition module is used to acquire data symbols after channel equalization, and cyclically shift the subcarrier index in the data symbol according to the timing error estimation value to obtain the subcarrier index after cyclic shift;

[0070] A data processing module is used to extract a phase tracking reference signal group in the data symbol according to the subcarrier index after cyclic shift, and conjugate-multiply the phase tracking reference signal group with the local phase tracking reference signal group to obtain a phase noise estimation value of the phase tracking reference signal group;

[0071] A data auxiliary processing module, used for selecting data subcarriers whose power values ​​exceed a power threshold from data symbols as auxiliary interpolation points, and calculating phase noise estimation values ​​of the auxiliary interpolation points;

[0072] A phase noise interpolation result calculation module is used to interpolate the phase noise estimation value of the phase tracking reference signal group and the phase noise estimation value of the auxiliary interpolation point to obtain the phase noise interpolation results of all subcarriers in the data symbol;

[0073] The phase noise compensation module is used to perform phase noise compensation on all subcarriers in the data symbol according to the phase noise interpolation result to obtain the phase noise compensated result.

[0074] In an optional embodiment of the present invention, the effect of the 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method provided in the above embodiment is verified by simulation experiments, specifically:

[0075] 1. Simulation conditions

[0076] The simulation experiment of this embodiment is carried out under MATLAB R2024b software. The PTRS sequence is generated using the Gold pseudo-random sequence specified in the 5GNR physical layer protocol, and the channel used in the simulation is an additive white Gaussian noise channel.

[0077] 2. Simulation content and result analysis

[0078] The system parameters are: subcarrier spacing of 120kHz, cyclic prefix type of normal, frequency domain resource occupancy of 256RB, orthogonal frequency division multiplexing (OFDM) modulation and demodulation points of 4096, modulation and coding scheme of MCS25 (64QAM modulation), number of PTRS groups of 8, number of PTRS symbols in each PTRS group of 4, and carrier frequency of 45GHz.

[0079] See also Figure 3 and Figure 4 , Figure 3 4 is a schematic diagram of a power spectral density curve of a phase noise model provided by an embodiment of the present invention under the condition that the carrier frequency is 45 GHz. Figure 4 FIG. 1 is a schematic diagram showing a comparison of the bit error rate of the method provided by the present invention and the conventional method under a signal-to-noise ratio of 17 to 28 dB provided by an embodiment of the present invention. Figure 3 The phase noise shown is affected by Figure 4 It can be seen that the method proposed in the present invention has significantly improved bit error rate performance compared with the traditional method, and can effectively improve the effect of phase noise estimation and compensation.

[0080] See also Figure 5 , Figure 5 FIG. 1 is a schematic diagram of a comparison between the minimum signal-to-noise ratio required for achieving a block error rate of 10% by the proposed method and the conventional method when the modulation coding scheme is MCS25 under different timing errors at the receiving end provided by the embodiment of the present invention, wherein Figure 5It can be seen that the signal-to-noise ratio threshold required by the traditional method in ideal synchronization is higher than that of the method proposed in the present invention, and the threshold changes significantly with the change of the timing error, while the signal-to-noise ratio threshold of the method proposed in the present invention remains basically unchanged under timing errors of various sizes. It can be seen that the phase noise compensation performance of the method proposed in the present invention under different timing errors is stable.

[0081] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply that there is any such actual relationship or order between these entities or operations. Moreover, the term "include", "comprise" or any other variant is intended to cover non-exclusive inclusion, so that the article or device including a series of elements includes not only those elements, but also other elements that are not explicitly listed. In the absence of more restrictions, the elements defined by the sentence "including one..." do not exclude the existence of other identical elements in the article or device including the elements. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. The orientation or position relationship indicated by "up", "down", "left", "right", etc. is based on the orientation or position relationship shown in the accompanying drawings, which is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation of the present invention.

[0082] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification.

[0083] The above contents are further detailed descriptions of the present invention in combination with specific preferred embodiments, and it cannot be determined that the specific implementation of the present invention is limited to these descriptions. For ordinary technicians in the technical field to which the present invention belongs, several simple deductions or substitutions can be made without departing from the concept of the present invention, which should be regarded as falling within the protection scope of the present invention.

Claims

1. A 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method, characterized in that: include: Acquire data symbols after channel equalization, and cyclically shift subcarrier indexes in the data symbols according to the timing error estimation value to obtain cyclically shifted subcarrier indexes; Extracting a phase tracking reference signal group in the data symbol according to the cyclically shifted subcarrier index, and conjugate-multiplying the phase tracking reference signal group with a local phase tracking reference signal group to obtain a phase noise estimation value of the phase tracking reference signal group; Selecting data subcarriers whose power values ​​exceed a power threshold from the data symbols as auxiliary interpolation points, and calculating phase noise estimation values ​​of the auxiliary interpolation points; Interpolate the phase noise estimation value of the phase tracking reference signal group and the phase noise estimation value of the auxiliary interpolation point to obtain the phase noise interpolation result of all subcarriers in the data symbol; According to the phase noise interpolation result, phase noise compensation is performed on all subcarriers in the data symbol to obtain a phase noise compensated result.

2. The 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method according to claim 1 is characterized in that: The cyclically shifted subcarrier index includes a cyclically shifted phase tracking reference signal subcarrier index and a cyclically shifted data subcarrier index.

3. The 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method according to claim 2 is characterized in that: Cyclic shifting the subcarrier index in the data symbol according to the timing error estimation value to obtain the cyclically shifted subcarrier index includes: According to the estimated value of the timing error, the cyclic shift length m is calculated, and its expression is: Among them, N STO represents the estimated value of the timing error, N SC Indicates the number of data symbol subcarriers, N FFT Indicates the number of Fourier transform points when demodulating data symbols, [·] indicates rounding to the nearest integer; According to the cyclic shift length m, the phase tracking reference signal subcarrier index after cyclic shift is calculated. and the cyclically shifted data subcarrier index The expressions are: Among them, k ptrs represents the subcarrier index of the phase tracking reference signal before cyclic shift, k data Represents the data subcarrier index before cyclic shift, rem() N It represents the remainder modulo N.

4. The 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method according to claim 1 is characterized in that: The phase noise estimation value θ of the phase tracking reference signal group ptrs The expression is: in, represents the subcarrier index of the i-th phase tracking reference signal group, N S represents the number of phase tracking reference signal symbols in each phase tracking reference signal group, r(n) represents the phase tracking reference signal group in the extracted data symbol, t(n) represents the local phase tracking reference signal group, * represents conjugate, and arg represents angle operation.

5. The 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method according to claim 4 is characterized in that: Also includes: Obtain the phase noise estimation value subcarrier index P of the phase tracking reference signal group ptrs , whose expression is: in, Indicates taking the minimum value of the subcarrier index in each phase tracking reference signal group, N G Indicates the number of phase tracking reference signal groups in each data symbol.

6. The 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method according to claim 1 is characterized in that: The method further comprises: selecting data whose power value exceeds a power threshold from the data symbols as auxiliary interpolation points, and calculating a phase noise estimation value of the auxiliary interpolation point, comprising: The power value S of the data subcarrier index in the data symbol is calculated, and its expression is: Among them, S( i ) represents the power value of the i-th data subcarrier index, and d(i) represents the data of the i-th data subcarrier; Select a data subcarrier whose power value exceeds the power threshold from the data symbol as an auxiliary interpolation point, and calculate the data subcarrier index P corresponding to the auxiliary interpolation point data , whose expression is: Among them, S th Indicates the power threshold; The data of the data subcarrier corresponding to the auxiliary interpolation point is conjugate-multiplied with the reference constellation point indexed by the data subcarrier corresponding to the auxiliary interpolation point to obtain the phase noise estimation value θ of the auxiliary interpolation point data , whose expression is: θ data ( i )=arg ( d(i)q * (and) ) ,i∈P data ; Among them, θ data ( i ) represents the phase noise estimation value of the i-th auxiliary interpolation point.

7. The 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method according to claim 6 is characterized in that: The process of acquiring the reference constellation point of the data subcarrier index corresponding to the auxiliary interpolation point includes: According to the positive and negative of the real and imaginary parts of the data subcarrier corresponding to the auxiliary interpolation point, the reference constellation point q of the data subcarrier index corresponding to the auxiliary interpolation point is calculated, and its expression is: Among them, q( i ) represents the reference constellation point of the i-th point in the data subcarrier index corresponding to the auxiliary interpolation point.

8. The 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method according to claim 1 is characterized in that: The expression of the phase noise interpolation result θ of all subcarriers in the data symbol is: Among them, θ( i ) represents the phase noise interpolation result of the i-th subcarrier in the data symbol, and P represents the phase noise estimation value subcarrier index P of the phase tracking reference signal group ptrs The data subcarrier index P corresponding to the auxiliary interpolation point data P1 and P2 represent the two subcarrier indices closest to the current auxiliary interpolation point in the set P, respectively. P1 and θ P2 They represent the phase noise estimates of P1 and P2 respectively.

9. The 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation method according to claim 1 is characterized in that: The phase noise compensation result is comp The expression is: Where r(i) represents the i-th received data on the data symbol, i = 0, ..., N SC -1.

10. A 5G non-terrestrial network millimeter wave uplink waveform phase noise estimation and compensation device, characterized in that: include: A data acquisition module, used to acquire data symbols after channel equalization, and cyclically shift the subcarrier index in the data symbols according to the timing error estimation value to obtain the subcarrier index after cyclic shift; A data processing module, configured to extract a phase tracking reference signal group in the data symbol according to the cyclically shifted subcarrier index, and conjugate-multiply the phase tracking reference signal group with a local phase tracking reference signal group to obtain a phase noise estimation value of the phase tracking reference signal group; A data auxiliary processing module, used to select data subcarriers whose power values ​​exceed a power threshold from the data symbols as auxiliary interpolation points, and calculate phase noise estimation values ​​of the auxiliary interpolation points; A phase noise interpolation result calculation module, used to interpolate the phase noise estimation value of the phase tracking reference signal group and the phase noise estimation value of the auxiliary interpolation point to obtain the phase noise interpolation results of all subcarriers in the data symbol; The phase noise compensation module is used to perform phase noise compensation on all subcarriers in the data symbol according to the phase noise interpolation result to obtain a phase noise compensated result.

Images