摘要
如今，随着卫星导航技术的飞速发展，位置信息服务已经融入到我们的日常生活中，导航目前被称为继移动互联网后第三大产业。卫星导航在维护国家的安全中也发挥着不可替代的作用。为了使导航系统不受干扰的影响，本文以北斗导航系统为平台，研究基于阵列天线的自适应抗干扰算法。
首先，文章就自适应抗干扰算法的原理和方法进行了系统介绍，并在 MATLAB 中建立阵列模型，对基于功率倒置算法的空域抗干扰算法和空时联合抗干扰算法进行性能仿真。然后根据系统的指标，确定了在 FPGA 中实现抗干扰算法的方案，包括数字下变频、权值计算、数据加权、数字上变频等模块。根据权值计算模块实现方式的不同，本文提供了两种抗干扰算法在 FPGA 中实现的方案：一种是基于 FPGA 嵌入式软核 NIOS II 的抗干扰实现，将权值计算的过程放在 NIOS II 软核中，用 C 语言进行实现；另一种是基于逻辑语言的抗干扰算法的实现，即用硬件描述语言 Verilog HDL 进行权值的计算。权值计算涉及到浮点数运算和 Hermite 矩阵求逆，本文给出了各模块的设计方法和仿真结果，并与 MATLAB 仿真结果进行对比。最后给出了两种实现方案的实测结果，表明两种实现方案都具备抗干扰性能，验证了抗干扰系统的正确性。
**关键词：**空域滤波，空时自适应处理，矩阵求逆，浮点运算，DDC，FPGA，NIOS II
6 北斗导航系统抗干扰算法实测验证
本章主要对空域和空时两种抗干扰算法进行实测验证。权值计算 FPGA 实现的过程中，两种算法的区别主要在于协方差矩阵的维数上。空域算法协方差矩阵维数低，采用分块求矩阵逆的方法，占用资源少，权值更新时间短，选用低成本、低功耗的 CYCLONE III 系列作为实现平台即可。空时算法协方差矩阵为高阶矩阵，分块求矩阵逆的方法不再适合。CYCLONE V 系列提供了实数矩阵求逆的 IP 核，在该 IP 核的基础上可以快速求出复数共轭对称矩阵的逆。
6.1 抗干扰算法的硬件实现平台和实验场景的搭建
图 6.1 是系统的数字板结构框图。数字板主要由一片 FPGA 芯片、一片时钟芯片、一片 4 通道 ADC 芯片、一片配置芯片、一片 DA 芯片组成。为了实现低功耗、小型化的设计要求，ADC 芯片选用 AD9253。ADC 采样后的信号传送到 FPGA。在 FPGA 芯片中完成数字下变频、空域滤波（空时二维联合）算法、数据加权，然后再数字上变频，交给射频通道进行模拟上变频后，最后把数据交给接收机进行解调。
空域算法的 FPGA 实验板主要是由 ALTERA 公司的 EP3C120F780 芯片和 ADI 公司的 AD9268 芯片组成，空时算法的 FPGA 实验板主要是由 ALTERA 公司的 5CGXFC7C6F23I7 芯片和 ADI 公司的 AD9253 芯片组成，两种算法使用的数字板的实物图如图 6.2 所示。
a) 空域算法数字板 b) 空时算法数字板
图 6.2 抗干扰算法的硬件实现平台
图 6.2（a）中 1 是 FPGA 芯片，2、3 是 AD 芯片，4 是时钟分配芯片，5 是 FLASH，6、7 是两个 DSP 芯片（本文没有用到），8 是 DA 芯片。数字板的尺寸为 95mm95mm。图 6.2（b）中 1 是 FPGA 芯片，2 是 AD 芯片，3 是 DA 芯片，5 是晶振，数字板的尺寸为 58mm45mm。可见少了 DSP 芯片后整个数字板的功耗、尺寸以及成本都大大降低。整个实验环境的搭建与仪器使用情况如图 6.3 所示。
图 6.3 北斗导航抗干扰系统实验环境
图 6.3 中，1 是抗干扰天线，2、3、5 是 3 个方位角俯仰角均不同的干扰源天线，4 是二维转台，6 是信号源，7 是便携式频谱仪，8 是接收机，9 是直流电源，10 是电脑，11 为八角天线。整个系统的工作流程为：
（1）搭建图 6.3 所示的实验环境。将 1 北斗导航抗干扰天线置于 4 二维转台上，可以测试天线各个方向上的抗干扰性能，并模拟干扰源运动的情况下抗干扰天线能否正常搜星定位。由于是 4 单元天线阵列，因此抗干扰天线最多能在 3 个干扰源下正常工作，将 3 个干扰源架在 3 个方位角不同的位置，同时调节干扰源的俯仰角角度。将 6 信号源与 3 个干扰天线相连作为干扰源，用 11 八角天线标定干扰源到达天线阵面的干扰信号功率的大小并调节干扰源到适当位置。将 1 导航天线与 8 接收机相连，并将接收机通过串口线 RS232 连接到 10 电脑，用搜星软件进行搜星定位。也可以用 7 频谱仪观察抗干扰天线输出的信号频谱，判断干扰是否被抑制。待整个环境搭建好之后，开启 9 直流电源，整个系统开始工作。
（2）调节 6 信号源的输出信号，包括信号的大小、类型（点频、宽带等）。
（3）启动转台，将天线 360 度旋转。
（4）在 10PC 端观察北斗导航抗干扰天线在有干扰下的搜星定位是否正常。
实验的判断指标主要包括：
（1）输出频谱是否稳定。
（2）搜星数量是否大于等于 4 个，是否会发生丢星现象。
（3）通过 signaltap 抓取抗干扰前的数据和抗干扰后的数据，在 MATLAB 中分析 FPGA 抗干扰的效果。
（4）将 signaltap 抓取的干扰数据在 MATLAB 中进行仿真，对比 MATLAB 与 FPGA 抗干扰之后的频谱图。
6.2.1 模数转换器 ADC（AD9253）
本系统选取的 ADC 是串行 LVDS 输出的 AD9253。AD9253 单个芯片包括四个通道，且采用的是 48 脚的 LFCSP 封装，大大简化了 PCB 布线工作，同时减少了在数字板上的使用面积。AD9253 的采样率可以达到 125MSPS，而且支持双倍速率（DDR）操作。AD9253 采用串行 LVDS 电平信号数据输出，拥有 14bit 分辨率，数据读取操作简单，接口时序如图 6.4 所示。其中 CLK 为差分采样时钟对，DCO 是输出比特数据同步时钟，FCO 为字节标志位，标志一个字节数据的开始。AD9253 的数据不是高低位依次输出的，因此需要在 FPGA 中根据标志位将数据重新组合，完成串并转换功能。AD9253 是 14 位分辨率，输出 16 位的时候，默认最低两位为 0。
