【阵列信号处理】DOA估计算法

DOA估计中的Capon、MUSIC算法

最近也在上阵列信号处理的课程，目前学习的主要内容包括阵列的信号模型、方向图、空时等效、Capon波束形成器、MMSE波束形成器、LCMV波束形成器、GSC波束形成器、传统测向法(比相法、CBF/Bartlett估计器)、MUSIC算法、ESPRIT算法以及空间平滑法。下面主要介绍DOA估计算法，主要内容有：MUSIC算法及改进、MUSIC及Capon法的比较、LS-ESPRIT以及TLS-ESPRIT。本文将介绍前两部分。

DOA估计

时域频谱表示信号在各个频率上的能量分布；空间谱表示信号在空间各个方向上的能量分布。所以如果能够得到信号的空间谱，就能够得到信号的波达方向 (Direction of arrival, DOA) 所以空间谱估计也被称为DOA估计。

1. DOA估计的分类主要有多信号源DOA、相干源DOA、宽带信号源DOA以及复杂环境下的DOA…
2. DOA估计的方法主要包括传统法、子空间法、最大似然法和综合法。传统法存在瑞利限；最大似然法高分辨最优但计算量大；子空间法属于次最优算法；综合法也会用到子空间法。过多的就不再展开赘述了，感兴趣者可自行查阅文献。

瑞利限: 一般的，在一定阵列长度下可达到的最小分辨率称为瑞利限，超过瑞利限的方法被称为 超分辨方法。

• 1979年，Schmidt提出MUSIC算法，1986年重新发表
• 1986年，Roy提出ESPRIT算法

以上两种方法属于早期经典的超分辨方法，今天大多数子空间法都是在这两种方法基础上改进的。

MUSIC算法

算法假设：

1.均匀线阵且阵元间距不大于处理最高频率信号波长的二分之一
2.处理器噪声为加性高斯分布，不同阵元间距噪声均为平稳随机过程，独立同分布，各阵元噪声方差相同
3.信号为零均值平稳随机过程且与噪声独立
4.阵元数大于信源数，信号源为窄带信号且入射信源之间弱相关或者不相关。
##窄带：信号通过天线阵列的时间远远小于信号带宽的倒数
1
2
3
4
5

MUSIC算法步骤

1.设置信源、阵列模型、接收信号
2.计算接收信号的协方差矩阵
3.协方差矩阵特征分解，得到特征值
4.特征值由大到小排序，求对应的信号子空间和噪声子空间
5.利用计算MUSIC谱，即噪声功率谱平方和的倒数
6.搜索MUSIC谱的峰值，找到最大的K个谱峰，对应的K个DOA估计的角度
1
2
3
4
5
6

• 计算MUSIC谱的公式如下

$$P_{music}=\frac{1}{\sum _{i=P+1}^N{\left|v_i^{H}a(\theta )\right|}^2}$$，其中N为阵元数，P为信号子空间数，$v_i$为特征向量。

• Capon方法的计算公式如下
  $$P_{capon}=\frac{1}{a^H(\theta )R_x^{-1}a(\theta )}$$
  其中$R_x^{-1}$为接收信号协方差的逆矩阵。

下面通过仿真来对比一下MUSIC与Capon的区别，参考的文章中给出了不同信噪比下以及不同角度间隔下两种方法的估计误差。
1. 信噪比
设置阵元数为10，阵元间隔为半波长，信源数为3(-10度,0度,20度)，快拍数为1024，下图为估计得到的信号谱，左侧信噪比设置为-8dB，右侧信噪比设置为10dB。

从上图可以看出，在信噪比较低时两种估计算法性能都急剧下降，但是MUSIC算法略优于Capon，而当信噪比较大时，两种算法基本一致。

关于谱峰强度的解释：

 MUSIC谱峰只反映阵列流形矢量与噪声子空间的正交性，与信噪比无关；
 Capon谱峰是真正的输出功率，与信噪比有关。
1
2

2.角度间隔/分辨率
设置阵元数为10，阵元间隔为半波长，信源数为3(-0.5度,0度,0.5度)，快拍数为1024，信噪比设置为10dB，下图为估计得到的信号谱，为方便观察进行了归一化。
可以看到这种情况下，MUSIC的分辨率是优于Capon法的。
上述仿真分析感觉也不是特别严谨，至于两种方法谁好谁坏，在处理具体问题的时候都试一下哪个好就用哪个吧