CN105553635B - 基于循环自相关的boc信号参数盲估计方法 - Google Patents
基于循环自相关的boc信号参数盲估计方法 Download PDFInfo
- Publication number
- CN105553635B CN105553635B CN201510881642.1A CN201510881642A CN105553635B CN 105553635 B CN105553635 B CN 105553635B CN 201510881642 A CN201510881642 A CN 201510881642A CN 105553635 B CN105553635 B CN 105553635B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- boc signal
- boc
- cycle
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L7/00—Arrangements for synchronising receiver with transmitter
- H04L7/0054—Detection of the synchronisation error by features other than the received signal transition
- H04L7/007—Detection of the synchronisation error by features other than the received signal transition detection of error based on maximum signal power, e.g. peak value, maximizing autocorrelation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/7073—Synchronisation aspects
- H04B1/7087—Carrier synchronisation aspects
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/69—Spread spectrum techniques
- H04B1/707—Spread spectrum techniques using direct sequence modulation
- H04B1/709—Correlator structure
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
- H04L1/0045—Arrangements at the receiver end
- H04L1/0046—Code rate detection or code type detection
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明请求保护一种基于循环自相关函数的BOC信号参数盲估计方法,属于信号处理技术领域。通过求取BOC信号的循环自相关函数,然后提取平均循环自相关函数延时τ=0的切面,搜索该切面的最大谱峰,最大谱峰所对应的循环频率即为所要估计的载频;提取平均循环自相关函数延时τ≠0的切面,搜索该切面距离零频最近的两个幅值大小不同的谱峰,根据它们所对应的循环频率可以实现伪码速率及副载波速率的估计。同时通过对多段信号的循环自相关函数进行累加平均,可以实现降低噪声和精确估计的目的。本方法可以在低信噪比下较准确地估计BOC信号的载频、副载波速率和伪码速率,从而对该信号的后续处理以及细微特征分析具有重要意义。
Description
技术领域
本发明涉及导航通信信号处理,具体为一种基于循环自相关的Binary OffsetCarrier(BOC)信号参数盲估计问题。
背景技术
Binary Offset Carrier(BOC)信号将原有的PSK(相位键控)调制信号的功率谱位于中心频率处的峰值搬移到中心频率两侧,使得原来一个中心主峰变成后来的位于两侧的两个主峰,从而消除了共频带信号之间的干扰,正是由于BOC信号的这种良好的频谱***特性,使得BOC信号广泛应用于各国的导航***中,同时由于该技术可以使***有更好的捕获跟踪性能以及更高的定位导航精度,因此受到各个导航大国的青睐。根据最新Galileo空间信号接口控制文档显示,10种Galileo导航信号中就有8种采用了BOC信号或者衍生型的BOC信号。与伽利略***的设计类似,美国的GPS导航***和我国的“北斗”导航***也采用了BOC调制方式,使得军用信号与民用信号分离传输,提高了军用信号的安全性。随着BOC调制信号在导航***中的逐渐应用,对其检测技术的研究也更加深入。
综上所述,BOC信号将会是未来导航通信甚至移动通信领域中非常重要的一种调制信号,因此对BOC信号的研究也越来越重要,BOC信号的研究主要包括伪码速率,副载波速率,伪码周期以及伪码的提取,这些研究对于信号解调、信息安全、电子对抗以及对信号进行捕获和跟踪都具有重要作用。
目前针对BOC信号的研究文献比较少,且集中在BOC信号的捕获与跟踪上,针对该信号的盲估计问题研究很少;文献“钱博.TDDM-BOC信号参数估计方法.信息与控制,2011”利用平方倍频法和自相关函数的多峰特性提出一种相关检测算法,但该方法中伪码速率和副载波速率的估计受解调性能的影响,且各参数的估计存在传递误差,估计效果不明显,尤其是在低信噪比情况下,估计性能不理想。文献“张天骐.基于谱相关的BOC信号参数估计.华中科技大学学报,2013”利用了循环平稳信号的谱相关理论对BOC信号进行了研究,实现了该信号的参数估计,虽然该方法的估计性能优异,但其计算量较大,很难适用于实际应用中。因此本发明提出基于循环自相关的BOC信号参数盲估计方法。
发明内容
本发明所要解决的技术问题,针对现有技术存在的BOC信号参数估计中传递误差大,低信噪比下估计性能差等缺陷,提出一种基于循环自相关的估计方法,解决了BOC信号参数盲估计的难题。该方法能够较精确的估计出BOC信号的伪码速率,副载波速率及载频。在信号循环自相关实现的过程中,采用累加平均的方法可以达到降低噪声和精确估计的目的。
本发明解决上述技术问题的技术方案是:一种基于循环自相关的BOC信号参数盲估计方法,其步骤在于,以采样频率fs对接收到的BOC信号进行采样,并将采样后的信号以一定长度L(一般L至少包含两个伪码周期的采样点数)进行分段;分别计算每段信号的循环自相关函数;将得到的循环自相关函数进行累加平均求取BOC信号的平均循环自相关函数;提取平均循环自相关函数延时τ=0的切面,搜索该切面的最大谱峰,最大谱峰所对应的循环频率即为所要估计的载频;提取平均循环自相关函数延时τ≠0的切面,搜索该切面距离零频最近的两个幅值大小不同的谱峰,根据它们所对应的循环频率可以实现伪码速率及副载波速率的估计。
BOC信号模型可以表示为:其中SB(t)为基带BOC信号,f0为载波频率,为载波初始相位,根据公式可以求出BOC信号的循环自相关函数通过对BOC信号进行循环自相关处理,可以发现在其循环自相关切面上,载频、副载波码片速率和伪码速率等参数会表现出较强的信号分量,从而可以实现BOC信号的载频,副载波速率和伪码速率的联合估计。
本发明运用循环自相关函数法对BOC信号的参数进行盲估计,分析推导了BOC信号的平均自相关函数,实现了载频、副载波速率和伪码速率的联合估计,克服了相关检测算法中存在的传递误差大,低信噪比下估计性能差等问题,同时循环自相关函数具有叠加特性,可以通过累加平均的方法提高该算法的抗噪声性能。
附图说明
图1本发明BOC信号参数盲估计方法流程图;
图2本发明BOC信号产生框图;
图3本发明BOC信号的循环自相关函数三维图;
图4本发明BOC信号的循环自相关τ=0的切面图;
图5本发明BOC信号的循环自相关τ≠0的切面图;
图6本发明载频估计的性能图;
图7本发明副载波速率估计的性能图;
图8本发明伪码速率估计的性能图;
具体实施方式
以下结合附图和具体实例,对本发明的实施作进一步的描述。
图1所示为本发明伪码速率估计方法流程图,具体步骤:对接收到的BOC信号进行采样,并将采样后的信号以一定长度(一般至少包含两个伪码周期的采样点数)进行分段;分别计算每段信号的循环自相关函数并取绝对值;将循环自相关函数进行累加平均求取BOC信号的平均循环自相关函数;提取平均循环自相关函数延时τ=0的切面,搜索该切面的最大谱峰,最大谱峰所对应的循环频率即为所要估计的载频;提取平均循环自相关函数延时τ≠0的切面,搜索该切面距离零频最近的两个幅值大小不同的谱峰,根据它们所对应的循环频率可以实现伪码速率及副载波速率的估计。
图2所示为BOC信号产生框图。信息数据先经伪码数据调制产生已扩序列,然后将已扩序列调制到方波副载波上,最后调制到载波上实现。
参照BOC信号产生框图,分析可得BOC信号的表达式为
式中,S(t)表示BOC信号;f0为载波频率;为载波初始相位;SB(t)为基带BOC信号,可表示为
其中,Sc(t)为正弦相位产生的方波副载波;{an}为经数据调制后的扩频序列;Tc为伪码码元宽度;uTc(t)是持续时间为Tc且幅度为1的矩形脉冲;可将基带BOC信号进一步表示为
其中,为扩频符号;Ts为副载波码片宽度,是持续时间为Ts且幅度为1的矩形脉冲。
综上所述,BOC信号可以表示为
式中,调制阶数为N=2Rs/Rc,伪码速率为Rc=1/Tc,副载波速率为Rs=1/(2Ts)。
对于循环平稳信号,可以通过延迟积的同步平均方法来提取这种周期性。假设x(t)为一周期平稳过程,其时变自相关函数定义为:
令N→∞(消除随机性),则有
其中,Rx(t,τ)是关于时间t且周期为T0的周期函数,将其做Fourier级数分解有
Fourier系数由下式给出
式中,α为循环频率,称为循环自相关函数,它是α和τ的二元函数。当循环频率α=0时,循环自相关函数即为传统的自相关函数。
由(1)和(8)式,可推出BOC信号的循环自相关函数的表达式为:
其中,为基带BOC信号的循环自相关函数。因此,求出就可以得到
根据式(6),分别令T0=Tc,T0=Ts,可得
由BOC信号的特征以及式(4)(10)和(11),可知基带BOC信号的时变自相关函数是关于时间t且周期分别为Ts和Tc的多周期函数。本文以BOC(1,1)信号(调制阶数N=2)为例,对循环自相关函数进行推导:
1.当周期为Ts时,为了便于分析可将基带BOC信号表示为
式中,am是信息码,扩频码,副载波三者调制后的序列。由其表达式可知其模型与PAM信号类似,而PAM信号的循环自相关函数易得为
其中
因此对于BOC信号,当周期为Ts时,可以得到
故此时基带BOC信号的循环自相关函数(τ≤Ts/2)为
将代入式(9)可得BOC信号的循环自相关函数为
其中,p为整数。
2.当周期为Tc时,基带BOC调制信号表示如(3)式,即
其中,对于BOC(1,1)来说,调制阶数N=2,即Tc=2Ts,则有
即此时
当时延|τ|≤Ts/2时,由式(13)和式(14)得
则
当时延Ts/2<|τ|≤Ts时,同理可得
则
其中,表示载波初始相位,t0表示BOC信号的起始时刻。由式(19)(24)和(26)可知BOC信号的循环自相关函数只在其循环频率处有非零值,且其循环频率为一族离散值,它们仅与信号的伪码速率Rc(1/Tc)、副载波速率Rs(1/2Ts)和载频f0有关,因此可以通过直接搜索循环自相关函数包络的一系列峰值来获取循环频率,再由其循环频率与信号的伪码速率、副载波速率和载频之间的关系对这三个参数进行联合估计。
利用仿真实验对本发明算法的理论推导进行验证,实验参数设置:BOC(1,1)信号,将BOC信号经过高斯白噪声信道,载频f0=4.092MHz,伪码速率为Rc=1.023MHz,伪码周期63,副载波速率Rs=1.023MHz,采样频率fs=16.368MHz。通过计算机仿真可得到BOC信号的循环自相关函数三维图如图3所示。
图4所示为延时τ=0的BOC信号循环自相关切面图,从图中可以看出该切面是关于零频正负轴对称的,且仅有2个幅度相同主峰,主峰所对应的循环频率即为载频的2倍。图5所示为延时τ≠0的BOC信号循环自相关切面图,搜索该切面距离零频最近的两个幅值大小不同的谱峰,离零频最近的谱峰对应的循环频率即为伪码速率,离零频第二近的谱峰对应的循环频率即为副载波码片速率,同时由于是BOC(1,1)信号调制阶数为2,则副载波码片速率为副载波速率的2倍,根据两个谱峰各自所对应的循环频率可以实现伪码速率及副载波速率的估计。
图6、图7和图8分别代表BOC信号的载频、副载波速率、伪码速率估计性能与数据长度的关系,此处估计性能用标准均方根误差(NRMSE)来衡量,从图中可以看出:在其他条件不变的情况下,随着数据长度成倍递增,载频、副载波速率和伪码速率的估计性能可以得到明显的改善。
本发明通过对BOC信号进行循环自相关函数处理以后,在载频,副载波码片速率和伪码速率等参数上表现出较强的信号分量,从该特征可以对BOC信号的载频,副载波速率和伪码速率进行联合估计。同时通过累加平均的方法可以进一步提高该算法的处理增益。该算法计算量相对较小,估计精度较好,克服了谱相关方法非线性运算复杂的问题。本方法可以在低信噪比下较准确地估计BOC信号的多个参数,从而对该信号的后续处理以及细微特征分析(副载波类型识别,信息序列估计,伪码序列估计)具有重要意义。
Claims (1)
1.一种基于循环自相关的BOC信号参数盲估计方法,其步骤在于,以采样频率fs对接收到的BOC信号进行采样,并将采样后的信号以一定长度L至少包含两个伪码周期的采样点数进行分段;分别计算每段信号的循环自相关函数;将得到的循环自相关函数进行累加平均求取BOC信号的平均循环自相关函数;提取平均循环自相关函数延时τ=0的切面,搜索该切面的最大谱峰,最大谱峰所对应的循环频率即为所要估计的载频;提取平均循环自相关函数延时τ≠0的切面,搜索该切面距离零频最近的两个幅值大小不同的谱峰,根据它们所对应的循环频率可以实现伪码速率及副载波速率的估计,具体特征在于,建立BOC信号循环自相关函数为
其中,为基带BOC信号的循环自相关函数,α表示循环频率;f表示频谱频率;f0表示载波频率;同时结合基带BOC信号的时变自相关函数是关于时间t且周期分别为Ts和Tc的多周期函数这一性质,可以求得BOC信号的循环自相关函数,其中Ts为副载波码片宽度,Tc为伪码码元宽度,当周期为Ts时,BOC信号循环自相关函数表示为
其中,表示载波初始相位,t0表示BOC信号的起始时刻;当周期为Tc时,BOC信号循环自相关函数与延时τ的取值范围有关;当时延|τ|≤Ts/2时,BOC信号循环自相关函数可表示为
当时延Ts/2<|τ|≤Ts时,BOC信号循环自相关函数可表示为
从上述分析可知BOC信号的循环自相关函数只在其循环频率处有非零值,且其循环频率为一族离散值,它们仅与信号的伪码速率Rc(1/Tc)、副载波速率Rs(1/2Ts)和载频f0有关,因此可以通过直接搜索循环自相关函数包络的一系列峰值来获取循环频率,再由其循环频率与信号的伪码速率、副载波速率和载频之间的关系对这三个参数进行联合估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510881642.1A CN105553635B (zh) | 2015-12-03 | 2015-12-03 | 基于循环自相关的boc信号参数盲估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510881642.1A CN105553635B (zh) | 2015-12-03 | 2015-12-03 | 基于循环自相关的boc信号参数盲估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105553635A CN105553635A (zh) | 2016-05-04 |
CN105553635B true CN105553635B (zh) | 2019-05-10 |
Family
ID=55832583
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510881642.1A Active CN105553635B (zh) | 2015-12-03 | 2015-12-03 | 基于循环自相关的boc信号参数盲估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105553635B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109714287A (zh) * | 2018-10-18 | 2019-05-03 | 重庆邮电大学 | 多径衰落信道下fbmc信号的符号周期盲估计 |
CN109361632B (zh) * | 2018-10-25 | 2021-03-30 | 重庆邮电大学 | 基于谱相关的多径boc信号的参数盲估计方法 |
CN109982253A (zh) * | 2019-01-29 | 2019-07-05 | 中国航空无线电电子研究所 | 基于无线电多元相关特性的波形加载平台 |
CN112803967B (zh) * | 2020-12-30 | 2022-07-12 | 湖南艾科诺维科技有限公司 | 用于非协同扩频信号的检测及参数估计方法及装置 |
CN114268341B (zh) * | 2021-12-21 | 2024-02-13 | 中国航天科工集团八五一一研究所 | 基于循环平稳特性的直接序列扩频信号伪码速率估计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104199064A (zh) * | 2014-09-12 | 2014-12-10 | 重庆邮电大学 | 基于相关波动的boc信号参数盲估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1598677B1 (en) * | 2004-05-17 | 2007-08-08 | University Technologies International Inc. | BOC signal acquisition and tracking method and apparatus |
-
2015
- 2015-12-03 CN CN201510881642.1A patent/CN105553635B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104199064A (zh) * | 2014-09-12 | 2014-12-10 | 重庆邮电大学 | 基于相关波动的boc信号参数盲估计方法 |
Non-Patent Citations (3)
Title |
---|
BOC信号的盲估计算法研究;阳锐;《中国优秀硕士学位论文全文数据库》;20150715;第44-62页 |
基于谱相关的BOC调制信号参数估计;张天琪等;《华中科技大学学报(自然科学版)》;20130930;第41卷(第9期);第1-3页 |
循环谱相关方法在调制信号参数估计中的应用;孟玲玲等;《燕山大学学报》;20100531;第34卷(第3期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN105553635A (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105553635B (zh) | 基于循环自相关的boc信号参数盲估计方法 | |
US10911178B2 (en) | System and method for blind detection of numerology | |
US9282525B2 (en) | Frequency-domain symbol and frame synchronization in multi-carrier systems | |
CN104618278B (zh) | 一种基于谱相关的多用户tddm‑boc信号的伪码速率估计方法 | |
CN105516045B (zh) | 一种ofdm训练序列构造及同步方法 | |
CN105445767B (zh) | 基于平均模糊函数的boc信号参数盲估计方法 | |
CN105553506B (zh) | 一种长码扩频信号的快速捕获方法及装置 | |
CN105743612A (zh) | 实时盲解调大频移短时突发信号的方法 | |
CN109587089A (zh) | 一种提升无人机信号识别准确度的方法 | |
CN104270167A (zh) | 一种基于多维特征神经网络的信号检测和估计方法 | |
CN109412644A (zh) | 一种直扩msk信号的多普勒频率估计方法 | |
CN105790788A (zh) | 一种直扩msk信号的伪码-多普勒联合捕获方法 | |
WO2021139163A1 (zh) | 射频指纹提取方法、装置、无线接收机及存储介质 | |
CN102158248A (zh) | 基于循环谱理论的ds/fh扩频信号参数估计方法 | |
CN104199064B (zh) | 基于相关波动的boc信号参数盲估计方法 | |
CN106526631B (zh) | 北斗b1频点卫星信号载波频率高精度估计方法 | |
CN107359902A (zh) | 一种基于频率槽和fft的pcma***时延和频差快速粗捕方法 | |
CN102707297A (zh) | 一种compass信号的快速捕获方法 | |
CN106027116B (zh) | 一种基于chirp信号的移动水声通信多普勒系数估计方法 | |
CN105656511B (zh) | 一种适应于有频偏和低信噪比环境下的差分相关捕获方法 | |
CN109660475B (zh) | 一种非合作相位编码水声通信信号自主辨识方法 | |
CN105306077B (zh) | 信号解码方法及装置 | |
CN110290084A (zh) | 一种基于数据频率能量峰值的短波信道盲符号同步方法 | |
CN109361632B (zh) | 基于谱相关的多径boc信号的参数盲估计方法 | |
CN107425904A (zh) | 一种基于频率槽和fft的pcma***时延和频差快速粗捕方法 |
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 |