基于贝叶斯估计的长码捕获方法与流程

本发明涉及卫星导航接收机的长码捕获方法。本发明可用于实现微弱卫星导航信号的快速长码捕获。

背景技术：

目前主流卫星导航接收机的长码捕获方法，主要是基于匹配滤波器的方式，但限于匹配滤波频率搜索范围较窄，解决的方式是采用部分匹配滤波级联fft的方式进行频率域搜索，这种方式受限于码多普勒，匹配滤波的长度不宜过长，长码捕获的灵敏度较低，通常的方式是增加非相干累积，而这又不可避免的延长捕获所需的搜索时间。

技术实现要素：

本发明的目的在于克服上述现有技术的不足，提供一种基于贝叶斯估计的长码捕获方法，在较少的硬件资源条件下实现较大多普勒频移微弱卫星导航信号环境下长码信号的快速捕获。

本发明的目的是这样实现的：基于贝叶斯估计的长码捕获方法，依次完成数字下变频、匹配滤波、相位旋转、门限判决、非相干累积以及二次门限判决，包括以下步骤：

(1)根据本地恢复载波对数字中频信号进行下变频并进行整码片累加变为两组累加信号：准时累加信号和延时累加信号；

(2)对步骤(1)得到的累加信号和延时累加信号分别进行多组并行的fft变换；

(3)对本地pn码进行多组并行的fft变换；

(4)对步骤(2)和步骤(3)得到的fft处理结果进行对位的共轭相乘；

(5)对步骤(4)得到的共轭相乘结果进行多组并行的ifft变换；

(6)对步骤(5)得到的ifft处理结果按组与不同的相位旋转系数相乘，并将结果在组间对位累积；

(7)对步骤(6)得到的累积结果进行门限判决，更换相位旋转系数，返回步骤(6)，直到遍历完设定的频率范围；

(8)对通过门限判决的码相位和多普勒频移进行非相干累加；

(9)对步骤(8)得到的非相干累加结果进行二次门限判决，如通过判决则捕获完成，否则返回步骤(1)搜索下一组数字中频信号。

其中，步骤(2)具体为：

对步骤(1)得到的准时累加信号s和延时累加信号sdelay分别取(2×m-1)×n/2个数据，按照50％数据重叠的方式分为n个为一组的2m个组，并行进行2m个n点的fft，准时累加信号s的fft结果为sfft、sfft-d两组信号，延时累加信号sdelay的fft结果为sdelay-fft、sdelay-fft-d两组信号，其中n＝2ⁿ，n为正整数，m为正整数。

其中，步骤(3)具体为：

对本地pn码序列取m×n个数据，分为n/2个为一组的2m个组，每组后面n/2点补零，进行n点的fft，产生pnfft信号。

其中，步骤(8)具体为：

根据通过门限判决的码相位选择准时累加信号s和延时累加信号sdelay之一，对选出的累加信号按通过门限判决的码相位进行偏移，用本地pn码序列与偏移后的累加信号对位相乘，然后分别与对应的相位旋转因子e^-jnω相乘、分组求和并取模，再将模值进行累加；其中n∈[1，n]，n＝2ⁿ，n为正整数。

本发明与现有技术相比的有益效果为：

(1)在微弱卫星导航信号环境下，能保证快速的完成长码捕获；

(2)在对信号捕获灵敏度、频率搜索范围及捕获时间均有要求的条件下，能够显著降低资源消耗。

附图说明

图1是本发明对采样信号进行累加分为准时累加信号/延时累加信号的示意图；

图2是本发明对累加信号进行fft处理的示意图；

图3是本发明对fft结果进行相位旋转的原理示意图；

图4是本发明采用相关器的方式进行验证的示意图。

具体实施方式

以下结合附图对本发明做进一步的说明。

基于贝叶斯估计的长码捕获方法，依次完成数字下变频、匹配滤波、相位旋转、门限判决、非相干累积以及二次门限判决，包括以下步骤：

(1)根据本地恢复载波对数字中频信号进行下变频，并进行整码片累加变为两组信号：准时累加信号s、延时累加信号sdelay，整码片累加的方式如图1所示；

(2)对步骤(1)得到的准时累加信号和延时累加信号分别取(2×m-1)×n/2个数据，按照50％数据重叠的方式分为n个为一组的2m个组，并行进行2m个n点的fft，准时累加信号s的fft结果为sfft、sfft-d两组信号，延时累加信号sdelay的fft结果为sdelay-fft、sdelay-fft-d两组信号，其中n＝2ⁿ，n为正整数，m为正整数；

50％数据重叠的fft分组方式如图2所示；

(3)对本地pn码序列取m×n个数据，分为n/2个为一组的2m个组，每组后面n/2点补零，进行n点的fft，产生pnfft信号；

(4)将sfft、sfft-d、sdelay-fft、sdelay-fft-d四组信号取共轭的结果分别与pnfft对位相乘，产生xfft、xfft-d、xdelay-fft、xdelay-fft-d四组信号；

(5)分别对xfft、xfft-d、xdelay-fft、xdelay-fft-d四组信号并行进行4m个n点的ifft，取每组ifft的前n/2点的结果；

(6)对ifft的结果进行相位旋转和累加，对于m组结果，分别与对应的相位旋转因子e^-jmω(m∈[1，m])相乘并求和；

(7)对累加结果的模值进行门限判决，记录通过判决的码相位以及多普勒频移；更换相位旋转因子的ω值并重复步骤(6)，直到遍历完所需的频率范围；

相位旋转和累加方式如图3所示；

(8)根据通过门限判决的码相位选择准时累加信号s和延时累加信号sdelay之一，对选出的累加信号按通过门限判决的码相位进行偏移p，用本地pn码序列与偏移后的累加信号对位相乘，然后分别与对应的相位旋转因子e^-jnω相乘、分组求和并取模，再将模值进行累加；其中n∈[1，n]，n＝2ⁿ，n为正整数。

(9)将模值进行累加进行二次门限判决，若通过二次门限判决，则捕获完成，若均未通过判决，则返回步骤(1)搜索下一组数字中频信号。用相关器进行验证的方式如图4所示。

技术特征：

技术总结
本发明公开了一种基于贝叶斯估计的长码捕获方法，它采用匹配滤波级联相位旋转的方法，对存在一定载波动态的长码信号进行相关峰捕获，同时考虑到弱信号的捕获，对出现的若干相关峰采用验证的方式有效的剔除虚警。该方法的特点是捕获时间快，适用于资源受限且对捕获时间有较高要求的情况。

技术研发人员：张丹;杜辉
受保护的技术使用者：中国电子科技集团公司第五十四研究所
技术研发日：2017.07.28
技术公布日：2017.11.24