CN103281282A - 一种dvb-t2***中频偏估计的方法 - Google Patents

一种dvb-t2***中频偏估计的方法 Download PDF

Info

Publication number
CN103281282A
CN103281282A CN2013102329004A CN201310232900A CN103281282A CN 103281282 A CN103281282 A CN 103281282A CN 2013102329004 A CN2013102329004 A CN 2013102329004A CN 201310232900 A CN201310232900 A CN 201310232900A CN 103281282 A CN103281282 A CN 103281282A
Authority
CN
China
Prior art keywords
symbol
frequency
sequence
frequency deviation
delta
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
Application number
CN2013102329004A
Other languages
English (en)
Other versions
CN103281282B (zh
Inventor
朱志敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan Goke Microelectronics Co Ltd
Original Assignee
Hunan Goke Microelectronics Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan Goke Microelectronics Co Ltd filed Critical Hunan Goke Microelectronics Co Ltd
Priority to CN201310232900.4A priority Critical patent/CN103281282B/zh
Publication of CN103281282A publication Critical patent/CN103281282A/zh
Application granted granted Critical
Publication of CN103281282B publication Critical patent/CN103281282B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

本发明公开了一种DVB-T2***中频偏估计的方法,使用接收机中重新生成的P1符号,和接收信号做相关,减小残留频偏的范围;使用P2符号导频的时域信息和接收信号做相关,估计整数频偏;使用P2符号导频的时域信息和接收信号做相关,估计小数频偏;多次计算小数频偏,求其均值提高准确度。本发明可以通过计算多个小数频偏然后求均值的方法提高估计的准确性;本发明避免了整数频偏16K,32K模式下复杂的相关运算,提高了整数频偏的检测速度和估计精度。

Description

一种DVB-T2***中频偏估计的方法
技术领域
本发明涉及一种适用于DVB-T2***、基于OFDM***的多载波***,例如DVB-T,DTMB等***的整数频偏和小数频偏的估计方法。
背景技术
DVB-T2***是一个以OFDM多载波技术为基础的地面数字电视广播***,其基带信号被分割成一个一个的超帧,每个超帧里面包含最多255个T2帧,每个T2帧以一个P1符号位开始,后面跟着N个P2符号,之后是多个数据符号,其中P2符号和数据符号都有相同的FFT长度和保护间隔(GI)长度。
P1符号的结构固定,长度为2048,图2中可以看到,A部分是由信令S1和S2经过编码之后,做1K IFFT得到的,C和B分别是A的前半部分和后半部分的加权复制,其生成流程见图3。
P2符号的特点是导频比较密集,导频间隔是3的倍数或者6的倍数,由FFT的长度和MISO的类型决定。
在OFDM***中,子载波之间是正交的,***对频偏很敏感。因此,如何准确快速的估计出当前***的频偏并且纠正之,是评价一个接收机性能好坏的重要指标。
在多载波***中,频偏估计一般分为两个部分,即整数频偏估计和小数频偏估计。下面就现有的技术分别对两个部分进行说明。
1)小数频偏的估计方法
在DVB-T***中,小数频偏是在进行符号同步以及模式检测的过程中计算出来的,根据ML算法,利用OFDM符号的循环前缀,进行滑动相关,得到相关峰值,峰值的位置就是OFDM符号的起始位置,而峰值的相位信息就包含了小数频偏的信息,原理描述如下:
定义相关函数:
R i ( n ) = 1 N i Σ j = 0 N i Q - 1 r ( n - j ) r * ( n - j - N i )
其中,Ni是FFT的长度,即子载波的个数,Q是GI的模式。
在找到相关峰值之后,对其相位求角度,即可得到小数频偏的估计值
Δ f ^ = - 1 2 π arg ( Λ ) = - 1 2 π tan - 1 [ Im ( Σ n = 0 N g - 1 r ( n ) r * ( n + N ) ) Re ( Σ n = 0 N g - 1 r ( n ) r * ( n + N ) ) ]
其中,Ng为GI的长度。
2)整数频偏的估计方法
在DVB-T***中,除了小数频偏之外,还存在着整数倍于子载波间隔的频偏,这个频偏在DVB-T***中是根据插在OFDM符号上的连续导频来进行估计的,即在FFT之后的频域进行的。
首先,用前后两个相邻的OFDM符号做相关,导频位置的相关值比较大,而非导频位置的相关近似于零。因此,可用频域最大似然估计理论估计出连续导频的位置,进而得到整数倍频偏。
定义整数倍频偏估计的相关函数:
C ( m ) = | Σ k = 0 N cp - 1 Y l - 1 * ( p k + m ) N Y l ( p k + m ) N |
其中,pk,Ncp分别为连续导频对应的子载波的位置和连续导频的数目,m为循环滑动量。
从上式可以看出,整数倍频偏为mt时,C(mt)可以认为是各连续导频点的相关累加和。因此,整数频偏的估计值可以由下式得出:
Δ f ^ I = arg | max [ C ( m ) ] | .
在DVB-T2***的帧结构中,连续导频只存在于data符号中,而且对于不同的离散导频模式,连续导频的位置也不一样。而离散导频的模式只有把信令解出来之后才能得到,但是,如果***中存在着整数频偏,就不可能把信令解出来,这就出现一个死循环。
因此,在DVB-T2***中,就不可能利用data符号的连续导频来估计整数频偏。
由于P2符号的导频位置不变,在1K~8K模式下,由于P2符号的个数大于2,前后两个P2符号可以做相关,因此可以将P2符号的导频当作连续导频来估计整数频偏。但是,在16K和32K模式下,只有一个P2符号,由于P2符号的导频密度比较大,其位置涵盖了后续的data符号的离散导频的位置。因此P2符号只能和后续的data符号做相关,然后提取data符号中离散导频位置的信息来处理,但是data符号的离散导频有八种,由于这种相关运算过程是在不知道哪一种离散导频的情况下来进行的,因此只能重复做8次这样的相关运算,从而使得计算量大大加大。
频偏估计模块是在帧同步和符号同步完成的基础上进行的。
在***启动之后,首先需要进行帧同步,即找到DVB-T2***的帧(Frame)的起始位置。在图1中的DVB-T2的帧结构中我们可以看到,每个T2帧的开头都有一个P1符号,这个P1符号具有很好的相关性。应用滑动自相关的方法,可以找到P1的起始位置,也就是T2帧的起始位置,关于帧同步的方法有很多关于这方面的专利,例如US2010/0284498,在参考文献[1]中也有这方面的介绍。需要注意的是,在帧同步的过程中,会检测到粗的频偏,会将***的残留频偏范围缩小到2个1K FFT的子载波频率间隔之内。
帧同步结束之后,由于存在一定的误差,只能找到T2帧的大概的起始位置。为了后续的正确接收信号,必须进一步的锁定帧起始的位置,这个过程就是符号同步。符号同步同样有很多种方法可以应用,在参考文献[1]中,也给出了符号同步的过程。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种DVB-T2***中频偏估计的方法,提高小数频偏估计的准确性,避免整数频偏16K,32K模式下复杂的相关运算,提高整数频偏的检测速度和估计精度。
为解决上述技术问题,本发明所采用的技术方案是:一种DVB-T2***中频偏估计的方法,该方法包括残留频偏估计、整数频偏估计和小数频偏估计三个部分:
所述残留频偏估计过程为:
1)根据符号同步的结果,取接收信号中的P1符号部分的信号,rp1(n),n=[0~1023];
2)根据帧同步的结果S1和S2,生成接收机在当前传输模式下的P1符号,sp1(n);
3)设置残留频偏的精度为1/M,相关次数为2M+1,其中M的取值范围为1
≤M≤32;
4)对接收信号中的P1符号部分的信号rp1(n)做频偏纠正,得到k个序列
Figure BDA00003338275500041
r ~ pl ( n , k ) = r pl ( n ) * exp ( - j * 2 π * k * 1 M * n * 1 1024 )
其中,n=0~1023,k=[-M,-(M-1),-(M-2)......-1,0,1,2,......M];
5)对sp1(n)和
Figure BDA00003338275500051
进行相关运算,得到运算结果c(k):
c ( k ) = Σ n = 0 1023 r ~ pl ( n , k ) × sp 1 * ( n ) ,
其中,
Figure BDA00003338275500053
为sp1(n)的共轭;
6)找到c(k)的2M+1个结果中的最大值cmax,cmax对应的角度就是粗频偏估计值
Figure BDA00003338275500054
Δ f ~ 0 = arg | msx [ c ( k ) ] | ;
所述整数频偏估计过程为:
1)根据符号同步的结果,取接收信号的第一个P2符号rp2(s),然后纠正上述粗频偏估计值
Figure BDA00003338275500056
得到序列
Figure BDA00003338275500057
r ~ p 2 ( s ) = r p 2 ( s ) * exp ( - j * 2 π * Δ f ~ 0 * s * 1 1024 ) ;
其中,s=0,1,2,…N-1,N为当前***FFT的长度
2)生成P2符号导频的时域信息sp2(s);
3)确定整数频偏的范围q:
q = N 1024 * 2 * 1 M
4)分别设置q个不同的频偏,对
Figure BDA000033382755000510
做频偏纠正,得到q个序列
Figure BDA000033382755000511
r - p 2 ( s , p ) = r ~ p 2 ( s ) * exp ( - j * 2 π * s * p * 1 N )
其中,s=0~N-1,p=[-q/2,q/2],
5)对P2序列sp2(s)和
Figure BDA000033382755000513
进行相关运算,得到运算结果c2(p):
c 2 ( p ) = Σ s = 0 N - 1 r - p 2 ( S , P ) * S * p 2 ( S ) ,
其中,s* p2(s)是sp2(s)的共轭;
6)找到c2(p)的q个结果的最大值c2max,c2max对应的p就是整数频偏估计值
Figure BDA000033382755000515
Δ f ~ 1 = arg | max [ c 2 ( p ) ] × 1024 N |
所述小数频偏估计过程为:
1)根据符号同步的结果,从第一个P2符号的保护间隔之内某一个位置t1开始,取一个完整的P2符号rp21(s),并利用粗频偏估计值
Figure BDA00003338275500061
和整数频偏估计值
Figure BDA00003338275500062
之和
Figure BDA00003338275500063
得到序列
Figure BDA00003338275500064
r ~ p 21 ( S ) = r p 21 ( S ) * exp ( - j * 2 π * ( Δ f ~ 0 + Δ f ~ 1 ) * S * 1 N )
其中所述保护间隔长度为GI*N,N为OFDM符号的长度,即当前***FFT的长度,GI为[1/4,1/8,1/16,1/32,1/128,19/128,19/256]中的一种;
2)对sp2(s)和
Figure BDA00003338275500066
进行相关运算,得到运算结果c31
c 31 = Σ s = 0 N - 1 r ~ p 21 ( s ) * s * p 2 ( s )
其中,s=0~N-1,s* p2(s)是sp2(s)的共轭;
3)根据符号同步的结果,从第一个P2符号的保护间隔之内另一个位置t2开始,取一个完整的P2符号rp22(s),然后纠正上述估计的频偏值
Figure BDA00003338275500068
得到序列
r ~ p 22 ( S ) = r p 22 ( S ) * exp ( - j * 2 π * ( Δ f ~ 0 + Δ f ~ 1 ) * S * 1 N )
4)对sp2(s)和
Figure BDA000033382755000611
进行相关运算,得到运算结果c32
c 32 = Σ s = 0 N - 1 r ~ p 22 ( s ) * s * p 2 ( s )
5)通过上述两个相关结果c31和c32计算小数频偏估计值
Figure BDA000033382755000613
Δ f ~ 2 = N 2 π ( t 2 - t 1 ) tan - 1 ( c 32 c 31 ) × 1024 N
最后,将上述粗频偏估计值
Figure BDA000033382755000615
整数频偏估计值和小数频偏估计值
Figure BDA000033382755000617
合并,得到最终的频偏
Figure BDA000033382755000618
Figure BDA000033382755000619
在DVB-T2***中,GI有七种可能,分别是[1/4,1/8,1/16,1/32,1/128,19/128,19/256],当符号同步和模式检测之后,GI的比例也就确定了,以下叙述中的GI的含义表示上述确定的比例,而保护间隔的长度为GI*N,N为OFDM符号的长度。
与现有技术相比,本发明所具有的有益效果为:本发明使用接收机中生成的P1符号,和接收信号做相关,减小残留频偏的范围;可以通过计算多个小数频偏然后求均值的方法提高估计的准确性;本发明避免了整数频偏16K,32K模式下复杂的相关运算,提高了整数频偏的检测速度和估计精度。
附图说明
图1为DVB-T2帧结构;
图2为DVB-T2P1符号结构;
图3为DVB-T2生成P1符号流程图;
图4为DVB-T2P2符号的导频图案(MISO模式);
图5为DVB-T2P2符号的导频图案(SISO模式);
图6为生成加扰序列的PRBS生成器示意图;
图7为生成p2_prbs0(v)的PRBS产生器;
图8为本发明一实施例残留频偏估计流程图;
图9为本发明一实施例整数频偏估计流程图。
具体实施方式
以下结合附图进一步说明本发明的原理。
1)残留频偏的估计方法
a)根据符号同步的结果,取接收信号的P1符号rp1(n)
b)根据帧同步的结果S1和S2,生成P1符号sp1(n)
P1符号的产生过程如下:
由于帧同步锁定之后,会提供两个参数S1和S2,这两个参数指示了当前***采用的FFT模式以及MISO模式,它决定了DVB-T2***的帧头P1符号的内容和结构。
表1S1和S2调制模式
Figure BDA00003338275500081
从上表可以看到,S1和S2分别是3bit和4bit的参数,分别对应着8和16个可能的序列。其中S1对应的序列长度为64bit,S2对应的序列长度为256bit。下面介绍P1序列的生成过程,详见参考文献[2]。
第一,根据S1和S2选择出调制序列:
{mss_seq0…mss_seq383}={cssS1,cssS2,cssS1}
={cssS1,0…cssS1,63,cssS2,0…cssS2,255,cssS1,0…cssS1,63}
其中,mss_seq是长度为384bit的序列,经过后续的变换,调制到1K FFT的相应的384个子载波位置上。它是由S1和S2序列产生的。S1是3bit的信息,有8种选择,分别对应着上面表格中的S1的每一行64个bit信息cssS1={cssS1,0…cssS1,63}。S2是4bit的信息,有16种选择,分别对应着上表中S2的每一行255bit的信息cssS2={cssS2,0…cssS2,255}。
第二,对上面的调制信息做DBPSK调制:
mss_diff=DBPSK(mss_seq)
mss _ diff i = mss _ diff i - 1 , mss _ seq i = 0 - mss _ diff i = 1 , mss _ seq i = 1
mss_diff-1=+1
其中,mss_diff是mss_seq经过DBPSK调制之后的长度为384bit的信息;mss_diff-1表示DBPSK的初始相位;mss_diffi表示当前正在调制的bit的调制序列;
第三,对上面的序列加扰:
mss_scr=SCRAMBLING{mss_diff}
加扰序列由图6的PRBS生成器产生,长度为384bit:
加扰PRBS生成器的生成多项式为:
1+x14+x15;其中x为0或1;
这样,经过加扰之后的序列可以用下面的公式表示:
mss_scri=mss_diffi×(1-2PRBSi);
其中,PRBSi是当前时刻PRBS生成器输出的bit(0或者1),对应着DBPSK结果中的当前正在做加扰处理的bit,i=[0~383]。
第四,将加扰序列mss_scri按照固定的位置分配到1K模式的子载波上。
第五,做IFFT,得到P1符号。
第六,对上面得到的P1符号的I/Q两路分别取符号位,得到本发明需
要的P1序列:
Figure BDA00003338275500095
其中
Figure BDA00003338275500092
表示取P1的实部,
Figure BDA00003338275500093
为P1的虚部。
c)由于残留的频偏在+1,-1之间,因此需要进一步的缩小范围,设置精
度为1/4,相关的次数为9,此处,我们设定精度为1/4,所以需要做9次相关。但是,精度并不局限于1/4。需要说明的是,这个精度范围表示在残留频偏估计阶段能够估计的最小的频偏的大小,更加精细的频偏估计是基于此来做的。如果此处的精度设置过大,那么后续的精细频偏估计需要搜索的范围就比较大。而在这个步骤中,相关序列的长度为1K,运算量相对小。而在精细频偏估计的步骤中,相关序列的长度最大可能为32768,运算量比较大。为了减小后续模块的运算量,在此处设置的精度要尽量的小。但是,由于1K的P1符号的分辨率有限,所以精度并不是越小越好,仿真显示,1/4是一个比较合适的数字。
d)若精度为1/4,那么分别设置9个不同的频偏分别是
对rp1(n)做频偏纠正,得到
r ~ pl ( n , k ) = r pl ( n ) * exp ( - j * 2 π * k * 1 4 * n * 1 1024 )
其中,n=0~1023,k=[-4,-3,-2,-1,0,1,2,3,4]
e)对P1序列sp1(n)和接收的P1序列
Figure BDA00003338275500104
进行相关运算,得到运算结果c(k)
c ( k ) = Σ n = 0 1023 r ~ pl ( n , k ) × sp 1 * ( n )
其中,n=0~1023,k=[-4,-3,-2,-1,0,1,2,3,4],
sp1 *(n)为sp1(n)的共轭;
f)比较c(k)的9个结果的最大值,其对应的频偏值就是粗频偏估计值
Figure BDA00003338275500106
Δ f ~ 0 = arg | msx [ c ( k ) ] |
这样,经过上述的步骤之后,可以将***的残留频偏控制在±1/4个1KFFT子载波的间隔之内。
2)整数频偏估计方法
a)根据符号同步的结果,取接收信号的第一个P2符号rp2(s),然后纠正上
述估计的频偏值得到序列
Figure BDA00003338275500112
r ~ p 2 ( s ) = r p 2 ( s ) * exp ( - j * 2 π * Δ f ~ 0 * s * 1 1024 )
s=[0,N-1]
b)生成P2符号导频的时域信息sp2(s)
i.根据图7所示的PRBS产生器生成长度为kmax个bit信息p2_prbs0(v)。kmax为当前***P2符号的有效子载波的个数;
图7所示的PRBS产生器的生成多项式为:
1+x2+x11
ii.将PRBS序列p2_prbs0(v)的0和1分别变成+1和-1
p 2 _ prbsl ( v ) = + 1 , p 2 _ prbs 0 ( v ) = 0 - 1 , p 2 _ prbs 0 ( v ) = - 1
其中,0≤v≤kmax,kmax为有效子载波的个数;
iii.将PRBS序列p2_prbs1(v)中序号(index)不是3的倍数的位置上的bit置成0:
Figure BDA00003338275500115
其中,0≤v≤kmax
iv.将上述的变换后的PRBS序列p2_prbs2(w)进行循环移位:
Figure BDA00003338275500116
其中,0≤w≤FFT_SIZE
v.将p2_prbs3(w)做IFFT变换:
p 2 _ prbs 4 ( w ) = 1 N Σ i = 0 N - 1 p 2 _ prbs 3 ( k ) × e j 2 πwi N ,
vi.对p2_prbs4(w)的I/Q两路分别取符号位,得到P2相关序列sp2(s):
Figure BDA00003338275500121
其中
Figure BDA00003338275500122
Figure BDA00003338275500123
表示取p2_prbs4的实部,
Figure BDA00003338275500124
为p2_prbs4的虚部。
c)确定整数频偏的范围
q = N 1024 * 2 * 1 M
其中N为当前***FFT的长度。
d)分别设置q个不同的频偏,对
Figure BDA00003338275500127
做频偏纠正,得到q个序列
Figure BDA00003338275500128
r - p 2 ( s , p ) = r ~ p 2 ( s ) * exp ( - j * 2 π * s * p * 1 N )
其中,s=0~N-1,p=[-q/2,q/2]
e)分别对P2序列sp2(s)和进行相关运算,得到q个运算结果c2(p)
c 2 ( p ) = Σ s = 0 N - 1 r ~ p 2 ( S , P ) * S * p 2 ( S ) ,
其中,s=0~N-1,p=[-q/2,q/2],s* p2(s)是sp2(s)的共轭
f)比较c2(p)的q个结果的最大值,其对应的频偏值就是整数频偏估计值
Figure BDA000033382755001217
Δ f ~ 1 = arg | max [ c 2 ( p ) ] × 1024 N |
3)小数频偏的估计方法
a)根据符号同步的结果,从第一个P2符号的GI之内某一个位置t1开始,取一个完整的P2符号rp21(n),然后纠正上述估计的频偏值
Figure BDA000033382755001213
得到序列
Figure BDA000033382755001214
r ~ p 21 ( S ) = r p 21 ( S ) * exp ( - j * 2 π * ( Δ f ~ 0 + Δ f ~ 1 ) * S * 1 N )
b)对sp2(s)和
Figure BDA000033382755001216
进行相关运算,得到运算结果c31
c 31 = Σ s = 0 N - 1 r ~ p 21 ( s ) * s * p 2 ( s )
其中,s=0~N-1
c)根据符号同步的结果,从第一个P2符号的GI之内另一个位置t2开始,取一个完整的P2符号rp22(s),然后纠正上述估计的频偏值
Figure BDA00003338275500132
得到序列
Figure BDA00003338275500133
d)对sp2(s)和
Figure BDA00003338275500134
进行相关运算,得到运算结果c32
c 31 = Σ s = 0 N - 1 r ~ p 22 ( s ) * s * p 2 ( s )
其中,s=0~N-1
e)通过上述两个相关结果计算小数频偏
Figure BDA00003338275500136
Δ f ~ 2 = N 2 π ( t 2 - t 1 ) tan - 1 ( c 32 c 31 ) × 1024 N
4)计算最终的频偏
Figure BDA00003338275500138
Δ f ~ = Δ f ~ 0 + Δ f ~ 1 + Δ f ~ 2 .
残留频偏估计流程图见图8。
整数频偏流程图见图9,本发明充分利用了P2导频信号的时域信息(简称为P2序列)来进行相关运算,从而得到整数频偏。
符号同步结束之后,我们可以找到非常精确的OFDM符号的起始位置。以此位置开始,取一个FFT窗口的接收数据,此时对应于第一个P2符号。将接收到的P2符号和本地存储的P2序列做互相关,因为在接收信号中存在着频偏,所以直接相关的结果中看不到峰值。为了估计残留的频偏,我们在相关之前先在接收信号中将之前计算的频偏都纠正掉,然后利用P2符号的相关特性通过两个窗口的相关运算推到计算出剩余的小数频偏。
本发明充分利用了P2的导频信号的时域信息,进行相关运算,从而得到小数频偏。
在符号同步结束之后,我们可以找到非常精确的OFDM符号的起始位置,以此位置开始取出第一个P2符号,在进行小数频偏估计之前,要将之前估计的所有频偏都纠正掉,然后再和本地存储的P2序列做互相关,其结果必然会是一个峰值,峰值的相位中就隐含了小数频偏的信息。
如果我们在循环前缀内的任一位置再取一个FFT窗口的信息,然后再和P2序列的相应循环移位序列做互相关,也会得到一个峰值。通过计算两个峰值的相位差就得到了小数频偏,如果我们重复多做几次这样的操作,然后将得到的小数频偏做加权平均的话,得到的小数频偏会更加准确。数学描述如下:
假设我们接收的信号可以表示成:
r ( u ) = { cons ( u ) + pilot ( u ) + noise ( u ) } × e j 2 πϵu N
其中ε是小数频偏,noise(u)是噪声分量。cons(u)是P2符号中星座图部分的时域信息,pilot(u)是P2序列。
r(u)和pilot(u)的互相关函数可以表示成:
corr ( z ) = Σ i = z 2 + N - 1 r ( i ) × pilot * ( i )
z是相关窗口的起点位置。
corr ( z ) = Σ i = z 2 + N - 1 r ( i ) × pilot * ( i )
= Σ i = z z + N - 1 { cons ( i ) + pilot ( i ) + noise ( i ) } × e j 2 πϵi N × pilot * ( i )
= Σ i = z z + N - 1 { | pilot ( i ) | 2 + noise ′ ( i ) } × e j 2 πϵi N
为了方便起见,我们只取z=0和z=GI*N-1两个窗口进行互相关。
corr ( 0 ) = Σ i = 0 N - 1 { | pilot ( i ) | 2 + noise ′ ( i ) } × e j 2 πϵi N
corr ( GI * N - 1 ) = Σ i = GI * N - 1 GI * N - 1 + N - 1 { | pilot ( i ) | 2 + noise ′ ( i ) } × e j 2 πϵi N
= corr ( 0 ) × e j 2 πϵGI * N N
因此我们可以通过这两个相关峰corr(0)和corr(GI*N-1)来得到小数频偏的估计值ε:
Figure BDA00003338275500154
ϵ = 1 2 π × GI × tan - 1 ( corr ( GI * N - 1 ) corr ( 0 ) )
其中GI为当前***参数中保护间隔的比例,有七种取值可能,分别是
Figure BDA00003338275500156
在这里,选取的两个窗口并不局限于z=0和z=GI*N-1。任意两个窗口都可适用做本发明所提出的小数频偏估计方法。
参考文献:
[1]Digital Video Broadcasting(DVB);Implementation guidelines for a secondgeneration digital terrestrial television broadcasting system,DVB Document A133,June2010
[2]ETSI EN302755V1.2.1,Digital Video Broadcasting(DVB);Frame structurechannel coding and modulation for a second generation digital terrestrial televisionbroadcasting system(DVB-T2).

