CN103095638A - 一种多径衰落信道下ofdm系统的采样频率偏移盲估算方法 - Google Patents
一种多径衰落信道下ofdm系统的采样频率偏移盲估算方法 Download PDFInfo
- Publication number
- CN103095638A CN103095638A CN2012105554480A CN201210555448A CN103095638A CN 103095638 A CN103095638 A CN 103095638A CN 2012105554480 A CN2012105554480 A CN 2012105554480A CN 201210555448 A CN201210555448 A CN 201210555448A CN 103095638 A CN103095638 A CN 103095638A
- Authority
- CN
- China
- Prior art keywords
- alpha
- sigma
- prime
- epsiv
- infin
- 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.)
- Granted
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Noise Elimination (AREA)
Abstract
一种多径衰落信道下OFDM系统的采样频率偏移盲估算方法。根据具有相关导频的OFDM信号的循环平稳特性,利用改进后的谱函数 克服采样频率偏移带来的衰减和多径衰落信道影响的数据估计出相关值点,再根据其估计值点处的循环谱值相位偏移量来估计采样频率偏移,可以克服频偏、高斯白噪声和多径衰落信道的影响,并利用跳变变换,充分利用数据,大大提高OFDM系统的采样频率偏移的估计性能。
Description
【技术领域】
本发明属于通信技术领域,具体涉及一种多径衰落信道下OFDM系统的采样频率偏移(SFO)盲估计方法,可用于非合作通信系统中时频同步之前的OFDM系统的采样偏移估计。
【背景技术】
以正交频分复用(OFDM)为代表的多载波技术已成功应用于数字视频广播、数字音频广播、数字用户线、无线接入网、电力线通信、卫星通信等领域,其具备很好的抗多径、抗窄带干扰性能和高效的频谱利用率,是宽带无线通信研究的热点。然而,由于收发两端的晶振的不完全匹配,及移动通信系统中的多普勒频移的影响,在通信系统的收发两端不可避免地存在采样频率偏差。尤其在非合作通信中,作为一种非授权接入通信模式,采样频率对于接收端是未知的,经前期的过采样率估计之后,必定还存在一定的采样频率偏差。在采样频率偏差存在的情况下,OFDM信号经FFT变换后将产生载波间干扰和时变的相位变化,且现有的时频同步方法基本上均是在无采样频率偏差条件下进行的。因此,研究非合作通信中多径衰落信道下时频同步之前的OFDM系统的采样频率偏移估计方法具有实际的工程意义。
近年来,已有学者对非合作通信中OFDM系统的采样频率偏移估计方法进行了一定的研究,但这些方法需要发射端配合或在频偏估计之后,且研究集中于利用冗余CP与对应数据之间的相关性,算法受信道和噪声影响较大。参见LiuGang,Li Bingbing,Pang Hengli.Blind Sampling Clock Offset Estimation Algorithmfor OFDM System[C]//Proc.of CISP’09,Tianjin,China.IEEE Press,2009:1-4。Liu Gang通过对各子载波加权,然后对接收信号相关函数的虚部进行变换,得到了一种非数据辅助的采样频率同步算法。但该方法需要发射端配合,且要求信道冲击响应在各子载波位置上能量相等。参见B.Ai,Y. Shen,Z.D.Zhong,B.H.Zhang.Enhanced Sampling Clock Offset Correction Based on Time DomainEstimation Scheme[J].IEEE Transactions on Consumer Electronics,2011,57(2):696-704。B.Ai提出一种基于保护间隔的时域SCO估计方法,并相应的加入一自适应模块,提高了系统的实时性,但此方法用到的CP长度要大于最大多经时延,且需要先估计频偏。参见Castillo-Sanchez E,Lopez-Martinez F.J,Martos-NayaE,etal.Joint Time,Frequency and Sampling Clock Synchronization for OFDM-basedsystems[C].Wireless Communications and Networking Conference2009:1-6。Castillo-Sanchez E在采用CP和OFDM符号尾部数据的相关进行符号定时及载波频率联合同步方法基础上,利用定时估计值的偏移来进行采样频率偏移估计,但该算法受信道和噪声影响较大。参见Arash Zahedi-Ghasabeh,Alireza Tarighat,andBabak Daneshrad.Spectrum Sensing of OFDM Waveforms Using Embedded Pilots inthe Presence of Impairments[J].IEEE Transactions onVehicular Technology:1208-1221。Arash Zahedi-Ghasabeh提出一种基于循环平稳的谱分析方法估计采样频率偏移,但是需要知道发射端导频信息,直接提取会受到采样频率偏移带来的衰减和多径衰落信道影响,且SFO估计范围很小,频偏对其性能也有一定影响。综上说明,这些研究均存在一定缺陷,不适合于实际的非合作通信系统中频率选择性多径衰落信道,况且在低信噪比情况下,时频同步之前估计性能不理想。因此,以上的方法不适合在实际的非合作通信系统中进行应用。
【发明内容】
本发明的目的是克服上述已有技术的不足,提供了一种多径衰落信道下时频同步之前进行估计的新方法,以提高实际的非合作通信系统中采样频率偏移(SFO)估计的精度。本发明选取子载波个数为N=64个,1/4循环前缀(CP)长度,符号周期Ts=10μs,具有一对相关导频的OFDM信号作为信号源。
实现本发明目的的技术方案,包括如下步骤:
(1)对接收到的信号y(t)采样得到y[n];
(5)根据S_n前后差值为正号的个数判断SFO符号为flag(flag值取1或-1),S_n乘以flag;
(6)从1到M,将an作为相位值变化因子,设定an初始值为-π,将S_n(m)转化为(an,an+2π]范围内相位S_a(m),当S_a(m)≥an+π时,an=an+π/2,计算S_a(m+1),以此进行类推;
(7)对所得S_a进行最小二乘直线拟合,求其斜率为k_a,再根据下式估算采样频率偏移值ε:
本发明与现有技术相比具有如下优点:
1)本发明根据具有相关导频的OFDM信号的循环平稳特性,利用改进后的谱函数 这一可以克服采样频率偏移带来的衰减和多径衰落信道影响的数据估计出相关值点,再根据其估计值点处的循环谱值相位偏移量来估计采样频率偏移,可以克服频偏、高斯白噪声和多径衰落信道的影响,并利用跳变变换,充分利用数据,大大提高OFDM系统的采样频率偏移的估计性能;
2)本发明可以适用于不同标准协议的频率选择性多径衰落信道模型,多径信道径数不受CP长度限制,并且在低信噪比情况下估计性能较好。
仿真结果表明,当存在频偏和定时误差时,在三种不同类型的多径衰落信道下,该发明方法具有较好性能和良好的稳健性;在相同的仿真实验环境、相同的信号参数设置及含有频偏和定时误差的情况下,本发明具有比现有的方法具有更好的性能。说明在多径衰落信道下,本发明更适合于非合作通信系统。
【附图说明】
图1为本发明一种多径衰落信道下OFDM系统的采样频率偏移估算方法步骤图;
图2中是本发明当存在频偏和定时误差时,在三种不同类型的多径衰落信道下对OFDM系统的采样频率偏移进行估计的性能图;
图3中是在相同的仿真实验环境、相同的信号参数设置及含有频偏和定时误差的情况下,本发明与现有方法的估计性能对比图。
【具体实施方式】
以下结合具体实施例,对本发明进行详细说明。
请参考图1,本发明的具体实现步骤如下:
步骤1,对接收到的信号y(t)采样得到y[n];
设基带接受信号y(t),存在采样频率偏移ε,频偏fo,加性高斯白噪声n(t)和多径衰落信道(径数为P)影响,表示为:
Y(f)=X(f-fo)·H(f-fo)+N(f)
其中Y(f)、X(f)、H(f)和N(f),分别为接收信号、发射信号、信道响应和加性高斯白噪声的傅里叶变换,
对y(t)以频率fs进行采样,并以L长度窗口作傅里叶变换,可得:
其中 为采样周期, 为相对频率, 为相对频偏,
直接将上式代入循环谱计算表达式,可得:
根据中心极限定理,可将近似成均值为0,方差为的高斯白噪声,其中E0=(|H(f0+α0/2)|2+|H*(f0-α0/2)|2)Ex2,值与其信噪比值均会发生衰减,尤其在某些点上,衰减将非常厉害,所以直接计算谱相关函数将无法提取相关点,同时,由于多径信道下窗口个数非无限长,又受到噪声和H(m,f-εf)的影响,其他非相关点值可能会高于相关点值,导致提取出的最大值点不是我们所要的相关值点,因此,我们对循环谱计算式进行了改进,得到如下表达式:
k为差值,一般为减小其他位置的相关性,以便于提取,k要取大于1的数,为尽量利用所得数据,我们选k=2。则:
当α=α0,f=f0时,
其中, 根据中心极限定理,可近似认为是均值为0,方差为的高斯白噪声,其中,EXH0=H(f0+α0/2)H*(f0-α0/2)E[X(f0+α0/2)X*f0-α0/2)], 即可近似认为是均值为0,方差为 的高斯白噪声,
当α≠α0或f≠f0时,
可近似认为是均值为 方差为 的高斯白噪声,其中 E′+=|H(f+α/2)|2Ex, E′-=|H(f-α2)|2Ex;
此方法可以消除频率偏移带来的衰减和多径影响,但受噪声影响较大,因此将其作为粗估计方法,将其搜索出的最大值位置作为下一步细估计频率范围的中心点。
其中,α∈[α1-range,α1+range],f∈[f1-range/2,f1+range/2],range为估计范围大小;
EXH=|H(f+α/2)H*(f-α/2)|E[|X(f+α/2)X*(f-α/2)|]
步骤5,根据S_n前后差值为正号的个数判断SFO符号为flag(flag值取1或-1),S_n乘以flag;
步骤6,从1到M,以an确定相位值范围,设定其初始值为-π,将S_n(m)转化为(an,an+2π]范围内相位S_a(m),当S_a(m)≥an+π时,an=an+π2,计算S_a(m+1),以此类推;
步骤7,对所得S_a进行最小二乘直线拟合,求其斜率为k_a,再根据下式估计采样频率偏移值ε:
仿真内容与结果:
为了验证本文方法的有效性,通过MATLAB仿真软件进行仿真实验,其所使用的仿真条件为:载波个数为N=64个,1/4循环前缀(CP)长度,符号周期Ts=10μs,具有一对相关导频的OFDM信号作为信号源,信道为SU13径信道、TU6径信和指数衰落9径信道三种多径信道,采样频率为8MHz,采样频率偏移为1×10-3,相对频率偏移为4.25,定时误差为20个符号,蒙特卡洛仿真次数为200次。
仿真当存在频偏和定时误差时,在三种不同类型的多径衰落信道下对OFDM系统的采样频率偏移进行估计,得到估计值的MSE曲线,其中从图2中可以看出,多径衰落信道对本发明方法有一定影响,但其影响不大,说明本方法具有良好的稳健性。
仿真在在相同的仿真实验环境、相同的信号参数设置及含有频偏和定时误差情况下,本发明方法与传统数据相关方法进行性能对比,其结果如图3所示。从图3中可以看出,本发明方法的性能有显著提高。由此可见,本发明方法明显优于现有的采样频率偏移方法。
Claims (3)
1.一种多径衰落信道下OFDM系统的采样频率偏移估算方法,其特征在于:包括如下步骤:
(1)对接收到的信号y(t)采样得到y[n];
(5)根据S_n前后差值为正号的个数判断SFO符号为flag(flag值取1或-1),S_n乘以flag;
(6)从1到M,将an作为相位值变化因子,设定an初始值为-π,将S_n(m)转化为(an,an+2π]范围内相位S_a(m),当S_a(m)≥an+π时,an=an+π/2,计算S_a(m+1),以此进行类推;
(7)对所得S_a进行最小二乘直线拟合,求其斜率为k_a,再根据下式估算采样频率偏移值ε:
设基带接受信号y(t),存在采样频率偏移ε,频偏fo,加性高斯白噪声n(t)和多径衰落信道(径数为P)影响,表示为:
Y(f)=X(f-fo)·H(f-fo)+N(f)
其中Y(f)、X(f)、H(f)和N(f),分别为接收信号、发射信号、信道响应和加性高斯白噪声的傅里叶变换,
对y(t)以频率fs进行采样,并以L长度窗口作傅里叶变换,可得:
其中 为采样周期, 为相对频率, 为相对频偏,
直接将上式代入循环谱计算表达式,可得:
根据中心极限定理,可将近似成均值为0,方差为的高斯白噪声,其中E0=(|H(f0+α0/2)|2+|H*(f0-α0/2)|2)Ex/2,由于循环谱不能直接提取相关点值,改进谱函数幅值计算式,改进后的谱函数幅值的计算式可表示如下:
k为差值,选k=2,则:
当α=α0,f=f0时,
其中, 根据中心极限定理,可近似认为是均值为0,方差为的高斯白噪声,其中,EXH0=H(f0+α0/2)H*(f0-α0/2)E[X(f0+α0/2)X*(f0-α0/2)], 即可近似认为是均值为0,方差为 的高斯白噪声,
可近似认为是均值为 方差为 的高斯白噪声,其中E′0+=|H(f0+α0/2)|2Ex, E′0-=|H(f0-α0/2)|2Ex,
当α≠α0或f≠f0时,
可近似认为是均值为 方差为 的高斯白噪声,其中 E′+=|H(f+α/2)|2Ex, E′-=|H(f-α2)|2Ex;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210555448.0A CN103095638B (zh) | 2012-12-19 | 2012-12-19 | 一种多径衰落信道下ofdm系统的采样频率偏移盲估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210555448.0A CN103095638B (zh) | 2012-12-19 | 2012-12-19 | 一种多径衰落信道下ofdm系统的采样频率偏移盲估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103095638A true CN103095638A (zh) | 2013-05-08 |
CN103095638B CN103095638B (zh) | 2016-06-08 |
Family
ID=48207783
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210555448.0A Expired - Fee Related CN103095638B (zh) | 2012-12-19 | 2012-12-19 | 一种多径衰落信道下ofdm系统的采样频率偏移盲估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103095638B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105516043A (zh) * | 2015-11-26 | 2016-04-20 | 中国电子科技集团公司第三十研究所 | 一种基于ofdm的多载波通信系统频偏估计方法及装置 |
CN106534033A (zh) * | 2016-12-06 | 2017-03-22 | 西安电子科技大学 | 一种多径信道下ofdm/oqam时频联合同步方法 |
CN108605031A (zh) * | 2016-02-15 | 2018-09-28 | 瑞典爱立信有限公司 | 以最小采样率操作的nb-iot接收器 |
CN109257128A (zh) * | 2018-11-01 | 2019-01-22 | 南京邮电大学 | 一种基于傅里叶级数拟合去噪的频谱信号识别方法及系统 |
CN110031675A (zh) * | 2019-04-19 | 2019-07-19 | 南京大学 | 一种数采系统实际采样频率的测量方法 |
CN110850385A (zh) * | 2019-11-20 | 2020-02-28 | 桂林电子科技大学 | 基于无源雷达与循环谱的无人机微动特征检测方法 |
CN113207167A (zh) * | 2021-05-10 | 2021-08-03 | 重庆邮电大学 | 基于序列最小二乘的一致性时钟同步频率偏移估计方法 |
US11239961B2 (en) | 2016-02-15 | 2022-02-01 | Telefonaktiebolaget Lm Ericsson (Publ) | Receiver circuit and methods |
US11245558B2 (en) | 2016-02-15 | 2022-02-08 | Telefonaktiebolaget Lm Ericsson (Publ) | Method and apparatus for generating OFDM signals |
CN114384968A (zh) * | 2021-12-29 | 2022-04-22 | 西安电子科技大学 | 利用特定频偏点可调控大小的相位噪声生成方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1435037A (zh) * | 1999-12-22 | 2003-08-06 | 汤姆森特许公司 | 在正交频分多路复用系统中采样频率偏移的纠正 |
CN1881823A (zh) * | 2005-06-17 | 2006-12-20 | 美国博通公司 | 校正通信系统中数据流的数据包的采样频率偏移的方法 |
CN1901527A (zh) * | 2005-07-19 | 2007-01-24 | 三星电子株式会社 | 用于正交频分复用系统的采样频率偏移估计设备和方法 |
US20090168641A1 (en) * | 2005-10-14 | 2009-07-02 | Via Technologies, Inc. | Method and circuit for frequency offset estimation in frequency domain in the orthogonal frequency division multiplexing baseband receiver for ieee 802.11a/g wireless lan standard |
US20100158144A1 (en) * | 2008-12-22 | 2010-06-24 | Samsung Electro-Mechanics Co., Ltd. | Sampling frequency offset estimation apparatus and method of ofdm system |
-
2012
- 2012-12-19 CN CN201210555448.0A patent/CN103095638B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1435037A (zh) * | 1999-12-22 | 2003-08-06 | 汤姆森特许公司 | 在正交频分多路复用系统中采样频率偏移的纠正 |
CN1881823A (zh) * | 2005-06-17 | 2006-12-20 | 美国博通公司 | 校正通信系统中数据流的数据包的采样频率偏移的方法 |
CN1901527A (zh) * | 2005-07-19 | 2007-01-24 | 三星电子株式会社 | 用于正交频分复用系统的采样频率偏移估计设备和方法 |
US20090168641A1 (en) * | 2005-10-14 | 2009-07-02 | Via Technologies, Inc. | Method and circuit for frequency offset estimation in frequency domain in the orthogonal frequency division multiplexing baseband receiver for ieee 802.11a/g wireless lan standard |
US20100158144A1 (en) * | 2008-12-22 | 2010-06-24 | Samsung Electro-Mechanics Co., Ltd. | Sampling frequency offset estimation apparatus and method of ofdm system |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105516043A (zh) * | 2015-11-26 | 2016-04-20 | 中国电子科技集团公司第三十研究所 | 一种基于ofdm的多载波通信系统频偏估计方法及装置 |
CN105516043B (zh) * | 2015-11-26 | 2018-11-16 | 中国电子科技集团公司第三十研究所 | 一种基于ofdm的多载波通信系统频偏估计方法及装置 |
US11245558B2 (en) | 2016-02-15 | 2022-02-08 | Telefonaktiebolaget Lm Ericsson (Publ) | Method and apparatus for generating OFDM signals |
CN108605031A (zh) * | 2016-02-15 | 2018-09-28 | 瑞典爱立信有限公司 | 以最小采样率操作的nb-iot接收器 |
US11239961B2 (en) | 2016-02-15 | 2022-02-01 | Telefonaktiebolaget Lm Ericsson (Publ) | Receiver circuit and methods |
CN106534033B (zh) * | 2016-12-06 | 2019-11-05 | 西安电子科技大学 | 一种多径信道下ofdm/oqam时频联合同步方法 |
CN106534033A (zh) * | 2016-12-06 | 2017-03-22 | 西安电子科技大学 | 一种多径信道下ofdm/oqam时频联合同步方法 |
CN109257128B (zh) * | 2018-11-01 | 2021-05-11 | 南京邮电大学 | 一种基于傅里叶级数拟合去噪的频谱信号识别方法及系统 |
CN109257128A (zh) * | 2018-11-01 | 2019-01-22 | 南京邮电大学 | 一种基于傅里叶级数拟合去噪的频谱信号识别方法及系统 |
CN110031675A (zh) * | 2019-04-19 | 2019-07-19 | 南京大学 | 一种数采系统实际采样频率的测量方法 |
CN110850385A (zh) * | 2019-11-20 | 2020-02-28 | 桂林电子科技大学 | 基于无源雷达与循环谱的无人机微动特征检测方法 |
CN113207167A (zh) * | 2021-05-10 | 2021-08-03 | 重庆邮电大学 | 基于序列最小二乘的一致性时钟同步频率偏移估计方法 |
CN113207167B (zh) * | 2021-05-10 | 2022-03-08 | 重庆邮电大学 | 基于序列最小二乘的一致性时钟同步频率偏移估计方法 |
CN114384968A (zh) * | 2021-12-29 | 2022-04-22 | 西安电子科技大学 | 利用特定频偏点可调控大小的相位噪声生成方法及系统 |
CN114384968B (zh) * | 2021-12-29 | 2023-09-29 | 西安电子科技大学 | 利用特定频偏点可调控大小的相位噪声生成方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103095638B (zh) | 2016-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103095638A (zh) | 一种多径衰落信道下ofdm系统的采样频率偏移盲估算方法 | |
US7961816B2 (en) | Device for and method of signal synchronization in a communication system | |
CN110290581B (zh) | 一种5g系统中的快速时频同步方法及终端 | |
CN101447970B (zh) | 利用训练序列进行lofdm系统定时和载波同步的方法 | |
CN102185822B (zh) | 一种ofdm/oqam系统及其时频同步方法 | |
CN101883412B (zh) | Lte在高速移动条件下初始小区搜索方法 | |
CN102404268B (zh) | 高速移动环境下莱斯信道中多普勒频偏估计与补偿方法 | |
US9065630B1 (en) | Systems and methods for detecting secondary synchronization signals in a wireless communication system | |
CN104410590A (zh) | 一种基于压缩感知的短波ofdm抑制干扰联合信道估计方法 | |
CN101651648A (zh) | 用于在无线通信系统中执行初始同步的系统和方法 | |
US20140036779A1 (en) | Robust downlink timing synchronization method in lte system | |
CN100518159C (zh) | 时域同步正交频分复用接收机及其信号处理方法 | |
CN101414990A (zh) | 一种单载波频域均衡系统中的载波频偏和时延的捕获方法 | |
WO2014063275A1 (zh) | 一种远距离同频干扰源的确定方法以及定位方法 | |
CN108989259B (zh) | 无线综测仪窄带物理上行共享信道的时偏估计方法及系统 | |
CN113438730B (zh) | 一种基于gfdm信号的无线定位方法 | |
CN105187352A (zh) | 一种基于ofdm前导的整数频偏估计方法 | |
CN111107028B (zh) | 一种5g系统基于pss与sss联合的频偏估计方法 | |
CN107070833A (zh) | 面向具有对称性的同步信号的检测方法、同步方法及终端 | |
CN102215184B (zh) | 一种上行定时偏差的估计方法及系统 | |
CN107707499B (zh) | 一种Alpha稳定分布噪声下OFDM信号调制参数估计方法 | |
CN101232315B (zh) | 一种WiMAX系统中多天线测距码的检测方法 | |
CN103078819B (zh) | 精符号同步方法及其装置 | |
US9369329B1 (en) | Low-complexity non-data-aided estimation of symbol time offset in OFDM systems | |
CN102487364B (zh) | 一种信道估计方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160608 Termination date: 20211219 |
|
CF01 | Termination of patent right due to non-payment of annual fee |