GNSS导航卫星信号抗干扰处理方法与流程

本发明涉及信号处理与导航领域，特别涉及一种民用领域卫星导航抗干扰信号处理技术。

背景技术：

全球导航卫星系统(gnss)距离地球约2万多公里，到达地球的信号功率很小，从而使得卫星导航信号接收系统很容易受到干扰。在民用领域，随着各种通信系统和无线数据传输系统的不断增加，电磁环境也愈加复杂，尽管这些系统本身可能不在gnss频段内，但信号的互调产品和带外传输都有可能对gnss信号的接收造成干扰，从而导致导航精度的降低或者接收机的完全失锁。

目前gnss接收机系统的抗干扰技术主要包括有:空域自适应滤波技术，时/频域滤波，空时/频滤波自适应滤波技术，直接p码捕获技术，抗多径技术，星上抗干扰技术，组合导航等抗干扰技术等。时域处理技术能处理多个窄带干扰，但对宽带干扰通常效果不佳；频域滤波对窄带干扰抑制可达35db以上，但对宽带噪音干扰及多个扫频瞄准式噪声干扰无效；空域滤波即是将自适应零陷技术应用于导航终端抗干扰，阵列自适应零陷是将包括多个阵元的天线阵所接收的信号用增益和相位可调的权值进行加权，从而在天线方向图中产生对着干扰方向的零点。在复杂的电磁环境下,采用单一的抗干扰方式已经无法满足接收机的需求。单纯的空域滤波技术是建立在窄带假设基础上，空域滤波算法只能有效抑制窄带干扰，而空时自适应处理(stap)是通过空时联合处理多阵元(空域)与多个时域接收到的数据，使干扰抑制在空时二维空间中进行。目前没有哪一种抗干扰技术能够解决所有的gnss干扰问题。针对繁多复杂的导航干扰环境，为了进一步提高gnss接收机生存能力及使用性能。未来gnss接收机抗干扰方案必然是多种抗干扰方法的综合使用。

随着我国北斗全球卫星导航系统(bds)的逐步建立，探索和研究北斗导航系统抗干扰能力提升的方法，具有重要应用和经济价值。

技术实现要素：

本发明所要解决的技术问题是：提升复杂电磁环境下的卫星导航信号接收系统的抗干扰能力，为此本发明提出一种基于自适应阵列信号处理的导航卫星信号抗干扰方法。

本发明的上述目的可以通过如下技术方案予以实现：一种gnss导航卫星信号抗干扰方法，具有如下技术特征：多阵元天线将接收到的gnss导航卫星射频信号，首先通过模数/数模(ad/da)处理芯片转换为数字中频采样信号，ad/da处理芯片将数字中频采样信号送入信号与数据处理模块，经过数字下变频处理后对i、q采样数据流进行控制和中频数据缓存抽取控制，其中缓存数据中前1ms用于数据抽取，后1ms缓存数据等待，通过对抽取的i、q采样数据和其共轭转置的乘积求解统计平均，得到自相关矩阵；再根据功率倒置(pi)准则，计算数字波束合成(dbf)最优权值系数；将dbf权值系数与对应时刻的i、q采样数据进行复数相乘，得到dbf合成处理的i、q数据；再将dbf合成后的数字中频信号输入ad/da处理芯片，进行模拟上变频处理，上变频为gnss射频信号，并进行功率匹配设置，得到抗干扰处理后的gnss射频信号送至接收机进行导航定位。

本发明相比于现有技术具有如下有益效果：本发明基于多阵元卫星导航天线，模数-数模转换处理芯片，信号与数据处理模块，gnss接收机等组成单元。通过ad-da处理芯片对阵列天线接收到的gnss射频信号中的一种频率信号的混频处理，同时获得两个频点信号的抗干扰能力；通过设计最大深度缓存队列缓存2ms的数据，权值计算完成后可以实现与对应时刻的采样数据的适配性处理。

附图说明

图1是本发明gnss导航卫星信号抗干扰信号处理流程示意图。

图2为gnss射频信号采样与下变频处理流程图。

图3为gnss抗干扰的中频数据抽取与缓存处理流程图。

具体实施方式

下面结合具体实施方式对本发明做进一步详细阐述。

参阅图1。根据本发明，阵列天线将接收到的gnss导航卫星射频信号，首先通过ad/da处理芯片转换为数字中频采样信号，ad/da处理芯片将数字中频采样信号送入信号与数据处理模块，并经过下变频处理后对i、q采样数据流进行缓存与抽取控制，并通过对选取的i、q采样数据和其共轭转置的乘积求解统计平均，得到自相关矩阵；再根据功率倒置(pi)准则，计算数字波束合成dbf最优权值系数；将dbf权值系数与对应时刻的i、q采样数据进行复数相乘，得到dbf合成处理的i、q数据；将dbf合成后的数字中频信号输入ad/da处理芯片，进行模拟上变频处理，上变频为gnss射频信，并进行功率匹配设置，得到抗干扰处理后的gnss射频信号送至接收机进行导航定位。其中，gnss接收机采用多阵元天线对视野内的导航卫星的信号进行接收。多阵元天线可以是四阵元无源天线，阵元数量包含但不仅限于四阵元阵列天线，还包括其他阵元数量的阵列天线。gnss接收机可以是普通民用接收机；ad-da处理芯片和信号与数据处理模块构成抗干扰处理单元。ad-da处理芯片内置混频处理单元，混频处理单元配置频率为1568mhz的单载波混频信号，对多阵元天线接收到的gnss射频信号进行混频处理,混频后的bds-b1信号频段可以覆盖gps-l1的频段。

