CN100355255C - 基于统计平均的大搜索范围ofdm系统载波同步方法 - Google Patents
基于统计平均的大搜索范围ofdm系统载波同步方法 Download PDFInfo
- Publication number
- CN100355255C CN100355255C CNB2006100766919A CN200610076691A CN100355255C CN 100355255 C CN100355255 C CN 100355255C CN B2006100766919 A CNB2006100766919 A CN B2006100766919A CN 200610076691 A CN200610076691 A CN 200610076691A CN 100355255 C CN100355255 C CN 100355255C
- Authority
- CN
- China
- Prior art keywords
- carrier
- phase
- carrier frequency
- ofdm
- integer part
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000001360 synchronised effect Effects 0.000 title claims description 6
- 125000004122 cyclic group Chemical group 0.000 claims description 20
- 241001269238 Data Species 0.000 claims description 3
- 230000021615 conjugation Effects 0.000 claims description 2
- 238000012935 Averaging Methods 0.000 abstract 1
- 230000005540 biological transmission Effects 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- ODKSFYDXXFIFQN-UHFFFAOYSA-M argininate Chemical compound [O-]C(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-M 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
一种基于统计平均的大搜索范围OFDM系统载波同步方法。由于数据符号是独立同分布的随机数据符号,如果子载波数足够大的话,可以认为由于载频偏移和差分调制所引起的对第k个子载波的相位差,是在(-π,π)上均匀分布的随机相位。假设子载波序号的取值区间为(-W,-W+1,…,-1,0,1,…,W-1,W),可以通过对所有子载波进行统计平均的方法来对载频偏移进行估计。为了提高z的估计精度,还可以采用对连续多个OFDM符号数据求平均的方法。
Description
所属技术领域
本发明涉及数字通信及雷达领域。
背景技术
在基于OFDM传输技术的数字地面电视广播系统中,由于信号采用了相互正交的多载波的调制技术,因此对时间及相位同步问题更为敏感,尤其是对载波频率及相位的同步要求与单载波系统相比也更为严格。
如果收发两端载波频率的失步,除了会对直达波信号造成一定衰减和相位失真外,还会造成子载波间的相互干扰,从而破坏OFDM系统内子载波之间的正交性,如果接收端不采取严格的同步措施,会对整个系统的性能带来严重的影响。
参照如图1所示的OFDM系统,由发射机输出的一个OFDM符号块可表示为:
其中:
为第k个子载波的频率。
当接收端存在频率偏差f0时,则接收到的经过载波解调后的信号可记为:
令:φ0=2πf0t,则上式可写为:
上式可看出,载波频率的偏差与载波相位存在固定的关系,随着载波频率的偏差的改变,子载波经过解调后的相位也会发生相应变化。
假设一个OFDM信号共有128个子载波,每个子载波均采用QPSK调制方式,QPSK星座为[1+j,1-j,-1+j,-1-j]。则当载波频率偏差f0为子载波间隔Δf的0.1时(即f0=0.1Δf),OFDM信号的星座如图2所示。当载波频率偏差f0为子载波间隔Δf的0.3时(即f0=0.3Δf),OFDM信号的星座如图3所示。图2及图3中,*表示载波频率没有偏差时星座的标准位置。
由图2及图3可以看出,当频率偏差较小时(f0=0.1Δf时),f0对OFDM星座图的影响还可忍受,但当频率偏差达到f0=0.3Δf时,星座图已基本无法辨认。
目前,利用循环前缀进行同步的方法主要有两种。其中,一种是在符号定时已知,并且载波频率偏移小于子载波间隔的1/2时,采用循环前缀与OFDM尾部符号进行相关的方法进行载频偏移ε估计的方法,即:
其中L要小于或等于保护间隔的长度。估计出的频偏可以用来对接收端的时钟进行校正,经过若干符号周期后,载波频偏便会下降。如果传输信道是AWGN,则保护间隔与OFDM符号的尾部除了由于频偏引起的相位差别外,其它完全相同,此时相关效果较好,如果传输信道是频率选择性信道,则在保护间隔内可能有些采样点受到损伤,但仍能保证很好的相关性,因此算法的估计精度较高。不过,其频偏的估计范围只有1/2个子载波间隔。
另外一种利用循环前缀进行同步估计的算法是假设在AWGN信道中,SNR已知,则循环前缀与OFDM尾部的相关系数为:
其中σs 2和σn 2分别为有用信号和加性白噪声的能量。
分别计算下面的相关方程和能量方程:
则定时偏差d及载波频偏ε估计方程为:
上面对载波频偏ε的估计只有当ε小于1/2子载波间隔时才有效。并且当传输信道不是AWGN时,(8)式便不能保证是最大似然估计。
发明内容
针对上述算法的不足,本发明提出了一种OFDM传输中的新的同步方法,扩大了载波频偏的估计范围。本发明主要内容为:
对于载频偏移的整数部分z所引起的对第k个子载波的相位差θz,其估计方法可用公式表示为:
同时,对于采用更多电平调制的系统,通过调整循环前缀的采样数L,来改变星座点的数量,使星座点的数量呈2的幂级数下降。
本发明采用的OFDM系统模型如图4所示,其中每个OFDM符号有用信息部分的采样点数为N,循环前缀的采样点数为L,子载波数为N,这样一个OFDM符号的总采样点数为N+L。首先考虑传输信道为AWGN型信道。这样,信号s(k)只受到加性白噪声n(k)的影响。
在上述模型中,考虑如下两种系统偏差:第一,OFDM符号定时偏差,该偏差会造成系统SNIR下降,并且还会对数据符号相位造成影响;第二,载波频率的偏移,即接收端本地时钟与发射端时钟偏差所导致的所有子载波频率发生偏移。第一种定时偏差可以用信道冲激响应时延δ(k-d)来表示,其中d表示未知的定时偏差量。第二种载波频率的偏移可表示为相应的相位偏移2πtf0,其中f0代表载波频率偏移,通常f0可用子载波间隔的倍数来表示,即f0=(z+y)/T,其中z为整数部分,y为小数部分,且规定|y|≤1/2。
首先对载频偏移的小数部分y进行估计。
定义OFDM符号定时估计器为:
由上可知,只有当d为循环前缀的起始时刻时,γ(d)才会出现最大,在此之前γ(d)近似为0,γ(d)达到最大后,逐渐变小直至再次变为0,为了检测γ(d)突变的峰值可以通过计算γ(d)前后两个时刻的商来实现:
因为γ(d)从近似为0突变到一个比较大的,在理想情况,当d为循环前缀的起始时刻时,M(d)可变为无穷大。
由11式得出符号定时d后,则对这一时刻的载波频率偏移小数部分y可按下式进行估计:
估计出载频偏移的小数部分y后,便可对其进行补偿,这可以通过在时域内给每个OFDM符号乘以补偿因子exp(-j2πyt/T)来完成,此时,载频偏移只剩下整数部分z。
假设两个连续OFDM符号的序号为n和n+1,对它们进行exp(-j2πyt/T)补偿后的DFT的输出分别记为:yn,k和yn+1,k,下标k表示第k个子载波。这样,yn,k和yn+1,k之间相位的差别便只包含载频偏移的整数部分。如果OFDM系统采用的是差分调制方案,那么,yn,k和yn+1,k之间的相位差别还包含了由于调制所引起的相位差别。
假设由于载频偏移的整数部分z和差分调制所引起的对第k个子载波的相位差分别为θz和θk,则有:
其中arg(·)为主值相位函数,
θk对应于相应调制数据符号在星座图中的矢量相位,由于数据符号是独立同分布的随机数据符号,如果子载波数足够大的话,可以认为θk是在(-π,π)上均匀分布的随机相位。
为了估计θz,我们计算yn,k和yn+1,k之间相位差的复对数,即:
Ln(exp(j(θz+θk)))=Ln(exp(jθz))+Ln(exp(jθk)) (15)
假设子载波序号的取值区间为(-W,-W+1,…,-1,0,1,…,W-1,W),对上式的所有子载波进行算术平均,则有:
由于θk是在(-π,π)上均匀分布的随机相位,则在W较大时有:
代入16式得:
于是对θz的估计方程为:
为了提高
的估计精度,可以采用对连续多个OFDM符号数据求平均的方法。由上式估计出
后,便可以通过ejθz的星座图来确定载频偏移的整数部分z。由14式知,在星座图中,ejθz的位置与OFDM符号的循环前缀Tg和有用信息T有关。通常在OFDM系统中,为了IDFT/DFT计算方便,有用信息T内的采样数N一般选取为2的幂级数,并且为保证接收机的采样率恒定,循环前缀长度Tg内采样数L也应为整数。图5给出了载波为32DPSK调制时,当N=32和L=6及N=32和L=5时,ejθz的星座图。
图中可以看出,当L取奇数(N=32和L=5)时,z在[-16,16]范围内都有唯一的星座点和其对应;当L取偶数(N=32和L=6)时,z在[-8,8]范围内也有唯一的星座点和其对应,但总的星座点数为L取奇数时的一半。假设OFDM信号的总带宽为B,那么此时本算法对载频偏移的搜索范围至少为
只有超出这个范围后,z的取值在星座图上才会出现模糊性。因此,算法极大地扩展了同步搜索的动态范围,也满足了在更恶劣的环境条件下,系统可靠工作的性能要求。
对于采用更多电平调制的系统,由于星座图中的点更加密集,就会使得查找z的精度下降。此时,可以通过调整循环前缀的采样数L,来改变星座点的数量,比如L值取为2的幂级数,可使星座点的数量呈2的幂级数下降。
附图说明
图1——OFDM系统的基本模型框图。
图2——OFDM星座图(f0=0.1Δf时)。
图3——OFDM星座图(f0=0.3Δf时)。
图4——本方法所采用的OFDM系统模型框图。
图5——载频偏移星座图(每个星座点旁边的数字为载频偏移的整数部分z)。
图6——频偏方差与循环前缀长度关系曲线图。
图7——频偏方差与SNR关系曲线图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
改进的基于循环前缀的载波频率盲同步方法的噪声性能可以通过分析接收信号r(k)来实现。根据前面定义的系统模型,在OFDM带宽范围内,接收信号r(k)实际上是发射信号s(k)与噪声信号n(k)的叠加,即:
r(k)=s(k)+n(k) (20)
注意上述信号实际上均可视为一组随机复数序列,假设s(k)与n(k)的实部和虚部的方差分别为:
由此可得系统的输入信噪比为:
由于在对载波频率偏移的估计时,频偏被分为整数部分和小数部分分别来估计,因此载波频偏的方差σf0 2也由载波频偏整数部分方差σz 2和载波频偏小数部分方差σy 2两部分组成。假设由于载波频偏小数部分引起的在一个OFDM符号周期内相位偏移记为exp(jθy),由于载波频偏整数部分引起的在一个OFDM符号周期内相位偏移记为exp(jθz)。
由12式可知,载波频偏小数部分的估算是通过计算一个OFDM符号周期内循环前缀与该符号尾部的相位差别来实现的,该相位差是通过计算循环前缀各采样值的共轭与相应符号尾部采样值的乘积来得到的,即:
式中第一项为没有噪声理想情况下的期望值。上式的近似是考虑到一般输入SNR大于1,因此,nk *nk+N的功率要远小于其它三项的功率,故可忽略不记。
上式的结果也可以表示成复矢量X+jY的形式,其中:
由21式和22式可得:
因此可得:
由于循环前缀的缘故,一般θ*较小,且σn与σs相比也很小,因此可近似的表示为右边的式子。
因此可得θy的方差为:
由于相位偏差与频率偏差的关系为:
式中L为循环前缀采样点数。因此可得载波频偏小数部分的方差σy 2为:
对于载波频偏整数部分z的估算是通过计算两个连续OFDM符号在时域上补偿后DFT输出的相位差来实现的,即首先需要估计出相位偏差,再对比星座图来查找z。同样按照上面的方法,可得载波频偏整数部分引起的相位偏差的方差为:
其中,N为IDFT/DFT长度。
由14式知:
因此,可得载波频偏整数部分的方差为:
由33式和36式可得总的载波频偏的方差σf0 2为:
如果SNR值对于接收到的载波频偏的方差太大,那么可以对12式和19式取多个符号进行平均,以增加有效的信噪比。对于K个符号的平均,σf0 2变为:
为了检验上述算法的性能,本发明做了大量的模拟仿真实验。在实验中,假定OFDM系统的子载波数N=512,调制采用DQPSK,每一次仿真都使用500个符号,图6给出了在SNR=8dB,10dB,16dB时,以循环前缀保护间隔长度L为函数的载波频偏的方差。图中的虚线为实际情况下的曲线,实线为采用5.3.34式近似后的曲线,这两条曲线的差别显示了采用30式简化后产生的影响。从图中可以看出,载波频偏的方差随着循环前缀的长度增加,而逐渐改善,但当循环前缀长度超过一定长度(在图中约为15)门限时,载波频偏的方差趋于稳定。另外,随着SNR的增加,两条曲线的差别也越来越小,对于本次仿真在大于16dB以后,其中的差别已经可以忽略不计。
图7给出了在L=9,18,52时,以信噪比SNR为函数的载波频偏的方差。图中的虚线仍为实际情况下的曲线,实线为采用30式近似后的曲线,这两条曲线的差别显示了采用30式简化后产生的影响。从图中可以看出,载波频偏的方差随着循环前缀的长度增加,而逐渐改善,并且随着SNR的增加,两条曲线的差别也越来越小,在SNR大于16dB以后,其中的差别已经可以忽略不计,这点与图6是一致的。
Claims (3)
1.一种基于统计平均的大搜索范围OFDM系统载波同步方法,载波频率的偏差与载波相位存在固定的关系,即载波频率的偏移可表示为相应的载波相位偏移,载波相位偏移又可以分解为整数部分偏移及小数部分偏移。当子载波数足够大时,由差分调制所引起的对第k个子载波的相位差θk可以认为是在(-π,π)上均匀分布的随机相位,在采用对连续多个OFDM符号数据求统计平均后,并通过星座图来确定载频相位偏移的整数部分z的方法,得到载波同步,从而扩展了同步搜索的动态范围,其特征在于得出符号定时d后,首先采用循环前缀与OFDM尾部符号进行相关的方法估计出载频偏移的小数部分y,然后在时域内给每个OFDM符号乘以补偿因子exp(-j2πyt/T),此时载频相位偏移只剩下整数部分,当OFDM系统采用差分调制时,相位差别还包含了由于调制所引起的相位差,计算连续多个OFDM符号之间相位差的复对数,进行统计平均,并结合相对应的星座图便可确定载频相位偏移的整数部分,获取载波同步的步骤如下:
(1)首先对载频偏移的小数部分y进行估算,估算是通过计算一个OFDM符号周期内循环前缀与该符号尾部的相位差别来实现的,该相位差是通过计算循环前缀各采样值的共轭与相应符号尾部采样值的乘积来得到;
(2)估计出载频偏移的小数部分y后,通过在时域内给每个OFDM符号乘以补偿因子exp(-j2πyt/T)来完成对载波偏移小数部分的补偿;
2.根据权利要求1,在估计出载频偏移的小数部分y,并在时域内对序号为n和n+1的OFDM符号乘以补偿因子exp(-j2πyt/T)后,估算出载频相位偏移的整数部分,其特征是通过对
所有子载波进行算术平均,从而得出载频相位偏移的整数部分的估计方程为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100766919A CN100355255C (zh) | 2006-04-29 | 2006-04-29 | 基于统计平均的大搜索范围ofdm系统载波同步方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100766919A CN100355255C (zh) | 2006-04-29 | 2006-04-29 | 基于统计平均的大搜索范围ofdm系统载波同步方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1835494A CN1835494A (zh) | 2006-09-20 |
CN100355255C true CN100355255C (zh) | 2007-12-12 |
Family
ID=37003097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100766919A Expired - Fee Related CN100355255C (zh) | 2006-04-29 | 2006-04-29 | 基于统计平均的大搜索范围ofdm系统载波同步方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100355255C (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101534272B (zh) * | 2008-03-14 | 2011-11-30 | 中兴通讯股份有限公司 | 一种正交频分复用系统中载波的重排方法和装置 |
CN102075460B (zh) * | 2009-11-20 | 2014-01-01 | 中兴通讯股份有限公司 | 一种基于数据的频偏估计方法和装置 |
CN102104439B (zh) * | 2009-12-21 | 2016-08-31 | 上海贝尔股份有限公司 | 一种传输方法及其设备 |
CN104683280B (zh) * | 2014-12-02 | 2018-04-24 | 北京星河亮点技术股份有限公司 | DFT-s-OFDM系统大频偏的精确估计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20020056986A (ko) * | 2000-12-30 | 2002-07-11 | 박태진 | 직교 주파수 분할 다중 통신 시스템의 프레임 구조 및분산 파일롯 부채널을 이용한 변조기 및 복조기 |
CN1516946A (zh) * | 2001-06-20 | 2004-07-28 | 皇家菲利浦电子有限公司 | 各正交频分复用码元中采用对称相位调整的频率校正 |
CN1705302A (zh) * | 2004-05-27 | 2005-12-07 | 华为技术有限公司 | 一种正交频分复用时频同步的方法 |
-
2006
- 2006-04-29 CN CNB2006100766919A patent/CN100355255C/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20020056986A (ko) * | 2000-12-30 | 2002-07-11 | 박태진 | 직교 주파수 분할 다중 통신 시스템의 프레임 구조 및분산 파일롯 부채널을 이용한 변조기 및 복조기 |
CN1516946A (zh) * | 2001-06-20 | 2004-07-28 | 皇家菲利浦电子有限公司 | 各正交频分复用码元中采用对称相位调整的频率校正 |
CN1705302A (zh) * | 2004-05-27 | 2005-12-07 | 华为技术有限公司 | 一种正交频分复用时频同步的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1835494A (zh) | 2006-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101425999B (zh) | 正交频分复用接收机的载频偏差同步的方法及装置 | |
CN109802912B (zh) | 宽带无线通信系统的同步方法、装置、设备及存储介质 | |
CN101588338B (zh) | 一种适用于分组传输的ofdm载波频偏估计方法 | |
CN108040028A (zh) | 一种基于本地序列互相关检测的ofdm系统抗干扰信号检测与同步方法 | |
CN102404268B (zh) | 高速移动环境下莱斯信道中多普勒频偏估计与补偿方法 | |
US7715484B2 (en) | Orthogonal frequency division multiplexing with PN-sequence | |
KR100729726B1 (ko) | 직교 주파수 분할 다중화 방식의 통신 시스템의 타이밍획득 및 반송파 주파수 오차 추정 장치 및 방법 | |
CN104125190A (zh) | 适于低信噪比信道环境的ofdm系统符号定时同步实现方法 | |
CN104836769B (zh) | 一种基于共轭结构前导的联合定时与频率同步方法 | |
Wang et al. | Optimized joint fine timing synchronization and channel estimation for MIMO systems | |
CN100389582C (zh) | 宽带无线接入系统中正交频分复用的同步方法 | |
CN104125188B (zh) | 一种基于Zadoff‑Chu序列的OFDM频率同步方法 | |
CN101207596B (zh) | 一种同步方法及接收端设备 | |
CN102215184B (zh) | 一种上行定时偏差的估计方法及系统 | |
CN100355255C (zh) | 基于统计平均的大搜索范围ofdm系统载波同步方法 | |
CN102594740B (zh) | 一种频偏估计方法及装置 | |
CN101958872B (zh) | 搜索最佳载波频率偏移校正值的方法 | |
Diaz et al. | A sign-bit auto-correlation architecture for fractional frequency offset estimation in OFDM | |
CN106100692A (zh) | Mimo‑ofdm水声通信系统多普勒扩展估计方法 | |
KR20100054987A (ko) | Ofdm 시스템에서의 주파수 추정을 위한 장치 및 방법 | |
CN103095627B (zh) | 一种正交频分复用技术系统同步方法和电子设备 | |
TWI396415B (zh) | 正交分頻多工系統的通道長度估測方法及其估測器 | |
CN105187351B (zh) | 一种多径信道下的ofdm定时同步检测方法 | |
Yang et al. | A Robust Time-Frequency Synchronization Method for Underwater Acoustic OFDM Communication Systems | |
CN111884979B (zh) | 一种基于ofdm智能电网抗脉冲噪声的符号同步方法 |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20071212 Termination date: 20100429 |