CN114915316B - 一种基于频域处理的带限直扩信号数字码跟踪方法 - Google Patents
一种基于频域处理的带限直扩信号数字码跟踪方法 Download PDFInfo
- Publication number
- CN114915316B CN114915316B CN202210345012.2A CN202210345012A CN114915316B CN 114915316 B CN114915316 B CN 114915316B CN 202210345012 A CN202210345012 A CN 202210345012A CN 114915316 B CN114915316 B CN 114915316B
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- signal
- spread spectrum
- band
- local
- 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 60
- 238000001228 spectrum Methods 0.000 title claims abstract description 26
- 238000012545 processing Methods 0.000 title claims abstract description 11
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 238000009825 accumulation Methods 0.000 claims description 11
- 230000010354 integration Effects 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 10
- 238000005314 correlation function Methods 0.000 claims description 8
- 108010003272 Hyaluronate lyase Proteins 0.000 claims description 5
- 230000007480 spreading Effects 0.000 claims description 4
- 238000005316 response function Methods 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- 230000001629 suppression Effects 0.000 abstract description 7
- 230000005764 inhibitory process Effects 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 230000003595 spectral effect Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical group [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- 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
-
- 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
-
- 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
- H04B2001/70706—Spread spectrum techniques using direct sequence modulation using a code tracking loop, e.g. a delay locked loop
-
- 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
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Noise Elimination (AREA)
Abstract
本发明提出了一种基于频域处理的带限直扩信号数字码跟踪方法,能够在不增加复杂度的条件下实现对窄带干扰的抑制。本发明通过将匹配滤波及时延补偿在频域实现,极大降低了计算复杂度,并且由于该方法由于引入了频域信号处理,能够实现在不增加计算复杂度的条件下实现对窄带干扰的抑制,能够满足数字扩频接收机对于性能和灵活性的需求。
Description
技术领域
本发明通信技术领域,具体涉及一种基于频域处理的带限直扩信号数字码跟踪方法。
背景技术
直接序列扩频系统在通信、导航、雷达等领域都有着广泛的应用。码同步是直接序列扩频接收机的重要组成部分,其中包括码捕获与码跟踪。码捕获部分将本地扩频信号与接收信号的码时延缩小到一个码片间隔内,码跟踪通过环路反馈系统进一步缩小本地信号与接收信号的码时延差。
传统的直接序列扩频系统一般采用矩形波形,其为时间有限波形,本地信号可以通过查表得方式形成,复杂度低,然而,矩形波形频谱效率较低,而且边带衰减较小,容易对其它频带信号形成干扰,因此,基于频带有限波形如升余弦波形的扩频系统,由于其具有较高的频谱效率和边带衰减,受到了广泛的研究。
针对带限扩频信号的码跟踪,现有的方法可以分为两种,一种是时域匹配滤波之后进行采样的半模拟半数字方式,另一种是以全数字方式实现匹配滤波及后续处理。数字化接收机是通信领域的发展方向,然而现存的全数字带限直扩信号码跟踪方法计算复杂度高,处理不够灵活,限制了数字带限扩频接收机的进一步发展。传统的数字带限直扩信号码跟踪方法将原半模拟半数字的解决方法推广到了全数字领域,降低了系统的复杂度,推动了全数字带限扩频接收机的发展,然而,随着设备的小型化和对降低计算复杂度、提高接收机在复杂环境下信号处理的灵活性的要求,现有的数字码跟踪方法逐渐不能满足任务的需求,主要体现在信号处理复杂度高,需要的计算资源较多,以及在频域干扰环境下,性能下降明显,需要添加独立干扰抑制部分。
发明内容
有鉴于此,本发明提出了一种基于频域处理的带限直扩信号数字码跟踪方法,能够在不增加复杂度的条件下实现对窄带干扰的抑制。
为实现上述目的,本发明的一种基于频域处理的带限直扩信号数字码跟踪方法,将匹配滤波及时延补偿在频域实现。
其中,频域实现具体方式为:采用移位反馈寄存器生成的伪码序列;整数时延补偿部分通过对接收信号进行移位实现。
其中,设采样周期Ts=NTc/2(N+1),令扩展因子M=N,考虑第i次迭代,经过整数时延补偿的输入信号为:
其中Tc为码片周期,ai为数据符号,|m|N表示m对N取模,N为伪码长度;gT为发射端成型脉冲,θ为未经同步的载波相位,服从[-π,π)上的均匀分布;ρ=2(N+1)/N,为整数部分;d表示时延;wn,i为低通等效的复高斯白噪声;
其中,采用了频域积分累加操作得到输出。
其中,所述积分累加的输出为:
其中P为积分累加点数,则积分累加输出的逆傅里叶变换为yn,i=IFFT[Yk,i],其中yn,i为第i次迭代输入信号与本地信号的相关函数。
其中,采用非相关误差鉴别器,则输出的误差信号为ei=|y1,i|2-|y-1,i|2,得到伪码跟踪环路的动态方程为:
有益效果:
本发明通过将匹配滤波及时延补偿在频域实现,极大降低了计算复杂度,并且由于该方法由于引入了频域信号处理,能够实现在不增加计算复杂度的条件下实现对窄带干扰的抑制,能够满足数字扩频接收机对于性能和灵活性的需求。
附图说明
图1为本发明方法基于的带限直扩信号频域码跟踪结构示意图。
图2为本发明频域跟踪方法环路时间误差理论与仿真对比示意图。
图3为本发明频域与时域跟踪方法环路时间误差性能对比示意图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明通过在频域完成匹配滤波及延迟补偿,实现了比原有时域方法更低的计算复杂度和相同的跟踪性能,能够在不增加复杂度的条件下实现对窄带干扰的抑制,将信号处理方法扩展到了频域,提高了方法的灵活性。与传统时域数字直扩信号码跟踪方法相比,本发明计算复杂度更低,在窄带干扰环境中跟踪性能更强。
本发明的信号模型中,发射端考虑采用根升余弦波形进行脉冲成型,则调制信号为:
其中,ai为数据符号,{m}M表示m/M的整数部分,M为扩展因子;ci=±1为扩频码,|m|N表示m对N取模,N为伪码长度;gT(t)为发射端成型脉冲,它的频率响应为其中Tc为码片周期,GN(f)为奈奎斯特升余弦滤波器的频域响应。由于调制信号s(t)为循环平稳随机过程,则s(t)的功率谱密度为Ss(f)=|GT(f)|2/Tc,则调制信号s(t)的平均功率为Ps=1。
考虑时延和高斯白噪声的影响,接收端采样之后的信号可以表示为
其中,Ts为采样周期,θ为未经同步的载波相位,服从[-π,π)上的均匀分布;τ为经过捕获过程之后剩余伪码时延;w(nTs)为低通等效的复高斯白噪声,功率谱密度为Sw(f)=N0/P,P为中频信号功率。
本发明提出的频域带限直扩信号码跟踪结构如图1所示,本发明采用移位反馈寄存器生成的伪码序列如m序列、Gold序列等,为了使FFT运算操作效率最大化的同时满足奈奎斯特采样定理,设采样周期Ts=NTc/2(N+1),为了简化分析,令扩展因子M=N,当扩展因子为伪码长度的整数倍时,能够自然地以相干累积的方式进行扩展。为了保证在FFT时间块内数据符号近似保持不变,整数时延补偿部分通过对接收信号进行移位实现。考虑第i次迭代,经过整数时延补偿的输入信号可以表示为:
输入信号经过FFT运算后为Xk,i=FFT[rn,i],本地初始频域扩频信号为其中gR(t)接收端匹配滤波脉冲,它的频率响应根据频域分数时延滤波器理论,归一化时延为D的2(N+1)点分数时延滤波器可以直接在频域实现,表达式为:
则第i次迭代的本地频域共轭信号为接收信号与本地信号的相关函数的频域表达式为Zk,i=Xk,iLk。因为在跟踪过程之前,伪码的粗时延同步已经完成,因此时域相关函数的有效值被限制在一个相较于伪码周期较小的范围之内,因此本方案采用了频域积分累加操作来提高运算效率,则积分累加的输出为:
其中P为积分累加点数,则积分累加输出的逆傅里叶变换为yn,i=IFFT[Yk,i],其中yn,i为第i次迭代输入信号与本地信号的相关函数。方案采用非相关误差鉴别器,则输出的误差信号为
ei=|y1,i|2-|y-1,i|2 (1.6)
由此,可以得到伪码跟踪环路的动态方程为
其中,*代表了线性卷积操作,Gd为环路增益系数,hi为环路滤波器的冲激响应函数。
对本发明方法进行误差鉴别分析:考虑方案的时域等效形式,则在频域实现的循环相关的时域表达式为:
经过化简,其中的功率谱密度为则相关输出的均值和方差分别为var(zp,i)=NN0/PρTs 3。由于超前支路与滞后支路的时间差为2Ts=2Tc/ρ≈Tc,则由zp的自相关函数零点特性可知z1与z-1不相关,同时在频域理想低通滤波情况下,有ei=|z1,i|2-|z-1,i|2,因此误差鉴别器输出的均值和方差分别为:
由公式(1.7)及ei的线性化,可以得到:
由于最终系统保持稳态,因为有:
对本发明方法进行计算复杂度分析:为了分析所提方法的计算复杂度并于时域方法进行对比,考虑单次跟踪的复数乘法次数,由于两种方法采用相同的环路滤波器,则分析主要包括误差鉴别部分和环路误差补偿部分。考虑实际应用,伪码长度为N,时域跟踪方法的匹配滤波器长度设为32点,延迟滤波器长度设为64点,频域跟踪方法的频域积分长度为P,干扰抑制时域方法设为增加前置FFT干扰抑制部分,频域方法的干扰抑制在FFT之后进行,则单次跟踪的计算复杂度可以由表1表示。
表1 时域与频域跟踪方法计算复杂度对比
在典型取值情况下,如N=1023,P=256,在无窄带干扰抑制时,时域方法的所需乘法次数为196416,频域方法所需的乘法次数为12299,为时域方法的0.0626;在有窄带干扰抑制时,时域方法所需的乘法次数约为2.18×105,频域方法所需乘法次数不变,为时域方法的0.0562。可见,本发明提出的频域跟踪方法较原时域跟踪方法计算复杂度降低,且能够在不增加复杂度的条件下实现窄带干扰抑制。
为了验证性能分析的正确性,设计了仿真程序。其中扩频序列采用长度为1023点的Gold序列,伪码速率为1/Tc=10.23MHz,成型函数为滚降系数为0.25的根升余弦函数,数据调制为BPSK符号,速率为1/T=10kHz,采样速率为1/Ts=20.48MHz,噪声为加性高斯白噪声,频域积分累加点数P=256,环路滤波器采用有源比例积分器,设置环路带宽分别分100Hz和10Hz,仿真曲线如图2所示。其中横坐标为中频接收信号信噪比,SNR=PTs/N0,纵坐标为码片的环路方差,理论曲线由公式(1.17)计算得出,σ2=E[|εi|2]/ρ2,仿真曲线由蒙特卡洛仿真得出,每个信噪比仿真200次求平均,每次仿真时间为1s。由图2可见,仿真结果与理论分析结果有较好的重合,验证了本发明分析方法的正确性。
为了对比本发明频域跟踪方法与之前提出的时域跟踪方法的性能,设计了蒙特卡洛仿真实验,同样采用长度为1023点Gold扩频序列,伪码速率为1/Tc=10.23MHz,BPSK符号速率为1/T=10kHz,时域跟踪方法采样率设为1/T′s=20.46MHz,频域跟踪方法的采样率设为1/Ts=20.48MHz,两者采用相同的环路滤波器及环路带宽,单次仿真时间都为1s,每个信噪比仿真200次求平均,得到的仿真曲线如图3所示。由图可知,在不同的环路带宽及信噪比条件下,本发明提出的频域跟踪方法与之前的时域方法相比,性能几乎没有损失,验证了本发明提出方法的有效性。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于频域处理的带限直扩信号数字码跟踪方法,其特征在于,将匹配滤波及时延补偿在频域实现;
频域实现具体方式为:采用移位反馈寄存器生成的伪码序列;整数时延补偿部分通过对接收信号进行移位实现;
设采样周期Ts=NTc/2(N+1),令扩展因子M=N,考虑第i次迭代,经过整数时延补偿的输入信号为:
其中Tc为码片周期,ai为数据符号,j为虚部符号,c为扩频码,|m|N表示m对N取模,N为伪码长度;gT为发射端成型脉冲,θ为未经同步的载波相位,服从[-π,π)上的均匀分布;ρ=2(N+1)/N,为整数部分;d表示时延;wn,i为低通等效的复高斯白噪声;
2.如权利要求1所述的方法,其特征在于,采用了频域积分累加操作得到输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210345012.2A CN114915316B (zh) | 2022-03-31 | 2022-03-31 | 一种基于频域处理的带限直扩信号数字码跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210345012.2A CN114915316B (zh) | 2022-03-31 | 2022-03-31 | 一种基于频域处理的带限直扩信号数字码跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114915316A CN114915316A (zh) | 2022-08-16 |
CN114915316B true CN114915316B (zh) | 2023-03-24 |
Family
ID=82762572
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210345012.2A Active CN114915316B (zh) | 2022-03-31 | 2022-03-31 | 一种基于频域处理的带限直扩信号数字码跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114915316B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116073855B (zh) * | 2022-10-08 | 2023-10-03 | 中国人民解放军战略支援部队航天工程大学 | 一种基于频域匹配滤波和时延调整的bl-dsss信号码跟踪方法 |
CN115940992B (zh) * | 2022-11-16 | 2023-10-03 | 中国人民解放军战略支援部队航天工程大学 | 一种基于频域子空间原理的bl-dsss信号码跟踪方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6975673B1 (en) * | 1998-07-14 | 2005-12-13 | Axonn, L.L.C. | Narrow-band interference rejecting spread spectrum radio system and method |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6058101A (en) * | 1997-06-11 | 2000-05-02 | Industrial Technology Research Institute | Synchronization method and system for a digital receiver |
CN102065043B (zh) * | 2010-12-09 | 2013-01-30 | 北京理工大学 | 一种高速通信系统中的频域并行解调方法 |
CN112804173A (zh) * | 2020-12-29 | 2021-05-14 | 重庆电子工程职业学院 | 高动态大频偏载波跟踪方法 |
US11201604B1 (en) * | 2021-02-24 | 2021-12-14 | North China Power Electric University | Resampling algorithm based on window function |
-
2022
- 2022-03-31 CN CN202210345012.2A patent/CN114915316B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6975673B1 (en) * | 1998-07-14 | 2005-12-13 | Axonn, L.L.C. | Narrow-band interference rejecting spread spectrum radio system and method |
Also Published As
Publication number | Publication date |
---|---|
CN114915316A (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114915316B (zh) | 一种基于频域处理的带限直扩信号数字码跟踪方法 | |
CN113238261B (zh) | 低轨卫星扩频通信体制信号捕获跟踪系统 | |
CN107659334B (zh) | 对多载波扩频信号使用滤波器组的方法和装置 | |
CN107171693B (zh) | 基于波形重构的宽带干扰抑制方法 | |
CN112260977A (zh) | 一种动态卫星通信系统信号跟踪方法及装置 | |
CN108377158B (zh) | 一种实现扩频信号的多频段分割和聚合方法 | |
CN109450488A (zh) | 窄带干扰下带脉冲成形扩频信号伪码周期盲估计 | |
CN104618278A (zh) | 一种基于谱相关的多用户tddm-boc信号的伪码速率估计方法 | |
CN112039816A (zh) | 一种窄带物联网系统下行同步方法 | |
CN113922843B (zh) | 通信系统信号同步方法及系统 | |
CN112803968B (zh) | 一种无人机机载测控方法 | |
Lu et al. | Efficiently all-digital code tracking for band-limited DSSS systems | |
CN112839001B (zh) | 一种无人机机载测控终端 | |
CN107621643B (zh) | 一种适用于导航信号质量评估的相关域参数精确解算方法 | |
CN103605139B (zh) | 适用于gnss接收机的载波频率和相位估计方法及系统 | |
US10581481B1 (en) | Communication device, spread-spectrum receiver, and related method using normalized matched filter for improving signal-to-noise ratio in harsh environments | |
Pan et al. | Analysis of an improved acquisition method for high-dynamic BOC signal | |
CN104767704B (zh) | 一种ofdm信道估计装置和方法 | |
CN115865127A (zh) | 一种直扩信号的参数估计和解调方法 | |
CN113630152B (zh) | 引导式数字抗截获抗干扰装置及方法 | |
Shen et al. | Improving symbol rate estimation accuracy by bandpass filter bank | |
Sun et al. | A Blind Estimation Method of QPSK/OQPSK Symbol Rate | |
Shen et al. | Blind symbol rate estimation based on BPF bank | |
CN110138382A (zh) | 一种基于等效扩频跟踪的自适应时频交叠信号分离方法 | |
CN116073855B (zh) | 一种基于频域匹配滤波和时延调整的bl-dsss信号码跟踪方法 |
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 |