CN113890563A - 一种直扩信号时频二维内插精估方法 - Google Patents
一种直扩信号时频二维内插精估方法 Download PDFInfo
- Publication number
- CN113890563A CN113890563A CN202111122632.1A CN202111122632A CN113890563A CN 113890563 A CN113890563 A CN 113890563A CN 202111122632 A CN202111122632 A CN 202111122632A CN 113890563 A CN113890563 A CN 113890563A
- Authority
- CN
- China
- Prior art keywords
- module
- frequency
- correlation
- pseudo code
- max
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000001228 spectrum Methods 0.000 title claims abstract description 15
- 238000009825 accumulation Methods 0.000 claims abstract description 57
- 238000005070 sampling Methods 0.000 claims abstract description 48
- 238000001514 detection method Methods 0.000 claims abstract description 19
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 230000001427 coherent effect Effects 0.000 claims abstract description 4
- 230000035508 accumulation Effects 0.000 claims description 55
- 230000000875 corresponding effect Effects 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 7
- 230000033001 locomotion Effects 0.000 claims description 5
- 230000000737 periodic effect Effects 0.000 claims description 5
- 230000002596 correlated effect Effects 0.000 claims description 4
- 230000021615 conjugation Effects 0.000 claims description 3
- 238000005314 correlation function Methods 0.000 claims description 2
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims 1
- 230000008569 process Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 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
- H04B1/7073—Synchronisation aspects
- H04B1/7075—Synchronisation aspects with code phase acquisition
- H04B1/708—Parallel implementation
-
- 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/7075—Synchronisation aspects with code phase acquisition
-
- 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/7075—Synchronisation aspects with code phase acquisition
- H04B1/70751—Synchronisation aspects with code phase acquisition using partial detection
- H04B1/70752—Partial correlation
-
- 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/7075—Synchronisation aspects with code phase acquisition
- H04B1/70751—Synchronisation aspects with code phase acquisition using partial detection
- H04B1/70753—Partial phase search
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种直扩信号时频二维内插精估装置,包括接收机天线、射频前端、存储器模块、并行频率搜索模块、FFT模块、伪码相关模块、非相干积累模块、峰值检测模块、相关积累存储模块以及三角包络拟合模块;射频前端包括下变频单元以及AD采样单元;伪码相关模块包括预存的本地伪码快速傅立叶变换、乘法单元以及逆快速傅立叶变换单元;并行频率搜索模块包含L个频率通道。本发明建立在传统码域并行搜索捕获结构的基础上,实现简单,无需加入额外的复杂结构,保证了捕获速度,在捕获时间受限的情况下能够快速、准确地完成时频二维精估。
Description
技术领域
本发明属于直接序列扩频通信以及伪码相位及多普勒频率精估技术领域,涉及一种基于三角包络拟合的直扩信号时频二维内插精估方法。
背景技术
直接扩频序列(Direct Sequences Spread Spectrum,DSSS)作为一种隐蔽性好、具有一定抗干扰、抗截获能力的信号处理方式,已被广泛应用于卫星通信、深空通信等环境较为恶劣的场景。其信号处理过程通常包括捕获、跟踪、位同步等步骤,其中,捕获过程作为对DSSS信号的第一步处理,其估计效果将直接影响后续跟踪和位同步是否能准确、稳定地工作,因此,捕获过程在DSSS信号处理过程中尤为重要。
DSSS信号的捕获过程需要进行时频二维估计,通过在较大范围内按照一定搜索步进进行搜索,实现载波多普勒频率及伪码相位的粗略估计,为了便于后续跟踪环路的快速、稳定锁定,通常捕获过程的估计精度越高越好。捕获过程中通常面临的主要问题分为两类:第一,由于发射信号的能量有限,且传输距离较远,接收信号的信噪比极低,导致捕获过程的搜索极为困难,捕获之后的估计精度无法满足后续跟踪的入锁要求;第二,由于捕获过程要求捕获时间极短,而较大相对运动带来的高动态将导致捕获搜索范围变大,若要在规定时间内完成搜索,只能降低捕获的估计精度。综上所述,DSSS信号的捕获过程亟需解决的问题为时频二维估计精度与高动态、低信噪比之间的矛盾。
发明内容
(一)发明目的
本发明的目的是:针对信噪比、高动态下传统捕获方法中存在载波多普勒频率和伪码相位估计精度不足的缺陷,提出一种直扩信号时频二维内插精估方法。
(二)技术方案
为了解决上述技术问题,本发明提供一种直扩信号时频二维内插精估装置;本装置包括接收机天线、射频前端、存储器模块、并行频率搜索模块、FFT模块、伪码相关模块、非相干积累模块、峰值检测模块、相关积累存储模块以及三角包络拟合模块;
其中,射频前端包括下变频单元以及AD采样单元;
伪码相关模块包括预存的本地伪码快速傅立叶变换、乘法单元以及逆快速傅立叶变换单元;
并行频率搜索模块包含L个频率通道;
本装置的各模块连接关系如下:
接收机天线与射频前端相连,射频前端与并行频率搜索模块相连,并行频率搜索模块与FFT模块相连,FFT模块和存储器模块均与伪码相关模块相连,伪码相关模块与非相干积累模块相连,非相干积累模块与峰值检测模块相连,峰值检测模块与相关积累存储模块相连,相关积累存储模块与三角包络拟合模块相连;
本装置的各模块功能如下:
接收机天线用于接收信号;
射频前端用于对接收信号进行下变频和采样;
存储器模块用于存储经过FFT变换后一个周期的本地伪码序列;
并行搜索模块用于按固定频率步进产生本地载波,再与接收信号相乘,以剥离多普勒频率,再将此剥离多普勒频率后的信号输出FFT模块;
FFT模块用于对并行搜索模块输出的每一组序列进行快速傅立叶变换,得到伪码相关模块的输入FFT序列;
伪码相关模块用于将FFT模块的输出序列与本地伪码信号的快速傅立叶变换序列相乘,再通过逆快速傅立叶变换得到相关结果,取模值,其中乘法单元调用预存的本地伪码FFT序列,与折叠相加后FFT模块的输出依次相乘,IFFT单元对乘法单元的输出进行逆快速傅立叶变换,取模单元对IFFT单元的输出进行取模运算;
非相干积累模块用于将伪码相关模块输出的结果按照相应位置进行非相干累加,提高输出信噪比;
峰值检测模块用于比较非相干积累模块输出的相关结果,找出其中的最大值,根据其所处的搜索通道位置得到载波多普勒频率估计值fmax,并由峰值检测模块输出的峰值位置得到码相位估计值τmax,作为捕获结果Pmax;
相关积累存储模块用于存储捕获结果Pmax以及其左右相邻各一搜索通道内的相关结果和并存储三个搜索通道对应的载波多普勒频率值fmax,fmid和fmin,同时存储捕获结果Pmax左右相邻各一码相位处的相关结果和并存储两个搜索位置对应的码相位值τmid和τmin;
本发明还提供一种直扩信号时频二维内插精估方法,包括以下步骤:
步骤一、接收机天线接收信号,再通过射频前端进行下变频及采样;
其中,下变频通过下变频单元实现,输出为带有数据调制的复信号;采样通过AD采样单元实现,采样率为fs,采样结果为带有数据调制的基带采样序列,可以用如下公式(1)表示:
其中,N为总体采样点数;ts=1/fs为时域采样间隔,D(nts)代表nts采样时刻的调制数据,为二进制数据,n代表第n个采样点;C(nts-τ0)表示接收信号的伪码,τ0代表伪码初始相位,表示载波多普勒频率,其中,c为光速,v为本装置相对于卫星的运动速度,fRF为射频载波频率,ξ=v/c,代表本装置相对于卫星的运动速度与光速的比值;exp代表以e为底的幂级数;j代表虚数单位;为射频载波初始相位;
步骤二、伪码相关模块以和AD采样单元的相同采样率对存储器模块产生的一个周期伪码序列进行采样,再输入FFT模块进行处理后存入存储器模块,具体为:
步骤2.1伪码相关模块以和AD采样单元的相同采样率对存储器模块产生的一个周期伪码序列,即对本地伪码进行采样,得到采样输出;
其中,表示带有估计相位的本地伪码;为本地伪码的估计相位,Nc=Tcfs是本地伪码一个伪码周期的采样点数,Tc=LTcp为本地伪码的伪码周期,Tcp=1/Rcp为一个本地伪码的码片持续时间,Rcp为本地伪码的码速率;
步骤2.2FFT模块对步骤2.1的采样输出进行FFT变换并取共轭;
步骤2.3将步骤2.2所得的结果存入存储器模块;
步骤三、并行频率搜索模块确定搜索频率,并进行载波多普勒剥离,具体为:
步骤3.1并行频率搜索模块根据信道先验信息估计载波多普勒频偏范围;
其中,所述的信道先验信息为可预知的最大载波多普勒频偏;所述的载波多普勒频偏范围,记为:(-fd,max,fd,max),fd,max即最大载波多普勒频偏;
步骤3.2将频偏范围按固定搜索步进划分成多个频率单元;
其中,固定搜索步进,记为Δfd,划分成L个频率单元,且第l个本地频率单元对应的搜索频率可以用如下公式(4)表示:
fl=-fd,max+lΔfd,l=0,1,2,...,L-1 (4)
步骤3.3并行频率搜索模块将本地载波与接收信号相乘完成载波多普勒剥离,具体为:
所述的载波多普勒剥离可以用如下公式(5)表示:
其中,Pl表示非相干积累模块中的累加次数;
步骤四、伪码相关模块计算步骤3.3中载波多普勒频率剥离后的信号与本地伪码的循环相关结果,并取模值,具体为:
步骤4.1先将第l条搜索通道的接收信号rl(n)以每行采样点长度为一个伪码周期Nc,分为Pl行,令i=0,1,...,Pl-1,经过划分后的rl(n)可写为:
步骤4.2伪码相关模块对步骤4.1划分后的信号按行做FFT,得到每行的FFT结果如下公式(7):
rl′(pl,k)=FFT[rl(pl,n)],k=0,1,2,...,Nc-1 (7)
步骤4.3伪码相关模块取步骤4.2的结果与步骤2.2的结果进行逐点相乘,并进行IFFT变换,得到第l条频率搜索通道的Pl次相关结果yl(pl,n)如下公式(8):
步骤五、非相干积累模块将步骤4.3所得的每个频率搜索通道内的Pl个相关结果进行非相干累加,得到长度为Nc的相关序列,如下公式(11):
步骤六、峰值检测模块对步骤五所得所有频率搜索通道内的相关积累结果进行幅值大小比较,设第l0条频率搜索通道内第τ0个位置出现最大幅值,则选取其结果定义为Pmax;
步骤七、相关积累存储模块存储峰值检测模块所检测得到的幅值最大结果Pmax,及其左右相邻各一频率通道内的相关积累结果,具体为:
步骤7.1相关积累存储模块存储峰值检测模块所检测得到的幅值最大结果Pmax,包括其对应的频率值fmax和相位值τmax;
步骤7.2相关积累存储模块同时存储与Pmax左右相邻各一频率通道l0-1和l0+1内,相同相位τmax处的相关积累结果,并比较三个相关积累结果的幅值大小,三个相关积累结果幅值可表示为:
步骤7.3相关积累存储模块同时存储与Pmax左右相邻各一码相位内,相同频率fmax处的相关积累结果,并比较三个相关积累结果的幅值大小,三个相关积累结果幅值可表示为:
其中,|τmax-τmid|=|τmid-τmin|=Δτ。
步骤8.6求取步骤8.4和步骤8.5中两条斜边的交点,即为三角包络的顶点,记作Pv τ,利用三角包络与伪码相关包络顶点的位置具有相似性可知,Pv τ所对应的伪码相位值即为精化后的伪码相位值;
至此,从步骤一到步骤八完成了一种直扩信号时频二维内插精估方法。
(三)有益效果
上述技术方案所提供的直扩信号时频二维内插精估方法,具有如下有益效果:
1.本发明建立在传统码域并行搜索捕获结构的基础上,实现简单,无需加入额外的复杂结构,保证了捕获速度,在捕获时间受限的情况下能够快速、准确地完成时频二维精估;
2.通过利用三角包络与sinc包络及伪码相关包络顶点位置的相似性,避免选取大量点数进行拟合,在低信噪比下可通过非相干积累有效获取信噪比增益,只需选取5个采样点便可完成时频二维精估。
附图说明
图1为本发明实施例1的示意图;
图2为本发明实施例1中载波多普勒频率精化结果仿真图;
图3为本发明实施例1中伪码相位精化结果仿真图。
具体实施方式
为使本发明的目的、内容和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
为了提高低信噪比、高动态下的捕获精度,本发明研究了一种基于三角包络拟合的时频二维精估算法,利用时频二维模糊函数(cross-ambiguity function,CAF)切面与三角包络的相似性,将精估过程简化为求取三角包络顶点位置的过程,在不提升捕获模块采样率和采样点数的情况下,有效改善低信噪比、高动态DSSS信号的估计精度,同时捕获时间没有明显增加。
实施例1
本实施例阐述了将本发明“一种直扩信号时频二维内插精估方法”应用于大频偏、低信噪比下直接序列扩频信号时频二维精估的结构框图。本实施例以射频频率为fRF=2.2GHz,伪码周期为1023,码片速率为Rcp=3.069Mcps,数据速率为RD=2kbps,取相关时间长度为1个伪码周期,载波多普勒频率为fd∈(-200kHz,200kHz)为例,由图1可以看出,接收机天线接收信号,再通过射频前端下变频和采样,再进行频率并行搜索,其结果输出给FFT模块,FFT模块处理后,输出值伪码相关模块并与存储器模块中存储的本地伪码序列FFT结果进行相关运算,再经非相干积累模块进行信噪比积累,之后将相关积累结果输出给峰值检测模块完成检测,利用相关积累存储模块对5个相关积累结果进行存储,最终经过三角包络拟合模块处理得出时频二维精估结果。
实施例2
本实施例按照实施例1所述的参数具体阐述了执行本发明步骤一到八后时频二维精估结果,频域结果如图2所示,时域结果如图3所示;
图2中,横坐标表示捕获之后的频差,其单位为Hz;纵坐标表示频率精化之后的频差,单位为Hz;
从图2中可以看出,本实施例中,捕获之后最大频差可达1kHz,经过本方法对载波多普勒频率进行精化后,残余频差最大不超过90Hz,极大地优化了捕获之后的频率误差,更有利于后续载波跟踪环路快速稳定地锁定信号。
图3中,横坐标表示捕获之后的码相位差,其单位为chip;纵坐标表示时域精化之后的码相位精化比例,定义码相位精化比例为经过精化处理后,码相位估计误差与处理前之比;
从图3中可以看出,本实施例中,精化算法在捕获后码相位误差越大的情况下,具有越为明显的精化效果,当码相位搜索步进为1码片时,捕获之后所得码相位误差相较于其他两组数值较大,但通过精化之后,码相位误差明显减小,完全满足后续伪码跟踪环路入锁需求。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (10)
1.一种直扩信号时频二维内插精估装置,其特征在于,包括接收机天线、射频前端、存储器模块、并行频率搜索模块、FFT模块、伪码相关模块、非相干积累模块、峰值检测模块、相关积累存储模块以及三角包络拟合模块;射频前端包括下变频单元以及AD采样单元;伪码相关模块包括预存的本地伪码快速傅立叶变换、乘法单元以及逆快速傅立叶变换单元;并行频率搜索模块包含L个频率通道;
接收机天线与射频前端相连,射频前端与并行频率搜索模块相连,并行频率搜索模块与FFT模块相连,FFT模块和存储器模块均与伪码相关模块相连,伪码相关模块与非相干积累模块相连,非相干积累模块与峰值检测模块相连,峰值检测模块与相关积累存储模块相连,相关积累存储模块与三角包络拟合模块相连;
接收机天线用于接收信号;射频前端用于对接收信号进行下变频和采样;存储器模块用于存储经过FFT变换后一个周期的本地伪码序列;并行搜索模块用于按固定频率步进产生本地载波,再与接收信号相乘,以剥离多普勒频率,再将此剥离多普勒频率后的信号输出FFT模块;FFT模块用于对并行搜索模块输出的每一组序列进行快速傅立叶变换,得到伪码相关模块的输入FFT序列;伪码相关模块用于将FFT模块的输出序列与本地伪码信号的快速傅立叶变换序列相乘,再通过逆快速傅立叶变换得到相关结果,取模值,其中乘法单元调用预存的本地伪码FFT序列,与折叠相加后FFT模块的输出依次相乘,IFFT单元对乘法单元的输出进行逆快速傅立叶变换,取模单元对IFFT单元的输出进行取模运算;非相干积累模块用于将伪码相关模块输出的结果按照相应位置进行非相干累加,提高输出信噪比;峰值检测模块用于比较非相干积累模块输出的相关结果,找出其中的最大值,根据其所处的搜索通道位置得到载波多普勒频率估计值fmax,并由峰值检测模块输出的峰值位置得到码相位估计值τmax,作为捕获结果Pmax;相关积累存储模块用于存储捕获结果Pmax以及其左右相邻各一搜索通道内的相关结果和并存储三个搜索通道对应的载波多普勒频率值fmax,fmid和fmin,同时存储捕获结果Pmax左右相邻各一码相位处的相关结果和并存储两个搜索位置对应的码相位值τmid和τmin;三角包络拟合模块用于计算相关积累存储模块中三个相关结果Pmax,和以及Pmax,和所构成的三角包络,将三角包络顶点和所对应的频率值fv和相位值τv作为载波多普勒频率和伪码相位的精估结果。
2.一种直扩信号时频二维内插精估方法,其特征在于,包括以下步骤:
步骤一、接收机天线接收信号,再通过射频前端进行下变频及采样;
步骤二、伪码相关模块以和AD采样单元的相同采样率对存储器模块产生的一个周期伪码序列进行采样,再输入FFT模块进行处理后存入存储器模块;
步骤三、并行频率搜索模块确定搜索频率,并进行载波多普勒剥离;
步骤四、伪码相关模块计算步骤3.3中载波多普勒频率剥离后的信号与本地伪码的循环相关结果,并取模值;
步骤五、非相干积累模块将步骤4.3所得的每个频率搜索通道内的Pl个相关结果进行非相干累加,得到长度为Nc的相关序列;
步骤六、峰值检测模块对步骤五所得所有频率搜索通道内的相关积累结果进行幅值大小比较,设第l0条频率搜索通道内第τ0个位置出现最大幅值,则选取其结果定义为Pmax;
3.如权利要求2所述的直扩信号时频二维内插精估方法,其特征在于,所述步骤一中,下变频通过下变频单元实现,输出为带有数据调制的复信号;采样通过AD采样单元实现,采样率为fs,采样结果为带有数据调制的基带采样序列,用如下公式(1)表示:
4.如权利要求3所述的直扩信号时频二维内插精估方法,其特征在于,所述步骤二包括以下子步骤:
步骤2.1伪码相关模块以和AD采样单元的相同采样率对存储器模块产生的一个周期伪码序列,即对本地伪码进行采样,得到采样输出;
步骤2.2FFT模块对步骤2.1的采样输出进行FFT变换并取共轭;
步骤2.3将步骤2.2所得的结果存入存储器模块。
6.如权利要求5所述的直扩信号时频二维内插精估方法,其特征在于,所述步骤三包括以下子步骤:
步骤3.1并行频率搜索模块根据信道先验信息估计载波多普勒频偏范围;
步骤3.2将频偏范围按固定搜索步进划分成多个频率单元;
步骤3.3并行频率搜索模块将本地载波与接收信号相乘完成载波多普勒剥离。
8.如权利要求7所述的直扩信号时频二维内插精估方法,其特征在于,所述步骤四包括以下子步骤:
步骤4.1先将第l条搜索通道的接收信号rl(n)以每行采样点长度为一个伪码周期Nc,分为Pl行,令i=0,1,...,Pl-1,经过划分后的rl(n)写为:
步骤4.2伪码相关模块对步骤4.1划分后的信号按行做FFT,得到每行的FFT结果如下公式(7):
rl′(pl,k)=FFT[rl(pl,n)],k=0,1,2,...,Nc-1 (7)
步骤4.3伪码相关模块取步骤4.2的结果与步骤2.2的结果进行逐点相乘,并进行IFFT变换,得到第l条频率搜索通道的Pl次相关结果yl(pl,n)如下公式(8):
步骤五中,长度为Nc的相关序列如下公式(11):
9.如权利要求8所述的直扩信号时频二维内插精估方法,其特征在于,所述步骤七包括以下子步骤:
步骤7.1相关积累存储模块存储峰值检测模块所检测得到的幅值最大结果Pmax,包括其对应的频率值fmax和相位值τmax;
步骤7.2相关积累存储模块同时存储与Pmax左右相邻各一频率通道l0-1和l0+1内,相同相位τmax处的相关积累结果,并比较三个相关积累结果的幅值大小,三个相关积累结果幅值表示为:
步骤7.3相关积累存储模块同时存储与Pmax左右相邻各一码相位内,相同频率fmax处的相关积累结果,并比较三个相关积累结果的幅值大小,三个相关积累结果幅值表示为:
其中,|τmax-τmid|=|τmid-τmin|=Δτ。
10.如权利要求9所述的直扩信号时频二维内插精估方法,其特征在于,所述步骤八包括以下子步骤:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111122632.1A CN113890563A (zh) | 2021-09-24 | 2021-09-24 | 一种直扩信号时频二维内插精估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111122632.1A CN113890563A (zh) | 2021-09-24 | 2021-09-24 | 一种直扩信号时频二维内插精估方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113890563A true CN113890563A (zh) | 2022-01-04 |
Family
ID=79006607
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111122632.1A Pending CN113890563A (zh) | 2021-09-24 | 2021-09-24 | 一种直扩信号时频二维内插精估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113890563A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115412125A (zh) * | 2022-08-26 | 2022-11-29 | 天津津航计算技术研究所 | 一种大频偏直扩信号一阶差分位同步方法及系统 |
CN115664914A (zh) * | 2022-10-27 | 2023-01-31 | 天津津航计算技术研究所 | 一种基于动态平均的高精度位同步方法 |
CN118174748A (zh) * | 2024-05-10 | 2024-06-11 | 中国人民解放军战略支援部队航天工程大学 | 一种低复杂度的非搜索一阶动态直扩信号捕获方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090037147A1 (en) * | 2007-08-03 | 2009-02-05 | Oracle International Corporation | Fast intrinsic mode decomposition of time series data with sawtooth transform |
CN103543456A (zh) * | 2013-10-24 | 2014-01-29 | 北京大学 | 一种基于分段相关结合fft运算的大频偏gnss信号捕获方法 |
CN109412644A (zh) * | 2018-09-13 | 2019-03-01 | 南京理工大学 | 一种直扩msk信号的多普勒频率估计方法 |
CN110943758A (zh) * | 2019-11-12 | 2020-03-31 | 天津津航计算技术研究所 | 基于多普勒频率变化率和频率联合搜索的二次捕获方法 |
CN110958035A (zh) * | 2019-11-12 | 2020-04-03 | 天津津航计算技术研究所 | 一种基于三角包络拟合的频率二次捕获装置及方法 |
-
2021
- 2021-09-24 CN CN202111122632.1A patent/CN113890563A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090037147A1 (en) * | 2007-08-03 | 2009-02-05 | Oracle International Corporation | Fast intrinsic mode decomposition of time series data with sawtooth transform |
CN103543456A (zh) * | 2013-10-24 | 2014-01-29 | 北京大学 | 一种基于分段相关结合fft运算的大频偏gnss信号捕获方法 |
CN109412644A (zh) * | 2018-09-13 | 2019-03-01 | 南京理工大学 | 一种直扩msk信号的多普勒频率估计方法 |
CN110943758A (zh) * | 2019-11-12 | 2020-03-31 | 天津津航计算技术研究所 | 基于多普勒频率变化率和频率联合搜索的二次捕获方法 |
CN110958035A (zh) * | 2019-11-12 | 2020-04-03 | 天津津航计算技术研究所 | 一种基于三角包络拟合的频率二次捕获装置及方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115412125A (zh) * | 2022-08-26 | 2022-11-29 | 天津津航计算技术研究所 | 一种大频偏直扩信号一阶差分位同步方法及系统 |
CN115412125B (zh) * | 2022-08-26 | 2023-07-21 | 天津津航计算技术研究所 | 一种大频偏直扩信号一阶差分位同步方法及系统 |
CN115664914A (zh) * | 2022-10-27 | 2023-01-31 | 天津津航计算技术研究所 | 一种基于动态平均的高精度位同步方法 |
CN118174748A (zh) * | 2024-05-10 | 2024-06-11 | 中国人民解放军战略支援部队航天工程大学 | 一种低复杂度的非搜索一阶动态直扩信号捕获方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113890563A (zh) | 一种直扩信号时频二维内插精估方法 | |
CN106230475B (zh) | 一种基于Tong检测器伪码多普勒补偿捕获方法及装置 | |
CN107135013B (zh) | 一种直接序列扩频信号快速捕获方法 | |
CN106646546B (zh) | 一种卫星信号多维快速捕获方法及系统 | |
CN103645483B (zh) | 一种弱信号环境下北斗信号捕获方法 | |
CN108401581B (zh) | 一种卫星扩频通信系统中的pn码快速捕获方法 | |
CN109655847B (zh) | 一种适于动态信号的快速捕获方法 | |
CN110943758A (zh) | 基于多普勒频率变化率和频率联合搜索的二次捕获方法 | |
CN105842713A (zh) | 一种基于sft的ins辅助bds信号快速捕获方法 | |
CN110501729A (zh) | 一种基于fpga分步码相位细化的gnss信号的捕获方法 | |
CN106526631B (zh) | 北斗b1频点卫星信号载波频率高精度估计方法 | |
CN103941269A (zh) | 用于卫星导航系统的pn码捕获方法 | |
CN110958035A (zh) | 一种基于三角包络拟合的频率二次捕获装置及方法 | |
CN110456393A (zh) | 北斗弱信号快速捕获方法 | |
CN103926604A (zh) | 基于重叠差分循环相干积分的弱信号捕获方法 | |
CN104035109A (zh) | 基于重叠1/5比特差分循环相干积分的弱信号捕获方法 | |
CN114114335A (zh) | 一种弱信号gnss接收机快速捕获方法 | |
Ahamed et al. | Fast acquisition of GPS signal using Radix-2 and Radix-4 FFT algorithms | |
Ahamed et al. | Fast acquisition of GPS signal using FFT decomposition | |
CN110018500B (zh) | 一种基于圆周移位的北斗卫星信号捕获方法 | |
CN103969664A (zh) | 一种基于差分循环相干积分的gps弱信号捕获方法 | |
CN115865127B (zh) | 一种直扩信号的参数估计和解调方法 | |
CN109633707B (zh) | 一种基于预平均处理的变系数匹配滤波的快速捕获方法 | |
CN101206259A (zh) | 基于数字匹配滤波器的多星座导航弱信号的捕获方法 | |
CN116338740A (zh) | 一种基于pcf重构的北斗b1c信号无模糊捕获计算方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20220104 |
|
RJ01 | Rejection of invention patent application after publication |