Claims (3)

1.一种DVB-T2***中频偏估计的方法,其特征在于,该方法包括残留频偏估计、整数频偏估计和小数频偏估计三个部分:
所述残留频偏估计过程为:
1)根据符号同步的结果,取接收信号中的P1符号部分的信号,rp1(n),n=[0~1023];
2)根据帧同步的结果S1和S2,生成接收机在当前传输模式下的P1符号,sp1(n);
3)设置残留频偏的精度为1/M,相关次数为2M+1,其中M的取值范围为1≤M≤32;
4)对接收信号中的P1符号部分的信号rp1(n)做频偏纠正,得到k个序列
Figure FDA00003338275400015
r ~ pl ( n , k ) = r pl ( n ) * exp ( - j * 2 π * k * 1 M * n * 1 1024 )
其中,n=0~1023,k=[-M,-(M-1),-(M-2)......-1,0,1,2,......M];
5)对sp1(n)和
Figure FDA00003338275400016
进行相关运算,得到运算结果c(k):
c ( k ) = Σ n = 0 1023 r ~ pl ( n , k ) × sp 1 * ( n ) ,
其中,sp1 *(n)为sp1(n)的共轭;
6)找到c(k)的2M+1个结果中的最大值cmax,cmax对应的角度就是粗频偏估计值
Figure FDA00003338275400013
Δ f ~ 0 = arg | msx [ c ( k ) ] | ;
所述整数频偏估计过程为:
1)根据符号同步的结果,取接收信号的第一个P2符号rp2(s),然后纠正上述粗频偏估计值
Figure FDA00003338275400021
得到序列
Figure FDA00003338275400022
r ~ p 2 ( s ) = r p 2 ( s ) * exp ( - j * 2 π * Δ f ~ 0 * s * 1 1024 ) ;
其中,s=0,1,2,…N-1,N为当前***FFT的长度;
2)生成P2符号导频的时域信息sp2(s);
3)确定整数频偏的范围q:
q = N 1024 * 2 * 1 M
4)分别设置q个不同的频偏,对
Figure FDA00003338275400025
做频偏纠正,得到q个序列
Figure FDA00003338275400026
r - p 2 ( s , p ) = r ~ p 2 ( s ) * exp ( - j * 2 π * s * p * 1 N )
其中,s=0~N-1,p=[-q/2,q/2],
5)对P2序列sp2(s)和
Figure FDA00003338275400028
进行相关运算,得到运算结果c2(p):
c 2 ( p ) = Σ s = 0 N - 1 r - p 2 ( s , p ) * s * p 2 ( s ) ,
其中,s* p2(s)是sp2(s)的共轭;
6)找到c2(p)的q个结果的最大值c2max,c2max对应的p就是整数频偏估计值
Δ f ~ 1 = arg | max [ c 2 ( p ) ] × 1024 N |
所述小数频偏估计过程为:
1)根据符号同步的结果,从第一个P2符号的保护间隔之内某一个位置t1开
始,取一个完整的P2符号rp21(s),并利用粗频偏估计值
Figure FDA000033382754000212
和整数频偏
估计值
Figure FDA000033382754000213
之和
Figure FDA000033382754000214
得到序列
Figure FDA000033382754000215
r ~ p 21 ( S ) = r p 21 ( S ) * exp ( - j * 2 π * ( Δ f ~ 0 + Δ f ~ 1 ) * S * 1 N )
其中所述保护间隔长度为GI*N,N为OFDM符号的长度,即当前***FFT的长度,GI为[1/4,1/8,1/16,1/32,1/128,19/128,19/256]中的一种;
2)对sp2(s)和
Figure FDA000033382754000217
进行相关运算,得到运算结果c31
c 31 = Σ s = 0 N - 1 r ~ p 21 ( s ) * s * p 2 ( s )
其中,s=0~N-1,s* p2(s)是sp2(s)的共轭;
3)根据符号同步的结果,从第一个P2符号的保护间隔之内另一个位置t2开始,取一个完整的P2符号rp22(s),然后纠正上述估计的频偏值
Figure FDA00003338275400032
得到序列
r ~ p 22 ( S ) = r p 22 ( S ) * exp ( - j * 2 π * ( Δ f ~ 0 + Δ f ~ 1 ) * S * 1 N )
4)对sp2(s)和
Figure FDA00003338275400034
进行相关运算,得到运算结果c32
c 32 = Σ s = 0 N - 1 r ~ p 22 ( s ) * s * p 2 ( s )
5)通过上述两个相关结果c31和c32计算小数频偏估计值
Figure FDA00003338275400036
Δ f ~ 2 = N 2 π ( t 2 - t 1 ) tan - 1 ( c 32 c 31 ) × 1024 N
最后,将上述粗频偏估计值
Figure FDA00003338275400038
整数频偏估计值
Figure FDA00003338275400039
和小数频偏估计值合并,得到最终的频偏
Figure FDA000033382754000311
Figure FDA000033382754000312
2.根据权利要求1所述的DVB-T2***中频偏估计的方法,其特征在于,所述残留频偏估计过程步骤3)中,M=4。
3.根据权利要求1所述的DVB-T2***中频偏估计的方法,其特征在于,所述整数频偏估计过程步骤2)中,生成P2符号导频的时域信息sp2(v)的方法为:
a)根据PRBS产生器生成长度为kmax个bit信息的序列p2_prbs0(v),kmax为当前***P2符号的有效子载波的个数,所述PRBS产生器的生成多项式为1+x2+x11,其中x为0或1;
b)将PRBS序列p2_prbs0(v)的0和1分别变成+1和-1:
p 2 _ prbsl ( v ) = + 1 , p 2 _ prbs 0 ( v ) = 0 - 1 , p 2 _ prbs 0 ( v ) = - 1 ,
其中,0≤v≤kmax;
c)将PRBS序列p2_prbs1(v)中序号不是3的倍数的位置上的bit置成0:
Figure FDA00003338275400041
其中,0≤v≤kmax;
d)将上述的变换后的PRBS序列p2_prbs2(w)进行转置,得到p2_prbs3(w):
Figure FDA00003338275400042
其中,0≤w≤N-1;
将p2_prbs3(w)做IFFT变换,得到p2_prbs4(w):
p 2 _ prbs 4 ( w ) = 1 N Σ i = 0 N - 1 p 2 _ prbs 3 ( k ) × e j 2 πwi N ,
其中,N为FFT的长度,w=[0~N-1];
e)对p2_prbs4(w)的I/Q两路分别取符号位,得到P2符号导频的时域信息
sp2(s):
Figure FDA00003338275400044
其中
Figure FDA00003338275400046
表示取p2_prbs4的实部,
Figure FDA00003338275400047
为p2_prbs4的虚部。
CN201310232900.4A 2013-06-13 2013-06-13 一种dvb-t2***中频偏估计的方法 Active CN103281282B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310232900.4A CN103281282B (zh) 2013-06-13 2013-06-13 一种dvb-t2***中频偏估计的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310232900.4A CN103281282B (zh) 2013-06-13 2013-06-13 一种dvb-t2***中频偏估计的方法

