CN102680988B - 高灵敏度导航卫星信号非线性捕获方法及装置 - Google Patents

高灵敏度导航卫星信号非线性捕获方法及装置 Download PDF

Info

Publication number
CN102680988B
CN102680988B CN 201210162393 CN201210162393A CN102680988B CN 102680988 B CN102680988 B CN 102680988B CN 201210162393 CN201210162393 CN 201210162393 CN 201210162393 A CN201210162393 A CN 201210162393A CN 102680988 B CN102680988 B CN 102680988B
Authority
CN
China
Prior art keywords
signal
satellite
maximal value
module
sampling
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.)
Expired - Fee Related
Application number
CN 201210162393
Other languages
English (en)
Other versions
CN102680988A (zh
Inventor
莫建文
张顺岚
张彤
首照宇
袁华
欧阳宁
陈利霞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guilin University of Electronic Technology
Original Assignee
Guilin University of Electronic Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guilin University of Electronic Technology filed Critical Guilin University of Electronic Technology
Priority to CN 201210162393 priority Critical patent/CN102680988B/zh
Publication of CN102680988A publication Critical patent/CN102680988A/zh
Application granted granted Critical
Publication of CN102680988B publication Critical patent/CN102680988B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明所要解决的技术问题是提供一种高灵敏度导航卫星信号非线性捕获方法及装置,该方法与装置采用DBZP技术来减小相关功率损失;用FFT技术减少了做相关所需时间;利用非线性双稳随机共振***的特性来提高***输出信噪比;捕获过程中引入了MTM(最大相关值与相关均值比值)阈值检测方法,提高了正确检测率。本发明能够最大限度地提高GPS信号捕获灵敏度;并且可以极大限度地减小捕获时间,为弱GPS信号的实时高灵敏度捕获提供技术支持。

Description

高灵敏度导航卫星信号非线性捕获方法及装置
技术领域
本发明适用于弱导航卫星信号捕获领域,具体涉及一种高灵敏度导航卫星信号非线性捕获方法及装置。
背景技术
近年来,全球定位***(Global Positioning System,GPS)已经得到越来越广泛的应用。随着这项定位技术发展的不断深入,用户对其使用要求也越来越高,无盲区、高灵敏度、实时定位与导航成为该技术发展的方向。然而,绝大多数使用手持定位设备的用户分布在高楼密集的城市地区,且大部分时间处于室内。由于受到的遮挡、多径和干扰等现象较严重,与普通环境相比,GPS信号在室内环境中,能量有更多的削弱和衰落,到达时间有更大的延迟,接收信噪比有更大程度的恶化,所以,其可用性和定位精度都会大大下降,GPS的总体性能严重恶化。在室内等低信噪比环境下,对GPS接收机最先启动的工作部分,即捕获功能模块的性能指标提出了更高的要求,必须能够快速高灵敏度捕获低信噪比环境下的信号。但是,目前的GPS捕获技术几乎很难满足这一要求,现有的GPS接收机在室内等环境下很难工作。常见的相干积分、非相干积分、差分相干积分等捕获方法能在一定程度上提高灵敏度,但均要求积累时间较长,无法满足室内弱信号高灵敏度及实时捕获的要求。如果GPS信号不能被快速高灵敏度捕获到,那么后续的所有功能都将受影响,比如跟踪、导航电文提取、位置解算等。
2011年04月《***工程与电子技术》第33卷第4期中提出了一种“基于随机共振的高灵敏度GPS信号捕获算法”,该算法首先用部分匹配滤波器对GPS信号进行分段相关预处理,然后利用随机共振提高预处理后信号的信噪比,实现在相对较短的相关累积时间内获得较高的捕获灵敏度。然而该算法在GPS信号多普勒频移较大且积累时间较长时,将造成码片速率发生较大变化,从而对码周期产生较大影响,当对数字中频信号与本地信号做相关时,会造成相关峰值衰减较大,从而影响捕获灵敏度;并且,上述算法中的部分匹配滤波是一小段一小段做相关,造成预处理过程耗时较大,难以达到快速高灵敏度GPS信号捕获的要求。
发明内容
本发明所要解决的技术问题是提供一种高灵敏度导航卫星信号非线性捕获方法及装置,该方法与装置能够有效减小大多普勒频率造成相关过程中相关功率的损失;利用非线性***-双稳随机共振***的特性,即在一定条件下,噪声能量会向信号能量转移的特性,最大限度地提高GPS信号捕获灵敏度;并且可以极大限度地减小捕获时间,为弱GPS信号的实时高灵敏度捕获提供技术支持。
为解决上述问题,本发明是通过以下方案实现的:
本发明一种高灵敏度导航卫星信号非线性捕获方法,包括如下步骤:
步骤一:对导航卫星发射的卫星模拟信号进行下变频及采样,并设定合适的多普勒频移变量初值及最大值,即设定多普勒频移搜索范围,同时设定多普勒频移搜索间隔、MTM(最大相关值与相关均值的比值)的检测阈值、双稳随机共振***参数、卫星号变量初值及卫星号最大值;
步骤二:根据待搜索的导航卫星号及多普勒频率偏移产生本地伪码卫星信号;
步骤三:分别对采样后的卫星信号和本地伪码卫星信号进行双块零拓展;
步骤四:对经双块零拓展后的卫星信号和本地伪码卫星信号利用快速傅里叶变换循环卷积进行相关运算;
步骤五:对相关结果进行处理,并保存相关结果的实部信号;
步骤六:对相关结果的实部信号进行二次采样随机共振;
步骤七:对步骤六中输出的信号做傅里叶变换并取模;
步骤八:对每一颗卫星,在频率方向和码延时方向上寻找步骤七中取模后的最大值;即对每一颗卫星,首先寻找初始频率对应的步骤七中取模后的最大值并保存;然后返回步骤二并找出下一多普勒频率对应的步骤七中取模后的最大值并保存;直到保存该颗卫星的所有频率对应的步骤七中取模后的最大值;最后再从上述所有频率对应的取模后的最大值中找出一个最大值作为该颗卫星取模后的最大值;
之后,将该颗卫星取模后的最大值同步骤一中设定的检测阈值进行比较,如果该最大值大于步骤一中设定的检测阈值,则该最大值所对应的频率和码延时即为要捕获的多普勒频率和码延时;如果该最大值未超过步骤一中设定的检测阈值,则重复步骤二搜索下一颗卫星,直到搜索完所有的卫星,则退出捕获,此时表示捕获不成功。
上述步骤六所述二次采样随机共振步骤包括二次采样和随机共振步骤,即首先对保存的相关结果的实部信号进行二次采样,使采样后信号的频率符合小参数随机共振的要求;然后让采样后的信号进入双稳随机共振***,并进行共振输出。
上述多普勒频移搜索范围最好为-10KHz~10KHz。
上述多普勒频移搜索间隔最好为5000Hz,这样只需搜索4次即可搜完1颗GPS卫星的-10KHz~10KHz多普勒频移范围,从而大大缩短了搜索时间,提高了捕获效率。
上述检测阈值最好介于3~7之间。
上述双稳随机共振***参数为a=1,b=1。
上述卫星号变量的初值为1,最大值为30。
本发明一种高灵敏度导航卫星信号非线性捕获装置,包括如下模块:
变频采样模块:对导航卫星发射的卫星模拟信号进行下变频及采样;
本地伪码生成器:根据待搜索的导航卫星号及多普勒频率偏移产生本地伪码卫星信号;
双块零拓展模块:分别对变频采样模块采样后的卫星信号和本地伪码生成器产生的本地伪码卫星信号进行双块零拓展;
快速傅里叶变换循环卷积相关模块:对经双块零拓展模块输出的卫星信号和本地伪码卫星信号利用快速傅里叶变换循环卷积进行相关运算;
取实部模块:对相关结果进行处理,并保存相关结果的实部信号;
二次采样随机共振模块:对取实部模块保存的相关结果的实部信号进行二次采样,使其通过双稳随机共振***,并进行共振输出;
傅里叶变换及取模模块:对二次采样随机共振模块输出的信号做傅里叶变换并取模;
MTM(最大相关值与相关均值的比值)阈值检测模块:对每颗卫星,在频率方向和码延时方向上寻找傅里叶变换及取模模块输出的最大值,即对每一颗卫星,首先寻找初始频率对应的取模后的最大值并保存;然后重复返回本地伪码生成器找出下一多普勒频率对应的取模后的最大值并保存,直到保存该颗卫星的所有频率对应的取模后的最大值;最后再从上述所有频率对应的取模后的最大值中找出一个最大值作为该颗卫星取模后的最大值;
如果该颗卫星取模后的最大值大于预先设定的检测阈值,则该最大值所对应的频率和码延时即为要捕获的多普勒频率和码延时;如果该最大值未超过预先设定的检测阈值,则返回本地伪码生成器搜索下一颗卫星,直到搜索完所有的卫星,则退出捕获,此时,表示捕获不成功。
上述二次采样随机共振模块最好包括二次采样单元和随机共振单元,其中,二次采样单元对取实部模块保存的相关结果的实部信号进行二次采样,使二次采样后信号的频率符合双稳随机共振的小参数信号要求;随机共振单元对二次采样后的信号进行共振输出。
与现有技术相比,本发明在尽量短的时间内,最大限度地提高了GPS信号捕获灵敏度。该发明采用DBZP技术来减小相关功率损失;用FFT技术减少了做相关所需时间;利用非线性双稳随机共振***的特性——强噪声干扰下的信号作用于双稳随机共振非线性***时,信号和噪声在非线性***的协同作用下,会发生噪声能量向信号能量的转移,产生类似于力学中的共振输出,从而极大地提高***输出信噪比;捕获过程中引入了MTM(最大相关值与相关均值比值)阈值检测方法,提高了正确检测率。
附图说明
图1为一种高灵敏度导航卫星信号捕获***原理图。
具体实施方式
一种高灵敏度导航卫星信号捕获方法,如图1所示,包括如下步骤:
步骤一:对导航卫星发射的卫星模拟信号进行下变频及采样,并设定合适的多普勒频移变量初值及最大值,即设定多普勒频移搜索范围,同时设定多普勒频移搜索间隔、MTM(最大相关值与相关均值的比值)的检测阈值、双稳随机共振***参数、卫星号变量初值及卫星号最大值。
在本发明中,所述多普勒频移搜索范围为-10KHz~10KHz。所述多普勒频移的搜索间隔为5000Hz。所述检测阈值是指MTM的检测阈值,其中MTM表示最大相关值与相关均值的比值,若计算所得MTM值超过检测阈值BT,即认为信号被准确捕获,在本发明中,所述检测阈值BT设定为3~7之间。所述双稳随机共振***参数a=1,b=1。根据目前GPS卫星的数目,在本发明中,所述卫星号变量的初值设定为1,最大值设定为30。
从卫星发射的L1频段GPS模拟信号下变频和采样后的数字中频信号为
r ( t j ) = Ad ( t j ) C ( ( t j - τ ) ( 1 + f d f L 1 ) ) cos ( 2 π · ( f IF + f d ) t j + φ ) + n ( t j ) - - - ( 1 )
式中,A为信号幅度,d(t)为导航数据信息,C(t)为C/A码,fd为多普勒频移,fL1为高频载波频率,fIF为中频载波频率,φ为初始相位,n(t)为加性高斯白噪声,tj为第j个采样时刻,τ为信号到达接收机的时延。
步骤二:根据待搜索的导航卫星号及多普勒频率偏移产生本地伪码卫星信号。
第k个时间段本地产生的伪码信号模型为
s k ( t i ) = C ( ( t i - τ ^ ) ( 1 + f ^ d f L 1 ) ) exp ( j 2 π ( f IF + f ^ d ) t i ) - - - ( 2 )
式中,sk(i)第k个时间段本地产生的伪码信号,C(t)为C/A码,
Figure GDA00003040545200043
为多普勒频移估计值,fL1为高频载波频率,fIF为中频载波频率,ti=(k·N+i)Ts为i个采样时刻,
Figure GDA00003040545200044
为时延的估计值。
步骤三:分别对采样后的卫星信号和本地伪码卫星信号进行双块零拓展。
由于多普勒频移对码相位的影响因子为
Figure GDA00003040545200045
当多普勒频移较大且积累时间较长时,造成码片速率发生较大变化,从而对码周期产生较大影响,直接对数字中频信号与本地信号做相关并进行差分累加,会造成相关峰值衰减较大,从而影响捕获灵敏度,DBZP技术可以改善相关过程中由大多普勒频移引起码片速率变化造成的相关功率损失。
拓展后的中频采样GPS信号和本地产生的伪码GPS信号分别表示为
r k &prime; ( i ) = r k ( i ) 0 &le; i < N r k + 1 ( i - N ) N &le; i < 2 N - - - ( 3 )
s k &prime; ( i ) = s k ( i ) 0 &le; i < N 0 N &le; i < 2 N - - - ( 4 )
式中,rk(i)为第k个相干时间段的数字中频信号;sk(i)第k个时间段本地产生的伪码信号。
步骤四:对经双块零拓展后的卫星信号和本地伪码卫星信号利用快速傅里叶变换循环卷积进行相关运算。
步骤五:对相关结果进行处理,并保存相关结果的实部信号。
将式(4)、(5)中的对应块利用FFT循环卷积进行相关运算,有
g ( k , &tau; ^ ) = Re [ IFFT ( FFT ( r k &prime; ( i ) ) FFT ( s k &prime; ( i ) ) ) * ) ]
= 1 2 Ad ( k ) R c ( &Delta;&tau; ) T a S a ( &pi;&Delta; f d T a ) cos ( ( 2 k - 1 ) &pi;&Delta; f d T a + &theta; ) + n ( k ) - - - ( 6 )
式中
Figure GDA00003040545200054
Figure GDA00003040545200055
下的不同码时延的相关结果,Rc(△τ)为C/A码循环卷积相关值,
Figure GDA00003040545200056
Figure GDA00003040545200057
为残余多普勒频率偏差,Sa(x)=sin(x)/x,n(k)为积分后的噪声项,其方差为
Figure GDA000030405452000511
改写式(6),有
g ( t k , &tau; ^ ) = A &prime; cos ( 2 &pi; ft k + &theta; ) + n ( k ) - - - ( 7 )
式中,
Figure GDA00003040545200059
为处理后信号幅度,f=△fd为要搜索的残余多普勒频率偏差,为第k段时间序列。
步骤六:对相关结果的实部信号进行二次采样随机共振。在本发明中,所述二次采样随机共振步骤包括二次采样和随机共振步骤,即首先,对保存的相关结果的实部信号进行二次采样,使采样后的信号频率符合小参数随机共振的要求,然后,让采样后的信号进入双稳随机共振***,并进行共振输出。
连续的双稳态随机共振一维动力学***的Langevin方程模型为
dx/dt=ax-bx3+Acos(2πf0t)+n(t)    (8)
式中A为被检测弱周期信号幅度,f0为信号频率,n(t)是强度为D、均值为0、方差为1的白噪声,a、b为***参数,x(t)为***输出信号,双稳***的势垒高度为△U=a2/4b。
由于双稳随机共振***要求的驱动频率f0很低,被限制在小参数(频率、幅度、噪声强度均远小于1)范围内。而式(7)中信号的频率f=△fd较大,一般为[-10KHz-+10KHz],远远超出小参数范围。
采用二次采样随机共振理论将式(7)中的f转换到小参数范围内,设输入信号频率为f,采样频率为fs,二次采样频率为fsr,二次采样后信号频率为fr,定义fr/fsr=f′0/fs,则fr=(f/fs)fsr,其中fr的数量级应控制在绝热近似理论的小参数范围内。二次采样后,(7)式可写为
g ( t k , &tau; ^ ) = A &prime; cos ( 2 &pi; f r t k + &theta; ) + n ( k ) - - - ( 9 )
以式(9)中
Figure GDA00003040545200062
为输入信号进入式(8)所表示的双稳态随机共振***,有
dx ( t k , &tau; ^ ) / dt = ax ( t k , &tau; ^ ) - b x 3 ( t k , &tau; ^ ) + g ( t k , &tau; ^ ) - - - ( 10 )
并利用四阶龙格-库塔方法求解方程(10),即
x ( t k + 1 , &tau; ^ ) = x ( t i , &tau; ^ ) + ( k 1 + 2 k 2 + 2 k 3 + k 4 ) / 6 - - - ( 11 )
k 1 = h &lsqb; ax ( t k , &tau; ^ ) - b x 3 ( t k , &tau; ^ ) + g ( t k , &tau; ^ ) &rsqb; - - - ( 12 )
k 2 = h { a [ x ( t k , &tau; ^ ) + k 1 / 2 ] - b [ x ( t k , &tau; ^ ) + k 1 / 2 ] 3 + g ( t k , &tau; ^ ) } - - - ( 13 )
k 3 = h { a [ x ( t k , &tau; ^ ) + k 2 / 2 ] - b [ x ( t k , &tau; ^ ) + k 2 / 2 ] 3 + g ( t k + 1 , &tau; ^ ) } - - - ( 14 )
k 4 = h { a [ x ( t k , &tau; ^ ) + k 3 ] - b [ x ( t k , &tau; ^ ) + k 3 ] 3 + g ( t k + 1 , &tau; ^ ) } - - - ( 15 )
其中h=1/fsr为时间步长。
步骤七:对步骤六中输出的信号做傅里叶变换并取模。
对双稳随机共振输出做傅里叶变换并取模,有
Y ( u , &tau; ^ ) = | FFT ( x ( t k , &tau; ^ ) ) | , u=0,1…N-1    (16)
式中,u为数字频率的序列号,
Figure GDA000030405452000610
为码延时。
步骤八:对每一颗卫星,在频率方向和码延时方向上寻找步骤七中取模后的最大值:即对每一颗卫星,首先寻找初始频率对应的步骤七中取模后的最大值并保存;然后返回步骤二并找出下一多普勒频率对应的步骤七中取模后的最大值并保存;直到保存该颗卫星的所有频率对应的步骤七中取模后的最大值;最后再从上述所有频率对应的取模后的最大值中找出一个最大值作为该颗卫星取模后的最大值。
之后,将该颗卫星取模后的最大值同步骤一中设定的检测阈值进行比较,如果该最大值大于步骤一中设定的检测阈值,则该最大值所对应的频率和码延时即为要捕获的多普勒频率和码延时;如果该最大值未超过步骤一中设定的检测阈值,则重复步骤二搜索下一颗卫星,直到搜索完所有的卫星,则退出捕获,此时表示捕获不成功。
上述过程即是,在u和
Figure GDA000030405452000611
方向上中寻找式(16)中
Figure GDA000030405452000612
的最大值
Figure GDA000030405452000613
如果该值大于步骤一中设定的检测阈值,则u及
Figure GDA000030405452000614
即为到要捕获的多普勒频率和码延时。
采用上述方法所设计的高灵敏度导航卫星信号非线性捕获装置,如图1所示,包括变频采样模块、本地伪码生成器、双块零拓展模块、快速傅里叶变换循环卷积相关模块、取实部模块、二次采样随机共振模块、傅里叶变换及取模模块和MTM(最大相关值与相关均值的比值)阈值检测模块。其中
变频采样模块:对导航卫星发射的卫星模拟信号进行下变频及采样。
本地伪码生成器:根据待搜索的导航卫星号及多普勒频率偏移产生本地伪码卫星信号。
双块零拓展模块:分别对变频采样模块采样后的卫星信号和本地伪码生成器产生的本地伪码卫星信号进行双块零拓展。
快速傅里叶变换循环卷积相关模块:对经双块零拓展模块输出的卫星信号和本地伪码卫星信号利用快速傅里叶变换循环卷积进行相关运算。
取实部模块:对相关结果进行处理,并保留相关结果的实部信号。
二次采样随机共振模块:对取实部模块保存的相关结果的实部信号进行二次采样,使其通过双稳随机共振***,并进行共振输出。本发明中所述二次采样随机共振模块包括二次采样单元和随机共振单元,其中二次采样单元对取实部模块保存的相关结果的实部信号进行二次采样,并使采样后信号频率符合双稳随机共振小参数信号的要求;随机共振单元对二次采样后的信号进行共振输出。
傅里叶变换及取模模块:对二次采样随机共振模块输出的信号做傅里叶变换并取模。
MTM(最大相关值与相关均值的比值)阈值检测模块:对每颗卫星,在频率方向和码延时方向上寻找傅里叶变换及取模模块输出的最大值,即对每一颗卫星,首先寻找初始频率对应的取模后的最大值并保存;然后重复返回本地伪码生成器找出下一多普勒频率对应的取模后的最大值并保存,直到保存该颗卫星的所有频率对应的取模后的最大值;最后再从上述所有频率对应的取模后的最大值中找出一个最大值作为该颗卫星取模后的最大值。如果该颗卫星取模后的最大值大于预先设定的检测阈值,则该最大值所对应的频率和码延时即为要捕获的多普勒频率和码延时;如果该最大值未超过预先设定的检测阈值,则返回本地伪码生成器搜索下一颗卫星,直到搜索完所有的卫星,则退出捕获,此时,表示捕获不成功。

