一种平面阵列多目标角度高分辨实现方法技术

本发明专利技术公开了一种平面阵列多目标角度高分辨实现方法，该方法包含如下步骤：S1，采用二维简化为一维的方式，分别对俯仰方向、方位方向阵列接收数据做波束域FFT，得到俯仰方向、方位方向目标大致方向；S2，根据得到的方位、俯仰角度进行二维角度遍历配对，并确定角度搜索区间；S3，对多目标二维角度的估计，确定目标个数和角度。本发明专利技术可以有效的估计多目标的角度，同时算法运算量远小于传统的平面阵列角分辨算法，有利于工程实现。

【技术实现步骤摘要】

本专利技术涉及相控阵雷达领域，特别涉及一种平面阵列多目标角度高分辨实现方 法。
技术介绍
波达方向（D0A)估计一直是阵列信号处理的重要研究内容之一，在通信、声纳、雷 达、海洋开发、生物医学工程等应用领域中占有十分重要的地位。 空间谱估计算法可突破阵列的角度分辨"瑞利限"限制，提高空间角度分辨能力。 以MUSIC算法为代表的子空间类算法由于其良好的分辨性能，受到研究者的广泛关注，是 D0A估计中最常用的方法。 在实际工程中，受条件限制，快拍数少、阵列误差的存在，将导致D0A估计算法性 能严重恶化；且传统的谱估计类算法当强弱信号同时存在时且相距较近时（强信号如主瓣 有源干扰），强信号会造成其附近的弱信号谱峰偏移甚至消失；另外实际工程中还需考虑 算法的实时实现，传统谱估计算法的全空域遍历谱峰搜索，很难满足实时性的要求。 现阶段出现了不少针对阵列误差、小快拍数，以及强弱信号、实时实现的空间谱估 计相关研究。但相关文献中的算法大部分都是针对其中一项或二项开展的研究，不利于工 程实现。本文综合考虑算法性能及工程实现性，提出了一种平面阵列多目标角度高分辨实 现方法。 侯颖妮等人在文献中针对均匀线阵提出一种通过通道ToePlitz化来实现幅相误 差条件下的D0A估计的方法，其避免了参数优化的庞大运算量。但未考虑二维阵列应用及 实时实现问题。李新潮等利用对角加载技术加强算法在快拍数少的情况下的稳定性，但该 方法需要根据实验来确定对角加载量的大小，且未考虑误差存在情况。Li Jian等提出的 Relax方法，采用减的思想消除强信号对弱信号D0A估计的影响，每估计一个新角度就重 新修整之前所有角度的估计，算法计算量较大；文献利用乘的思想抑制强信号，设计阻塞矩 阵，阻塞矩阵与强信号导向矢量正交，作用在于阻止强功率信号通过，进而接收信号经过阻 塞矩阵输出后不存在强信号分量。文献中利用投影变换将阵列接收数据投影到信源大致方 向区间，一方面通过对接收数据的预处理来抑制模型误差，增强算法的稳健性；另一方面使 数据维数降低算法运算量，但是该方法在方向区间内构建投影变换矩阵及求投影矩阵会带 来新的大运算量，且对于二维阵面和多个目标的情况下算法实时性更差。采用预先存储的 方式提高实时性需遍历所有的二维角度，将对系统的硬件存储带来很大的压力。
技术实现思路
本专利技术的目的是提供，快拍数少、阵列 误差以及强弱信号相近的情况下，可以有效的估计多目标的角度，同时算法运算量远小于 传统的平面阵列角分辨算法，有利于工程实现。 为了实现以上目的，本专利技术是通过以下技术方案实现的： -种平面阵列多目标角度高分辨实现方法，其特点是，该方法包含如下步骤： S1，采用二维简化为一维的方式，分别对俯仰方向、方位方向阵列接收数据做波束 域FFT，得到俯仰方向、方位方向目标大致方向； S2,根据得到的方位、俯仰角度进行二维角度遍历配对，并确定角度搜索区间； S3,对多目标二维角度的估计，确定目标个数和角度。 所述的步骤S1具体为：采用二维简化为一维的方式，先对俯仰方向阵列接收数据 做波束域FFT，得到俯仰方向目标大致方向Θpi= 1，2. . .I，再对方位方向阵列接收数据 做波束域FFT，得到方位方向目标大致方向βj，j= 1，2, ...J其中I，J分别为俯仰方向和 方位方向初步确定的信号个数。 所述的步骤S2具体为： S2. 1，判断方位方向和俯仰方向得到的角度是否均为1，若是，则执行步骤 S2. 1. 1，若否，则执行步骤S2. 2; S2. 1. 1，计算平面阵列方位角，，其中，W=arcsin(siii/?/co#>、且二维角度搜索 区域设为（θ± 93dB，);S2. 2,对方位、俯仰角度遍历配对，得到相应的二维阵列方位角，每个二维阵列方 位角采用如下公式：P=arcsin(siii/?/cos6?)，同时对每个配对后的角度设定相应的二维角 度搜索区域（θρ±Θ3dB，(pp±q)3dB)，p= 1，2...，P,为配对后角度编号。 所述的步骤S3具体包含： 步骤S3. 1，通过特征值分解得到L个大特征值对应的特征向量θι，i= 1，2,. ..，L; 并得到噪声子空间UN= : 式中ε为小值常数；x(t)为阵列接收到的信号；#为协方差矩阵；λλ2，…λΜ 为特征值分解得到的特征向量％，e2，…eM为相应的特征向量;Μ为二维阵列阵元个数； 步骤S3. 2,计算搜索空间角度区域每个空间二维角度上方向性矢量在噪声子空间 上的投影值： 式中，UN= [eL+1，eL+2, ."ej为噪声子空间，糾，％>为角度0，心对 应的二维方向性矢量。 步骤S3. 3,计算搜索空间角度区域每个角度上方向性矢量在每个大特征值对应的 特征向量上的投影值： 式中，θι，e2，…^为大特征值对应的特征向量；为角度他％)对应的二维 方向性矢量； 步骤S3. 4,通过峰值搜索确定目标的个数和角度： 式中，fs(n，i，j)为空间角度区域每个角度上方向性矢量在每个大特征值对应的 特征向量上的投影值；fN(i，j)为搜索空间角度区域每个空间二维角度上方向性矢量在噪 声子空间上的投影值。 本专利技术与现有技术相比，具有以下优点： 1、本专利技术将二维阵列简化为方位、俯仰一维，通过波束域FFT得到搜索的大致角 度区域，然后再利用整个阵面的接收数据来进一步获得更准确的目标个数和角度。采用 此方式可大大降低谱估计算法的角度遍历搜索运算量，提高算法的实时性。当目标个数 为2个时，阵列方位、俯仰3dB波束宽度均为3°时，角度搜索运算量可降低到传统算法的 4/225。 2、本专利技术提出改进的基于投影理论的二维MUSIC算法实现多目标二维角度的估 计，确定目标个数和角度，该方法利用大特征值对应的特征向量与信号之间的相关性以及 信号与噪声空间的正交性，提高算法对小快拍数、阵列误差、强弱信号相近情况下算法的有 效性。【附图说明】 图1为本专利技术一种平面相控雷达天线阵元排列示意图； 图2为本专利技术流程图； 图3为传统谱估计算法空间谱图； 图4A为最大特征值对应特征向量的空间谱图；图4B为第二大特值对应特征向量的空间谱图； 图4C为第三大特值对应特征向量的空间谱图。【具体实施方式】 以下结合附图，通过详细说明一个较佳的具体实施例，对本专利技术做进一步阐述。 如图2所示，，该方法包含如下步骤： S1，采用二维简化为一维的方式，分别对俯仰方向、方位方向阵列接收数据做波束 域FFT，得到俯仰方向、方位方向目标大致方向； S2,根据得到的方位、俯当前第1页1 2 本文档来自技高网...

【技术保护点】
一种平面阵列多目标角度高分辨实现方法，其特征在于，该方法包含如下步骤：S1，采用二维简化为一维的方式，分别对俯仰方向、方位方向阵列接收数据做波束域FFT，得到俯仰方向、方位方向目标大致方向；S2，根据得到的方位、俯仰角度进行二维角度遍历配对，并确定角度搜索区间；S3，对多目标二维角度的估计，确定目标个数和角度。

【技术特征摘要】

【专利技术属性】
技术研发人员：黄飞，李毓琦，王佳，朱剑，
申请(专利权)人：上海无线电设备研究所，
类型：发明
国别省市：上海;31