Publications (2)

Publication Number Publication Date
CN103281282A true CN103281282A (zh) 2013-09-04
CN103281282B CN103281282B (zh) 2016-04-06

Family

ID=49063734

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310232900.4A Active CN103281282B (zh) 2013-06-13 2013-06-13 一种dvb-t2***中频偏估计的方法

Country Status (1)

Country Link
CN (1) CN103281282B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104469236A (zh) * 2013-09-25 2015-03-25 扬智科技股份有限公司 采样频偏的纠正装置与纠正方法
CN105282080A (zh) * 2014-06-26 2016-01-27 无锡明波微电子技术有限公司 Ofdm符号同步位置搜索方法及其同步装置
CN108270712A (zh) * 2017-12-08 2018-07-10 中国电子科技集团公司第三十研究所 一种改进的电力载波同步检测方法
CN110505172A (zh) * 2019-07-12 2019-11-26 四川安迪科技实业有限公司 一种基于离散导频辅助的频偏估计方法及装置
CN110661590A (zh) * 2015-03-09 2020-01-07 第一媒体有限责任公司 ***发现与信令
CN113179233A (zh) * 2021-04-01 2021-07-27 成都星联芯通科技有限公司 一种对抗低轨宽带卫星***大频偏大时延的联合估计方法
US11923966B2 (en) 2014-08-25 2024-03-05 ONE Media, LLC Dynamic configuration of a flexible orthogonal frequency division multiplexing PHY transport data frame preamble

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070019763A1 (en) * 2005-07-19 2007-01-25 Samsung Electronics Co., Ltd. Sampling frequency offset estimation apparatus and method for OFDM system
CN101212440A (zh) * 2006-12-29 2008-07-02 大唐移动通信设备有限公司 一种载波频率偏移的估计方法及装置
EP2134044A1 (en) * 2008-06-04 2009-12-16 Sony Corporation New frame and signalling pattern structure for multi-carrier systems
CN102377711A (zh) * 2010-08-11 2012-03-14 中国科学院微电子研究所 一种ofdm无线终端的同步方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070019763A1 (en) * 2005-07-19 2007-01-25 Samsung Electronics Co., Ltd. Sampling frequency offset estimation apparatus and method for OFDM system
CN101212440A (zh) * 2006-12-29 2008-07-02 大唐移动通信设备有限公司 一种载波频率偏移的估计方法及装置
EP2134044A1 (en) * 2008-06-04 2009-12-16 Sony Corporation New frame and signalling pattern structure for multi-carrier systems
CN102377711A (zh) * 2010-08-11 2012-03-14 中国科学院微电子研究所 一种ofdm无线终端的同步方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张博等: "DVB - T2 中一种基于P1 符号的改进同步算法", 《信息技术》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104469236B (zh) * 2013-09-25 2017-11-07 扬智科技股份有限公司 采样频偏的纠正装置与纠正方法
CN104469236A (zh) * 2013-09-25 2015-03-25 扬智科技股份有限公司 采样频偏的纠正装置与纠正方法
CN105282080A (zh) * 2014-06-26 2016-01-27 无锡明波微电子技术有限公司 Ofdm符号同步位置搜索方法及其同步装置
CN105282080B (zh) * 2014-06-26 2018-08-07 无锡明波微电子技术有限公司 Ofdm符号同步位置搜索方法及其同步装置
US11923966B2 (en) 2014-08-25 2024-03-05 ONE Media, LLC Dynamic configuration of a flexible orthogonal frequency division multiplexing PHY transport data frame preamble
US11627030B2 (en) 2015-03-09 2023-04-11 ONE Media, LLC System discovery and signaling
CN110661590A (zh) * 2015-03-09 2020-01-07 第一媒体有限责任公司 ***发现与信令
CN110661590B (zh) * 2015-03-09 2021-04-20 第一媒体有限责任公司 ***发现与信令
US11012282B2 (en) 2015-03-09 2021-05-18 ONE Media, LLC System discovery and signaling
CN108270712A (zh) * 2017-12-08 2018-07-10 中国电子科技集团公司第三十研究所 一种改进的电力载波同步检测方法
CN110505172A (zh) * 2019-07-12 2019-11-26 四川安迪科技实业有限公司 一种基于离散导频辅助的频偏估计方法及装置
CN113179233B (zh) * 2021-04-01 2022-11-04 成都星联芯通科技有限公司 一种对抗低轨宽带卫星***大频偏大时延的联合估计方法
CN113179233A (zh) * 2021-04-01 2021-07-27 成都星联芯通科技有限公司 一种对抗低轨宽带卫星***大频偏大时延的联合估计方法