Claims (9)

1.高灵敏度导航卫星信号非线性捕获方法,其特征是包括如下步骤:
步骤一:对导航卫星发射的卫星模拟信号进行下变频及采样,并设定合适的多普勒频移变量初值及最大值,即设定多普勒频移搜索范围,同时设定多普勒频移搜索间隔、最大相关值与相关均值的比值的检测阈值、双稳随机共振***参数、卫星号变量初值及卫星号最大值;
步骤二:根据待搜索的导航卫星号及多普勒频率偏移产生本地伪码卫星信号;
步骤三:分别对采样后的卫星信号和本地伪码卫星信号进行双块零拓展;
步骤四:对经双块零拓展后的卫星信号和本地伪码卫星信号利用快速傅里叶变换循环卷积进行相关运算;
步骤五:对相关结果进行处理,并保留相关结果的实部信号;
步骤六:对保存的相关结果的实部信号进行二次采样随机共振;
步骤七:对步骤六中输出的信号做傅里叶变换并取模;
步骤八:对每一颗卫星,在频率方向和码延时方向上寻找步骤七中取模后的最大值;即对每一颗卫星,首先寻找初始频率对应的步骤七中取模后的最大值并保存;然后返回步骤二并找出下一多普勒频率对应的步骤七中取模后的最大值并保存;直到保存该颗卫星的所有频率对应的步骤七中取模后的最大值;最后再从上述所有频率对应的取模后的最大值中找出一个最大值作为该颗卫星取模后的最大值;
之后,将该颗卫星取模后的最大值同步骤一中设定的检测阈值进行比较,如果该最大值大于步骤一中设定的检测阈值,则该最大值所对应的频率和码延时即为要捕获的多普勒频率和码延时;如果该最大值未超过步骤一中设定的检测阈值,则重复步骤二搜索下一颗卫星,直到搜索完所有的卫星,则退出捕获,此时表示捕获不成功。
2.根据权利要求1所述的高灵敏度导航卫星信号非线性捕获方法,其特征是:
步骤六所述二次采样随机共振步骤包括二次采样步骤和随机共振步骤,即首先,对保存的相关结果的实部信号进行二次采样,使采样后的信号频率符合小参数随机共振的要求,然后,让采样后的信号进入双稳随机共振***,并进行随机共振。
3.根据权利要求1或2所述高灵敏度导航卫星信号非线性捕获方法,其特征是:所述多普勒频移搜索范围为-10KHz~10KHz。
4.根据权利要求3所述高灵敏度导航卫星信号非线性捕获方法,其特征是:所述多普勒频移搜索间隔为5000Hz。
5.根据权利要求1或2所述高灵敏度导航卫星信号非线性捕获方法,其特征是,所述检测阈值介于3~7之间。
6.根据权利要求1或2所述高灵敏度导航卫星信号非线性捕获方法,其特征是:双稳随机共振***参数初值为a=1,b=1。
7.根据权利要求1或2所述高灵敏度导航卫星信号非线性捕获方法,其特征是:卫星号变量的初值为1,最大值为30。
8.高灵敏度导航卫星信号非线性捕获装置,其特征是包括如下模块:
变频采样模块:对导航卫星发射的卫星模拟信号进行下变频及采样;
本地伪码生成器:根据待搜索的导航卫星号及多普勒频率偏移产生本地伪码卫星信号;
双块零拓展模块:分别对变频采样模块采样后的卫星信号和本地伪码生成器产生的本地伪码卫星信号进行双块零拓展;
快速傅里叶变换循环卷积相关模块:对经双块零拓展模块输出的卫星信号和本地伪码卫星信号利用快速傅里叶变换循环卷积进行相关运算;
取实部模块:对相关结果进行处理,并保留相关结果的实部信号;
二次采样随机共振模块:对取实部模块保存的相关结果的实部信号进行二次采样,使其通过双稳随机共振***,并进行共振输出;
傅里叶变换及取模模块:对二次采样随机共振模块输出的信号做傅里叶变换并取模;
最大相关值与相关均值的比值阈值检测模块:对每颗卫星,在频率方向和码延时方向上寻找傅里叶变换及取模模块输出的最大值;如果该颗卫星取模后的最大值大于预先设定的检测阈值,则该最大值所对应的频率和码延时即为要捕获的多普勒频率和码延时;如果该最大值未超过预先设定的检测阈值,则返回本地伪码生成器搜索下一颗卫星,直到搜索完所有的卫星,则退出捕获,此时,表示捕获不成功。
9.根据权利要求8所述的高灵敏度导航卫星信号非线性捕获装置,其特征是:所述二次采样随机共振模块包括二次采样单元和随机共振单元,其中
二次采样单元,对取实部模块保存的相关结果的实部信号进行二次采样,使采样后信号频率符合双稳随机共振小参数信号的要求;
随机共振单元,对二次采样后的信号进行共振输出。
CN 201210162393 2012-05-23 2012-05-23 高灵敏度导航卫星信号非线性捕获方法及装置 Expired - Fee Related CN102680988B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210162393 CN102680988B (zh) 2012-05-23 2012-05-23 高灵敏度导航卫星信号非线性捕获方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210162393 CN102680988B (zh) 2012-05-23 2012-05-23 高灵敏度导航卫星信号非线性捕获方法及装置

