本发明属于卫星导航抗干扰领域，涉及一种简化的卫星导航干扰信号测向方法。
背景技术：
在卫星导航抗干扰信号处理过程中，确定干扰信号的空间位置，即波达方向(doa)估计是一个基本问题。最早的超分辨doa估计方法是著名的music(多重信号分类：multiplesignalclassification)方法。
music方法属于特征结构的子空间方法，它建立在这样一个基本观察之上：若传感器个数比信源个数多，则阵列数据的信号分量一定位于一个低质的子空间，在一定条件下，这个子空间将唯一确定信号的波达方向，并且可以使用数值稳定的奇异值分解精确确定波达方向。由于在频域中进行doa估计，music算法对波达方向估计的性能改进是明显的。但是，该方法需要进行fft变换，运算量大，占用硬件资源多，硬件集成化程度低。
技术实现要素：
为了克服现有技术的不足，本发明提供一种改进的卫星导航干扰测向方法，能够简化测向过程，易于工程实现。
本发明解决其技术问题所采用的技术方案包括以下步骤：
第一步，对单天线阵元输出矩阵进行单比特量化，单个阵元在n时刻接收到的射频数字信号
第二步，对均匀线阵的单比特量化阵列输出矩阵进行fft变换，得到
第三步，测频并获取频域样本，对阵元1的单比特量化阵列fft后的结果搜索谱峰，然后记录谱峰值和相应的频点坐标，其他m-1个阵元取出该频点坐标位置的频域复数值，对应m个阵元得到m个峰值，称为频域峰值快拍矢量x(f)，该矢量是一个m×p的矩阵，其中m是阵元个数，p是信源个数，f为信源的频点位置标记；
第四步，对频域峰值快拍矢量x(f)的快拍数进行统计平均，得到m×m维的自相关估计
第五步，对自相关估计
第六步，搜索方向矢量到噪声子空间投影的峰值，得到波达方向的估计；搜索函数
本发明的有益效果是：采用单比特方法，使得数据位宽降低；由于数据位宽少，只有1位，极大的简化了fft变换过程，省去了乘法运算，只剩下加法运算，运算量大大降低，使得算法实现时占用的硬件资源减少70％以上，可以在一个芯片上实现整个设计，极大的简化了测向过程，易于工程实现。
附图说明
图1是本发明的方法流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施例。
本发明的算法流程如图1所示，实现步骤如下：
第一步：对单天线阵元输出矩阵进行单比特量化，量化过程如下式：
其中，x(n)为单个阵元在n时刻接收到的射频数字信号；x(t)为射频模拟信号。
第二步：对均匀线阵的单比特量化阵列输出矩阵x(n)进行fft变换，变换公式如下：
其中n是fft变换的采样点数；
第三步：频域检测，测频并获取频域样本。对阵元1的单比特量化阵列fft后的结果搜索谱峰，然后记录谱峰值和相应的频点坐标，其他m-1个阵元取出该频点坐标位置的频域复数值，对应m个阵元得到m个峰值，称为频域峰值快拍矢量x(f)，该矢量是一个m×p的矩阵，其中m是阵元个数，p是信源个数，f为信源的频点位置标记。
第四步：对频域峰值快拍矢量x(f)的快拍数进行统计平均，得到自相关估计，计算表达式如下：
其中，p为信源个数；上标h为共轭转置运算符；
第五步：对自相关估计
第六步：搜索方向矢量到噪声子空间投影的峰值，得到波达方向的估计。搜索函数表达式如下：
其中a(w)为方向矢量，g为噪声子空间矩阵，上标h为共轭转置运算符。
本发明适用于gps、bds和glonass卫星导航系统干扰来向检测。以8单元均匀线阵检测2个不同方向干扰源来向为例阐述本发明的具体实施方式。
步骤1：对8单元均匀线阵各阵元输出矩阵依次进行单比特量化，数据采样快拍数n＝1024，计算过程如下：
步骤2：对均匀线阵的单比特量化阵列输出矩阵x(n)进行fft变换，变换公式如下：
步骤3：频域检测，测频并获取频域样本。对阵元1的单比特量化阵列fft后的结果搜索谱峰，然后记录谱峰值和相应的频点坐标，其他7个阵元根据频点坐标取出相应位置的频域复数值，对应8个阵元得到8个峰值，称为频域峰值快拍矢量x(f)，该矢量是一个8×2的矩阵。
步骤4：对频域峰值快拍矢量x(f)进行统计平均，得到自相关估计，计算表达式如下：
其中，2是信源个数；上标h为共轭转置运算符；
步骤5：依次对自相关估计
步骤6：搜索方向矢量到噪声子空间投影的峰值，得到波达方向的估计。搜索函数表达式如下：
其中a(w)为方向矢量，g为噪声子空间矩阵，上标h为共轭转置运算符。
至此，完成对8单元均匀线阵检测2个不同方向干扰源来向检测。