在计算自相关矩阵之前，信号与数据处理模块对i、q数据流进行2ms时间长度的缓存，其中缓存数据中前1ms用于数据抽取并计算自相关矩阵并进一步计算权值系数，后1ms缓存数据等待。这种设计可以保证在dbf合成时，根据缓存的i、q采样数据的输出时刻，将dbf权值系数与之进行时间对齐进行复数相乘。

gnss导航卫星信号抗干扰具体处理流程包括：阵列天线gnss信号接收，gnss信号采样与下变频处理，中频数据缓存与抽取，自相关矩阵计算，数字波束合成权值矩阵计算，波束合成以及信号上变频处理。

步骤1：gnss信号采样与下变频处理。采用四阵元天线接收gnss射频信号，并对四路天线信号进行独立处理，以某一天线阵元的处理流程为例进行说明，其他三路处理完全一致。

参阅图2。ad-da处理芯片内置混频处理单元，混频处理单元配置频率为1568mhz的单载波混频信号，对多阵元天线接收到的gnss射频信号进行混频处理。经过混频后，gnss信号中gps-l1(频段：1575.42±1.023mhz)信号混频后的中心频点为7.42mhz±1.023mhz，bds-b1(频段：1561.098±2.046mhz)信号混频后的中心频点为6.902mhz±2.046mhz，混频后的bds-b1信号频段可以覆盖gps-l1的频段。ad-da处理芯片使用62mhz模数转换a/d进行数字采样，得到数字中频采样信号并送入信号与数据处理模块；信号与数据处理模块配置数控振荡器nco和频率为6.902mhz的直接数字频率合成器(dds)产生的正、余弦信号，对数字中频采样信号进行数字下变频，并设计带宽为±2.046mhz的有限长单位冲激响应(fir)低通滤波器进行滤波处理，从而获得零中频且带宽为±2.046mhz，并包含gps-l1和bds-b1两个频点的i路和q路信号，之后用于dbf权值系数的计算。

上述信号采样与下变频处理方法的优点在于，通过对一种频率信号的处理，可以同时获得gps-l1与bds-b1两个频点信号的抗干扰能力。

参阅图3。步骤2：中频数据抽取与缓存：i、q路信号的62mhz数字采样在1ms内的采样点个数为62000，采用最大深度缓存队列缓存2ms，124000组数据，以进行数据流控制，数据流中前1ms，62000组数据用于数据抽取，后1ms缓存数据等待。在1ms的62000组采样数据中，抽取的间隔点数为62，共抽取1000组i、q采样数据，用于dbf权值系数中自相关矩阵的计算。通过上述最大缓存队列的设计，可以实现权值计算完成后可以与计算权值的采样数据的时间同步适配性处理。中频数据抽取与缓存包含但不仅限于2ms数据缓存，还包括其他时间长度的数据缓存。

[步骤3：自相关矩阵计算

在自相关矩阵计算中，信号与数据处理模块对选取的1000组采样数据和其共轭转置的乘积求解统计平均，得到如下自相关矩阵r:

x(n)＝[x1(n),x2(n),x3(n),x4(n)]^t，对于第i行第j列的元素的自相关矩阵其计算公式为：

式中，表示阵列天线中阵元j所对应的第n个i、q采样数据的共轭，n为采样数据组数，n为采样数据序号，t表示矩阵转置运算，xi(n)表示阵列线中第i个阵元所对应的第n个i、q采样数据。

步骤4：dbf权值系数计算

信号与数据处理模块根据功率倒置(pi)准则计算最优权值系数wpi：然后对上述最优权值系数wpi计算公式进行简化得到：

其中，r^-1为阵列输入信号自相关矩阵的逆矩阵，b为信号的导向矢量，在这里根据pi准则，b＝[1,0,0,0]^t。

步骤5：信号与数据处理模块在dbf合成中，根据缓存的i、q采样数据的输出时刻，将dbf权值系数与之进行时间对齐，实现dbf权值系数，与自相关矩阵求解时所用i、q采样数据的匹配。

信号与数据处理模块将dbf权值系数与对应的i、q采样数据进行复数相乘，获得对应于m采样时刻的dbf合成处理的i、q数据：

其中，xi(m)表示i天线第m个i、q采样数据，为wpi第i行系数的共轭，m＝1,2,…,62000。

步骤6：上变频

上变频处理分为数字和模拟两个部分，数字上变频部分采用6.902mhz的dds与dbf合成后的i、q数据进行复数相乘，再经过fir滤波器处理，获得6.902mhz为中心频率且±2.046mhz为带宽的数字信号，之后，将该数字信号输入ad-da处理芯片，进行模拟上变频处理，混合1568mhz的单载波信号后，进行功率匹配设置，处理后的gnss射频信号传输给普通接收机使用。根据以上步骤的实施，从而实现对gps-l1与bds-b1信号的抗干扰处理。

本领域技术人员将认识到，对以上描述做出众多变通修改是可能的，所以本实施例仅是用来描述一个或多个特定的实施方式。

尽管已经描述和叙述了被看作本发明的实施范例，本领域技术人员将会明白，可以对其做出各种变通，而不会脱离本发明的精髓。另外，也可以做出许多修改以将特定情况适配到本发明的教义，而不会脱离在此描述的本发明中心概念。所以，本发明不受限于在此披露的特定实施例，但本发明可能还包括属于本发明范围的所有实施例及其等同物。