一种无人机阵列幅相误差与信号DOA联合估计方法

摘要

本发明公开了一种无人机阵列幅相误差与信号DOA联合估计方法，包括：在每架无人机上搭载单个天线阵元，多架无人机组成的无人机群构成接收阵列以接收信源信号；在观测基线不变的情况下，通过无人机运动使阵列结构发生改变，每次阵列结构变化后对信源信号进行再次采集，得到多个信源信号；针对每个采集到的信源信号，计算其协方差矩阵，通过特征值分解得到相应的噪声子空间；通过噪声子空间和方向向量重构二次优化问题，构造代价函数，通过谱峰搜索得到幅相误差和DOA联合估计。本发明能够实现DOA和阵列幅相误差联合估计，校正幅相误差，提高无源定位的精准性。

说明书

技术领域

本发明涉及阵列信号处理技术领域，具体而言涉及一种基于无人机群体阵列的波达方向（Direction of Arrival, DOA）和幅相误差联合估计方法。

背景技术

阵列信号处理作为现代信号处理的一个重要分支，是近几十年来发展十分迅速的一个技术领域，在雷达、声呐、无线通讯等军事及民用领域得到广泛的应用。波达方向估计是阵列信号处理技术的关键技术之一，其主要目的是对空间信源位置进行估计。DOA估计技术发展迅速，相关的理论和技术仍处于不断完善之中。

阵列误差的普遍存在是空间谱估计技术应用到实际工程中困难的一个重要原因。一般情况下，几乎所有的DOA估计算法都是建立在精确已知阵列流型的前提上，为了得到良好的算法估计效果，必须保证实际使用的阵列和理论研究中的标准阵列模型完全一致。但在实际应用中无论器件自身的因素还是实际环境因素都可能导致阵列出现误差。当使用理想的阵列流型进行实际的DOA估计时，就不可避免地得到误差较大的测向结果，或者测向结果根本就是无效的。大多数阵列误差的影响最终都能归结为阵列幅相误差。因此，研究阵列幅相误差存在条件下的DOA估计算法，对于空间谱估计技术的实用化具有重要意义。

发明内容

本发明针对现有技术中的不足，提供一种无人机阵列幅相误差与信号DOA联合估计方法，基于无人机群实现阵列位置变化，能够实现DOA和阵列幅相误差联合估计，校正幅相误差，提高无源定位的精准性。

为实现上述目的，本发明采用以下技术方案：

一种无人机阵列幅相误差与信号DOA联合估计方法，所述方法包括以下步骤：

S1，在每架无人机上搭载单个天线阵元，多架无人机组成的无人机群构成接收阵列以接收信源信号；

S2，在观测基线不变的情况下，通过无人机运动使阵列结构发生改变，每次阵列结构变化后对信源信号进行再次采集，得到多个信源信号；

S3，针对每个采集到的信源信号，计算其协方差矩阵，通过特征值分解得到相应的噪声子空间；

S4，通过噪声子空间和方向向量重构二次优化问题，构造代价函数，通过谱峰搜索得到幅相误差和DOA联合估计。

为优化上述技术方案，采取的具体措施还包括：

进一步地，步骤S1中，在每架无人机上搭载单个天线阵元，多架无人机组成的无人机群构成接收阵列以接收信源信号的过程包括以下步骤：

S11，使M架无人机均匀排列，每架无人机上搭载一个天线阵元，阵元间距为单位 间隔

S12，假设有K个平行的平面波从

其中，

式中，

进一步地，步骤S3中，针对每个采集到的信源信号，计算其协方差矩阵，通过特征值分解得到相应的噪声子空间的过程包括以下步骤：

S21，在基线不变的情况下，通过无人机的运动，使相应的阵元位置发生改变以构成新的阵列，对于每个新构成的阵列，采集相应的信源信号；

S22，对于第i个采集的信源信号

式中，L表示数据的快拍数；

S23，对

其中，

进一步地，步骤S4中，通过噪声子空间和方向向量重构二次优化问题，构造代价函数，通过谱峰搜索得到幅相误差和DOA联合估计的过程包括以下步骤：

S41，构建二次优化问题：

其中，

式中，

S42，构造代价函数：

式中，

S43，对

其中，

S44，获得角度与幅相误差的估计值：

S45，将

进一步地，所述方法还包括以下步骤：

采用均方根误差作为性能估计指标对估计结果的有效性进行评估；根据下述公式计算得到相应的均方根误差：

其中，N表示蒙特卡洛仿真次数，

本发明的有益效果是：

与现有技术相比，本发明突破了现有技术中DOA估计对幅相误差的局限，能够获得精确的角度估计值，具有更准确的定位性能；在幅相误差存在的情况下，本发明无需辅助校准信源、无需辅助校准阵元、无需迭代求解就可以估计出幅相误差的值并进行校正，能获得高分辨率估计。

附图说明

图1是本发明的无人机阵列幅相误差与信号DOA联合估计方法流程图。

图2是本发明的无人机群构成的接收阵列场景图。

图3是未校正幅相误差与使用本发明方法校正后的谱峰比较图。

图4是本发明在不同信噪比下的DOA估计性能比较图。

图5是本发明在不同信噪比下的幅相误差实部与虚部估计性能比较图。

具体实施方式

现在结合附图对本发明作进一步详细的说明。

需要注意的是，发明中所引用的如“上”、“下”、“左”、“右”、“前”、“后”等的用语，亦仅为便于叙述的明了，而非用以限定本发明可实施的范围，其相对关系的改变或调整，在无实质变更技术内容下，当亦视为本发明可实施的范畴。

为了便于描述，本实施例中的符号表示含义如下：

图1是本发明的无人机阵列幅相误差与信号DOA联合估计方法流程图。该方法用多架无人机作为平台搭载阵元，多个阵元构成一个阵列接收信号。该方法包括以下步骤：

S1，在每架无人机上搭载单个天线阵元，多架无人机组成的无人机群构成接收阵列以接收信源信号。

S2，在观测基线不变的情况下，通过无人机运动使阵列结构发生改变，每次阵列结构变化后对信源信号进行再次采集，得到多个信源信号。

S3，针对每个采集到的信源信号，计算其协方差矩阵，通过特征值分解得到相应的噪声子空间。

S4，通过噪声子空间和方向向量重构二次优化问题，构造代价函数，通过谱峰搜索得到幅相误差和DOA联合估计。

具体实现如下：

步骤1：信号接收

其中

式中，

根据数据模型可以得到接收信号信息，计算协方差矩阵：

式中，L表示数据的快拍数。对

其中，

步骤2：获得多个噪声子空间

在基线不变的情况下(即信号入射角不变)，通过无人机的运动，相应的阵元位置 发生改变，构成了一个新的阵列。此时再次接收信号并进行同样的处理，可以得到噪声子空 间

步骤3：DOA与幅相误差联合估计

存在幅相误差时，MUSIC函数变为：

其中

令

其中

对

将

为验证本发明所述算法的有效性，下面通过MATLAB仿真分析进行证明，性能估计指标为均方根误差（root mean square error, RMSE），定义为：

其中N表示蒙特卡洛仿真次数，

如图2所示为本发明的场景图，仿真中选取的无人机数目为M=6，P=2，

如图3所示为幅相误差在利用本发明方法校正前后的谱峰对比图。仿真中假设信噪比SNR设置为20dB，快拍数设置为J=500。从谱峰图可以看出幅相误差会影响空间谱谱峰的位置和幅度大小。当幅相误差存在时MUSIC算法估计性能将下降，如果误差较大甚至会导致算法失效。

图4为本发明方法与未校正方法在不同信噪比下的DOA估计性能比较图。从仿真结果可以看出，本发明方法具有很高的估计精度，并且可以通过适当提高信噪比的方法来提高本发明所述的基于无人机群体阵列的联合幅相误差与信号DOA估计方法的估计精度。

图5为本发明方法在不同信噪比下的幅相误差估计性能图。从仿真结果可以看出本发明方法可以准确地估计出实部与虚部。

本发明公开了一种基于无人机群体阵列的联合幅相误差与信号DOA估计方法，属于阵列信号处理技术领域。本发明通过每架无人机搭载单个天线阵元，由多架无人机组成的无人机群构成了接收阵列，通过无人机群内个体位置的改变，对应阵元位置也会发生变化，由此变换阵列结构。对多次阵列接收到的信号计算协方差矩阵进行特征值分解，得到多个信号噪声子空间。通过噪声子空间和方向向量重构二次优化问题，构造代价函数，最后通过谱峰搜索确定DOA，并同时获得幅相误差估计。本发明突破了传统的无人机群协同感知中DOA估计精度受限于无人机间幅相误差的藩篱，且无需辅助信源和阵元，也无需迭代求解，可以获得高精度的DOA和幅相误差联合估计，具有重要的应用价值。

以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。