CN115955380B - Phase noise estimation and compensation method - Google Patents
Phase noise estimation and compensation method Download PDFInfo
- Publication number
- CN115955380B CN115955380B CN202211722484.1A CN202211722484A CN115955380B CN 115955380 B CN115955380 B CN 115955380B CN 202211722484 A CN202211722484 A CN 202211722484A CN 115955380 B CN115955380 B CN 115955380B
- Authority
- CN
- China
- Prior art keywords
- phase angle
- phase noise
- phase
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000012937 correction Methods 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 238000004891 communication Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000003780 insertion Methods 0.000 description 2
- 230000037431 insertion Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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 is characterized in that a phase noise initial estimation value is determined 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 carrying out grouping calculation on the frequency domain data to obtain a phase angle correction factor; reversely updating the phase angle estimated value and the phase noise estimated value by using the phase angle correction factor; and performing phase noise compensation on the received signal by using the phase noise estimated value. According to the phase angle correction factor calculation 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 correction factor is used for updating the phase angle estimation value, so that the problem of large phase angle estimation error under the condition of low signal-to-noise ratio can be effectively solved, and the accuracy of phase noise estimation is improved.
Description
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 rates and spectral efficiency of communication systems continue to increase, communication systems will employ higher order quadrature amplitude modulation (Quadrature Amplitude Modulation, QAM) and greater system bandwidth to meet these requirements. Higher order QAM modulated signals are very sensitive to phase noise due to the increased degree of constellation point density, and the lower intensity of phase noise can severely degrade system performance. In addition, due to the trend of the communication system to have a higher operating frequency band and a larger system bandwidth, the effect of phase noise on the system will become more serious. How to suppress the influence of phase noise is critical to the performance of the high-order quadrature amplitude modulation communication system.
The prior art mainly suppresses the influence of phase noise on the system performance by estimating and compensating the phase noise. While the estimation compensation method can be generally classified into a pilot data-aided method and a decision data-directed method. The performance of the pilot data auxiliary method is closely related to the density of pilot insertion, the number of the inserted pilot and the like, and the increase of the pilot insertion density and the number can effectively improve the accuracy of phase noise estimation, but greatly reduce the transmission efficiency of the system. The method for judging data guiding can effectively improve the transmission efficiency of the system, but the error propagation problem of the method can seriously influence the performance of a high-order QAM system.
Patent CN112688891a provides a method for phase noise estimation and compensation using PTRS signals in the time domain. Patent CN110460385B provides a method of estimating phase noise from pilot signals, which has no targeted improvement effect on the increase of phase noise estimation error due to 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 initially estimated phase noise.
The aim of the invention is realized by the following technical scheme:
Assume that L pilot signals are configured in total in the time domain, where the pilot signals are PTRS, K REs are configured for each PTRS signal, and frequency domain data of the PTRS on each signal is divided into M groups, that is, 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 a phase noise initial 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 carrying out grouping calculation on the frequency domain data to obtain a phase angle correction factor;
S4: reversely updating the phase angle estimated value and the phase noise estimated value by using the phase angle correction factor;
s5: and performing phase noise compensation on the received signal by using the phase noise estimated value.
Further, the specific calculation process in step S1 is as follows:
Wherein Y ptrs( K, L) is the received pilot signal, X ptrs (K, L) is the local pilot signal, H temp (K, L) is the initial phase noise estimation, L is the range of [1, L ], L is the number of pilot signals configured in the time domain, K is the range of [1, K ], and K is the number of REs configured for each pilot signal.
Further, the specific calculation process in step S2 is as follows:
wherein, For the phase angle estimation value, H temp (K, L) is the initial phase noise estimation value, the value range of L is [1, L ], L is the number of pilot signals configured on the time domain, the value range of K is [1, K ], and K is the number of REs configured on each pilot signal.
Further, the step S3 of traversing the phase angle estimation value to determine the flag (k, l) specifically includes: traversing phase angle estimatesIf it isFlag (k, l) =1; if it isFlag (k, l) =0.
Further, in the step S3, the grouping of the frequency domain data is specifically: the corresponding flag (k, l) of the first signal is divided into M groups according to the frequency domain configuration of the PTRS, and the N flags (k, l) in each group are subjected to exclusive OR operation:
Wherein F (m, l) represents the mth group of phase noise estimation reliability flags of the mth pilot symbol, the range of N is [1, N ], the range of m is [1, m ], if the values of N flags (k, l) in each group are the same, F (m, l) =0, otherwise F (m, l) =1.
Further, the phase angle correction factor calculated in the step S3 is specifically:
wherein F (M, L) represents the M-th group of phase noise estimation reliability marks of the first pilot symbol, M is an index, the value range is [1, M ], M is the number of frequency domain data packets, L is an index, the value range is [1, L ], L is the number of pilot signals configured on the time domain, and beta is a phase angle correction factor.
Further, the step S4 specifically includes:
Calculating a phase angle estimate for each signal based on the phase angle correction factor
Wherein, beta is a phase angle correction factor, K ptrs is the Kth RE of the pilot signal, K is an index, l is an index, the value range is [1, L ], and H ptrs (K, l) is a phase noise estimated value;
Obtaining a final estimated value H ptrs (k, l) of the phase noise:
the invention has the beneficial effects that:
(1) The phase noise correction factors are calculated in a mode of grouping the frequency domain data, so that the reliability of the phase angle correction factors is improved;
(2) The phase angle correction factor is used for updating the phase angle estimation value, so that the problem of large phase angle estimation error under the condition of 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 that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to the structures shown in these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of the present invention;
FIG. 2 is a graph comparing the performance of the present invention with prior art implementations.
Detailed Description
It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to be limiting. The present invention is defined.
The following description of the embodiments of the present invention will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention. In addition, the technical solutions of the embodiments may be combined with each other, but it is necessary to base that the technical solutions can be realized 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 is not within the scope of protection claimed by the present invention.
As shown in fig. 1, a phase noise estimation and compensation method includes the following specific steps:
s1: determining a phase noise initial 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 carrying out grouping calculation on the frequency domain data to obtain a phase angle correction factor;
S4: reversely updating the phase angle estimated value and the phase noise estimated value by using the phase angle correction factor;
s5: and performing phase noise compensation on the received signal by using the phase noise estimated value.
Assuming that 2 PTRS signals are received in total, each PTRS signal is configured with 2 REs, the received frequency domain data of PTRS is divided into one group, i.e., the frequency domain data of every two REs is one group.
(1) Phase noise estimation H temp (k, l) is obtained from the received PTRS signal Y ptrs (k, l) and the local PTRS signal X ptrs (k, l), k represents a subcarrier index, and l represents a signal index.
Wherein,
(2) Calculating phase angle estimation from phase noise estimation H temp (k, l)
(3) Traversing phase angle estimatesIf it isFlag (k, l) =1; if it is
Flag (k, l) =0;
(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.
Wherein, the value range of n is [1,2].
(5) Calculating the phase angle reliability factor beta
(6) Calculating a phase angle estimation value according to the phase angle reliability factor
(7) Obtaining a final estimate H ptrs (l) of the phase noise of each signal
(8) And performing phase noise compensation by using the estimated phase noise H ptrs (l).
As shown in fig. 2, in order to compare the effect graphs before and after implementing the present embodiment, it can be seen from fig. 2 that the SNR point corresponding to the block error rate 1*e-03 after implementing the present embodiment is about 13.5dB, and the SNR point corresponding to the block error rate 1*e-03 before implementing the present embodiment is about 14.4dB, which is about 0.9dB higher than the performance after implementing the present embodiment, that is, the demodulation performance corresponding to the present embodiment is improved.
The foregoing is merely a preferred embodiment of the present invention, and it is to be understood that the invention is not limited to the forms disclosed herein and is not to be construed as excluding other embodiments, but rather: various other combinations, modifications, and environments are possible within the scope of the concepts described herein, through the above teachings or through other variations or modifications within the skill or knowledge of the relevant art. And that modifications and variations which do not depart from the spirit and scope of the invention are intended to be within the scope of the appended claims.
Claims (1)
1. The phase noise estimation and compensation method is characterized by comprising the following specific steps:
S1: the initial phase noise estimation value is determined through the received pilot signal and the local pilot signal, and the specific calculation process is as follows:
Wherein Y ptrs (K, L) is the received pilot signal, X ptrs (K, L) is the local pilot signal, H temp (K, L) is the initial phase noise estimation value, the value range of L is [1, L ], L is the number of pilot signals configured in the time domain, the value range of K is [1, K ], and K is the number of REs configured for each pilot signal;
S2: the phase angle estimation value is further determined according to the initial phase noise estimation value, and the specific calculation process is as follows:
wherein, H temp (K, L) is an initial phase noise estimation value, the value range of L is [1, L ], L is the number of pilot signals configured on the time domain, the value range of K is [1, K ], and K is the number of REs configured on each pilot signal;
s3: traversing the phase angle estimation value to determine a flag (k, l), and simultaneously carrying out grouping calculation on the frequency domain data to obtain a phase angle correction factor;
The step of determining the flag (k, l) by the traversal phase angle estimation value is specifically as follows: traversing phase angle estimates If it isFlag (k, l) =1; if it isFlag (k, l) =0;
The grouping of the frequency domain data specifically includes: the corresponding flag (k, l) of the first signal is divided into M groups according to the frequency domain configuration of the PTRS, and the N flags (k, l) in each group are subjected to exclusive OR operation:
Wherein F (m, l) represents an mth group of phase noise estimation reliability flags of the mth 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 identical, F (m, l) =0, otherwise F (m, l) =1;
the calculated phase angle correction factor is specifically:
Wherein F (M, L) represents an mth group phase noise estimation reliability mark of the ith pilot symbol, M is an index, the value range is [1, M ], M is the number of frequency domain data packets, L is the index, the value range is [1, L ], L is the number of pilot signals configured on the time domain, and beta is a phase angle correction factor;
S4: the phase angle correction factor is used for reversely updating the phase angle estimated value and the phase noise estimated value, specifically:
Calculating a phase angle estimate for each signal based on the phase angle correction factor
Wherein, beta is a phase angle correction factor, K ptrs is the Kth RE of the pilot signal, K is an index, the value range is [1, K ], and K is the RE quantity configured by each pilot signal; l is an index, the value range is [1, L ], and L is the number of pilot signals configured on the time domain; h ptrs (k, l) is a phase noise estimate;
Obtaining a final estimated value H ptrs (k, l) of the phase noise:
s5: and performing phase noise compensation on the received signal by using the phase noise estimated value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211722484.1A CN115955380B (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 CN115955380B (en) | 2022-12-30 | 2022-12-30 | Phase noise estimation and compensation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115955380A CN115955380A (en) | 2023-04-11 |
| CN115955380B true CN115955380B (en) | 2024-08-27 |
Family
ID=87282346
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211722484.1A Active CN115955380B (en) | 2022-12-30 | 2022-12-30 | Phase noise estimation and compensation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115955380B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101394381A (en) * | 2007-09-19 | 2009-03-25 | 中国科学院微电子研究所 | Self-adaptive equalizer for phase noise compensation |
| CN111786925A (en) * | 2015-09-10 | 2020-10-16 | Lg 电子株式会社 | Broadcast signal transmitting apparatus and method and broadcast signal receiving apparatus and method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8897412B2 (en) * | 2012-12-14 | 2014-11-25 | Intel Corporation | Method and apparatus for phase noise mitigation |
| CN113472712B (en) * | 2021-06-30 | 2023-05-19 | 中铁二院工程集团有限责任公司 | Phase noise suppression method |
| CN114039829B (en) * | 2021-11-09 | 2022-08-09 | 北京邮电大学 | Phase recovery method and system based on pilot frequency assistance |
-
2022
- 2022-12-30 CN CN202211722484.1A patent/CN115955380B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101394381A (en) * | 2007-09-19 | 2009-03-25 | 中国科学院微电子研究所 | Self-adaptive equalizer for phase noise compensation |
| CN111786925A (en) * | 2015-09-10 | 2020-10-16 | Lg 电子株式会社 | Broadcast signal transmitting apparatus and method and broadcast signal receiving apparatus and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115955380A (en) | 2023-04-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100553544B1 (en) | Frequency Offset Estimation and Channel Equalization in Burst Orthogonal Frequency Division Multiplexing Systems | |
| KR100794979B1 (en) | Method and system for determining signal-to-interference ratio with reduced bias effect | |
| CN108683620B (en) | An Anti-Phase Noise Blind Carrier Recovery Method Suitable for Higher Order Modulation | |
| WO2001020831A1 (en) | Ofdm communication device and detecting method | |
| CN115882923B (en) | Symbol-level residual phase deviation compensation method for broadband satellite communication | |
| US20040203465A1 (en) | Mode adaptation in wireless systems | |
| US8422354B2 (en) | Channel estimation method and apparatus using data channel | |
| US20170324595A1 (en) | Circuits and methods for frequency offset estimation in fsk communications | |
| US8223825B2 (en) | Method and computer program for estimating signal quality value, and signal quality estimator, receiver and communication apparatus | |
| CN115955380B (en) | Phase noise estimation and compensation method | |
| CN117040976A (en) | Blind equalizer under non-cooperative communication condition and dual-mode joint blind equalization method | |
| CN103188176B (en) | The method of single-frequency or arrowband interference effect is reduced in ofdm system | |
| CN110830398B (en) | An Intra-Symbol Frequency Domain Average Channel Estimation Method Applied in Optical Fiber DMT Systems | |
| US20130294545A1 (en) | Method and apparatus for eliminating direct current offset | |
| CN119544081A (en) | An adaptive CKF phase noise compensation method in coherent optical OFDM systems | |
| CN114900267B (en) | Optimization method for eliminating abrupt change of frequency domain data symbol power by NR-5G | |
| US7864905B2 (en) | Interference alleviation equalizing apparatus of multi-carrier communication system and method thereof | |
| CN115913851B (en) | Carrier phase estimation method based on cubic spline interpolation | |
| CN114389710A (en) | A Novel CKF Phase Noise Compensation Method in Coherent Optical OFDM Systems | |
| CN112260973B (en) | Combined channel estimation method for power line carrier communication | |
| CN113746577A (en) | Subcarrier signal-to-noise ratio estimation method of OFDM power line communication system | |
| CN113794661A (en) | Method and system for optimizing receiving performance based on phase noise | |
| KR102627059B1 (en) | Method and apparatus for estimating channel based on phase difference | |
| CN101132389B (en) | A Compensation Method for Signal-to-Interference Ratio Estimation | |
| CN121173639B (en) | Signal processing method, device and receiver |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |