CN115955380A - Phase noise estimation and compensation method - Google Patents

Phase noise estimation and compensation method Download PDF

Info

Publication number
CN115955380A
CN115955380A CN202211722484.1A CN202211722484A CN115955380A CN 115955380 A CN115955380 A CN 115955380A CN 202211722484 A CN202211722484 A CN 202211722484A CN 115955380 A CN115955380 A CN 115955380A
Authority
CN
China
Prior art keywords
phase noise
phase
phase angle
value
correction factor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211722484.1A
Other languages
Chinese (zh)
Inventor
刘伟利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Innogence Technology Co Ltd
Original Assignee
Sichuan Innogence Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sichuan Innogence Technology Co Ltd filed Critical Sichuan Innogence Technology Co Ltd
Priority to CN202211722484.1A priority Critical patent/CN115955380A/en
Publication of CN115955380A publication Critical patent/CN115955380A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Noise Elimination (AREA)

Abstract

The invention discloses a phase noise estimation and compensation method, which comprises the steps of determining an initial value of phase noise through a received pilot signal and a local pilot signal; further determining a phase angle estimation value according to the initial phase noise estimation value; traversing the phase angle estimation value to determine a flag (k, l), and simultaneously performing grouping calculation on the frequency domain data to obtain a phase angle correction factor; reversely updating the phase angle estimation value and the phase noise estimation value by using the phase angle correction factor; and performing phase noise compensation on the received signal by using the phase noise estimation value. According to the method, the phase noise correction factor is calculated in a mode of grouping frequency domain data, so that the reliability of the phase angle correction factor is improved; the phase angle estimation value is updated by using the phase angle correction factor, so that the problem of large phase angle estimation error under low signal-to-noise ratio can be effectively solved, and the accuracy of phase noise estimation is improved.

Description

