CN101567868A - 一种cmmb制式手机电视小数频偏估计方法 - Google Patents

一种cmmb制式手机电视小数频偏估计方法 Download PDF

Info

Publication number
CN101567868A
CN101567868A CNA2008100667895A CN200810066789A CN101567868A CN 101567868 A CN101567868 A CN 101567868A CN A2008100667895 A CNA2008100667895 A CN A2008100667895A CN 200810066789 A CN200810066789 A CN 200810066789A CN 101567868 A CN101567868 A CN 101567868A
Authority
CN
China
Prior art keywords
mrow
msub
msubsup
math
preprocessing
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.)
Pending
Application number
CNA2008100667895A
Other languages
English (en)
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.)
ZTE Corp
Original Assignee
ZTE Corp
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 ZTE Corp filed Critical ZTE Corp
Priority to CNA2008100667895A priority Critical patent/CN101567868A/zh
Publication of CN101567868A publication Critical patent/CN101567868A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Noise Elimination (AREA)

Abstract

本发明公开了一种CMMB制式手机电视小数频偏估计方法。该方法将用于运算的数据先通过预处理步骤进行处理后,对其进行相关累加,使用相关累加值的相位求取小数频偏估计的结果。所述预处理为r(n)-r(n-M)形式的差分运算,L为重复符号长度,n=0,1,…,L-1,M为小于L-1的正整数;或者为由FIR形式的频域特性为高通性质的滤波器处理,或者为由FIR形式的频域特性为带通性质的直流分量不在其通带内的滤波器处理。本发明增加了对进行小数频偏估计的数据进行预处理的步骤,预处理后再使用现有方法进行小数频偏估计,由于预处理的作用本发明的估计结果不受数据中所含直流分量的影响,较现有方法大幅提高了小数频偏估计的鲁棒性。

Description

