基于MUSIC算法的测向性能分析

基于MUSIC算法的测向性能仿真2013年1月16日大作业--基于MUSIC算法的测向性能仿真I摘要随着移动通信技术的飞速发展，智能天线技术研究的不断深入，来波方向(DOA)估计技术逐渐成为研究的热点之一，而MUSIC算法是智能天线技术的典型算法。本文在对MUSIC算法进行分析的基础上，设计了MUSIC算法的仿真程序，对不同情况下该算法的性能进行了仿真分析。仿真结果表明该算法在不同阵列结构、信号入射角度时具有不同的性能。关键词：智能天线；DOA；MUSIC；阵元大作业--基于MUSIC算法的测向性能仿真II目录摘要...............................................................................................................................................I引言..............................................................................................................................................1一、MUSIC算法介绍.....................................................................................................................11.1MUSIC算法的提出...........................................................................................................11.2波达方向估计问题中的阵列信号数学模型.....................................................................21.3阵列协方差矩阵的特征分解.............................................................................................41.4MUSIC算法的原理及实现...............................................................................................61.5MUSIC算法的实现步骤：...............................................................................................8二、MUSIC算法的DOA估计仿真..............................................................................................82.1MUSIC算法的基本仿真....................................................................................................82.2MUSIC算法DOA估计与阵元数的关系.........................................................................92.3MUSIC算法DOA估计与阵元间距的关系...................................................................102.4MUSIC算法DOA估计与快拍数的关系.......................................................................112.5MUSIC算法DOA估计与信噪比的关系.......................................................................122.6MUSIC算法DOA估计与信号入射角度差的关系.......................................................13三、MUSIC算法性能分析小结...................................................................................................15参考文献.................................................................................................................................15附录.....................................................................................................................................16附录一：MUSIC算法的基本仿真源代码...........................................................................16附录二：MUSIC算法DOA估计与不同阵元数关系仿真源代码.....................................17附录三：MUSIC算法DOA估计与阵元间距的关系仿真源代码...................................18附录四：MUSIC算法DOA估计与快拍数的关系仿真源代码.........................................21附录五：MUSIC算法DOA估计与信噪比的关系仿真源代码.......................................22附录六：MUSIC算法DOA估计与信号入射角度差的关系仿真源代码.........................24图目录图1-1等距线阵与远场信号..................................................................................................2图2-1MUSIC算法的DOA估计谱......................................................................................9图2-2阵元数不同时MUSIC算法的DOA估计谱............................................................10图2-3阵元间距不同时MUSIC算法的DOA估计谱.......................................................11图2-4快拍数不同时MUSIC算法的DOA估计谱...........................................................12图2-5信噪比不同时MUSIC算法的DOA估计谱...........................................................13图2-6角度间隔不同时MUSIC算法的DOA估计谱.......................................................14大作业--基于MUSIC算法的测向性能仿真1引言智能天线技术是当前无线移动通信领域颇为关注和研究的热点领域之一，可将无线电的信号导向到具体的方向上，产生空间定向波束，使天线主波束对准用户信号到达方向，旁瓣或零陷对准干扰信号的到达方向，起到充分高效利用移动用户信号并删除或抑制干扰信号的目的。而波束形成的关键是要准确知道信号的到达方向，即波达方向，所以波达角估计(DOA)是波束形成的基础。本文着重分析了用于DOA估计的典型算法--MUSIC(MultipleSignalClassification)算法，然后对不同的条件下MUSIC算法的性能进行了Matlab的仿真和分析。一、MUSIC算法介绍1.1MUSIC算法的提出多重信号分类（MUSIC）算法是Schmidt等人在1979年提出的。这一算法的提出开创了空间谱估计算法研究的新时代，促进了特征结构类算法的兴起和发展，该算法已成为空间谱估计理论体系中的标志性算法。此算法提出之前的有关算法都是针对阵列接收数据协方差矩阵进行直接处理，而MUSIC算法的基本思想则是对任意阵列输出数据的协方差矩阵进行特征分解，从而得到与信号分类相对应的信号子空间和与信号分量相正交的噪声子空间，然后利用这两个子空间的正交性构造空间谱函数，通过谱峰搜索，检测信号的DOA。正是由于MUSIC算法在特定的条件下具有很高的分辨力、估计精度及稳定性，从而吸引了大量的学者对其进行深入的研究和分析。总的来说，它用于阵列的波达方向估计有以下一些突出的优点：(1)多信号同时测向能力(2)高精度测向(3)对天线波束内的信号的高分辨测向(4)可适用于短数据情况大作业--基于MUSIC算法的测向性能仿真2(5)采用高速处理技术后可实现实时处理1.2波达方向估计问题中的阵列信号数学模型为了分析推导的方便，现将波达方向估计问题中的数学模型作理想状态的假设如下：(1)各待测信号源具有相同的极化、且互不相关的。一般考虑信号源为窄带的，且各信号源具有相同的中心频率0。待测信号源的个数为D。(2)天线阵列是由M(MD)个阵元组成的等间距直线阵，各阵元特性相同，各向同性，阵元间隔为d，并且阵元间隔不大于最高频率信号半波长。(3)天线阵列处于各信号源的远场中，即天线阵列接收从各信号源传来的信号为平面波。(4)各阵元上有互不相关，与各待测信号也不相关，方差为2的零均值高斯白噪声)(tnm。(5)各接收支路具有完全相同的特性。图1-1等距线阵与远场信号设由第k（k=1，2，…D）个信号源辐射到天线阵列的波前信号为)(tSk，前k1d23M大作业--基于MUSIC算法的测向性能仿真3面已假设)(tSk为窄带信号，则)(tSk可以表示为以下形式：}exp{)()(tjtstSkkk（1.1）式中)(tsk是)(tSk的复包络，k是信号)(tSk的角频率。前面已经假设D个信号具有相同的中心频率，所以有：ck20（1.2）式中c是电磁波波速，是公用的信号波长。设电磁波通过天线阵列尺寸所需的时间为1t，则根据窄带假设，有如下近似：)()(1tSttSkk（1.3）故延迟后的波前信号为：)](exp[)()](exp[)()(~101011ttjtSttjttSttSkkk（1.4）所以，若以第一个阵元为参考点，则t时刻等间距直线阵中的第m(m=1，2，…M)个阵元对第k个信号源的感应信号为：]sin2)1(exp[)(kkkdmjtSa（1.5）其中，ka为第m个阵元对第k个信号源的影响，前面以假设各阵元无方向性，所以可取1ka。k为第k个信号源的方位角，cdmksin)1(表示由第m个阵元与第1个阵元间的波程差所引起的信号相位差。计及测量噪声和所有信号源来波，第m个阵元的输出信号为：)(]sin2)1(exp[)()(1tndmjtStxmDkkkm