Phase noise estimation and compensation method
Technical Field
The invention relates to the technical field of communication, in particular to a phase noise estimation and compensation method.
Background
As the requirements for transmission rate and spectrum efficiency of communication systems are increasing, the communication systems will use higher order Quadrature Amplitude Modulation (QAM) and larger system bandwidth to meet these requirements. The high-order QAM modulated signal is sensitive to phase noise due to the increased density of modulation constellation points and the small phase error tolerance, and the system performance is severely deteriorated by the phase noise with small intensity. In addition, due to the trend of higher operating frequency band and larger system bandwidth of the communication system, the influence of phase noise on the system will become more serious. Therefore, how to suppress the influence of the phase noise is crucial to the performance of the high-order quadrature amplitude modulation communication system.
In the prior art, the influence of phase noise on the system performance is mainly suppressed by estimating and compensating the phase noise. And the estimation compensation method can be generally classified into a pilot data aided method and a decision data oriented method. The performance of the pilot frequency data auxiliary method is closely related to the density of pilot frequency insertion, the quantity of the inserted pilot frequency and the like, and the accuracy of phase noise estimation can be effectively improved by increasing the density and the quantity of the pilot frequency insertion, but the transmission efficiency of the system is greatly reduced. The decision data-oriented method can effectively improve the transmission efficiency of the system, but the error propagation problem can seriously affect the performance of the high-order QAM system.
Patent CN112688891A provides a method for phase noise estimation and compensation using PTRS signals in time domain. Patent CN110460385B provides a method for estimating phase noise from pilot signals, which has no targeted improvement effect on the increase of phase noise estimation error caused by noise.
Disclosure of Invention
The technical scheme of the invention solves the problem of performance loss caused by inaccurate phase noise estimation due to noise under low signal-to-noise ratio, and improves the performance of phase noise estimation and compensation by judging the reliability of the preliminarily estimated phase noise.
The purpose of the invention is realized by the following technical scheme:
assuming that L pilot signals are configured in the time domain, here, the pilot signal is referred to as PTRS, each PTRS signal is configured with K REs, and the frequency domain data of the PTRS on each signal is divided into M groups, that is, the frequency domain data of every N = K/M REs is a group.
A phase noise estimation and compensation method comprises the following specific steps:
s1: determining an initial noise estimation value through the received pilot signal and a local pilot signal;
s2: further determining a phase angle estimation value according to the initial phase noise estimation value;
s3: traversing the phase angle estimation value to determine a flag (k, l), and simultaneously performing grouping calculation on the frequency domain data to obtain a phase angle correction factor;
s4: reversely updating the phase angle estimation value and the phase noise estimation value by using the phase angle correction factor;
s5: and performing phase noise compensation on the received signal by using the phase noise estimation value.
Further, the specific calculation process in step S1 is as follows:
Figure BDA0004030100570000021
/>
wherein, Y ptrs( k, l) is the received pilot signal, X ptrs (k, l) is a local pilot signal, H temp (k, l) is the initial estimate of phase noise, and l is in the range of [1]L is the number of pilot signals configured in the time domain, and k has a value range of [1, K ]]K is the number of REs configured for each pilot signal.
Further, the specific calculation process in step S2 is as follows:
Figure BDA0004030100570000022
wherein,
Figure BDA0004030100570000023
for the phase angle estimation, H temp (k, l) is the initial estimate of phase noise, and l is in the range of [1]L is the number of pilot signals configured in the time domain, and k has a value range of [1, K ]]And K is the number of REs configured for each pilot signal.
Further, the step S3 of traversing the estimated phase angle value to determine flag (k, l) specifically includes: ergodic phase angle estimation
Figure BDA0004030100570000024
If>
Figure BDA0004030100570000025
Then flag (k, l) =1; if/or>
Figure BDA0004030100570000026
Then flag (k, l) =0.
Further, the grouping of the frequency domain data in step S3 specifically includes: dividing the flag (k, l) corresponding to the ith signal into M groups according to the frequency domain configuration of the PTRS, and carrying out XOR operation on N flags (k, l) in each group:
Figure BDA0004030100570000027
wherein, F (m, l) represents the mth group phase noise estimation reliability flag of the ith pilot symbol, the range of N is [1, N ], the range of m is [1, m ], if N flag (k, l) values in each group are the same, F (m, l) =0, otherwise, F (m, l) =1.
Further, the phase angle correction factor calculated in step S3 is specifically:
Figure BDA0004030100570000031
wherein, F (M, L) represents the mth group phase noise estimation reliability flag of the mth pilot symbol, M is an index, the value range is [1, M ], M is the frequency domain data grouping number, L is an index, the value range is [1, L ], L is the pilot signal number configured on the time domain, and β is the phase angle correction factor.
Further, step S4 specifically includes:
calculating a phase angle estimate for each signal based on the phase angle correction factor
Figure BDA0004030100570000032
Figure BDA0004030100570000033
Wherein beta is a phase angle correction factor, K ptrs Is the Kth RE of the pilot signal, K is the index, l is the index, the value range is [1, L],H ptrs (k, l) is the phase noise estimate;
obtaining the final estimated value H of the phase noise ptrs (k,l):
Figure BDA0004030100570000034
The invention has the beneficial effects that:
(1) The phase noise correction factor is calculated in a mode of grouping frequency domain data, so that the reliability of the phase angle correction factor is improved;
(2) The phase angle estimation value is updated by using the phase angle correction factor, so that the problem of large phase angle estimation error under low signal-to-noise ratio can be effectively solved, and the accuracy of phase noise estimation is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the structures shown in the drawings without creative efforts.
FIG. 1 is a schematic flow diagram of the present invention;
figure 2 is a graph comparing the performance of the present invention with the prior art implementation.
Detailed Description
It should be understood that the specific embodiments described herein are merely illustrative of the invention and do not: the present invention is defined.
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention. In addition, the technical solutions in the embodiments may be combined with each other, but must be based on the realization of the technical solutions by those skilled in the art, and when the technical solutions are contradictory or cannot be realized, the technical solutions should be considered that the combination does not exist, and the technical solutions are not within the protection scope of the present invention.
As shown in fig. 1, a method for estimating and compensating phase noise includes the following steps:
s1: determining an initial noise estimation value through the received pilot signal and a local pilot signal;
s2: further determining a phase angle estimation value according to the initial phase noise estimation value;
s3: traversing the phase angle estimation value to determine a flag (k, l), and simultaneously performing grouping calculation on the frequency domain data to obtain a phase angle correction factor;
s4: reversely updating the phase angle estimation value and the phase noise estimation value by using the phase angle correction factor;
s5: and performing phase noise compensation on the received signal by using the phase noise estimation value.
Assuming that 2 PTRS signals are received in total, each PTRS signal configures 2 REs, and the received frequency domain data of the PTRS is divided into one group, i.e., the frequency domain data of every two REs is a group.
(1) From the received PTRS signal Y ptrs (k, l) and local PTRS signal X ptrs (k, l) obtaining an estimated phase noise value H temp (k, l), k denotes a subcarrier index, and l denotes a signal index.
Figure BDA0004030100570000041
Wherein,
Figure BDA0004030100570000042
Figure BDA0004030100570000043
Figure BDA0004030100570000044
/>
(2) Based on the phase noise estimate H temp (k, l) calculating a phase angle estimate
Figure BDA0004030100570000051
Figure BDA0004030100570000052
(3) Ergodic phase angle estimation
Figure BDA0004030100570000053
If/or>
Figure BDA0004030100570000054
Then flag (k, l) =1; if/or>
Figure BDA0004030100570000055
Then flag (k, l) =0;
Figure BDA0004030100570000056
(4) The flag (k, l) is divided into 1 group according to the frequency domain configuration of PTRS, 2 flags (k, l) in each group
The following operation is performed.
Figure BDA0004030100570000057
Figure BDA0004030100570000058
Wherein, the value range of n is [1,2].
(5) Calculating a phase angle reliability factor beta
Figure BDA0004030100570000059
(6) Calculating a phase angle estimation value according to the phase angle reliability factor
Figure BDA00040301005700000510
Figure BDA00040301005700000511
Figure BDA00040301005700000512
(7) Obtaining a final estimation value H of the phase noise of each signal ptrs (l)
Figure BDA00040301005700000513
(8) Using estimated phase noise H ptrs (l) And performing phase noise compensation.
As shown in fig. 2, which is a graph comparing the effect before and after the implementation of the present embodiment, it can be seen from fig. 2 that the SNR point corresponding to the block error rate of 1 × e-03 after the implementation of the present embodiment is about 13.5dB, and the SNR point corresponding to the block error rate of 1 × e-03 before the implementation of the present embodiment is about 14.4dB, which is about 0.9dB higher than the performance after the implementation of the present embodiment before the implementation of the present embodiment, that is, the demodulation performance after the implementation of the present embodiment is improved.
The foregoing is merely a preferred embodiment of the invention, it being understood that the invention is not limited to the forms disclosed herein, but is not to be construed as excluding other embodiments, and: various other combinations, modifications, and environments and can be made within the scope of the concepts described herein, as indicated by the above teachings or by the skill or knowledge of the relevant art. And that modifications and variations may be effected by those skilled in the art without departing from the spirit and scope of the invention as defined by the appended claims.

