CN104022981B - 一种正交幅度调制信号的盲载波频偏估计方法 - Google Patents
一种正交幅度调制信号的盲载波频偏估计方法 Download PDFInfo
- Publication number
- CN104022981B CN104022981B CN201410298701.8A CN201410298701A CN104022981B CN 104022981 B CN104022981 B CN 104022981B CN 201410298701 A CN201410298701 A CN 201410298701A CN 104022981 B CN104022981 B CN 104022981B
- Authority
- CN
- China
- Prior art keywords
- value
- estimation
- algorithms
- wph
- carrier frequency
- 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
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明公开了一种正交幅度调制信号的盲载波频偏估计方法,属于移动通信技术领域。在该方法中,QAM信号解调***接收的射频信号经过匹配滤波、模数转换以及定时恢复等处理后得到含有载波频偏f0的样值rn,对于rn中含有的频偏成分,不依赖其他任何的数据信息,联合FFT算法和WPH算法,通过多级估计逐步校正估计误差得到较为精确的频偏估计值。该方法结合了FFT算法和WPH算法,规避了这两种算法中存在的不足,结合了各自的优势,通过多级估计的方式计算频偏,并对信号进行补偿,为下一阶段进行的符号判决提供了必要的先决条件。该方法在一定程度上缓解了算法中存在的计算精度和计算复杂度的有效性和可靠性之间的矛盾。
Description
技术领域
本发明属于移动通信技术领域,涉及一种正交幅度调制信号的盲载波频偏估计方法,特别是一种对数字调相信号进行非数据辅助的多普勒频偏估计的方法。
背景技术
在数字通信中,由于接收机振荡器不稳定和多普勒效应等因素,使得接收信号的载波与本地载波之间不能完全同步,存在一定的频率偏差,因此在经过相干解调和下变频后得到的基带信号中,通常含有载波频率偏移。在突发模式传输通信***的相干检测中频偏估计具有关键作用,估计的准确性直接影响接收信号的误码率。
假设数据辅助的前导信息已知或者星座图类型和增益已知,那么在高信噪比条件下,可以较容易的估算接近理想值的频率偏移值。但是,额外的辅助数据需要占用频谱资源,降低了***传输的有效性。
非数据辅助同步的频偏估计方面研究主要包括快似然函数估计算法、快速傅里叶变换(FFT)算法、直方图算法等。其中文献(Blind carrier frequency offset estimationfor QAM signals based on weighted 4th power of signal samples)中的似然估计算法的均方误差接近理想值,然而此算法没有闭式解,因此计算量较大,难以进行实时处理。而文献(Non-data-aided wide-range frequency offset estimation for QAM opticalcoherent receivers)中的FFT算法由于需要较长数据长度才能使频谱更接近连续模型,因此当有较短长度数据时,估计精度不高。相位直方图估计法如文献(Frequency offsetestimation for unknown QAM constellations),采用了粗-细估计,估计精度不够好,而且运算复杂度较大。因此有必要找到一种既能保证频率估计有较高的估计精度,又有较好运算效率的非数据辅助估计算法。
发明内容
有鉴于此,本发明的目的在于提供一种正交幅度调制信号的盲载波频偏估计方法,该方法利用WPH算法和FFT算法联合估计载波频偏,利用多级迭代保证估计精度,并对估计过程进行简化,减少运算量。
为达到上述目的,本发明提供如下技术方案:
一种正交幅度调制信号的盲载波频偏估计方法,在该方法中,QAM信号解调***接收的射频信号经过匹配滤波、模数转换以及定时恢复等处理后得到含有载波频偏f0的样值rn,对于rn中含有的频偏成分,不依赖其他任何的数据信息,联合FFT算法和WPH算法,通过多级估计逐步校正估计误差得到较为精确的频偏估计值。
进一步,所述多级估计的第一级采用FFT算法模型,缩小频偏估计范围,以下各级利用WPH算法模型迭代运算。
进一步,所述FFT算法,包括利用对四次方后的接收信号样值进行FFT变换,经过测量最大谱线的信息得到载波频偏值。
进一步,所述WPH算法包括,样值rn经过fe频率补偿后yn的相位设为fe为一个变量,取值范围为[-1/8,1/8],n=(1,2,···,N)表示第n个数据样值,N表示总数据长度;通过统计在[0,2π]的K等分相位间隔ψk中的特性得到yn的相位直方图,即K为常数,可由设计者根据具体情况自行设定,一般来说,高阶QAM时K值设定较大,低阶QAM时K值设定较小;k=(0,1,···,K-1);rect(·)表示矩形运算,ψk=π(2k+1)/K;进而经过测量该相位直方图加权后的函数的二阶矩的峰值,得到载波频率偏移值。
进一步,所述的计算相位直方图的方法如下:直方图ξk的计算从逐步统计相位的位置解析为由公式直接判决确定其位置,其中arg(Y)表示求yn的角度;表示取整,如果,,其余ξ值为0;否则ξk=0。
本发明的有益效果在于:本发明所述方法结合了FFT算法和WPH算法,规避了这两种算法中存在的不足,结合了各自的优势,通过多级估计的方式计算频偏,并对信号进行补偿,为下一阶段进行的符号判决提供了必要的先决条件。该方法在一定程度上缓解了算法中存在的计算精度和计算复杂度的有效性和可靠性之间的矛盾。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为接收端QAM信号解调模型;
图2为WPH算法流程图;
图3为FFT算法流程图;
图4为QAM信号盲载波频偏估计的WPH和FFT联合算法流程图。
具体实施方式
本发明提供的盲载波频偏估计方法,采用以下方式实现:
对于QAM信号解调***,接收的射频信号经过匹配滤波、模数转换以及定时恢复等处理后得到含有载波频偏f0的rn,WPH算法利用fe对rn频率进行四倍补偿,其中fe的取值范围为[-1/8,1/8],即
arg(·)表示求角度。我们将[0,2π]相位范围K等分,设为ψk,判断此时yn的相位相对于ψk的矩形运算,得到式(2)所示的加权相位直方图
N为rn的样值个数,rect(·)代表矩形运算。运算需要依赖逐步统计判断相位的位置,复杂度较高,可以进行如下简化:采用直接计算确定其位置,其中arg(Y)表示求yn的角度;代表取整。当时,,其余ξ值为0;否则ξk=0。此时测量直方图a二阶矩的峰值得到载波频偏值,如式(3)所示
本发明的载波频偏估计采用多级迭代估计算法,然而第一级估计时估计范围为[-1/8,1/8],当设定fe取值时,就应为fe=-1/8+i·△f,其中△f为步进值,i=(0,1,2,···,I-1),I=(1/8-(-1/8))/△f,由分析可知,这样的估计范围使得fe取值较多,相应的增加了计算复杂度,因此第一步采用FFT算法,确定一个粗略的频偏值,具体方法是
利用式(4)先求得rn 4的频谱,此时峰值谱线对应的频率信息,即为所求频偏,但是由于此时的FFT算法采用有限长度的样值数据进行计算,估计值偏差较大,因此,要联合WPH算法进行迭代计算以纠正误差。
下面将结合附图,对本发明的优选实施例进行详细的描述。
QAM信号接收端解调如图1所示,S(t)为发送信号,N(t)为噪声,接收端经过相干解调和下变频后得到基带信号,由于接收机振荡器不稳定和多普勒效应等因素,使得接收信号的载波与本地载波之间不能完全同步,存在一定的频率偏差,即载波频率偏移。因此要经过算法估计频偏值然后进行,以达到正确解调的目的。
在本发明的算法中,包含了两种估计方法,即WPH算法和FFT算法。图2为WPH算法的流程图。具体步骤包括如下:
步骤21,对接收数据进行四倍频率补偿,补偿值为fe,因此得到补偿后的样值yn;
步骤22,利用判断此时yn的相位在等间隔相位槽ψk中对应的位置,然后对yn求取加权相位直方图a;
步骤23,求a的二阶矩,定义为p;
步骤24,根据二阶矩p测量其峰值,得到频偏估计值;
相应的,图3为FFT计算流程,具体过程如下:
步骤31,对接收样值求四次方;
步骤32,对四次方后的样值进行FFT变换;
步骤33,测量变换后的离散谱谱峰值;
步骤34,根据峰值信息计算得到频偏估计值,但是此值误差较大,为粗估计值。
由于本发明采用多级估计,并且第一步粗估计采用FFT算法,然后根据WPH算法进行迭代计算以矫正误差。总体计算流程可以归纳为图4所示的流程图。
步骤41,根据接收样值yn采用FFT算法,计算粗频偏估计值σ1;
步骤42,设定此次搜索步进值△f;
步骤43,以前一步计算的频偏估计值σ为中心确定fe的检索范围[σ_left,σ_right]
步骤44,依据步骤42、步骤43确定的参数值,得到fe=σ_left+i·△f,i=(0,1,2,···,I-1),其中I=(σ_right-σ_left)/△f,代入WPH算法中进行计算得到估计值σ=σ_left+imax·Δf;
步骤45,判定迭代次数是否达到预设值。否,则返回步骤42,重新设定更精细的步进值,并确定此时检索区域,进行校正估计;是,则输出估算值即为。
最后说明的是,以上优选实施例仅用以说明本发明的技术方案而非限制,尽管通过上述优选实施例已经对本发明进行了详细的描述,但本领域技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离本发明权利要求书所限定的范围。
Claims (2)
1.一种正交幅度调制信号的盲载波频偏估计方法,其特征在于:QAM信号解调***接收的射频信号经过匹配滤波、模数转换以及定时恢复处理后得到含有载波频偏f0的样值rn,对于rn中含有的频偏成分,不依赖其他任何的数据信息,联合FFT算法和WPH算法,通过多级估计逐步校正估计误差得到较为精确的频偏估计值;
所述多级估计的第一级采用FFT算法模型,缩小频偏估计范围,以下各级利用WPH算法模型迭代运算;
所述FFT算法,包括利用对四次方后的接收信号样值进行FFT变换,经过测量最大谱线的信息得到载波频偏值;
所述WPH算法包括,样值rn经过fe频率补偿后的样值yn的相位设为fe为一个变量,取值范围为[-1/8,1/8],n=(1,2,...,N)表示第n个数据样值,N表示总数据长度;通过统计在[0,2π]的K等分相位间隔ψk中的特性得到补偿后的样值yn的相位直方图,即K为常数,由设计者根据具体情况自行设定,高阶QAM时K值设定较大,低阶QAM时K值设定较小;k=(0,1,...,K-1);rect(·)表示矩形运算,ψk=π(2k+1)/K;进而经过测量该相位直方图加权后的函数的二阶矩的峰值,得到载波频率偏移值。
2.根据权利要求1所述的一种正交幅度调制信号的盲载波频偏估计方法,其特征在于:计算所述的相位直方图的方法如下:直方图ξk的计算从逐步统计相位的位置解析为由公式直接判决确定其位置,其中arg(Y)表示求yn的角度;表示取整,如果ki∈[0,K-1],其余ξ值为0;否则ξk=0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298701.8A CN104022981B (zh) | 2014-06-27 | 2014-06-27 | 一种正交幅度调制信号的盲载波频偏估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298701.8A CN104022981B (zh) | 2014-06-27 | 2014-06-27 | 一种正交幅度调制信号的盲载波频偏估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104022981A CN104022981A (zh) | 2014-09-03 |
CN104022981B true CN104022981B (zh) | 2018-02-16 |
Family
ID=51439554
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410298701.8A Active CN104022981B (zh) | 2014-06-27 | 2014-06-27 | 一种正交幅度调制信号的盲载波频偏估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104022981B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104486288B (zh) * | 2014-12-01 | 2017-09-15 | 北京理工大学 | 一种适用于pcm/fm遥测接收机的载波频偏抑制方法 |
CN105790848B (zh) * | 2014-12-23 | 2018-06-08 | 深圳市中兴微电子技术有限公司 | 一种载波相位估计方法和装置 |
CN106059984A (zh) * | 2015-12-10 | 2016-10-26 | 国网山东省电力公司烟台供电公司 | 一种数字调相信号载波相偏估计方法 |
CN106992950A (zh) * | 2016-01-20 | 2017-07-28 | 晨星半导体股份有限公司 | 载波频偏估测装置与载波频偏估测方法 |
TWI623233B (zh) * | 2016-03-30 | 2018-05-01 | 晨星半導體股份有限公司 | 訊號檢測方法及訊號檢測裝置 |
CN105933265B (zh) * | 2016-03-31 | 2019-04-23 | 泉州装备制造研究所 | 一种对qam信号进行非数据辅助的相位噪声盲估计方法 |
CN106027454A (zh) * | 2016-04-02 | 2016-10-12 | 上海大学 | 基于频率补偿的qam信号的频偏估计方法 |
CN110247874A (zh) * | 2018-03-09 | 2019-09-17 | 中国电力科学研究院有限公司 | 基于星座图的电力线宽带载波性能分析方法及装置 |
CN109005139A (zh) * | 2018-06-30 | 2018-12-14 | 华南理工大学 | 一种正交幅度调制解调方法及其星座图优化方法 |
CN109088839A (zh) * | 2018-09-10 | 2018-12-25 | 金陵科技学院 | 一种适用于π/4-DQPSK的频率捕获方法 |
CN116346215A (zh) * | 2021-03-31 | 2023-06-27 | 北京百度网讯科技有限公司 | 信号处理方法及信号处理*** |
CN114448760B (zh) * | 2022-01-14 | 2023-12-19 | 山东大学 | 一种基于人工蜂群和网格聚类的数字非恒包络调制信号频偏盲估计方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101222473A (zh) * | 2008-02-01 | 2008-07-16 | 北京天碁科技有限公司 | 一种频偏估计方法及装置 |
CN102195905A (zh) * | 2011-03-31 | 2011-09-21 | 重庆重邮信科通信技术有限公司 | 一种频偏估计方法及装置 |
CN102497339A (zh) * | 2011-12-12 | 2012-06-13 | 安徽省菲特科技股份有限公司 | 一种qam通信***中的频偏估计方法 |
US8446974B2 (en) * | 2010-01-27 | 2013-05-21 | National Instruments Corporation | Blind mechanism for the joint estimation of frequency offset and phase offset for QAM modulated signals |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8908809B2 (en) * | 2010-12-15 | 2014-12-09 | At&T Intellectual Property I, L.P. | Complexity reduced feed forward carrier recovery methods for M-QAM modulation formats |
-
2014
- 2014-06-27 CN CN201410298701.8A patent/CN104022981B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101222473A (zh) * | 2008-02-01 | 2008-07-16 | 北京天碁科技有限公司 | 一种频偏估计方法及装置 |
US8446974B2 (en) * | 2010-01-27 | 2013-05-21 | National Instruments Corporation | Blind mechanism for the joint estimation of frequency offset and phase offset for QAM modulated signals |
CN102195905A (zh) * | 2011-03-31 | 2011-09-21 | 重庆重邮信科通信技术有限公司 | 一种频偏估计方法及装置 |
CN102497339A (zh) * | 2011-12-12 | 2012-06-13 | 安徽省菲特科技股份有限公司 | 一种qam通信***中的频偏估计方法 |
Non-Patent Citations (2)
Title |
---|
Frequency Offset Estimation for Unknown QAM Constellations;Stefania Colonnese等;《IEEE Transactions on Communications》;20120210;正文第Ⅱ-Ⅲ部分 * |
Non-Data-Aided Wide-Range Frequency Offset Estimator for QAM Optical Coherent Receivers;Tadao Nakagawa等;《Optical Fiber Communication Conference and Exposition (OFC/NFOEC), 2011 and National Fiber Optic Engineers Conference》;20110613;正文第2部分 * |
Also Published As
Publication number | Publication date |
---|---|
CN104022981A (zh) | 2014-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104022981B (zh) | 一种正交幅度调制信号的盲载波频偏估计方法 | |
CN109495410B (zh) | 高动态pcm/fm信号载波频率精确估计方法 | |
CN105245303B (zh) | 一种高速突发解调同步*** | |
US20030220084A1 (en) | Frequency offset estimator | |
CN111683034B (zh) | 一种基于ofdm的大多普勒无线通信时频同步方法 | |
CN102546500A (zh) | 基于导频和软信息联合辅助的soqpsk载波同步方法 | |
CN108512791A (zh) | 基于定时频偏补偿的星载ais解调方法 | |
CN108055058A (zh) | 载波多普勒及其变化率的高精度测量方法 | |
EP2984802A2 (en) | Digital radio transmissions | |
CN106059984A (zh) | 一种数字调相信号载波相偏估计方法 | |
CN109561035A (zh) | 一种基于太赫兹的超宽带单通道信号盲检测方法 | |
CN113037663B (zh) | 一种改进的适用于非恒包络信号的码元速率估计方法 | |
CN110691051A (zh) | 一种基于fft的gmsk信号频偏估计算法 | |
CN104363194A (zh) | 基于波形变换的psk调制识别方法 | |
CN108347396B (zh) | 追踪并补偿调制信号的定时偏移和频率偏移的***和方法 | |
CN113472483A (zh) | 数字调制信号码元速率及码元转换时刻盲估计方法 | |
US20120140858A1 (en) | Fine symbol timing estimation | |
CN114531329A (zh) | 一种多路msk信号载波频率估计方法、***及应用 | |
WO2021255124A1 (en) | Demodulating modulated signals | |
CN105933265B (zh) | 一种对qam信号进行非数据辅助的相位噪声盲估计方法 | |
CN105721075B (zh) | 数据辅助同频复用卫星干扰信号强度测量方法及装置 | |
CN112039591A (zh) | 一种基于二分法的载波相位估计算法 | |
CN105262505B (zh) | 跳频通信干扰条件下的信噪比估计方法 | |
CN112671684B (zh) | 一种短时突发bpsk信号的自适应解调方法 | |
CN108353066B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |