CN114268341A - 基于循环平稳特性的直接序列扩频信号伪码速率估计方法 - Google Patents
基于循环平稳特性的直接序列扩频信号伪码速率估计方法 Download PDFInfo
- Publication number
- CN114268341A CN114268341A CN202111576310.4A CN202111576310A CN114268341A CN 114268341 A CN114268341 A CN 114268341A CN 202111576310 A CN202111576310 A CN 202111576310A CN 114268341 A CN114268341 A CN 114268341A
- Authority
- CN
- China
- Prior art keywords
- frequency
- spread spectrum
- direct sequence
- sequence spread
- spectrum signal
- 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
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 53
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 20
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 2
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
本发明公开了一种基于循环平稳特性的直接序列扩频信号伪码速率估计方法,包括以下步骤:S10、获取直接序列扩频信号数据;S20、估计直接序列扩频信号的载频,得到估计载频;S30、利用估计载频计算直接序列扩频信号频率等于载频的循环谱截面;S40、根据估计载频,得到直接序列扩频信号数据的伪码符号频率的预估值,利用伪码符号频率预估值,并结合循环谱截面上循环频率大于零部分的幅度特征,估计出直接序列扩频信号的伪码速率。该方法解决了低信噪比条件下,以及由于循环频率分辨率和频率分辨率设置不合理导致循环谱截面出现虚假线谱时,难以估计出真实的符号频率对应的线谱位置问题,适合工程上对直接序列扩频信号伪码速率进行快速稳健地估计。
Description
技术领域
本发明属于信号处理领域,具体涉及一种基于循环平稳特性的直接序列扩频信号伪码速率估计方法。
背景技术
直接序列扩频信号具有伪随机编码调制和信号相关处理两大特点,使其存在许多优点,如抗噪声、抗干扰、抗多径衰落等。直接序列扩频信号在诸多领域,如保密通信、多址通信、卫星导航定位中得到了广泛的应用。伪码速率是直接序列扩频信号的一个重要参数,伪码速率的稳定准确估计,是正确解调直接序列扩频信号的一个重要前提,因此直接序列扩频信号的伪码估计问题在军事侦察中具有重要的实际意义。
目前工程上常用的信号码速率估计方法主要有:瞬时频率法、延迟相乘法、小波变换法以及循环平稳法等。其中瞬时频率法以及延迟相乘法运算量小,计算速度快,但是抗噪声能力较差,估计精度一般。基于小波变换法估计精度较高,抗噪声能力较强,但是选用不同的尺度因子以及小波函数时估计性能差异有时较大,如何选取合适的尺度因子以及小波函数至今尚未解决,因此该方法鲁棒性较差。基于循环平稳的方法,抗噪声能力以及鲁棒性较强,估计精度高,因此得到了较广泛的应用,但是在噪声为非白噪声且信噪比较低,以及在循环平稳的方法中计算信号的循环谱时,循环频率以及频率分辨率选取不当出现虚假线谱时,基于循环平稳的方法估计性能将严重下降,有时甚至失效。
发明内容
本发明的目的在于提供一种基于循环平稳特性的直接序列扩频信号伪码速率估计方法,该方法解决了低信噪比条件下,以及由于循环频率分辨率和频率分辨率设置不合理导致循环谱截面出现虚假线谱时,难以估计出真实的符号频率对应的线谱位置问题,且运算速度快,鲁棒性强。
实现本发明的技术解决方案为:一种基于循环平稳特性的直接序列扩频信号伪码速率估计方法,包括以下步骤:
S10、获取直接序列扩频信号数据。
S20、估计直接序列扩频信号的载频,得到估计载频。
S30、利用估计载频计算直接序列扩频信号频率等于载频的循环谱截面。
S40、根据估计载频,得到直接序列扩频信号数据的伪码符号频率的预估值,利用伪码符号频率预估值,并结合循环谱截面上循环频率大于零部分的幅度特征,估计出直接序列扩频信号的伪码速率。
本发明与现有技术相比,其显著优点在于:本发明解决了低信噪比条件下,以及由于循环频率分辨率和频率分辨率设置不合理导致循环谱截面出现虚假线谱时,难以估计出真实的符号频率对应的线谱位置问题,且运算速度快,鲁棒性强,同时不需要存储中间数据,适合工程上对直接序列扩频信号伪码序列进行快速稳健估计。
附图说明
图1为本发明的方法流程图。
图2为根据本发明实施例的直接序列扩频信号实际的海试数据。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
另外,本发明各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在,也不在本发明要求的保护范围指内。
下面将结合本设计实例对具体实施方式、以及本次发明的技术难点、发明点进行进一步介绍。
目前工程上常用的信号码速率估计方法主要有:瞬时频率法、延迟相乘法、小波变换法以及循环平稳法等。其中瞬时频率法以及延迟相乘法运算量小,计算速度快,但是抗噪声能力较差,估计精度一般。基于小波变换法估计精度较高,抗噪声能力较强,但是选用不同的尺度因子以及小波函数时估计性能差异有时较大,如何选取合适的尺度因子以及小波函数至今尚未解决,因此该方法鲁棒性较差。基于循环平稳的方法,抗噪声能力以及鲁棒性较强,估计精度高,因此得到了较广泛的应用,但是在噪声为非白噪声且信噪比较低,以及在循环平稳的方法中计算信号的循环谱时,循环频率以及频率分辨率选取不当出现虚假线谱时,基于循环平稳的方法估计性能将严重下降,有时甚至失效。而本发明的直接序列扩频信号伪码速率估计算法中,使用预估符号频率,并结合循环谱截面的幅度特征,解决了低信噪比条件下,以及由于循环频率分辨率和频率分辨率设置不合理导致循环谱截面出现虚假线谱时,难以估计出真实的符号频率对应的线谱位置问题,同时本发明方法不需要估计噪声的循环谱背景,适合工程上对直接序列扩频信号伪码速率进行快速稳健地估计。
结合图1,本发明所述的一种基于循环平稳特性的直接序列扩频信号伪码速率估计方法,包括以下步骤:
S10、获取直接序列扩频信号数据:
从传感器接收N个采样点的实时采集数据,或者从存储器中提取从检测到信号时刻起始的N个采样点的数据,作为待处理的直接序列扩频信号数据x(n),其中,n=0,1,…,N-1,数据采样频率记为fs,直接序列扩频信号为DSSS-BPSK信号。
S20、估计直接序列扩频信号的载频,得到估计载频,其具体步骤如下:
S21、对直接序列扩频信号数据x(n)进行平方,并去平均得到x1(n),n=0,1,…,N-1即:
x1(n)=x2(n)-mean{x(n)},n=0,1,…,N-1
这里mean{·},表示求平均计算。
S22、对x1(n)做离散傅里叶变换,并计算得到其频谱X1(l),计算过程如下:
这里功率谱只取其非负频率部分,k为P(k)的离散频率索引,|·|为求幅度运算。
S23、搜索P(k)最大值所对应的离散频率索引k0,其计算过程如下:
S24、根据Rife插值算法,估计x1(n)的频率f1,其计算过程如下:
γ表示插值方向;当|X(k0+1)|>|X(k0-1)|时,γ=1,否则γ=-1。
S25、根据x1(n)的频率f1,计算出直接序列扩频信号序列x(n)的估计载频fc,计算过程如下:
S30、利用估计的载频,计算直接序列扩频信号频率等于载频的循环谱截面,其具体步骤如下:
S31、计算直接序列扩频信号数据x(n)的傅里叶变换,即
其中l′为X(l′)的离散频率索引。
S40、根据估计载频,得到直接序列扩频信号数据的伪码符号频率的预估值,利用伪码符号频率预估值,并结合循环谱截面上循环频率大于零部分的幅度特征,估计出直接序列扩频信号的伪码速率,具体如下:
S41、根据直接序列扩频信号序列x(n)的估计载频fc,计算出x(n)的伪码的符号频率的预估值f2(u),其计算过程如下:
其中,u为设置的直接序列扩频信号序列x(n)的伪码序列的阶数,U为设置的直接序列扩频信号序列x(n)的伪码序列的最大阶数,根据实际的应用环境选取;
S43、寻找f2(u)中满足如下公式的值,记为f2s(i),i=0,1,…I:
且满足,
这里频率间隔Δf=ρf/N,需要平滑的频点个数ρ为正整数,一般取10,临时变量ε也为正整数一般取5,设置的临时门限ξ为正数,一般也取10,I为满足上述公式的f2(u)的个数,round{·}为四舍五入运算,同时求取f2s(i),i=0,1,…I的最小值f2smin,即:
f2smin=min{f2s(i)},i=0,1,…I
其中min{·}表示求最小值运算。
S44、若I=0,则直接序列扩频信号序列x(n)的伪码符号频率即为f2max,此时x(n)的伪码速率即为1/f2max,否则x(n)的伪码符号频率为f2smin,此时x(n)的伪码速率即为1/f2smin。
仿真中实际海试数据的信号参数为:采样频率为4000Hz,直接序列扩频信号的真实伪码符号频率为37.5Hz,载频为300Hz;
图2所示为直接序列扩频信号实际的海试时域数据;图3为计算得出的图2中实际直接序列扩频海试数据的循环谱截面。图4为根据图3中的循环谱截面谱峰特征,估计的20帧直接序列扩频信号的伪码符号频率值,其倒数即为直接序列扩频信号的伪码速率值。
从实施例的结果可以看出,本发明估计方法在信噪比低以及噪声背景复杂的海试数据中,仍能正确地估计出直接序列扩频信号的伪码速率,适用于快速稳健估计直接序列扩频信号伪码速率的场合。
Claims (6)
1.一种基于循环平稳特性的直接序列扩频信号伪码速率估计方法,其特征在于,包括以下步骤:
S10、获取直接序列扩频信号数据,转入步骤S20;
S20、估计直接序列扩频信号的载频,得到估计载频,转入步骤S30;
S30、利用估计载频计算直接序列扩频信号频率等于载频的循环谱截面,转入步骤S40;
S40、根据估计载频,得到直接序列扩频信号数据的伪码符号频率的预估值,利用伪码符号频率预估值,并结合循环谱截面上循环频率大于零部分的幅度特征,估计出直接序列扩频信号的伪码速率。
2.根据权利要求1所述的一种基于循环平稳特性的直接序列扩频信号伪码速率估计方法,其特征在于,所述步骤S10中,获取直接序列扩频信号数据,具体如下:
从传感器接收N个采样点的实时采集数据,或从存储器中提取从检测到信号时刻起始的N个采样点的数据,作为待处理的直接序列扩频信号数据x(n),其中,n=0,1,…,N-1,数据采样频率记为fs,直接序列扩频信号为DSSS-BPSK信号。
3.根据权利要求2所述的基于循环平稳特性的直接序列扩频信号伪码速率估计方法,其特征在于,所述步骤S20中,估计直接序列扩频信号的载频,得到估计载频,具体如下:
S21、对直接序列扩频信号数据x(n)进行平方,并去平均得到x1(n),n=0,1,…,N-1,即:
x1(n)=x2(n)-mean{x(n)}
这里mean{·},表示求平均计算;
S22、对x1(n)做离散傅里叶变换,并计算得到其频谱X1(l),计算过程如下:
这里功率谱只取其非负频率部分,k为P(k)的离散频率索引,|·|为求幅度运算;
S23、搜索P(k)最大值所对应的离散频率索引k0,其计算过程如下:
S24、根据Rife插值算法,估计x1(n)的频率f1,其计算过程如下:
γ表示插值方向;当|X(k0+1)|>|X(k0-1)|时,γ=1,否则γ=-1;
S25、根据x1(n)的频率f1,计算出直接序列扩频信号序列x(n)的估计载频fc,计算过程如下:
5.根据权利要求4所述的一种基于循环平稳特性的直接序列扩频信号伪码速率估计方法,其特征在于,所述步骤S40中,根据估计载频,得到直接序列扩频信号数据的伪码符号频率的预估值,利用伪码符号频率预估值,并结合循环谱截面上循环频率大于零部分的幅度特征,估计出直接序列扩频信号的伪码速率,具体如下:
S41、根据直接序列扩频信号序列x(n)的估计载频fc,计算出x(n)的伪码的符号频率的预估值f2(u),其计算过程如下:
其中,u为设置的直接序列扩频信号序列x(n)的伪码序列的阶数,U为设置的直接序列扩频信号序列x(n)的伪码序列的最大阶数,根据实际的应用环境选取;
S43、寻找f2(u)中满足如下公式的值,记为f2s(i),i=0,1,…I:
f2(u)<f2max,u=1,2,…,U
且满足,
这里频率间隔Δf=ρf/N,ρ为需要平滑的频点个数,ε为临时变量,取值为正整数,设置的临时门限ξ为正数,I为满足上述公式的f2(u)的个数,round{·}为四舍五入运算,同时求取f2s(i),i=0,1,…I的最小值f2smin,即:
f2smin=min{f2s(i)},i=0,1,…I
其中min{·}表示求最小值运算;
S44、若I=0,则直接序列扩频信号序列x(n)的伪码符号频率即为f2max,此时x(n)的伪码速率即为1/f2max,否则x(n)的伪码符号频率为f2smin,此时x(n)的伪码速率即为1/f2smin。
6.根据权利要求4所述的一种基于循环平稳特性的直接序列扩频信号伪码速率估计方法,其特征在于:步骤S43中,ρ取10,ε取5,ξ取10。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111576310.4A CN114268341B (zh) | 2021-12-21 | 2021-12-21 | 基于循环平稳特性的直接序列扩频信号伪码速率估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111576310.4A CN114268341B (zh) | 2021-12-21 | 2021-12-21 | 基于循环平稳特性的直接序列扩频信号伪码速率估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114268341A true CN114268341A (zh) | 2022-04-01 |
CN114268341B CN114268341B (zh) | 2024-02-13 |
Family
ID=80828540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111576310.4A Active CN114268341B (zh) | 2021-12-21 | 2021-12-21 | 基于循环平稳特性的直接序列扩频信号伪码速率估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114268341B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070091987A1 (en) * | 2005-10-26 | 2007-04-26 | Samsung Electro-Mechanics Co., Ltd. | Apparatus and method for detecting code of direct sequence spread spectrum signal |
CN101917239A (zh) * | 2010-09-08 | 2010-12-15 | 中国人民解放军信息工程大学 | 基于频域平滑的获取循环谱α截面的方法和系统 |
CN102158248A (zh) * | 2011-03-04 | 2011-08-17 | 西安电子科技大学 | 基于循环谱理论的ds/fh扩频信号参数估计方法 |
CN105553635A (zh) * | 2015-12-03 | 2016-05-04 | 重庆邮电大学 | 基于循环自相关的boc信号参数盲估计方法 |
CN106443178A (zh) * | 2016-09-08 | 2017-02-22 | 东南大学 | 一种基于IQuinn‑Rife综合的正弦信号频率估计方法 |
US20170317712A1 (en) * | 2016-04-28 | 2017-11-02 | Raytheon Company | Rate line suppression using chaotic spreading codes |
CN109150236A (zh) * | 2018-08-01 | 2019-01-04 | 东南大学 | 一种基于变步长leap神经网络的直扩信号伪码序列估计方法 |
CN109361632A (zh) * | 2018-10-25 | 2019-02-19 | 重庆邮电大学 | 基于谱相关的多径boc信号的参数盲估计方法 |
CN109379310A (zh) * | 2018-09-13 | 2019-02-22 | 东南大学 | 一种基于Rife-Quinn综合的MPSK信号载频估计方法 |
-
2021
- 2021-12-21 CN CN202111576310.4A patent/CN114268341B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070091987A1 (en) * | 2005-10-26 | 2007-04-26 | Samsung Electro-Mechanics Co., Ltd. | Apparatus and method for detecting code of direct sequence spread spectrum signal |
CN101917239A (zh) * | 2010-09-08 | 2010-12-15 | 中国人民解放军信息工程大学 | 基于频域平滑的获取循环谱α截面的方法和系统 |
CN102158248A (zh) * | 2011-03-04 | 2011-08-17 | 西安电子科技大学 | 基于循环谱理论的ds/fh扩频信号参数估计方法 |
CN105553635A (zh) * | 2015-12-03 | 2016-05-04 | 重庆邮电大学 | 基于循环自相关的boc信号参数盲估计方法 |
US20170317712A1 (en) * | 2016-04-28 | 2017-11-02 | Raytheon Company | Rate line suppression using chaotic spreading codes |
CN106443178A (zh) * | 2016-09-08 | 2017-02-22 | 东南大学 | 一种基于IQuinn‑Rife综合的正弦信号频率估计方法 |
CN109150236A (zh) * | 2018-08-01 | 2019-01-04 | 东南大学 | 一种基于变步长leap神经网络的直扩信号伪码序列估计方法 |
CN109379310A (zh) * | 2018-09-13 | 2019-02-22 | 东南大学 | 一种基于Rife-Quinn综合的MPSK信号载频估计方法 |
CN109361632A (zh) * | 2018-10-25 | 2019-02-19 | 重庆邮电大学 | 基于谱相关的多径boc信号的参数盲估计方法 |
Non-Patent Citations (10)
Title |
---|
余敬东,郭伟,李牧: "直接序列扩频通信的循环谱检测和参数估计技术", no. 3 * |
史建锋;沈辉;: "基于循环谱包络的BPSK信号码元速率改进估计算法", no. 05 * |
史建锋;王可人;: "基于循环谱包络的多径直扩信号参数估计", no. 08 * |
周钱;臧国珍;宋慧颖;: "基于循环谱的DSSS信号检测方法研究", no. 11 * |
汪赵华;陈昊;郭立;: "基于频域平滑循环周期图法的直接序列扩频信号的参数估计", no. 05 * |
王利萍;邵高平;殷小龙;孙军明;: "一种基于Hilbert-Huang变换和循环谱的直扩信号参数估计方法", no. 06 * |
王钢;江东;: "基于谱相关的直扩信号载频及伪码速率估计技术", no. 13 * |
谢皓;陈光;: "基于循环谱的GPS信号快速捕获算法", no. 11 * |
郑鹏;张鑫;刘锋;: "直扩信号载频、码元速率循环谱估计改进算法", no. 09 * |
郭小芳;李懿;李锋;: "基于改进循环谱的直扩信号检测方法研究", no. 32 * |
Also Published As
Publication number | Publication date |
---|---|
CN114268341B (zh) | 2024-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109802912B (zh) | 宽带无线通信系统的同步方法、装置、设备及存储介质 | |
CN109379310B (zh) | 一种基于Rife-Quinn综合的MPSK信号载频估计方法 | |
CN108845973B (zh) | 一种基于改进Quinn算法的多普勒频率估计方法 | |
CN110068727B (zh) | 一种基于Candan-Rife综合内插的单频信号频率估计方法 | |
CN101702701A (zh) | 极低信噪比下频率偏移的估计与补偿方法 | |
CN105871413B (zh) | 低信噪比直接序列扩频信号检测方法 | |
CN114280366B (zh) | 基于改进型频率插值算法的正弦信号频率估计方法 | |
Phukan et al. | An algorithm for blind symbol rate estimation using second order cyclostationarity | |
CN106878213B (zh) | 一种lte上行频偏估计的方法 | |
CN109104215B (zh) | 一种基于小波变换的跳频信号码元速率盲估计方法 | |
CN111654308B (zh) | 一种突发扩频弱信号的精密频偏估计方法 | |
CN111158026B (zh) | 一种boc信号的边峰消除方法 | |
Pan et al. | Analysis of an improved acquisition method for high-dynamic BOC signal | |
CN114268341A (zh) | 基于循环平稳特性的直接序列扩频信号伪码速率估计方法 | |
CN110018500B (zh) | 一种基于圆周移位的北斗卫星信号捕获方法 | |
CN116032709B (zh) | 无先验知识fsk信号盲解调和调制特征解析方法及装置 | |
CN105656511B (zh) | 一种适应于有频偏和低信噪比环境下的差分相关捕获方法 | |
CN110048976B (zh) | 面向中频的短码直扩信号数组盲解扩方法及装置 | |
CN111539323A (zh) | 一种循环前缀线性调频信号的频率估计方法与装置 | |
CN115051899A (zh) | 频偏估计方法及装置、计算机可读存储介质 | |
CN115643139A (zh) | 一种抗频偏的突发信号检测系统及检测方法 | |
CN105681240B (zh) | 一种适应于低信噪比环境下的高精度频偏估计方法 | |
CN117176518B (zh) | 频偏估计方法、装置、电子设备和存储介质 | |
Wang et al. | Assisted GPS signal acquisition algorithm based on the joint code-frequency search | |
CN116055262B (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 |