Also Published As

Publication number Publication date
CN103281282B (zh) 2016-04-06

Similar Documents

Publication Publication Date Title
CN103281282B (zh) 一种dvb-t2***中频偏估计的方法
CN107154908B (zh) 前导符号的生成方法
US8422595B2 (en) Channel estimation for communication systems with multiple transmit antennas
US8548101B2 (en) Demodulation method and demodulator for orthogonal frequency multiplexing—multiple input multiple output system
CN105282076B (zh) 前导符号的生成方法及频域ofdm符号的生成方法
CN112039810B (zh) 频域ofdm符号的生成方法及前导符号的生成方法
CN103532898A (zh) 基于cazac序列的ofdm训练符号生成及同步方法
KR20160064823A (ko) 수신 신호의 주파수 오프셋을 추정하는 방법 및 장치
Nasir et al. Performance of coarse and fine timing synchronization in OFDM receivers
US20170265202A1 (en) Time domain pilot of single-carrier mimo system and synchronization method thereof
CN105024952A (zh) 频域ofdm符号的生成方法及前导符号的生成方法
CN101534287A (zh) 移动通信***中进行载波频率偏移校正的方法和装置
CN1787509A (zh) 一种用于WiMAX***基站接收端的整数倍频偏校正方法
CN100355255C (zh) 基于统计平均的大搜索范围ofdm***载波同步方法
CN104901918A (zh) 基于Chirp信号产生OFDM数据序列的方法及同步方法
CN103501284A (zh) 一种dvb-t2***中精细定时同步的方法
US7830990B2 (en) Method for estimating and compensating frequency offset and frequency offset estimation module
CN103546416A (zh) Ofdm符号同步方法
CN103346990B (zh) 一种dvb-t2***中符号定时校正的方法
CN102137056B (zh) 伪噪声序列相位检测的方法
Wang et al. Estimation of carrier frequency offset and channel state information of generalize frequency division multiplexing systems by using a Zadoff–Chu sequence
CN107835141A (zh) 自相关与互相关结合的多段重复序列ofdm同步算法
CN114401178B (zh) 频偏估计方法、装置、设备及存储介质
CN102868663B (zh) 多载波数字信息传输***中采样时钟同步方法及***
CN105262710B (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
CB02 Change of applicant information

Address after: 410125 Hunan, Changsha economic and Technological Development Zone, the east side of the south section of the No. ten road, Tong Tong Street, No.

Applicant after: GOKE MICROELECTRONICS CO., LTD.

Address before: 410125 Hunan province Changsha City Economic and Technological Development Zone Li Xiang Road No. 98

Applicant before: Hunan Guoke Microelectronics Co., Ltd.

COR Change of bibliographic data
C14 Grant of patent or utility model
GR01 Patent grant