Publications (2)

Publication Number Publication Date
CN102680988A CN102680988A (zh) 2012-09-19
CN102680988B true CN102680988B (zh) 2013-11-06

Family

ID=46813197

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210162393 Expired - Fee Related CN102680988B (zh) 2012-05-23 2012-05-23 高灵敏度导航卫星信号非线性捕获方法及装置

Country Status (1)

Country Link
CN (1) CN102680988B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103728632B (zh) * 2013-12-02 2017-10-03 西安合众思壮导航技术有限公司 卫星数据捕获装置和捕获方法
CN103809192A (zh) * 2014-02-25 2014-05-21 浙江理工大学 一种gnss接收机的动态校正算法
CN105093243A (zh) * 2014-05-08 2015-11-25 北京大学 一种基于随机共振算法的gnss载波环路跟踪方法
CN106021181B (zh) * 2016-05-17 2019-06-25 中国电子科技集团公司第四十一研究所 一种对fft数据实施的相邻点相关性均值降噪方法
CN111623785B (zh) * 2020-07-17 2022-02-15 中南大学 基于星间时间延迟量测的深空探测器自主导航方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4605140B2 (ja) * 2006-10-26 2011-01-05 セイコーエプソン株式会社 測位装置、電子機器及びプログラム
CN101718874A (zh) * 2009-12-14 2010-06-02 浙江大学 高灵敏度导航卫星接收机中的信号快速捕获方法及装置
CN101975960B (zh) * 2010-09-17 2013-06-05 桂林电子科技大学 利用gps卫星导航信号实现目标探测和定位的方法及装置

