一维DOA估计之Capon算法(含程序)

一、数据模型

   考虑一个由M个传感器构成的阵列k个窄带信号源激励。那么M×1维传感器阵列输出向量X可用下式表示：

                                                   X(t)=As(t)+e(t)

其中，s(t)是在一定参考点测量的K×1维源信号向量，e(t)是加性噪声，

A=[a(θ1)，a(θ2)........,a(θk)]       

a(θk) =[1,exp(-j*2*pi*(d/λ)*sin(θk),.........,exp(-j*2*pi*(M-1)*(d/λ)*sin(θk)]T

其中的0：d：（M-1）*d是指均匀线阵排列的传感器的位置，λ指信号波长

二、Capon算法

Capon算法的DOA估计之由以下函数取极小值时的角度所决定

                              

其中，θ表示DOA变量,R是协方差矩阵

                                      R=(X*X')/N

N是指快拍数

算法步骤：

1.先计算协方差矩阵

2，对1/(a’*R*a）进行谱峰搜索，其峰值即为DOA估计

在这上面打公式太过费劲，下次还是WORD打完截图

MATLAB程序：https://download.csdn.net/download/shaoyede/19520454?spm=1001.2014.3001.5503