本发明涉及阵列信号处理、无源定位技术、doa估计、无人机故障检测与隔离领域，特别涉及一种基于无人机机载阵列的干扰源定位方法及系统。
背景技术：
1、无人机具有造价较低、飞行环境要求低、安全性高、便于携带等优点，因此无人机在定位、电力巡检、影视拍摄、农业喷洒等领域有着广泛的应用。近年来，电磁干扰事件经常发生，对卫星通信与导航造成干扰。针对日益增加的监测区域要求，地基设备灵活性较差，而无人机具有较强的灵活性，且观测视距的影响也远小于地面设备。因此基于无人机平台的定位具有一定现实意义。
2、到达角度(doa)估计在无源定位中有着至关重要的地位，经典的超分辨率doa估计算法主要有music算法，esprit算法，这些都是基于特征结构的子空间算法。与传统技术不同，子空间算法利用接收信号的特征结构，而不是利用收集数据的统计特征，突破了“瑞利限”的限制，这使得估计分辨力显著提高。许多天线阵列都可以用来进行doa估计，最流行的平面阵列是均匀矩形阵列(ura)，其中，元件间的间距需要在半波长以内以避免空间混叠。然而，传统ura的自由度受到传感器数量的严格限制。此外，还受到阵列元件之间严重的相互耦合效应的影响。
3、稀疏阵列可以有效解决上述问题。互质阵列是稀疏阵列的一种，和传统ura相比，互质阵列能够增加自由度，同时提高算法的估计性能。虽然与嵌套阵列相比，互质阵列的空间自由度较低，但是它减少了相互耦合的影响。
技术实现思路
1、为解决上述技术问题，本发明提供了一种基于无人机载天线阵列的doa角度定位方法，该方法在确保无人机位置准确的前提下，引入稀疏表示方法，建立凸优化模型对干扰源进行角度估计，实现高精度定位。
2、本发明还提供了一种基于无人机载天线阵列的doa角度定位系统。
3、为达到上述目的，本发明的技术方案如下：
4、一种基于无人机载天线阵列的卫星导航干扰源定位方法，包括：
5、步骤1，数个搭载有天线阵列的无人机按照固定轨迹运动，采集无人机在不同的飞行状态下的多种传感器数据信号，多种传感器数据信号包括三轴加速度，三轴角速度，滚转角、俯仰角、偏航角，对采集到的传感器数据信号进行去噪和归一化处理，提取出信号的有效特征；
6、步骤2，根据采集数据的种类和模式设置神经网络的参数，通过误差函数确定预测误差，当误差满足需要的精度时结束训练，利用训练好的网络进行无人机故障检测；
7、步骤3，无人机上搭载的天线阵列接收多个干扰信号，得到互质阵列接收信号模型；
8、步骤4，计算所接收的干扰信号的协方差矩阵并将其向量化；将机载互质阵列虚拟化得到虚拟差分阵列，分别计算虚拟阵列正负两部分的协方差矩阵并将其向量化；
9、步骤5，通过空间平滑构建原始阵列和虚拟阵列的综合协方差矩阵；
10、步骤6，利用稀疏表示，构造整个监测空间内的字典，进而构造多测向向量稀疏表示模型；
11、步骤7，将多测向向量稀疏表示模型转化为凸优化问题，构建二阶锥优化问题，通过谱峰搜索得到干扰源的doa角度估计。
12、根据本发明优选的，步骤1中，对采集到的传感器数据信号进行去噪和归一化处理，提取出信号的有效特征；包括：
13、使用小波阈值降噪对采集到的传感器数据信号抑制噪声，假设含有噪声的信号公式如式(1)所示：
14、f(t)＝x(t)+n(t)(1)
15、式(1)中，f(t)表示含有噪声的信号，x(t)表示目标信号，n(t)表示噪声信号；
16、选择软硬阈值折衷的方法分离出噪声，如式(2)所示：
17、
18、式(2)中，表示阈值化后的小波系数，wj,k表示原始的小波系数，sign表示符号函数，α表示变化因子，t表示阈值；
19、对信号进行归一化处理，归一化公式如式(4)所示：
20、
21、式(4)中，xi为输入的第i个值，xmin表示输入的最小值，xmax表示输入的最大值，表示归一化之后的结果。
22、根据本发明优选的，步骤2中，根据采集数据的种类和模式设置神经网络的参数，通过误差函数确定预测误差，当误差满足需要的精度时结束训练，利用训练好的网络进行无人机故障检测；包括：
23、bp神经网络包括输入层、隐含层和输出层；设x1、x2、…xn为输入，i为输入的节点数，y1、y2、…ym为输出，k为输出节点，隐含层为单隐层，则隐含层第j个节点的输出pj公式如式(5)所示：
24、
25、式(5)中，f1表示激活函数，wji表示输入层和隐含层之间的权值，aj表示隐含层的阈值；则输出层第k个节点的输出yk公式如式(6)所示：
26、
27、式(6)中，f2表示激活函数，wkj表示输入层和隐含层之间的权值，bk表示隐含层的阈值；设n为输入样本的数量，则网络误差公式即误差函数如式(7)、式(8)所示：
28、
29、
30、式(7)、式(8)中，p表示样本的数量，m表示输出节点的数量，dk表示期望输出的结果，yk表示实际输出的结果，ep表示第p个样本的正向传播结果与期望结果的误差，表示第p个样本的第k个输出节点的期望输出结果，表示第p个样本的第k个输出节点的实际输出结果，综合起来，误差函数ep衡量了每个样本的正向传播结果与期望输出之间的差异，而总误差函数e则是所有样本误差的累加；
31、输出层的权值变化量公式如式(9)所示：
32、
33、式(9)中，△wjk表示输出层的权值变化量，η表示学习率；输出层的阈值变化量△bk、隐含层的权值变化量△wji、隐含层的阈值变化量△aj公式分别如式(10)、式(11)、式(12)所示：
34、
35、
36、
37、权值和阈值的更新公式如式(13)、式(14)、式(15)、式(16)所示：
38、wji(n+1)＝wji(n)+△wji(13)
39、wkj(n+1)＝wkj(n)+△wkj(14)
40、aj(n+1)＝aj(n)+△aj(15)
41、bk(n+1)＝bk(n)+△bk(16)
42、式(13)、式(14)、式(15)、式(16)中，wji(n)、wkj(n)、aj(n)、bk(n)分别表示前面一层的权值和阈值，而wji(n+1)、wkj(n+1)、aj(n+1)、bk(n+1)表示更新之后的权值和阈值；
43、利用训练好的网络进行无人机故障检测：输入层的数据为一个9维的向量，代表x轴加速度、y轴加速度、z轴加速度、x轴角速度、y轴角速度、z轴角速度、滚转角、俯仰角以及偏航角的信号特征，采集到的传感器数据经过小波去噪得到去噪后的信号，再进行归一化处理作为训练好的网络的输入变量，输出层的数据为一个3维向量，代表三种状态：正常状态、慢变故障、突变故障。
44、根据本发明优选的，步骤3中，无人机上搭载的天线阵列接收多个干扰信号，并得到接收信号模型；包括：
45、空间中有k个干扰源sk(t)，k＝1,2,…,k，其来波方向为θ＝[θ1,θ2,…,θk]；空间中有数个无人机，机上各载有一个包含m＝2线阵和n＝3线阵的互质阵列；该互质阵列的传感器分布位置为{0,md,2md,…,(n-1)md}和{0,nd,2nd,…,(m-1)nd}，d＝λ2为单位阵元间距，λ是信号波长；两个线阵的第一个元素重叠，互质阵列的阵元数为m+n-1，互质阵列的位置集合如式(17)所示：
46、
47、每个互质阵列的两个线阵接收到的信号如式(18)、式(19)所示：
48、
49、
50、式(18)、式(19)中，x1(t)和x2(t)是在t时刻无人机上的互质阵列的两个线阵分别接收到的信号，τ1(θk)和τ2(θk)表示来波角度为θk的信号到达两个线阵的时间延迟，n1(t)和n2(t)为观测噪声，且该噪声均符合零均值高斯分布并和信号sk(i)不相关；
51、对x1(t)和x2(t)进行采样得到离散形式并作l点dft变换，其频域模型如式(20)、式(21)所示：
52、x1＝a1s+n1(20)
53、x2＝a2s+n2(21)
54、其中，a1＝[α1(θ1),α1(θ2),…,α1(θk)]；a2＝[α2(θ1),α2(θ2),…,α2(θk)]；s＝[s1(t),s2(t),…,sk(t)]t；n1＝[n1,0,n1,1,…n1,m-1]t；n2＝[n2,0,n2,1,…n2,n-1]t；a1和a2表示两个线阵的流形矩阵，和表示流形矩阵各列指向向量，s是信号源向量，n1和n2是两个线阵的噪声向量；将式(20)、式(21)中的两个信号合并为一个信号，得到互质阵列接收信号模型x；如式(22)所示：
55、
56、互质阵列接收信号模型x的阵列流形矩阵为
57、根据本发明优选的，计算所接收的干扰信号的协方差矩阵并将其向量化；将机载互质阵列虚拟化得到虚拟差分阵列，分别计算虚拟阵列正负两部分的协方差矩阵并将其向量化；包括：
58、计算互质阵列接收信号模型的协方差矩阵；如式(23)所示：
59、cxx＝e[xxh]＝acssah+w (23)
60、式(23)中，w＝e[nnh]＝σ2i；为信号功率，σ2为噪声功率，i是大小为(m+n-1)×(m+n-1)的单位矩阵；
61、将协方差矩阵向量化，得到向量化协方差矩阵，如式(24)所示：
62、v1＝d1p+μ(24)
63、式(24)中，d1＝a*⊙a，为新的阵列流形矩阵，⊙表示矩阵的khatri-rao积，μ为矩阵w的向量化形式；
64、式(24)看成一个阵元位置为集合{xi-xj,0≤i,j≤m+n-1}的虚拟线性阵列，其中，xi表示第i个阵元的位置，d1即为虚拟阵列流形矩阵；利用其稀疏性，对两线阵进行差分操作，获得虚拟阵列；虚拟阵列的阵元位置由集合给出；对虚拟阵列的正、负两部分分别作与信号模型相同的操作；即：先计算出正部分与负部分的协方差矩阵，再将协方差矩阵向量化；故对于正位置部分，如式(25)所示：
65、c2＝e[xxt]＝asstat+w (25)
66、将其向量化得到式(26)：
67、v2＝d2p2+μ (26)
68、其中，d2是正位置阵元的流形矩阵，信号功率为实值，故p2＝p；
69、再对虚拟阵列的负位置部分求协方差矩阵，如式(27)所示：
70、c3＝e[xxh]＝asshah+w (27)
71、将其向量化得到式(28)：
72、v3＝d3p3+μ (28)
73、其中，d3是负位置阵元的流形矩阵，p3＝p。
74、根据本发明优选的，通过空间平滑构建原始阵列和虚拟阵列的综合协方差矩阵；包括：将式(24)、式(26)、式(28)三个向量模型合并，如式(29)所示：
75、
76、相应的阵列流形矩阵为
77、采用空间平滑技术，首先，从流形矩阵d中抽取一个维度是(2mn+2m+2n-1)×k的矩阵，将虚拟阵列中的2mn+2m+2n-1个阵元分成mn+m+n个子阵，接收信号模型式(29)重构为式(30)：
78、
79、且信号模型的综合协方差矩阵为式(31)：
80、
81、其中，ω＝σ2e，ω＝σ2e，e是(2mn+2m+2n-1)×1的列向量，是(mn+m+n)×k的流形矩阵。
82、根据本发明优选的，步骤(6)中，利用稀疏表示，构造整个监测空间内的字典，进而构造多测向向量稀疏表示模型；包括：
83、基于阵列和式(31)的稀疏性，将观测的空域角度区间按照一定的角度间隔分成q个角度值；q＞＞k，定义集合如式(32)所示：
84、
85、式(32)，θ表示整个观测空域的搜索网格；表示整个观测空间的第q个网格的角度值，q＝1,2,…q，流形矩阵表示为
86、那么式(31)的第j列如式(33)所示：
87、
88、式(33)中，1≤j≤mn+m+n，gj∈rq×1是字典的稀疏表示系数；故将mn+m+n个第j列合并为一个矩阵形式，得到多测向向量稀疏表示模型，如式(34)所示：
89、
90、式(34)中，g＝[g1,g2,…,gmn+m+n]是一个稀疏系数矩阵。
91、根据本发明优选的，步骤(7)中，将多测向向量稀疏表示模型转化为凸优化问题，构建二阶锥优化问题，通过谱峰搜索得到干扰源的doa角度估计；包括：
92、定义γ(g)表示g中的非零行索引集合，|γ(g)|＝k表示非零行个数，即稀疏度，那么对应的doa估计问题表示为约束优化问题，如式(35)所示：
93、
94、定义为矩阵g的估计，其行之间具有稀疏性，列之间不具有稀疏性；定义b＝[b1,b2,…,bq]t，其每一个元素代表中对应行向量的2范数，即而b的1范数||b||1代表在空域θ上的稀疏约束；因此，优化问题更改为式(36)：
95、
96、噪声功率用sr的最小特征值作为其估计，定义为υ，故将式(36)正则化得到正则化模型，如式(37)所示：
97、
98、其中，是sr的估计，δ是正则化参数，代表frobenius范数；正则化模型为凸优化问题，采用二阶锥优化(socp)对其进行求解，如式(38)所示：
99、
100、式(38)中，辅助变量ε和η来线性化目标函数，i为q×1的单位向量，ξ为1×q的向量，ξq是ξ的第q个元素；通过式(38)求出稀疏向量b后，通过谱峰搜索即得到干扰源的doa角度估计；根据多个阵列的doa估计值，通过交叉定位得到目标源的位置。
101、一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现基于无人机载天线阵列的卫星导航干扰源定位方法的步骤。
102、一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现基于无人机载天线阵列的卫星导航干扰源定位方法的步骤。
103、一种基于无人机载天线阵列的卫星导航干扰源定位系统，包括：
104、信号有效特征提取模块，被配置为：数个载有天线阵列的无人机按照固定轨迹运动，采集无人机在不同的飞行状态下的多种传感器数据信号，对采集到的传感器数据信号进行去噪和归一化处理，从而尽可能的提取出信号的有效特征；
105、无人机故障检测模块，被配置为：根据采集数据的种类和模式设置神经网络的参数，通过误差函数确定预测误差，当误差满足需要的精度时结束训练，利用训练好的网络进行无人机故障检测；
106、接收信号模型获取模块，被配置为：无人机上搭载的天线阵列接收多个干扰信号，并得到接收信号模型；
107、多测向向量稀疏表示模型构造模块，被配置为：计算所接收的干扰信号的协方差矩阵并将其向量化；将机载互质阵列虚拟化得到虚拟差分阵列，分别计算虚拟阵列正负两部分的协方差矩阵并将其向量化；通过空间平滑构建原始阵列和虚拟阵列的综合协方差矩阵；利用稀疏表示，构造整个监测空间内的字典，进而构造多测向向量稀疏表示模型；
108、doa角度估计模块，被配置为：将多测向向量稀疏表示模型转化为凸优化问题，构建二阶锥优化问题，通过谱峰搜索得到干扰源的doa角度估计。
109、通过上述技术方案，本发明提供的一种基于机载天线阵列的卫星导航干扰源定位方法具有如下
110、有益效果：
111、1、本发明基于机载天线阵列的卫星导航干扰源定位方法，采用无人机机载阵列，灵活度高，定位前首先对无人机进行故障检测，若有无人机出现故障或者无人机位置不在预定位置则删除该无人机的接收信号，避免了接收无人机的位置误差引起的定位误差；
112、2、本发明基于机载天线阵列的卫星导航干扰源定位方法，采用的机载阵列为互质阵列，该阵列有较高的自由度，受阵元间的相互耦合的影响较低，且可以通过虚拟化扩展阵元，利用较少的阵元数量完成多干扰源估计；
113、3、本发明基于机载天线阵列的卫星导航干扰源定位方法，引入稀疏表示，构造凸优化模型，相比于music等传统doa估计算法，拥有更高的估计精度和更高的谱分辨率。
1.一种基于无人机载天线阵列的卫星导航干扰源定位方法，其特征在于，包括：
2.根据权利要求1所述的一种基于无人机载天线阵列的卫星导航干扰源定位方法，其特征在于，步骤1中，对采集到的传感器数据信号进行去噪和归一化处理，提取出信号的有效特征；包括：
3.根据权利要求1所述的一种基于无人机载天线阵列的卫星导航干扰源定位方法，其特征在于，步骤2中，根据采集数据的种类和模式设置神经网络的参数，通过误差函数确定预测误差，当误差满足需要的精度时结束训练，利用训练好的网络进行无人机故障检测；包括：
4.根据权利要求1所述的一种基于无人机载天线阵列的卫星导航干扰源定位方法，其特征在于，步骤3中，无人机上搭载的天线阵列接收多个干扰信号，并得到接收信号模型；包括：
5.根据权利要求4所述的一种基于无人机载天线阵列的卫星导航干扰源定位方法，其特征在于，通过空间平滑构建原始阵列和虚拟阵列的综合协方差矩阵；包括：将式(24)、式(26)、式(28)三个向量模型合并，如式(29)所示：
6.根据权利要求1所述的一种基于无人机载天线阵列的卫星导航干扰源定位方法，其特征在于，步骤(6)中，利用稀疏表示，构造整个监测空间内的字典，进而构造多测向向量稀疏表示模型；包括：
7.根据权利要求1所述的一种基于无人机载天线阵列的卫星导航干扰源定位方法，其特征在于，步骤(7)中，将多测向向量稀疏表示模型转化为凸优化问题，构建二阶锥优化问题，通过谱峰搜索得到干扰源的doa角度估计；包括：
8.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1-7任一所述的基于无人机载天线阵列的卫星导航干扰源定位方法的步骤。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-7任一所述的基于无人机载天线阵列的卫星导航干扰源定位方法的步骤。
10.一种基于无人机载天线阵列的卫星导航干扰源定位系统，其特征在于，包括：
