谱估计ODA，MUSIC，ESPRIT等ＤＯＡ算法

内含5个matlab文件，都是自己编写，注释很清楚，适合小白学习理解。内含capon_MVM算法，ESPRIT算法，MUSIC，MUSIC平滑，以及性能比较

%capon和CBF波束形成
clear all;close all;clc
M=16;         %均匀线阵列数目
P=2;          %信号源数目
f0=2500;f1=1000;%信号频率Hz
nn=2000;      %快拍数
angle1=-10;angle2=30;%信号入射角度
th=[angle1;angle2]‘;
SN1=-40;SN2=-40;%信噪比
sn=[SN1;SN2];
degrad=pi/180;
tt=0:1/2000:（1-1/2000）;
x0=exp（-j*2*pi*f0*tt）; %2个信号x0、x1
x1=exp（-j*2*pi*f1*tt）; 
t=1:nn;
S=[x0（t）;x1（t）];
nr=randn（Mnn）;
ni=randn（Mnn）;
u=nr+j*ni;      %复高斯白噪声
pn=2 ;          %复高斯噪声功率
tmp=sqrt（pn*10.^（sn/10））;%信号功率
S2=diag（tmp）*S+u（[1 2]:）; %计算信号幅值，满足所设信噪比
AA1=[exp（-pi*j*[0:M-1]*sin（angle1*pi/180））]‘;%阵列流型;系数的设置对结果影响非常大%
AA2=[exp（-pi*j*[0:M-1]*sin（angle2*pi/180））]‘;%估计是涉及到周期的问题            %
A=[AA1AA2];
X=A*S2;  %接收到的信号
Rxx=X*X‘./nn;  %相关矩阵
invRxx=inv（Rxx）;
%信号搜索
th2=[-pi/2:pi/180:pi/2];
for n=1:length（th2）
    a（n）=th2（n）;
    A2=[exp（-1.1*pi*j*[0:M-1]*sin（a（n）））]‘;%搜索空间的方向向量表示                              
    den=A2‘*invRxx*A2;
    P_capon（n）=abs（1./den）;%capon波束形成
    P_cbf（n）=A2‘*Rxx*A2;%CBF波束形成
end;
semilogy（th2*180/piP_capon‘r‘）;
hold on;
semilogy（th2*180/piP_cbf‘b‘）;
axis（[-90 90 1e-4 1e5]）;
grid on;
xlabel（‘角度‘）;ylabel（‘P功率谱‘）;
title（‘SNR1=-40 SNR2=20 th1=-10 th2=30‘）;
legend（‘capon‘‘CBF‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1325  2002-01-27 21:17  capon_MVM.m
     文件         756  2015-05-03 23:10  esprit.m
     文件        1010  2015-04-21 20:02  music.m
     文件        2096  2015-03-31 21:02  music_capon_practice.m
     文件         919  2015-04-21 22:38  music_pinghua.m