基于协方差矩阵重构的改进型MUSIC算法.pdf

第 32 卷第 3 期淮阴工学院学报Vol.32 No.32023年 6 月Journal of Huaiyin Institute of TechnologyJun.2023基于协方差矩阵重构的改进型MUSIC算法屠亚杰1，方遒1，2，李艳玲1（1.厦门理工学院 福建省客车先进设计与制造重点实验室，福建 厦门 361024；2.厦门大学 航空航天学院，福建 厦门 361102）摘要：针对经典MUSIC算法在信源相干条件下算法性能失效问题，提出一种基于数据协方差矩阵二次重构并对噪声子空间进行加权处理的改进型MUSIC算法。首先对接收到的数据协方差矩阵进行重构处理，然后根据特征值分解得到噪声特征值的最大值与最小值的比值关系，使用合适的加权系数得到新的噪声子空间并推导得到改进后的空间谱函数，最后结合谱峰搜索完成对信源目标的DOA估计。仿真结果表明，改进型MUSIC算法有效消除了信号的相干性，不造成阵列孔径损失，在低信噪比、小快拍数以及信号入射间隔较小的情况下，能有效估计出相干信号的波达方向。关键词：DOA估计；MUSIC算法；相干信号；协方差矩阵重构中图分类号：TN911文献标志码：A文章编号：1009-7961（2023）03-0036-06DOA Estimation of Improved MUSIC Algorithm Based on CovarianceMatrix ReconstructionTU Yajie1,FANG Qiu1,2,LI Yanling1（1.Fujian Key Laboratory of Bus Advanced Design and Manufacture,Xiamen University ofTechnology,Xiamen Fujian 361024,China;2.School of Aerospace Engineering,Xiamen University,Xiamen Fujian 361102,China）Abstract：An improved MUSIC algorithm based on quadratic reconstruction of data covariance matrixand weighting of noise subspaces is presented to overcome the performance failure of the classical MUSIC algorithm under the condition of source coherence.First,the covariance matrix of the received data is reconstructed,then the ratio between the maximum and minimum noise eigenvalues is obtainedby eigenvalue decomposition,a new noise subspace is obtained by using appropriate weighting coefficients,and an improved spatial spectral function is derived.Finally,the DOA estimation of the sourcetarget is completed by combining the peak search.The simulation results show that the improved MUSIC algorithm effectively eliminates the signal coherence,does not cause array aperture loss,and caneffectively estimate the direction of arrival of the coherent signal at low signal-to-noise ratio,smallnumber of snapshots and small signal incidence interval.Key words：DOA estimation;MUSIC algorithm;coherent signal;covariance matrix reconstruction收稿日期：2023-01-11基金项目：福建省自然科学基金项目（2022J011247）作者简介：屠亚杰（1997-），男，江苏淮安人，在读硕士，主要从事智能驾驶技术、车载毫米波雷达研究。通讯作者：方遒（1967-），男，福建漳州人，教授，硕士，主要从事汽车电子技术、新能源汽车研究。波达方向（DOA）估计作为阵列信号处理的研究内容之一，在雷达、航空、声纳、日常通信等诸多领域都有着广泛应用1-2。DOA估计的角度分辨率取决于阵列孔径，阵列孔径确定后，可以计算得到对应的角度分辨率，称之为瑞利限3。为突破瑞利限提出了一系列方法称为超分辨率方法，在众多超分辨率算法中多重信号分类4（Multiple Signal Classification，MUSIC）算法最为经典，它利用噪声特征向量和信号向量的正交关系得到阵列空间谱函数，在空域内进行谱峰搜索求出信源方向。基于MUSIC算法容易实第 3 期现、测角精度高等优点，学术界基于MUSIC算法展开了一系列研究。同非等3提出一种基于空间平滑的改进MUSIC算法，该算法解决MUSIC算法对相干信号DOA估计失效的难题，但会损失阵列孔径。石要武等5提出一种基于特征空间的MUSIC算法，该算法不仅可以对相干信号进行DOA估计，还可以对信号源的功率进行估计，提高对小能量目标信号的估计成功率，但该算法依然存在阵列孔径损失问题。刘晓志等6提出一种基于伪数据相关矩阵二次重构的DOA估计算法，该算法可以实现对相干信号的DOA估计且不损失阵列的孔径，能对强弱邻近信号作出准确估计，但在低信噪比环境下对相同信号强度的邻近相干信号的检测成功率下降。张贺勇等7提出一种基于空间平滑的单次快拍DOA估计算法，该算法应用空间平滑算法可以对相干和非相干信号进行比较准确的DOA估计，可以满足工程应用中对算法易于实现的需求，有着较高的应用价值，但该算法仅使用单快拍数据参与运算，在低信噪比环境下DOA估计结果与实际结果偏差大。针对上述方法的限制与不足，本文提出一种基于数据协方差矩阵二次重构并对噪声子空间进行加权处理的改进型MUSIC算法，通过重构数据协方差矩阵，对得到的噪声特征值进行加权处理，得到新的噪声子空间和改进后的空间谱函数，最后通过谱峰搜索可以实现相干信号入射角度邻近条件下的DOA估计，并对新算法的有效性进行仿真验证。1阵列模型和相干信号考虑有P个窄带远场目标信号从不同入射角=1,2,p入射到以下由M个接收阵元组成的均匀线阵，相邻阵元间的间距d=/2，为载波波长，信号与均匀线阵法线之间的夹角/2,/2。均匀线阵结构如图1所示。图1均匀线阵模型均匀线阵在t时刻接收到的数据向量为：y()t=y1()t,y2()t,yM()tT=x()t+n()t=As()t+n(t)（1）式（1）中A=a(1),a(p)为MP维方向矩阵，a()i=1,j2dsini,j2(M1)dsiniT为入射角度i方向的M1维导向向量，其中包含角度信息。s()t=s1()t,sP(t)T为P1维的入射目标信号向量，n()t=n1(t),n2(t),nM(t)T为M1维噪声向量。噪声服从零均值、方差为2的高斯分布，噪声向量相互独立，且与信号向量独立。T表示转置运算。对于阵列接收到的多个入射信号，入射信号之间的关联程度可以根据相关系数来划分为不相关、部分相关和相干3种情况。对于两个平稳信号Si(t)和Sj(t)，其相关系数定义为：ij=ESi()t ESi(t)Sj()t ESj(t)ESi()t ESi(t)2ESj()t ESj(t)2（2）根据公式（2），相关系数满足|ij1。当ij=0时，Si(t)和Sj(t)不 相 关，相 互 独 立；当0|iji+1=M=2（4）由于信号向量和噪声向量相互独立，可进一步得到R=UUH=USSUSH+UNNUNH，其中S为前P个较大的信号特征值组成的对角阵，N为后MP个较小的噪声特征值组成的对角阵；US为前P个较大的信号特征值对应的信号特征向量张成的信号子空间，UN为后MP个较小噪声特征值对应的噪声特征向量张成的噪声子空间。当入射目标信号源之间相互独立时，方向矩阵A中各列的导向向量a()i与噪声子空间UN正交，可得UHNa()i=0，则MUSIC算法的空间谱函数为：PMUSIC()i=1aH()iUNUHNa()i（5）根据公式（5）得到的空间谱函数，通过谱峰搜索可以估计出入射信号的DOA。但在实际应用中协方差矩阵R无法直接获得，一般情况下采用样本协方差矩阵R来代替R进行特征值分解8。R=1LLi=1y(t)y(t)H（6）在实际情况下，由于存在噪声的干扰导致导向向量a()i与噪声子空间UN不完全正交，aH()iUNUHNa()i结果不严格为0，是一个趋近于0的极小值，此时MUSIC算法的空间谱函数可以表示为：PMUSIC()i=1aH()iUNUHNa()i（7）对以上过程进行总结，得到经典MUSIC算法的主要步骤为：步骤1：由公式（3）得到均匀线阵接收到的信号协方差矩阵R；步骤2：对信号协方差矩阵R进行特征值分解，将得到的特征值按从大到小排序；步骤3：提取与入射信号个数相等的P个较大特征值对应特征向量张成的信号子空间US，提取剩下的MP个较小特征值对应特征向量张成的噪声子空间UN；步骤4：使变化，根据公式（5）计算空间谱函数并进行谱峰搜索，得到入射信号的DOA估计值。3改进型MUSIC算法根据理论分析，经典MUSIC在大快拍数、高信噪比、阵列阵元个数足够的条件下，MUSIC算法能充分发挥作用，可以达到任意高的分辨率。在信号弱相干或不相干条件下，接收到信号的协方差矩阵R满秩，经过特征值分解得到的信号子空间US和噪声子空间UN正交，根据公式（5）的空间谱函数进行谱峰搜索可以达到较高的分辨率。但是实际应用中MUSIC算法的估计性能在某些情况下会失效，一方面可能因为入射信号存在多径传播和同频干扰导致信号相干情况发生，部分入射信号的能量会散发到噪声子空间，造成协方差矩阵R不满秩9；另一方面可能在目标附近存在强反射点时，协方差矩阵R经特征分解后得到的信号子空间中目标信号的信息被强反射点遮盖，导致弱信号强度的目标难以估计10。针对以上问题，在经典MUSIC算法基础上提出一种基于协方差矩阵重构的改进型MUSIC算法。首先为解决信号相干问题，参考空间平滑算法的原理，重构处理接收到信号协方差矩阵，使信号 协 方 差 矩 阵 的 秩 达 到 满 秩 状 态，即rank()R=M；接着为加强对弱信号强度目标的识别，对噪声子空间根据实际的DOA估计场景进行加权修正，得出新的噪声子空间和改进后的空间谱函数；最后进行谱峰搜索得到入射目标信号的DOA估计。假设入射信号目标的个数已知，本文基于经典MUSIC算法提出的改进型MUSIC算法的具体步骤为：步骤1：重构入射信号协方差矩阵：将阵列接收到的数据向量y()t变换为Y()t：Y()t=Jy()t（8）在式（8）中，Y()t为y()t的复共轭矩阵；J为MM维反单位矩阵，即J=001010100MM（9）步骤2：求Y()t的协方差矩阵：R1=EY()t YH()t=JRJ（10）公式（10）中，R为R的复共轭矩阵。步骤3：根据R1和R的算数平均值重构协方差矩阵RY：RY=R1+R2=JRJ+R2（11）步骤4：根据矩阵运算原理，矩阵R，R1和RY具有相同的信号子空间和噪声子空间。对重构后的协方差矩阵RY进行特征值分解：RY=UUH（12）38第 3 期式（12）中，是由全部特征值组成的对角矩阵。步骤5：对噪声子空间进行加权修正，改进后的噪声特征值定义为：i=i+0.1i,i=P+1,P+2,M（13）式（13）中，i为初始噪声的特征值，i为改进后的噪声特征值，0.1为校正值。校正值0.1=0时，校正值不对特征值起到校正作用，改进后的特征值等于经典MUSIC算法的特征值；校正值0.1过大时，会导致噪声子空间的特征值接近信号子空间的特征值，会引起欠估计问题，算法分辨率下降；校正值0.1过小时，无法起到特征值校正作用，噪声子空间的特征值与信号子空间的特征值之间的差值会更大，导致弱信号强度的目标难以估计。因此，需要找到合适的校正值。根据信息论准则并参考束强11的研究得到：噪声特征值的最大值与最小值之比应小于2，此时经典MUSIC算法的性能达到最佳，即：P+1M2（14）将改进后的特征值i=i+0.1i代入式（14），可得：P+1M=P+1+0.1P+1M+0.1M2（15）根据式（15）计算出符合要求的，并取符合条件的最小整数计算得到最终校正值0.1。步骤6：根据最终校正值得到改进后的噪声子空间为：UN=P+1vP+1,P+2vP+2,MvM（16）式（16）中，P+1,P+2,M为改进后的噪声特征值，vP+1,vP+2,vM为经典MUSIC算法得到的噪声特征值对应的特征向量。步骤7：改进后的空间谱函数为：P()i=1aH()iUNUHNa()i（17）步骤8：根据改进后的空间谱函数进行谱峰搜索，得到入射信号的DOA估计。根据理论分析，本文提出的改进型MUSIC算法相比较前向空间平滑（forward spatial smoothing，FSS）算法、前后向空间平滑（forward backward spatial smoothing，FBSS）算法这类空间平滑算法可以在不牺牲阵列孔径的基础上，实现对相干信号的DOA估计。当阵列阵元个数为M时，本文提出的算法最多可以对M1个相干信源目标进行DOA估计。因此本文算法具有更大的阵列孔径，能提高DOA估计的精度。此外，相比较空间平滑类算法需要对构造出的多个子阵的数据协方差进行特征分解，本文提出的改进型MUSIC算法只需要进行两次特征分解，算法复杂程度降低，整个算法的运行时间减少，有着更高的应用价值。4仿真实验与分析为对比经典MUSIC算法、FSS算法、FBSS算法和本文提出的改进型MUSIC算法的性能，在阵元个数M=10、相邻阵元间距d=/2的均匀线阵，噪声为理想高斯白噪声条件下进行仿真实验。对于需要进行空间平滑处理过程的仿真实验，将整个均匀线阵划分为4个子阵，每个子阵包含6个阵元。FSS算法和FBSS算法的具体实现步骤参考张权等12、姚昕彤等13的方法。4.1算法对相距很近相干信源的估计性能有两个相同功率的相干信源分别从0和6方向入射到M=10均匀线阵上，信噪比为5 dB，快拍数为200，仿真结果如图2所示。/KAE#图2近距离相干信源的DOA估计由图2可知，当两个相干目标以6间隔入射到均匀线阵时，经典MUSIC算法已经失效，无法对相干目标进行识别；FSS算法和FBSS算法虽然有两个峰值点，且FBSS算法得到的尖峰比FSS算法得到的尖峰更尖锐，但这两种算法得到的DOA估计值明显与真实入射角度值之间存在偏差；本文提出的算法形成的尖峰更加尖锐，DOA估计值与真实入射角度值之间的差值更小，说明本文算法对入射目标的角度分辨能力更强，有着更好的DOA估计性能。屠亚杰，方遒，李艳玲：基于协方差矩阵重构的改进型MUSIC算法39淮阴工学院学报2023 年4.2算法对多个相干信源的估计性能由于经典MUSIC算法无法对相干信号进行识别，因此仅对比FSS算法、FBSS算法和本文提出的改进型MUSIC算法在多个相干信源条件下的DOA估计性能。有 6 个相同功率的相干信号分别从-6，0，6，18，30，40方向入射到M=10的均匀线阵，信噪比均为10 dB，快拍数为200，仿真结果如图3所示。图3多相干信源条件下的DOA估计唐晓杰等14研究表明，在阵元个数为M的均匀线阵条件下，FSS算法最多只可以对M/2个相干信源进行DOA估计；FBSS算法因为额外利用了后向子阵列，增加了平滑次数，最多可以对2M/3个相干信源目标进行DOA估计。因此，在本研究设置的仿真条件下，FSS算法最多可以对5个相干信源进行DOA估计；FBSS算法最多可以对6个相干信源进行DOA估计。从图3仿真结果可以看到，6个相干信号入射时，FSS算法只产生5个峰值点，无法完成对全部6个入射角度的DOA估计，FSS算法失效；FBSS算法虽然产生6个尖锐的峰值点，但生成的6个峰值点中只有4个峰值点在真实峰值点附近，分别位于-6、18、30、40这4个真实峰值点附近，其它两个峰值点为产生的虚假峰值点，FBSS算法也失效。本文提出的改进型MUSIC算法在阵元个数为10时，最多可以对9个相干信源目标进行DOA估计，所以该算法可以对这6个相干信号进行DOA估计。仿真结果也表明，本文算法可以完成对这6个入射目标的DOA估计，且峰值点都在真实角度附近，产生的峰值点与真实峰值点之间的误差不超过0.5。验证了本文提出的改进型MUSIC算法相较FSS算法和FBSS算法在相同快照数、信噪比条件下，在保证DOA估计精度的同时能利用更少的阵元个数完成更多的入射信号的DOA估计。4.3不同信噪比条件下的算法性能比较已知均方根误差（RMSE）公式为：RMSE=1WPPi=1Ww=1()i,wi2。（18）式（18）中，W表示蒙特卡洛实验重复的次数，P是信号源的个数，i,w是在第w次蒙特卡洛实验中对第i个信号源i的DOA估计值。有两个相干信号分别从 0、10方向入射到M=10的均匀线阵上，快拍数为 200。信噪比在1010 dB 范围内变化，步长为 2 dB，根据式（18）在不同的信噪比条件下进行100次蒙特卡洛实验，得出FSS算法、FBSS算法和本文提出算法的DOA估计均方根误差随着信噪比的变化曲线，如图4所示。图4算法RMSE随信噪比变化曲线根据图4可知，在信噪比大于0 dB时，三种算法的DOA估计均方根误差十分接近于0，说明此时这3种算法的DOA估计性能都很好，3种算法都可以用于相干信源目标的DOA估计；当信噪比小于0 dB时，特别是SNR=10 dB时可以明显发现本文提出算法的均方根误差要远小于FSS算法和FBSS算法的均方根误差。说明在低信噪比条件下，本文提出算法的均方根误差更小，本文提出的算法有着更好的DOA估计性能。4.4不同快拍数条件下的算法性能比较有两个相干信号分别从 0、10方向入射到M=10的均匀线阵上，信噪比设置为8 dB。快拍数在20200范围内变化，步长为20，根据式（18）在不同的快拍数条件下进行100次蒙特卡洛实验，得出FSS算法、FBSS算法和本文提出算法的DOA40第 3 期估计均方根误差随着快拍数的变化曲线，如图5所示。图5算法RMSE随快拍数变化曲线由图5可知，FSS算法、FBSS算法和本文提出算法的DOA估计均方根误差均随着快拍数的增加而逐渐降低，即这3种算法的估计精度会随着快拍数的增加而提高，并且随着快拍数增加到足够大时，这3种算法的估计精度将会保持同一水平，即在大快拍数条件下这3种算法都能对入射信号进行准确的DOA估计。在快拍数小于140的小快拍条件下，本文提出算法的估计精度优于FSS算法和FBSS算法。当快拍数为60时，FSS算法和FBSS算法的均方根误差均大于 1，分别为 1.10和1.18；此时本文提出算法的均方根误差仅为0.93。可以说明在小快拍数条件下，本文提出的算法相比FSS算法和FBSS算法有着更高的角度估计精度；并且在实验过程中还发现，本文提出的算法在小快拍条件下对相干信号DOA估计的成功率高于FSS算法和FBSS算法，具备更好的DOA估计性能。5结论针对经典 MUSIC算法无法对相干信源进行DOA估计的问题，提出了一种基于协方差矩阵重构的改进型MUSIC算法。该算法在对信号协方差矩阵进行重构处理的基础上，对其噪声特征值进行加权处理，得到新的噪声子空间和改进后的空间谱函数，结合谱峰搜索实现同等信号强度信源邻近情况下相干信源的DOA估计。经过仿真验证，本文提出算法可以在不牺牲阵列有效孔径的基础上完成对相干信号的DOA估计，并且在低信噪比、小快拍数条件下，相比FSS算法和FBSS算法有着更低的均方根误差，有着更好的DOA估计性能和最少的计算量，更能满足实际工程应用中对算法高准确性、易于实现的要求。参考文献：1刘琳.稀疏阵列的DOA估计研究 D.南京：南京理工大学，2019.2PATOLE S，TORLAK M，WANG D，et al.Automotive radars a review of signal processing techniquesJ.IEEE Signal Processing Magazine，2017，34（2）：22-35.3同非，郭磊，连豪.基于MUSIC及其改进算法的空间信号到达方向估计方法研究 J.火控雷达技术，2018，47（2）：58-62.4 张小飞.阵列信号处理及MATLAB实现 M.北京：电子工业出版社，2019.5石要武，陈淼，单泽涛，等.基于特征空间MUSIC算法的相干信号波达方向空间平滑估计 J.吉林大学学报（工学版），2017，47（1）：268-273.6刘晓志，宋牧野，李鸿儒.基于伪数据相关矩阵二次重构的DOA估计新算法 J.东北大学学报（自然科学版），2018，39（6）：761-765，770.7张贺勇，王新宇，杨杰，等.基于空间平滑的单次快拍DOA 估计方法 J.火力与指挥控制，2021，46（6）：105-109，114.8黎昕婷，钟舜聪，钟剑锋.基于改进MUSIC算法的宽带信号 DOA 估计 J.计算机工程，2022，48（11）：201-206.9HU B，XUE J L，CHONG P.Denoised modified incoherent signal subspace method for DOA of coherent signalsC/IEEE Asia-Pacific Conference onAntennas and Propagation，2018：539-540.10蒋繁.基于压缩感知的探地雷达信号时延估计 J.舰船电子工程，2020，40（4）：140-143.11束强.基于软件无线电的无线电测向算法研究 D.成都：电子科技大学，2019.12张权，李思敏，唐智灵.基于混合型MUSIC算法对相干信源 DOA 估计 J.计算机应用研究，2020，37（5）：1536-1540.13姚昕彤，王玉文，刘奇，等.基于MUSIC及其改进算法的 DOA 估计研究 J.通信技术，2021，54（6）：1363-1369.14唐晓杰，赵迪，何明浩，等.一种基于前后向Toeplitz矩阵重构的相干信号DOA估计算法 J.中国电子科学研究院学报，2019，14（8）：823-829.（责任编辑：王萍）屠亚杰，方遒，李艳玲：基于协方差矩阵重构的改进型MUSIC算法41