Also Published As

Publication number Publication date
CN102680988A (zh) 2012-09-19

Similar Documents

Publication Publication Date Title
CN102662183B (zh) Gps信号捕获方法与***
CN102680988B (zh) 高灵敏度导航卫星信号非线性捕获方法及装置
CN103645483B (zh) 一种弱信号环境下北斗信号捕获方法
CN102520423B (zh) 卫星导航接收机长周期伪随机扩频码的快速捕获方法
CN106855628A (zh) 一种高动态卫星导航信号的快速捕获和跟踪***和方法
CN100553168C (zh) 一种长周期扩频码的频率折叠时频并行搜索方法
CN105553506B (zh) 一种长码扩频信号的快速捕获方法及装置
CN102914782B (zh) 一种适用于gps弱信号的快速捕获方法
CN106125053A (zh) 脉冲多普勒雷达极化抗干扰方法
CN104833988A (zh) 大多普勒动态下高数据速率直扩信号的捕获***及方法
CN104936287A (zh) 基于矩阵补全的传感网室内指纹定位方法
CN104765050A (zh) 一种新的北斗信号二次捕获算法
CN105005057A (zh) 一种北斗导航***d1导航电文的捕获方法
CN106230475A (zh) 一种基于Tong检测器伪码多普勒补偿捕获方法及装置
CN104483684A (zh) 一种快速捕获北斗d1卫星导航***弱信号的方法
CN103217695A (zh) 一种gps的l1、l2c信号的捕获方法
CN103176189A (zh) 高灵敏度卫星导航接收机的远近效应抑制器及其方法
CN101893698A (zh) 噪声源测试分析方法及其装置
CN103926604A (zh) 基于重叠差分循环相干积分的弱信号捕获方法
CN104270167A (zh) 一种基于多维特征神经网络的信号检测和估计方法
CN103885071A (zh) 一种卫星信号捕获的方法和设备
CN106291610A (zh) 一种用于gnss信号压缩捕获处理装置的压缩并行相关模块及其实现方法
CN101441265B (zh) 利用混沌***进行导航卫星信号捕获的方法
CN103809192A (zh) 一种gnss接收机的动态校正算法
CN101520505A (zh) 全球卫星导航***弱信号下的自适应相干累积联合捕获方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20131106

Termination date: 20150523

EXPY Termination of patent right or utility model