Claims (7)

1. A phase noise estimation and compensation method is characterized by comprising the following specific steps:
s1: determining an initial noise estimation value through the received pilot signal and the local pilot signal;
s2: further determining a phase angle estimation value according to the initial phase noise estimation value;
s3: traversing the phase angle estimation value to determine a flag (k, l), and simultaneously performing grouping calculation on the frequency domain data to obtain a phase angle correction factor;
s4: reversely updating the phase angle estimation value and the phase noise estimation value by using the phase angle correction factor;
s5: and performing phase noise compensation on the received signal by using the phase noise estimation value.
2. The method for estimating and compensating phase noise according to claim 1, wherein the step S1 specifically comprises:
Figure FDA0004030100560000011
wherein, Y ptrs (k, l) is the received pilot signal, X ptrs (k, l) is a local pilot signal, H temp (k, l) is the initial estimate of phase noise, and l is in the range of [1]L is the number of pilot signals configured in the time domain, and k has a value range of [1, K ]]And K is the number of REs configured for each pilot signal.
3. The method for estimating and compensating phase noise according to claim 1, wherein the step S2 specifically comprises the following steps:
Figure FDA0004030100560000012
wherein,
Figure FDA0004030100560000013
is a phase angleEstimate of H temp (k, l) is the initial estimate of phase noise, and l has a value range of [1, L ]]L is the number of pilot signals configured in the time domain, and k has a value range of [1, K ]]And K is the number of REs configured for each pilot signal.
4. The method for estimating and compensating phase noise according to claim 1, wherein the step S3 of determining flag (k, l) by traversing the estimated phase angle value is specifically as follows: ergodic phase angle estimation
Figure FDA0004030100560000014
If/or>
Figure FDA0004030100560000015
Then flag (k, l) =1; if/or>
Figure FDA0004030100560000016
Then flag (k, l) =0.
5. The method for phase noise estimation and compensation according to claim 1, wherein the step S3 of grouping the frequency domain data specifically comprises: dividing the flag (k, l) corresponding to the ith signal into M groups according to the frequency domain configuration of the PTRS, and carrying out XOR operation on N flags (k, l) in each group:
Figure FDA0004030100560000017
wherein, F (m, l) represents the mth group phase noise estimation reliability flag of the ith pilot symbol, the range of N is [1, N ], the range of m is [1, m ], if N flag (k, l) values in each group are the same, F (m, l) =0, otherwise, F (m, l) =1.
6. The method according to claim 1, wherein the phase angle correction factor calculated in step S3 is specifically:
Figure FDA0004030100560000021
wherein, F (M, L) represents the mth group phase noise estimation reliability flag of the lth pilot symbol, M is an index, and has a value range of [1, M ], M is the frequency domain data grouping number, L is an index, and has a value range of [1, L ], L is the pilot signal number configured on the time domain, and β is a phase angle correction factor.
7. The method for estimating and compensating phase noise according to claim 1, wherein the step S4 specifically comprises:
calculating a phase angle estimate for each signal based on the phase angle correction factor
Figure FDA0004030100560000022
Figure FDA0004030100560000023
Wherein beta is a phase angle correction factor, K ptrs Is the Kth RE of the pilot signal, K is the index, l is the index, the value range is [1, L],H ptrs (k, l) is the phase noise estimate;
obtaining the final estimated value H of the phase noise ptrs (k,l):
Figure FDA0004030100560000024
/>
CN202211722484.1A 2022-12-30 2022-12-30 Phase noise estimation and compensation method Pending CN115955380A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211722484.1A CN115955380A (en) 2022-12-30 2022-12-30 Phase noise estimation and compensation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211722484.1A CN115955380A (en) 2022-12-30 2022-12-30 Phase noise estimation and compensation method

