基于改进的MUSIC算法的波达方向估计方法技术

本发明专利技术属于通信技术领域，公开了基于改进的MUSIC算法的波达方向估计方法。该方法包括：建立对称阵列模型，对各天线阵元接收的射频信号进行采样，正交下变频和滤波得到复信号，对复信号构成的阵列输出矩阵做线性变换得到复矩阵，进而处理得到实对称矩阵，对实对称矩阵进行特征值分解，构建噪声子空间矩阵，构造阵列方向矩阵，利用噪声子空间和阵列方向矩阵构造空间谱函数，进而进行谱峰局部极值搜索，得到函数极大值对应的角度值即为波达方向。本发明专利技术能够在不损失算法性能的情况下将特征值分解由复矩阵运算转换为实矩阵运算，并利用实方向矩阵进行峰值搜索，从而能有效地估计出信号的波达方向，提高了波达方向估计的准确性并且降低了计算的复杂度。

【技术实现步骤摘要】
基于改进的MUSIC算法的波达方向估计方法
本专利技术涉及通信
，尤其涉及基于改进的MUSIC算法的波达方向估计方法。
技术介绍
在无线电通信、导航、雷达、声纳、电子侦察等领域，获得信号的来波方向极其重要。在军事无线通信领域，获得高分辨率的来波方向不仅有助于对我方目标的导航、跟踪，还可以有助于我方对敌方目标进行准确定位以及打击。在民用无线通信领域，准确的来波方向估计能够提升系统容量和减少多址干扰等。空间谱表示信号在各个方向上的能量分布，若能够估计出空间谱就可以得到信号源的波达方向，所以空间谱估计也常称为波达方向(英文全称：DirectionOfArrival，英文缩写：DOA)估计，DOA估计技术的关键在于利用空间不同位置的天线阵列，接收来自不同方向的多个信号源发出的信号，运用现代信号处理技术估计出信号源的来波方向。在多种空间谱估计算法中，多重信号分类(英文全称：MultipleSignalClassification，英文缩写：MUSIC)算法是一种基于子空间分解的算法，其特点为测向误差小、精度高、分辨率高、灵敏度高，可以对入射信号数目、到达方向及波形的强度做出无偏估计。因此研究MUSIC算法对于促进更广泛的应用有着十分重要的意义。但是，传统的MUSIC算法在协方差矩阵计算及特征分解时涉及大量的复数运算，导致算法的复杂度高，不利于硬件实现。因此，如何在降低计算复杂度的同时又不损失算法性能就显得尤为重要。
技术实现思路
本专利技术的实施例提供基于改进的MUSIC算法的波达方向估计方法，能够在不损失算法性能的情况下将特征值分解由复矩阵运算转换为实矩阵运算，并利用实方向矩阵进行峰值搜索，从而能有效地估计出信号的波达方向。此专利技术提高了波达方向估计的准确性并且降低了计算的复杂度。为达到上述目的，本专利技术的实施例采用如下技术方案：步骤1、建立对称阵列模型，天线阵列为均匀线阵或者平面均匀圆阵，阵元个数为偶数；当天线阵元为均匀线阵时，以阵列的对称轴为基准，向右将阵列对称轴右侧天线阵元依次编号为向左将阵列对称轴左侧天线阵元依次编号为当天线阵元为平面均匀圆阵时，以其中任一阵元为起点，按顺时针方向依次编号为1，2...M，阵元间隔为d；步骤2、对各天线阵元接收的射频模拟信号进行K次采样，其中，第m个阵元天线的第k次采样得到的采样信号为sm(k)，m＝1，2，…M，k＝1，2，…K；K为快拍数，K＞2fs/Δf，fs为采样频率，Δf为需要分辨的最小频率，M为天线阵元个数，M为偶数；步骤3、对各天线阵元的K个采样信号做I、Q正交下变频和滤波，得到各天线阵元的K个复信号，并对第一阵元的K个复信号做快速傅里叶变换，判断快速傅里叶变换的最大值是否超过预设的门限值，如果是，则执行步骤4；步骤4、利用各天线阵元的K个复信号构成阵列输出矩阵X，对阵列输出矩阵X做线性变换得到复矩阵Y，利用复矩阵Y进一步计算得到复共轭对称矩阵R，利用复共轭对称矩阵R进一步计算得到实对称矩阵RX；步骤5、对实对称矩阵RX进行特征值分解，得到实对称矩阵RX的特征值和对应的特征向量，并将特征值按照从大到小的顺序进行排序；步骤6、基于最小描述长度MDL准则，利用排序后的实对称矩阵RX的特征值估计信号源个数，进而根据信号源个数，利用实对称矩阵RX的特征向量得到噪声子空间矩阵；步骤7、基于对称阵列模型构造阵列方向矩阵，利用噪声子空间和阵列方向矩阵构造空间谱函数，进而利用空间谱函数进行谱峰局部极值搜索，得到空间谱函数极大值对应的角度值，角度值即为信号源的入射方向。本专利技术的实施例提供基于改进的MUSIC算法的波达方向估计方法，首先建立天线的对称阵列模型，通过对中心对称阵列结构的数学分析，利用采样数据及其共轭信息得到了一种线性变换方法，能够在不损失算法性能的情况下将特征值分解由复矩阵运算转换为实矩阵运算，并利用实方向矩阵进行峰值搜索，从而能有效地估计出信号的波达方向。此专利技术提高了波达方向估计的准确性并且降低了计算的复杂度。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术实施例提供的基于改进的MUSIC算法的波达方向估计方法的流程示意图；图2为四阵元均匀直线阵的阵列结构示意图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。图1为本专利技术实施例提供的基于改进的MUSIC算法的波达方向估计方法的流程示意图，如图1所示，步骤1、建立对称阵列模型，天线阵列为均匀线阵或者平面均匀圆阵，阵元个数为偶数；当天线阵元为均匀线阵时，以阵列的对称轴为基准，向右将阵列对称轴右侧天线阵元依次编号为向左将阵列对称轴左侧天线阵元依次编号为当天线阵元为平面均匀圆阵时，以其中任一阵元为起点，按顺时针方向依次编号为1，2...M，阵元间隔为d；步骤2、对各天线阵元接收的射频模拟信号进行K次采样，其中，第m个阵元天线的第k次采样得到的采样信号为sm(k)，m＝1，2，…M，k＝1，2，…K；K为快拍数，K＞2fs/Δf，fs为采样频率，Δf为需要分辨的最小频率，M为天线阵元个数，M为偶数；步骤3、对各天线阵元的K个采样信号做I、Q正交下变频和滤波，得到各天线阵元的K个复信号，并对第一阵元的K个复信号做快速傅里叶变换，判断快速傅里叶变换的最大值是否超过预设的门限值，如果是，则执行步骤4；进一步的，利用公式对第一阵元的K个复信号做快速傅里叶变换；其中，第m个天线阵元第k次采样对应的复信号为xm(k)＝Im(k)+jQm(k)，m＝1，2，…M，k＝1，2，…K；N为快速傅里叶变换的点数，步骤4、利用各天线阵元的K个复信号构成阵列输出矩阵X，对阵列输出矩阵X做线性变换得到复矩阵Y，利用复矩阵Y进一步计算得到复共轭对称矩阵R，利用复共轭对称矩阵R进一步计算得到实对称矩阵RX；进一步的，步骤4具体包括：(4.1)利用各天线阵元的K个复信号构成阵列输出矩阵X，X＝[X1(k)，X2(k)，...，Xm(k)，...，XM(k)]T，Xm(k)＝[xm(1)，xm(2)，...xm(k)...，xm(K)]，m＝1，2，…M；(4.2)构造线性变换矩阵根据线性变换矩阵T和阵列输出矩阵X，计算得到复矩阵Y＝TX；其中，M为表示阵元个数，M为偶数；(4.3)利用复矩阵Y，计算得到复共轭对称矩阵R＝Y×YH，其中，上标H表示矩阵的共轭转置；(4.4)取复共轭对称矩阵R的实部，得到实对称矩阵RX＝Re{R}，其中，Re{·}表示取矩阵的实部。步骤5、对实对称矩阵RX进行特征值分解，得到实对称矩阵RX的特征值和对应的特征向量，并将特征值按照从大到小的顺序进行排序；进一步的，步骤5具体包括：(5.1)利用预设公式：RX＝UAUH，对实对称矩阵RX进行特征值分解；其中，Λ为对角矩阵，其对角本文档来自技高网...

【技术保护点】
1.一种基于改进的MUSIC算法的波达方向估计方法，其特征在于，包括以下步骤：步骤1、建立对称阵列模型，天线阵列为均匀线阵或者平面均匀圆阵，阵元个数为偶数；当天线阵元为均匀线阵时，以阵列的对称轴为基准，向右将阵列对称轴右侧天线阵元依次编号为

【技术特征摘要】
1.一种基于改进的MUSIC算法的波达方向估计方法，其特征在于，包括以下步骤：步骤1、建立对称阵列模型，天线阵列为均匀线阵或者平面均匀圆阵，阵元个数为偶数；当天线阵元为均匀线阵时，以阵列的对称轴为基准，向右将阵列对称轴右侧天线阵元依次编号为向左将阵列对称轴左侧天线阵元依次编号为当天线阵元为平面均匀圆阵时，以其中任一阵元为起点，按顺时针方向依次编号为1，2...M，阵元间隔为d；步骤2、对各天线阵元接收的射频模拟信号进行K次采样，其中，第m个阵元天线的第k次采样得到的采样信号为sm(k)，m＝1，2，…M，k＝1，2，…K；K为快拍数，K＞2fs/Δf，fs为采样频率，Δf为需要分辨的最小频率，M为天线阵元个数，M为偶数；步骤3、对所述各天线阵元的K个采样信号做I、Q正交下变频和滤波，得到所述各天线阵元的K个复信号，并对第一阵元的K个复信号做快速傅里叶变换，判断所述快速傅里叶变换的最大值是否超过预设的门限值，如果是，则执行步骤4；步骤4、利用所述各天线阵元的K个复信号构成阵列输出矩阵X，对阵列输出矩阵X做线性变换得到复矩阵Y，利用复矩阵Y进一步计算得到复共轭对称矩阵R，利用复共轭对称矩阵R进一步计算得到实对称矩阵RX；步骤5、对实对称矩阵RX进行特征值分解，得到实对称矩阵RX的特征值和对应的特征向量，并将所述特征值按照从大到小的顺序进行排序；步骤6、基于最小描述长度MDL准则，利用排序后的实对称矩阵RX的特征值估计信号源个数，进而根据信号源个数，利用实对称矩阵RX的特征向量得到噪声子空间矩阵；步骤7、基于对称阵列模型构造阵列方向矩阵，利用噪声子空间和阵列方向矩阵构造空间谱函数，进而利用空间谱函数进行谱峰局部极值搜索，得到所述空间谱函数极大值对应的角度值，所述角度值即为所述信号源的入射方向。2.根据权利要求1所述的方法，其特征在于，在步骤3中，利用公式对第一阵元的K个复信号做快速傅里叶变换；其中，第m个天线阵元第k次采样对应的复信号为xm(k)＝Im(k)+jQm(k)，m＝1，2，…M，k＝1，...

【专利技术属性】
技术研发人员：李晓婷，胡山峰，王君方，张宝强，张凡，李瑞，赵微，张忠波，刘洋，
申请(专利权)人：西安烽火电子科技有限责任公司，
类型：发明
国别省市：陕西,61