CN114760178B - 一种基于多径能量的频偏估计方法 - Google Patents
一种基于多径能量的频偏估计方法 Download PDFInfo
- Publication number
- CN114760178B CN114760178B CN202210380713.XA CN202210380713A CN114760178B CN 114760178 B CN114760178 B CN 114760178B CN 202210380713 A CN202210380713 A CN 202210380713A CN 114760178 B CN114760178 B CN 114760178B
- Authority
- CN
- China
- Prior art keywords
- frequency offset
- offset estimation
- points
- multipath
- point
- 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 60
- 238000005070 sampling Methods 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000001360 synchronised effect Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2657—Carrier synchronisation
- H04L27/2659—Coarse or integer frequency offset determination and synchronisation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2657—Carrier synchronisation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2655—Synchronisation arrangements
- H04L27/2657—Carrier synchronisation
- H04L27/266—Fine or fractional frequency offset determination and synchronisation
-
- 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
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Noise Elimination (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明涉及一种基于多径能量的频偏估计方法,属于频偏估计技术领域。本发明在进行频偏补偿后再进行相关求和,计算窗口内所有多径的能量值。通过上述类似的过程获得五个特定频点的多径能量和,然后再通过这五个点基于插值的LJEA的方法计算最大值所在的位置,即得到频偏估计结果。本发明提升了在低信噪比下的稳定性,提升了估计精度;在多径的条件下提升了信息的利用率。
Description
技术领域
本发明属于频偏估计技术领域,具体涉及一种基于多径能量的频偏估计方法。
背景技术
随着信息传输速率的提升,以及对时延和传输可靠性的要求提高。同步技术是无线信号解调的基础,其中载波同步技术(频偏估计)在低信噪比、选择性衰落信道下进行稳定、高精度的频偏估计可以有效提升数据解调灵敏度。如果存在频偏如图1所示,会造成误码率的提升,降低了解调灵敏度。
传统的单载波频偏估计方法大多基于特殊同步头设计的差分方法,如图2所示典型的同步头结构,特别是在低信噪比条件下,无法利用数据部分进行频偏估计。如果能够利用数据部分的有效信息可以进一步的提升频偏估计精度。
如果利用数据部分的信息进行频偏估计需要使用基于频域的频偏估计方法,在盲信号解调中常用此方法,但是往往只能选择一条最强径进行频偏估计。而其他的多径中也包含着频偏的信息,因此本发明通过多径能量和与频偏的关系进行频偏估计。同时,可以通过LJEA频域插值技术降低复杂度。
在实际的应用系统中会存在多普勒或者本地晶振频偏,特别是跳频通信中,需要准确的频偏估计。在突发信号传输系统中,信道环境复杂,由于实际应单位时间内跳数和带宽不断增加,那么对解调性能的提升或者维持提出了更高的要求。其中,频偏估计就是一个重要的过程,由于同步头比较短,如果只使用传统的差分方法无法提升性能。那么可以利用数据部分信息进行进一步的频偏修正。进而,能够在较低的信噪比下达到较高的灵敏度。
如果利用数据部分进行频偏估计,就破坏了同步头格式,无法进行差分估计,那么可以通过频率谱线估计方法。而传统通过频率谱线的估计方式是选择最强径进行频偏估计,如图3中的”A”径。本文利用窗口中”A~F”多径的能量之和进行频偏估计,进而提升低信噪比下的估计精度。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是如何提供一种基于多径能量的频偏估计方法,以解决通过选择最强径进行频偏估计在低信噪比下精度不高的问题。
(二)技术方案
为了解决上述技术问题,本发明提供了一种基于多径能量的频偏估计方法,该方法包括:假设经过经过粗频偏估计后,剩余频偏小于预设值,则在进行频偏补偿后再进行相关求和,计算窗口内所有多径的能量值,获得五个特定频点的多径能量和,然后再通过这五个点基于插值的LJEA的方法计算最大值所在的位置,即得到频偏估计结果。
进一步地,预设值为FFT变换的分辨率。
进一步地,分辨率为数据采样率/采用数据的个数。
进一步地,通过这五个点基于插值的LJEA的方法计算最大值所在的位置,即得到频偏估计结果包括:依据归一化频域位置取0和±1,±1/2五个点的幅度ψ(0),ψ(±1),ψ(±1/2),通过LJEA方法对频偏Δf=fs·Δξ进行计算,fs表示采样率,Δξ具体计算公式如下:
其中,|ψ(Ns-1)|表示傅里叶变换后的第Ns-1个点的幅度值,表示傅里叶变换后的第/>个点的幅度值,Ns表示采用傅里叶变换的点数。
本发明还提供一种基于多径能量的频偏估计方法,该方法包括如下步骤:
第一步,假设采用计算的码片为N,每个码片一个采样点,那么接收的信号表示为:x(k)=s(k)e2πjΔfk+n(k),其中s(k)为信息码片序列,n(k)为接收的噪声,e2πjΔfk表示发生频偏的载波,Δf表示频偏量,k表示第k个采样点;假设信号的频偏小于预设值,则通过信道到达检测算法获取信号的同步头位置,然后获得带有频偏Δf的信号x(k);
第二步,对x(k)进行频偏补偿,假设补偿的频偏为f,通过载波进行补偿,计算相关求和结果/> 其中/>表示第n个起始点的采样序列,对应的就是第n个多径采样序列;其中,fs表示采样率,N表示参与计算的所有采样点的总个数;
第三步,利用第二步中的幅度值计算所有n个多径的能量和,即:其中width表示窗口宽度,表示选择计算总的多径个数;
第四步,选择f分别为频点补偿位置计算R(f),通过这五个点基于插值的LJEA的方法计算最大值所在的位置,即得到频偏估计结果。
进一步地,预设值为FFT变换的分辨率。
进一步地,分辨率为fs表示采样率,N表示参与计算的所有采样点的总个数。
进一步地,频点补偿位置对应归一化频域位置取0和±1,±1/2的位置,即幅度最高点、幅度最低点和幅度为最高点1/2的位置。
进一步地,依据归一化频域位置取0和±1,±1/2五个点的幅度ψ(0),ψ(±1),ψ(±1/2),通过LJEA方法对频偏Δf=fs·Δξ进行计算,Δξ具体计算公式如下:
其中,|ψ(Ns-1)|表示傅里叶变换后的第Ns-1个点的幅度值,表示傅里叶变换后的第/>个点的幅度值,Ns表示采用傅里叶变换的点数。
进一步地,Δξ计算公式通过傅里叶变换公式在零点的泰勒展开,留取低阶项化简得到插值结果。
(三)有益效果
本发明提供一种基于多径能量的频偏估计方法,本发明提出的改进方法具有两个效果:1)提升了在低信噪比下的稳定性,提升了估计精度;2)在多径的条件下提升了信息的利用率。
本发明基于多径条件下频域的频偏估计做了进一步改进,能够在低信噪比下获得更加准确的频偏估计结果。利用多径能量求和的方法,在低信噪比情况下使得频偏估计更加准确。如图4所示,当信噪比小于2dB时本文的方法优于基于差分的方法,如果再利用捕获部分的数据将会有更好的效果。同时在整个信噪比区间内,基于多径能量和的频偏估计方法要优于基于最强径的频偏估计方法。同时计算量没有明显的增加,在可接受的范围内。
附图说明
图1为存在频偏情况下的单载波信号示意图;
图2为同步数据帧结构示意图;
图3为多径示意图;
图4为估计精度对比结果;
图5为通过LJEA方法频偏估计方法计算说明示意图;
图6为基于多径能量的频偏估计方法示意。
具体实施方式
为使本发明的目的、内容、和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
本发明的基于多径能量的频偏估计方法,包括:
假设经过经过粗频偏估计后,剩余频偏小于预设值,即小于FFT变换的分辨率(数据采样率/采用数据的个数)。进行频偏补偿后再进行相关求和,计算窗口内所有多径的能量值。通过上述类似的过程获得五个特定频点的多径能量和,然后再通过这五个点基于插值的LJEA的方法计算最大值所在的位置,即得到频偏估计结果。具体过程如下所示。
实施例1:
本发明的基于多径能量的频偏估计方法,包括:
假设采用计算的码片为N,每个码片一个采样点,那么接收的信号可以表示为:其中s(k)为信息码片序列,n(k)为接收的噪声,e2πjΔfk表示发生频偏的载波,Δf表示频偏量,单位是Hz,k表示第k个采样点。假设信号的频偏小于预设值,预设值为/>fs表示采样率,N表示参与计算的所有采样点的总个数。
第一步,通过信道到达检测算法获取信号的同步头位置,然后获得带有频偏Δf的信号x(k);
第二步,对x(k)进行频偏补偿,假设补偿的频偏为f,通过载波进行补偿,计算相关求和结果/> 其中/>表示第n个起始点的采样序列,对应的就是第n个多径采样序列;
第三步,利用第二步中的幅度值计算所有n个多径的能量和,即:其中width表示窗口宽度,表示选择计算总的多径个数;
第四步,如上图5所示,f分别为频点补偿位置计算R(f),对应图5中归一化频域位置取0和±1,±1/2五个点(幅度最高点、幅度最低点和幅度为最高点1/2的位置)的幅度ψ(0),ψ(±1),ψ(±1/2)。通过LJEA方法对频偏Δf=fs·Δξ进行计算,Δξ具体计算公式如下:
其中,|ψ(Ns-1)|表示傅里叶变换后的第Ns-1个点的幅度值,表示傅里叶变换后的第/>个点的幅度值,Ns表示采用傅里叶变换的点数。
上述公式通过傅里叶变换公式在零点的泰勒展开,留取低阶项化简得到上述插值结果。这是一个成熟的方法具体参考相应的文献。
本发明基于多径条件下频域的频偏估计做了进一步改进,能够在低信噪比下获得更加准确的频偏估计结果。利用多径能量求和的方法,在低信噪比情况下使得频偏估计更加准确。如图4所示,当信噪比小于2dB时本文的方法优于基于差分的方法,如果再利用捕获部分的数据将会有更好的效果。同时在整个信噪比区间内,基于多径能量和的频偏估计方法要优于基于最强径的频偏估计方法。同时计算量没有明显的增加,在可接受的范围内。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (9)
1.一种基于多径能量的频偏估计方法,其特征在于,该方法包括:假设经过粗频偏估计后,剩余频偏小于预设值,则在进行频偏补偿后再进行相关求和,计算窗口内所有多径的能量值,获得五个频点的多径能量和,然后再通过这五个点基于插值的低复杂度联合频偏估计LJEA的方法计算最大值所在的位置,即得到频偏估计结果;
其中,
通过这五个点基于插值的LJEA的方法计算最大值所在的位置,即得到频偏估计结果包括:依据归一化频域位置取0和±1,±1/2五个点的幅度ψ(0),ψ(±1),ψ(±1/2),通过LJEA方法对频偏Δf=fs·Δξ进行计算,fs表示采样率,Δξ具体计算公式如下:
其中,|ψ(Ns-1)|表示傅里叶变换后的第Ns-1个点的幅度值,表示傅里叶变换后的第/>个点的幅度值,Ns表示采用傅里叶变换的点数。
2.如权利要求1所述的基于多径能量的频偏估计方法,其特征在于,预设值为FFT变换的分辨率。
3.如权利要求2所述的基于多径能量的频偏估计方法,其特征在于,分辨率为数据采样率/采用数据的个数。
4.一种基于多径能量的频偏估计方法,其特征在于,该方法包括如下步骤:
第一步,假设采用计算的码片为N,每个码片一个采样点,那么接收的信号表示为:x(k)=s(k)e2πjΔfk+n(k),其中s(k)为信息码片序列,n(k)为接收的噪声,表示发生频偏的载波,Δf表示频偏量,k表示第k个采样点;假设信号的频偏小于预设值,则通过信道到达检测算法获取信号的同步头位置,然后获得带有频偏Δf的信号x(k);
第二步,对x(k)进行频偏补偿,假设补偿的频偏为f,通过载波进行补偿,计算相关求和结果/> 其中/>表示第n个起始点的采样序列,对应的就是第n个多径采样序列;其中,fs表示采样率,N表示参与计算的所有采样点的总个数;
第三步,利用第二步中的幅度值计算所有n个多径的能量和,即:其中width表示窗口宽度,表示选择计算总的多径个数;
第四步,选择f分别为 频点补偿位置计算R(f),通过这五个点基于插值的低复杂度联合频偏估计LJEA的方法计算最大值所在的位置,即得到频偏估计结果。
5.如权利要求4所述的基于多径能量的频偏估计方法,其特征在于,预设值为FFT变换的分辨率。
6.如权利要求5所述的基于多径能量的频偏估计方法,其特征在于,分辨率为fs表示采样率,N表示参与计算的所有采样点的总个数。
7.如权利要求4-6任一项所述的基于多径能量的频偏估计方法,其特征在于,频点补偿位置对应归一化频域位置取0和±1,±1/2的位置,即幅度最高点、幅度最低点和幅度为最高点1/2的位置。
8.权利要求7所述的基于多径能量的频偏估计方法,其特征在于,依据归一化频域位置取0和±1,±1/2五个点的幅度ψ(0),ψ(±1),ψ(±1/2),通过LJEA方法对频偏Δf=fs·Δξ进行计算,Δξ具体计算公式如下:
其中,|ψ(Ns-1)|表示傅里叶变换后的第Ns-1个点的幅度值,表示傅里叶变换后的第/>个点的幅度值,Ns表示采用傅里叶变换的点数。
9.如权利要求8所述的基于多径能量的频偏估计方法,其特征在于,Δξ计算公式通过傅里叶变换公式在零点的泰勒展开,留取低阶项化简得到插值结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210380713.XA CN114760178B (zh) | 2022-04-08 | 2022-04-08 | 一种基于多径能量的频偏估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210380713.XA CN114760178B (zh) | 2022-04-08 | 2022-04-08 | 一种基于多径能量的频偏估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114760178A CN114760178A (zh) | 2022-07-15 |
CN114760178B true CN114760178B (zh) | 2024-03-15 |
Family
ID=82329918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210380713.XA Active CN114760178B (zh) | 2022-04-08 | 2022-04-08 | 一种基于多径能量的频偏估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114760178B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1972137A (zh) * | 2006-10-25 | 2007-05-30 | 北京北方烽火科技有限公司 | 一种用于wcdma的接收信号频偏估计与补偿方法 |
CN101924723A (zh) * | 2009-06-09 | 2010-12-22 | 中兴通讯股份有限公司 | Ofdm信号解调方法和装置 |
CN102611671A (zh) * | 2012-04-13 | 2012-07-25 | 豪威科技(上海)有限公司 | 一种载波频偏估计方法、装置和系统 |
CN110691051A (zh) * | 2019-09-29 | 2020-01-14 | 天津大学 | 一种基于fft的gmsk信号频偏估计算法 |
CN111092964A (zh) * | 2019-12-19 | 2020-05-01 | 中国北方车辆研究所 | 车载网络中的设备标识方法 |
CN111654308A (zh) * | 2020-04-30 | 2020-09-11 | 中国科学院上海微系统与信息技术研究所 | 一种突发扩频弱信号的精密频偏估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8125885B2 (en) * | 2008-07-11 | 2012-02-28 | Texas Instruments Incorporated | Frequency offset estimation in orthogonal frequency division multiple access wireless networks |
-
2022
- 2022-04-08 CN CN202210380713.XA patent/CN114760178B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1972137A (zh) * | 2006-10-25 | 2007-05-30 | 北京北方烽火科技有限公司 | 一种用于wcdma的接收信号频偏估计与补偿方法 |
CN101924723A (zh) * | 2009-06-09 | 2010-12-22 | 中兴通讯股份有限公司 | Ofdm信号解调方法和装置 |
CN102611671A (zh) * | 2012-04-13 | 2012-07-25 | 豪威科技(上海)有限公司 | 一种载波频偏估计方法、装置和系统 |
CN110691051A (zh) * | 2019-09-29 | 2020-01-14 | 天津大学 | 一种基于fft的gmsk信号频偏估计算法 |
CN111092964A (zh) * | 2019-12-19 | 2020-05-01 | 中国北方车辆研究所 | 车载网络中的设备标识方法 |
CN111654308A (zh) * | 2020-04-30 | 2020-09-11 | 中国科学院上海微系统与信息技术研究所 | 一种突发扩频弱信号的精密频偏估计方法 |
Non-Patent Citations (3)
Title |
---|
Low-complexity_joint_estimation_of_frequency_offset_and_carrier_phase_for_M-ary_PSK;Dae-Ki Hong et al;《Proceedings of 2004 International Symposium on Intelligent Signal PROCESSING and Communication SYstems,ISPACS 2004》;20041119;第639-644页 * |
一种基于频域插值的非同频采样信号的补偿方法;邓轶亮;于宏毅;沈彩耀;;信息工程大学学报;20121215(第06期);全文 * |
卫星移动通信信号处理同步技术研究;陈志星;杨金孝;纪鹏;;计算机测量与控制;20180625(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114760178A (zh) | 2022-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11108612B2 (en) | Anti-interference signal detection and synchronization method for wireless broadband communication system | |
CN110071756B (zh) | 一种旋翼遮挡下的突发信号检测和载波同步方法 | |
KR100457987B1 (ko) | 심볼처리가향상된전송시스템및수신기 | |
EP2530894B1 (en) | Method and apparatus for implementing automatic frequency control | |
CN101449534B (zh) | 具有一个或多个接收器的无线通信系统中的联合分组检测 | |
CN107911329B (zh) | 一种信号分析仪ofdm信号解调方法 | |
US20030043947A1 (en) | GFSK receiver | |
CN111683034B (zh) | 一种基于ofdm的大多普勒无线通信时频同步方法 | |
US20120114031A1 (en) | OFDM signal demodulation method and device thereof | |
JP2001520484A (ja) | 無線通信のための同期方法およびシステム | |
CN101414990A (zh) | 一种单载波频域均衡系统中的载波频偏和时延的捕获方法 | |
CN108881089B (zh) | 遥测系统中soqpsk调制信号的导频检测方法 | |
CN113259291B (zh) | 利用水声连续信号动态多普勒跟踪实现的相位补偿方法 | |
EP4322489A1 (en) | M-fsk modulation-based receiver and reception method therefor | |
CN111147406B (zh) | 基于恒虚警信号检测及频偏校正的asm体制解调方法 | |
CN113253313A (zh) | 一种改进pmf-fft系统多普勒频率估计精度的方法 | |
CN111654308B (zh) | 一种突发扩频弱信号的精密频偏估计方法 | |
CN112738000B (zh) | 一种pss分块互相关检测方法 | |
CN114760178B (zh) | 一种基于多径能量的频偏估计方法 | |
TW201412062A (zh) | 數位資料輔助頻率偏移估計 | |
CN110290084B (zh) | 一种基于数据频率能量峰值的短波信道盲符号同步方法 | |
CN112968850A (zh) | 微弱bpsk信号载波捕获方法 | |
CN115499909B (zh) | 一种用于超宽带的多级相位跟踪方法及装置 | |
US7346098B2 (en) | Communication receiver | |
CN114465691A (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 |