基于L阵列的目标波达方向估计算法改进(共25页).doc

精选优质文档-----倾情为你奉上基于L型阵列的目标波达方向估计算法改进摘 要：从目标波达方向估计的工程应用角度分析，普通的阵列布局已不能满足现代人们的需求本文主要阐述了一种采用次最小冗余阵的L型阵列的波达方向估计，与传统L型阵列相比，采用次最小冗余阵能够减少阵元的使用同时提高角度分辨率本文采用9个阵元的L型阵列，通过采用次最小冗余阵能够实现只使用9个阵元就能达到19个阵元的效果实验表明，结合传统的MUSIC算法，在进行DOA判断时比普通的L型阵列的分辨率具有明显的提高关键词：目标波达方向；次最小冗余阵；阵元；MUSIC专心---专注---专业Research on directionofarrival estimation algorithm based on L - type arrayAbstract：From the perspective of engineering application of reach direction estimation, the general array layout can no longer meet the needs of modern people. In this paper, a method for estimating the direction of arrival of wave in an l-type array with subminimum redundancy matrix is presented. In this paper, the l-type array with 9 elements is adopted, and the effect of 19 elements can be achieved by using only 9 elements with sub-minimal redundant matrix. The experimental results show that the resolution of DOA with traditional MUSIC algorithm is significantly higher than that of ordinary l-type arrayKey words: directionofarrival; Subminimum redundant matrix;Pixelarray; MUSIC目 录1 引言1.1 研究背景及意义技术不断革新的当代生活，人们对信息化以及智能化的设备需求也在不断提升，表现形式也愈加丰富多彩，波达方向（Direction of arrival, DOA）估计主要用于判断信源方向的一种手段。

随着近年来关于阵列信号理论和研究的不断发展与进步，加之越来越多的地方应用到智能化的处理，DOA技术又被研究学者们推上研究的重心阵列信号也就是将一或多组的阵元（传感器）按照一定的规律进行排布组合，构成一个阵列，形成多个阵元接收信号的系统通过系统不断的在空间内进行对信号的搜寻以及发射接收等动作，将有用的信号保留并放大，对无用的信号或者噪声等信号进行抑制甚至消除，将目标信号中的特征信息也就是有用信息进行提取[1]对于传统技术而言，DOA估计算法就必须是符合空间理论的满阵，也就是平常所说的阵元之间的距离必须小于等于入射信号的半波长，以此用来避免产生测量时产生的模糊误差但是在现代所有产品都趋向微型化和智能化的社会，且在现实生活中，各种干扰因素过多，信号复杂且交叠，参数范围变化幅度愈加不明显，这对传统的满阵DOA方法充满了挑战，导致传统满阵的DOA系统显然已经逐渐被淘汰，所以如何实现空间方位识别率更准，鲁棒性更高，实时性更快的DOA估计方法成为全球相关科研人员及研究机构研究的重点内容[2,3]通过查阅资料了解到，阵列的孔径直接对阵列对空间的空间分别率起着重要的作用，且成正比关系，也就是孔径越大则系统的空间分辨率越低。

同时，阵元增加不仅增加了硬件成本，对多通道信号的的处理也提出了更高的要求，不适合工程应用[4]稀疏阵列指的是阵元之间的距离不小于入射信源半波长的阵列排布和传统的满阵相比，稀疏阵列孔径变得更大，自由度提高，在方向测量精准度、空间分辨率和多信源监测时拥有更好的表现本系统采用最小冗余阵和MUSIC算法结合的方法来提高系统分辨率1.2 研究历史和现状在初期的研究中，最常用的估计方向的方法现在被称为波束形成算法，由于多种因素的原因，该算法在很多应用环境下受限严重，虽然现在技术已经很发达了，但是还是很难在该算法上做出新的突破，影响最严重的是因为该算法受到瑞利的制约，然而根据现代科研学者的研究，发现通过增大天线的孔径可以实现对精度的提高，但是如何合理的增大天线孔径又成了一个难题[5]知道20世纪60年代，关于如何提高分辨率的问题再次出现在人们的视野，其中以最大熵法和最小方差法在这方面应用最为突出，但是这两种算法的局限性也很强，使用这两种算法在对线性函数进行预测时准确率很高[6]，但针对非线性预测时却效果很不理想而且由于在实际进行DOA估计时，条件非常苛刻，理想环境下的情况几乎不存在，两种线性算法在这时使用的效果很不理想，所以在DOA估计方面很大程度限制了这两种线性算法的发展[7]。

后续出现了更多的优秀算例，例如有些研究学者将时频分布和小波分析加入进了超分辨算法当中综上所述，窄带信号不是该领域的局限，而当代研究学者的研究方向更倾向于再实际情况下使用，因此对宽带信号得研究逐渐被挖掘出来[9,10]信源的数量估计有着不可替代的关键地位，在目前已知的算法中其大部分需要事先知道入射信号源的数量，如果对信源数量的估计不正确的话，不仅是虚警和漏报的问题，同时也会影响子空间之间的正交特性信源数量估计的错误还会对子空间和信号维度的判断造成影响，致使对方位进行判断时产生的角度差过大在这类算法的初期阶段，通过增加惩罚函数的方法避免了主观门限的弊端上述准则的局限性也非常大，仅仅适应在白噪声和非相干信号的前提下，而且这种情况对快拍数的要求也非常的高，然而后面通过Wu等人针对于盖氏圆盘定理的信源数目估算方法，该算法在对信源数量进行估计时不需用对协方差矩阵的特征值，而是利用噪声和信号圆盘半径差距的特性通过上述分析可以知道，超分辨测向技术存在这非常多的缺点，最知名的缺点就是该技术对环境要求非常高，而且对阵列矢量的精度要求非常高，必须在没有任何误差的情况下，才能最大限度的保持其精准度，这样就限制了其在实际生活中的应用，因为在正常的实际使用中，实际的阵列形状会受到多种外界环境的影响，会对测量的结果形成偏差和扰动，从而对方向的判定效果和性能造成极大影响。

因此，在实际的生活应用中，阵列误差必须被考虑在内，并且必须要通过手段来对其进行调整，确保其始终保持在一定的范围内，如何保证让误差始终保持在一定的范围内，这个问题一直到上世纪90年代末才被研究的一种新方法解决，该方法通过将阵列校正通过一种参数估计的方法，使得其波动范围始终在一定的范围内，这样有效的解决了上述的问题，而且，使用该方法硬件成本低，设备兼容性高，同时还能够针对不同的环境进行自适应，从而将该项技术又向前推进了一步2 DOA基本概念2.1 DOA估计阵列数学模型在进行DOA空间谱估的数学建模时，可在理想状态下记性假设，其假设如下：⑴ 每个待检测的信源极化相同且互不相关，通常我们认为信源为窄带信号，且每个信源有着一样的中心频率，假设待检测的信源数量为D⑵ 天线阵列看成是由M（M>D）个阵元构成的同距线阵，阵元特性相同、各项同性，设d为阵元之间的距离，根据阵元间距规则阵元之间的距离不能大于最高频率波长信号的一半⑶ 在每个远场的信源中必须要放置天线，因为对方向的估计就是通过信源中天线的波动来实现断定的⑷ 对于除了主接收线路外，其他各支路在理论和实际应用中的特征完全一样图1为DOA估计的一般模型图。

图 21 DOA估计的一般模型假使使用第个信源用于对天线阵列的辐射，其中辐射信号为，此时我们令作为这次辐射信号中的窄带信号，这样我们可以将表示为 通过一系列的变化可以得到如式2.2所示： 式中，为光速；为信号波长在我们常规的测量中使用来表示经过天线阵列时电磁波所需用的时间，同时使用窄带假设理论对其进行计算，得到如下的结果： 因此延迟后的波前信号如下所示： 因此，假设选择的参考点为第一个阵元，那么在t时刻同距线阵中的第个阵元能够对第k个信号源进行感应，感应的信号为 通过上式，我们可以知道是第k个信源受第m个阵元影响产生的上面已经对阵元方向进行了假设，即所有阵元均无方向，所以当时第k个信源的方位角我们使用来进行表示式2.6是我们通过计算得到的第m个阵元信号的输出方程，想要得到该方程就必须要对所有的声源进行监测，通过要去除所采集的声源中噪音和无用的部分 式中，为测量噪声标号即代表着第多少个阵元设 通过矩阵运算，我们可以变换出更加简单明了的表达式： 式中，在一般的信号处理中，一段真实的语音信号可以分解成若干信号的叠加，各种各样的声音通过叠加形成了我们采集到的声音，在这些叠加的声音中，有我们需要的特征声音，当然也包含了许多杂乱的无用信号，如何去除这些无用信号也是语音信号处理的关键要解决的问题，在DOA估计上，我们可以通过对X进行阵列输出的相关处理，从而得到 式中，表示矩阵共轭转置。

2.2 DOA估计方法波束形成和电子引导作为DOA估计的传统方法的理论概念，这两种方法并非是单一的对信号进行接收，然后通过信号的特征进行识别然后建模而实现的估计，在这里我们需要用这两种方法对信号进行采集，与此同时对阵列进行一定程度上的引导，通过原始的DOA估计方法结合现有的技术电子引导，这样可以使我们在进行DOA估计时任意的调节我们需要的方向，并且能够针对峰值信息进行寻找下面我们将对传统的DOA估计方法进行描述和讨论2.2.1 延迟—相加法延时-相加法其采用的是最经典的波束形成和傅里叶变换的方法图2给出了经典窄带波束形成器的结构图 22 窄带波束形成器构造输出信号y(t)可以看做是传感器阵元的输出的线性加权和，即 式中，中针对对个方面的因素进行了包含和容纳原始的波束形成器的总输出功率可以用下式进行表示 式中， 是关于的导引向量；是阵列输入端的噪声向量；和分别是信号和噪声的功率可以看出，当时，输出功率最大2.2.2 Capon最小方差法Capon最小方差法的初衷是为了解决延时法针对分辨率过低这一缺点而提出的一种新的方法，这种方法既然被提取出来使因为其算法的优越性主要体现在仅使用一个波束就可对多方面的问题进行监测和估计，使用这种方法可以将输出功率降低，达到增益观测方向保持为零数的前提下降低非期望干扰的贡献，增益观测方向通常为1。

其约束条件为利用拉格朗日乘子，权向量解为 利用此方法，可以通过Capon空间谱得到输出功率关于DOA的函数，即 2.2.3 MUSIC算法MUSIC算法和其他算法相比起来具有很强的适应性，可以对一维和二维的阵列进行DOA估计，在知道阵列分布的情况下即可得 (2.19)根据上式，得到矩阵为满秩矩，非奇异，所以有逆存在所以上式通过变化可转换为，这种变化足以说明中噪声子空间和各个列矢量之间存在正交关系，所以有。