Publications (1)

Publication Number Publication Date
CN115955380A true CN115955380A (en) 2023-04-11

Family

ID=87282346

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211722484.1A Pending CN115955380A (en) 2022-12-30 2022-12-30 Phase noise estimation and compensation method

Country Status (1)

Country Link
CN (1) CN115955380A (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101394381A (en) * 2007-09-19 2009-03-25 中国科学院微电子研究所 Self-adaptive equalizer for phase noise compensation
US20140169434A1 (en) * 2012-12-14 2014-06-19 Thushara Hewavithana Method and apparatus for phase noise mitigation
CN111786925A (en) * 2015-09-10 2020-10-16 Lg 电子株式会社 Broadcast signal transmitting apparatus and method, and broadcast signal receiving apparatus and method
CN113472712A (en) * 2021-06-30 2021-10-01 中铁二院工程集团有限责任公司 Phase noise suppression method
CN114039829A (en) * 2021-11-09 2022-02-11 北京邮电大学 Phase recovery method and system based on pilot frequency assistance

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101394381A (en) * 2007-09-19 2009-03-25 中国科学院微电子研究所 Self-adaptive equalizer for phase noise compensation
US20140169434A1 (en) * 2012-12-14 2014-06-19 Thushara Hewavithana Method and apparatus for phase noise mitigation
CN111786925A (en) * 2015-09-10 2020-10-16 Lg 电子株式会社 Broadcast signal transmitting apparatus and method, and broadcast signal receiving apparatus and method
CN113472712A (en) * 2021-06-30 2021-10-01 中铁二院工程集团有限责任公司 Phase noise suppression method
CN114039829A (en) * 2021-11-09 2022-02-11 北京邮电大学 Phase recovery method and system based on pilot frequency assistance

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JEAN-CHRISTOPHE SIBEL: ""Tracking the phase noise in sub-THz bands"", 《2022 IEEE WIRELESS COMMUNICATIONS AND NETWORKING CONFERENCE (WCNC)》, 16 May 2022 (2022-05-16) *
吴继文: ""TDS-OFDM***相位噪声补偿技术研究与实现"", 《中国优秀硕士学位论文全文数据库信息科技辑》, 15 August 2009 (2009-08-15) *

Similar Documents

Publication Publication Date Title
CN103095615B (en) The method and device that a kind of phase noise is estimated, compensated
US7483497B2 (en) Method and apparatus for calculating log-likelihood ratio for decoding in a receiver for a mobile communication system
CN102307175A (en) Multicarrier system soft-decision method
US8532216B2 (en) MIMO detection method and system based on channel correlation
CN108683620B (en) Phase noise resistant blind carrier recovery method suitable for high-order modulation mode
CN109862545A (en) Frequency bias compensation method, device, computer equipment and the storage medium of Bluetooth signal
KR101461837B1 (en) Channel-adaptive likelihood determination
CN101120559A (en) Method, apparatus and computer program product for channel estimation of transmission channels with memory in digital telecommunications system by estimating the RMS delay spread
CN106452652B (en) A kind of MPI suppression method based on chaos wireless communication system
US9071401B2 (en) Pilot-less noise estimation
CN110830398B (en) Frequency domain average channel estimation method in symbol applied in optical fiber DMT system
CN115955380A (en) Phase noise estimation and compensation method
US7315587B2 (en) Demodulation method and apparatus based on differential detection system for π/4 shifted QPSK modulated wave
CN112039591A (en) Carrier phase estimation algorithm based on dichotomy
CN101599772B (en) Method and device for demodulating signal
CN113411279B (en) Reusable Q power demapping method and system based on DVB-S2 system
CN113794661B (en) Method and system for optimizing receiving performance based on phase noise
US7684522B2 (en) Method and system for determining a log-likelihood ratio (LLR) corresponding to each bit of a symbol
CN101577689B (en) Method for detecting and correcting residual frequency offset of national standard of digital terrestrial broadcasting
CN109981503B (en) Layered k-best spherical decoding method and device
CN113746577A (en) Subcarrier signal-to-noise ratio estimation method of OFDM power line communication system
CN107508777B (en) Channel estimation method based on enhanced self-adaptive polarization linear interpolation
CN115426234B (en) Millimeter wave FBMC system phase noise compensation method and application thereof
CN112260973B (en) Combined channel estimation method for power line carrier communication
CN101316143B (en) Signal-to-noise ratio estimation device, system and method based on star map measurement

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination