CN107623647B - 一种基于离散导频辅助的载波同步方法 - Google Patents
一种基于离散导频辅助的载波同步方法 Download PDFInfo
- Publication number
- CN107623647B CN107623647B CN201710694666.5A CN201710694666A CN107623647B CN 107623647 B CN107623647 B CN 107623647B CN 201710694666 A CN201710694666 A CN 201710694666A CN 107623647 B CN107623647 B CN 107623647B
- Authority
- CN
- China
- Prior art keywords
- frequency offset
- carrier
- pilot
- symbol
- change rate
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000009825 accumulation Methods 0.000 claims abstract description 32
- 230000001427 coherent effect Effects 0.000 claims abstract description 18
- 230000003595 spectral effect Effects 0.000 claims description 50
- 238000005311 autocorrelation function Methods 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000004891 communication Methods 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
Images
Landscapes
- Synchronisation In Digital Transmission Systems (AREA)
Abstract
本发明涉及一种基于离散导频辅助的载波同步方法。该方法要求输入信号具有非连续、等间隔分布的离散导频点的结构特点,方法实施步骤为:(1)从采集的数据符号中抽取导频符号;(2)根据抽取的导频符号获取载波频偏变化率,并根据所述载波频偏变化率对包含有所述导频符号的数据符号进行频偏变化率补偿;(3)根据频偏变化率补偿后的导频符号获取载波频偏,并根据所述载波频偏对包含有所述导频符号的数据符号进行频偏补偿;(4)将所述频偏补偿后的导频符号进行相干累加,得到载波初始相位;(5)通过载波初始相位对进行频偏补偿后的各数据符号进行相位校正,完成载波同步。
Description
技术领域
本发明涉及猝发通信技术领域,特别涉及一种基于离散导频辅助的载波同步方法。
背景技术
在TDMA网络***中,数据帧的长短决定了组网节点数目上限,进而决定着通信网络的容量。为保证网络具有较大容量并且兼顾数据帧的传输效率,帧长设计通常为10ms级别,并且在一些特殊应用领域如深远海通信,出于保密的需要,通信方式多选择猝发通信体制。由于数据帧持续时间较短,故在信噪比较低,多普勒动态较大场景下,传统的闭环跟踪方式受限于环路入锁时间,并不适用,因此基于导频辅助的开环跟踪方式是猝发通信的重要应用方向。
目前,在进行导频辅助的短时数据帧同步算法时,导频符号非连续、均匀分布在整个信号帧内,基于此通信格式,首先估计出导频符号的相位值,然后利用插值的方法,插出相邻导频符号之间的所有数据符号值,然而如果导频符号相位估计错误,则数据符号的相位也都校正错误,造成载波同步误差较大的问题。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于离散导频辅助的载波同步方法,利用导频符号进行频偏变化率估计,并对数据符号进行频偏变化率剥离,然后对频偏变化率剥离后的数据符号进行载波频偏估计,并对数据符号进行载波频偏剥离,最后进行初始相位校正,完成载波同步。可以解决现有载波同步误差较大,无法适用于大频率动态应用场景下的载波同步。
本发明的上述目的是通过如下技术方案予以实现:
一种基于离散导频辅助的载波同步方法,包括:
(1)从采集的数据符号中抽取导频符号;
(2)根据抽取的导频符号获取载波频偏变化率,并根据所述载波频偏变化率对包含有所述导频符号的数据符号进行频偏变化率补偿;
(3)根据频偏变化率补偿后的导频符号获取载波频偏,并根据所述载波频偏对包含有所述导频符号的数据符号进行频偏补偿;
(4)将所述频偏补偿后的导频符号进行相干累加,得到载波初始相位;
(5)通过载波初始相位对进行频偏补偿后的各数据符号进行相位校正,完成载波同步。
进一步地,所述数据符号的长度为N,其中,N=K×(P+L),K表示将长度为N的数据符号均匀化分为K段,L为每一段数据符号中包含的数据符号个数,P为每一段数据符号中导频符号的个数。
进一步地,所述步骤(2)中获取载波频偏变化率的方法为:
将每一段数据符号中的P个导频符号进行相干累加,得到K个累加结果p1,p2,…pK;
进一步地,所述步骤(3)中获取载波频偏的方法为:
将进行频偏变化率补偿后的每一段数据符号中的P个导频符号进行相干累加,得到K个累加结果p1,p2,…pK;
对累加结果p1,p2,…pK进行FFT运算,得到所述自相关函数序列对应的FFT谱线;
进一步地,所述步骤(4)中获取载波初始相位的方法为:
进一步地,所述步骤(5)中通过载波初始相位对进行频偏补偿后的各数据符号进行相位校正,完成载波同步的方法为:
本发明与现有技术相比具有如下有益效果:
(1)本发明提供的基于离散导频辅助的载波同步方法,在估计多普勒频率变化率和多普勒频偏时,均采用Jacobsen估计算法而非传统的FFT估计算法,即获取到FFT谱线后,对峰值谱线频率进行修正,并采用修正后的峰值谱线频率获取载波频偏变化率和载波频偏,从而可以保证该同步方法在低信噪比条件下,通过极小的FFT运算规模获得比传统FFT运算方案更高的估计精度,从而获得更理想的载波同步效果。
(2)本发明在频偏变化率估计、频偏估计以及最终的载波初始相位校正上,均采用相干处理方式,利用多个导频符号的相干累加提升判决信噪比,降低了错误判决概率,从而保证本发明提供的方法可以工作于低信噪比条件下。
(3)本发明提供的基于离散导频辅助的载波同步方法,首先利用导频符号进行多普勒频偏变化率估计,在多普勒频偏变化率补偿之后,再进行多普勒频偏估计。从而可以保证本发明提供的方法可以应用在更大的多普勒动态条件下,进而可以适用于大频率动态应用场景下的载波同步。
附图说明
图1为本发明的离散导频辅助的信号格式;
图2为本发明的信号处理流程图。
具体实施方式
下面结合附图2本发明的实现过程进行进一步详细阐述:
(1)从采集的数据符号中抽取导频符号。
对于本发明,要求输入的数据符号具有非连续、等间隔导频符号分布的特点,即每隔L个数据符号,有P个导频符号,信号结构可以如图1所示。则采集的数据符号符号长度N=K×(P+L),其中K表示将长度为N的数据符号均匀化分为K段,K的取值应满足累计导频符号信噪比的要求,通常需要满足其中,表示用dB表示的单个导频符号信噪比。并且预置输入信号已经完成粗同步,残余的多普勒频偏在导频图案频率估计的范围内。
(2)根据抽取的导频符号获取载波频偏变化率,并根据所述载波频偏变化率对包含有所述导频符号的数据符号进行频偏变化率补偿。
本发明中的频偏变化率估计方法采用延迟自相关FFT算法,算法原理描述如下:
设基带接收信号为:
其中,fd为载波频偏,f′d为载波频偏变化率,则信号自相关函数为:
其中,(·)*表示(·)共轭函数。由此可见自相关函数实际上为频率为f′dτ的单载波信号,因此频偏变化率估计的问题实际为频偏估计的问题。
本发明的频偏估计算法采用Jacobsen估计算法,具体为:
a、将每一段数据符号中的P个导频符号进行相干累加,得到K个累加结果p1,p2,…pK。
对于本发明,此处进行频偏变化率估计,是为了解决大多普勒频偏变化率条件下,频偏无法直接估计的问题,同时,此处采用Jacobsen估计算法进行频率估计,进而计算出频偏变化率,是因为Jacobsen估计算法较传统的FFT算法相比,具有更小的FFT规模,并且在相同信噪比条件下,获得更高的估计精度。
(3)根据频偏变化率补偿后的导频符号获取载波频偏,并根据所述载波频偏对包含有所述导频符号的数据符号进行频偏补偿。
本发明中载波频偏算法采用Jacobsen估计算法,具体步骤为:
a、将进行频偏变化率补偿后的每一段数据符号中的P个导频符号进行相干累加,得到K个累加结果p1,p2,…pK。
b、对累加结果p1,p2,…pK进行FFT运算,得到所述自相关函数序列对应的FFT谱线。
(4)将所述频偏补偿后的导频符号进行相干累加,得到载波初始相位。
具体地,根据公式进行计算,其中,为初始相位,sp,0为导频符号的初始值,通常sp,0=1+j,S'p为频偏补偿后的导频符号进行相干累加后的结果。其中,s′p,j第j个导频符号,PK为存储的N个接收信号中,导频符号的总个数。
(5)通过载波初始相位对进行频偏补偿后的各数据符号进行相位校正,完成载波同步。
实施例
假设接收信号具有图1所示结构,其中导频符号初始值为s0=1+j,业务数据符号采用QPSK调制方式,则接收信号可用下式描述:
其中,fd为捕获之后残余的多普勒频偏,f′d为多普勒频偏变化率,Ts为符号周期,sn=In+jQn为发送的第n个数据符号,n=1,2,…,N。特别地,对于导频符号,sn=1+j。
按照步骤(1)提取数据符号中的所有导频符号,得到序列集合{s′1,1,s'2,1,…,s'P,1,s1',2,s'2,2,…,s'P,K}。
按照步骤(2)利用抽取的导频符号,通过Jacobsen估计算法计算出载波频偏变化率:
a、将抽取出来的导频符号按照P个为一组进行相干累加,得到K个累加结果p1,p2,…pK,其中m=1,2,…K。一般情况下,残余多普勒在连续P个导频符号内累计相位忽略不计,则其中,T′s=(P+L)Ts表示导频累加后样点序列p1,p2,…pK的周期。为方便讨论,省略T′s,于是m=1,2,…K。
d、找出FFT的峰值谱线k其前后相邻两根谱线k-1,k+1所对应的FFT变换值Xk,Xk-1,Xk+1,然后按照下式得到峰值位置的修正项δ,δ可以为正数或者负数。
其中,Re[·]表示对[·]取实部。
e、修正后的峰值谱线kpeak及其对应的估计频偏分别为下式所示:
kpeak=k+δ
此时,对式1的数据符号进行频偏变化率补偿,补偿后的接收符号序列为
其中,θ0为补偿后的初始相位。
按照步骤(3),利用经频偏变化率补偿后的导频符号,通过Jacobsen估计算法计算出载波频偏,具体地:
C、获取FFT的峰值谱线k其前后相邻两根谱线k-1,k+1所对应的FFT变换值Xk,Xk-1,Xk+1,然后按照下式得到峰值位置的修正项δ,δ可以为正数或者负数。
其中,Re[·]表示对[·]取实部。
D、修正后的峰值谱线kpeak及其对应的估计频偏分别为下式所示。
kpeak=k+δ
此时,对式2的数据符号进行频偏补偿,补偿后的接收符号序列为
按照步骤(5),对式3数据符号进行初始相位校正,补偿后的接收符号序列为
s””(n)=sn (式4)
s””(n)进行解调、译码等操作。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (6)
1.一种基于离散导频辅助的载波同步方法,其特征在于,包括:
(1)从采集的数据符号中抽取导频符号;
(2)根据抽取的导频符号获取载波频偏变化率,并根据所述载波频偏变化率对包含有所述导频符号的数据符号进行载波频偏变化率补偿;
(3)根据载波频偏变化率补偿后的导频符号获取载波频偏,并根据所述载波频偏对包含有所述导频符号的数据符号进行频偏补偿;
(4)将所述频偏补偿后的导频符号进行相干累加,得到载波初始相位;
(5)通过载波初始相位对进行频偏补偿后的各数据符号进行相位校正,完成载波同步。
2.根据权利要求1所述的基于离散导频辅助的载波同步方法,其特征在于,所述数据符号的长度为N,其中,N=K×(P+L),K表示将长度为N的数据符号均匀划分为K段;L为每一段数据符号中包含的数据符号个数,P为每一段数据符号中导频符号的个数。
3.根据权利要求2所述的一种基于离散导频辅助的载波同步方法,其特征在于,所述步骤(2)中获取载波频偏变化率的方法为:
将每一段数据符号中的P个导频符号进行相干累加,得到K个累加结果p1,p2,…pK;
4.根据权利要求2所述的一种基于离散导频辅助的载波同步方法,其特征在于,所述步骤(3)中获取载波频偏的方法为:
将进行载波频偏变化率补偿后的每一段数据符号中的P个导频符号进行相干累加,得到K个累加结果p1',p2',…pK';
对累加结果p1',p2',…pK'进行FFT运算,得到自相关函数序列对应的FFT谱线;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710694666.5A CN107623647B (zh) | 2017-08-15 | 2017-08-15 | 一种基于离散导频辅助的载波同步方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710694666.5A CN107623647B (zh) | 2017-08-15 | 2017-08-15 | 一种基于离散导频辅助的载波同步方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107623647A CN107623647A (zh) | 2018-01-23 |
CN107623647B true CN107623647B (zh) | 2020-07-14 |
Family
ID=61088047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710694666.5A Active CN107623647B (zh) | 2017-08-15 | 2017-08-15 | 一种基于离散导频辅助的载波同步方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107623647B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109391572A (zh) * | 2018-11-22 | 2019-02-26 | 中国电子科技集团公司第五十四研究所 | 一种基于相位增量的载波频偏估计方法 |
CN109600331B (zh) * | 2019-01-30 | 2020-09-15 | 北京智芯微电子科技有限公司 | 无线通信网络中的频偏同步方法 |
CN110505172B (zh) * | 2019-07-12 | 2022-04-26 | 四川安迪科技实业有限公司 | 一种基于离散导频辅助的频偏估计方法及装置 |
CN110768919B (zh) * | 2019-09-11 | 2022-09-16 | 武汉船舶通信研究所(中国船舶重工集团公司第七二二研究所) | 用于卫星通信***的载波同步方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102624419A (zh) * | 2012-04-23 | 2012-08-01 | 西安电子科技大学 | 突发直接序列扩频***的载波同步方法 |
CN103401661A (zh) * | 2013-08-08 | 2013-11-20 | 电子科技大学 | 一种基于mimo雷达通信的一体化编解码方法 |
CN104852875A (zh) * | 2015-03-26 | 2015-08-19 | 中国人民解放军理工大学通信工程学院卫星通信军队重点实验室 | 高动态大频偏突发信号的频偏估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI412257B (zh) * | 2009-03-20 | 2013-10-11 | Univ Nat Central | Combined Carrier Synchronization and Channel Equalization Method for Orthogonal Frequency Division Multiplexing System |
-
2017
- 2017-08-15 CN CN201710694666.5A patent/CN107623647B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102624419A (zh) * | 2012-04-23 | 2012-08-01 | 西安电子科技大学 | 突发直接序列扩频***的载波同步方法 |
CN103401661A (zh) * | 2013-08-08 | 2013-11-20 | 电子科技大学 | 一种基于mimo雷达通信的一体化编解码方法 |
CN104852875A (zh) * | 2015-03-26 | 2015-08-19 | 中国人民解放军理工大学通信工程学院卫星通信军队重点实验室 | 高动态大频偏突发信号的频偏估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107623647A (zh) | 2018-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107623647B (zh) | 一种基于离散导频辅助的载波同步方法 | |
US10778496B2 (en) | OFDM system synchronization tracking method and apparatus | |
CN100448184C (zh) | 正交频分多路复用***的接收机 | |
CN111683034B (zh) | 一种基于ofdm的大多普勒无线通信时频同步方法 | |
CN110300079B (zh) | 一种msk信号相干解调方法及*** | |
CN101827052B (zh) | 一种用于lte***时间同步和频率同步的方法及装置 | |
CN107241286A (zh) | 一种结合频域和时域估计的载波同步方法 | |
US8422614B2 (en) | Methods and apparatus for determining timing in a wireless communication system | |
CN106059973B (zh) | 频偏估计方法和*** | |
CN109361634B (zh) | 接收机载波频偏的补偿方法及*** | |
CN105763500A (zh) | 一种连续相位调制信号的频偏、时延和相偏的联合同步方法 | |
CN111147123A (zh) | 一种低轨卫星宽带ofdm通信***的载波同步方法 | |
EP2736208B1 (en) | Frequency offset estimation | |
JP5219528B2 (ja) | 無線受信装置および無線受信方法 | |
CN111131114A (zh) | 前导符号的接收方法及装置 | |
CN105610755B (zh) | 一种突发信号的频偏估计方法及装置 | |
CN101534287A (zh) | 移动通信***中进行载波频率偏移校正的方法和装置 | |
CN114039829B (zh) | 一种基于导频辅助的相位恢复方法及*** | |
CN103795676A (zh) | 一种基于循环累积的多普勒频偏估计方法 | |
CN107707499B (zh) | 一种Alpha稳定分布噪声下OFDM信号调制参数估计方法 | |
CN111131110B (zh) | 适合突发通信逆序组合的判决引导载波恢复***及方法 | |
Gudovskiy et al. | A novel nondata-aided synchronization algorithm for MSK-type-modulated signals | |
CN111064687B (zh) | 一种频率和相位补偿的方法及装置 | |
CN110943949B (zh) | 一种基于sc-fde传输体制的宽带角跟踪方法 | |
CN110445738B (zh) | 一种基于离散导频辅助的相位估计方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |