本发明专利技术公开了一种基于稀疏重构的互质阵列波达方向角估计方法，主要解决现有技术运算量大，需要先验知识多，识别信源数目少，无源定位估计误差大的问题，实现步骤是：天线接收机形成互质线阵；对空间信号采样获得观测数据；由观测数据得到虚拟阵元接收数据矢量；空域网格划分，构造超完备基；用虚拟阵元接收数据矢量与超完备基稀疏关系，建稀疏约束方程；采用凸优化法求解稀疏约束方程，得最稀疏解；以稀疏解与空间角度对应关系，绘制幅度谱图，获得波达方向角度值。本发明专利技术在低先验知识的情况下提高了无源测向的精度和运算速度，相同阵元数目，增加了阵列可识别的信源数目，提高了低信噪比下对信号方向角估计精度，可用于目标侦察和无源定位。
【技术实现步骤摘要】
本专利技术属于信号处理
，特别涉及声学信号、电磁信号的阵列信号处理技 术，具体是一种，可用于目标侦察与无源定 位。
技术介绍
信号的波达方向角DOA估计是阵列信号处理领域的一个重要分支，它是指利用天 线阵列对空间声学信号、电磁信号进行感应接收，再运用现代信号处理方法快速准确的估 计出信号源的方向，在雷达、声纳、无线通信等领域具有重要应用价值。随着科技的不断进 步，对信号波达方向估计的精确度和和分辨率也有越来越高的要求。 过去的DOA估计方法，通常使用的是典型线性均匀阵列。对于一个L阵元的典型 线性阵，传统的MUSIC类计算方法可检测的信源数目是L-1。为了在少的阵元条件下得到 尽量大的角度自由度，检测更多的信源，一些新的阵列结构被提出，比较有代表性的是嵌 套阵列以及互质阵列。对于这些阵列的DOA估计技术，通常的方法是平滑MUSIC子空间方 法（S-MUSIC)，但S-MUSIC方法需要较多的先验知识：需要预先知道信号数目，知道准确的 噪声分布方差，或者花费很多运算来估计噪声分布方差。近年来，一些新颖的观点：基于 Khatri-Rao积子空间协方差稀疏重构的DOA估计方法，也被应用到DOA估计中。但已有的 方法在计算互质阵列中虚拟阵列接收数据时需要计算自相关矩阵、互相关矩阵，之后抽取 排序，且所选用字典并非最优，从而不能估计最多的信源数目，在稀疏重构的过程中也没有 给出一个很好的误差分布参量，造成较大的D0A估计误差。而在实际应用中，信源数目和噪 声分布一般都是未知的，上述传统方法中存在运算量大和估计误差大等缺陷，会对目标侦 查、无源定位的速度和精度带来较大的影响：不能实时快速的检测目标，目标个数很多时甚 至无法识别，导致目标捕获失败。
技术实现思路
本专利技术的目的在于针对上述已有技术中需要已知信源数目和同等条件下识别信 源数量不足等问题，提出一种，在降低运算 量的情况下，提高目标侦察准确性，提高了无源定位在低信噪比、低快拍数、低先验知识条 件下的D0A估计精度，避免因角度估计误差引起的目标侦察失误，同时在一定阵元数条件 下增加了可识别的信源数目，在侦查目标多于阵元数目时表现尤为突出。 为实现上述目的，本专利技术的实现步骤包括如下： 1)用2M+N-1个天线接收机形成互质阵列，每个天线接收机称为一个阵元，并假设 有K个目标声学或电磁信号入射到该互质阵列，K多1。 la)用N个天线接收机形成均匀线性阵列1，简称线阵1，用2M-1个天线接收机形 成均匀线性阵列2,简称线阵2,线性阵列1的阵元间距为Md，线性阵列2的阵元间距为Nd， 其中N>M多2且M、N互质，0〈d<A/2,A为入射到互质阵列的窄带信号波长，定义线阵1 的第一个阵元为互质阵列的阵元0 ; lb)组合线性阵列1和线性阵列2为互质阵列：将线阵2的第一个阵元放置于阵元 0后相距为Nd的位置，线阵2的所有阵元依次插于线阵1中，线阵1和线阵2在同一条线上， 从互质阵列阵元0开始，由头至尾命名各个阵元依次为阵元0,阵元1，……，阵元2M+N-2。 2)由互质阵列天线接收机对空间目标的声学或电磁信号进行采样，得到阵列输出 信号Y(t)，定义Ym(t)为第m个阵元的输出信号，m= 0, 1，. . .，2M+N-2。 3)根据互质阵列输出信号Y(t)计算虚拟阵列接收数据矢量y。 4)对探测空域波达方向角观测空间进行网格划分，构造超完备基0(0)，并定义 一个空域稀疏向量s。 5)将探测空域波达方向角估计转化为求解如下稀疏约束方程：min||s| |!subjectto| |y-? ( 9 )s| 1n,s^0 其中s是一个QX1维的未知矢量，q是需要预先估计的误差分布参量，II?lli 是矩阵的一阶范数，11 ? 112是矩阵的二阶范数。6)采用凸优化方法求解稀疏约束方程，得到未知矢量s的最稀疏解5。 7)以波达方向角范围0 = 的值为x轴坐标，以5向 量的幅度值为y轴坐标，绘制幅度谱图，从该幅度谱图中按照从高到低的顺序寻找幅值较 大的前K个谱峰，这些谱峰的峰值点所对应的x轴坐标即为目标的波达方向角度值。 本专利技术与现有技术相比具有以下优点： 1)本专利技术将互质阵列应用到D0A估计中，将实际阵元接收信号转化为虚拟阵列接 收信号，通过选择无冗余虚拟接收信号，降低了信号数据量，减小了运算量，增加了阵元数 目一定条件下，阵列可识别的信源数目。 2)本专利技术采用稀疏表示技术将波达方向角的估计转化为稀疏信号的重构，是新理 论技术与传统问题的结合，利用入射信号源的空域稀疏特性进行建模，突破了阵列分辨率 的瑞利限，提高目标侦察和无源定位在低信噪比、低快拍数、低先验知识条件下的D0A估计 精度，避免因角度估计误差引起的目标侦察失误。 3)本专利技术利用信号源的空域稀疏性得到的稀疏解中大系数对应的角度，为信号源 的波达方向角，无目标的方向角对应的系数约等于〇,因此无需预先知道目标的数目；计算 虚拟阵列接收信号时采用阵列的互相关信息，噪声的自相关信息没有或者很少量的被带入 模型所选取的虚拟阵列接收信号，因此无需已知或估计噪声分布方差，在现实环境中具有 更广泛的应用价值。【附图说明】 图1是本专利技术的实现流程图； 图2是本专利技术互质阵列结构示意图； 图3是本专利技术中M= 2,N= 3时的虚拟阵元所有可能的位置集合信息图； 图4是本专利技术与现有一种传统互质阵列波达方向角估计方法的估计效果比较图； 图5是本专利技术对估计信源数目扩展能力的效果图。【具体实施方式】 以下参照附图，对本专利技术的技术方案和效果作进一步的详细说明。 实施例1 在军事侦察中，比如在雷达领域，探测者预先不知道目标的个数和波达方向角等 信息，而这些信息又尤为重要，为获取这些信息，通常的做法是先估计出目标的个数，然后 通过子空间等方法得到目标信号的空间谱，传统空间谱估计一般都是建立在目标个数已知 的情况下，而且传统方案识别的目标数目要少于阵元数，假如所要捕获的目标数目要多于 阵元数目，传统的做法是增加阵列接收机的阵元数目，这在一定程度上增加了成本或者是 现场无法实现，传统方法中运算量大和估计误差等缺陷将对目标侦查、无源定位的速度和 精度带来较大的影响：不能实时快速的检测目标，目标个数很多时无法识别，目标捕获失败 等。 为此，本专利技术经过创新研宄，提供了一种基于稀疏重构的互质阵列波达方向角估 计方法，参照图1，本专利技术的实现步骤如下： 1)用2M+N-1个天线接收机形成互质阵列，每个天线接收机称为一个阵元，并假设 有K个窄带信号入射到该互质阵列，1。 la)参见图2,每隔Md放置一个天线接收机，共放置N个天线接收机，形成均匀线 性阵列1，简称线阵1 ;再每隔Nd放置一个天线接收机，共放置2M-1个天线接收机，形成均 匀线性阵列2,简称线阵2,线性阵列1的阵元间距为Md，线性阵列2的阵元间距为Nd，其其 N>M多2且M、N互互，假设有个K远场窄带信号入射到该线阵上，且信号在传播过程中加入 了均值为〇的复高斯白噪声，〇〈d<X/2,X为入射到互质阵列的窄带信号波长，定义线阵 1的第一个阵元为互质阵列的阵元0。 lb)组合线性阵列1和线性阵列2为互质阵列：将线阵2的第一个阵元放置于阵 元0后相距为Nd本文档来自技高网...
【技术保护点】
一种基于稀疏重构的互质阵列波达方向角估计方法，其特征在于，包括以下步骤：1)用2M+N‑1个天线接收机形成互质阵列，每个天线接收机称为一个阵元，并假设有K个信号入射到该互质阵列，K≥1；1a)用N个天线接收机形成均匀线性阵列1，简称线阵1，用2M‑1个天线接收机形成均匀线性阵列2，简称线阵2，线性阵列1的阵元间距为Md，线性阵列2的阵元间距为Nd，其中N>M≥2且M、N互质，0<d≤λ/2，λ为入射到互质阵列的窄带信号波长，定义线阵1的第一个阵元为互质阵列的阵元0；1b)组合线性阵列1和线性阵列2为互质阵列：将线阵2的第一个阵元放置于与阵元0相距为Nd的位置，线阵2的所有阵元依次插于线阵1中，从互质阵列阵元0开始，由头至尾命名各个阵元依次为阵元0，阵元1，……，阵元2M+N‑2；2)由互质阵列天线接收机对空间目标信号进行采样，得到阵列输出信号Y(t)，定义Ym(t)为第m个阵元的输出信号，m＝0,1,...,2M+N‑2；3)根据互质阵列输出信号Y(t)计算虚拟阵列接收数据矢量y；4)对探测空域波达方向角观测空间进行网格划分，构造超完备基Φ(θ)，并定义一个空域稀疏向量s；5)将探测空域波达方向角估计转化为求解如下稀疏约束方程：min||s||1subject to ||y‑Φ(θ)s||2≤η,s≥0其中s是一个Q×1维的未知矢量，η是需要预先估计的误差分布参量，||·||1表示矩阵的一阶范数，||·||2表示矩阵的二阶范数；6)采用凸优化方法求解稀疏约束方程，得到未知矢量s的最稀疏解7)以波达方向角范围θ＝[θ1,θ2,...,θq,...,θQ]的值为x轴坐标，以向量的幅度值为y轴坐标，绘制幅度谱图，从该幅度谱图中按照从高到低的顺序寻找幅值较大的前K个谱峰，这些谱峰的峰值点所对应的x轴坐标即为目标的波达方向角度值。...
【技术特征摘要】
【专利技术属性】
技术研发人员：蔡晶晶，鲍丹，武斌，刘高高，秦国栋，李鹏，马亚东，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。