一种CMMB制式手机电视小数频偏估计方法
技术领域
本发明属于通信与信息技术中载波频偏估计和补偿领域,特别涉及一种针对CMMB(China Mobile Multimedia Broadcasting,中国移动数字多媒体广播)制式手机电视标准终端设备的小数载波频偏的估计方法。
背景技术
国家广电总局于06年10月颁布了中国移动多媒体广播行业标准,确定采用我国自主研发的移动电视接收标准STiMi,该标准从06年11月1日起实施。它是国内自主研发的第一套面向手机、PDA(Personal Digital Assistant,个人数字助理)、MP3、MP4、数码相机、笔记本电脑多种移动终端的系统,利用S波段卫星信号实现“天地”一体覆盖、全国漫游,支持25套电视节目和30套广播节目。CMMB规定了在广播业务频率范围内,移动多媒体广播系统广播信道传输信号的帧结构、信道编码和调制,该标准适用于30MHz到3000MHz频率范围内的广播业务频率,通过卫星和/或地面无线发射电视、广播、数据信息等多媒体信号的广播系统,可以实现全国漫游。
CMMB标准采用了OFDM(Orthogonal Frequency Division Multiplexing,正交频分复用)技术,它是一种公知的多载波调制技术,其主要原理是:将信道分成若干正交子信道,将高速数据转换成并行的低速子数据流,调制到每个子信道上进行传输。正交信号可以通过在接收端采用相关技术来分开,这样可以减少子信道之间的相互干扰。每个子信道上的信号带宽小于信道的相关带宽,因此每个子信道可以看成平坦性衰落,从而可以消除符号间干扰。而且由于每个子信道带宽仅仅是原信道带宽的一小部分,信道均衡变得相对容易。OFDM目前已被用于数种无线系统标准中,譬如欧洲数字音频和数字视频广播系统(DAB、DVB-T、DVB-H)、5GHz高数据数率无线LAN(IEEE802.11a,HiperLan2,MMAC)系统等。
由于OFDM各子信道带宽较小,对载波频率偏差的敏感程度很高,需要非常精确的频率同步。根据频率偏差和子载波间隔之间的关系可以将其分为小数部分和整数部分,整数部分仅使信息信号在子信道上平移,并不破坏各子载波间的正交性;而小数部分会造成子信道干扰,破坏子载波间的正交性,导致系统的误码率性能严重下降。因此,在接收端先进行时域小数频偏估计和补偿,消除小数部分引起的ICI(子载波间干扰),随后在频域上进行平移即可消除载波频偏的影响。
一、现有小数频偏估计方法
现有小数频偏估计方法是利用训练符号信息基于时域相关的方法。
如图1所示,它是相隔Nd个样值的两个时域重复训练符号示意图。假设发送的两个时域同步符号x1(n)和x2(n)之间有Nd个样值的时延,重复符号长度为L。
当存在Δfc
Figure A20081006678900041
时,可得接收到的两个训练符号可分别表示为
Figure A20081006678900042
Figure A20081006678900043
n=0,1,…,L-1
在接收端定义中间变量(两个训练符号的时域相关)
R t = Σ n = 0 L - 1 r 1 ( n ) r 2 * ( n )
Figure A20081006678900045
= exp ( j 2 πΔ f c N d T s ) Σ n = 0 L - 1 x 1 ( n ) x 2 * ( n ) + η
其中,
在不考虑噪声的情况下,并注意到x1(n)=x2(n),n=0,1,2,…,L-1,有
R t = exp ( j 2 πΔ f c N d T s ) Σ n = 0 L - 1 | x 1 ( n ) | 2
相关值Rt的相位为arg(Rt)=2πΔfcNdTs
可见这种情况下,所得相关值的相位与载波相位偏差
Figure A20081006678900052
没有关系,所以有载波频率偏差
Δ f c = f ^ c - f c = arg ( R t ) 2 π N d T s = arg ( R t ) NΔF 2 π N d
式中载波间隔ΔF=1/NTs
1)估计范围。对于相关值Rt的相位而言,arg(Rt)的变化范围是[-π,π),所以可估计载波频率范围是 | Δf c | ≤ N N d ΔF 2 .
2)估计的准确度。Jan等人已经证明,在AWGN信道中,这种估计是最大似然估计。
另外,同样是在AWGN信道中,Schimdl等人分析了该算法的性能,并给出了载波偏差Δfc估计值的方差 σ Δ f c 2 ∝ 1 LSNR .
现有估计方法的实施步骤如下:
步骤1:将接收到的信号按照对应顺序分别存入数组r1(n)和r2(n);
步骤2:使用 R t = Σ n = 0 L - 1 r 1 ( n ) r 2 * ( n ) 对r1(n)和r2 *(n)进行乘累加运算求得相关值Rt
步骤3:使用 Δf c = arg ( R t ) NΔF 2 π N d 得出小数频偏估计结果。
二、直流分量对现有小数频偏估计方法的影响
当存在直流分量α=αxy·j时,可得接收到的两个训练符号可分别表示为
Figure A20081006678900058
Figure A20081006678900059
n=0,1,…,L-1
在接收端定义中间变量(两个训练符号的时域相关)
R t = Σ n = 0 L - 1 r 1 ( n ) r 2 * ( n )
Figure A20081006678900062
= exp ( j 2 πΔ f c N d T s ) Σ n = 0 L - 1 x 1 ( n ) x 2 * ( n ) + L × | α | 2 + η
其中
Figure A20081006678900064
的期望为0。
相关值Rt的相位估计期望为 arg ( R t ) = arg ( exp ( j 2 πΔ f c N d T s ) Σ n = 0 L - 1 | x 1 ( n ) | 2 + L × | α | 2 )
由此可见,直流分量的引入不仅使得现有小数频偏方法从无偏估计退化为有偏估计,而且估计值的方差也明显增大。
图2是直流分量对现有小数频偏估计方法的影响示意图,它是小数频偏为1/4个子载波且存在直流分量时,CMMB手机电视仿真平台小数频偏估计精度的仿真结果。测试环境为AWGN信道,横坐标是信噪比,纵坐标为估计结果和真实值间的标准偏差。可见,当不存在直流分量时,标准偏差随着信噪比的提高单调快速下降,在SNR=-4dB达到0.01,SNR=10dB时接近千分之一。但在接收端引入和有用信号功率相等(0dB)的直流分量后,直接导致现有频偏估计方法崩溃,估计结果错误,且对信噪比变化不敏感。即使直流分量降低到相对于有用信号-10dB,估计偏差也非常大,且性能盆底劣于百分之一,达不到对OFDM高质量解调的要求。仅在直流分量下降到-20dB时,估计性能才能在大部分信噪比区间内满足要求。
发明内容
本发明的目的在于克服现有技术的不足,提供一种CMMB制式手机电视小数频偏估计方法,该方法大幅度降低了直流分量对估计精度的影响。
为解决上述技术问题,本发明是通过以下技术方案实现的:
一种CMMB制式手机电视小数频偏估计方法,包括如下步骤:
A:对接收到的基带数字信号进行预处理;
B:将经过预处理后的信号按照对应顺序分别存入数组r′1(n)和r′2(n);
C:使用 R t ′ = Σ n = M L - 1 r 1 ′ ( n ) r 2 ′ * ( n ) 对r′1(n)和r′2 *(n)进行乘累加运算求得相关值R′t
D:使用 Δf c = arg ( R t ′ ) NΔF 2 π N d 得出小数频偏估计结果。
所述步骤A中预处理为r(n)-r(n-M)形式的差分运算,L为重复符号长度,n=0,1,…,L-1,M为小于L-1的正整数。
所述步骤A中预处理为由FIR(Finite Impulse Response,有限长单位脉冲响应)形式的频域特性为高通性质的滤波器处理。
所述步骤A中预处理为由FIR形式的频域特性为带通性质的直流分量不在其通带内的滤波器处理。
本发明与现有方法相比具有以下有益效果:本发明增加了对进行小数频偏估计的数据进行预处理的步骤,预处理后再使用现有方法进行小数频偏估计,由于预处理的作用本发明的估计结果不受数据中所含直流分量的影响,较现有方法大幅提高了小数频偏估计的鲁棒性。
附图说明
图1是相隔Nd个样值的两个时域重复训练符号示意图;
图2是直流分量对现有小数频偏估计方法的影响示意图;
图3是存在直流分量情况下本发明和现有方法间的性能比较示意图;
图4是M=2时实现差分运算的硬件装置示意图;
图5是1-Z2的幅频响应示意图;
图6是用频域特性为高通性质的滤波器实现预处理的示意图;
图7是图6中滤波器的幅频响应示意图;
图8是用频域特性为带通性质的滤波器实现预处理的示意图;
图9是图8中滤波器的幅频响应示意图。
具体实施方式
为了更好地理解本发明,下面结合附图和具体实施例对本发明作进一步地描述。
本发明的CMMB制式手机电视小数频偏估计方法,增加了对进行小数频偏估计的数据进行预处理的步骤,预处理后再使用现有方法进行小数频偏估计,其步骤如下:
A:对接收到的基带数字信号进行预处理;
B:将经过预处理后的信号按照对应顺序分别存入数组r′1(n)和r′2(n);
C:使用 R t ′ = Σ n = M L - 1 r 1 ′ ( n ) r 2 ′ * ( n ) 对r′1(n)和r′2 *(n)进行乘累加运算求得相关值R′t
D:使用 Δf c = arg ( R t ′ ) NΔF 2 π N d 得出小数频偏估计结果。
其步骤A中预处理为r(n)-r(n-M)形式的差分运算,L为重复符号长度,n=0,1,…,L-1,M为小于L-1的正整数。
求取r′1(n)和r′2(n)为该差分运算预处理的过程,定义:
r 1 ′ ( n ) = r ( n ) - r ( n - M )
Figure A20081006678900084
Figure A20081006678900085
Figure A20081006678900086
r 2 ′ ( n ) = r ( n + N d ) - r ( n + N d - M )
Figure A20081006678900088
Figure A20081006678900089
Figure A200810066789000810
其中,n=0,1,…,L-1,M为小于L-1的正整数,两个时域同步符号之间有Nd个样值的时延,重复符号长度为L。
两个信号的时域相关为:
R t ′ = Σ n = M L - 1 r 1 ′ ( n ) r 2 ′ * ( n )
Figure A200810066789000812
= exp ( j 2 πΔ f c N d T s ) Σ n = M L - 1 [ x 2 ( n ) - x 2 ( n - M ) exp ( j 2 πΔ f c M T s ) ] * + η ′ ] [ [ x 1 ( n ) - x 1 ( n - M ) exp ( j 2 πΔ f c M T s ) ] ·
其中,
Figure A20081006678900091
在不考虑噪声的情况下,并注意到x1(n)=x2(n),n=0,1,2,…,L-1,有
R t ′ = exp ( j 2 πΔ f c N d T s ) Σ n = 0 L - 1 | x 1 ( n ) - x 1 ( n - M ) exp ( j 2 πΔ f c M T s ) | 2
相关值Rt的相位为arg(Rt)=2πΔfcNdTs
可见,本发明改进后的估计算法完全消除了直流分量对频偏估计的影响。请参阅图3,其仿真了M=2时针对CMMB手机电视系统现有方法和本发明方法的性能比较,仿真结果同样证实了改进算法对直流分量并不敏感,在不同程度的直流分量下都可以得到较为准确的估计结果。
上述差分运算预处理步骤(M=2时)可以采用如图4所示的硬件装置实现,该硬件装置使用了延迟器和加法器这两种器件。两级串联的延迟器401、402对输入数据进行延迟产生r(n-2)数据,加法器403配置为减运算方式,被减数输入端接输入数据r(n),减数输入端接延迟后的数据r(n-2),差也即输出端为差分运算结果,连接到现有小数频偏估计装置的输入端。此外,该差分运算预处理步骤同样可以用软件实现。
在实现复杂程度上,该差分运算预处理用r(n)=r(n)-r(n-M)取代了r(n),每个同步符号样点仅增加了一次减法运算,非常适合软硬件实现,且本发明对后续解调用到的接收信号未做任何变动,不会附加引入新的误差,由于后续模块均工作于频域,直流分量对后续模块的影响很小。
上述利用了差分结果依然满足延迟相关条件这一特性来消除直流分量的影响。除了差分解释外,本发明还可以理解为将同步符号通过1-ZM的滤波器进行预处理,由于该滤波器在Z=1有一个零点,因此经过预处理后的信号抵消了直流分量。此外,由于同步符号只利用了2048个子载波中频率较低的1536个子载波,高频段没有有用信号,在Z=-1处引入另一个零点,可以在一定程度上抑制带外噪声对估计性能的影响,却不用额外增加任何实现复杂度。如图5所示,它是1-Z2的幅频响应示意图,从图中可以看出该差分运算满足带通滤波器的特性,低、高频段无法通过,只有中间频段才能通过。
步骤A中预处理也可以为由FIR形式的频域特性为高通性质的滤波器进行处理。如图6所示,它是用频域特性为高通性质的滤波器实现预处理的示意图,该滤波器硬件装置使用了延迟器、乘法器和加法器这三种器件,滤波器为10阶,滤波器系数分别为
C0=C10=-0.036128128207098237;C1=C9=-0.060475430582165868;
C2=C8=-0.085696027762971985;C3=C7=-0.10775116816567994;
C4=C6=-0.12281831866385055;C5=0.87182999472165412。
图7为该滤波器的幅频响应示意图,低频段无法通过,高频段通过。预处理时,基带数字信号连接到滤波器的输入端,滤波器的输出端连接到现有小数频偏估计装置的输入端即可。
步骤A中预处理还可以为由FIR形式的频域特性为带通性质的滤波器进行处理,该带通性质的滤波器其直流分量不在其通带内。如图8所示,它是用频域特性为带通性质的滤波器实现预处理的示意图,该滤波器硬件装置使用了延迟器、乘法器和加法器这三种器件,滤波器为9阶,滤波器系数分别为C0=C4=-0.117765767330552;C1=C3=-0.21761846811752925;C2=0.73889497419310113。图9为该滤波器的幅频响应示意图,从图中可以看出低、高频段无法通过,只有中间频段才能通过。预处理时,基带数字信号连接到滤波器的输入端,输出端连接到现有小数频偏估计装置的输入端即可。
本方法对用于运算的数据先通过预处理步骤进行处理后,对其进行相关累加,使用相关累加值的相位求取小数频偏估计的结果,大幅度降低了直流分量对估计精度的影响,大幅提高了小数频偏估计的鲁棒性。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,应当指出,对于本领域的普通技术人员来说,凡是本发明的精神和原则之内所作的任何修改、等同替换或改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种CMMB制式手机电视小数频偏估计方法,其特征在于,包括如下步骤:
A:对接收到的基带数字信号进行预处理;
B:将经过预处理后的信号按照对应顺序分别存入数组r′1(n)和r′2(n);
C:使用 R t ′ = Σ n = M L - 1 r 1 ′ ( n ) r 2 ′ * ( n ) 对r′1(n)和r′2 *(n)进行乘累加运算求得相关值R′t
D:使用 Δ f c = arg ( R t ′ ) NΔF 2 π N d 得出小数频偏估计结果。
2.根据权利要求1所述的CMMB制式手机电视小数频偏估计方法,其特征在于,所述步骤A中预处理为r(n)-r(n-M)形式的差分运算,L为重复符号长度,n=0,1,…,L-1,M为小于L-1的正整数。
3.根据权利要求1所述的CMMB制式手机电视小数频偏估计方法,其特征在于,所述步骤A中预处理为由FIR形式的频域特性为高通性质的滤波器处理。
4.根据权利要求1所述的CMMB制式手机电视小数频偏估计方法,其特征在于,所述步骤A中预处理为由FIR形式的频域特性为带通性质的直流分量不在其通带内的滤波器处理。
CNA2008100667895A 2008-04-25 2008-04-25 一种cmmb制式手机电视小数频偏估计方法 Pending CN101567868A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNA2008100667895A CN101567868A (zh) 2008-04-25 2008-04-25 一种cmmb制式手机电视小数频偏估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNA2008100667895A CN101567868A (zh) 2008-04-25 2008-04-25 一种cmmb制式手机电视小数频偏估计方法

Publications (1)

Publication Number Publication Date
CN101567868A true CN101567868A (zh) 2009-10-28

Family

ID=41283821

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2008100667895A Pending CN101567868A (zh) 2008-04-25 2008-04-25 一种cmmb制式手机电视小数频偏估计方法

Country Status (1)

Country Link
CN (1) CN101567868A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103401537A (zh) * 2013-07-31 2013-11-20 中科院微电子研究所昆山分所 一种判决电路及接收机
CN104735011A (zh) * 2013-12-24 2015-06-24 卓胜微电子(上海)有限公司 用于复杂信道环境下鲁棒估计采样频率偏差的装置及方法
CN110932770A (zh) * 2019-11-13 2020-03-27 北京邮电大学 一种低轨卫星星地链路同步序列设计及频偏估计方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103401537A (zh) * 2013-07-31 2013-11-20 中科院微电子研究所昆山分所 一种判决电路及接收机
CN103401537B (zh) * 2013-07-31 2016-06-22 中科院微电子研究所昆山分所 一种判决电路及接收机
CN104735011A (zh) * 2013-12-24 2015-06-24 卓胜微电子(上海)有限公司 用于复杂信道环境下鲁棒估计采样频率偏差的装置及方法
CN104735011B (zh) * 2013-12-24 2019-01-29 江苏卓胜微电子股份有限公司 用于复杂信道环境下鲁棒估计采样频率偏差的装置及方法
CN110932770A (zh) * 2019-11-13 2020-03-27 北京邮电大学 一种低轨卫星星地链路同步序列设计及频偏估计方法
CN110932770B (zh) * 2019-11-13 2020-10-16 北京邮电大学 一种低轨卫星星地链路同步序列设计及频偏估计方法

Similar Documents

Publication Publication Date Title
Yeh et al. Channel estimation using pilot tones in OFDM systems
US7535820B2 (en) OFDM signal receiving apparatus and OFDM signal receiving method
EP2362596B1 (en) Method and device for filtering orthogonal frequency division multiplexing channel estimate result
US7693039B2 (en) Apparatus and method for carrier frequency synchronization in an OFDM system
USRE44622E1 (en) Coarse carrier frequency offset estimation for a receiver
EP1460814B1 (en) Coarse frequency synchronization for multicarrier receivers
US20100166050A1 (en) Time error estimation for data symbols
JP2004032748A (ja) 直交周波数分割多重化通信システムのためのマルチパス等化技術
US20080118006A1 (en) Symbol synchronization using phase discontinuity for dvb-t systems in awgn channels
EP0838928B1 (en) Equalisation of multicarrier signals
US20060158366A1 (en) Receiver and method thereof
CN101068232B (zh) 得到信道时域响应方法及装置、ofdm符号精同步方法及装置
US8238481B2 (en) Blind channel estimation for PSK and D-PSK modulated multicarrier communication systems
US20110167321A1 (en) Apparatus and method of transmitting and receiving data for improving transmission rate
CN109194594A (zh) 一种基于连续载波聚合的相位噪声抑制方法
US20100074346A1 (en) Channel estimation in ofdm receivers
CN101902426B (zh) 一种实现小数频偏估计的方法及装置
CN101888363B (zh) 一种ofdm接收机中的信号解调方法及ofdm接收机
US20090097596A1 (en) Methods and Systems for Impulse Noise Compensation for OFDM Systems
CN101567868A (zh) 一种cmmb制式手机电视小数频偏估计方法
CN100550998C (zh) 多阶结构的载波恢复装置
Yang et al. Channel estimation of DMB-T
EP1331783A2 (en) Apparatus and method for recovering symbol timing in OFDM receiver
CN107026806B (zh) 接收设备、接收方法和程序
US9369329B1 (en) Low-complexity non-data-aided estimation of symbol time offset in OFDM systems

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20091028