CN102043158A - 微弱卫星导航信号捕获中的信号检测判决方法及其装置 - Google Patents
微弱卫星导航信号捕获中的信号检测判决方法及其装置 Download PDFInfo
- Publication number
- CN102043158A CN102043158A CN2010105922021A CN201010592202A CN102043158A CN 102043158 A CN102043158 A CN 102043158A CN 2010105922021 A CN2010105922021 A CN 2010105922021A CN 201010592202 A CN201010592202 A CN 201010592202A CN 102043158 A CN102043158 A CN 102043158A
- Authority
- CN
- China
- Prior art keywords
- module
- peak
- signal
- accumulation results
- simple crosscorrelation
- 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
Images
Abstract
本发明公开了一种微弱卫星导航信号捕获中的信号检测判决方法及其装置。在弱信号环境下,高灵敏度卫星导航接收机需要对伪码的相关结果进行长时间累积,以达到捕获的灵敏度增益要求。本发明旨在使用基于信号的峰均比判决方法,减少相关结果检测判决过程带来的信噪比损失,并给出了一种判决门限的自适应处理方法。同时本发明考虑远近效应的影响,在初步判决信号存在之后提供排除互相关的判决模块,以此削弱C/A码互相关干扰的影响。控制模块对捕获器进行的二维搜索,采用检测最大值下降趋势的方法判决当前信号存在的多普勒频偏,从而缩短频率维的搜索时间,提高了捕获效率。本发明适用于GPS和BD-2等系统,冷启动捕获灵敏度可达-145dbm。
Description
技术领域
本发明涉及一种微弱卫星导航信号捕获中的信号检测判决方法及其装置。
背景技术
全球卫星导航系统GNSS(Global Navigation Satellite System),作为一种星基导航系统,它利用用户到位置已知的卫星的距离或方向角等信息进行定位和导航. 近年来,随着各大GNSS系统的不断完善,卫星导航技术不仅深入军事、测绘、通信、遥感等军用领域,其在汽车定位系统、个人手持通信终端等民用领域的应用也日益广泛。因其重要的战略地位和无限的市场潜力,卫星定位技术已成为各国研究和发展的重点和热点。
卫星导航接收机通过从足够多的卫星信号中,获得信号传播时间信息和卫星位置信息,根据这些信息计算得到接收机的位置、时间和速度等。为了进行导航定位解算,导航接收机至少需要从不同的4颗卫星得到信息。接收机需要先与接收到的卫星信号完成粗同步,即扩频信号的捕获过程,得到伪随机码的相位偏移和信号载波多普勒偏移。再由跟踪过程得到伪码相位、载波频率/相位的精确估计值,对导航信号进行实时跟踪,解扩后得到该卫星信号中调制的信息。
当接收机位于室内、森林、城市等复杂环境下,导航信号受到严重的衰弱,信号功率将远远低于一般卫星导航接收机的工作范围。微弱信号条件大大限制了卫星导航定位的应用领域。在弱信号环境下,高灵敏度卫星导航接收机需要对伪码的相关结果进行长时间累积,以达到捕获的灵敏度增益要求。无论是非相干累积还是差分相干累积,积分带来的灵敏度增益都是有限的,而基于累积结果幅值的判决过程又会引入灵敏度损失。对于捕获方案本身来说,其所处的弱信号使用环境是不可设定的,因而信号捕获中的噪声基底是难以统计其分布的。另外,接收机使用不同的射频前端也会影响信号捕获结果中噪声基底的大小。因此可以明显看出,基于高斯分布的固定门限判决并不是最佳选择。因此,微弱信号条件下的卫星导航信号捕获中的信号检测判决技术,已成为研究高灵敏度捕获技术不可或缺的组成部分。
发明内容
本发明的目的是克服现有技术的不足,提供一种微弱卫星导航信号捕获中的信号检测判决方法及其装置。
微弱卫星导航信号捕获中的信号检测判决方法的步骤如下:
1)将从捕获器获得的一个扩频码周期的卫星信号累积结果送入检测模块,检测模块检测到数据同步信号之后,同步启动峰值查找模块和均值计算模块,并将该周期内所有累积结果存入存储模块;
2)存储模块、峰值查找模块和均值计算模块将同步结束工作,然后由存储模块将存储的数据导出,连同均值计算模块计算的均值和峰值查找模块找到的峰值,一起送入峰均比计算模块和互相关比例因子计算模块,计算该组捕获累积结果的峰均比和互相关比例因子;
3)控制模块将当前频点的捕获累积结果峰均比与所设比例门限进行比较,判决信号是否存在;对于判决为信号存在的各组累积结果,控制模块同时记录捕获器提供的频点信息,通过比较各频点下各组累积结果峰值的大小,找到目标频点;
4)将目标频点对应的捕获累积结果,送入互相关判决模块,通过互相关比例因子与所设门限的比较,排除远近效率引起的互相关干扰,确认信号是否存在。
所述的将从捕获器获得的一个扩频码周期的卫星信号累积结果送入检测模块。检测模块检测到数据同步信号之后,同步启动峰值查找模块和均值计算模块,并将该周期内所有累积结果存入存储模块步骤,包括:
1)捕获器在完成所需次数的非相干累积或者差分相干累积之后,将卫星信号在一个扩频码周期的累积结果依次输出,同时给出数据有效信号,输出信号为I,Q两路信号取模运算的结果,采样率为系统采样时钟,一个扩频码周期的数据长度为1ms;
2)检测模块检测到数据同步信号的同时,启动峰值查找模块,遍历一个扩频码周期的M个累积结果,找到这M个结果中的最大值即峰值max以及峰值在M个值中所处的位置pos;
3)检测模块检测到数据同步信号的同时,启动均值计算模块,在接收M个累积结果的同时进行M个累积结果求和操作,得到所有M个累积结果的代数和sum,最后将此代数和除以累积结果个数M,同步得到累积结果的均值mean;
4)接收到的数据一边进行峰值查找和均值计算,一边存入存储模块用以进行后续的信号处理,存储模块的存储数据量即为一个扩频码周期累积结果的个数M,接收下一个周期的累积结果时存储模块空间可被覆盖并重复使用。
所述的存储模块、峰值查找模块和均值计算模块将同步结束工作,然后由存储模块将存储的数据导出,连同均值计算模块计算的均值和峰值查找模块找到的峰值,一起送入峰均比计算模块和互相关比例因子计算模块,计算该组捕获累积结果的峰均比和互相关比例因子步骤,包括:
1)存储模块、峰值查找模块和均值计算模块将同步结束工作,即得到M个累积结果的均值、M个累积结果的峰值及其所处位置信息,同时存储模块存满M个数据;
2)存储模块存储数据量达到要求时,打开存储模块读取使能,重新将M个累积结果依序导出;
3)检测到存储模块读使能有效后,开启峰均比计算模块的使能。直至存储模块将所存M个数据全部输出之后,方可得到峰均比的计算结果;
4)开启峰均比计算模块使能的同时开启互相关比例因子计算模块,同样需等存储模块所有M个数据全部输出之后,方可得到该组累积结果的互相关比例因子;
5)峰均比和互相关比例因子计算完成之后,向检测模块发送运算结束标志。
所述的检测到存储模块读使能有效后,开启峰均比计算模块的使能。直至存储模块将所存M个数据全部输出之后,方可得到峰均比的计算结果步骤,包括:
1)将存储模块同步输出的每个累积结果,先减去均值mean,再逐个进行平方累加;
2)完成M次平方累加之后,将累加值除以数据总数M得到方差var,作为近似的统计噪声功率;
3)将找到的峰值max减去均值mean后作平方运算,即消除直流分量影响后统计信号的功率。然后除以统计噪声功率var,得到该组累积结果的峰均比。
4) 计算结束后,送出峰均比结果同时给出计算完成标志。
所述的开启峰均比计算模块使能的同时开启互相关比例因子计算模块,同样需等存储模块所有M个数据全部输出之后,方可得到该组累积结果的互相关比例因子步骤,包括:
1)确定次大值查找的范围,所谓次大值是指除去信号之后噪声采样点的最大值,由于信号的输出会扩展到它邻近的采样点,因此最大值位置pos左右5个采样点都应该从查找范围里剔除;
2)存在3种状态的查找范围,它们分别对应3种剔除方式:
①当峰值位置5<pos<M-1-5,剔除pos±5范围位置的采样点;
②当峰值位置pos≦5,剔除[0,pos+5]和(M-1-(5-pos),M-1]范围位置的采样点;
③当峰值位置pos≧M-1-5, 剔除[pos-5,M-1]和[0,5-(M-1-pos))范围位置的采样点;
3) 在指定范围内找到次大值sub_max,将其与峰值max相减后作平方,然后除以近似噪声功率var,得到互相关比例因子。
所述的控制模块将当前频点的捕获累积结果峰均比与所设比例门限进行比较,判决信号是否存在;对于判决为信号存在的各组累积结果,控制模块同时记录捕获器提供的频点信息,通过比较各频点下各组累积结果峰值的大小,找到目标频点步骤,包括:
1)将当前频点的捕获累积结果峰均比与所设比例门限进行比较,若大于门限则判决信号存在,否则判决信号不存在,与固定门限判决相比,比例门限判决可以适应各种强度条件的信号,等效于给出了一种判决门限的自适应解决方法;
2)一旦判决某个频点上的信号存在,则进入检测峰值下降趋势的等待状态;
3)控制模块记录捕获器提供的频点信息,一旦检测到峰值的下降趋势,则判决发现下降趋势当前频点的前一频点为目标频点。
所述的将目标频点对应的捕获累积结果,送入互相关判决模块。通过互相关比例因子与所设门限的比较,排除远近效率引起的互相关干扰,确认信号是否存在步骤,包括:
1)找到目标频点后,将该频点对应的互相关比例因子送入互相关判决模块;
2)若互相关比例因子大于互相关判决门限,则认为属于自相关结果,确认信号存在;否则认为该结果为互相关干扰下的误捕,重新判决信号不存在;
3)确认信号存在则输出伪码相位和多普勒频点;重新判决信号不存在,则令捕获器进入下一个捕获状态。
微弱卫星导航信号捕获中的信号检测判决装置包括控制模块、检测模块、存储模块、峰值查找模块、均值计算模块、峰均比计算模块、互相关比例因子计算模块、互相关判决模块;捕获器输出的累积结果送入检测模块,检测模块与存储模块,峰值查找模块和均值计算模块相连;存储模块,峰值查找模块和均值计算模块的计算结果送入峰均比计算模块和互相关比例因子计算模块;峰均比计算模块和互相关比例因子计算模块的结果前馈进入检测模块,进行信号存在性判决;控制器完成各模块的接口控制和使能控制,控制不同频点累积结果峰值的下降趋势检测流程以及进入互相关判决模块的时间;互相关判决模块与控制模块相连,输出结果可能前馈到前级捕获器也可能送到后级跟踪通道。
所述的峰均比计算模块,互相关比例因子计算模块都包括方差运算器、平方乘法器和除法器模块,除法器模块分别与方差运算器和平方乘法器相连接。
所述的方差运算器包括顺次相连接的平方乘法器、加法器、除法器和存放累加结果的寄存器dat_mean2_acc。
本发明旨在使用基于信号的峰均比判决方法,减少相关结果检测判决过程带来的信噪比损失,并给出了一种判决门限的自适应处理方法。同时本发明考虑远近效应的影响,在初步判决信号存在之后提供排除互相关的判决模块,以此削弱C/A码互相关干扰的影响。控制模块对捕获器进行的二维搜索,采用检测最大值下降趋势的方法判决当前信号存在的多普勒频偏,从而缩短频率维的搜索时间,提高了捕获效率。本发明适用于GPS和BD-2等系统,冷启动捕获灵敏度可达-145dbm。
附图说明
图1是微弱卫星导航信号捕获中的信号检测判决装置的电路框图;
图2是本发明中峰均比计算模块和互相关比例因子计算模块的电路框图;
图3是本发明中方差运算器的电路框图;
图4是本发明中存储模块、峰值查找模块和均值计算模块完成时刻时序图;
图5是本发明中峰均比计算模块、互相关比例因子计算模块完成时刻时序图。
具体实施方式
本发明实现了微弱卫星导航信号捕获中的信号检测判决装置。从捕获器获得一个扩频码周期的卫星信号累积结果,送入检测模块找到该组累积结果的峰值及其所处的位置标号并计算均值。再利用所得峰值和均值,计算信号功率与噪声基底之间的峰均比和互相关比例因子。然后在控制模块中判决当前频点的信号是否存在,同时通过比较各频点下各组累积结果最大值的大小,找到目标频点。最后对目标频点对应的捕获累积结果进行互相关判决,排除互相关干扰,最终确认信号是否存在。在实施例中,系统的采样时钟为4.092MHz,对于GPS信号,扩频码的码率为1.023MHz;对于BD-2信号则码率为2.046MHz。4.092MHz的采样率对两者都满足奈奎斯特采样定律,两者的扩频码周期均为1ms,这样在一个扩频码周期内累积结果的采样点数为4092个。
微弱卫星导航信号捕获中的信号检测判决方法的步骤如下:
1)将从捕获器获得的一个扩频码周期的卫星信号累积结果送入检测模块,检测模块检测到数据同步信号之后,同步启动峰值查找模块和均值计算模块,并将该周期内所有累积结果存入存储模块;
2)存储模块、峰值查找模块和均值计算模块将同步结束工作,然后由存储模块将存储的数据导出,连同均值计算模块计算的均值和峰值查找模块找到的峰值,一起送入峰均比计算模块和互相关比例因子计算模块,计算该组捕获累积结果的峰均比和互相关比例因子;
3)控制模块将当前频点的捕获累积结果峰均比与所设比例门限进行比较,判决信号是否存在;对于判决为信号存在的各组累积结果,控制模块同时记录捕获器提供的频点信息,通过比较各频点下各组累积结果峰值的大小,找到目标频点;
4)将目标频点对应的捕获累积结果,送入互相关判决模块,通过互相关比例因子与所设门限的比较,排除远近效率引起的互相关干扰,确认信号是否存在。
所述的将从捕获器获得的一个扩频码周期的卫星信号累积结果送入检测模块。检测模块检测到数据同步信号之后,同步启动峰值查找模块和均值计算模块,并将该周期内所有累积结果存入存储模块步骤,包括:
1)捕获器在完成所需次数的非相干累积或者差分相干累积之后,将卫星信号在一个扩频码周期的累积结果依次输出,同时给出数据有效信号,输出信号为I,Q两路信号取模运算的结果,采样率为系统采样时钟,一个扩频码周期的数据长度为1ms;
在实施例中,一个扩频码周期的累积结果为I,Q两路信号取模运算的结果data,共有4092个值。数据有效信号data_valid为高电平有效,即在输出第一个结果的同时输出为高电平,输出最后一个结果后恢复低电平。图1中所示的捕获器并不是本发明装置的组成部分,它只描述输入信号。
2)检测模块检测到数据同步信号的同时,启动峰值查找模块,遍历一个扩频码周期的M个累积结果,找到这M个结果中的最大值即峰值max以及峰值在M个值中所处的位置pos;
实施例中的具体实施方式如下:当检测到数据有效信号data_valid高电平后,启动峰值查找模块。该模块设有一个计数器count,被启动后count信号从0开始计数。另有两个寄存器max和pos,分别寄存当前为止的最大值以及最大值所处位置。当每个采样时钟上升沿到来,则比较当前datain与max的大小。一旦datain大于max,则将max的值更新为datain,同时将pos更新为当前的count值;否则max和pos均保持不变。遍历所有4092个数据之后,宣布查找完成给出done信号有效脉冲。max和pos初始化值为0。
3)检测模块检测到数据同步信号的同时,启动均值计算模块,在接收M个累积结果的同时进行M个累积结果求和操作,得到所有M个累积结果的代数和sum,最后将此代数和除以累积结果个数M,同步得到累积结果的均值mean;
具体实施方式如下:
实施例中M=4092。均值计算模块中有一个寄存器sum存放当前累加和,当每个采样时钟上升沿到来,把sum中已有的结果与当前datain相加后继续存放在sum中。该模块中包括一个除法器,sum直接连接除法器的分子,分母为常数4092。遍历所有4092个数据之后,sum完成所有累加,同时将除法器输出赋给mean。由于所用除法器为组合电路实现,不存在计算时延。
4)接收到的数据一边进行峰值查找和均值计算,一边存入存储模块用以进行后续的信号处理,存储模块的存储数据量即为一个扩频码周期累积结果的个数M,接收下一个周期的累积结果时存储模块空间可被覆盖并重复使用。
实施例中的存储模块采用双口RAM实现,可读可写。对取模后的捕获累积结果使用14bit进行存储,存储模块的大小为4092×14bit=57288bit。存储模块当数据有效信号为高时,从地址0到地址4091进行顺序存储。接收下一组累积结果时重新从地址0开始存储,如此重复使用。
所述的存储模块、峰值查找模块和均值计算模块将同步结束工作,然后由存储模块将存储的数据导出,连同均值计算模块计算的均值和峰值查找模块找到的峰值,一起送入峰均比计算模块和互相关比例因子计算模块,计算该组捕获累积结果的峰均比和互相关比例因子步骤,包括:
1)存储模块、峰值查找模块和均值计算模块将同步结束工作,即得到M个累积结果的均值、M个累积结果的峰值及其所处位置信息,同时存储模块存满M个数据;
如图4时序图所示,峰值查找模块的完成标志done信号脉冲有效的同时,均值计算模块的输出mean有效,存储模块完成存储,写使能置低电平。在实施例中,找到的峰值max为289,所处位置pos为1566,计算得到的均值为42。
2)存储模块存储数据量达到要求时,打开存储模块读取使能,重新将M个累积结果依序导出;
当存储模块存满后,读取使能有效,开始从地址0到地址4091进行顺序读取。每个采样时钟上升沿到来时输出一个ramdata。
3)检测到存储模块读使能有效后,开启峰均比计算模块的使能。直至存储模块将所存M个数据全部输出之后,方可得到峰均比的计算结果;
由于实施例中的存储模块用双口RAM实现,从读使能有效到输出数据有效存在两个系统时钟clk的时延,因此峰均比计算模块的使能因在读取使能有效两个系统时钟clk后开启。
4)开启峰均比计算模块使能的同时开启互相关比例因子计算模块,同样需等存储模块所有M个数据全部输出之后,方可得到该组累积结果的互相关比例因子;
5)峰均比和互相关比例因子计算完成之后,向检测模块发送运算结束标志。
如图5时序图所示,当存储模块的读取使能关闭三个系统时钟clk后,运算结束标志finish发出有效脉冲,峰均比和互相关比例因子计算完成。其中三个系统时钟clk的延时有两个来自RAM的读取延时,另外一个是计算延时。在实施例中,如图5计算得到的峰均比pm_ratio为347,互相关比例因子un_ratio为125。
所述的检测到存储模块读使能有效后,开启峰均比计算模块的使能。直至存储模块将所存M个数据全部输出之后,方可得到峰均比的计算结果步骤,包括:
1)将存储模块同步输出的每个累积结果,先减去均值mean,再逐个进行平方累加;
具体实施方式如下:
实施例中M=4092。如图3所示的方差运算器中设有一个寄存器dat_mean2直接与一个乘法器的输出相连。该乘法器的两个输入均为寄存器dat_mean,即该乘法器完成平方运算。当每个采样时钟上升沿到来时,将当前数据ramdata减去之前mean的结果存放在dat_mean中。由于乘法器用组合电路实现,同步可以得到dat_mean2=dat_mean*dat_mean。有另一寄存器dat_mean2_acc,用来存放当前计算的平方和,把dat_mean2_acc中已有的结果与当前dat_mean2相加后继续存放在dat_mean2_acc中。
2)完成M次平方累加之后,将累加值除以数据总数M得到方差var,作为近似的统计噪声功率;
实施例中M=4092。同如图3方差运算器所示,dat_mean2_acc直接与一个除法器相连作为分子。除法器的分母为常量4092,除法器的输出为变量divout2。当存储模块的读取使能关闭三个系统时钟clk后,除法器的输出有效并赋值给变量var。
3)将找到的峰值max减去均值mean后作平方运算,即消除直流分量影响后统计信号的功率。然后除以统计噪声功率var,得到该组累积结果的峰均比。
实施例中如图2所示,方差运算器的结果作为峰均比计算模块中的除法器的分母。峰值max减去均值mean后作为图2中平方乘法器的输入,该平方乘法器的结果直接接入图2中的除法器作为分子。除法器的输出变量即为峰均比的计算结果pm_ratio。由于图2中各部分均由组合电路实现,因此同样当存储模块的读取使能关闭的三个系统时钟clk后,pm_ratio值有效,计算完成。
4) 计算结束后,送出峰均比结果同时给出计算完成标志。
所述的开启峰均比计算模块使能的同时开启互相关比例因子计算模块,同样需等存储模块所有M个数据全部输出之后,方可得到该组累积结果的互相关比例因子步骤,包括:
1)确定次大值查找的范围,所谓次大值是指除去信号之后噪声采样点的最大值,由于信号的输出会扩展到它邻近的采样点,因此最大值位置pos左右5个采样点都应该从查找范围里剔除;
2)存在3种状态的查找范围,它们分别对应3种剔除方式:
①当峰值位置5<pos<M-1-5,剔除pos±5范围位置的采样点;
②当峰值位置pos≦5,剔除[0,pos+5]和(M-1-(5-pos),M-1]范围位置的采样点;
③当峰值位置pos≧M-1-5, 剔除[pos-5,M-1]和[0,5-(M-1-pos))范围位置的采样点;
在实施例中,互相关比例因子计算模块中设置两个变量excludeRangeIndex1和excludeRangeIndex2,其中excludeRangeIndex1固定赋值为pos-5,excludeRangeIndex2固定赋值为pos+5。最大值位置pos的范围从0到4091。另有两个变量RangeReg1和RangeReg2,表示次大值的查找范围从RangeReg1到RangeReg2。变量Range_State分别对应以上三种剔除方式:
①当pos∈(5,4086)时,RangeReg1= excludeRangeIndex1,RangeReg2= excludeRangeIndex2,Range_State=11。
②当pos∈[0,5]即excludeRangeIndex2超过4091时,则Range_State=10,RangeReg1= excludeRangeIndex2-4091,RangeReg2=excludeRangeIndex1。
③当pos∈[4086,4091]即excludeRangeIndex1为负,则Range_State=01,
RangeReg1= excludeRangeIndex2,RangeReg2=4091+excludeRangeIndex1。
3) 在指定范围内找到次大值sub_max,将其与峰值max相减后作平方,然后除以近似噪声功率var,得到互相关比例因子。
实施例中在RangeReg1和RangeReg2之间的范围内查找次大值,具体实施步骤与峰值查找模块中的实施方法一致,这里不再赘述。当存储模块的读取使能关闭的两个clk后可以找到次大值sub_max。将峰值max减去sub_max后作为如图2中平方乘法器的输入。平方乘法器的输出直接接入除法器的分子,方差运算器的结果接入除法器作为分母,则除法器的输出即为互相关比例因子un_ratio。由于减法操作又引入一个系统时钟clk延时,相当于在存储模块的读取使能关闭的三个系统时钟clk后,un_ratio值有效,与峰均比计算模块同时完成计算。
所述的控制模块将当前频点的捕获累积结果峰均比与所设比例门限进行比较,判决信号是否存在;对于判决为信号存在的各组累积结果,控制模块同时记录捕获器提供的频点信息,通过比较各频点下各组累积结果峰值的大小,找到目标频点步骤,包括:
1)将当前频点的捕获累积结果峰均比与所设比例门限进行比较,若大于门限则判决信号存在,否则判决信号不存在,与固定门限判决相比,比例门限判决可以适应各种强度条件的信号,等效于给出了一种判决门限的自适应解决方法;
在实施例中,针对捕获器采用分段捕获方式的情况,根据不同载噪比分段,使用不同的峰均比门限如下:
对于GPS信号,载噪比从47db/hz到38db/hz,峰均比门限为169;从37db/hz到33db/hz,峰均比门限为121;从32db/hz到28db/hz,峰均比门限为121;从27db/hz到25db/hz,峰均比门限为64。
对于BD-2信号,载噪比从47db/hz到40db/hz,峰均比门限为169;从39db/hz到36db/hz,峰均比门限为64;从35db/hz到31db/hz,峰均比门限为64;从30db/hz到26db/hz,峰均比门限为36。
以上所使用的峰均比门限都是根据仿真得到的,并且经过实测验证。
2)一旦判决某个频点上的信号存在,则进入检测峰值下降趋势的等待状态;
3)控制模块记录捕获器提供的频点信息,一旦检测到峰值的下降趋势,则判决发现下降趋势当前频点的前一频点为目标频点。
所述的将目标频点对应的捕获累积结果,送入互相关判决模块。通过互相关比例因子与所设门限的比较,排除远近效率引起的互相关干扰,确认信号是否存在步骤,包括:
1)找到目标频点后,将该频点对应的互相关比例因子送入互相关判决模块;
2)若互相关比例因子大于互相关判决门限,则认为属于自相关结果,确认信号存在;否则认为该结果为互相关干扰下的误捕,重新判决信号不存在;
在实施例中,针对捕获器采用分段捕获方式的情况,根据不同载噪比分段,使用的不同的互相关判决门限如下:
对于GPS信号,载噪比从47db/hz到27db/hz,互相关判决门限为14;从27db/hz到25db/hz,互相关判决门限为6。
对于BD-2信号,载噪比从47db/hz到31db/hz,互相关判决门限为14;从30db/hz到26db/hz,互相关判决门限为6。
以上所使用的互相关判决门限都是根据仿真得到的,并且经过实测验证。
3)确认信号存在则输出伪码相位和多普勒频点;重新判决信号不存在,则令捕获器进入下一个捕获状态。
微弱卫星导航信号捕获中的信号检测判决装置包括控制模块、检测模块、存储模块、峰值查找模块、均值计算模块、峰均比计算模块、互相关比例因子计算模块、互相关判决模块;捕获器输出的累积结果送入检测模块,检测模块与存储模块,峰值查找模块和均值计算模块相连;存储模块,峰值查找模块和均值计算模块的计算结果送入峰均比计算模块和互相关比例因子计算模块;峰均比计算模块和互相关比例因子计算模块的结果前馈进入检测模块,进行信号存在性判决;控制器完成各模块的接口控制和使能控制,控制不同频点累积结果峰值的下降趋势检测流程以及进入互相关判决模块的时间;互相关判决模块与控制模块相连,输出结果可能前馈到前级捕获器也可能送到后级跟踪通道。
所述的峰均比计算模块,互相关比例因子计算模块都包括方差运算器、平方乘法器和除法器模块,除法器模块分别与方差运算器和平方乘法器相连接。
所述的方差运算器包括顺次相连接的平方乘法器、加法器、除法器和存放累加结果的寄存器dat_mean2_acc。
以上所述实施例所使用的峰均比门限和互相关判决门限,为在不同载噪比条件下仿真得到,进过实测验证的。任何熟悉本领域的技术人员可以根据所采用的不同捕获器,修改本发明所用门限,实现其他场合的应用。本发明仅由后附权利要求书及其等效技术方案来限定,我们要求将这些作为本发明来保护。
Claims (10)
1. 一种微弱卫星导航信号捕获中的信号检测判决方法,其特征在于它的步骤如下:
1)将从捕获器获得的一个扩频码周期的卫星信号累积结果送入检测模块,检测模块检测到数据同步信号之后,同步启动峰值查找模块和均值计算模块,并将该周期内所有累积结果存入存储模块;
2)存储模块、峰值查找模块和均值计算模块将同步结束工作,然后由存储模块将存储的数据导出,连同均值计算模块计算的均值和峰值查找模块找到的峰值,一起送入峰均比计算模块和互相关比例因子计算模块,计算该组捕获累积结果的峰均比和互相关比例因子;
3)控制模块将当前频点的捕获累积结果峰均比与所设比例门限进行比较,判决信号是否存在;对于判决为信号存在的各组累积结果,控制模块同时记录捕获器提供的频点信息,通过比较各频点下各组累积结果峰值的大小,找到目标频点;
4)将目标频点对应的捕获累积结果,送入互相关判决模块,通过互相关比例因子与所设门限的比较,排除远近效率引起的互相关干扰,确认信号是否存在。
2.根据权利要求1所述的一种微弱卫星导航信号捕获中的信号检测判决方法,其特征在于,所述的将从捕获器获得的一个扩频码周期的卫星信号累积结果送入检测模块;检测模块检测到数据同步信号之后,同步启动峰值查找模块和均值计算模块,并将该周期内所有累积结果存入存储模块步骤,包括:
1)捕获器在完成所需次数的非相干累积或者差分相干累积之后,将卫星信号在一个扩频码周期的累积结果依次输出,同时给出数据有效信号,输出信号为I,Q两路信号取模运算的结果,采样率为系统采样时钟,一个扩频码周期的数据长度为1ms;
2)检测模块检测到数据同步信号的同时,启动峰值查找模块,遍历一个扩频码周期的M个累积结果,找到这M个结果中的最大值即峰值max以及峰值在M个值中所处的位置pos;
3)检测模块检测到数据同步信号的同时,启动均值计算模块,在接收M个累积结果的同时进行M个累积结果求和操作,得到所有M个累积结果的代数和sum,最后将此代数和除以累积结果个数M,同步得到累积结果的均值mean;
4)接收到的数据一边进行峰值查找和均值计算,一边存入存储模块用以进行后续的信号处理,存储模块的存储数据量即为一个扩频码周期累积结果的个数M,接收下一个周期的累积结果时存储模块空间可被覆盖并重复使用。
3. 根据权利要求1所述的一种微弱卫星导航信号捕获中的信号检测判决方法,其特征在于,所述的存储模块、峰值查找模块和均值计算模块将同步结束工作,然后由存储模块将存储的数据导出,连同均值计算模块计算的均值和峰值查找模块找到的峰值,一起送入峰均比计算模块和互相关比例因子计算模块,计算该组捕获累积结果的峰均比和互相关比例因子步骤,包括:
1)存储模块、峰值查找模块和均值计算模块将同步结束工作,即得到M个累积结果的均值、M个累积结果的峰值及其所处位置信息,同时存储模块存满M个数据;
2)存储模块存储数据量达到要求时,打开存储模块读取使能,重新将M个累积结果依序导出;
3)检测到存储模块读使能有效后,开启峰均比计算模块的使能,直至存储模块将所存M个数据全部输出之后,方可得到峰均比的计算结果;
4)开启峰均比计算模块使能的同时开启互相关比例因子计算模块,同样需等存储模块所有M个数据全部输出之后,方可得到该组累积结果的互相关比例因子;
5)峰均比和互相关比例因子计算完成之后,向检测模块发送运算结束标志。
4. 根据权利要求3所述的一种微弱卫星导航信号捕获中的信号检测判决方法,其特征在于,所述的检测到存储模块读使能有效后,开启峰均比计算模块的使能;
直至存储模块将所存M个数据全部输出之后,方可得到峰均比的计算结果步骤,包括:
1)将存储模块同步输出的每个累积结果,先减去均值mean,再逐个进行平方累加;
2)完成M次平方累加之后,将累加值除以数据总数M得到方差var,作为近似的统计噪声功率;
3)将找到的峰值max减去均值mean后作平方运算,即消除直流分量影响后统计信号的功率;
然后除以统计噪声功率var,得到该组累积结果的峰均比;
4) 计算结束后,送出峰均比结果同时给出计算完成标志;
根据权利要求3所述的一种微弱卫星导航信号捕获中的信号检测判决方法,其特征在于,所述的开启峰均比计算模块使能的同时开启互相关比例因子计算模块,同样需等存储模块所有M个数据全部输出之后,方可得到该组累积结果的互相关比例因子步骤,包括:
1)确定次大值查找的范围,所谓次大值是指除去信号之后噪声采样点的最大值,由于信号的输出会扩展到它邻近的采样点,因此最大值位置pos左右5个采样点都应该从查找范围里剔除;
2)存在3种状态的查找范围,它们分别对应3种剔除方式:
①当峰值位置5<pos<M-1-5,剔除pos±5范围位置的采样点;
②当峰值位置pos≦5,剔除[0,pos+5]和(M-1-(5-pos),M-1]范围位置的采样点;
③当峰值位置pos≧M-1-5, 剔除[pos-5,M-1]和[0,5-(M-1-pos))范围位置的采样点;
3) 在指定范围内找到次大值sub_max,将其与峰值max相减后作平方,然后除以近似噪声功率var,得到互相关比例因子。
5. 根据权利要求1所述的一种微弱卫星导航信号捕获中的信号检测判决方法,其特征在于,所述的控制模块将当前频点的捕获累积结果峰均比与所设比例门限进行比较,判决信号是否存在;对于判决为信号存在的各组累积结果,控制模块同时记录捕获器提供的频点信息,通过比较各频点下各组累积结果峰值的大小,找到目标频点步骤,包括:
1)将当前频点的捕获累积结果峰均比与所设比例门限进行比较,若大于门限则判决信号存在,否则判决信号不存在,与固定门限判决相比,比例门限判决可以适应各种强度条件的信号,等效于给出了一种判决门限的自适应解决方法;
2)一旦判决某个频点上的信号存在,则进入检测峰值下降趋势的等待状态;
3)控制模块记录捕获器提供的频点信息,一旦检测到峰值的下降趋势,则判决发现下降趋势当前频点的前一频点为目标频点。
6. 根据权利要求1所述的一种微弱卫星导航信号捕获中的信号检测判决方法,其特征在于,所述的将目标频点对应的捕获累积结果,送入互相关判决模块。
7.通过互相关比例因子与所设门限的比较,排除远近效率引起的互相关干扰,确认信号是否存在步骤,包括:
1)找到目标频点后,将该频点对应的互相关比例因子送入互相关判决模块;
2)若互相关比例因子大于互相关判决门限,则认为属于自相关结果,确认信号存在;否则认为该结果为互相关干扰下的误捕,重新判决信号不存在;
3)确认信号存在则输出伪码相位和多普勒频点;重新判决信号不存在,则令捕获器进入下一个捕获状态。
8. 一种如权利要求1所述方法设计的微弱卫星导航信号捕获中的信号检测判决装置,其特征在于,包括控制模块、检测模块、存储模块、峰值查找模块、均值计算模块、峰均比计算模块、互相关比例因子计算模块、互相关判决模块;捕获器输出的累积结果送入检测模块,检测模块与存储模块,峰值查找模块和均值计算模块相连;存储模块,峰值查找模块和均值计算模块的计算结果送入峰均比计算模块和互相关比例因子计算模块;峰均比计算模块和互相关比例因子计算模块的结果前馈进入检测模块,进行信号存在性判决;控制器完成各模块的接口控制和使能控制,控制不同频点累积结果峰值的下降趋势检测流程以及进入互相关判决模块的时间;互相关判决模块与控制模块相连,输出结果可能前馈到前级捕获器也可能送到后级跟踪通道。
9.根据权利要求8所述的一种微弱卫星导航信号捕获中的信号检测判决装置,其特征在于,所述的峰均比计算模块,互相关比例因子计算模块都包括方差运算器、平方乘法器和除法器模块,除法器模块分别与方差运算器和平方乘法器相连接。
10.根据权利要求8所述的一种微弱卫星导航信号捕获中的信号检测判决装置,其特征在于,所述的方差运算器包括顺次相连接的平方乘法器、加法器、除法器和存放累加结果的寄存器dat_mean2_acc。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010592202 CN102043158B (zh) | 2010-12-17 | 2010-12-17 | 微弱卫星导航信号捕获中的信号检测判决方法及其装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010592202 CN102043158B (zh) | 2010-12-17 | 2010-12-17 | 微弱卫星导航信号捕获中的信号检测判决方法及其装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102043158A true CN102043158A (zh) | 2011-05-04 |
CN102043158B CN102043158B (zh) | 2013-04-24 |
Family
ID=43909484
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010592202 Active CN102043158B (zh) | 2010-12-17 | 2010-12-17 | 微弱卫星导航信号捕获中的信号检测判决方法及其装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102043158B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102761503A (zh) * | 2012-07-03 | 2012-10-31 | 华为技术有限公司 | 信道速度估计方法及装置 |
CN103176189A (zh) * | 2013-03-08 | 2013-06-26 | 浙江大学 | 高灵敏度卫星导航接收机的远近效应抑制器及其方法 |
CN104101883A (zh) * | 2013-04-11 | 2014-10-15 | 中国科学院微电子研究所 | 一种全球卫星导航信号的快速捕获方法及装置 |
CN104101884A (zh) * | 2014-07-31 | 2014-10-15 | 豪芯微电子科技(上海)有限公司 | 卫星信号的捕获方法及接收机 |
CN104730337A (zh) * | 2015-03-06 | 2015-06-24 | 华东师范大学 | 一种基于频谱图的信号检测方法 |
CN105656511A (zh) * | 2016-01-20 | 2016-06-08 | 上海物联网有限公司 | 一种适应于有频偏和低信噪比环境下的差分相关捕获方法 |
CN105699990A (zh) * | 2015-12-31 | 2016-06-22 | 西安航天华迅科技有限公司 | 一种gnss微弱信号的载波跟踪环路及跟踪方法 |
CN106597492A (zh) * | 2017-02-06 | 2017-04-26 | 桂林电子科技大学 | 卫星导航接收机及其抗远近效应的方法和室内定位方法 |
CN106646548A (zh) * | 2016-10-27 | 2017-05-10 | 湖南易通星云科技发展有限公司 | 一种多模的超快速、低功耗的定位接收机系统及方法 |
CN107015253A (zh) * | 2017-03-09 | 2017-08-04 | 沈阳理工大学 | 一种导航接收机中可动态威胁预测的搜索检测方法 |
CN108663694A (zh) * | 2018-03-23 | 2018-10-16 | 交通运输部水运科学研究所 | 基带芯片性能测试方法和装置 |
CN109150407A (zh) * | 2018-09-25 | 2019-01-04 | 中国科学院国家授时中心 | 多路复用信号伪随机码盲估方法、误码纠错方法及装置 |
CN110291421A (zh) * | 2017-02-13 | 2019-09-27 | 株式会社久保田 | 用于作业车的卫星电波灵敏度分布管理系统及方法 |
CN111343118A (zh) * | 2019-10-31 | 2020-06-26 | 北京昂瑞微电子技术有限公司 | 一种同步电路及同步检测方法 |
CN112073990A (zh) * | 2020-09-09 | 2020-12-11 | 重庆重邮汇测电子技术研究院有限公司 | 一种串行数据流中峰均比检测实现方法 |
CN113608242A (zh) * | 2021-06-18 | 2021-11-05 | 西安空间无线电技术研究所 | 一种基于码周期扩频码认证的导航信号安全性增强方法 |
CN113890562A (zh) * | 2021-09-10 | 2022-01-04 | 深圳市力合微电子股份有限公司 | 一种低功耗微功率无线通信的扩频信号捕获方法 |
CN115291258A (zh) * | 2022-10-08 | 2022-11-04 | 成都星航时空科技有限公司 | 一种gnss基带捕获方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070085736A1 (en) * | 2005-10-14 | 2007-04-19 | Accord Software & Systems Pvt. Ltd. | Weak signal acquisition |
US20090129449A1 (en) * | 2007-11-20 | 2009-05-21 | Seiko Epson Corporation | Received signal determination method, positioning device, and electronic instrument |
CN101581776A (zh) * | 2009-06-19 | 2009-11-18 | 哈尔滨工业大学 | 一种微弱gnss信号的差分相干累积捕获方法 |
CN101793968A (zh) * | 2010-03-10 | 2010-08-04 | 北京航空航天大学 | 一种适用于微弱卫星信号捕获的双门限检测方法 |
-
2010
- 2010-12-17 CN CN 201010592202 patent/CN102043158B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070085736A1 (en) * | 2005-10-14 | 2007-04-19 | Accord Software & Systems Pvt. Ltd. | Weak signal acquisition |
US20090129449A1 (en) * | 2007-11-20 | 2009-05-21 | Seiko Epson Corporation | Received signal determination method, positioning device, and electronic instrument |
CN101581776A (zh) * | 2009-06-19 | 2009-11-18 | 哈尔滨工业大学 | 一种微弱gnss信号的差分相干累积捕获方法 |
CN101793968A (zh) * | 2010-03-10 | 2010-08-04 | 北京航空航天大学 | 一种适用于微弱卫星信号捕获的双门限检测方法 |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102761503B (zh) * | 2012-07-03 | 2015-11-25 | 华为技术有限公司 | 信道速度估计方法及装置 |
CN102761503A (zh) * | 2012-07-03 | 2012-10-31 | 华为技术有限公司 | 信道速度估计方法及装置 |
CN103176189A (zh) * | 2013-03-08 | 2013-06-26 | 浙江大学 | 高灵敏度卫星导航接收机的远近效应抑制器及其方法 |
CN104101883A (zh) * | 2013-04-11 | 2014-10-15 | 中国科学院微电子研究所 | 一种全球卫星导航信号的快速捕获方法及装置 |
CN104101883B (zh) * | 2013-04-11 | 2016-08-03 | 中国科学院微电子研究所 | 一种全球卫星导航信号的快速捕获方法及装置 |
CN104101884B (zh) * | 2014-07-31 | 2017-02-15 | 上海北伽导航科技有限公司 | 卫星信号的捕获方法及接收机 |
CN104101884A (zh) * | 2014-07-31 | 2014-10-15 | 豪芯微电子科技(上海)有限公司 | 卫星信号的捕获方法及接收机 |
CN104730337A (zh) * | 2015-03-06 | 2015-06-24 | 华东师范大学 | 一种基于频谱图的信号检测方法 |
CN104730337B (zh) * | 2015-03-06 | 2017-07-25 | 华东师范大学 | 一种基于频谱图的信号检测方法 |
CN105699990B (zh) * | 2015-12-31 | 2018-09-28 | 西安航天华迅科技有限公司 | 一种gnss微弱信号的载波跟踪环路及跟踪方法 |
CN105699990A (zh) * | 2015-12-31 | 2016-06-22 | 西安航天华迅科技有限公司 | 一种gnss微弱信号的载波跟踪环路及跟踪方法 |
CN105656511A (zh) * | 2016-01-20 | 2016-06-08 | 上海物联网有限公司 | 一种适应于有频偏和低信噪比环境下的差分相关捕获方法 |
CN105656511B (zh) * | 2016-01-20 | 2020-04-07 | 上海物联网有限公司 | 一种适应于有频偏和低信噪比环境下的差分相关捕获方法 |
CN106646548A (zh) * | 2016-10-27 | 2017-05-10 | 湖南易通星云科技发展有限公司 | 一种多模的超快速、低功耗的定位接收机系统及方法 |
CN106597492B (zh) * | 2017-02-06 | 2023-09-05 | 桂林电子科技大学 | 卫星导航接收机及其抗远近效应的方法和室内定位方法 |
CN106597492A (zh) * | 2017-02-06 | 2017-04-26 | 桂林电子科技大学 | 卫星导航接收机及其抗远近效应的方法和室内定位方法 |
CN110291421B (zh) * | 2017-02-13 | 2023-07-28 | 株式会社久保田 | 用于作业车的卫星电波灵敏度分布管理系统及方法 |
CN110291421A (zh) * | 2017-02-13 | 2019-09-27 | 株式会社久保田 | 用于作业车的卫星电波灵敏度分布管理系统及方法 |
CN107015253B (zh) * | 2017-03-09 | 2021-06-18 | 沈阳理工大学 | 一种导航接收机中可动态威胁预测的搜索检测方法 |
CN107015253A (zh) * | 2017-03-09 | 2017-08-04 | 沈阳理工大学 | 一种导航接收机中可动态威胁预测的搜索检测方法 |
CN108663694A (zh) * | 2018-03-23 | 2018-10-16 | 交通运输部水运科学研究所 | 基带芯片性能测试方法和装置 |
CN109150407A (zh) * | 2018-09-25 | 2019-01-04 | 中国科学院国家授时中心 | 多路复用信号伪随机码盲估方法、误码纠错方法及装置 |
CN109150407B (zh) * | 2018-09-25 | 2021-12-17 | 中国科学院国家授时中心 | 多路复用信号伪随机码盲估方法、误码纠错方法及装置 |
CN111343118A (zh) * | 2019-10-31 | 2020-06-26 | 北京昂瑞微电子技术有限公司 | 一种同步电路及同步检测方法 |
CN112073990A (zh) * | 2020-09-09 | 2020-12-11 | 重庆重邮汇测电子技术研究院有限公司 | 一种串行数据流中峰均比检测实现方法 |
CN112073990B (zh) * | 2020-09-09 | 2023-03-03 | 重庆重邮汇测电子技术研究院有限公司 | 一种串行数据流中峰均比检测实现方法 |
CN113608242B (zh) * | 2021-06-18 | 2023-08-11 | 西安空间无线电技术研究所 | 一种基于码周期扩频码认证的导航信号安全性增强方法 |
CN113608242A (zh) * | 2021-06-18 | 2021-11-05 | 西安空间无线电技术研究所 | 一种基于码周期扩频码认证的导航信号安全性增强方法 |
CN113890562A (zh) * | 2021-09-10 | 2022-01-04 | 深圳市力合微电子股份有限公司 | 一种低功耗微功率无线通信的扩频信号捕获方法 |
CN115291258A (zh) * | 2022-10-08 | 2022-11-04 | 成都星航时空科技有限公司 | 一种gnss基带捕获方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102043158B (zh) | 2013-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102043158B (zh) | 微弱卫星导航信号捕获中的信号检测判决方法及其装置 | |
CN102608626B (zh) | 一种高灵敏度卫星导航信号捕获方法及系统 | |
CN101029923B (zh) | 全球定位系统室内定位与导航技术的软件接收机 | |
CN103885057A (zh) | 自适应变滑窗多目标跟踪方法 | |
CN103439697B (zh) | 基于动态规划的目标检测方法 | |
CN107741231B (zh) | 一种基于机器视觉的多运动目标快速测距方法 | |
CN104749594B (zh) | 一种gps双频非差周跳探测与修复方法及装置 | |
CN105487060B (zh) | 一种双通道四斜率调制的多目标提取方法 | |
CN102722706B (zh) | 基于粒子滤波的红外弱小目标检测与跟踪方法及装置 | |
CN104880717B (zh) | 一种卫星测量值多路径误差检测装置及其算法 | |
CN102155955A (zh) | 立体视觉里程计及测量方法 | |
CN106796295B (zh) | 用于在多路径误差方面评估全球导航卫星系统中的卫星信号的方法、用于全球导航卫星系统的接收器和机动车 | |
CN104833988A (zh) | 大多普勒动态下高数据速率直扩信号的捕获系统及方法 | |
CN104730537A (zh) | 基于多尺度模型的红外/激光雷达数据融合目标跟踪方法 | |
CN111983655B (zh) | 城市峡谷区域的确定方法、装置、电子设备及存储介质 | |
CN103176189A (zh) | 高灵敏度卫星导航接收机的远近效应抑制器及其方法 | |
CN111999735B (zh) | 一种基于径向速度和目标跟踪的动静目标分离方法 | |
CN108919311A (zh) | 用于北斗导航系统的抗干扰方法 | |
CN104931980B (zh) | 载波相位测量半周模糊解除方法 | |
CN104407345B (zh) | 一种用于弱目标检测跟踪中的改进型动态规划方法 | |
CN103901427A (zh) | 一种测速雷达多目标跟踪的方法和装置 | |
CN106054221A (zh) | 一种位同步方法及其装置 | |
CN110488273A (zh) | 一种基于雷达的车辆跟踪检测方法及装置 | |
CN103235325A (zh) | 基于线段二分法的镜面反射点估计算法 | |
CN105204044B (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 |