CN116016049B - 一种适合高动态低信噪比的频偏估计方法 - Google Patents
一种适合高动态低信噪比的频偏估计方法 Download PDFInfo
- Publication number
- CN116016049B CN116016049B CN202211688178.0A CN202211688178A CN116016049B CN 116016049 B CN116016049 B CN 116016049B CN 202211688178 A CN202211688178 A CN 202211688178A CN 116016049 B CN116016049 B CN 116016049B
- Authority
- CN
- China
- Prior art keywords
- phase
- frequency offset
- module
- cnt
- noise ratio
- 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 37
- 238000009825 accumulation Methods 0.000 claims abstract description 9
- 238000004891 communication Methods 0.000 claims abstract description 6
- 238000005311 autocorrelation function Methods 0.000 claims description 21
- 238000012937 correction Methods 0.000 claims description 13
- 238000013507 mapping Methods 0.000 claims description 10
- 239000000654 additive Substances 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 230000000996 additive effect Effects 0.000 claims description 4
- 230000003631 expected effect Effects 0.000 abstract description 3
- 238000013459 approach Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 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
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明属于通信技术领域,具体提供一种适合高动态低信噪比的频偏估计方法,用以解决传统时域频偏估计方法估计范围与估计精度无法统一的问题,实现低信噪比、高动态和高估计精度的频偏估计。本发明针对传统时域频偏估计方法出现的相位翻转情况进行补偿,使得时域频偏估计方法不再局限于估计范围,且在高信噪比下的估计精度逼近CRLB;并且,利用系统中的速度参数,对相位补偿进行判决,能够解决补偿过多的问题;同时,通过相位累加模块能够将多点能量累加在一起,在低信噪比条件下具有良好的估计性能,符合预期效果;综上所述,本发明提出的适合高动态低信噪比的频偏估计方法能够实现低信噪比、高动态和高估计精度的频偏估计。
Description
技术领域
本发明属于通信技术领域,具体提供一种适合高动态低信噪比的频偏估计方法。
背景技术
在通信中,收发端之间的相对运动会使得载波频率发生变化,产生预期之外的频偏差;因此,接收端为了能够正常工作,通常需要进行频偏估计;但由于通信的数据容量有限,在频偏估计中的估计精度、估计范围和计算复杂度等无法全部满足。现有的频偏估计技术分为频域估计和时域估计,频域估计方法的估计范围大,但存在着估计精度受FFT点数影响的问题,要想得到精确的估计,则需要不断增加导频长度;时域估计方法的估计精度可以达到很高,但会存在估计范围偏小的问题,如Fitz算法,该算法估计精度高,但估计范围很窄,无法适应高动态情况。
发明内容
本发明的目的在于针对上述现有技术的缺陷,提供一种适合高动态低信噪比的频偏估计方法;本发明能够解决传统时域频偏估计方法估计范围与估计精度无法统一的问题,实现低信噪比、高动态和高估计精度的频偏估计。
为实现上述目的,本发明采用的技术方案为:
一种适合高动态低信噪比的频偏估计方法,其特征在于,包括以下步骤:
S1:采样得到接收信号r(k):
其中,c(k)为调制符号,fd为载波频偏值,T为符号间隔,θ为初始相位,N(k)为离散化后的加性高斯白噪声,L为导频序列长度,j为虚数单位;
S2:对接收信号r(k)进行星座逆映射,得到消除导频调制信息后接收信号z(k):
S3:计算接收信号z(k)的自相关函数R(m),并对自相关函数R(m)求相位,得到自相关函数中每一点的相位A(m);
S4:将自相关函数中每n点相位进行累加,得到合并后相位A′(p):
其中,n为累加点数;
S5:对合并后相位进行补偿,具体过程为:
S5.1:给定初始矫正相位Q(1)=A′(1);
S5.2:计算相位不会发生翻转的最大相位点数cntmax,进而设置判决点cnt=cntmax-2;
S5.3:对合并后相位中判决点cnt及其前的每个点(即1≤p≤cnt)进行2πM的相位补偿,M为:
Q(p)=A′(p)+2πM
其中,Q(p)为矫正相位;
S5.4:计算判决点cnt及其前的所有点的频偏fbk:
若频偏fbk未超过系统最大多普勒频偏fdmax,则Q(1:p)不变;否则,矫正模块执行:Q(1:p)=A′(1:p);
S5.5:对合并后相位中判决点cnt后每个点(即cnt<p≤P)进行2πM的相位补偿;
S6:计算多普勒频偏估计值f:
进一步的,步骤S3中,自相关函数R(m)为:
其中,N′(m)为自相关后加性高斯白噪声,N为自相关函数长度。
进一步的,步骤S5.2中,最大相位点数cntmax具体为:
其中,fs为采样频率;fdmax为系统最大多普勒频偏。
进一步的,所述频偏估计方法应用于无线通信的接发收系统中,所述接发收系统包括:发射端、接收端与信道;发射端包括:星座映射模块与加载波发送模块,信道为高斯白噪声信道,接收端包括:去载波模块、星座逆映射模块、求相关模块、相位累加模块、相位补偿模块、判决模块与矫正模块;所述去载波模块执行步骤S1,所述星座逆映射模块执行步骤S2,所述求相关模块执行步骤S3,所述相位累加模块执行步骤S4,所述相位补偿模块执行步骤S5.1、S5.2、S5.3、S5.5,所述判决模块与矫正模块执行步骤S5.4。
基于上述技术方案,本发明的有益效果在于:
本发明提供一种适合高动态低信噪比的频偏估计方法,能够解决传统时域频偏估计方法估计范围与估计精度无法统一的问题;本发明针对传统时域频偏估计方法出现的相位翻转情况进行补偿,使得时域频偏估计方法不再局限于估计范围,且在高信噪比下的估计精度逼近CRLB。然而,由于在信噪比较低时会出现补偿过多的情况,使得计算出的频率极大或极小,导致在低信噪比情况下估计精度会出现门限的情况,门限以下性能会急速恶化;针对此问题,本发明利用系统中的速度参数,对相位补偿进行判决,能够解决补偿过多的问题;同时,通过相位累加模块能够将多点能量累加在一起,在低信噪比条件下具有良好的估计性能,符合预期效果;综上所述,本发明能够实现低信噪比、高动态和高估计精度的频偏估计。
附图说明
图1为本发明中适合高动态低信噪比的频偏估计方法的流程示意图。
图2为本发明实施例中相位过补偿与正确补偿的测试结果图。
图3为本发明实施例中时域算法自相关函数长度与估计误差的测试结果图。
具体实施方式
为使本发明的目的、技术方案与有益效果更加清楚明白,下面结合附图和实施例对本发明进行进一步的说明。
本实施例提供一种适合高动态低信噪比的频偏估计方法,所述频偏估计方法应用于无线通信的接发收系统中,如图1所示;所述接发收系统包括:发射端、接收端与信道;发射端包括:星座映射模块与加载波发送模块,信道为高斯白噪声信道,接收端包括:去载波模块、星座逆映射模块、相位累加模块、求相关模块、相位补偿模块、判决模块与矫正模块。
发射端上电时,执行以下步骤:
产生系统载波频率fc=60GHz,最大运动速度Vmax=20马赫;
接收端开始工作时,去载波模块执行以下步骤:
S1:采样后接收信号表示为:
其中,r(k)为接收信号,c(k)为调制符号,fd为载波频偏值,T为符号间隔(T=1/(10MHz)),θ为初始相位,N(k)为离散化后的加性高斯白噪声,L为导频序列长度,j为虚数单位;
针对接收信号,星座逆映射模块执行以下步骤:
S2:PSK调制系统中满足表示c(k)的共轭),对接收信号进行星座逆映射消除导频调制信息后信号为:
其中,z(k)为消除导频调制信息后的接收信号;
去调制信息后的信号输入求相关模块后,求相关模块执行以下步骤:
S3.1:计算消除导频调制信息后的接收信号的自相关函数为:
其中,R(m)为自相关函数,N为自相关函数长度;
上式可带入z(k)表达式求得自相关函数R(m)(包含频偏信息)的另一种表达式:
其中,N′(m)为相关后的加性高斯白噪声;
S3.2:对自相关函数求相位,得到自相关函数中每一点的相位:
A(m)=arg(R(m))
其中,arg为求相位符号;
求得相位信息后,相位累加模块执行以下步骤:
S4:将自相关函数中每n点相位进行累加,得到合并后的相位:
其中,n为累加点数,在高信噪比系统下令n=1;
对合并后的相位进行补偿,相位补偿模块执行以下步骤:
S5.1:给定初始矫正相位Q(1)=A′(1);
S5.2:计算相位不会发生翻转的最大相位点数cntmax,进而设置判决点cnt=cntmax-2;最大相位点数cntmax具体为:
其中,fs为采样频率;fdmax为系统最大多普勒频偏:fdmax=(Vmax×fc)/c,c为光速;
S5.3:对合并后的相位中判决点cnt及其前的每个点(即1≤p≤cnt)进行2πM的相位补偿,M的计算如下:
Q(p)=A′(p)+2πM
再对补偿进行判决,判决模块执行以下步骤:
S5.4:计算判决点cnt及其前的所有点的频偏fbk:
若频偏fbk未超过系统最大多普勒频偏fdmax,则Q(1:p)不变;否则,矫正模块执行:Q(1:p)=A′(1:p);Q(1:p)表示第1至第p点的矫正相位,A′(1:p)表示第1至第p点的合并后相位;
经判决后,相位补偿模块继续执行以下步骤:
S5.5:对合并后的相位中判决点cnt后每个点(即cnt<p≤P)进行2πM的相位补偿,M的计算过程与步骤S5.3相同;
完成所有点的相位补偿后,相位补偿模块执行以下步骤:
S6:计算系统的多普勒频偏估计f:
基于本实施例中频偏估计方法,能够解决时域频偏估计方法估计范围与估计精度无法统一的问题;本实施例针对其出现的相位翻转情况进行补偿,使得时域频偏估计方法不再局限于估计范围,且在高信噪比下的估计精度逼近CRLB;但在低信噪比情况下估计精度会出现门限的情况,门限以下性能会急速恶化,这是由于在信噪比较低时会出现补偿过多的情况,使得计算出的频率极大或极小;在此基础上,本实施例利用系统中的速度参数,对相位补偿进行判决,能够解决补偿过多的问题;同时,通过相位累加模块能够将多点能量累加在一起,在低信噪比条件下具有良好的估计性能,符合预期效果;本实施例中相位过补偿与正确补偿的测试结果图如图2所示,时域算法自相关函数长度与估计误差的测试结果图如图3所示;由图可见,本实施例能够实现低信噪比、高动态和高估计精度的频偏估计。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (2)
1.一种适合高动态低信噪比的频偏估计方法,其特征在于,包括以下步骤:
S1:采样得到接收信号r(k):
其中,c(k)为调制符号,fd为载波频偏值,T为符号间隔,θ为初始相位,N(k)为离散化后的加性高斯白噪声,L为导频序列长度,j为虚数单位;
S2:对接收信号r(k)进行星座逆映射,得到消除导频调制信息后接收信号z(k):
S3:计算接收信号z(k)的自相关函数R(m),并对自相关函数R(m)求相位,得到自相关函数中每一点的相位A(m);自相关函数R(m)为:
其中,N′(m)为自相关后加性高斯白噪声,N为自相关函数长度;
S4:将自相关函数中每n点相位进行累加,得到合并后相位A′(p):
其中,n为累加点数;
S5:对合并后相位进行补偿,具体过程为:
S5.1:给定初始矫正相位Q(1)=A′(1);
S5.2:计算相位不会发生翻转的最大相位点数cntmax,进而设置判决点cnt=cntmax-2;最大相位点数cntmax具体为:
其中,fs为采样频率;fdmax为系统最大多普勒频偏;
S5.3:对合并后相位中判决点cnt及其前的每个点(即1≤p≤cnt)进行2πM的相位补偿,M为:
Q(p)=A′(p)+2πM
其中,Q(p)为矫正相位;
S5.4:计算判决点cnt及其前的所有点的频偏fbk:
若频偏fbk未超过系统最大多普勒频偏fdmax,则Q(1:p)不变;否则,执行:Q(1:p)=A′(1:p);
S5.5:对合并后相位中判决点cnt后每个点(即cnt<p≤P)进行2πM的相位补偿;
S6:计算多普勒频偏估计值f:
2.按权利要求1所述适合高动态低信噪比的频偏估计方法,其特征在于,所述频偏估计方法应用于无线通信的接发收系统中,所述接发收系统包括:发射端、接收端与信道;发射端包括:星座映射模块与加载波发送模块,信道为高斯白噪声信道,接收端包括:去载波模块、星座逆映射模块、求相关模块、相位累加模块、相位补偿模块、判决模块与矫正模块;所述去载波模块执行步骤S1,所述星座逆映射模块执行步骤S2,所述求相关模块执行步骤S3,所述相位累加模块执行步骤S4,所述相位补偿模块执行步骤S5.1、S5.2、S5.3、S5.5,所述判决模块与矫正模块执行步骤S5.4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211688178.0A CN116016049B (zh) | 2022-12-27 | 2022-12-27 | 一种适合高动态低信噪比的频偏估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211688178.0A CN116016049B (zh) | 2022-12-27 | 2022-12-27 | 一种适合高动态低信噪比的频偏估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116016049A CN116016049A (zh) | 2023-04-25 |
CN116016049B true CN116016049B (zh) | 2024-04-19 |
Family
ID=86027787
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211688178.0A Active CN116016049B (zh) | 2022-12-27 | 2022-12-27 | 一种适合高动态低信噪比的频偏估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116016049B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003073717A1 (en) * | 2002-02-22 | 2003-09-04 | Intersil Americas Inc. | Method and apparatus for timing synchronization in a multicarrier system |
CN101404633A (zh) * | 2008-11-20 | 2009-04-08 | 北京韦加航通科技有限责任公司 | 基于块传输的单载波系统的载波跟踪方法 |
CN105681245A (zh) * | 2016-03-29 | 2016-06-15 | 成都国腾电子集团有限公司 | 一种基于连续分段滑窗估计频偏的频率跟踪方法 |
CN106789764A (zh) * | 2016-11-18 | 2017-05-31 | 杭州电子科技大学 | 联合加权阈值去噪与均衡判决的变换域二次估计方法 |
CN109633709A (zh) * | 2018-12-28 | 2019-04-16 | 四川安迪科技实业有限公司 | 一种卫星通信系统中实用高效的频偏估计方法 |
CN112039818A (zh) * | 2020-08-31 | 2020-12-04 | 中原工学院 | 基于频域移位平均周期图和神经网络译码的载波同步方法 |
CN113364716A (zh) * | 2021-06-04 | 2021-09-07 | 深圳智微电子科技有限公司 | 一种ofdm系统中采样频偏的相位补偿方法 |
CN115242587A (zh) * | 2022-07-27 | 2022-10-25 | 西安电子科技大学 | 一种低信噪比环境下基于数据辅助的载波频偏估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10382244B2 (en) * | 2017-02-14 | 2019-08-13 | Samsung Electronics Co., Ltd. | System and method for providing time offset and frequency offset estimation for vehicle to everything communication system |
-
2022
- 2022-12-27 CN CN202211688178.0A patent/CN116016049B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003073717A1 (en) * | 2002-02-22 | 2003-09-04 | Intersil Americas Inc. | Method and apparatus for timing synchronization in a multicarrier system |
CN101404633A (zh) * | 2008-11-20 | 2009-04-08 | 北京韦加航通科技有限责任公司 | 基于块传输的单载波系统的载波跟踪方法 |
CN105681245A (zh) * | 2016-03-29 | 2016-06-15 | 成都国腾电子集团有限公司 | 一种基于连续分段滑窗估计频偏的频率跟踪方法 |
CN106789764A (zh) * | 2016-11-18 | 2017-05-31 | 杭州电子科技大学 | 联合加权阈值去噪与均衡判决的变换域二次估计方法 |
CN109633709A (zh) * | 2018-12-28 | 2019-04-16 | 四川安迪科技实业有限公司 | 一种卫星通信系统中实用高效的频偏估计方法 |
CN112039818A (zh) * | 2020-08-31 | 2020-12-04 | 中原工学院 | 基于频域移位平均周期图和神经网络译码的载波同步方法 |
CN113364716A (zh) * | 2021-06-04 | 2021-09-07 | 深圳智微电子科技有限公司 | 一种ofdm系统中采样频偏的相位补偿方法 |
CN115242587A (zh) * | 2022-07-27 | 2022-10-25 | 西安电子科技大学 | 一种低信噪比环境下基于数据辅助的载波频偏估计方法 |
Non-Patent Citations (3)
Title |
---|
I.M. Ngebani ; .Analysis and Optimization of Auto-Correlation Based Frequency Offset Estimation.SAIEE Africa Research Journal.2015,全文. * |
LTE中一种改进的基于探测参考信号的定时估计算法;田浩;杨霖;李少谦;;电讯技术;20131120(11);全文 * |
基于5GNR的LEO卫星通信大尺度时变频偏估计算法研究;田丁;信息科技辑;20210515;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116016049A (zh) | 2023-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11108612B2 (en) | Anti-interference signal detection and synchronization method for wireless broadband communication system | |
US10778496B2 (en) | OFDM system synchronization tracking method and apparatus | |
JP3214159B2 (ja) | キャリア検出器 | |
CN110518936B (zh) | 高超声速飞行器大动态多普勒快速捕获方法及通信系统 | |
CN105245303A (zh) | 一种高速突发解调同步系统 | |
CN104852875A (zh) | 高动态大频偏突发信号的频偏估计方法 | |
JP2006352746A (ja) | 直交周波数分割多重伝送用受信機 | |
CN109687912B (zh) | 相干光ofdm系统中一种时域容积卡尔曼相位噪声补偿方法 | |
CN111884973B (zh) | 一种用于单载波频域均衡系统接收端的数据接收方法 | |
CN110943758A (zh) | 基于多普勒频率变化率和频率联合搜索的二次捕获方法 | |
CN107342960A (zh) | 一种适合幅度相移键控的非数据辅助频偏估计方法 | |
CN111131098B (zh) | 基于水声ofdm信号自相关输出峰插值的多普勒因子估计方法 | |
CN103117980A (zh) | 用于正交频分复用接收机的快速数字自动频率控制方法 | |
CN112787730A (zh) | 基于联合波形设计的多普勒估计方法 | |
CN106656454B (zh) | 基于宽带跳频系统的频偏估计与补偿方法 | |
CN111654308B (zh) | 一种突发扩频弱信号的精密频偏估计方法 | |
CN109756435B (zh) | 一种对信号的频偏估计方法 | |
CN116016049B (zh) | 一种适合高动态低信噪比的频偏估计方法 | |
CN111832708B (zh) | 一种基于长短期记忆网络的窄带信号时差估计系统 | |
CN113115430A (zh) | 一种高速突发数字解调系统 | |
CN115643139A (zh) | 一种抗频偏的突发信号检测系统及检测方法 | |
US7428272B2 (en) | Method and apparatus for joint phase offset and frequency offset estimator for MPSK transmission | |
CN100355255C (zh) | 基于统计平均的大搜索范围ofdm系统载波同步方法 | |
CN110943949B (zh) | 一种基于sc-fde传输体制的宽带角跟踪方法 | |
CN114285713A (zh) | 一种低轨宽带卫星时频偏估计方法及系统 |
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 |